-
Can massive stars form in low mass clouds?
Authors:
Jamie D. Smith,
Sarah E. Jaffa,
Martin G. H. Krause
Abstract:
The conditions required for massive star formation are debated, particularly whether massive stars must form in conjunction with massive clusters. Some authors have advanced the view that stars of any mass (below the total cluster mass) can form in clusters of any mass with some probability (random sampling). Others pointed out that the scatter in the determinations of the most massive star mass f…
▽ More
The conditions required for massive star formation are debated, particularly whether massive stars must form in conjunction with massive clusters. Some authors have advanced the view that stars of any mass (below the total cluster mass) can form in clusters of any mass with some probability (random sampling). Others pointed out that the scatter in the determinations of the most massive star mass for a given cluster mass was consistent with the measurement error, such that the mass of the most massive star was determined by the total cluster mass (optimal sampling). Here we investigate the relation between cluster mass (M\textsubscript{ecl}) and the maximum stellar mass (M\textsubscript{max}) using a suite of SPH simulations. Varying cloud mass and turbulence random seed results in a range of cluster masses which we compare with their respective maximum star masses. We find that more massive clusters will have, on average, higher mass stars with this trend being steeper at lower cluster masses ($M\textsubscript{max} \propto M\textsubscript{ecl}^{0.31}$ for $M\textsubscript{ecl}<500M\,_{\odot}$) and flattening at higher cluster masses ($M\textsubscript{max} \propto M\textsubscript{ecl}^{0.11}$ for $M\textsubscript{ecl}>500M\,_{\odot}$). This rules out purely stochastic star formation in our simulations. Significant scatter in the maximum masses with identical initial conditions also rules out the possibility that the relation is purely deterministic (that is that a given cluster mass will result in a specific maximum stellar mass). In conclusion our simulations disagree with both random and optimal sampling of the initial mass function.
△ Less
Submitted 28 September, 2023; v1 submitted 11 September, 2023;
originally announced September 2023.
-
RJ-plots: An improved method to classify structures objectively
Authors:
Seamus D. Clarke,
Sarah E. Jaffa,
Anthony P. Whitworth
Abstract:
The interstellar medium is highly structured, presenting a range of morphologies across spatial scales. The large data sets resulting from observational surveys and state-of-the-art simulations studying these hierarchical structures means that identification and classification must be done in an automated fashion to be efficient. Here we present RJ-plots, an improved version of the automated morph…
▽ More
The interstellar medium is highly structured, presenting a range of morphologies across spatial scales. The large data sets resulting from observational surveys and state-of-the-art simulations studying these hierarchical structures means that identification and classification must be done in an automated fashion to be efficient. Here we present RJ-plots, an improved version of the automated morphological classification technique J-plots developed by Jaffa et al. (2018). This method allows clear distinctions between quasi-circular/elongated structures and centrally over/under-dense structures. We use the recent morphological SEDIGISM catalogue of Neralwar et al. (2022) to show the improvement in classification resulting from RJ-plots, especially for ring-like and concentrated cloud types. We also find a strong correlation between the central concentration of a structure and its star formation efficiency and dense gas fraction, as well as a lack of correlation with elongation. Furthermore, we use the accreting filament simulations of Clarke et al. (2020) to highlight a multi-scale application of RJ-plots, finding that while spherical structures become more common at smaller scales they are never the dominant structure down to $r\sim0.03$ pc.
△ Less
Submitted 15 August, 2022;
originally announced August 2022.
-
Star Cluster Formation in Clouds with Externally Driven Turbulence
Authors:
Jamie D. Smith,
James E. Dale,
Sarah E. Jaffa,
Martin G. H. Krause
Abstract:
Star clusters are known to be formed in turbulent molecular clouds. How turbulence is driven in molecular clouds and what effect this has on star formation is still unclear. We compare a simulation setup with turbulent driving everywhere in a periodic box with a setup where turbulence is only driven around the outside of the box. We analyse the resulting gas distribution, kinematics, and the popul…
▽ More
Star clusters are known to be formed in turbulent molecular clouds. How turbulence is driven in molecular clouds and what effect this has on star formation is still unclear. We compare a simulation setup with turbulent driving everywhere in a periodic box with a setup where turbulence is only driven around the outside of the box. We analyse the resulting gas distribution, kinematics, and the population of stars that are formed from the cloud. Both setups successfully produce a turbulent velocity field with a power law structure function, the externally driven cloud has a more central, monolithic, clump, while the fully driven cloud has many smaller, more dispersed, clumps. The star formation follows the cloud morphology producing large clusters, with high star forming efficiency in the externally driven simulations and sparse individual star formation with much lower star formation efficiency in the fully driven case. We conclude that the externally driven method, which resembles a Global Hierarchical Collapse (GHC) scenario, produces star clusters that more closely match with observations.
△ Less
Submitted 10 August, 2022;
originally announced August 2022.
-
Grouped star formation: converting sink particles to stars in hydrodynamical simulations
Authors:
K. Y. Liow,
S. Rieder,
C. L. Dobbs,
S. E. Jaffa
Abstract:
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. How…
▽ More
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.
△ Less
Submitted 9 December, 2021;
originally announced December 2021.
-
A simple approach to CO cooling in molecular clouds
Authors:
A. P. Whitworth,
S. E. Jaffa
Abstract:
CO plays an important role in interstellar molecular clouds, both as a coolant, and as a diagnostic molecule. However, a proper evaluation of the cooling rate due to CO requires a determination of the populations of many levels, the spontaneous and stimulated radiative de-excitation rates between these levels, and the transfer of the emitted multi-line radiation; additionally, this must be done fo…
▽ More
CO plays an important role in interstellar molecular clouds, both as a coolant, and as a diagnostic molecule. However, a proper evaluation of the cooling rate due to CO requires a determination of the populations of many levels, the spontaneous and stimulated radiative de-excitation rates between these levels, and the transfer of the emitted multi-line radiation; additionally, this must be done for three isotopologues. It would be useful to have a simple analytic formulation that avoided these complications; this could then be used in situations where CO plays an important role as a coolant, but the details of this role are not the main concern. We derive such a formulation here, by first considering the two asymptotic forms that obtain in the limits of (a) low volume-density and optical depth, and (b) high volume-density and optical depth. These forms are then combined in such a way as to fit the detailed numerical results from Goldsmith \& Langer (1978), which cover low temperatures, and a range of physical conditions where the interplay of thermal and sub-thermal excitation, optical-depth effects, and the contributions from rare isotopologues, are all important. The fit is obtained using the Metropolis-Hastings method, and reproduces the results of GL78 well. It is a purely local and analytic function of state --- specifically a function of the density, $ρ$, isothermal sound speed, $a$, CO abundance, $X_{_{\rm CO}}$, and velocity divergence, $\nabla\cdot{\boldsymbol\upsilon}$. As an application, we consider the cooling layer following a slow steady non-magnetic planar J-shock. We show that, if the post-shock cooling is dominated by CO and its isotopologues, the thickness of the post-shock cooling layer is very small and approximately independent of the pre-shock velocity, $\upsilon_o$, or pre-shock isothermal sound speed, $a_o$.
△ Less
Submitted 16 November, 2018;
originally announced November 2018.
-
Multifractal Analysis of the Interstellar Medium. First application to Hi-GAL Observations
Authors:
D. Elia,
F. Strafella,
S. Dib,
N. Schneider,
P. Hennebelle,
S. Pezzuto,
S. Molinari,
E. Schisano,
S. E. Jaffa
Abstract:
The multifractal geometry remains an under-exploited approach to describe and quantify the large-scale structure of interstellar clouds. In this paper, the typical tools of multifractal analysis are applied to Herschel far-infrared (70-500 $μ$m) dust continuum maps, which represent an ideal case of study. This dust component is a relatively optically thin tracer at these wavelengths and the size i…
▽ More
The multifractal geometry remains an under-exploited approach to describe and quantify the large-scale structure of interstellar clouds. In this paper, the typical tools of multifractal analysis are applied to Herschel far-infrared (70-500 $μ$m) dust continuum maps, which represent an ideal case of study. This dust component is a relatively optically thin tracer at these wavelengths and the size in pixel of the maps is well suitable for this statistical analysis. We investigate the so-called multifractal spectrum and generalised fractal dimension of six Hi-GAL regions in the third Galactic quadrant. We measure systematic variations of the spectrum at increasing wavelength, which generally correspond to an increasing complexity of the image, and we observe peculiar behaviours of the investigated fields, strictly related to the presence of high-emission regions, which in turn are connected to star formation activity. The same analysis is applied to synthetic column density maps, generated from numerical turbulent molecular cloud models and from fractal Brownian motion (fBm), allowing for the confrontation of the observations with models with well controlled physical parameters. This comparison shows that, in terms of multifractal descriptors of the structure, fBm images exhibit a very different, quite unrealistic behaviour when compared with Hi-GAL observations, whereas the numerical simulations appear in some cases (depending on the specific model) more similar to the observations. Finally, the link between mono-fractal parameters (commonly found in the literature) and multifractal indicators is investigated: the former appear to be only partially connected with the latter, with the multifractal analysis offering a more extended and complete characterization of the cloud structure.
△ Less
Submitted 8 August, 2018;
originally announced August 2018.
-
Synthetic C$^{18}$O observations of fibrous filaments: the problems of mapping from PPV to PPP
Authors:
S. D. Clarke,
A. P. Whitworth,
R. L. Spowage,
A. Duarte-Cabral,
S. T. Suri,
S. E. Jaffa,
S. Walch,
P. C. Clark
Abstract:
Molecular-line observations of filaments in star-forming regions have revealed the existence of elongated coherent features within the filaments; these features are termed fibres. Here we caution that, since fibres are traced in PPV space, there is no guarantee that they represent coherent features in PPP space. We illustrate this contention using simulations of the growth of a filament from a tur…
▽ More
Molecular-line observations of filaments in star-forming regions have revealed the existence of elongated coherent features within the filaments; these features are termed fibres. Here we caution that, since fibres are traced in PPV space, there is no guarantee that they represent coherent features in PPP space. We illustrate this contention using simulations of the growth of a filament from a turbulent medium. Synthetic C$^{18}$O observations of the simulated filaments reveal the existence of fibres very similar to the observed ones, i.e. elongated coherent features in the resulting PPV data-cubes. Analysis of the PPP data-cubes (i.e. 3D density fields) also reveals elongated coherent features, which we term sub-filaments. Unfortunately there is very poor correspondence between the fibres and the sub-filaments in the simulations. Both fibres and sub-filaments derive from inhomogeneities in the turbulent accretion flow onto the main filament. As a consequence, fibres are often affected by line-of-sight confusion. Similarly, sub-filaments are often affected by large velocity gradients, and even velocity discontinuities. These results suggest that extreme care should be taken when using velocity coherent features to constrain the underlying substructure within a filament.
△ Less
Submitted 25 June, 2018; v1 submitted 22 June, 2018;
originally announced June 2018.
-
J plots: a new method for characterizing structures in the interstellar medium
Authors:
S. E. Jaffa,
A. P. Whitworth,
S. D. Clarke,
A. D. P. Howard
Abstract:
Large scale surveys have brought about a revolution in astronomy. To analyse the resulting wealth of data, we need automated tools to identify, classify, and quantify the important underlying structures. We present here a method for classifying and quantifying a pixelated structure, based on its principal moments of inertia. The method enables us to automatically detect, and objectively compare, c…
▽ More
Large scale surveys have brought about a revolution in astronomy. To analyse the resulting wealth of data, we need automated tools to identify, classify, and quantify the important underlying structures. We present here a method for classifying and quantifying a pixelated structure, based on its principal moments of inertia. The method enables us to automatically detect, and objectively compare, centrally condensed cores, elongated filaments and hollow rings. We illustrate the method by applying it to (i) observations of surface-density from Hi-GAL, and (ii) simulations of filament growth in a turbulent medium. We limit the discussion here to 2D data; in a future paper we will extend the method to 3D data.
△ Less
Submitted 5 March, 2018;
originally announced March 2018.
-
Properties of Hi-GAL clumps in the inner Galaxy]{The Hi-GAL compact source catalogue. I. The physical properties of the clumps in the inner Galaxy ($-71.0^{\circ}< \ell < 67.0^{\circ}$)
Authors:
D. Elia,
S. Molinari,
E. Schisano,
M. Pestalozzi,
S. Pezzuto,
M. Merello,
A. Noriega-Crespo,
T. J. T. Moore,
D. Russeil,
J. C. Mottram,
R. Paladini,
F. Strafella,
M. Benedettini,
J. P. Bernard,
A. Di Giorgio,
D. J. Eden,
Y. Fukui,
R. Plume,
J. Bally,
P. G. Martin,
S. E. Ragan,
S. E. Jaffa,
F. Motte,
L. Olmi,
N. Schneider
, et al. (61 additional authors not shown)
Abstract:
Hi-GAL is a large-scale survey of the Galactic plane, performed with Herschel in five infrared continuum bands between 70 and 500 $μ$m. We present a band-merged catalogue of spatially matched sources and their properties derived from fits to the spectral energy distributions (SEDs) and heliocentric distances, based on the photometric catalogs presented in Molinari et al. (2016a), covering the port…
▽ More
Hi-GAL is a large-scale survey of the Galactic plane, performed with Herschel in five infrared continuum bands between 70 and 500 $μ$m. We present a band-merged catalogue of spatially matched sources and their properties derived from fits to the spectral energy distributions (SEDs) and heliocentric distances, based on the photometric catalogs presented in Molinari et al. (2016a), covering the portion of Galactic plane $-71.0^{\circ}< \ell < 67.0^{\circ}$. The band-merged catalogue contains 100922 sources with a regular SED, 24584 of which show a 70 $μ$m counterpart and are thus considered proto-stellar, while the remainder are considered starless. Thanks to this huge number of sources, we are able to carry out a preliminary analysis of early stages of star formation, identifying the conditions that characterise different evolutionary phases on a statistically significant basis. We calculate surface densities to investigate the gravitational stability of clumps and their potential to form massive stars. We also explore evolutionary status metrics such as the dust temperature, luminosity and bolometric temperature, finding that these are higher in proto-stellar sources compared to pre-stellar ones. The surface density of sources follows an increasing trend as they evolve from pre-stellar to proto-stellar, but then it is found to decrease again in the majority of the most evolved clumps. Finally, we study the physical parameters of sources with respect to Galactic longitude and the association with spiral arms, finding only minor or no differences between the average evolutionary status of sources in the fourth and first Galactic quadrants, or between "on-arm" and "inter-arm" positions.
△ Less
Submitted 4 June, 2017;
originally announced June 2017.
-
$\mathcal{Q}^{+}$: Characterising the structure of young star clusters
Authors:
S. E. Jaffa,
A. P. Whitworth,
O. Lomax
Abstract:
Many young star clusters appear to be fractal, i.e. they appear to be concentrated in a nested hierarchy of clusters within clusters. We present a new algorithm for statistically analysing the distribution of stars to quantify the level of sub-structure. We suggest that, even at the simplest level, the internal structure of a fractal cluster requires the specification of three parameters. (i) The…
▽ More
Many young star clusters appear to be fractal, i.e. they appear to be concentrated in a nested hierarchy of clusters within clusters. We present a new algorithm for statistically analysing the distribution of stars to quantify the level of sub-structure. We suggest that, even at the simplest level, the internal structure of a fractal cluster requires the specification of three parameters. (i) The 3D fractal dimension, $\mathcal{D}$, measures the extent to which the clusters on one level of the nested hierarchy fill the volume of their parent cluster. (ii) The number of levels, $\mathcal{L}$, reflects the finite ratio between the linear size of the large root-cluster at the top of the hierarchy, and the smallest leaf-clusters at the bottom of the hierarchy. (iii) The volume-density scaling exponent, $\mathcal{C}=-\textrm{d}\ln[δn]/\textrm{d}\ln[L]$ measures the factor by which the excess density, $δn$, in a structure of scale $L$, exceeds that of the background formed by larger structures; it is similar, but not exactly equivalent, to the exponent in Larson's scaling relation between density and size for molecular clouds. We describe an algorithm which can be used to constrain the values of $({\cal D},{\cal L},{\cal C})$ and apply this method to artificial and observed clusters. We show that this algorithm is able to reliably describe the three dimensional structure of an artificial star cluster from the two dimensional projection, and quantify the varied structures observed in real and simulated clusters.
△ Less
Submitted 7 December, 2016; v1 submitted 30 November, 2016;
originally announced November 2016.
-
Star Formation triggered by cloud-cloud collisions
Authors:
S. K. Balfour,
A. P. Whitworth,
D. A. Hubber,
S. E. Jaffa
Abstract:
We present the results of SPH simulations in which two clouds, each having mass $M_{_{\rm{o}}}\!=\!500\,{\rm M}_{_\odot}$ and radius $R_{_{\rm{o}}}\!=\!2\,{\rm pc}$, collide head-on at relative velocities of $Δv_{_{\rm{o}}} =2.4,\;2.8,\;3.2,\;3.6\;{\rm and}\;4.0\,{\rm km}\,{\rm s}^{-1}$. There is a clear trend with increasing $Δv_{_{\rm{o}}}$. At low $Δv_{_{\rm{o}}}$, star formation starts later,…
▽ More
We present the results of SPH simulations in which two clouds, each having mass $M_{_{\rm{o}}}\!=\!500\,{\rm M}_{_\odot}$ and radius $R_{_{\rm{o}}}\!=\!2\,{\rm pc}$, collide head-on at relative velocities of $Δv_{_{\rm{o}}} =2.4,\;2.8,\;3.2,\;3.6\;{\rm and}\;4.0\,{\rm km}\,{\rm s}^{-1}$. There is a clear trend with increasing $Δv_{_{\rm{o}}}$. At low $Δv_{_{\rm{o}}}$, star formation starts later, and the shock-compressed layer breaks up into an array of predominantly radial filaments; stars condense out of these filaments and fall, together with residual gas, towards the centre of the layer, to form a single large-$N$ cluster, which then evolves by competitive accretion, producing one or two very massive protostars and a diaspora of ejected (mainly low-mass) protostars; the pattern of filaments is reminiscent of the hub and spokes systems identified recently by observers. At high $Δv_{_{\rm{o}}}$, star formation occurs sooner and the shock-compressed layer breaks up into a network of filaments; the pattern of filaments here is more like a spider's web, with several small-$N$ clusters forming independently of one another, in cores at the intersections of filaments, and since each core only spawns a small number of protostars, there are fewer ejections of protostars. As the relative velocity is increased, the {\it mean} protostellar mass increases, but the {\it maximum} protostellar mass and the width of the mass function both decrease. We use a Minimal Spanning Tree to analyse the spatial distributions of protostars formed at different relative velocities.
△ Less
Submitted 17 September, 2015;
originally announced September 2015.