-
NEATH III: a molecular line survey of a simulated star-forming cloud
Authors:
F. D. Priestley,
P. C. Clark,
S. C. O. Glover,
S. E. Ragan,
O. Fehér,
L. R. Prole,
R. S. Klessen
Abstract:
We present synthetic line observations of a simulated molecular cloud, utilising a self-consistent treatment of the dynamics and time-dependent chemical evolution. We investigate line emission from the three most common CO isotopologues ($^{12}$CO, $^{13}$CO, C$^{18}$O) and six supposed tracers of dense gas (NH$_3$, HCN, N$_2$H$^+$, HCO$^+$, CS, HNC). Our simulation produces a range of line intens…
▽ More
We present synthetic line observations of a simulated molecular cloud, utilising a self-consistent treatment of the dynamics and time-dependent chemical evolution. We investigate line emission from the three most common CO isotopologues ($^{12}$CO, $^{13}$CO, C$^{18}$O) and six supposed tracers of dense gas (NH$_3$, HCN, N$_2$H$^+$, HCO$^+$, CS, HNC). Our simulation produces a range of line intensities consistent with that observed in real molecular clouds. The HCN-to-CO intensity ratio is relatively invariant with column density, making HCN (and chemically-similar species such as CS) a poor tracer of high-density material in the cloud. The ratio of N$_2$H$^+$ to HCN or CO, on the other hand, is highly selective of regions with densities above $10^{22} \, {\rm cm^{-2}}$, and the N$_2$H$^+$ line is a very good tracer of the dynamics of high volume density ($>10^4 \, {\rm cm^{-3}}$) material. Focusing on cores formed within the simulated cloud, we find good agreement with the line intensities of an observational sample of prestellar cores, including reproducing observed CS line intensities with an undepleted elemental abundance of sulphur. However, agreement between cores formed in the simulation, and models of isolated cores which have otherwise-comparable properties, is poor. The formation from and interaction with the large-scale environment has a significant impact on the line emission properties of the cores, making isolated models unsuitable for interpreting observational data.
△ Less
Submitted 10 June, 2024;
originally announced June 2024.
-
Self-consistent modelling of the Milky Way structure using live potentials
Authors:
Eva Durán-Camacho,
Ana Duarte-Cabral,
Alex R. Pettitt,
Robin G. Treß,
Paul C. Clark,
Ralf S. Klessen,
Kamran R. J. Bogue,
Rowan J. Smith,
Mattia C. Sormani
Abstract:
To advance our understanding of the evolution of the interstellar medium (ISM) of our Galaxy, numerical models of Milky Way (MW) type galaxies are widely used. However, most models only vaguely resemble the MW (e.g. in total mass), and often use imposed analytic potentials (which cannot evolve dynamically). This poses a problem in asserting their applicability for the interpretation of observation…
▽ More
To advance our understanding of the evolution of the interstellar medium (ISM) of our Galaxy, numerical models of Milky Way (MW) type galaxies are widely used. However, most models only vaguely resemble the MW (e.g. in total mass), and often use imposed analytic potentials (which cannot evolve dynamically). This poses a problem in asserting their applicability for the interpretation of observations of our own Galaxy. The goal of this work is to identify a numerical model that is not only a MW-type galaxy, but one that can mimic some of the main observed structures of our Galaxy, using dynamically evolving potentials, so that it can be used as a base model to study the ISM cycle in a galaxy like our own. This paper introduces a suite of 15 MW-type galaxy models developed using the {\sc arepo} numerical code, that are compared to Galactic observations of $^{12}$CO and \ion{H}{I} emission via longitude-velocity plots, from where we extract and compare the skeletons of major galactic features and the terminal gas velocities. We found that our best-fitting model to the overall structure, also reproduces some of the more specific observed features of the MW, including a bar with a pattern speed of $30.0 \pm 0.2$ km\,s$^{-1}$\,kpc$^{-1}$, a bar half-length of $3.2 \pm 0.8$\,kpc. Our model shows large streaming motions around spiral arms, and strong radial motions well beyond the inner bar. This model highlights the complex motions of a dynamic MW-type galaxy and has the potential to offer valuable insight into how our Galaxy regulates the ISM and star formation.
△ Less
Submitted 11 June, 2024; v1 submitted 15 May, 2024;
originally announced May 2024.
-
Assessing the accuracy of the star formation rate measurements by direct star count in molecular clouds
Authors:
Sami Dib,
Jian Wen Zhou,
Sébastien Comerón,
Luis E. Garduño,
Valery V. Kravtsov,
Paul C. Clark,
Guang-Xing Li,
Maritza A. Lara-López,
Tie Liu,
Mohsen Shadmehri,
James R. Doughty
Abstract:
Star formation estimates based on the counting of YSOs is commonly applied to nearby star-forming regions in the Galaxy. With this method, the SFRs are measured using the counts of YSOs in a particular protostellar Class, a typical protostellar mass, and the lifetime associated with this Class. However, the assumptions underlying the validity of the method such as that of a constant star formation…
▽ More
Star formation estimates based on the counting of YSOs is commonly applied to nearby star-forming regions in the Galaxy. With this method, the SFRs are measured using the counts of YSOs in a particular protostellar Class, a typical protostellar mass, and the lifetime associated with this Class. However, the assumptions underlying the validity of the method such as that of a constant star formation history (SFH) and whether the method is valid for all protostellar Classes has never been fully tested. In this work, we use Monte Carlo models to test the validity of the method. We build synthetic clusters in which stars form at times that are randomly drawn from a specified SFH. The latter is either constant or time-dependent with a burst like behavior. The masses of the protostars are randomly drawn from an IMF which can be either similar to that of the Milky Way field or be variable . For each star in every cluster, the lifetimes associated with the different protostellar classes are also randomly drawn from Gaussian distribution functions centered around their most likely value as suggested by the observations. We find that only the SFR derived using the Class 0 population can reproduce the true SFR at all epochs, and this is true irrespective of the shape of the SFH. For a constant SFH, the SFR derived using the more evolved populations of protostars (Classes I, F, II, and III) reproduce the real SFR only at later epochs which correspond to epochs at which their numbers have reached a steady state. For a time-dependent burst-like SFH, all SFR estimates based on the number counts of the evolved populations fail to reproduce the true SFR. We also show how the offsets between Class I and Class II based SFRs and the true SFR plotted as a function of the number ratios of Class I and Class II versus Class III YSOs can be used in order to constrain the SFH of observed molecular clouds.
△ Less
Submitted 30 April, 2024;
originally announced May 2024.
-
N$_2$H$^+$(1-0) as a tracer of dense gas in and between spiral arms
Authors:
O. Feher,
S. E. Ragan,
F. D. Priestley,
P. C. Clark,
T. J. T. Moore
Abstract:
Recent advances in identifying giant molecular filaments in galactic surveys allow us to study the interstellar material and its dense, potentially star forming phase on scales comparable to resolved extragalactic clouds. Two large filaments detected in the CHIMPS $^{13}$CO(3-2) survey, one in the Sagittarius-arm and one in an inter-arm region, were mapped with dense gas tracers inside a 0.06 deg…
▽ More
Recent advances in identifying giant molecular filaments in galactic surveys allow us to study the interstellar material and its dense, potentially star forming phase on scales comparable to resolved extragalactic clouds. Two large filaments detected in the CHIMPS $^{13}$CO(3-2) survey, one in the Sagittarius-arm and one in an inter-arm region, were mapped with dense gas tracers inside a 0.06 deg$^2$ area and with a spatial resolution of around 0.4 and 0.65 pc at the distance of the targets using the IRAM 30m telescope, to investigate the environmental dependence of the dense gas fraction. The N$_2$H$^+$(1-0) transition, an excellent tracer of the dense gas, was detected in parsec-scale, elliptical clumps and with a filling factor of around 8.5% in our maps. The N$_2$H$^+$-emitting areas appear to have higher dense gas fraction (e.g. the ratio of N$_2$H$^+$ and $^{13}$CO emission) in the inter-arm than in the arm which is opposite to the behaviour found by previous studies, using dust emission rather than N$_2$H$^+$ as a tracer of dense gas. However, the arm filament is brighter in $^{13}$CO and the infrared emission of dust, and the dense gas fraction determined as above is governed by the $^{13}$CO brightness. We caution that measurements regarding the distribution and fraction of dense gas on these scales may be influenced by many scale- and environment-dependent factors, as well as the chemistry and excitation of the particular tracers, then consider several scenarios that can reproduce the observed effect.
△ Less
Submitted 28 March, 2024;
originally announced March 2024.
-
Heavy Black Hole Seed Formation in High-z Atomic Cooling Halos
Authors:
Lewis R. Prole,
John A. Regan,
Simon C. O. Glover,
Ralf S. Klessen,
Felix D. Priestley,
Paul C. Clark
Abstract:
Halos with masses in excess of the atomic limit are believed to be ideal environments in which to form heavy black hole seeds with masses above 10^3 Msun. In cases where the H_2 fraction is suppressed this is expected to lead to reduced fragmentation of the gas and the generation of a top heavy initial mass function. In extreme cases this can result in the formation of massive black hole seeds. Re…
▽ More
Halos with masses in excess of the atomic limit are believed to be ideal environments in which to form heavy black hole seeds with masses above 10^3 Msun. In cases where the H_2 fraction is suppressed this is expected to lead to reduced fragmentation of the gas and the generation of a top heavy initial mass function. In extreme cases this can result in the formation of massive black hole seeds. Resolving the initial fragmentation scale and the resulting protostellar masses has, until now, not been robustly tested. Cosmological simulations were performed with the moving mesh code Arepo using a primordial chemistry network until z = 11. Three haloes with masses in excess of the atomic cooling mass were then selected for detailed examination via zoom-ins. The highest resolution simulations resolve densities up to 10^-6 g cm^-3 (10^18 cm^-3) and capture a further 100 yr of fragmentation behaviour at the center of the halo. Our simulations show intense fragmentation in the central region of the halos, leading to a large number of near-solar mass protostars. Despite the increased fragmentation the halos produce a protostellar mass spectrum that peaks at higher masses relative to standard Population III star forming halos. The most massive protostars have accretion rates of 10^-3-10^-1 Msun yr^-1 after the first 100 years of evolution, while the total mass of the central region grows at 1 Msun yr^-1. Lower resolution zoom-ins show that the total mass of the system continues to accrete at 1 Msun yr^-1 for at least 10^4 yr, although how this mass is distributed amongst the rapidly growing number of protostars is unclear. However, assuming that a fraction of stars can continue to accrete rapidly the formation of a sub-population of stars with masses in excess of 10^3 Msun is likely in these halos.
△ Less
Submitted 11 December, 2023;
originally announced December 2023.
-
Population III star formation: multiple gas phases prevent the use of an equation of state at high densities
Authors:
Lewis R. Prole,
Paul C. Clark,
Felix D. Priestley,
Simon C. O. Glover,
John A. Regan
Abstract:
Advanced primordial chemistry networks have been developed to model the collapse of metal-free baryonic gas within the gravitational well of dark matter (DM) halos and its subsequent collapse into Population III stars. At the low densities of 10^-26-10^-21 g cm-3 (10-3-10^2 cm-3) the collapse is dependent on H2 production, which is a function of the compressional heating provided by the DM potenti…
▽ More
Advanced primordial chemistry networks have been developed to model the collapse of metal-free baryonic gas within the gravitational well of dark matter (DM) halos and its subsequent collapse into Population III stars. At the low densities of 10^-26-10^-21 g cm-3 (10-3-10^2 cm-3) the collapse is dependent on H2 production, which is a function of the compressional heating provided by the DM potential. Once the gas decouples from the DM, the temperature-density relationship follows a well established path dictated by various chemical reactions until the formation of the protostar at 10^-4 g cm-3 (10^19 cm-3). Here we explore the feasibility of replacing the chemical network (CN) with a barotropic equation of state (EoS) just before the formation of the first protostar, to reduce the computational load of simulating the further fragmentation, evolution and characteristics of the very high density gas. We find a significant reduction in fragmentation when using the EoS. The EoS method produces a protostellar mass distribution that peaks at higher masses when compared to CN runs. The change in fragmentation behaviour is due to a lack of cold gas falling in through the disc around the first protostar when using an EoS. Despite this, the total mass accreted across all sinks was invariant to the switch to an EoS, hence the star formation rate (Msun yr^-1) is accurately predicted using an EoS. The EoS routine is approximately 4000 times faster than the CN, however this numerical gain is offset by the lack of accuracy in modelling secondary protostar formation and hence its use must be considered carefully.
△ Less
Submitted 19 January, 2024; v1 submitted 16 October, 2023;
originally announced October 2023.
-
NEATH II: N$_2$H$^+$ as a tracer of imminent star formation in quiescent high-density gas
Authors:
F. D. Priestley,
P. C. Clark,
S. C. O. Glover,
S. E. Ragan,
O. Fehér,
L. R. Prole,
R. S. Klessen
Abstract:
Star formation activity in molecular clouds is often found to be correlated with the amount of material above a column density threshold of $\sim 10^{22} \, {\rm cm^{-2}}$. Attempts to connect this column density threshold to a ${\it volume}$ density above which star formation can occur are limited by the fact that the volume density of gas is difficult to reliably measure from observations. We po…
▽ More
Star formation activity in molecular clouds is often found to be correlated with the amount of material above a column density threshold of $\sim 10^{22} \, {\rm cm^{-2}}$. Attempts to connect this column density threshold to a ${\it volume}$ density above which star formation can occur are limited by the fact that the volume density of gas is difficult to reliably measure from observations. We post-process hydrodynamical simulations of molecular clouds with a time-dependent chemical network, and investigate the connection between commonly-observed molecular species and star formation activity. We find that many molecules widely assumed to specifically trace the dense, star-forming component of molecular clouds (e.g. HCN, HCO$^+$, CS) actually also exist in substantial quantities in material only transiently enhanced in density, which will eventually return to a more diffuse state without forming any stars. By contrast, N$_2$H$^+$ only exists in detectable quantities above a volume density of $10^4 \, {\rm cm^{-3}}$, the point at which CO, which reacts destructively with N$_2$H$^+$, begins to deplete out of the gas phase onto grain surfaces. This density threshold for detectable quantities of N$_2$H$^+$ corresponds very closely to the volume density at which gas becomes irreversibly gravitationally bound in the simulations: the material traced by N$_2$H$^+$ never reverts to lower densities, and quiescent regions of molecular clouds with visible N$_2$H$^+$ emission are destined to eventually form stars. The N$_2$H$^+$ line intensity is likely to directly correlate with the star formation rate averaged over timescales of around a Myr.
△ Less
Submitted 9 October, 2023;
originally announced October 2023.
-
Non-Equilibrium Abundances Treated Holistically (NEATH): the molecular composition of star-forming clouds
Authors:
F. D. Priestley,
P. C. Clark,
S. C. O. Glover,
S. E. Ragan,
O. Fehér,
L. R. Prole,
R. S. Klessen
Abstract:
Much of what we know about molecular clouds, and by extension star formation, comes from molecular line observations. Interpreting these correctly requires knowledge of the underlying molecular abundances. Simulations of molecular clouds typically only model species that are important for the gas thermodynamics, which tend to be poor tracers of the denser material where stars form. We construct a…
▽ More
Much of what we know about molecular clouds, and by extension star formation, comes from molecular line observations. Interpreting these correctly requires knowledge of the underlying molecular abundances. Simulations of molecular clouds typically only model species that are important for the gas thermodynamics, which tend to be poor tracers of the denser material where stars form. We construct a framework for post-processing these simulations with a full time-dependent chemical network, allowing us to model the behaviour of observationally-important species not present in the reduced network used for the thermodynamics. We use this to investigate the chemical evolution of molecular gas under realistic physical conditions. We find that molecules can be divided into those which reach peak abundances at moderate densities ($10^3 \, {\rm cm^{-3}}$) and decline sharply thereafter (such as CO and HCN), and those which peak at higher densities and then remain roughly constant (e.g. NH$_3$, N$_2$H$^+$). Evolving the chemistry with physical properties held constant at their final values results in a significant overestimation of gas-phase abundances for all molecules, and does not capture the drastic variations in abundance caused by different evolutionary histories. The dynamical evolution of molecular gas cannot be neglected when modelling its chemistry.
△ Less
Submitted 24 July, 2023;
originally announced July 2023.
-
Do simulated molecular clouds look like real ones?
Authors:
F. D. Priestley,
P. C. Clark,
A. P Whitworth
Abstract:
Simulations of molecular clouds often begin from highly idealised initial conditions, such as a uniform-density sphere with an artificially imposed turbulent velocity field. While the resulting structures may appear qualitatively similar to those detected in continuum and line observations, it is unclear whether they are genuinely representative of real molecular clouds. Recent observational work…
▽ More
Simulations of molecular clouds often begin from highly idealised initial conditions, such as a uniform-density sphere with an artificially imposed turbulent velocity field. While the resulting structures may appear qualitatively similar to those detected in continuum and line observations, it is unclear whether they are genuinely representative of real molecular clouds. Recent observational work has discovered a tight, often close-to-linear relationship between the integrated intensity of molecular lines and the total column density of the cloud material. We combine magnetohydrodynamical simulations, time-dependent chemistry, and radiative transfer to produce synthetic molecular line observations of model clouds. We find similarly tight correlations between line intensity and column density to those observed, although the linear behaviour is only seen in isolated (as opposed to colliding) model clouds. This linear relationship is not due to optically thin emission; all lines investigated have high optical depths, and the increase in integrated intensity with column density is due to higher velocity dispersion along the line of sight. Overall, the idealised models commonly used in the literature appear to be reasonably accurate representations of real molecular clouds.
△ Less
Submitted 12 January, 2023;
originally announced January 2023.
-
From dark matter halos to pre-stellar cores: High resolution follow-up of cosmological Lyman-Werner simulations
Authors:
Lewis R. Prole,
Anna T. P. Schauer,
Paul C. Clark,
Simon C. O. Glover,
Felix D. Priestley,
Ralf S. Klessen
Abstract:
Molecular hydrogen allows cooling in primordial gas, facilitating its collapse into Population III stars within primordial halos. Lyman-Werner (LW) radiation from these stars can escape the halo and delay further star formation by destroying H$_2$ in other halos. As cosmological simulations show that increasing the background LW field strength increases the average halo mass required for star form…
▽ More
Molecular hydrogen allows cooling in primordial gas, facilitating its collapse into Population III stars within primordial halos. Lyman-Werner (LW) radiation from these stars can escape the halo and delay further star formation by destroying H$_2$ in other halos. As cosmological simulations show that increasing the background LW field strength increases the average halo mass required for star formation, we perform follow-up simulations of selected halos to investigate the knock-on effects this has on the Population III IMF. We follow 5 halos for each of the $J_{21}$ = 0, 0.01 and 0.1 LW field strengths, resolving the pre-stellar core density of $10^{-6}$ g cm$^{-3}$ (10$^{18}$ cm$^{-3}$) before inserting sink particles and following the fragmentation behaviour for hundreds of years further. We find that the mass accreted onto sinks by the end of the simulations is proportional to the mass within the $\sim 10^{-2}$ pc molecular core, which is not correlated to the initial mass of the halo. As such, the IMFs for masses above the brown dwarf limit show little dependence on the LW strength, although they do show variance in the number of low-mass clumps formed. As the range of background LW field strengths tested here covers the most likely values from literature, we conclude that the IMF for so-called Pop III.2 stars is not significantly different from the initial population of Pop III.1 stars. The primordial IMF therefore likely remains unchanged until the formation of the next generation of Population II stars.
△ Less
Submitted 19 January, 2023; v1 submitted 2 January, 2023;
originally announced January 2023.
-
Fragmentation induced starvation in Population III star formation: a resolution study
Authors:
Lewis R. Prole,
Paul C. Clark,
Ralf S. Klessen,
Simon C. O. Glover
Abstract:
The Population III initial mass function (IMF) is currently unknown, but recent studies agree that fragmentation of primordial gas gives a broader IMF than the initially suggested singular star per halo. In this study we introduce sink particle mergers into Arepo, to perform the first resolution study for primordial star formation simulations and present the first Population III simulations to run…
▽ More
The Population III initial mass function (IMF) is currently unknown, but recent studies agree that fragmentation of primordial gas gives a broader IMF than the initially suggested singular star per halo. In this study we introduce sink particle mergers into Arepo, to perform the first resolution study for primordial star formation simulations and present the first Population III simulations to run up to densities of 10-6g cm-3 for hundreds of years after the formation of sink particles. The total number of sinks formed increases with increasing sink particle creation density, without achieving numerical convergence. The total mass in sinks remains invariant to the maximum resolution and is safe to estimate using low resolution studies. This results in an IMF that shifts towards lower masses with increasing resolution. Greater numbers of sinks cause increased fragmentation-induced starvation of the most massive sink, yielding lower accretion rates, masses and ionising photons emitted per second. The lack of convergence up to densities 2 orders of magnitudes higher than all relevant chemical reactions suggests that the number of sinks will continue to grow with increasing resolution until H2 is fully dissociated and the collapse becomes almost adiabatic at 10-4g cm-3. These results imply that many Population III studies utilising sink particles have produced IMFs which have overestimated the masses of primordial stars, and underestimated the number of stars formed. In the highest resolution runs, sinks with masses capable of surviving until the present day had an ejection fraction of 0.21.
△ Less
Submitted 20 December, 2021;
originally announced December 2021.
-
On the density regime probed by HCN emission
Authors:
Gerwyn H. Jones,
Paul C. Clark,
Simon C. O. Glover,
Alvaro Hacar
Abstract:
HCN J$\, =\,$1$\, -\,$0 emission is commonly used as a dense gas tracer, thought to mainly arise from gas with densities $\mathrm{\sim 10^4\ -\ 10^5\ cm^{-3}}$. This has made it a popular tracer in star formation studies. However, there is increasing evidence from observational surveys of `resolved' molecular clouds that HCN can trace more diffuse gas. We investigate the relationship between gas d…
▽ More
HCN J$\, =\,$1$\, -\,$0 emission is commonly used as a dense gas tracer, thought to mainly arise from gas with densities $\mathrm{\sim 10^4\ -\ 10^5\ cm^{-3}}$. This has made it a popular tracer in star formation studies. However, there is increasing evidence from observational surveys of `resolved' molecular clouds that HCN can trace more diffuse gas. We investigate the relationship between gas density and HCN emission through post-processing of high resolution magnetohydrodynamical simulations of cloud-cloud collisions. We find that HCN emission traces gas with a mean volumetric density of $\mathrm{\sim 3 \times 10^3\ cm^{-3}}$ and a median visual extinction of $\mathrm{\sim 5\ mag}$. We therefore predict a characteristic density that is an order of magnitude less than the "standard" characteristic density of $\mathrm{n \sim 3 \times 10^4\ cm^{-3}}$. Indeed, we find in some cases that there is clear HCN emission from the cloud even though there is no gas denser than this standard critical density. We derive luminosity-to-mass conversion factors for the amount of gas at $A_{\rm V} > 8$ or at densities $n > 2.85 \times 10^{3} \: {\rm cm^{-3}}$ or $n > 3 \times 10^{4} \: {\rm cm^{-3}}$, finding values of $α_{\rm HCN} = 6.79, 8.62$ and $27.98 \: {\rm M_{\odot}} ({\rm K \, km \, s^{-1} \, pc^{2}})$, respectively. In some cases, the luminosity to mass conversion factor predicted mass in regions where in actuality there contains no mass.
△ Less
Submitted 18 January, 2023; v1 submitted 10 December, 2021;
originally announced December 2021.
-
Towards the impact of GMC collisions on the star formation rate
Authors:
Glen H. Hunter,
Paul C. Clark,
Simon C. O. Glover,
Ralf S. Klessen
Abstract:
Collisions between giant molecular clouds (GMCs) are one of the pathways for massive star formation, due to the high densities created. However the enhancement of the star formation rate (SFR) is not well constrained. In this study we perform a parameter study of cloud-cloud collisions, and investigate how the resulting SFR depends on the details of set-up. Our parameter study explores variations…
▽ More
Collisions between giant molecular clouds (GMCs) are one of the pathways for massive star formation, due to the high densities created. However the enhancement of the star formation rate (SFR) is not well constrained. In this study we perform a parameter study of cloud-cloud collisions, and investigate how the resulting SFR depends on the details of set-up. Our parameter study explores variations in: collision speed; magnetic field inclination (with respect to the collisional axis); and resolution, as defined by the number of cells per Jeans length. In all our collision simulations we find a factor of 2-3 increase in the SFR compared to our no collision simulation, with star formation beginning sooner with a) high collisional velocities, b) parallel orientation between the magnetic field and collision axis, c) and lower resolution. The mean virial parameter of high density (and thus possible star-forming) gas increases with collisional velocity, but has little variation with magnetic field inclination. The alignment of the velocity and magnetic field remains uniform in low density environments but becomes more perpendicular with increasing density, indicating the compression of the magnetic field by collapsing gas. Comparing the trends in the SFR with other GMC collision studies, we find good agreement with studies that account for the gravitational boundedness of the gas in their star formation algorithm, but not with those that simply form stars above a prescribed density threshold. This suggests that the latter approach should be used with caution when modelling star formation on resolved cloud scales.
△ Less
Submitted 11 January, 2023; v1 submitted 13 September, 2021;
originally announced September 2021.
-
Extreme adaptive optics astrometry of R136. Searching for high proper motion stars
Authors:
Zeinab Khorrami,
M. Langlois,
F. Vakili,
P. C. Clark,
A. S. M. Buckner,
M. Gonzalez,
P. Crowther,
R. Wunsch,
J. Palous,
A. Boccaletti,
S. Lumsden,
E. Moraux
Abstract:
We compared high-contrast near-infrared images of the core of R136 taken by VLT/SPHERE, in two epochs separated by 3.06 years. For the first time we monitored the dynamics of the detected sources in the core of R136 from a ground-based telescope with adaptive optics. The aim of these observations was to search for High prOper Motion cAndidates (HOMAs) in the central region of R136 (r<6") where it…
▽ More
We compared high-contrast near-infrared images of the core of R136 taken by VLT/SPHERE, in two epochs separated by 3.06 years. For the first time we monitored the dynamics of the detected sources in the core of R136 from a ground-based telescope with adaptive optics. The aim of these observations was to search for High prOper Motion cAndidates (HOMAs) in the central region of R136 (r<6") where it has been challenging for other instruments. Two bright sources (K<15mag and V<16mag) are located near R136a1 and R136c (massive WR stars) and have been identified as potential HOMAs. These sources have significantly shifted in the images with respect to the mean shift of all reliable detected sources and their neighbours, and six times their own astrometric errors. We calculate their proper motions to be 1.36\pm0.22 mas/yr (321\pm52 km/s) and 1.15\pm0.11 mas/yr (273\pm26 km/s). We discuss different possible scenarios to explain the magnitude of such extreme proper motions, and argue for the necessity to conduct future observations to conclude on the nature of HOMAs in the core of R136.
△ Less
Submitted 26 April, 2021;
originally announced April 2021.
-
High contrast and resolution near infrared photometry of the core of R136
Authors:
Zeinab Khorrami,
Maud Langlois,
Paul C. Clark,
Farrokh Vakili,
Anne S. M. Buckner,
Marta Gonzalez,
Paul Crowther,
Richard Wunsch,
Jan Palous,
Stuart Lumsden,
Estelle Moraux
Abstract:
We present the sharpest and deepest near infrared photometric analysis of the core of R136, a newly formed massive star cluster at the centre of the 30 Doradus star forming region in the Large Magellanic Cloud. We used the extreme adaptive optics of the SPHERE focal instrument implemented on the ESO Very Large Telescope and operated in its IRDIS imaging mode, for the second time with longer exposu…
▽ More
We present the sharpest and deepest near infrared photometric analysis of the core of R136, a newly formed massive star cluster at the centre of the 30 Doradus star forming region in the Large Magellanic Cloud. We used the extreme adaptive optics of the SPHERE focal instrument implemented on the ESO Very Large Telescope and operated in its IRDIS imaging mode, for the second time with longer exposure time in the H- and K filters. Our aim was to (i) increase the number of resolved sources in the core of R136, and (ii) to compare with the first epoch to classify the properties of the detected common sources between the two epochs. Within the field of view (FOV) of 10.8"x12.1" (2.7pc x3.0pc), we detected 1499 sources in both H and K filters, for which 76% of these sources have visual companions closer than 0.2". The larger number of detected sources, enabled us to better sample the mass function (MF). The MF slopes are estimated at ages of 1, 1.5 and 2 Myr, at different radii, and for different mass ranges. The MF slopes for the mass range of 10-300 solar-mass are about 0.3 dex steeper than the mass range of 3-300 solar-mass, for the whole FOV and different radii. Comparing the JHK colours of 790 sources common in between the two epochs, 67% of detected sources in the outer region (r >3") are not consistent with evolutionary models at 1-2 Myr and with extinctions similar to the average cluster value, suggesting an origin from ongoing star formation within 30 Doradus, unrelated to R136.
△ Less
Submitted 11 February, 2021;
originally announced February 2021.
-
The Cloud Factory II: Gravoturbulent Kinematics of Resolved Molecular Clouds in a Galactic Potential
Authors:
Andres F. Izquierdo,
Rowan J. Smith,
Simon C. O. Glover,
Ralf S. Klessen,
Robin G. Tress,
Mattia C. Sormani,
Paul C. Clark,
Ana Duarte-Cabral,
Catherine Zucker
Abstract:
We present a statistical analysis of the gravoturbulent velocity fluctuations in molecular cloud complexes extracted from our "Cloud Factory" galactic-scale ISM simulation suite. For this purpose, we produce non-LTE $^{12}$CO J=1-0 synthetic observations and apply the Principal Component Analysis (PCA) reduction technique on a representative sample of cloud complexes. The velocity fluctuations are…
▽ More
We present a statistical analysis of the gravoturbulent velocity fluctuations in molecular cloud complexes extracted from our "Cloud Factory" galactic-scale ISM simulation suite. For this purpose, we produce non-LTE $^{12}$CO J=1-0 synthetic observations and apply the Principal Component Analysis (PCA) reduction technique on a representative sample of cloud complexes. The velocity fluctuations are self-consistently generated by different physical mechanisms at play in our simulations, which include galactic-scale forces, gas self-gravity, and supernova feedback. The statistical analysis suggests that, even though purely gravitational effects are necessary to reproduce standard observational laws, they are not sufficient in most cases. We show that the extra injection of energy from supernova explosions plays a key role in establishing the global turbulent field and the local dynamics and morphology of molecular clouds. Additionally, we characterise structure function scaling parameters as a result of cloud environmental conditions: some of the complexes are immersed in diffuse (inter-arm) or dense (spiral-arm) environments, and others are influenced by embedded or external supernovae. In quiescent regions, we obtain time-evolving trajectories of scaling parameters driven by gravitational collapse and supersonic turbulent flows. Our findings suggests that a PCA-based statistical study is a robust method to diagnose the physical mechanisms that drive the gravoturbulent properties of molecular clouds. Also, we present a new open source module, the PCAFACTORY, which smartly performs PCA to extract velocity structure functions from simulated or real data of the ISM in a user-friendly way. Software DOI: 10.5281/zenodo.3822718
△ Less
Submitted 4 November, 2020;
originally announced November 2020.
-
CHIMPS2: Survey description and $^{12}$CO emission in the Galactic Centre
Authors:
D. J. Eden,
T. J. T. Moore,
M. J. Currie,
A. J. Rigby,
E. Rosolowsky,
Y. Su,
Kee-Tae Kim,
H. Parsons,
O. Morata,
H. -R. Chen,
T. Minamidani,
Geumsook Park,
S. E. Ragan,
J. S. Urquhart,
R. Rani,
K. Tahani,
S. J. Billington,
S. Deb,
C. Figura,
T. Fujiyoshi,
G. Joncas,
L. W. Liao,
T. Liu,
H. Ma,
P. Tuan-Anh
, et al. (81 additional authors not shown)
Abstract:
The latest generation of Galactic-plane surveys is enhancing our ability to study the effects of galactic environment upon the process of star formation. We present the first data from CO Heterodyne Inner Milky Way Plane Survey 2 (CHIMPS2). CHIMPS2 is a survey that will observe the Inner Galaxy, the Central Molecular Zone (CMZ), and a section of the Outer Galaxy in $^{12}$CO, $^{13}$CO, and C…
▽ More
The latest generation of Galactic-plane surveys is enhancing our ability to study the effects of galactic environment upon the process of star formation. We present the first data from CO Heterodyne Inner Milky Way Plane Survey 2 (CHIMPS2). CHIMPS2 is a survey that will observe the Inner Galaxy, the Central Molecular Zone (CMZ), and a section of the Outer Galaxy in $^{12}$CO, $^{13}$CO, and C$^{18}$O $(J = 3\rightarrow2)$ emission with the Heterodyne Array Receiver Program on the James Clerk Maxwell Telescope (JCMT). The first CHIMPS2 data presented here are a first look towards the CMZ in $^{12}$CO J = 3$\rightarrow$2 and cover $-3^{\circ}\leq\,\ell\,\leq\,5^{\circ}$ and $\mid$b$\mid \leq 0.5^{\circ}$ with angular resolution of 15 arcsec, velocity resolution of 1 km s$^{-1}$, and rms $ΔT_A ^\ast =$ 0.58 K at these resolutions. Such high-resolution observations of the CMZ will be a valuable data set for future studies, whilst complementing the existing Galactic Plane surveys, such as SEDIGISM, the Herschel infrared Galactic Plane Survey, and ATLASGAL. In this paper, we discuss the survey plan, the current observations and data, as well as presenting position-position maps of the region. The position-velocity maps detect foreground spiral arms in both absorption and emission.
△ Less
Submitted 10 September, 2020;
originally announced September 2020.
-
On the emergent System Mass Function: the contest between accretion and fragmentation
Authors:
Paul C. Clark,
Anthony P. Whitworth
Abstract:
We propose a new model for the evolution of a star cluster's System Mass Function (SMF). The model involves both turbulent fragmentation and competitive accretion. Turbulent fragmentation creates low-mass seed proto-systems (i.e. single and multiple protostars). Some of these low-mass seed proto-systems then grow by competitive accretion to produce the high-mass power-law tail of the SMF. Turbulen…
▽ More
We propose a new model for the evolution of a star cluster's System Mass Function (SMF). The model involves both turbulent fragmentation and competitive accretion. Turbulent fragmentation creates low-mass seed proto-systems (i.e. single and multiple protostars). Some of these low-mass seed proto-systems then grow by competitive accretion to produce the high-mass power-law tail of the SMF. Turbulent fragmentation is relatively inefficient, in the sense that the creation of low-mass seed proto-systems only consumes a fraction, $\sim 23\%$ (at most $\sim 50\%$), of the mass available for star formation. The remaining mass is consumed by competitive accretion. Provided the accretion rate onto a proto-system is approximately proportional to its mass ($dm/dt \propto m$), the SMF develops a power-law tail at high masses with the Salpeter slope ($\sim -2.3$). If the rate of supply of mass accelerates, the rate of proto-system formation also accelerates, as appears to be observed in many clusters. However, even if the rate of supply of mass decreases, or ceases and then resumes, the SMF evolves homologously, retaining the same overall shape, and the high-mass power-law tail simply extends to ever higher masses until the supply of gas runs out completely. The Chabrier SMF can be reproduced very accurately if the seed proto-systems have an approximately log-normal mass distribution with median mass $\sim 0.11 {\rm M}_{_\odot}$ and logarithmic standard deviation $σ_{\log_{10}(M/{\rm M}_\odot)}\sim 0.47$).
△ Less
Submitted 22 August, 2020;
originally announced August 2020.
-
Characteristic scale of star formation. I. Clump formation efficiency on local scales
Authors:
D. J. Eden,
T. J. T. Moore,
R. Plume,
A. J. Rigby,
J. S. Urquhart,
K. A. Marsh,
C. H. Peñaloza,
P. C. Clark,
M. W. L. Smith,
K. Tahani,
S. E. Ragan,
M. A. Thompson,
D. Johnstone,
H. Parsons,
R. Rani
Abstract:
We have used the ratio of column densities (CDR) derived independently from the 850-$μ$m continuum JCMT Plane Survey (JPS) and the $^{13}$CO/C$^{18}$O $(J=3-2)$ Heterodyne Inner Milky Way Plane Survey (CHIMPS) to produce maps of the dense-gas mass fraction (DGMF) in two slices of the Galactic Plane centred at $\ell$=30$^{\circ}$ and $\ell$=40$^{\circ}$. The observed DGMF is a metric for the instan…
▽ More
We have used the ratio of column densities (CDR) derived independently from the 850-$μ$m continuum JCMT Plane Survey (JPS) and the $^{13}$CO/C$^{18}$O $(J=3-2)$ Heterodyne Inner Milky Way Plane Survey (CHIMPS) to produce maps of the dense-gas mass fraction (DGMF) in two slices of the Galactic Plane centred at $\ell$=30$^{\circ}$ and $\ell$=40$^{\circ}$. The observed DGMF is a metric for the instantaneous clump-formation efficiency (CFE) in the molecular gas. We split the two fields into velocity components corresponding to the spiral arms that cross them, and a two-dimensional power-spectrum analysis of the spiral arm DGMF maps reveals a break in slope at the approximate size scale of molecular clouds. We interpret this as the characteristic scale of the amplitude of variations in the CFE and a constraint on the dominant mechanism regulating the CFE and, hence, the star-formation efficiency in CO-traced clouds.
△ Less
Submitted 12 October, 2020; v1 submitted 30 June, 2020;
originally announced July 2020.
-
Simulations of the Milky Way's central molecular zone -- II. Star formation
Authors:
Mattia C. Sormani,
Robin G. Tress,
Simon C. O. Glover,
Ralf S. Klessen,
Cara D. Battersby,
Paul C. Clark,
H Perry Hatchfield,
Rowan J. Smith
Abstract:
The Milky Way's central molecular zone (CMZ) has emerged in recent years as a unique laboratory for the study of star formation. Here we use the simulations presented in Tress et al. 2020 to investigate star formation in the CMZ. These simulations resolve the structure of the interstellar medium at sub-parsec resolution while also including the large-scale flow in which the CMZ is embedded. Our ma…
▽ More
The Milky Way's central molecular zone (CMZ) has emerged in recent years as a unique laboratory for the study of star formation. Here we use the simulations presented in Tress et al. 2020 to investigate star formation in the CMZ. These simulations resolve the structure of the interstellar medium at sub-parsec resolution while also including the large-scale flow in which the CMZ is embedded. Our main findings are as follows. (1) While most of the star formation happens in the CMZ ring at $R\gtrsim100 {\, \rm pc}$, a significant amount also occurs closer to SgrA* at $R \lesssim 10{\, \rm pc}$. (2) Most of the star formation in the CMZ happens downstream of the apocentres, consistent with the "pearls-on-a-string" scenario, and in contrast to the notion that an absolute evolutionary timeline of star formation is triggered by pericentre passage. (3) Within the timescale of our simulations ($\sim100$ Myr), the depletion time of the CMZ is constant within a factor of $\sim2$. This suggests that variations in the star formation rate are primarily driven by variations in the mass of the CMZ, caused for example by AGN feedback or externally-induced changes in the bar-driven inflow rate, and not by variations in the depletion time. (4) We study the trajectories of newly born stars in our simulations. We find several examples that have age and 3D velocity compatible with those of the Arches and Quintuplet clusters. Our simulations suggest that these prominent clusters originated near the collision sites where the bar-driven inflow accretes onto the CMZ, at symmetrical locations with respect to the Galactic centre, and that they have already decoupled from the gas in which they were born.
△ Less
Submitted 7 October, 2020; v1 submitted 14 April, 2020;
originally announced April 2020.
-
Simulations of the Milky Way's central molecular zone -- I. Gas dynamics
Authors:
Robin G. Tress,
Mattia C. Sormani,
Simon C. O. Glover,
Ralf S. Klessen,
Cara D. Battersby,
Paul C. Clark,
H Perry Hatchfield,
Rowan J. Smith
Abstract:
We use hydrodynamical simulations to study the Milky Way's central molecular zone (CMZ). The simulations include a non-equilibrium chemical network, the gas self-gravity, star formation and supernova feedback. We resolve the structure of the interstellar medium at sub-parsec resolution while also capturing the interaction between the CMZ and the bar-driven large-scale flow out to $R\sim 5\kpc$. Ou…
▽ More
We use hydrodynamical simulations to study the Milky Way's central molecular zone (CMZ). The simulations include a non-equilibrium chemical network, the gas self-gravity, star formation and supernova feedback. We resolve the structure of the interstellar medium at sub-parsec resolution while also capturing the interaction between the CMZ and the bar-driven large-scale flow out to $R\sim 5\kpc$. Our main findings are as follows: (1) The distinction between inner ($R\lesssim120$~pc) and outer ($120\lesssim R\lesssim450$~pc) CMZ that is sometimes proposed in the literature is unnecessary. Instead, the CMZ is best described as single structure, namely a star-forming ring with outer radius $R\simeq 200$~pc which includes the 1.3$^\circ$ complex and which is directly interacting with the dust lanes that mediate the bar-driven inflow. (2) This accretion can induce a significant tilt of the CMZ out of the plane. A tilted CMZ might provide an alternative explanation to the $\infty$-shaped structure identified in Herschel data by Molinari et al. 2011. (3) The bar in our simulation efficiently drives an inflow from the Galactic disc ($R\simeq 3$~kpc) down to the CMZ ($R\simeq200$~pc) of the order of $1\rm\,M_\odot\,yr^{-1}$, consistent with observational determinations. (4) Supernova feedback can drive an inflow from the CMZ inwards towards the circumnuclear disc of the order of $\sim0.03\,\rm M_\odot\,yr^{-1}$. (5) We give a new interpretation for the 3D placement of the 20 and 50 km s$^{-1}$ clouds, according to which they are close ($R\lesssim30$~pc) to the Galactic centre, but are also connected to the larger-scale streams at $R\gtrsim100$~pc.
△ Less
Submitted 7 October, 2020; v1 submitted 14 April, 2020;
originally announced April 2020.
-
Symmetry-breaking and spin-blockage effects on carrier dynamics in single-layer tungsten diselenide
Authors:
Ro-Ya Liu,
Meng-Kai Lin,
Peng Chen,
Takeshi Suzuki,
Pip C. J. Clark,
Nathan K. Lewis,
Cephise Cacho,
Emma Springate,
Chia-Seng Chang,
Kozo Okazaki,
Wendy Flavell,
Iwao Matsuda,
Tai-Chang Chiang
Abstract:
Understanding carrier creation and evolution in materials initiated by pulsed optical excitation is central to developing ultrafast optoelectronics. We demonstrate herein that the dynamic response of a system can be drastically modified when its physical dimension is reduced to the atomic scale, the ultimate limit of device miniaturization. A comparative study of single-layer (SL) tungsten diselen…
▽ More
Understanding carrier creation and evolution in materials initiated by pulsed optical excitation is central to developing ultrafast optoelectronics. We demonstrate herein that the dynamic response of a system can be drastically modified when its physical dimension is reduced to the atomic scale, the ultimate limit of device miniaturization. A comparative study of single-layer (SL) tungsten diselenide(WSe2) relative to bulk WSe2 shows substantial differences in the transient response as measured by time- and angle-resolved photoemission spectroscopy (TRARPES). The conduction-band minimum in bulk WSe2, populated by optical pumping, decays promptly. The corresponding decay for SL WSe2 is much slower and exhibits two time constants. The results indicate the presence of two distinct decay channels in the SL that are correlated with the breaking of space inversion symmetry in the two-dimensional limit. This symmetry breaking lifts the spin degeneracy of the bands, which in turn causes the blockage of decay for one spin channel. The stark contrast between the single layer and the bulk illustrates the basic carrier scattering processes operating at different timescales that can be substantially modified by dimensional and symmetry-reduction effects.
△ Less
Submitted 22 December, 2019;
originally announced December 2019.
-
Formation sites of Population III star formation: The effects of different levels of rotation and turbulence on the fragmentation behavior of primordial gas
Authors:
Katharina M. J. Wollenberg,
Simon C. O. Glover,
Paul C. Clark,
Ralf S. Klessen
Abstract:
We use the moving-mesh code AREPO to investigate the effects of different levels of rotation and turbulence on the fragmentation of primordial gas and the formation of Population III stars. We consider 9 different combinations of turbulence and rotation and carry out 5 different realizations of each setup, yielding one of the largest sets of simulations of Population III star formation ever perfor…
▽ More
We use the moving-mesh code AREPO to investigate the effects of different levels of rotation and turbulence on the fragmentation of primordial gas and the formation of Population III stars. We consider 9 different combinations of turbulence and rotation and carry out 5 different realizations of each setup, yielding one of the largest sets of simulations of Population III star formation ever performed. We find that fragmentation in Population III star-forming systems is a highly chaotic process and show that the outcomes of individual realizations of the same initial conditions often vary significantly. However, some general trends are apparent. Increasing the turbulent energy promotes fragmentation, while increasing the rotational energy inhibits fragmentation. Within the 1000 yr period that we simulate, runs including turbulence yield flat protostellar mass functions while purely rotational runs show a more top-heavy distribution. The masses of the individual protostars are distributed over a wide range from a few $10^{-3} \, {\rm M_{\odot}}$ to several tens of ${\rm M_\odot}$. The total mass growth rate of the stellar systems remains high throughout the simulations and depends only weakly on the degree of rotation and turbulence. Mergers between protostars are common, but predictions of the merger fraction are highly sensitive to the criterion used to decide whether two protostars should merge. Previous studies of Population III star formation have often considered only one realization per set of initial conditions. However, our results demonstrate that robust trends can only be reliably identified by considering averages over a larger sample of runs.
△ Less
Submitted 29 January, 2020; v1 submitted 13 December, 2019;
originally announced December 2019.
-
The Cloud Factory I: Generating resolved filamentary molecular clouds from galactic-scale forces
Authors:
Rowan J. Smith,
Robin G. Treß,
Mattia C. Sormani,
Simon C. O. Glover,
Ralf S. Klessen,
Paul C. Clark,
Andrés F. Izquierdo,
Ana Duarte Cabral,
Catherine Zucker
Abstract:
We introduce a new suite of simulations, "The Cloud Factory", which self-consistently forms molecular cloud complexes at high enough resolution to resolve internal substructure (up to 0.25 Msol in mass) all while including galactic-scale forces. We use a version of the Arepo code modified to include a detailed treatment of the physics of the cold molecular ISM, and an analytical galactic gravitati…
▽ More
We introduce a new suite of simulations, "The Cloud Factory", which self-consistently forms molecular cloud complexes at high enough resolution to resolve internal substructure (up to 0.25 Msol in mass) all while including galactic-scale forces. We use a version of the Arepo code modified to include a detailed treatment of the physics of the cold molecular ISM, and an analytical galactic gravitational potential for computational efficiency. The simulations have nested levels of resolution, with the lowest layer tied to tracer particles injected into individual cloud complexes. These tracer refinement regions are embedded in the larger simulation so continue to experience forces from outside the cloud. This allows the simulations to act as a laboratory for testing the effect of galactic environment on star formation. Here we introduce our method and investigate the effect of galactic environment on filamentary clouds. We find that cloud complexes formed after a clustered burst of feedback, have shorter lengths and are less likely to fragment compared to quiescent clouds (e.g. the Musca filament) or those dominated by the galactic potential (e.g. Nessie). Spiral arms and differential rotation preferentially align filaments, but strong feedback randomises them. Long filaments formed within the cloud complexes are necessarily coherent with low internal velocity gradients, which has implications for the formation of filamentary star-clusters. Cloud complexes formed in regions dominated by supernova feedback have fewer star-forming cores, and these are more widely distributed. These differences show galactic-scale forces can have a significant impact on star formation within molecular clouds.
△ Less
Submitted 18 November, 2019; v1 submitted 13 November, 2019;
originally announced November 2019.
-
Simulations of the star-forming molecular gas in an interacting M51-like galaxy
Authors:
Robin G. Tress,
Rowan J. Smith,
Mattia C. Sormani,
Simon C. O. Glover,
Ralf S. Klessen,
Mordecai-Mark Mac Low,
Paul C. Clark
Abstract:
We present here the first of a series of papers aimed at better understanding the evolution and properties of giant molecular clouds (GMCs) in a galactic context. We perform high resolution, three-dimensional {\sc arepo} simulations of an interacting galaxy inspired by the well-observed M51 galaxy. Our fiducial simulations include a non-equilibrium, time-dependent, chemical network that follows th…
▽ More
We present here the first of a series of papers aimed at better understanding the evolution and properties of giant molecular clouds (GMCs) in a galactic context. We perform high resolution, three-dimensional {\sc arepo} simulations of an interacting galaxy inspired by the well-observed M51 galaxy. Our fiducial simulations include a non-equilibrium, time-dependent, chemical network that follows the evolution of atomic and molecular hydrogen as well as carbon and oxygen self-consistently. Our calculations also treat gas self-gravity and subsequent star formation (described by sink particles), and coupled supernova feedback. In the densest parts of the simulated interstellar medium (ISM) we reach sub-parsec resolution, granting us the ability to resolve individual GMCs and their formation and destruction self-consistently throughout the galaxy. In this initial work we focus on the general properties of the ISM with a particular focus on the cold star-forming gas. We discuss the role of the interaction with the companion galaxy in generating cold molecular gas and controlling stellar birth. We find that while the interaction drives large-scale gas flows and induces spiral arms in the galaxy, it is of secondary importance in determining gas fractions in the different ISM phases and the overall star-formation rate. The behaviour of the gas on small GMC scales instead is mostly controlled by the self-regulating property of the ISM driven by coupled feedback.
△ Less
Submitted 23 September, 2019;
originally announced September 2019.
-
The geometry of the gas surrounding the Central Molecular Zone: on the origin of localised molecular clouds with extreme velocity dispersions
Authors:
Mattia C. Sormani,
Robin G. Tress,
Simon C. O. Glover,
Ralf S. Klessen,
Ashley T. Barnes,
Cara D. Battersby,
Paul C. Clark,
H Perry Hatchfield,
Rowan J. Smith
Abstract:
Observations of molecular gas near the Galactic centre ($| l | < 10^\circ$, $| b | < 1^\circ$) reveal the presence of a distinct population of enigmatic compact clouds which are characterised by extreme velocity dispersions ($Δv > 100\, \rm km/s$). These Extended Velocity Features (EVFs) are very prominent in the datacubes and dominate the kinematics of molecular gas just outside the Central Molec…
▽ More
Observations of molecular gas near the Galactic centre ($| l | < 10^\circ$, $| b | < 1^\circ$) reveal the presence of a distinct population of enigmatic compact clouds which are characterised by extreme velocity dispersions ($Δv > 100\, \rm km/s$). These Extended Velocity Features (EVFs) are very prominent in the datacubes and dominate the kinematics of molecular gas just outside the Central Molecular Zone (CMZ). The prototypical example of such a cloud is Bania Clump 2. We show that similar features are naturally produced in simulations of gas flow in a realistic barred potential. We analyse the structure of the features obtained in the simulations and use this to interpret the observations. We find that the features arise from collisions between material that has been infalling rapidly along the dust lanes of the Milky Way bar and material that belongs to one of the following two categories: (i) material that has `overshot' after falling down the dust lanes on the opposite side; (ii) material which is part of the CMZ. Both types of collisions involve gas with large differences in the line-of-sight velocities, which is what produces the observed extreme velocity dispersions. Examples of both categories can be identified in the observations. If our interpretation is correct, we are directly witnessing (a) collisions of clouds with relative speeds of $\sim 200\, \rm km/s$ and (b) the process of accretion of fresh gas onto the CMZ.
△ Less
Submitted 24 June, 2019;
originally announced June 2019.
-
A Code to Make Your Own Synthetic ObservaTIonS (MYOSOTIS)
Authors:
Zeinab Khorrami,
Pouria Khalaj,
Anne S. M. Buckner,
Paul C. Clark,
Estelle Moraux,
Stuart Lumsden,
Isabelle Joncour,
Rene D. Oudmaijer,
Ignacio de la Calle,
Jose M. Herrera-Fernandez,
Frederique Motte,
Jose Manuel Blanco,
Luis Valero-Martin
Abstract:
We introduce our new code MYOSOTIS (Make Your Own Synthetic ObservaTIonS) which is designed to produce synthetic observations from simulated clusters. The code can synthesise observations from both ground- and spaced-based observatories, for a range of different filters, observational conditions and angular/spectral resolution. In this paper, we highlight some of the features of MYOSOTIS, creating…
▽ More
We introduce our new code MYOSOTIS (Make Your Own Synthetic ObservaTIonS) which is designed to produce synthetic observations from simulated clusters. The code can synthesise observations from both ground- and spaced-based observatories, for a range of different filters, observational conditions and angular/spectral resolution. In this paper, we highlight some of the features of MYOSOTIS, creating synthetic observations from young massive star clusters. Our model clusters are simulated using nbody6 code and have different total masses, half-mass radii, and binary fractions. The synthetic observations are made at the age of 2 Myr with Solar metallicity and under different extinction conditions. For each cluster, we create synthetic images of the Hubble Space Telescope (HST) in the visible (WFPC2/F555W) as well as Very Large Telescopes (VLT) in the nearIR (SPHERE/IRDIS/Ks). We show how MYOSOTIS can be used to look at mass function (MF) determinations. For this aim we re-estimate stellar masses using a photometric analysis on the synthetic images. The synthetic MF slopes are compared to their actual values. Our photometric analysis demonstrate that depending on the adopted filter, extinction, angular resolution and pixel sampling of the instruments, the power-law index of the underlying MFs can be shallower than the observed ones by at least 0.25 dex which is in agreement with the observed discrepancies reported in the literature, specially for young star clusters.
△ Less
Submitted 14 February, 2019;
originally announced February 2019.
-
Synthetic [CII] emission maps of a simulated molecular cloud in formation
Authors:
A. Franeck,
S. Walch,
D. Seifried,
S. D. Clarke,
V. Ossenkopf-Okada,
S. C. O. Glover,
R. S. Klessen,
P. Girichidis,
T. Naab,
R. Wünsch,
P. C. Clark,
E. Pellegrini,
T. Peters
Abstract:
The C$^{+}$ ion is an important coolant of interstellar gas, and so the [CII] fine structure line is frequently observed in the interstellar medium. However, the physical and chemical properties of the [CII]-emitting gas are still unclear. We carry out non-LTE radiative transfer simulations with RADMC-3D to study the [CII] line emission from a young, turbulent molecular cloud before the onset of s…
▽ More
The C$^{+}$ ion is an important coolant of interstellar gas, and so the [CII] fine structure line is frequently observed in the interstellar medium. However, the physical and chemical properties of the [CII]-emitting gas are still unclear. We carry out non-LTE radiative transfer simulations with RADMC-3D to study the [CII] line emission from a young, turbulent molecular cloud before the onset of star formation, using data from the SILCC-Zoom project. The [CII] emission is optically thick over 40% of the observable area with $I_{[\textrm{CII}]} > 0.5$ K km s$^{-1}$. To determine the physical properties of the [CII] emitting gas, we treat the [CII] emission as optically thin. We find that the [CII] emission originates primarily from cold, moderate density gas ($40 \lesssim T \lesssim 65$ K and $50 \lesssim n \lesssim 440$ cm$^{-3}$), composed mainly of atomic hydrogen and with an effective visual extinction between $\sim 0.50$ and $\sim 0.91$. Gas dominated by molecular hydrogen contributes only $\lesssim$20% of the total [CII] line emission. Thus, [CII] is not a good tracer for CO-dark H$_2$ at this early phase in the cloud's lifetime. We also find that the total gas, H and C$^+$ column densities are all correlated with the integrated [CII] line emission, with power law slopes ranging from 0.5 to 0.7. Further, the median ratio between the total column density and the [CII] line emission is $Y_{\rm CII}\approx 1.1 \times 10^{21}$ cm$^{-2}$ (K km s$^{-1}$)$^{-1}$, and $Y_{\rm CII}$ scales with $I_{[\textrm{CII}]}^{-0.3}$. We expect $Y_{\rm CII}$ to change in environments with a lower or higher radiation field than simulated here.
△ Less
Submitted 27 September, 2018;
originally announced September 2018.
-
Histogram of oriented gradients: a technique for the study of molecular cloud formation
Authors:
J. D. Soler,
H. Beuther,
M. Rugel,
Y. Wang,
P. C. Clark,
S. C. O. Glover,
P. F. Goldsmith,
M. Heyer,
L. D. Anderson,
A. Goodman,
Th. Henning,
J. Kainulainen,
R. S. Klessen,
S. N. Longmore,
N. M. McClure-Griffiths,
K. M. Menten,
J. C. Mottram,
J. Ott,
S. E. Ragan,
R. J. Smith,
J. S. Urquhart,
F. Bigiel,
P. Hennebelle,
N. Roy,
P. Schilke
Abstract:
We introduce the histogram of oriented gradients (HOG), a tool developed for machine vision that we propose as a new metric for the systematic characterization of observations of atomic and molecular gas and the study of molecular cloud formation models. In essence, the HOG technique takes as input extended spectral-line observations from two tracers and provides an estimate of their spatial corre…
▽ More
We introduce the histogram of oriented gradients (HOG), a tool developed for machine vision that we propose as a new metric for the systematic characterization of observations of atomic and molecular gas and the study of molecular cloud formation models. In essence, the HOG technique takes as input extended spectral-line observations from two tracers and provides an estimate of their spatial correlation across velocity channels.
We characterize HOG using synthetic observations of HI and $^{13}$CO(J=1-0) emission from numerical simulations of MHD turbulence leading to the formation of molecular gas after the collision of two atomic clouds. We find a significant spatial correlation between the two tracers in velocity channels where $v_{HI}\approx v_{^{13}CO}$, independent of the orientation of the collision with respect to the line of sight.
We use HOG to investigate the spatial correlation of the HI, from the THOR survey, and the $^{13}$CO(J=1-0) emission, from the GRS, toward the portion of the Galactic plane 33.75$\lt l\lt$35.25$^{o}$ and $|b|\lt$1.25$^{o}$. We find a significant spatial correlation between the tracers in extended portions of the studied region. Although some of the regions with high spatial correlation are associated with HI self-absorption features, suggesting that it is produced by the cold atomic gas, the correlation is not exclusive to this kind of region. The HOG results also indicate significant differences between individual regions: some show spatial correlation in channels around $v_{HI}\approx v_{^{13}CO}$ while others present this correlation in velocity channels separated by a few km/s. We associate these velocity offsets to the effect of feedback and to the presence of physical conditions that are not included in the atomic-cloud-collision simulations, such as more general magnetic field configurations, shear, and global gas infall.
△ Less
Submitted 2 January, 2019; v1 submitted 21 September, 2018;
originally announced September 2018.
-
On the Origin of the Two-Dimensional Electron Gas at the CdO (100) Surface
Authors:
Pip C. J. Clark,
Andrew I. Williamson,
Nathan K. Lewis,
Ruben Ahumada-Lazo,
Mathieu Silly,
James J. Mudd,
Chris F. McConville,
Wendy R. Flavell
Abstract:
Synchrotron-radiation angle-resolved and core-level photoemission spectroscopy are used together to investigate the origin of the two-dimensional electron gas on the surface of single crystal CdO (100) films. A reduction in the two-dimensional electron density of the surface state is observed under the synchrotron beam during ARPES, which is shown to be accompanied by a concomitant reduction in th…
▽ More
Synchrotron-radiation angle-resolved and core-level photoemission spectroscopy are used together to investigate the origin of the two-dimensional electron gas on the surface of single crystal CdO (100) films. A reduction in the two-dimensional electron density of the surface state is observed under the synchrotron beam during ARPES, which is shown to be accompanied by a concomitant reduction in the surface-adsorbed species (monitored through the O 1s core level signal). This shows that surface adsorbates donate electrons into the surface accumulation layer. When the surface is cleaned, the surface conduction band state empties. A surface doped with atomic H is also studied. Here, interstitial H increases the two-dimensional electron density at the surface. This demonstrates that reversible donor doping is possible. The surface band bending profiles, 2D electron densities, and effective masses are calculated from subband dispersion simulations.
△ Less
Submitted 21 September, 2018;
originally announced September 2018.
-
Tracing the formation of molecular clouds via [CII], [CI] and CO emission
Authors:
Paul C. Clark,
Simon C. O. Glover,
Sarah E. Ragan,
Ana Duarte-Cabral
Abstract:
Our understanding of how molecular clouds form in the interstellar medium (ISM) would be greatly helped if we had a reliable observational tracer of the gas flows responsible for forming the clouds. Fine structure emission from singly ionised and neutral carbon ([CII], [CI]) and rotational line emission from CO are all observed to be associated with molecular clouds. However, it remains unclear wh…
▽ More
Our understanding of how molecular clouds form in the interstellar medium (ISM) would be greatly helped if we had a reliable observational tracer of the gas flows responsible for forming the clouds. Fine structure emission from singly ionised and neutral carbon ([CII], [CI]) and rotational line emission from CO are all observed to be associated with molecular clouds. However, it remains unclear whether any of these tracers can be used to study the inflow of gas onto an assembling cloud, or whether they primarily trace the cloud once it has already assembled. In this paper, we address this issue with the help of high resolution simulations of molecular cloud formation that include a sophisticated treatment of the chemistry and thermal physics of the ISM. Our simulations demonstrate that both [CI] and CO emission trace gas that is predominantly molecular, with a density n ~ 500-1000 cm$^{-3}$, much larger than the density of the inflowing gas. [CII] traces lower density material (n ~ 100 cm$^{-3}$) that is mainly atomic at early times. A large fraction of the [CII] emission traces the same structures as the [CI] or CO emission, but some arises in the inflowing gas. Unfortunately, this emission is very faint and will be difficult to detect with current observational facilities, even for clouds situated in regions with an elevated interstellar radiation field.
△ Less
Submitted 12 May, 2019; v1 submitted 3 September, 2018;
originally announced September 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.
-
CO line ratios in molecular clouds: the impact of environment
Authors:
Camilo H. Peñaloza,
Paul C. Clark,
Simon C. O. Glover,
Ralf S. Klessen
Abstract:
Line emission is strongly dependent on the local environmental conditions in which the emitting tracers reside. In this work, we focus on modelling the CO emission from simulated giant molecular clouds (GMCs), and study the variations in the resulting line ratios arising from the emission from the $J=1-0$, $J=2-1$ and $J=3-2$ transitions. We perform a set of smoothed particle hydrodynamics (SPH) s…
▽ More
Line emission is strongly dependent on the local environmental conditions in which the emitting tracers reside. In this work, we focus on modelling the CO emission from simulated giant molecular clouds (GMCs), and study the variations in the resulting line ratios arising from the emission from the $J=1-0$, $J=2-1$ and $J=3-2$ transitions. We perform a set of smoothed particle hydrodynamics (SPH) simulations with time-dependent chemistry, in which environmental conditions -- including total cloud mass, density, size, velocity dispersion, metallicity, interstellar radiation field (ISRF) and the cosmic ray ionisation rate (CRIR) -- were systematically varied. The simulations were then post-processed using radiative transfer to produce synthetic emission maps in the 3 transitions quoted above. We find that the cloud-averaged values of the line ratios can vary by up to $\pm 0.3$ dex, triggered by changes in the environmental conditions. Changes in the ISRF and/or in the CRIR have the largest impact on line ratios since they directly affect the abundance, temperature and distribution of CO-rich gas within the clouds. We show that the standard methods used to convert CO emission to H$_2$ column density can underestimate the total H$_2$ molecular gas in GMCs by factors of 2 or 3, depending on the environmental conditions in the clouds.
△ Less
Submitted 15 December, 2017; v1 submitted 3 November, 2017;
originally announced November 2017.
-
New ALMA constraints on the star-forming ISM at low metallicity: A 50 pc view of the blue compact dwarf galaxy SBS0335-052
Authors:
D. Cormier,
G. J. Bendo,
S. Hony,
V. Lebouteiller,
S. C. Madden,
F. Galliano,
S. C. O. Glover,
R. S. Klessen,
N. P. Abel,
F. Bigiel,
P. C. Clark
Abstract:
Properties of the cold interstellar medium of low-metallicity galaxies are not well-known due to the faintness and extremely small scale on which emission is expected. We present deep ALMA band 6 (230GHz) observations of the nearby, low-metallicity (12 + log(O/H) = 7.25) blue compact dwarf galaxy SBS0335-052 at an unprecedented resolution of 0.2 arcsec (52 pc). The 12CO J=2-1 line is not detected…
▽ More
Properties of the cold interstellar medium of low-metallicity galaxies are not well-known due to the faintness and extremely small scale on which emission is expected. We present deep ALMA band 6 (230GHz) observations of the nearby, low-metallicity (12 + log(O/H) = 7.25) blue compact dwarf galaxy SBS0335-052 at an unprecedented resolution of 0.2 arcsec (52 pc). The 12CO J=2-1 line is not detected and we report a 3-sigma upper limit of LCO(2-1) = 3.6x10^4 K km/s pc^2. Assuming that molecular gas is converted into stars with a given depletion time, ranging from 0.02 to 2 Gyr, we find lower limits on the CO-to-H2 conversion factor alpha_CO in the range 10^2-10^4 Msun pc^-2 (K km/s)^-1. The continuum emission is detected and resolved over the two main super star clusters. Re-analysis of the IR-radio spectral energy distribution suggests that the mm-fluxes are not only free-free emission but are most likely also associated with a cold dust component coincident with the position of the brightest cluster. With standard dust properties, we estimate its mass to be as large as 10^5 Msun. Both line and continuum results suggest the presence of a large cold gas reservoir unseen in CO even with ALMA.
△ Less
Submitted 7 April, 2017; v1 submitted 27 February, 2017;
originally announced February 2017.
-
Variable interstellar radiation fields in simulated dwarf galaxies: supernovae versus photoelectric heating
Authors:
Chia-Yu Hu,
Thorsten Naab,
Simon C. O. Glover,
Stefanie Walch,
Paul C. Clark
Abstract:
We present high-resolution hydrodynamical simulations of isolated dwarf galaxies including self-gravity, non-equilibrium cooling and chemistry, interstellar radiation fields (IRSF) and shielding, star formation, and stellar feedback. This includes spatially and temporally varying photoelectric (PE) heating, photoionization, resolved supernova (SN) blast waves and metal enrichment. A new flexible m…
▽ More
We present high-resolution hydrodynamical simulations of isolated dwarf galaxies including self-gravity, non-equilibrium cooling and chemistry, interstellar radiation fields (IRSF) and shielding, star formation, and stellar feedback. This includes spatially and temporally varying photoelectric (PE) heating, photoionization, resolved supernova (SN) blast waves and metal enrichment. A new flexible method to sample the stellar initial mass function allows us to follow the contribution to the ISRF, the metal output and the SN delay times of individual massive stars. We find that SNe play the dominant role in regulating the global star formation rate, shaping the multi-phase interstellar medium (ISM) and driving galactic outflows. Outflow rates (with mass-loading factors of a few) and hot gas fractions of the ISM increase with the number of SNe exploding in low-density environments where radiative energy losses are low. While PE heating alone can suppress star formation slightly more (a factor of a few) than SNe alone can do, it is unable to drive outflows and reproduce the multi-phase ISM that emerges naturally when SNe are included. These results are in conflict with recent results of Forbes et al. who concluded that PE heating is the dominant process suppressing star formation in dwarfs, about an order of magnitude more efficient than SNe. Potential origins for this discrepancy are discussed. In the absence of SNe and photoionization (mechanisms to disperse dense clouds), the impact of PE heating is highly overestimated owing to the (unrealistic) proximity of dense gas to the radiation sources. This leads to a substantial boost of the infrared continuum emission from the UV-irradiated dust and a far infrared line-to-continuum ratio too low compared to observations. Though sub-dominant in regulating star formation, the ISRF controls the abundance of molecular hydrogen via photodissociation.
△ Less
Submitted 30 January, 2017;
originally announced January 2017.
-
Using CO line ratios to trace the physical properties of molecular clouds
Authors:
Camilo H. Peñaloza,
Paul C. Clark,
Simon C. O. Glover,
Rahul Shetty,
Ralf S. Klessen
Abstract:
The carbon monoxide (CO) rotational transition lines are the most common tracers of molecular gas within giant molecular clouds (MCs). We study the ratio ($R_{2-1/1-0}$) between CO's first two emission lines and examine what information it provides about the physical properties of the cloud. To study $R_{2-1/1-0}$ we perform smooth particle hydrodynamic simulations with time dependent chemistry (u…
▽ More
The carbon monoxide (CO) rotational transition lines are the most common tracers of molecular gas within giant molecular clouds (MCs). We study the ratio ($R_{2-1/1-0}$) between CO's first two emission lines and examine what information it provides about the physical properties of the cloud. To study $R_{2-1/1-0}$ we perform smooth particle hydrodynamic simulations with time dependent chemistry (using GADGET-2), along with post-process radiative transfer calculations on an adaptive grid (using RADMC-3D) to create synthetic emission maps of a MC. $R_{2-1/1-0}$ has a bimodal distribution that is a consequence of the excitation properties of each line, given that $J=1$ reaches local thermal equilibrium (LTE) while $J=2$ is still sub-thermally excited in the considered clouds. The bimodality of $R_{2-1/1-0}$ serves as a tracer of the physical properties of different regions of the cloud and it helps constrain local temperatures, densities and opacities. Additionally this bimodal structure shows an important portion of the CO emission comes from diffuse regions of the cloud, suggesting that the commonly used conversion factor of $R_{2-1/1-0}\sim 0.7$ between both lines may need to be studied further.
△ Less
Submitted 7 November, 2016;
originally announced November 2016.
-
The turbulent life of dust grains in the supernova-driven, multi-phase interstellar medium
Authors:
Thomas Peters,
Svitlana Zhukovska,
Thorsten Naab,
Philipp Girichidis,
Stefanie Walch,
Simon C. O. Glover,
Ralf S. Klessen,
Paul C. Clark,
Daniel Seifried
Abstract:
Dust grains are an important component of the interstellar medium (ISM) of galaxies. We present the first direct measurement of the residence times of interstellar dust in the different ISM phases, and of the transition rates between these phases, in realistic hydrodynamical simulations of the multi-phase ISM. Our simulations include a time-dependent chemical network that follows the abundances of…
▽ More
Dust grains are an important component of the interstellar medium (ISM) of galaxies. We present the first direct measurement of the residence times of interstellar dust in the different ISM phases, and of the transition rates between these phases, in realistic hydrodynamical simulations of the multi-phase ISM. Our simulations include a time-dependent chemical network that follows the abundances of H^+, H, H_2, C^+ and CO and take into account self-shielding by gas and dust using a tree-based radiation transfer method. Supernova explosions are injected either at random locations, at density peaks, or as a mixture of the two. For each simulation, we investigate how matter circulates between the ISM phases and find more sizeable transitions than considered in simple mass exchange schemes in the literature. The derived residence times in the ISM phases are characterised by broad distributions, in particular for the molecular, warm and hot medium. The most realistic simulations with random and mixed driving have median residence times in the molecular, cold, warm and hot phase around 17, 7, 44 and 1 Myr, respectively. The transition rates measured in the random driving run are in good agreement with observations of Ti gas-phase depletion in the warm and cold phases in a simple depletion model, although the depletion in the molecular phase is under-predicted. ISM phase definitions based on chemical abundance rather than temperature cuts are physically more meaningful, but lead to significantly different transition rates and residence times because there is no direct correspondence between the two definitions.
△ Less
Submitted 20 October, 2016;
originally announced October 2016.
-
The IMF as a function of supersonic turbulence
Authors:
Clio Bertelli Motta,
Paul C. Clark,
Simon C. O. Glover,
Ralf S. Klessen,
Anna Pasquali
Abstract:
Recent studies seem to suggest that the stellar initial mass function (IMF) in early-type galaxies might be different from a classical Kroupa or Chabrier IMF, i.e. contain a larger fraction of the total mass in low-mass stars. From a theoretical point of view, supersonic turbulence has been the subject of interest in many analytical theories proposing a strong correlation with the characteristic m…
▽ More
Recent studies seem to suggest that the stellar initial mass function (IMF) in early-type galaxies might be different from a classical Kroupa or Chabrier IMF, i.e. contain a larger fraction of the total mass in low-mass stars. From a theoretical point of view, supersonic turbulence has been the subject of interest in many analytical theories proposing a strong correlation with the characteristic mass of the core mass function (CMF) in star forming regions, and as a consequence with the stellar IMF. Performing two suites of smoothed particles hydrodynamics (SPH) simulations with different mass resolutions, we aim at testing the effects of variations in the turbulent properties of a dense, star forming molecular cloud on the shape of the system mass function in different density regimes. While analytical theories predict a shift of the peak of the CMF towards lower masses with increasing velocity dispersion of the cloud, we observe in the low-density regime the opposite trend, with high Mach numbers giving rise to a top-heavy mass distribution. For the high-density regime we do not find any trend correlating the Mach number with the characteristic mass of the resulting IMF, implying that the dynamics of protostellar accretion discs and fragmentation on small scales is not strongly affected by turbulence driven at the scale of the cloud. Furthermore, we suggest that a significant fraction of dense cores are disrupted by turbulence before stars can be formed in their interior through gravitational collapse. Although this particular study has limitations in its numerical resolution, we suggest that our results, along with those from other studies, cast doubt on the turbulent fragmentation models on the IMF that simply map the CMF to the IMF.
△ Less
Submitted 3 August, 2016;
originally announced August 2016.
-
The SILCC project: III. Regulation of star formation and outflows by stellar winds and supernovae
Authors:
A. Gatto,
S. Walch,
T. Naab,
P. Girichidis,
R. Wünsch,
S. C. O. Glover,
R. S. Klessen,
P. C. Clark,
T. Peters,
D. Derigs,
C. Baczynski,
J. Puls
Abstract:
We study the impact of stellar winds and supernovae on the multi-phase interstellar medium using three-dimensional hydrodynamical simulations carried out with FLASH. The selected galactic disc region has a size of (500 pc)$^2$ x $\pm$ 5 kpc and a gas surface density of 10 M$_{\odot}$/pc$^2$. The simulations include an external stellar potential and gas self-gravity, radiative cooling and diffuse h…
▽ More
We study the impact of stellar winds and supernovae on the multi-phase interstellar medium using three-dimensional hydrodynamical simulations carried out with FLASH. The selected galactic disc region has a size of (500 pc)$^2$ x $\pm$ 5 kpc and a gas surface density of 10 M$_{\odot}$/pc$^2$. The simulations include an external stellar potential and gas self-gravity, radiative cooling and diffuse heating, sink particles representing star clusters, stellar winds from these clusters which combine the winds from indi- vidual massive stars by following their evolution tracks, and subsequent supernova explosions. Dust and gas (self-)shielding is followed to compute the chemical state of the gas with a chemical network. We find that stellar winds can regulate star (cluster) formation. Since the winds suppress the accretion of fresh gas soon after the cluster has formed, they lead to clusters which have lower average masses (10$^2$ - 10$^{4.3}$ M$_{\odot}$) and form on shorter timescales (10$^{-3}$ - 10 Myr). In particular we find an anti-correlation of cluster mass and accretion time scale. Without winds the star clusters easily grow to larger masses for ~5 Myr until the first supernova explodes. Overall the most massive stars provide the most wind energy input, while objects beginning their evolution as B-type stars contribute most of the supernova energy input. A significant outflow from the disk (mass loading $\gtrsim$ 1 at 1 kpc) can be launched by thermal gas pressure if more than 50% of the volume near the disc mid-plane can be heated to T > 3x10$^5$ K. Stellar winds alone cannot create a hot volume-filling phase. The models which are in best agreement with observed star formation rates drive either no outflows or weak outflows.
△ Less
Submitted 16 June, 2016;
originally announced June 2016.
-
The role of cosmic ray pressure in accelerating galactic outflows
Authors:
Christine M. Simpson,
Ruediger Pakmor,
Federico Marinacci,
Christoph Pfrommer,
Volker Springel,
Simon C. O. Glover,
Paul C. Clark,
Rowan J. Smith
Abstract:
We study the formation of galactic outflows from supernova explosions (SNe) with the moving-mesh code AREPO in a stratified column of gas with a surface density similar to the Milky Way disk at the solar circle. We compare different simulation models for SNe placement and energy feedback, including cosmic rays (CR), and find that models that place SNe in dense gas and account for CR diffusion are…
▽ More
We study the formation of galactic outflows from supernova explosions (SNe) with the moving-mesh code AREPO in a stratified column of gas with a surface density similar to the Milky Way disk at the solar circle. We compare different simulation models for SNe placement and energy feedback, including cosmic rays (CR), and find that models that place SNe in dense gas and account for CR diffusion are able to drive outflows with similar mass loading as obtained from a random placement of SNe with no CRs. Despite this similarity, CR-driven outflows differ in several other key properties including their overall clumpiness and velocity. Moreover, the forces driving these outflows originate in different sources of pressure, with the CR diffusion model relying on non-thermal pressure gradients to create an outflow driven by internal pressure and the random-placement model depending on kinetic pressure gradients to propel a ballistic outflow. CRs therefore appear to be non-negligible physics in the formation of outflows from the interstellar medium.
△ Less
Submitted 16 August, 2016; v1 submitted 7 June, 2016;
originally announced June 2016.
-
Evidence that widespread star formation may be underway in G0.253+016, "The Brick"
Authors:
K. A. Marsh,
S. E. Ragan,
A. P. Whitworth,
P. C. Clark
Abstract:
Image cubes of differential column density as a function of dust temperature are constructed for Galactic Centre molecular cloud G0.253+0.016 ("The Brick") using the recently described PPMAP procedure. The input data consist of continuum images from the Herschel Space Telescope in the wavelength range 70-500 $μ$m, supplemented by previously published interferometric data at 1.3 mm wavelength. Whil…
▽ More
Image cubes of differential column density as a function of dust temperature are constructed for Galactic Centre molecular cloud G0.253+0.016 ("The Brick") using the recently described PPMAP procedure. The input data consist of continuum images from the Herschel Space Telescope in the wavelength range 70-500 $μ$m, supplemented by previously published interferometric data at 1.3 mm wavelength. While the bulk of the dust in the molecular cloud is consistent with being heated externally by the local interstellar radiation field, our image cube shows the presence, near one edge of the cloud, of a filamentary structure whose temperature profile suggests internal heating. The structure appears as a cool ($\sim 14$ K) tadpole-like feature, $\sim 6$ pc in length, in which is embedded a thin spine of much hotter ($\sim$ 40-50 K) material. We interpret these findings in terms of a cool filament whose hot central region is undergoing gravitational collapse and fragmentation to form a line of protostars. If confirmed, this would represent the first evidence of widespread star formation having started within this cloud.
△ Less
Submitted 26 April, 2016;
originally announced April 2016.
-
Hunting for the signatures of molecular cloud formation
Authors:
Simon C. O. Glover,
Paul C. Clark
Abstract:
In order to understand how molecular clouds form in the Galactic interstellar medium, we would like to be able to map the structure and kinematics of the gas flows responsible for forming them. However, doing so is observationally challenging. CO, the workhorse molecule for studies of molecular clouds, traces only relatively dense gas and hence only allows us to study those portions of the clouds…
▽ More
In order to understand how molecular clouds form in the Galactic interstellar medium, we would like to be able to map the structure and kinematics of the gas flows responsible for forming them. However, doing so is observationally challenging. CO, the workhorse molecule for studies of molecular clouds, traces only relatively dense gas and hence only allows us to study those portions of the clouds that have already assembled. Numerical simulations suggest that the inflowing gas that forms these clouds is largely composed of CO-dark H2. These same simulations allow us to explore the usefulness of different tracers of this CO-dark molecular material, and we use them here to show that the [CII] fine structure line is potentially a very powerful tracer of this gas and should be readily detectable using modern instrumentation.
△ Less
Submitted 15 November, 2015;
originally announced November 2015.
-
Impact of supernova and cosmic-ray driving on the surface brightness of the galactic halo in soft X-rays
Authors:
Thomas Peters,
Philipp Girichidis,
Andrea Gatto,
Thorsten Naab,
Stefanie Walch,
Richard Wünsch,
Simon C. O. Glover,
Paul C. Clark,
Ralf S. Klessen,
Christian Baczynski
Abstract:
The halo of the Milky Way contains a hot plasma with a surface brightness in soft X-rays of the order $10^{-12}$erg cm$^{-2}$ s$^{-1}$ deg$^{-2}$. The origin of this gas is unclear, but so far numerical models of galactic star formation have failed to reproduce such a large surface brightness by several orders of magnitude. In this paper, we analyze simulations of the turbulent, magnetized, multi-…
▽ More
The halo of the Milky Way contains a hot plasma with a surface brightness in soft X-rays of the order $10^{-12}$erg cm$^{-2}$ s$^{-1}$ deg$^{-2}$. The origin of this gas is unclear, but so far numerical models of galactic star formation have failed to reproduce such a large surface brightness by several orders of magnitude. In this paper, we analyze simulations of the turbulent, magnetized, multi-phase interstellar medium including thermal feedback by supernova explosions as well as cosmic-ray feedback. We include a time-dependent chemical network, self-shielding by gas and dust, and self-gravity. Pure thermal feedback alone is sufficient to produce the observed surface brightness, although it is very sensitive to the supernova rate. Cosmic rays suppress this sensitivity and reduce the surface brightness because they drive cooler outflows. Self-gravity has by far the largest effect because it accumulates the diffuse gas in the disk in dense clumps and filaments, so that supernovae exploding in voids can eject a large amount of hot gas into the halo. This can boost the surface brightness by several orders of magnitude. Although our simulations do not reach a steady state, all simulations produce surface brightness values of the same order of magnitude as the observations, with the exact value depending sensitively on the simulation parameters. We conclude that star formation feedback alone is sufficient to explain the origin of the hot halo gas, but measurements of the surface brightness alone do not provide useful diagnostics for the study of galactic star formation.
△ Less
Submitted 22 October, 2015;
originally announced October 2015.
-
Star formation and molecular hydrogen in dwarf galaxies: a non-equilibrium view
Authors:
Chia-Yu Hu,
Thorsten Naab,
Stefanie Walch,
Simon C. O. Glover,
Paul C. Clark
Abstract:
We study the connection of star formation to atomic (HI) and molecular hydrogen (H$_2$) in isolated, low metallicity dwarf galaxies with high-resolution ($m_{\rm gas}$ = 4 M$_\odot$, $N_{\rm ngb}$ = 100) SPH simulations. The model includes self-gravity, non-equilibrium cooling, shielding from an interstellar radiation field, the chemistry of H$_2$ formation, H$_2$-independent star formation, super…
▽ More
We study the connection of star formation to atomic (HI) and molecular hydrogen (H$_2$) in isolated, low metallicity dwarf galaxies with high-resolution ($m_{\rm gas}$ = 4 M$_\odot$, $N_{\rm ngb}$ = 100) SPH simulations. The model includes self-gravity, non-equilibrium cooling, shielding from an interstellar radiation field, the chemistry of H$_2$ formation, H$_2$-independent star formation, supernova feedback and metal enrichment. We find that the H$_2$ mass fraction is sensitive to the adopted dust-to-gas ratio and the strength of the interstellar radiation field, while the star formation rate is not. Star formation is regulated by stellar feedback, keeping the gas out of thermal equilibrium for densities $n <$ 1 cm$^{-3}$. Because of the long chemical timescales, the H$_2$ mass remains out of chemical equilibrium throughout the simulation. Star formation is well-correlated with cold ( T $\leqslant$ 100 K ) gas, but this dense and cold gas - the reservoir for star formation - is dominated by HI, not H$_2$. In addition, a significant fraction of H$_2$ resides in a diffuse, warm phase, which is not star-forming. The ISM is dominated by warm gas (100 K $<$ T $\leqslant 3\times 10^4$ K) both in mass and in volume. The scale height of the gaseous disc increases with radius while the cold gas is always confined to a thin layer in the mid-plane. The cold gas fraction is regulated by feedback at small radii and by the assumed radiation field at large radii. The decreasing cold gas fractions result in a rapid increase in depletion time (up to 100 Gyrs) for total gas surface densities $Σ_{\rm HI+H_2} \lesssim$ 10 M$_\odot$pc$^{-2}$, in agreement with observations of dwarf galaxies in the Kennicutt-Schmidt plane.
△ Less
Submitted 7 March, 2016; v1 submitted 19 October, 2015;
originally announced October 2015.
-
Launching cosmic-ray-driven outflows from the magnetized interstellar medium
Authors:
Philipp Girichidis,
Thorsten Naab,
Stefanie Walch,
Michal Hanasz,
Mordecai-Mark Mac Low,
Jeremiah P. Ostriker,
Andrea Gatto,
Thomas Peters,
Richard Wünsch,
Simon C. O. Glover,
Ralf S. Klessen,
Paul C. Clark,
Christian Baczynski
Abstract:
We present a hydrodynamical simulation of the turbulent, magnetized, supernova (SN)-driven interstellar medium (ISM) in a stratified box that dynamically couples the injection and evolution of cosmic rays (CRs) and a self-consistent evolution of the chemical composition. CRs are treated as a relativistic fluid in the advection-diffusion approximation. The thermodynamic evolution of the gas is comp…
▽ More
We present a hydrodynamical simulation of the turbulent, magnetized, supernova (SN)-driven interstellar medium (ISM) in a stratified box that dynamically couples the injection and evolution of cosmic rays (CRs) and a self-consistent evolution of the chemical composition. CRs are treated as a relativistic fluid in the advection-diffusion approximation. The thermodynamic evolution of the gas is computed using a chemical network that follows the abundances of H+, H, H2, CO, C+, and free electrons and includes (self-)shielding of the gas and dust. We find that CRs perceptibly thicken the disk with the heights of 90% (70%) enclosed mass reaching ~1.5 kpc (~0.2 kpc). The simulations indicate that CRs alone can launch and sustain strong outflows of atomic and ionized gas with mass loading factors of order unity, even in solar neighborhood conditions and with a CR energy injection per SN of 10^50 erg, 10% of the fiducial thermal energy of an SN. The CR-driven outflows have moderate launching velocities close to the midplane (~100 km/s) and are denser (ρ~1e-24 - 1e-26 g/cm^3), smoother, and colder than the (thermal) SN-driven winds. The simulations support the importance of CRs for setting the vertical structure of the disk as well as the driving of winds.
△ Less
Submitted 6 January, 2016; v1 submitted 24 September, 2015;
originally announced September 2015.
-
Synthetic observations of molecular clouds in a galactic center environment: I. Studying maps of column density and integrated intensity
Authors:
Erik Bertram,
Simon C. O. Glover,
Paul C. Clark,
Sarah E. Ragan,
Ralf S. Klessen
Abstract:
We run numerical simulations of molecular clouds (MCs), adopting properties similar to those found in the Central Molecular Zone (CMZ) of the Milky Way. For this, we employ the moving mesh code Arepo and perform simulations which account for a simplified treatment of time-dependent chemistry and the non-isothermal nature of gas and dust. We perform simulations using an initial density of n_0 = 10^…
▽ More
We run numerical simulations of molecular clouds (MCs), adopting properties similar to those found in the Central Molecular Zone (CMZ) of the Milky Way. For this, we employ the moving mesh code Arepo and perform simulations which account for a simplified treatment of time-dependent chemistry and the non-isothermal nature of gas and dust. We perform simulations using an initial density of n_0 = 10^3 cm^{-3} and a mass of 1.3x10^5 M_sun. Furthermore, we vary the virial parameter, defined as the ratio of kinetic and potential energy, alpha = E_{kin} / |E_{pot}|, by adjusting the velocity dispersion. We set it to alpha = 0.5, 2.0 and 8.0, in order to analyze the impact of the kinetic energy on our results. We account for the extreme conditions in the CMZ and increase both the interstellar radiation field (ISRF) and the cosmic-ray flux (CRF) by a factor of 1000 compared to the values found in the solar neighbourhood. We use the radiative transfer code RADMC-3D to compute synthetic images in various diagnostic lines. These are [CII] at 158 micron, [OI] (145 micron), [OI] (63 micron), 12CO (J = 1 -> 0) and 13CO (J = 1 -> 0) at 2600 micron and 2720 micron, respectively. When alpha is large, the turbulence disperses much of the gas in the cloud, reducing its mean density and allowing the ISRF to penetrate more deeply into the cloud's interior. This significantly alters the chemical composition of the cloud, leading to the dissociation of a significant amount of the molecular gas. On the other hand, when alpha is small, the cloud remains compact, allowing more of the molecular gas to survive. We show that in each case the atomic tracers accurately reflect most of the physical properties of both the H2 and the total gas of the cloud and that they provide a useful alternative to molecular lines when studying the ISM in the CMZ.
△ Less
Submitted 4 November, 2015; v1 submitted 21 September, 2015;
originally announced September 2015.
-
Is atomic carbon a good tracer of molecular gas in metal-poor galaxies?
Authors:
Simon C. O. Glover,
Paul C. Clark
Abstract:
Carbon monoxide (CO) is widely used as a tracer of molecular hydrogen (H2) in metal-rich galaxies, but is known to become ineffective in low metallicity dwarf galaxies. Atomic carbon has been suggested as a superior tracer of H2 in these metal-poor systems, but its suitability remains unproven. To help us to assess how well atomic carbon traces H2 at low metallicity, we have performed a series of…
▽ More
Carbon monoxide (CO) is widely used as a tracer of molecular hydrogen (H2) in metal-rich galaxies, but is known to become ineffective in low metallicity dwarf galaxies. Atomic carbon has been suggested as a superior tracer of H2 in these metal-poor systems, but its suitability remains unproven. To help us to assess how well atomic carbon traces H2 at low metallicity, we have performed a series of numerical simulations of turbulent molecular clouds that cover a wide range of different metallicities. Our simulations demonstrate that in star-forming clouds, the conversion factor between [CI] emission and H2 mass, $X_{\rm CI}$, scales approximately as $X_{\rm CI} \propto Z^{-1}$. We recover a similar scaling for the CO-to-H2 conversion factor, $X_{\rm CO}$, but find that at this point in the evolution of the clouds, $X_{\rm CO}$ is consistently smaller than $X_{\rm CI}$, by a factor of a few or more. We have also examined how $X_{\rm CI}$ and $X_{\rm CO}$ evolve with time. We find that $X_{\rm CI}$ does not vary strongly with time, demonstrating that atomic carbon remains a good tracer of H2 in metal-poor systems even at times significantly before the onset of star formation. On the other hand, $X_{\rm CO}$ varies very strongly with time in metal-poor clouds, showing that CO does not trace H2 well in starless clouds at low metallicity.
△ Less
Submitted 4 December, 2015; v1 submitted 7 September, 2015;
originally announced September 2015.
-
The SILCC (SImulating the LifeCycle of molecular Clouds) project - II. Dynamical evolution of the supernova-driven ISM and the launching of outflows
Authors:
Philipp Girichidis,
Stefanie Walch,
Thorsten Naab,
Andrea Gatto,
Richard Wünsch,
Simon C. O. Glover,
Ralf S. Klessen,
Paul C. Clark,
Thomas Peters,
Dominik Derigs,
Christian Baczynski
Abstract:
The SILCC project (SImulating the Life-Cycle of molecular Clouds) aims at a more self-consistent understanding of the interstellar medium (ISM) on small scales and its link to galaxy evolution. We present three-dimensional (magneto)hydrodynamic simulations of the ISM in a vertically stratified box including self-gravity, an external potential due to the stellar component of the galactic disc, and…
▽ More
The SILCC project (SImulating the Life-Cycle of molecular Clouds) aims at a more self-consistent understanding of the interstellar medium (ISM) on small scales and its link to galaxy evolution. We present three-dimensional (magneto)hydrodynamic simulations of the ISM in a vertically stratified box including self-gravity, an external potential due to the stellar component of the galactic disc, and stellar feedback in the form of an interstellar radiation field and supernovae (SNe). The cooling of the gas is based on a chemical network that follows the abundances of H+, H, H2, C+, and CO and takes shielding into account consistently. We vary the SN feedback by comparing different SN rates, clustering and different positioning, in particular SNe in density peaks and at random positions, which has a major impact on the dynamics. Only for random SN positions the energy is injected in sufficiently low-density environments to reduce energy losses and enhance the effective kinetic coupling of the SNe with the gas. This leads to more realistic velocity dispersions (σ_HI ~ 0.8σ_(300-8000K) ~ 10-20km/s, σ_Hα~ 0.6σ_(8000-3e5K) ~ 20-30km/s), and strong outflows with mass loading factors of up to 10 even for solar neighbourhood conditions. Clustered SNe abet the onset of outflows compared to individual SNe but do not influence the net outflow rate. The outflows do not contain any molecular gas and are mainly composed of atomic hydrogen. The bulk of the outflowing mass is dense (ρ~ 1e-25-1e-24g/cc) and slow (v ~ 20-40km/s) but there is a high-velocity tail of up to v ~ 500km/s with ρ~ 1e-28-1e-27g/cc.
△ Less
Submitted 11 December, 2015; v1 submitted 26 August, 2015;
originally announced August 2015.
-
Does the CO-to-H2 conversion factor depend on the star formation rate?
Authors:
Paul C. Clark,
Simon C. O. Glover
Abstract:
We present a series of numerical simulations that explore how the `X-factor', $X_{CO}$ -- the conversion factor between the observed integrated CO emission and the column density of molecular hydrogen -- varies with the environmental conditions in which a molecular cloud is placed. Our investigation is centred around two environmental conditions in particular: the cosmic ray ionisation rate (CRIR)…
▽ More
We present a series of numerical simulations that explore how the `X-factor', $X_{CO}$ -- the conversion factor between the observed integrated CO emission and the column density of molecular hydrogen -- varies with the environmental conditions in which a molecular cloud is placed. Our investigation is centred around two environmental conditions in particular: the cosmic ray ionisation rate (CRIR) and the strength of the interstellar radiation field (ISRF). Since both these properties of the interstellar medium have their origins in massive stars, we make the assumption in this paper that both the strength of the ISRF and the CRIR scale linearly with the local star formation rate (SFR). The cloud modelling in this study first involves running numerical simulations that capture the cloud dynamics, as well as the time-dependent chemistry, and ISM heating and cooling. These simulations are then post-processed with a line radiative transfer code to create synthetic 12CO (1-0) emission maps from which $X_{CO}$ can be calculated. We find that for 1e4 solar mass virialised clouds with mean density 100 cm$^{-3}$, $X_{CO}$ is only weakly dependent on the local SFR, varying by a factor of a few over two orders of magnitude in SFR. In contrast, we find that for similar clouds but with masses of 1e5 solar masses, the X-factor will vary by an order of magnitude over the same range in SFR, implying that extra-galactic star formation laws should be viewed with caution. However, for denser ($10^4$ cm$^{-3}$), super-virial clouds such as those found at the centre of the Milky Way, the X-factor is once again independent of the local SFR.
△ Less
Submitted 22 June, 2015;
originally announced June 2015.
-
Star formation efficiencies of molecular clouds in a galactic center environment
Authors:
Erik Bertram,
Simon C. O. Glover,
Paul C. Clark,
Ralf S. Klessen
Abstract:
We use the Arepo moving mesh code to simulate the evolution of molecular clouds exposed to a harsh environment similar to that found in the galactic center (GC), in an effort to understand why the star formation efficiency (SFE) of clouds in this environment is so small. Our simulations include a simplified treatment of time-dependent chemistry and account for the highly non-isothermal nature of t…
▽ More
We use the Arepo moving mesh code to simulate the evolution of molecular clouds exposed to a harsh environment similar to that found in the galactic center (GC), in an effort to understand why the star formation efficiency (SFE) of clouds in this environment is so small. Our simulations include a simplified treatment of time-dependent chemistry and account for the highly non-isothermal nature of the gas and the dust. We model clouds with a total mass of 1.3x10^5 M_{sun} and explore the effects of varying the mean cloud density and the virial parameter, alpha = E_{kin}/|E_{pot}|. We vary the latter from alpha = 0.5 to alpha = 8.0, and so many of the clouds that we simulate are gravitationally unbound. We expose our model clouds to an interstellar radiation field (ISRF) and cosmic ray flux (CRF) that are both a factor of 1000 higher than the values found in the solar neighbourhood. As a reference, we also run simulations with local solar neighbourhood values of the ISRF and the CRF in order to better constrain the effects of the extreme conditions in the GC on the SFE. Despite the harsh environment and the large turbulent velocity dispersions adopted, we find that all of the simulated clouds form stars within less than a gravitational free-fall time. Increasing the virial parameter from alpha = 0.5 to alpha = 8.0 decreases the SFE by a factor ~4-10, while increasing the ISRF/CRF by a factor of 1000 decreases the SFE again by a factor ~2-6. However, even in our most unbound clouds, the SFE remains higher than that inferred for real GC clouds. We therefore conclude that high levels of turbulence and strong external heating are not enough by themselves to lead to a persistently low SFE at the center of the Galaxy.
△ Less
Submitted 1 June, 2015; v1 submitted 17 May, 2015;
originally announced May 2015.