Showing 1–50 of 83 results for author: Dale, J

Signatures of mass segregation from competitive accretion and monolithic collapse

Authors: Richard J. Parker, Emily J. Pinson, Hayley L. Alcock, James E. Dale

Abstract: The two main competing theories proposed to explain the formation of massive ($>10$M$_\odot$) stars -- competitive accretion and monolithic core collapse -- make different observable predictions for the environment of the massive stars during, and immediately after, their formation. Proponents of competitive accretion have long predicted that the most massive stars should have a different spatial… ▽ More The two main competing theories proposed to explain the formation of massive ($>10$M$_\odot$) stars -- competitive accretion and monolithic core collapse -- make different observable predictions for the environment of the massive stars during, and immediately after, their formation. Proponents of competitive accretion have long predicted that the most massive stars should have a different spatial distribution to lower-mass stars, either through the stars being mass segregated, or being in areas of higher relative densities, or sitting deeper in gravitational potential wells. We test these predictions by analysing a suite of SPH simulations where star clusters form massive stars via competitive accretion with and without feedback. We find that the most massive stars have higher relative densities, and sit in deeper potential wells, only in simulations in which feedback is not present. When feedback is included, only half of the simulations have the massive stars residing in deeper potential wells, and there are no other distinguishing signals in their spatial distributions. Intriguingly, in our simple models for monolithic core collapse, the massive stars may also end up in deeper potential wells, because if massive cores fragment the stars are still massive, and dominate their local environs. We find no robust diagnostic test in the spatial distributions of massive stars that can distinguish their formation mechanisms, and so other predictions for distinguishing between competitive accretion and monolithic collapse are required. △ Less

Submitted 7 August, 2024; originally announced August 2024.

Comments: 18 pages, 7 figures, with an Appendix containing a further 18 figures showing all the SPH simulation results. Accepted for publication in ApJ

Star Cluster Formation in Clouds with Externally Driven Turbulence

Authors: Jamie D. Smith, James E. Dale, Sarah E. Jaffa, Martin G. H. Krause

Abstract: Star clusters are known to be formed in turbulent molecular clouds. How turbulence is driven in molecular clouds and what effect this has on star formation is still unclear. We compare a simulation setup with turbulent driving everywhere in a periodic box with a setup where turbulence is only driven around the outside of the box. We analyse the resulting gas distribution, kinematics, and the popul… ▽ More Star clusters are known to be formed in turbulent molecular clouds. How turbulence is driven in molecular clouds and what effect this has on star formation is still unclear. We compare a simulation setup with turbulent driving everywhere in a periodic box with a setup where turbulence is only driven around the outside of the box. We analyse the resulting gas distribution, kinematics, and the population of stars that are formed from the cloud. Both setups successfully produce a turbulent velocity field with a power law structure function, the externally driven cloud has a more central, monolithic, clump, while the fully driven cloud has many smaller, more dispersed, clumps. The star formation follows the cloud morphology producing large clusters, with high star forming efficiency in the externally driven simulations and sparse individual star formation with much lower star formation efficiency in the fully driven case. We conclude that the externally driven method, which resembles a Global Hierarchical Collapse (GHC) scenario, produces star clusters that more closely match with observations. △ Less

Submitted 10 August, 2022; originally announced August 2022.

Comments: 8 pages, 6 figures, accepted by MNRAS

Collision Between Molecular Clouds-III: The effects of cloud initial density profile on head-on collisions

Authors: Tabassum S Tanvir, James E Dale

Abstract: In this third paper in our cloud-collision series, we present the results from simulations of head-on collisions with a strongly centrally-condensed initial density profile of $\rm ρ\propto R^{-2}$. We investigate the impact of these density profiles on the overall evolution of the simulations: the structures formed, their dynamical evolution, and their star formation activity. We consider clouds… ▽ More In this third paper in our cloud-collision series, we present the results from simulations of head-on collisions with a strongly centrally-condensed initial density profile of $\rm ρ\propto R^{-2}$. We investigate the impact of these density profiles on the overall evolution of the simulations: the structures formed, their dynamical evolution, and their star formation activity. We consider clouds which are globally bound and globally unbound, leading to three different scenarios -- the collision of a bound cloud with a bound cloud, the collision of two unbound clouds, or the collision of one cloud of each type. In all the simulations dense star clusters form before the collisions occur, and we find that star formation remains confined to these systems and is little affected by the collisions. If the clouds' are both initially bound, the collision forms a filamentary structure, but otherwise, this does not occur. We observe that rotating structures form around the clusters, but they also form in our non-colliding control simulations, so are not a consequence of the collisions. Dissipation of kinetic energy in these simulations is inefficient because of the substructure created in the clouds by turbulence before the collisions. As a result, although some gas is left bound in the COM frame, the star clusters formed in the simulations do not become bound to each other. △ Less

Submitted 10 May, 2021; originally announced May 2021.

Comments: 16 Pages, 21 Figures; Accepted by MNRAS

The complex multi-scale structure in simulated and observed emission maps of the proto-cluster cloud G0.253+0.016 (`the Brick')

Authors: Maya A. Petkova, J. M. Diederik Kruijssen, A. Louise Kluge, Simon C. O. Glover, Daniel L. Walker, Steven N. Longmore, Jonathan D. Henshaw, Stefan Reissl, James E. Dale

Abstract: The Central Molecular Zone (CMZ; the central ~ 500 pc of the Milky Way) hosts molecular clouds in an extreme environment of strong shear, high gas pressure and density, and complex chemistry. G0.253+0.016, also known as `the Brick', is the densest, most compact and quiescent of these clouds. High-resolution observations with the Atacama Large Millimeter/submillimeter Array (ALMA) have revealed its… ▽ More The Central Molecular Zone (CMZ; the central ~ 500 pc of the Milky Way) hosts molecular clouds in an extreme environment of strong shear, high gas pressure and density, and complex chemistry. G0.253+0.016, also known as `the Brick', is the densest, most compact and quiescent of these clouds. High-resolution observations with the Atacama Large Millimeter/submillimeter Array (ALMA) have revealed its complex, hierarchical structure. In this paper we compare the properties of recent hydrodynamical simulations of the Brick to those of the ALMA observations. To facilitate the comparison, we post-process the simulations and create synthetic ALMA maps of molecular line emission from eight molecules. We correlate the line emission maps to each other and to the mass column density, and find that HNCO is the best mass tracer of the eight emission lines within the simulations. Additionally, we characterise the spatial structure of the observed and simulated cloud using the density probability distribution function (PDF), spatial power spectrum, fractal dimension, and moments of inertia. While we find good agreement between the observed and simulated data in terms of power spectra and fractal dimensions, there are key differences in the density PDFs and moments of inertia, which we attribute to the omission of magnetic fields in the simulations. This demonstrates that the presence of the Galactic potential can reproduce many cloud properties, but additional physical processes are needed to fully explain the gas structure. △ Less

Submitted 1 August, 2023; v1 submitted 19 April, 2021; originally announced April 2021.

Comments: 27 pages; 21 figures; accepted by MNRAS

Journal ref: MNRAS, 520, 2245-2268 (2023)

Galactic 26Al traces metal loss through hot chimneys

Authors: Martin G. H. Krause, Donna Rodgers-Lee, James E. Dale, Roland Diehl, Chiaki Kobayashi

Abstract: Radioactive 26Al is an excellent tracer for metal ejection in the Milky Way, and can provide a direct constraint on the modelling of supernova feedback in galaxy evolution. Gamma-ray observations of the 26Al decay line have found high velocities and hence require a significant fraction of the Galactic 26Al in the hot component. At the same time, meteoritic data combined with simulation results sug… ▽ More Radioactive 26Al is an excellent tracer for metal ejection in the Milky Way, and can provide a direct constraint on the modelling of supernova feedback in galaxy evolution. Gamma-ray observations of the 26Al decay line have found high velocities and hence require a significant fraction of the Galactic 26Al in the hot component. At the same time, meteoritic data combined with simulation results suggest that a significant amount of 26Al makes its way into stars before decay. We investigated the distribution into hot and cold channels with a simulation of a Milky-Way-like galaxy with massive-star feedback in superbubbles and with ejecta traced by 26Al. About 30-40 per cent of the ejecta remain hot, with typical cooling times of the order Gyr. 26Al traces the footpoints of a chimney-fed outflow that mixes metals turbulently into the halo of the model galaxy on a scale of at least 50 kpc. The rest diffuses into cold gas with temperatures less than about 10,000 K, and may therefore be quickly available for star formation. We discuss the robustness of the result by comparison to a simulation with a different global flow pattern. The branching ratio into hot and cold components is comparable to that of longer term average results from chemical evolution modelling of galaxies, clusters and the intracluster medium. △ Less

Submitted 17 November, 2020; originally announced November 2020.

Comments: 10 pages, 10 figures, accepted for MNRAS

Which feedback mechanisms dominate in the high-pressure environment of the Central Molecular Zone?

Authors: Ashley T. Barnes, Steven N. Longmore, James E. Dale, Mark R. Krumholz, J. M. Diederik Kruijssen, Frank Bigiel

Abstract: Supernovae (SNe) dominate the energy and momentum budget of stellar feedback, but the efficiency with which they couple to the interstellar medium (ISM) depends strongly on how effectively early, pre-SN feedback clears dense gas from star-forming regions. There are observational constraints on the magnitudes and timescales of early stellar feedback in low ISM pressure environments, yet no such con… ▽ More Supernovae (SNe) dominate the energy and momentum budget of stellar feedback, but the efficiency with which they couple to the interstellar medium (ISM) depends strongly on how effectively early, pre-SN feedback clears dense gas from star-forming regions. There are observational constraints on the magnitudes and timescales of early stellar feedback in low ISM pressure environments, yet no such constraints exist for more cosmologically typical high ISM pressure environments. In this paper, we determine the mechanisms dominating the expansion of HII regions as a function of size-scale and evolutionary time within the high-pressure ($P/k_\rm{B}$~$10^{7-8}$K cm$^{-3}$) environment in the inner 100pc of the Milky Way. We calculate the thermal pressure from the warm ionised ($P_\rm{HII}$; 10$^{4}$K) gas, direct radiation pressure ($P_\rm{dir}$), and dust processed radiation pressure ($P_\rm{IR}$). We find that (1) $P_\rm{dir}$ dominates the expansion on small scales and at early times (0.01-0.1pc; $<$0.1Myr); (2) the expansion is driven by $P_\rm{HII}$ on large scales at later evolutionary stages ($>0.1$pc; $>1$Myr); (3) during the first ~1Myr of growth, but not thereafter, either $P_{\rm IR}$ or stellar wind pressure likely make a comparable contribution. Despite the high confining pressure of the environment, natal star-forming gas is efficiently cleared to radii of several pc within ~2Myr, i.e. before the first SNe explode. This `pre-processing' means that subsequent SNe will explode into low density gas, so their energy and momentum will efficiently couple to the ISM. We find the HII regions expand to a radius of 3pc, at which point they have internal pressures equal with the surrounding external pressure. A comparison with HII regions in lower pressure environments shows that the maximum size of all HII regions is set by pressure equilibrium with the ambient ISM. △ Less

Submitted 8 September, 2020; originally announced September 2020.

Comments: 20+1 pages, 8 figures, 6 tables. A machine-readable version of Table A1 is available online. Accepted for publication in MNRAS

Collision between Molecular Clouds I. The effect of the cloud virial ratio in head-on collisions

Authors: Tabassum S Tanvir, James E Dale

Abstract: In a series of papers we investigate the effect of collisions between turbulent molecular clouds on their structure, evolution and star formation activity. In this paper we look into the role of the clouds' initial virial ratios. Three different scenarios were examined: both clouds initially bound, one cloud bound and one unbound, and both clouds initially unbound. Models in which one or both clou… ▽ More In a series of papers we investigate the effect of collisions between turbulent molecular clouds on their structure, evolution and star formation activity. In this paper we look into the role of the clouds' initial virial ratios. Three different scenarios were examined: both clouds initially bound, one cloud bound and one unbound, and both clouds initially unbound. Models in which one or both clouds are bound generate filamentary structures aligned along the collision axis and discernible in position-position and position--velocity space. If neither cloud is bound, no filaments result. Unlike in previous simulations of collisions between smooth clouds, owing to the substructure created in the clouds by turbulence before the collisions, dissipation of kinetic energy by the collision is very inefficient and in none of our simulations is sufficient bulk kinetic energy lost to render the clouds bound. Simulations where both clouds are bound created twice as much stellar mass than the bound-unbound model, and both these scenarios produced much more stellar mass than the simulation in which both clouds are unbound. Each simulation was also compared with a control run in which the clouds do not collide. We find the bound-bound collision increases the overall star formation efficiency by a factor of approximately two relative to the control, but that the bound-unbound collision produces a much smaller increase, and the collision has very little effect on the unbound--unbound cloud collision. △ Less

Submitted 12 March, 2020; v1 submitted 10 March, 2020; originally announced March 2020.

Comments: 13 pages, 15 figures; Accepted by MNRAS

Supernova feedback and the energy deposition in molecular clouds

Authors: William E. Lucas, Ian A. Bonnell, James E. Dale

Abstract: Feedback from supernovae is often invoked as an important process in limiting star formation, removing gas from galaxies and hence as a determining process in galaxy formation. Here we report on numerical simulations investigating the interaction between supernova explosions and the natal molecular cloud. We also consider the cases with and without previous feedback from the high-mass star in the… ▽ More Feedback from supernovae is often invoked as an important process in limiting star formation, removing gas from galaxies and hence as a determining process in galaxy formation. Here we report on numerical simulations investigating the interaction between supernova explosions and the natal molecular cloud. We also consider the cases with and without previous feedback from the high-mass star in the form of ionising radiation and stellar winds. The supernova is able to find weak points in the cloud and create channels through which it can escape, leaving much of the well shielded cloud largely unaffected. This effect is increased when the channels are pre-existing due to the effects of previous stellar feedback. The expanding supernova deposits its energy in the gas that is in these exposed channels, and hence sweeps up less mass when feedback has already occurred, resulting in faster outflows with less radiative losses. The full impact of the supernova explosion is then able to impact the larger scale of the galaxy in which it abides. We conclude that supernova explosions only have moderate effects on their dense natal environments but that with pre-existing feedback, the energetic effects of the supernova are able to escape and affect the wider scale medium of the galaxy. △ Less

Submitted 12 February, 2020; originally announced February 2020.

Comments: 12 pages, 10 figures, MNRAS in press

VVV-WIT-01: highly obscured classical nova or protostellar collision?

Authors: P. W. Lucas, D. Minniti, A. Kamble, D. L. Kaplan, N. Cross, I. Dekany, V. D. Ivanov, R. Kurtev, R. K. Saito, L. C. Smith, M. Catelan, N. Masetti, I. Toledo, M. Hempel, M. A. Thompson, C. Contreras Peña, J. Forbrich, M. Krause, J. Dale, J. Borissova, J. Emerson

Abstract: A search of the first Data Release of the VISTA Variables in the Via Lactea (VVV) Survey discovered the exceptionally red transient VVV-WIT-01 ($H-K_s=5.2$). It peaked before March 2010, then faded by $\sim$9.5 mag over the following two years. The 1.6--22 $μ$m spectral energy distribution in March 2010 was well fit by a highly obscured black body with $T \sim 1000$ K and $A_{K_s} \sim 6.6$ mag. T… ▽ More A search of the first Data Release of the VISTA Variables in the Via Lactea (VVV) Survey discovered the exceptionally red transient VVV-WIT-01 ($H-K_s=5.2$). It peaked before March 2010, then faded by $\sim$9.5 mag over the following two years. The 1.6--22 $μ$m spectral energy distribution in March 2010 was well fit by a highly obscured black body with $T \sim 1000$ K and $A_{K_s} \sim 6.6$ mag. The source is projected against the Infrared Dark Cloud (IRDC) SDC G331.062$-$0.294. The chance projection probability is small for any single event ($p \approx 0.01$ to 0.02) which suggests a physical association, e.g. a collision between low mass protostars. However, black body emission at $T \sim 1000$ K is common in classical novae (especially CO novae) at the infrared peak in the light curve, due to condensation of dust $\sim$30--60 days after the explosion. Radio follow up with the Australia Telescope Compact Array (ATCA) detected a fading continuum source with properties consistent with a classical nova but probably inconsistent with colliding protostars. Considering all VVV transients that could have been projected against a catalogued IRDC raises the probability of a chance association to $p=0.13$ to 0.24. After weighing several options, it appears likely that VVV-WIT-01 was a classical nova event located behind an IRDC. △ Less

Submitted 15 January, 2020; originally announced January 2020.

Comments: 13 pages, 7 figures. Accepted by MNRAS

Carina's Pillars of Destruction: the view from ALMA

Authors: P. D. Klaassen, M. R. Reiter, A. F. McLeod, J. C. Mottram, J. E. Dale, M. Gritschneder

Abstract: Forming high-mass stars have a significant effect on their natal environment. Their feedback pathways, including winds, outflows, and ionising radiation, shape the evolution of their surroundings which impacts the formation of the next generation of stars. They create or reveal dense pillars of gas and dust towards the edges of the cavities they clear. They are modelled in feedback simulations, an… ▽ More Forming high-mass stars have a significant effect on their natal environment. Their feedback pathways, including winds, outflows, and ionising radiation, shape the evolution of their surroundings which impacts the formation of the next generation of stars. They create or reveal dense pillars of gas and dust towards the edges of the cavities they clear. They are modelled in feedback simulations, and the sizes and shapes of the pillars produced are consistent with those observed. However, these models predict measurably different kinematics which provides testable discriminants. Here we present the first ALMA Compact Array (ACA) survey of 13 pillars in Carina, observed in $^{12}$CO, $^{13}$CO and C$^{18}$O J=2-1, and the 230 GHz continuum. The pillars in this survey were chosen to cover a wide range in properties relating to the amount and direction of incident radiation, proximity to nearby irradiating clusters and cloud rims, and whether they are detached from the cloud. With these data, we are able to discriminate between models. We generally find pillar velocity dispersions of $<$ 1 km s$^{-1}$ and that the outer few layers of molecular emission in these pillars show no significant offsets from each other, suggesting little bulk internal motions within the pillars. There are instances where the pillars are offset in velocity from their parental cloud rim, and some with no offset, hinting at a stochastic development of these motions. △ Less

Submitted 21 October, 2019; originally announced October 2019.

Comments: 24 Pages, 19 figures. Accepted to MNRAS

Illuminating the Tadpole's metamorphosis I: MUSE observations of a small globule in a sea of ionizing photons

Authors: Megan Reiter, Anna F. McLeod, Pamela D. Klaassen, Andrés E. Guzmán, J. E. Dale, Joseph C. Mottram, Guido Garay

Abstract: We present new MUSE/VLT observations of a small globule in the Carina H II region that hosts the HH 900 jet+outflow system. Data were obtained with the GALACSI ground-layer adaptive optics system in wide-field mode, providing spatially-resolved maps of diagnostic emission lines. These allow us to measure the variation of the physical properties in the globule and jet+outflow system. We find high t… ▽ More We present new MUSE/VLT observations of a small globule in the Carina H II region that hosts the HH 900 jet+outflow system. Data were obtained with the GALACSI ground-layer adaptive optics system in wide-field mode, providing spatially-resolved maps of diagnostic emission lines. These allow us to measure the variation of the physical properties in the globule and jet+outflow system. We find high temperatures ($T_e \approx 10^4$ K), modest extinction ($A_V \approx 2.5$ mag), and modest electron densities ($n_e \approx 200$ cm$^{-3}$) in the ionized gas. Higher excitation lines trace the ionized outflow; both the excitation and ionization in the outflow increase with distance from the opaque globule. In contrast, lower excitation lines that are collisionally de-excited at densities $\gtrsim 10^4$ cm$^{-3}$ trace the highly collimated protostellar jet. Assuming the globule is an isothermal sphere confined by the pressure of the ionization front, we compute a Bonnor-Ebert mass of $\sim 3.7$ M$_{\odot}$. This is two orders of magnitude higher than previous mass estimates, calling into question whether small globules like the Tadpole contribute to the bottom of the IMF. Derived globule properties are consistent with a cloud that has been and/or will be compressed by the ionization front on its surface. At the estimated globule photoevaporation rate of $\sim 5 \times 10^{-7}$ M$_{\odot}$ yr$^{-1}$, the globule will be completely ablated in $\sim 7$ Myr. Stars that form in globules like the Tadpole will emerge into the H II later and may help resolve some of the temporal tension between disk survival and enrichment. △ Less

Submitted 27 September, 2019; originally announced September 2019.

Comments: MNRAS, accepted

Synthetic $^{26}$Al emission from galactic-scale superbubble simulations

Authors: D. Rodgers-Lee, M. G. H. Krause, J. Dale, R. Diehl

Abstract: Emission from the radioactive trace element $^{26}$Al has been observed throughout the Milky Way with the COMPTEL and INTEGRAL satellites. In particular the Doppler shifts measured with INTEGRAL connect $^{26}$Al with superbubbles, which may guide $^{26}$Al flows off spiral arms in the direction of Galactic rotation. In order to test this paradigm, we have performed galaxy-scale simulations of sup… ▽ More Emission from the radioactive trace element $^{26}$Al has been observed throughout the Milky Way with the COMPTEL and INTEGRAL satellites. In particular the Doppler shifts measured with INTEGRAL connect $^{26}$Al with superbubbles, which may guide $^{26}$Al flows off spiral arms in the direction of Galactic rotation. In order to test this paradigm, we have performed galaxy-scale simulations of superbubbles with $^{26}$Al injection in a Milky Way-type galaxy. We produce all-sky synthetic $γ-$ray emission maps of the simulated galaxies. We find that the 1809keV emission from the radioactive decay of $^{26}$Al is highly variable with time and the observer's position. This allows us to estimate an additional systematic variability of 0.2dex for a star formation rate derived from $^{26}$Al for different times and measurement locations in Milky Way-type galaxies. High-latitude morphological features indicate nearby emission with correspondingly high integrated $γ-$ray intensities. We demonstrate that the $^{26}$Al scale height from our simulated galaxies depends on the assumed halo gas density. We present the first synthetic 1809keV longitude-velocity diagrams from 3D hydrodynamic simulations. The line-of-sight velocities for $^{26}$Al can be significantly different from the line-of-sight velocities associated with the cold gas. Over time, $^{26}$Al velocities consistent with the INTEGRAL observations, within uncertainties, appear at any given longitude, broadly supporting previous suggestions that $^{26}$Al injected into expanding superbubbles by massive stars may be responsible for the high velocities found in the INTEGRAL observations. We discuss the effect of systematically varying the location of the superbubbles relative to the spiral arms. △ Less

Submitted 24 September, 2019; originally announced September 2019.

Comments: 20 pages, 12 figures, accepted for publication in MNRAS

The dynamical evolution of molecular clouds near the Galactic Centre -- III. Tidally--induced star formation in protocluster clouds

Authors: James E. Dale, J. M. Diederik Kruijssen, S. N. Longmore

Abstract: As part of a series of papers aimed at understanding the evolution of the Milky Way's Central Molecular Zone (CMZ), we present hydrodynamical simulations of turbulent molecular clouds orbiting in an accurate model of the gravitational potential extant there. We consider two sets of model clouds differing in the energy content of their velocity fields. In the first, self--virialised set, the turbul… ▽ More As part of a series of papers aimed at understanding the evolution of the Milky Way's Central Molecular Zone (CMZ), we present hydrodynamical simulations of turbulent molecular clouds orbiting in an accurate model of the gravitational potential extant there. We consider two sets of model clouds differing in the energy content of their velocity fields. In the first, self--virialised set, the turbulent kinetic energies are chosen to be close in magnitude to the clouds' self--gravitational potential energies. Comparison with isolated clouds evolving without an external potential shows that the self--virialised clouds are unable to withstand the compressive tidal field of the CMZ and rapidly collapse, forming stars much faster and reaching gas exhaustion after a small fraction of a Galactocentric orbit. In the second, tidally--virialised, set of simulations, the clouds' turbulent kinetic energies are in equilibrium with the external tidal field. These models are better supported against the field and the stronger turbulence suppresses star formation. Our results strongly support the inference that anomalously low star formation rates in the CMZ are due primarily to high velocity dispersions in the molecular gas. The clouds follow open, eccentric orbits oscillating in all three spatial coordinates. We examine the consequences of the orbital dynamics, particularly pericentre passage, by performing companion simulations of clouds on circular orbits. The increased tidal forces at pericentre produce transient accelerations in star formation rates of at most a factor of 2.7. Our results demonstrate that modelling star formation in galactic centres requires the inclusion of tidal forces. △ Less

Submitted 25 March, 2019; originally announced March 2019.

Comments: 21 pages, 18 figures, accepted by MNRAS

'The Brick' is not a brick: A comprehensive study of the structure and dynamics of the Central Molecular Zone cloud G0.253+0.016

Authors: J. D. Henshaw, A. Ginsburg, T. J. Haworth, S. N. Longmore, J. M. D. Kruijssen, E. A. C. Mills, V. Sokolov, D. L. Walker, A. T. Barnes, Y. Contreras, J. Bally, C. Battersby, H. Beuther, N. Butterfield, J. E. Dale, T. Henning, J. M. Jackson, J. Kauffmann, T. Pillai, S. Ragan, M. Riener, Q. Zhang

Abstract: In this paper we provide a comprehensive description of the internal dynamics of G0.253+0.016 (a.k.a. 'the Brick'); one of the most massive and dense molecular clouds in the Galaxy to lack signatures of widespread star formation. As a potential host to a future generation of high-mass stars, understanding largely quiescent molecular clouds like G0.253+0.016 is of critical importance. In this paper… ▽ More In this paper we provide a comprehensive description of the internal dynamics of G0.253+0.016 (a.k.a. 'the Brick'); one of the most massive and dense molecular clouds in the Galaxy to lack signatures of widespread star formation. As a potential host to a future generation of high-mass stars, understanding largely quiescent molecular clouds like G0.253+0.016 is of critical importance. In this paper, we reanalyse Atacama Large Millimeter Array cycle 0 HNCO $J=4(0,4)-3(0,3)$ data at 3 mm, using two new pieces of software which we make available to the community. First, scousepy, a Python implementation of the spectral line fitting algorithm scouse. Secondly, acorns (Agglomerative Clustering for ORganising Nested Structures), a hierarchical n-dimensional clustering algorithm designed for use with discrete spectroscopic data. Together, these tools provide an unbiased measurement of the line of sight velocity dispersion in this cloud, $σ_{v_{los}, {\rm 1D}}=4.4\pm2.1$ kms$^{-1}$, which is somewhat larger than predicted by velocity dispersion-size relations for the Central Molecular Zone (CMZ). The dispersion of centroid velocities in the plane of the sky are comparable, yielding $σ_{v_{los}, {\rm 1D}}/σ_{v_{pos}, {\rm 1D}}\sim1.2\pm0.3$. This isotropy may indicate that the line-of-sight extent of the cloud is approximately equivalent to that in the plane of the sky. Combining our kinematic decomposition with radiative transfer modelling we conclude that G0.253+0.016 is not a single, coherent, and centrally-condensed molecular cloud; 'the Brick' is not a \emph{brick}. Instead, G0.253+0.016 is a dynamically complex and hierarchically-structured molecular cloud whose morphology is consistent with the influence of the orbital dynamics and shear in the CMZ. △ Less

Submitted 7 February, 2019; originally announced February 2019.

Comments: 29 pages, 19 figures, 1 table (including appendix). Accepted for publication in MNRAS (February 4, 2019). Scousepy is available here: https://github.com/jdhenshaw/scousepy. Acorns is available here: https://github.com/jdhenshaw/acorns

The dynamical evolution of molecular clouds near the Galactic Centre - II. Spatial structure and kinematics of simulated clouds

Authors: J. M. D. Kruijssen, J. E. Dale, S. N. Longmore, D. L. Walker, J. D. Henshaw, S. M. R. Jeffreson, M. A. Petkova, A. Ginsburg, A. T. Barnes, C. D. Battersby, K. Immer, J. M. Jackson, E. R. Keto, N. Krieger, E. A. C. Mills, Á. Sánchez-Monge, A. Schmiedeke, S. T. Suri, Q. Zhang

Abstract: The evolution of molecular clouds in galactic centres is thought to differ from that in galactic discs due to a significant influence of the external gravitational potential. We present a set of numerical simulations of molecular clouds orbiting on the 100-pc stream of the Central Molecular Zone (the central $\sim500$ pc of the Galaxy) and characterise their morphological and kinematic evolution i… ▽ More The evolution of molecular clouds in galactic centres is thought to differ from that in galactic discs due to a significant influence of the external gravitational potential. We present a set of numerical simulations of molecular clouds orbiting on the 100-pc stream of the Central Molecular Zone (the central $\sim500$ pc of the Galaxy) and characterise their morphological and kinematic evolution in response to the background potential and eccentric orbital motion. We find that the clouds are shaped by strong shear and torques, by tidal and geometric deformation, and by their passage through the orbital pericentre. Within our simulations, these mechanisms control cloud sizes, aspect ratios, position angles, filamentary structure, column densities, velocity dispersions, line-of-sight velocity gradients, spin angular momenta, and kinematic complexity. By comparing these predictions to observations of clouds on the Galactic Centre 'dust ridge', we find that our simulations naturally reproduce a broad range of key observed morphological and kinematic features, which can be explained in terms of well-understood physical mechanisms. We argue that the accretion of gas clouds onto the central regions of galaxies, where the rotation curve turns over and the tidal field is fully compressive, is accompanied by transformative dynamical changes to the clouds, leading to collapse and star formation. This can generate an evolutionary progression of cloud collapse with a common starting point, which either marks the time of accretion onto the tidally-compressive region or of the most recent pericentre passage. Together, these processes may naturally produce the synchronised starbursts observed in numerous (extra)galactic nuclei. △ Less

Submitted 7 February, 2019; v1 submitted 5 February, 2019; originally announced February 2019.

Comments: 21 pages, 8 figures, 3 tables; accepted by MNRAS (February 5, 2019); Figures 2, 3, and 4 show the main results of the paper. An animated version of Figure 2 is available as an ancillary file in the arXiv source and will be included as online Supporting Information with the MNRAS publication

Feedback from massive stars at low metallicities: MUSE observations of N44 and N180 in the Large Magellanic Cloud

Authors: A. F. McLeod, J. E. Dale, C. J. Evans, A. Ginsburg, J. M. D. Kruijssen, E. W. Pellegrini, S. K. Ramsay, L. Testi

Abstract: We present MUSE integral field data of two HII region complexes in the Large Magellanic Cloud (LMC), N44 and N180. Both regions consist of a main superbubble and a number of smaller, more compact HII regions that formed on the edge of the superbubble. For a total of 11 HII regions, we systematically analyse the radiative and mechanical feedback from the massive O-type stars on the surrounding gas.… ▽ More We present MUSE integral field data of two HII region complexes in the Large Magellanic Cloud (LMC), N44 and N180. Both regions consist of a main superbubble and a number of smaller, more compact HII regions that formed on the edge of the superbubble. For a total of 11 HII regions, we systematically analyse the radiative and mechanical feedback from the massive O-type stars on the surrounding gas. We exploit the integral field property of the data and the coverage of the HeII$λ$5412 line to identify and classify the feedback-driving massive stars, and from the estimated spectral types and luminosity classes we determine the stellar radiative output in terms of the ionising photon flux $Q_{0}$. We characterise the HII regions in terms of their sizes, morphologies, ionisation structure, luminosity and kinematics, and derive oxygen abundances via emission line ratios. We analyse the role of different stellar feedback mechanisms for each region by measuring the direct radiation pressure, the pressure of the ionised gas, and the pressure of the shock-heated winds. We find that stellar winds and ionised gas are the main drivers of HII region expansion in our sample, while the direct radiation pressure is up to three orders of magnitude lower than the other terms. We relate the total pressure to the star formation rate per unit area, $Σ_{SFR}$, for each region and find that stellar feedback has a negative effect on star formation, and sets an upper limit to $Σ_{SFR}$ as a function of increasing pressure. △ Less

Submitted 2 October, 2018; originally announced October 2018.

Comments: Accepted for publication in MNRAS, 27 pages, 21 figures

Synthetic observations of star formation and the interstellar medium

Authors: Thomas J. Haworth, Simon C. O. Glover, Christine M. Koepferl, Thomas G. Bisbas, James E. Dale

Abstract: Synthetic observations are playing an increasingly important role across astrophysics, both for interpreting real observations and also for making meaningful predictions from models. In this review, we provide an overview of methods and tools used for generating, manipulating and analysing synthetic observations and their application to problems involving star formation and the interstellar medium… ▽ More Synthetic observations are playing an increasingly important role across astrophysics, both for interpreting real observations and also for making meaningful predictions from models. In this review, we provide an overview of methods and tools used for generating, manipulating and analysing synthetic observations and their application to problems involving star formation and the interstellar medium. We also discuss some possible directions for future research using synthetic observations. △ Less

Submitted 9 November, 2018; v1 submitted 14 November, 2017; originally announced November 2017.

Comments: 71 Pages, 16 figures. Revised following reviewers comments

No preferential spatial distribution for massive stars expected from their formation

Authors: Richard J. Parker, James E. Dale

Abstract: We analyse N-body and Smoothed Particle Hydrodynamic (SPH) simulations of young star-forming regions to search for differences in the spatial distributions of massive stars compared to lower-mass stars. The competitive accretion theory of massive star formation posits that the most massive stars should sit in deeper potential wells than lower-mass stars. This may be observable in the relative surf… ▽ More We analyse N-body and Smoothed Particle Hydrodynamic (SPH) simulations of young star-forming regions to search for differences in the spatial distributions of massive stars compared to lower-mass stars. The competitive accretion theory of massive star formation posits that the most massive stars should sit in deeper potential wells than lower-mass stars. This may be observable in the relative surface density or spatial concentration of the most massive stars compared to other, lower-mass stars. Massive stars in cool--collapse N-body models do end up in significantly deeper potentials, and are mass segregated. However, in models of warm (expanding) star-forming regions, whilst the massive stars do come to be in deeper potentials than average stars, they are not mass segregated. In the purely hydrodynamical SPH simulations, the massive stars do come to reside in deeper potentials, which is due to their runaway growth. However, when photoionisation and stellar winds are implemented in the simulations, these feedback mechanisms regulate the mass of the stars and disrupt the inflow of gas into the clouds' potential wells. This generally makes the potential wells shallower than in the control runs, and prevents the massive stars from occupying deeper potentials. This in turn results in the most massive stars having a very similar spatial concentration and surface density distribution to lower-mass stars. Whilst massive stars do form via competitive accretion in our simulations, this rarely translates to a different spatial distribution and so any lack of primordial mass segregation in an observed star-forming region does not preclude competitive accretion as a viable formation mechanism for massive stars. △ Less

Submitted 12 May, 2017; originally announced May 2017.

Comments: 12 pages, 6 figures, accepted for publication in MNRAS (12.05.2017)

Cloud structure of three Galactic infrared dark star-forming regions from combining ground and space based bolometric observations

Authors: Yuxin Lin, Hauyu Baobab Liu, James E. Dale, Di Li, Gemma Busquet, Zhi-Yu Zhang, Adam Ginsburg, Roberto Galvan-Madrid, Attila Kovacs, Eric Koch, Lei Qian, Ke Wang, Steve Longmore, Huei-Ru Chen, Daniel Walker

Abstract: We have modified the iterative procedure introduced by Lin et al. (2016), to systematically combine the submm images taken from ground based (e.g., CSO, JCMT, APEX) and space (e.g., Herschel, Planck) telescopes. We applied the updated procedure to observations of three well studied Infrared Dark Clouds (IRDCs): G11.11-0.12, G14.225-0.506 and G28.34+0.06, and then performed single-component, modifi… ▽ More We have modified the iterative procedure introduced by Lin et al. (2016), to systematically combine the submm images taken from ground based (e.g., CSO, JCMT, APEX) and space (e.g., Herschel, Planck) telescopes. We applied the updated procedure to observations of three well studied Infrared Dark Clouds (IRDCs): G11.11-0.12, G14.225-0.506 and G28.34+0.06, and then performed single-component, modified black-body fits to derive $\sim$10$"$ resolution dust temperature and column density maps. The derived column density maps show that these three IRDCs exhibit complex filamentary structures embedding with rich clumps/cores. We compared the column density probability distribution functions (N-PDFs) and two-point correlation (2PT) functions of the column density field between these IRDCs with several OB cluster-forming regions. Based on the observed correlation and measurements, and complementary hydrodynamical simulations for a 10$^{4}$ $\rm M_{\odot}$ molecular cloud, we hypothesize that cloud evolution can be better characterized by the evolution of the (column) density distribution function and the relative power of dense structures as a function of spatial scales, rather than merely based on the presence of star-forming activity. Based on the small analyzed sample, we propose four evolutionary stages, namely: {\it cloud integration, stellar assembly, cloud pre-dispersal and dispersed-cloud.} The initial {\it cloud integration} stage and the final {\it dispersed cloud} stage may be distinguished from the two intermediate stages by a steeper than $-$4 power-law index of the N-PDF. The {\it cloud integration} stage and the subsequent {\it stellar assembly} stage are further distinguished from each other by the larger luminosity-to-mass ratio ($>$40 $\rm L_{\odot}/M_{\odot}$) of the latter. △ Less

Submitted 21 April, 2017; originally announced April 2017.

Comments: 14 pages, 13 figures; accepted for publication in ApJ

Thermal Feedback in the high-mass star and cluster forming region W51

Authors: Adam Ginsburg, Ciriaco Goddi, J. M. Diederik Kruijssen, John Bally, Rowan Smith, Roberto Galván-Madrid, Elisabeth A. C. Mills, Ke Wang, James E. Dale, Jeremy Darling, Erik Rosolowsky, Robert Loughnane, Leonardo Testi, Nate Bastian

Abstract: High-mass stars have generally been assumed to accrete most of their mass while already contracted onto the main sequence, but this hypothesis has not been observationally tested. We present ALMA observations of a 3 x 1.5 pc area in the W51 high-mass star-forming complex. We identify dust continuum sources and measure the gas and dust temperature through both rotational diagram modeling of CH3OH a… ▽ More High-mass stars have generally been assumed to accrete most of their mass while already contracted onto the main sequence, but this hypothesis has not been observationally tested. We present ALMA observations of a 3 x 1.5 pc area in the W51 high-mass star-forming complex. We identify dust continuum sources and measure the gas and dust temperature through both rotational diagram modeling of CH3OH and brightness-temperature-based limits. The observed region contains three high-mass YSOs that appear to be at the earliest stages of their formation, with no signs of ionizing radiation from their central sources. The data reveal high gas and dust temperatures (T > 100 K) extending out to about 5000 AU from each of these sources. There are no clear signs of disks or rotating structures down to our 1000 AU resolution. The extended warm gas provides evidence that, during the process of forming, these high-mass stars heat a large volume and correspondingly large mass of gas in their surroundings, inhibiting fragmentation and therefore keeping a large reservoir available to feed from. By contrast, the more mature massive stars that illuminate compact H ii regions have little effect on their surrounding dense gas, suggesting that these main sequence stars have completed most or all of their accretion. The high luminosity of the massive protostars (L > 10^4 solar luminosities), combined with a lack of centimeter continuum emission from these sources, implies that they are not on the main sequence while they accrete the majority of their mass; instead, they may be bloated and cool. △ Less

Submitted 17 May, 2017; v1 submitted 5 April, 2017; originally announced April 2017.

Comments: Accepted in ApJ. Data are publicly available at http://dx.doi.org/10.7910/DVN/8QJT3K . Repository of analysis scripts is at https://github.com/adamginsburg/W51_ALMA_2013.1.00308.S

The effect of the virial state of molecular clouds on the influence of feedback from massive stars

Authors: James E. Dale

Abstract: A set of Smoothed Particle Hydrodynamics simulations of the influence of photoionising radiation and stellar winds on a series of 10$^{4}$M$_{\odot}$ turbulent molecular clouds with initial virial ratios of 0.7, 1.1, 1.5, 1.9 and 2.3 and initial mean densities of 136, 1135 and 9096\,cm$^{-3}$ are presented. Reductions in star formation efficiency rates are found to be modest, in the range… ▽ More A set of Smoothed Particle Hydrodynamics simulations of the influence of photoionising radiation and stellar winds on a series of 10$^{4}$M$_{\odot}$ turbulent molecular clouds with initial virial ratios of 0.7, 1.1, 1.5, 1.9 and 2.3 and initial mean densities of 136, 1135 and 9096\,cm$^{-3}$ are presented. Reductions in star formation efficiency rates are found to be modest, in the range $30\%-50\%$ and to not vary greatly across the parameter space. In no case was star formation entirely terminated over the $\approx3$\,Myr duration of the simulations. The fractions of material unbound by feedback are in the range $20-60\%$, clouds with the lowest escape velocities being the most strongly affected.\\ Leakage of ionised gas leads to the HII regions rapidly becoming underpressured. The destructive effects of ionisation are thus largely not due to thermally--driven expansion of the HII regions, but to momentum transfer by photoevaporation of fresh material. Our simulations have similar global ionisation rates and we show that the effects of feedback upon them can be adequately modelled as a steady injection of momentum, resembling a momentum--conserving wind. △ Less

Submitted 13 January, 2017; originally announced January 2017.

Comments: 18 pages, 14 figures (mostly mangled to get under size limit), accepted by MNRAS

Connecting the dots: a correlation between ionising radiation and cloud mass-loss rate traced by optical integral field spectroscopy

Authors: A. F. McLeod, M. Gritschneder, J. E. Dale, A. Ginsburg, P. D. Klaassen, J. C. Mottram, T. Preibisch, S. Ramsay, M. Reiter, L. Testi

Abstract: We present an analysis of the effect of feedback from O- and B-type stars with data from the integral field spectrograph MUSE mounted on the Very Large Telescope of pillar-like structures in the Carina Nebular Complex, one of the most massive star-forming regions in the Galaxy. For the observed pillars, we compute gas electron densities and temperatures maps, produce integrated line and velocity m… ▽ More We present an analysis of the effect of feedback from O- and B-type stars with data from the integral field spectrograph MUSE mounted on the Very Large Telescope of pillar-like structures in the Carina Nebular Complex, one of the most massive star-forming regions in the Galaxy. For the observed pillars, we compute gas electron densities and temperatures maps, produce integrated line and velocity maps of the ionised gas, study the ionisation fronts at the pillar tips, analyse the properties of the single regions, and detect two ionised jets originating from two distinct pillar tips. For each pillar tip we determine the incident ionising photon flux $Q_\mathrm{0,pil}$ originating from the nearby massive O- and B-type stars and compute the mass-loss rate $\dot{M}$ of the pillar tips due to photo-evaporation caused by the incident ionising radiation. We combine the results of the Carina data set with archival MUSE data of a pillar in NGC 3603 and with previously published MUSE data of the Pillars of Creation in M16, and with a total of 10 analysed pillars, find tight correlations between the ionising photon flux and the electron density, the electron density and the distance from the ionising sources, and the ionising photon flux and the mass-loss rate. The combined MUSE data sets of pillars in regions with different physical conditions and stellar content therefore yield an empirical quantification of the feedback effects of ionising radiation. In agreement with models, we find that $\dot{M}\propto Q_\mathrm{0,pil}^{1/2}$. △ Less

Submitted 29 July, 2016; originally announced August 2016.

Comments: Accepted for publication in MNRAS, 36 pages, 30 figures

Insights from Synthetic Star-forming Regions: III. Calibration of Measurement Techniques of Star-formation Rates

Authors: Christine M. Koepferl, Thomas P. Robitaille, James E. Dale

Abstract: Through an extensive set of realistic synthetic observations (produced in Paper I), we assess in this part of the paper series (Paper III) how the choice of observational techniques affects the measurement of star-formation rates (SFRs) in star-forming regions. We test the accuracy of commonly used techniques and construct new methods to extract the SFR, so that these findings can be applied to me… ▽ More Through an extensive set of realistic synthetic observations (produced in Paper I), we assess in this part of the paper series (Paper III) how the choice of observational techniques affects the measurement of star-formation rates (SFRs) in star-forming regions. We test the accuracy of commonly used techniques and construct new methods to extract the SFR, so that these findings can be applied to measure the SFR in real regions throughout the Milky Way. We investigate diffuse infrared SFR tracers such as those using 24 μm, 70 μm and total infrared emission, which have been previously calibrated for global galaxy scales. We set up a toy model of a galaxy and show that the infrared emission is consistent with the intrinsic SFR using extra-galactic calibrated laws (although the consistency does not prove their reliability). For local scales, we show that these techniques produce completely unreliable results for single star-forming regions, which are governed by different characteristic timescales. We show how calibration of these techniques can be improved for single star-forming regions by adjusting the characteristic timescale and the scaling factor and give suggestions of new calibrations of the diffuse star-formation tracers. We show that star-forming regions that are dominated by high-mass stellar feedback experience a rapid drop in infrared emission once high-mass stellar feedback is turned on, which implies different characteristic timescales. Moreover, we explore the measured SFRs calculated directly from the observed young stellar population. We find that the measured point sources follow the evolutionary pace of star formation more directly than diffuse star-formation tracers. △ Less

Submitted 28 June, 2016; originally announced June 2016.

Comments: Accepted by the Astrophysical Journal (ApJ), online material

Insights from Synthetic Star-forming Regions: II. Verifying Dust Surface Density, Dust Temperature & Gas Mass Measurements with Modified Blackbody Fitting

Authors: Christine M. Koepferl, Thomas P. Robitaille, James E. Dale

Abstract: We use a large data-set of realistic synthetic observations (PaperI) to assess how observational techniques affect the measurement of physical properties of star-forming regions. In this paper (PaperII), we explore the reliability of the measured total gas mass, dust surface density and dust temperature maps derived from modified blackbody fitting of synthetic Herschel observations. We found from… ▽ More We use a large data-set of realistic synthetic observations (PaperI) to assess how observational techniques affect the measurement of physical properties of star-forming regions. In this paper (PaperII), we explore the reliability of the measured total gas mass, dust surface density and dust temperature maps derived from modified blackbody fitting of synthetic Herschel observations. We found from our pixel-by-pixel analysis of the measured dust surface density and dust temperature a worrisome error spread especially close to star-formation sites and low-density regions, where for those "contaminated" pixels the surface densities can be under/overestimated by up to three orders of magnitude. In light of this, we recommend to treat the pixel-based results from this technique with caution in regions with active star formation. In regions of high background typical in the inner Galactic plane, we are not able to recover reliable surface density maps of individual synthetic regions, since low-mass regions are lost in the FIR background. When measuring the total gas mass of regions in moderate background, we find that modified blackbody fitting works well (absolute error:+9%;-13%) up to 10kpc distance (errors increase with distance). Commonly, the initial images are convolved to the largest common beam-size, which smears contaminated pixels over large areas. The resulting information loss makes this commonly-used technique less verifiable as now chi^2-values cannot be used as a quality indicator of a fitted pixel. Our control measurements of the total gas mass (without the step of convolution to the largest common beam size) produce similar results (absolute error:+20%;-7%) while having much lower median errors especially for the high-mass stellar feedback phase. In upcoming papers (III&IV) we test the reliability of measured star-formation rate with direct and indirect techniques. △ Less

Submitted 27 June, 2016; originally announced June 2016.

Comments: Accepted by the Astrophysical Journal (ApJ), 7 page plots, 12 figures, extensive online material

Cloud Structure of Galactic OB Cluster Forming Regions from Combining Ground and Space Based Bolometric Observations

Authors: Yuxin Lin, Hauyu Baobab Liu, Di Li, Zhiyu Zhang, Adam Ginsburg, Jaime E. Pineda, Lei Qian, Roberto Galván-Madrid, Anna Faye McLeod, Erik Rosolowsky, James E. Dale, Katharina Immer, Eric Koch, Steve Longmore, Daniel Walker, Leonardo Testi

Abstract: We have developed an iterative procedure to systematically combine the millimeter and submillimeter images of OB cluster-forming molecular clouds, which were taken by ground based (CSO, JCMT, APEX, IRAM-30m) and space telescopes (Herschel, Planck). For the seven luminous ($L$$>$10$^{6}$ $L_{\odot}$) Galactic OB cluster-forming molecular clouds selected for our analyses, namely W49A, W43-Main, W43-… ▽ More We have developed an iterative procedure to systematically combine the millimeter and submillimeter images of OB cluster-forming molecular clouds, which were taken by ground based (CSO, JCMT, APEX, IRAM-30m) and space telescopes (Herschel, Planck). For the seven luminous ($L$$>$10$^{6}$ $L_{\odot}$) Galactic OB cluster-forming molecular clouds selected for our analyses, namely W49A, W43-Main, W43-South, W33, G10.6-0.4, G10.2-0.3, G10.3-0.1, we have performed single-component, modified black-body fits to each pixel of the combined (sub)millimeter images, and the Herschel PACS and SPIRE images at shorter wavelengths. The $\sim$10$"$ resolution dust column density and temperature maps of these sources revealed dramatically different morphologies, indicating very different modes of OB cluster-formation, or parent molecular cloud structures in different evolutionary stages. The molecular clouds W49A, W33, and G10.6-0.4 show centrally concentrated massive molecular clumps that are connected with approximately radially orientated molecular gas filaments. The W43-Main and W43-South molecular cloud complexes, which are located at the intersection of the Galactic near 3-kpc (or Scutum) arm and the Galactic bar, show a widely scattered distribution of dense molecular clumps/cores over the observed $\sim$10 pc spatial scale. The relatively evolved sources G10.2-0.3 and G10.3-0.1 appear to be affected by stellar feedback, and show a complicated cloud morphology embedded with abundant dense molecular clumps/cores. We find that with the high angular resolution we achieved, our visual classification of cloud morphology can be linked to the systematically derived statistical quantities (i.e., the enclosed mass profile, the column density probability distribution function, the two-point correlation function of column density, and the probability distribution function of clump/core separations). △ Less

Submitted 26 June, 2016; v1 submitted 24 June, 2016; originally announced June 2016.

Comments: ApJ accepted, Some figures in appendix are omitted

Toward gas exhaustion in the W51 high-mass protoclusters

Authors: Adam Ginsburg, W. M. Goss, Ciriaco Goddi, Roberto Galván-Madrid, James E. Dale, John Bally, Cara D. Battersby, Allison Youngblood, Ravi Sankrit, Rowan Smith, Jeremy Darling, J. M. Diederik Kruijssen, Hauyu Baobab Liu

Abstract: We present new JVLA observations of the high-mass cluster-forming region W51A from 2 to 16 GHz with resolution $θ_{fwhm} \approx$ 0.3 - 0.5". The data reveal a wealth of observational results: (1) Currently-forming, very massive (proto-O) stars are traced by o-H2CO $2_{1,1}-2_{1,2}$ emission, suggesting that this line can be used efficiently as a massive protostar tracer. (2) There is a spatially… ▽ More We present new JVLA observations of the high-mass cluster-forming region W51A from 2 to 16 GHz with resolution $θ_{fwhm} \approx$ 0.3 - 0.5". The data reveal a wealth of observational results: (1) Currently-forming, very massive (proto-O) stars are traced by o-H2CO $2_{1,1}-2_{1,2}$ emission, suggesting that this line can be used efficiently as a massive protostar tracer. (2) There is a spatially distributed population of $\sim$mJy continuum sources, including hypercompact H ii regions and candidate colliding wind binaries, in and around the W51 proto-clusters. (3) There are two clearly detected protoclusters, W51e and W51 IRS2, that are gas-rich but may have most of their mass in stars within their inner $\sim$ 0.05 pc. The majority of the bolometric luminosity in W51 most likely comes from a third population of OB stars between these clusters. The presence of a substantial population of exposed O-stars coincident with a population of still-forming massive stars, along with a direct measurement of the low mass loss rate via ionized gas outflow from W51 IRS2, together imply that feedback is ineffective at halting star formation in massive protoclusters. Instead, feedback may shut off the large-scale accretion of diffuse gas onto the W51 protoclusters, implying that they are evolving towards a state of gas exhaustion rather than gas expulsion. Recent theoretical models predict gas exhaustion to be a necessary step in the formation of gravitationally bound stellar clusters, and our results provide an observational validation of this process. △ Less

Submitted 19 June, 2016; v1 submitted 30 May, 2016; originally announced May 2016.

Comments: Accepted to A&A. Associated code at https://github.com/adamginsburg/paper_w51_evla

Journal ref: A&A 595, A27 (2016)

Insights from Synthetic Star-forming Regions: I. Reliable Mock Observations from SPH Simulations

Authors: Christine M. Koepferl, Thomas P. Robitaille, James E. Dale, Francesco Biscani

Abstract: Through synthetic observations of a hydrodynamical simulation of an evolving star-forming region, we assess how the choice of observational techniques affects the measurements of properties which trace star formation. Testing and calibrating observational measurements requires synthetic observations which are as realistic as possible. In this part of the paper series (Paper I), we explore differen… ▽ More Through synthetic observations of a hydrodynamical simulation of an evolving star-forming region, we assess how the choice of observational techniques affects the measurements of properties which trace star formation. Testing and calibrating observational measurements requires synthetic observations which are as realistic as possible. In this part of the paper series (Paper I), we explore different techniques for how to map the distributions of densities and temperatures from the particle-based simulations onto a Voronoi mesh suitable for radiative transfer and consequently explore their accuracy. We further test different ways to set up the radiative transfer in order to produce realistic synthetic observations. We give a detailed description of all methods and ultimately recommend techniques. We have found that the flux around 20 microns is strongly overestimated when blindly coupling the dust radiative transfer temperature with the hydrodynamical gas temperature. We find that when instead assuming a constant background dust temperature in addition to the radiative transfer heating, the recovered flux is consistent with actual observations. We present around 5800 realistic synthetic observations for Spitzer and Herschel bands, at different evolutionary time-steps, distances and orientations. In the upcoming papers of this series (Paper II, Paper III and Paper IV), we will test and calibrate measurements of the star-formation rate (SFR), gas mass and the star-formation efficiency (SFE) using our realistic synthetic observations. △ Less

Submitted 7 March, 2016; originally announced March 2016.

Comments: 17 figures; online data, 5800 synthetic observations (open source), algorithms to create realistic synthetic observations

Molecular gas kinematics within the central 250 pc of the Milky Way

Authors: J. D. Henshaw, S. N. Longmore, J. M. D. Kruijssen, B. Davies, J. Bally, A. Barnes, C. Battersby, M. Burton, M. R. Cunningham, J. E. Dale, A. Ginsburg, K. Immer, P. A. Jones, S. Kendrew, E. A. C. Mills, S. Molinari, T. J. T. Moore, J. Ott, T. Pillai, J. Rathborne, P. Schilke, A. Schmiedeke, L. Testi, D. Walker, A. Walsh , et al. (1 additional authors not shown)

Abstract: Using spectral-line observations of HNCO, N2H+, and HNC, we investigate the kinematics of dense gas in the central ~250 pc of the Galaxy. We present SCOUSE (Semi-automated multi-COmponent Universal Spectral-line fitting Engine), a line fitting algorithm designed to analyse large volumes of spectral-line data efficiently and systematically. Unlike techniques which do not account for complex line pr… ▽ More Using spectral-line observations of HNCO, N2H+, and HNC, we investigate the kinematics of dense gas in the central ~250 pc of the Galaxy. We present SCOUSE (Semi-automated multi-COmponent Universal Spectral-line fitting Engine), a line fitting algorithm designed to analyse large volumes of spectral-line data efficiently and systematically. Unlike techniques which do not account for complex line profiles, SCOUSE accurately describes the {l, b, v_LSR} distribution of CMZ gas, which is asymmetric about Sgr A* in both position and velocity. Velocity dispersions range from 2.6 km/s<σ<53.1 km/s. A median dispersion of 9.8 km/s, translates to a Mach number, M_3D>28. The gas is distributed throughout several "streams", with projected lengths ~100-250 pc. We link the streams to individual clouds and sub-regions, including Sgr C, the 20 and 50 km/s clouds, the dust ridge, and Sgr B2. Shell-like emission features can be explained by the projection of independent molecular clouds in Sgr C and the newly identified conical profile of Sgr B2 in {l ,b, v_LSR} space. These features have previously invoked supernova-driven shells and cloud-cloud collisions as explanations. We instead caution against structure identification in velocity-integrated emission maps. Three geometries describing the 3-D structure of the CMZ are investigated: i) two spiral arms; ii) a closed elliptical orbit; iii) an open stream. While two spiral arms and an open stream qualitatively reproduce the gas distribution, the most recent parameterisation of the closed elliptical orbit does not. Finally, we discuss how proper motion measurements of masers can distinguish between these geometries, and suggest that this effort should be focused on the 20 km/s and 50 km/s clouds and Sgr C. △ Less

Submitted 10 February, 2016; v1 submitted 14 January, 2016; originally announced January 2016.

Comments: Accepted for publication in MNRAS. 30 pages (including appendices), 22 figures, 3 tables. This updated version includes minor corrections (including typos and updated references). SCOUSE is available for download at https://github.com/jdhenshaw/SCOUSE. A selection of PPV movies can be downloaded at https://github.com/jdhenshaw/PPV_movies

Using young massive star clusters to understand star formation and feedback in high-redshift-like environments

Authors: Steven Longmore, Ashley Barnes, Cara Battersby, John Bally, J. M. Diederik Kruijssen, James Dale, Jonathan Henshaw, Daniel Walker, Jill Rathborne, Leonardo Testi, Juergen Ott, Adam Ginsburg

Abstract: The formation environment of stars in massive stellar clusters is similar to the environment of stars forming in galaxies at a redshift of 1 - 3, at the peak star formation rate density of the Universe. As massive clusters are still forming at the present day at a fraction of the distance to high-redshift galaxies they offer an opportunity to understand the processes controlling star formation and… ▽ More The formation environment of stars in massive stellar clusters is similar to the environment of stars forming in galaxies at a redshift of 1 - 3, at the peak star formation rate density of the Universe. As massive clusters are still forming at the present day at a fraction of the distance to high-redshift galaxies they offer an opportunity to understand the processes controlling star formation and feedback in conditions similar to those in which most stars in the Universe formed. Here we describe a system of massive clusters and their progenitor gas clouds in the centre of the Milky Way, and outline how detailed observations of this system may be able to: (i) help answer some of the fundamental open questions in star formation and (ii) quantify how stellar feedback couples to the surrounding interstellar medium in this high-pressure, high-redshift analogue environment. △ Less

Submitted 11 January, 2016; originally announced January 2016.

Comments: 6 pages, 2 figures, to appear in proceedings of The 6th Zermatt ISM Symposium: Conditions and Impact of Star Formation From Lab to Space, eds. R. Simon, M. Rollig

Did the Solar System form in a sequential triggered star formation event?

Authors: Richard J. Parker, James E. Dale

Abstract: The presence and abundance of the short-lived radioisotopes (SLRs) $^{26}$Al and $^{60}$Fe during the formation of the Solar System is difficult to explain unless the Sun formed in the vicinity of one or more massive star(s) that exploded as supernovae. Two different scenarios have been proposed to explain the delivery of SLRs to the protosolar nebula: (i) direct pollution of the protosolar disc b… ▽ More The presence and abundance of the short-lived radioisotopes (SLRs) $^{26}$Al and $^{60}$Fe during the formation of the Solar System is difficult to explain unless the Sun formed in the vicinity of one or more massive star(s) that exploded as supernovae. Two different scenarios have been proposed to explain the delivery of SLRs to the protosolar nebula: (i) direct pollution of the protosolar disc by supernova ejecta and (ii) the formation of the Sun in a sequential star formation event in which supernovae shockwaves trigger further star formation which is enriched in SLRs. The sequentially triggered model has been suggested as being more astrophysically likely than the direct pollution scenario. In this paper we investigate this claim by analysing a combination of $N$-body and SPH simulations of star formation. We find that sequential star formation would result in large age spreads (or even bi-modal age distributions for spatially coincident events) due to the dynamical relaxation of the first star-formation event(s). Secondly, we discuss the probability of triggering spatially and temporally discrete populations of stars and find this to be only possible in very contrived situations. Taken together, these results suggest that the formation of the Solar System in a triggered star formation event is as improbable, if not more so, than the direct pollution of the protosolar disc by a supernova. △ Less

Submitted 20 November, 2015; originally announced November 2015.

Comments: 8 pages, 3 figures, accepted for publication in MNRAS

Observing gas and dust in simulations of star formation with Monte Carlo radiation transport on Voronoi meshes

Authors: D. A. Hubber, B. Ercolano, J. Dale

Abstract: Ionising feedback from massive stars dramatically affects the interstellar medium local to star forming regions. Numerical simulations are now starting to include enough complexity to produce morphologies and gas properties that are not too dissimilar from observations. The comparison between the density fields produced by hydrodynamical simulations and observations at given wavelengths relies how… ▽ More Ionising feedback from massive stars dramatically affects the interstellar medium local to star forming regions. Numerical simulations are now starting to include enough complexity to produce morphologies and gas properties that are not too dissimilar from observations. The comparison between the density fields produced by hydrodynamical simulations and observations at given wavelengths relies however on photoionisation/chemistry and radiative transfer calculations. We present here an implementation of Monte Carlo radiation transport through a Voronoi tessellation in the photoionisation and dust radiative transfer code MOCASSIN. We show for the first time a synthetic spectrum and synthetic emission line maps of an hydrodynamical simulation of a molecular cloud affected by massive stellar feedback. We show that the approach on which previous work is based, which remapped hydrodynamical density fields onto Cartesian grids before performing radiative transfer/photoionisation calculations, results in significant errors in the temperature and ionisation structure of the region. Furthermore, we describe the mathematical process of tracing photon energy packets through a Voronoi tessellation, including optimisations, treating problematic cases and boundary conditions. We perform various benchmarks using both the original version of MOCASSIN and the modified version using the Voronoi tessellation. We show that for uniform grids, or equivalently a cubic lattice of cell generating points, the new Voronoi version gives the same results as the original Cartesian-grid version of MOCASSIN for all benchmarks. For non-uniform initial conditions, such as using snapshots from Smoothed Particle Hydrodynamics simulations, we show that the Voronoi version performs better than the Cartesian grid version, resulting in much better resolution in dense regions. △ Less

Submitted 16 November, 2015; originally announced November 2015.

Comments: 12 pages, 4 figures, 3 tables. Accepted for publication by MNRAS

A nebular analysis of the central Orion Nebula with MUSE

Authors: A. F. Mc Leod, P. M. Weilbacher, A. Ginsburg, J. E. Dale, S. Ramsay, L. Testi

Abstract: A nebular analysis of the central Orion Nebula and its main structures is presented. We exploit MUSE integral field observations in the wavelength range 4595-9366 Å to produce the first O, S and N ionic and total abundance maps of a region spanning 6' x 5' with a spatial resolution of 0.2". We use the S$_{23}$ ( = ([SII]$λ$6717,31+[SIII]$λ$9068)/H$β$) parameter, together with [OII]/[OIII] as an in… ▽ More A nebular analysis of the central Orion Nebula and its main structures is presented. We exploit MUSE integral field observations in the wavelength range 4595-9366 Å to produce the first O, S and N ionic and total abundance maps of a region spanning 6' x 5' with a spatial resolution of 0.2". We use the S$_{23}$ ( = ([SII]$λ$6717,31+[SIII]$λ$9068)/H$β$) parameter, together with [OII]/[OIII] as an indicator of the degree of ionisation, to distinguish between the various small-scale structures. The only Orion Bullet covered by MUSE is HH 201, which shows a double component in the [FeII]$λ$8617 line throughout indicating an expansion, and we discuss a scenario in which this object is undergoing a disruptive event. We separate the proplyds located south of the Bright Bar into four categories depending on their S$_{23}$ values, propose the utility of the S$_{23}$ parameter as an indicator of the shock-contribution to the excitation of line-emitting atoms, and show that the MUSE data is able to identify the proplyds associated with disks and microjets. We compute the second order structure function for the H$α$, [OIII]$λ$5007, [SII]$λ$6731 and [OI]$λ$6300 emission lines to analyse the turbulent velocity field of the region covered with MUSE. We find that the spectral and spatial resolution of MUSE is not able to faithfully reproduce the structure functions of previous works. △ Less

Submitted 9 December, 2015; v1 submitted 5 November, 2015; originally announced November 2015.

Comments: 34 pages, 31 figures, published in MNRAS, added missing citation

Journal ref: MNRAS, 455, 4057 (2016)

A network of filaments detected by Herschel in the Serpens Core: A laboratory to test simulations of low-mass star formation

Authors: V. Roccatagliata, J. E. Dale, T. Ratzka, L. Testi, A. Burkert, C. Koepferl, A. Sicilia-Aguilar, C. Eiroa, B. Gaczkowski

Abstract: Filaments represent a key structure during the early stages of the star formation process. Simulations show filamentary structure commonly formed before and during the formation of cores. Aims. The Serpens Core represents an ideal laboratory to test the state-of-the-art of simulations of turbulent Giant Molecular Clouds. We use Herschel observations of the Serpens Core to compute temperature and c… ▽ More Filaments represent a key structure during the early stages of the star formation process. Simulations show filamentary structure commonly formed before and during the formation of cores. Aims. The Serpens Core represents an ideal laboratory to test the state-of-the-art of simulations of turbulent Giant Molecular Clouds. We use Herschel observations of the Serpens Core to compute temperature and column density maps of the region. Among the simulations of Dale et al. (2012), we select the early stages of their Run I, before stellar feedback is initiated, with similar total mass and physical size as the Serpens Core. We derive temperature and column density maps also from the simulations. The observed distribution of column densities of the filaments has been analysed first including and then masking the cores. The same analysis has been performed on the simulations as well. A radial network of filaments has been detected in the Serpens Core. The analysed simulation shows a striking morphological resemblance to the observed structures. The column density distribution of simulated filaments without cores shows only a log-normal distribution, while the observed filaments show a power-law tail. The power-law tail becomes evident in the simulation if one focuses just on the column density distribution of the cores. In contrast, the observed cores show a flat distribution. Even though the simulated and observed filaments are subjectively similar-looking, we find that they behave in very different ways. The simulated filaments are turbulence-dominated regions, the observed filaments are instead self-gravitating structures that will probably fragment into cores. △ Less

Submitted 21 September, 2015; originally announced September 2015.

Comments: Accepted for publication on A&A. Main text: 12 pages, 7 figures

Dense gas in the Galactic central molecular zone is warm and heated by turbulence

Authors: Adam Ginsburg, Christian Henkel, Yiping Ao, Denise Riquelme, Jens Kauffmann, Thushara Pillai, Elisabeth A. C. Mills, Miguel A. Requena-Torres, Katharina Immer, Leonardo Testi, Juergen Ott, John Bally, Cara Battersby, Jeremy Darling, Susanne Aalto, Thomas Stanke, Sarah Kendrew, J. M. Diederik Kruijssen, Steven Longmore, James Dale, Rolf Guesten, Karl M. Menten

Abstract: The Galactic center is the closest region in which we can study star formation under extreme physical conditions like those in high-redshift galaxies. We measure the temperature of the dense gas in the central molecular zone (CMZ) and examine what drives it. We mapped the inner 300 pc of the CMZ in the temperature-sensitive J = 3-2 para-formaldehyde (p-H$_2$CO) transitions. We used the… ▽ More The Galactic center is the closest region in which we can study star formation under extreme physical conditions like those in high-redshift galaxies. We measure the temperature of the dense gas in the central molecular zone (CMZ) and examine what drives it. We mapped the inner 300 pc of the CMZ in the temperature-sensitive J = 3-2 para-formaldehyde (p-H$_2$CO) transitions. We used the $3_{2,1} - 2_{2,0} / 3_{0,3} - 2_{0,2}$ line ratio to determine the gas temperature in $n \sim 10^4 - 10^5 $cm$^{-3}$ gas. We have produced temperature maps and cubes with 30" and 1 km/s resolution and published all data in FITS form. Dense gas temperatures in the Galactic center range from ~60 K to > 100 K in selected regions. The highest gas temperatures T_G > 100 K are observed around the Sgr B2 cores, in the extended Sgr B2 cloud, the 20 km/s and 50 km/s clouds, and in "The Brick" (G0.253+0.016). We infer an upper limit on the cosmic ray ionization rate $ζ_{CR} < 10^{-14}$ 1/s. The dense molecular gas temperature of the region around our Galactic center is similar to values found in the central regions of other galaxies, in particular starburst systems. The gas temperature is uniformly higher than the dust temperature, confirming that dust is a coolant in the dense gas. Turbulent heating can readily explain the observed temperatures given the observed line widths. Cosmic rays cannot explain the observed variation in gas temperatures, so CMZ dense gas temperatures are not dominated by cosmic ray heating. The gas temperatures previously observed to be high in the inner ~75 pc are confirmed to be high in the entire CMZ. △ Less

Submitted 4 September, 2015; originally announced September 2015.

Comments: Accepted to A&A. The data (raw & processed) are available from http://doi.org/10.7910/DVN/27601 . The processing pipeline, analysis code, and paper source are available at https://github.com/adamginsburg/APEX_CMZ_H2CO

Journal ref: A&A 586, A50 (2016)

Isolating signatures of major cloud-cloud collisions II: The lifetimes of broad bridge features

Authors: T. J. Haworth, K. Shima, E. J. Tasker, Y. Fukui, K. Torii, J. E. Dale, K. Takahira, A. Habe

Abstract: We investigate the longevity of broad bridge features in position-velocity diagrams that appear as a result of cloud-cloud collisions. Broad bridges will have a finite lifetime due to the action of feedback, conversion of gas into stars and the timescale of the collision. We make a series of analytic arguments with which to estimate these lifetimes. Our simple analytic arguments suggest that for c… ▽ More We investigate the longevity of broad bridge features in position-velocity diagrams that appear as a result of cloud-cloud collisions. Broad bridges will have a finite lifetime due to the action of feedback, conversion of gas into stars and the timescale of the collision. We make a series of analytic arguments with which to estimate these lifetimes. Our simple analytic arguments suggest that for collisions between clouds larger than R~10 pc the lifetime of the broad bridge is more likely to be determined by the lifetime of the collision rather than the radiative or wind feedback disruption timescale. However for smaller clouds feedback becomes much more effective. This is because the radiative feedback timescale scales with the ionising flux Nly as R^{7/4}Nly^{-1/4} so a reduction in cloud size requires a relatively large decrease in ionising photons to maintain a given timescale. We find that our analytic arguments are consistent with new synthetic observations of numerical simulations of cloud-cloud collisions (including star formation and radiative feedback). We also argue that if the number of observable broad bridges remains ~ constant, then the disruption timescale must be roughly equivalent to the collision rate. If this is the case our analytic arguments also provide collision rate estimates, which we find are readily consistent with previous theoretical models at the scales they consider (clouds larger than about 10 pc) but are much higher for smaller clouds. △ Less

Submitted 2 September, 2015; originally announced September 2015.

Comments: 11 pages, 10 figures, accepted for publication in the MNRAS main journal

The modelling of feedback in star formation simulations

Authors: James E. Dale

Abstract: I review the current state of numerical simulations of stellar feedback in the context of star formation at scales ranging from the formation of individual stars to models of galaxy formation including cosmic reionisation. I survey the wealth of algorithms developed recently to solve the radiative transfer problem and to simulate stellar winds, supernovae and protostellar jets. I discuss the resul… ▽ More I review the current state of numerical simulations of stellar feedback in the context of star formation at scales ranging from the formation of individual stars to models of galaxy formation including cosmic reionisation. I survey the wealth of algorithms developed recently to solve the radiative transfer problem and to simulate stellar winds, supernovae and protostellar jets. I discuss the results of these simulations with regard to star formation in molecular clouds, the interaction of different feedback mechanisms with each other and with magnetic fields, and in the wider context of galactic-- and cosmological--scale simulations. △ Less

Submitted 25 August, 2015; originally announced August 2015.

Comments: In press at New Astronomy Reviews

Journal ref: New Astronomy Reviews (2015), pp. 1-33

StarBench: The D-type expansion of an HII region

Authors: T. G. Bisbas, T. J. Haworth, R. J. R. Williams, J. Mackey, P. Tremblin, A. C. Raga, S. J. Arthur, C. Baczynski, J. E. Dale, T. Frostholm, S. Geen, T. Haugboelle, D. Hubber, I. T. Iliev, R. Kuiper, J. Rosdahl, D. Sullivan, S. Walch, R. Wuensch

Abstract: StarBench is a project focused on benchmarking and validating different star-formation and stellar feedback codes. In this first StarBench paper we perform a comparison study of the D-type expansion of an HII region. The aim of this work is to understand the differences observed between the twelve participating numerical codes against the various analytical expressions examining the D-type phase o… ▽ More StarBench is a project focused on benchmarking and validating different star-formation and stellar feedback codes. In this first StarBench paper we perform a comparison study of the D-type expansion of an HII region. The aim of this work is to understand the differences observed between the twelve participating numerical codes against the various analytical expressions examining the D-type phase of HII region expansion. To do this, we propose two well-defined tests which are tackled by 1D and 3D grid- and SPH- based codes. The first test examines the `early phase' D-type scenario during which the mechanical pressure driving the expansion is significantly larger than the thermal pressure of the neutral medium. The second test examines the `late phase' D-type scenario during which the system relaxes to pressure equilibrium with the external medium. Although they are mutually in excellent agreement, all twelve participating codes follow a modified expansion law that deviates significantly from the classical Spitzer solution in both scenarios. We present a semi-empirical formula combining the two different solutions appropriate to both early and late phases that agrees with high-resolution simulations to $\lesssim2\%$. This formula provides a much better benchmark solution for code validation than the Spitzer solution. The present comparison has validated the participating codes and through this project we provide a dataset for calibrating the treatment of ionizing radiation hydrodynamics codes. △ Less

Submitted 20 July, 2015; originally announced July 2015.

Comments: 20 pages, 12 Figures, 4 Tables. Accepted for publication in MNRAS. Comments are welcome. Participation in future StarBench tests is also welcome

On the spatial distributions of stars and gas in numerical simulations of molecular clouds

Authors: Richard J. Parker, James E. Dale

Abstract: We compare the spatial distribution of stars which form in hydrodynamical simulations to the spatial distribution of the gas, using the $\mathcal{Q}$-parameter. The $\mathcal{Q}$-parameter enables a self-consistent comparison between the stars and gas because it uses a pixelated image of the gas as a distribution of points, in the same way that the stars (sink particles in the simulations) are a d… ▽ More We compare the spatial distribution of stars which form in hydrodynamical simulations to the spatial distribution of the gas, using the $\mathcal{Q}$-parameter. The $\mathcal{Q}$-parameter enables a self-consistent comparison between the stars and gas because it uses a pixelated image of the gas as a distribution of points, in the same way that the stars (sink particles in the simulations) are a distribution of points. We find that, whereas the stars have a substructured, or hierarchical spatial distribution ($\mathcal{Q} \sim 0.4 - 0.7$), the gas is dominated by a smooth, concentrated component and typically has $\mathcal{Q} \sim 0.9$. We also find no statistical difference between the structure of the gas in simulations that form with feedback, and those that form without, despite these two processes producing visually different distributions. These results suggest that the link between the spatial distributions of gas, and the stars which form from them, is non-trivial. △ Less

Submitted 1 June, 2015; originally announced June 2015.

Comments: 8 pages, 7 figures, accepted for publication in MNRAS

Early evolution of embedded clusters

Authors: J. E. Dale, B. Ercolano, I. A. Bonnell

Abstract: We examine the combined effects of winds and photoionizing radiation from O--type stars on embedded stellar clusters formed in model turbulent molecular clouds covering a range of masses and radii. We find that feedback is able to increase the quantities of dense gas present, but decreases the rate and efficiency of the conversion of gas to stars relative to control simulations in which feedback i… ▽ More We examine the combined effects of winds and photoionizing radiation from O--type stars on embedded stellar clusters formed in model turbulent molecular clouds covering a range of masses and radii. We find that feedback is able to increase the quantities of dense gas present, but decreases the rate and efficiency of the conversion of gas to stars relative to control simulations in which feedback is absent. Star formation in these calculations often proceeds at a rate substantially slower than the freefall rate in the dense gas. This decoupling is due to the weakening of, and expulsion of gas from, the deepest parts of the clouds' potential wells where most of the star formation occurs in the control simulations. This results in large fractions of the stellar populations in the feedback simulation becoming dissociated from dense gas. However, where star formation \emph{does} occur in both control and feedback simulations, it does so in dense gas, so the correlation between star formation activity and dense gas is preserved. The overall dynamical effects of feedback on the \emph{clusters} are minimal, with only small fraction of stars becoming unbound, despite large quantities of gas being expelled from some clouds. This owes to the settling of the stars into virialised and stellar--dominated configurations before the onset of feedback. By contrast, the effects of feedback on the observable properties of the clusters -- their U--, B-- and V--band magnitudes -- are strong and sudden. The timescales on which the clusters become visible and unobscured are short compared with the timescales which the clouds are actually destroyed. △ Less

Submitted 22 April, 2015; originally announced April 2015.

Comments: 19 pages, 17 figures (many degraded to get below submission size limit), accepted by MNRAS

The Pillars of Creation revisited with MUSE: gas kinematics and high-mass stellar feedback traced by optical spectroscopy

Authors: A. F. Mc Leod, J. E. Dale, A. Ginsburg, B. Ercolano, M. Gritschneder, S. Ramsay, L. Testi

Abstract: Integral field unit (IFU) data of the iconic Pillars of Creation in M16 are presented. The ionisation structure of the pillars was studied in great detail over almost the entire visible wavelength range, and maps of the relevant physical parameters, e.g. extinction, electron density, electron temperature, line-of-sight velocity of the ionised and neutral gas are shown. In agreement with previous a… ▽ More Integral field unit (IFU) data of the iconic Pillars of Creation in M16 are presented. The ionisation structure of the pillars was studied in great detail over almost the entire visible wavelength range, and maps of the relevant physical parameters, e.g. extinction, electron density, electron temperature, line-of-sight velocity of the ionised and neutral gas are shown. In agreement with previous authors, we find that the pillar tips are being ionised and photo-evaporated by the massive members of the nearby cluster NGC 6611. They display a stratified ionisation structure where the emission lines peak in a descending order according to their ionisation energies. The IFU data allowed us to analyse the kinematics of the photo-evaporative flow in terms of the stratified ionisation structure, and we find that, in agreement with simulations, the photo-evaporative flow is traced by a blueshift in the position-velocity profile. The gas kinematics and ionisation structure have allowed us to produce a sketch of the 3D geometry of the Pillars, positioning the pillars with respect to the ionising cluster stars. We use a novel method to detect a previously unknown bipolar outflow at the tip of the middle pillar and suggest that it has an embedded protostar as its driving source. Furthermore we identify a candidate outflow in the leftmost pillar. With the derived physical parameters and ionic abundances, we estimate a mass loss rate due to the photo-evaporative flow of 70 M$_{\odot}$ Myr$^{-1}$ which yields an expected lifetime of approximately 3 Myr. △ Less

Submitted 13 April, 2015; originally announced April 2015.

Comments: 25 pages, 21 figures, appendix 10 pages and 10 figures, MNRAS in press

Isolating signatures of major cloud-cloud collisions using position-velocity diagrams

Authors: T. J. Haworth, E. J. Tasker, Y. Fukui, K. Torii, J. E. Dale, K. Shima, K. Takahira, A. Habe, K. Hasegawa

Abstract: Collisions between giant molecular clouds are a potential mechanism for triggering the formation of massive stars, or even super star clusters. The trouble is identifying this process observationally and distinguishing it from other mechanisms. We produce synthetic position-velocity diagrams from models of: cloud-cloud collisions, non-interacting clouds along the line of sight, clouds with interna… ▽ More Collisions between giant molecular clouds are a potential mechanism for triggering the formation of massive stars, or even super star clusters. The trouble is identifying this process observationally and distinguishing it from other mechanisms. We produce synthetic position-velocity diagrams from models of: cloud-cloud collisions, non-interacting clouds along the line of sight, clouds with internal radiative feedback and a more complex cloud evolving in a galactic disc, to try and identify unique signatures of collision. We find that a broad bridge feature connecting two intensity peaks, spatially correlated but separated in velocity, is a signature of a high velocity cloud-cloud collision. We show that the broad bridge feature is resilient to the effects of radiative feedback, at least to around 2.5Myr after the formation of the first massive (ionising) star. However for a head on 10km/s collision we find that this will only be observable from 20-30 per cent of viewing angles. Such broad-bridge features have been identified towards M20, a very young region of massive star formation that was concluded to be a site of cloud-cloud collision by Torii et al (2011), and also towards star formation in the outer Milky Way by Izumi et al (2014). △ Less

Submitted 23 March, 2015; originally announced March 2015.

Comments: 12 pages. Accepted for publication in the MNRAS main journal

The dangers of being trigger--happy

Authors: J. E. Dale, T. J. Haworth, E. Bressert

Abstract: We examine the evidence offered for triggered star formation against the backdrop provided by recent numerical simulations of feedback from massive stars at or below giant molecular cloud sizescales. We compile a catalogue of sixty--seven observational papers, mostly published over the last decade, and examine the signposts most commonly used to infer the presence of triggered star formation. We t… ▽ More We examine the evidence offered for triggered star formation against the backdrop provided by recent numerical simulations of feedback from massive stars at or below giant molecular cloud sizescales. We compile a catalogue of sixty--seven observational papers, mostly published over the last decade, and examine the signposts most commonly used to infer the presence of triggered star formation. We then determine how well these signposts perform in a recent suite of hydrodynamic simulations of star formation including feedback from O--type stars performed by Dale et al (2012a, b, 2013a, b, 2014). We find that none of the observational markers improve the chances of correctly identifying a given star as triggered by more than factors of two at most. This limits the fidelity of these techniques in interpreting star formation histories. We therefore urge caution in interpreting observations of star formation near feedback--driven structures in terms of triggering. △ Less

Submitted 20 February, 2015; originally announced February 2015.

Comments: 16 pages, 6 figures, accepted by MNRAS

Studying the YMC population of M83: how long clusters remain embedded, their interaction with the ISM and implications for GC formation theories

Authors: Katherine Hollyhead, Nate Bastian, Angela Adamo, Esteban Silva-Villa, Jim Dale, Jenna Ryon, Zack Gazak

Abstract: The study of young massive clusters can provide key information for the formation of globular clusters, as they are often considered analogues. A currently unanswered question in this field is how long these massive clusters remain embedded in their natal gas, with important implications for the formation of multiple populations that have been used to explain phenomena observed in globular cluster… ▽ More The study of young massive clusters can provide key information for the formation of globular clusters, as they are often considered analogues. A currently unanswered question in this field is how long these massive clusters remain embedded in their natal gas, with important implications for the formation of multiple populations that have been used to explain phenomena observed in globular clusters. We present an analysis of ages and masses of the young massive cluster population of M83. Through visual inspection of the clusters, and comparison of their SEDs and position in colour-colour space, the clusters are all exposed (no longer embedded) by < 4 Myr, most likely less, indicating that current proposed age spreads within older clusters are unlikely. We also present several methods of constraining the ages of very young massive clusters. This can often be difficult using SED fitting due to a lack of information to disentangle age-extinction degeneracies and possible inaccurate assumptions in the models used for the fitting. The individual morphology of the Halpha around each cluster has a significant effect on the measured fluxes, which contributes to inaccuracies in the age estimates for clusters younger than 10 Myr using SED fitting. This is due to model uncertainties and aperture effects. Our methods to help constrain ages of young clusters include using the near-infrared and spectral features, such as Wolf-Rayet stars. △ Less

Submitted 12 February, 2015; originally announced February 2015.

Comments: 14 pages, accepted for publication in MNRAS

Photoionising feedback and the star formation rates in galaxies

Authors: J. M. MacLachlan, I. A. Bonnell, K. Wood, J. E. Dale

Abstract: Aims. We investigate the effects of ionising photons on accretion and stellar mass growth in a young star forming region, using a Monte Carlo radiation transfer code coupled to a smoothed particle hydrodynamics (SPH) simulation. Methods. We introduce the framework with which we correct stellar cluster masses for the effects of photoionising (PI) feedback and compare to the results of a full ionisa… ▽ More Aims. We investigate the effects of ionising photons on accretion and stellar mass growth in a young star forming region, using a Monte Carlo radiation transfer code coupled to a smoothed particle hydrodynamics (SPH) simulation. Methods. We introduce the framework with which we correct stellar cluster masses for the effects of photoionising (PI) feedback and compare to the results of a full ionisation hydrodynamics code. Results. We present results of our simulations of star formation in the spiral arm of a disk galaxy, including the effects of photoionising radiation from high mass stars. We find that PI feedback reduces the total mass accreted onto stellar clusters by approximately 23 per cent over the course of the simulation and reduces the number of high mass clusters, as well as the maximum mass attained by a stellar cluster. Mean star formation rates (SFRs) drop from 0.042 solar masses per year in our control run to 0.032 solar masses per year after the inclusion of PI feedback with a final instantaneous SFR reduction of 62 per cent. The overall cluster mass distribution appears to be affected little by PI feedback. Conclusions. We compare our results to the observed extra-galactic Schmidt-Kennicutt relation and the observed properties of local star forming regions in the Milky Way and find that internal photoionising (PI) feedback is unlikely to reduce star formation rates by more than a factor of approximately 2 and thus may play only a minor role in regulating star formation. △ Less

Submitted 22 December, 2014; originally announced December 2014.

Comments: 10 pages, 11 figures, A&A, in press

The dynamical evolution of molecular clouds near the Galactic Centre - I. Orbital structure and evolutionary timeline

Authors: J. M. Diederik Kruijssen, James E. Dale, Steven N. Longmore

Abstract: We recently proposed that the star-forming potential of dense molecular clouds in the Central Molecular Zone (CMZ, i.e. the central few 100 pc) of the Milky Way is linked to their orbital dynamics, potentially giving rise to an absolute-time sequence of star-forming clouds. In this paper, we present an orbital model for the gas stream(s) observed in the CMZ. The model is obtained by integrating or… ▽ More We recently proposed that the star-forming potential of dense molecular clouds in the Central Molecular Zone (CMZ, i.e. the central few 100 pc) of the Milky Way is linked to their orbital dynamics, potentially giving rise to an absolute-time sequence of star-forming clouds. In this paper, we present an orbital model for the gas stream(s) observed in the CMZ. The model is obtained by integrating orbits in the observed gravitational potential and represents a good fit to the distribution of dense gas, reproducing all of its key properties. The orbit is also consistent with observational constraints not included in the fitting process, such as the velocities of Sgr B2 and the Arches and Quintuplet clusters. It differs from previous models: (1) the orbit is open rather than closed due to the extended mass distribution in the CMZ, (2) its orbital velocity is twice as high as in previous models, and (3) Sgr A$^*$ coincides with the focus of the (eccentric) orbit rather than being offset. Our orbital solution supports the scenario in which the dust ridge between G0.253+0.016 ('the Brick') and Sgr B2 represents an absolute-time sequence of star-forming clouds, triggered by the tidal compression during their recent pericentre passage. We position the clouds on a common timeline and find that their pericentre passages occurred 0.30-0.74 Myr ago. Given their short free-fall times (0.3-0.4 Myr), the quiescent cloud G0.253+0.016 and the vigorously star-forming complex Sgr B2 are separated by a single free-fall time of evolution, implying that star formation proceeds rapidly once collapse has been initiated. We provide several quantitative predictions of our model and conclude with a discussion of the model in the Galactic context, highlighting its relation to large-scale gas accretion, the dynamics of the bar, the $x_2$ orbital family, and the origin of the Arches and Quintuplet clusters. (Abridged) △ Less

Submitted 1 December, 2014; originally announced December 2014.

Comments: 22 pages (including appendices), 11 figures, 4 tables; accepted by MNRAS (27 November 2014). A high-resolution, machine-readable version of Table C1 will be available in the Supporting Information of the published paper; for earlier access, please contact the first author

Turbulence in Giant Molecular Clouds: The effect of photoionisation feedback

Authors: D. M. Boneberg, J. E. Dale, P. Girichidis, B. Ercolano

Abstract: Giant Molecular Clouds (GMCs) are observed to be turbulent, but theory shows that without a driving mechanism turbulence should quickly decay. The question arises by which mechanisms turbulence is driven or sustained. It has been shown that photoionising feedback from massive stars has an impact on the surrounding GMC and can for example create vast HII bubbles. We therefore address the question o… ▽ More Giant Molecular Clouds (GMCs) are observed to be turbulent, but theory shows that without a driving mechanism turbulence should quickly decay. The question arises by which mechanisms turbulence is driven or sustained. It has been shown that photoionising feedback from massive stars has an impact on the surrounding GMC and can for example create vast HII bubbles. We therefore address the question of whether turbulence is a consequence of this effect of feedback on the cloud. To investigate this, we analyse the velocity field of simulations of high mass star forming regions by studying velocity structure functions and power spectra. We find that clouds whose morphology is strongly affected by photoionising feedback also show evidence of driving of turbulence by preserving or recovering a Kolmogorov-type velocity field. On the contrary, control run simulations without photoionising feedback have a velocity distribution that bears the signature of gravitational collapse and of the dissipation of energy, where the initial Kolmogorov-type structure function is erased. △ Less

Submitted 28 November, 2014; originally announced November 2014.

Comments: accepted by MNRAS, 13 pages, 12 figures

Primordial mass segregation in simulations of star formation?

Authors: Richard J. Parker, James E. Dale, Barbara Ercolano

Abstract: We take the end result of smoothed particle hydrodynamics (SPH) simulations of star formation which include feedback from photoionisation and stellar winds and evolve them for a further 10Myr using $N$-body simulations. We compare the evolution of each simulation to a control run without feedback, and to a run with photoionisation feedback only. In common with previous work, we find that the prese… ▽ More We take the end result of smoothed particle hydrodynamics (SPH) simulations of star formation which include feedback from photoionisation and stellar winds and evolve them for a further 10Myr using $N$-body simulations. We compare the evolution of each simulation to a control run without feedback, and to a run with photoionisation feedback only. In common with previous work, we find that the presence of feedback prevents the runaway growth of massive stars, and the resulting star-forming regions are less dense, and preserve their initial substructure for longer. The addition of stellar winds to the feedback produces only marginal differences compared to the simulations with just photoionisation feedback. We search for mass segregation at different stages in the simulations; before feedback is switched on in the SPH runs, at the end of the SPH runs (before $N$-body integration) and during the $N$-body evolution. Whether a simulation is primordially mass segregated (i.e. before dynamical evolution) depends extensively on how mass segregation is defined, and different methods for measuring mass segregation give apparently contradictory results. Primordial mass segregation is also less common in the simulations when star formation occurs under the influence of feedback. Further dynamical mass segregation can also take place during the subsequent (gas-free) dynamical evolution. Taken together, our results suggest that extreme caution should be exercised when interpreting the spatial distribution of massive stars relative to low-mass stars in simulations. △ Less

Submitted 11 November, 2014; originally announced November 2014.

Comments: 14 pages, 11 figures, accepted for publication in MNRAS

First Investigation of the Combined Impact of Ionizing Radiation and Momentum Winds from a Massive Star on a Self-Gravitating Core

Authors: Judith Ngoumou, David Hubber, James E. Dale, Andreas Burkert

Abstract: Massive stars shape the surrounding ISM by emitting ionizing photons and ejecting material through stellar winds. To study the impact of the momentum from the wind of a massive star on the surrounding neutral or ionized material, we implemented a new HEALPix-based momentum conserving wind scheme in the Smoothed Particle Hydrodynamics (SPH) code SEREN. A qualitative study of the impact of the feedb… ▽ More Massive stars shape the surrounding ISM by emitting ionizing photons and ejecting material through stellar winds. To study the impact of the momentum from the wind of a massive star on the surrounding neutral or ionized material, we implemented a new HEALPix-based momentum conserving wind scheme in the Smoothed Particle Hydrodynamics (SPH) code SEREN. A qualitative study of the impact of the feedback from an O7.5-like star on a self gravitating sphere shows that, on its own, the transfer of momentum from a wind onto cold surrounding gas has both a compressing and dispersing effect. It mostly affects gas at low and intermediate densities. When combined with a stellar source's ionizing UV radiation, we find the momentum driven wind to have little direct effect on the gas. We conclude that, during a massive star's main sequence, the UV ionizing radiation is the main feedback mechanism shaping and compressing the cold gas. Overall, the wind's effects on the dense gas dynamics and on the triggering of star formation are very modest. The structures formed in the ionization-only simulation and in the combined feedback simulation are remarkably similar. However, in the combined feedback case, different SPH particles end up being compressed. This indicates that the microphysics of gas mixing differ between the two feedback simulations and that the winds can contribute to the localized redistribution and reshuffling of gas. △ Less

Submitted 12 December, 2014; v1 submitted 20 October, 2014; originally announced October 2014.

Comments: Accepted for publication in ApJ (12 pages, 15 figures)

Old pre-main-sequence Stars: Disc reformation by Bondi-Hoyle accretion

Authors: P. Scicluna, G. Rosotti, J. E. Dale, L. Testi

Abstract: Young stars show evidence of accretion discs which evolve quickly and disperse with an e-folding time of $\sim$ 3Myr. This is in striking contrast with recent observations that suggest evidence for numerous $>30$ Myr old stars with an accretion disc in large star-forming complexes. We consider whether these observations of apparently old accretors could be explained by invoking Bondi-Hoyle accreti… ▽ More Young stars show evidence of accretion discs which evolve quickly and disperse with an e-folding time of $\sim$ 3Myr. This is in striking contrast with recent observations that suggest evidence for numerous $>30$ Myr old stars with an accretion disc in large star-forming complexes. We consider whether these observations of apparently old accretors could be explained by invoking Bondi-Hoyle accretion to rebuild a new disc around these stars during passage through a clumpy molecular cloud. We combine a simple Monte Carlo model to explore the capture of mass by such systems with a viscous evolution model to infer the levels of accretion that would be observed. We find that a significant fraction of stars may capture enough material via the Bondi-Hoyle mechanism to rebuild a disc of mass $\gtrsim$ 1 minimum-mass solar nebula, and $\lesssim 10\%$ accrete at observable levels at any given time. A significant fraction of the observed old accretors may be explained with our proposed mechanism. Such accretion may provide a chance for a second epoch of planet formation, and have unpredictable consequences for planetary evolution. △ Less

Submitted 23 May, 2014; originally announced May 2014.

Comments: 5 pages, 2 figures. Accepted by A&A

Before the first supernova: combined effects of HII regions and winds on molecular clouds

Authors: J. E. Dale, J. Ngoumou, B. Ercolano, I. A. Bonnell

Abstract: We model the combined effects of photoionization and momentum--driven winds from O--stars on molecular clouds spanning a parameter space of initial conditions. The dynamical effects of the winds are very modest. However, in the lower--mass clouds, they influence the morphologies of the HII regions by creating 10pc--scale central cavities.\\ The inhomogeneous structures of the model GMCs make them… ▽ More We model the combined effects of photoionization and momentum--driven winds from O--stars on molecular clouds spanning a parameter space of initial conditions. The dynamical effects of the winds are very modest. However, in the lower--mass clouds, they influence the morphologies of the HII regions by creating 10pc--scale central cavities.\\ The inhomogeneous structures of the model GMCs make them highly permeable to photons, ionized gas and supernova ejecta, and the leaking of ionized gas in particular strongly affects their evolution, reducing the effectiveness of feedback. Nevertheless, feedback is able to expel large fractions of the mass of the lower escape--velocity clouds. Its impact on star formation is more modest, decreasing final star formation efficiencies by 10--20$\%$, and the rate of change of the star formation efficiency per freefall time by about one third. However, the clouds still form stars substantially faster than observed GMCs. △ Less

Submitted 24 April, 2014; originally announced April 2014.

Comments: 21 pages, 19 figures, accepted for publication in MNRAS