-
Self-consistent modelling of the Milky Way structure using live potentials
Authors:
Eva Durán-Camacho,
Ana Duarte-Cabral,
Alex R. Pettitt,
Robin G. Treß,
Paul C. Clark,
Ralf S. Klessen,
Kamran R. J. Bogue,
Rowan J. Smith,
Mattia C. Sormani
Abstract:
To advance our understanding of the evolution of the interstellar medium (ISM) of our Galaxy, numerical models of Milky Way (MW) type galaxies are widely used. However, most models only vaguely resemble the MW (e.g. in total mass), and often use imposed analytic potentials (which cannot evolve dynamically). This poses a problem in asserting their applicability for the interpretation of observation…
▽ More
To advance our understanding of the evolution of the interstellar medium (ISM) of our Galaxy, numerical models of Milky Way (MW) type galaxies are widely used. However, most models only vaguely resemble the MW (e.g. in total mass), and often use imposed analytic potentials (which cannot evolve dynamically). This poses a problem in asserting their applicability for the interpretation of observations of our own Galaxy. The goal of this work is to identify a numerical model that is not only a MW-type galaxy, but one that can mimic some of the main observed structures of our Galaxy, using dynamically evolving potentials, so that it can be used as a base model to study the ISM cycle in a galaxy like our own. This paper introduces a suite of 15 MW-type galaxy models developed using the {\sc arepo} numerical code, that are compared to Galactic observations of $^{12}$CO and \ion{H}{I} emission via longitude-velocity plots, from where we extract and compare the skeletons of major galactic features and the terminal gas velocities. We found that our best-fitting model to the overall structure, also reproduces some of the more specific observed features of the MW, including a bar with a pattern speed of $30.0 \pm 0.2$ km\,s$^{-1}$\,kpc$^{-1}$, a bar half-length of $3.2 \pm 0.8$\,kpc. Our model shows large streaming motions around spiral arms, and strong radial motions well beyond the inner bar. This model highlights the complex motions of a dynamic MW-type galaxy and has the potential to offer valuable insight into how our Galaxy regulates the ISM and star formation.
△ Less
Submitted 11 June, 2024; v1 submitted 15 May, 2024;
originally announced May 2024.
-
Magnetic field morphology and evolution in the Central Molecular Zone and its effect on gas dynamics
Authors:
R. G. Tress,
M. C. Sormani,
P. Girichidis,
S. C. O. Glover,
R. S. Klessen,
R. J. Smith,
E. Sobacchi,
L. Armillotta,
A. T. Barnes,
C. Battersby,
K. R. J. Bogue,
N. Brucy,
L. Colzi,
C. Federrath,
P. García,
A. Ginsburg,
J. Göller,
H P. Hatchfield,
C. Henkel,
P. Hennebelle,
J. D. Henshaw,
M. Hirschmann,
Y. Hu,
J. Kauffmann,
J. M. D. Kruijssen
, et al. (12 additional authors not shown)
Abstract:
The interstellar medium in the Milky Way's Central Molecular Zone (CMZ) is known to be strongly magnetised, but its large-scale morphology and impact on the gas dynamics are not well understood. We explore the impact and properties of magnetic fields in the CMZ using three-dimensional non-self gravitating magnetohydrodynamical simulations of gas flow in an external Milky Way barred potential. We f…
▽ More
The interstellar medium in the Milky Way's Central Molecular Zone (CMZ) is known to be strongly magnetised, but its large-scale morphology and impact on the gas dynamics are not well understood. We explore the impact and properties of magnetic fields in the CMZ using three-dimensional non-self gravitating magnetohydrodynamical simulations of gas flow in an external Milky Way barred potential. We find that: (1) The magnetic field is conveniently decomposed into a regular time-averaged component and an irregular turbulent component. The regular component aligns well with the velocity vectors of the gas everywhere, including within the bar lanes. (2) The field geometry transitions from parallel to the Galactic plane near $z=0$ to poloidal away from the plane. (3) The magneto-rotational instability (MRI) causes an in-plane inflow of matter from the CMZ gas ring towards the central few parsecs of $0.01-0.1$ M$_\odot$ yr$^{-1}$ that is absent in the unmagnetised simulations. However, the magnetic fields have no significant effect on the larger-scale bar-driven inflow that brings the gas from the Galactic disc into the CMZ. (4) A combination of bar inflow and MRI-driven turbulence can sustain a turbulent vertical velocity dispersion of $σ_z \simeq 5$ km s$^{-1}$ on scales of $20$ pc in the CMZ ring. The MRI alone sustains a velocity dispersion of $σ_z \simeq 3$ km s$^{-1}$. Both these numbers are lower than the observed velocity dispersion of gas in the CMZ, suggesting that other processes such as stellar feedback are necessary to explain the observations. (5) Dynamo action driven by differential rotation and the MRI amplifies the magnetic fields in the CMZ ring until they saturate at a value that scales with the average local density as $B \simeq 102 (n/10^3 cm^{-3})^{0.33}$ $μ$G. Finally, we discuss the implications of our results within the observational context in the CMZ.
△ Less
Submitted 19 March, 2024;
originally announced March 2024.
-
Cloud properties across spatial scales in simulations of the interstellar medium
Authors:
Tine Colman,
Noé Brucy,
Philipp Girichidis,
Simon C. O Glover,
Milena Benedettini,
Juan D. Soler,
Robin G. Tress,
Alessio Traficante,
Patrick Hennebelle,
Ralf S. Klessen,
Sergio Molinari,
Marc-Antoine Miville-Deschênes
Abstract:
Molecular clouds (MC) are structures of dense gas in the interstellar medium (ISM), that extend from ten to a few hundred parsecs and form the main gas reservoir available for star formation. Hydrodynamical simulations of varying complexity are a promising way to investigate MC evolution and their properties. However, each simulation typically has a limited range in resolution and different cloud…
▽ More
Molecular clouds (MC) are structures of dense gas in the interstellar medium (ISM), that extend from ten to a few hundred parsecs and form the main gas reservoir available for star formation. Hydrodynamical simulations of varying complexity are a promising way to investigate MC evolution and their properties. However, each simulation typically has a limited range in resolution and different cloud extraction algorithms are used, which complicates the comparison between simulations. In this work, we aim to extract clouds from different simulations covering a wide range of spatial scales. We compare their properties, such as size, shape, mass, internal velocity dispersion and virial state. We apply the Hop cloud detection algorithm on (M)HD numerical simulations of stratified ISM boxes and isolated galactic disk simulations that were produced using Flash Ramses and Arepo We find that the extracted clouds are complex in shape ranging from round objects to complex filamentary networks in all setups. Despite the wide range of scales, resolution, and sub-grid physics, we observe surprisingly robust trends in the investigated metrics. The mass spectrum matches in the overlap between simulations without rescaling and with a high-mass slope of $\mathrm{d} N/\mathrm{d}\ln M\propto-1$ in accordance with theoretical predictions. The internal velocity dispersion scales with the size of the cloud as $σ\propto R^{0.75}$ for large clouds ($R\gtrsim3\,\mathrm{pc}$). For small clouds we find larger sigma compared to the power-law scaling, as seen in observations, which is due to supernova-driven turbulence. Almost all clouds are gravitationally unbound with the virial parameter scaling as $α_\mathrm{vir}\propto M^{-0.4}$, which is slightly flatter compared to observed scaling, but in agreement given the large scatter.
△ Less
Submitted 1 March, 2024;
originally announced March 2024.
-
Hidden Gems on a Ring: Infant Massive Clusters and Their Formation Timeline Unveiled by ALMA, HST, and JWST in NGC 3351
Authors:
Jiayi Sun,
Hao He,
Kyle Batschkun,
Rebecca C. Levy,
Kimberly Emig,
M. Jimena Rodriguez,
Hamid Hassani,
Adam K. Leroy,
Eva Schinnerer,
Eve C. Ostriker,
Christine D. Wilson,
Alberto D. Bolatto,
Elisabeth A. C. Mills,
Erik Rosolowsky,
Janice C. Lee,
Daniel A. Dale,
Kirsten L. Larson,
David A. Thilker,
Leonardo Ubeda,
Bradley C. Whitmore,
Thomas G. Williams,
Ashley. T. Barnes,
Frank Bigiel,
Melanie Chevance,
Simon C. O. Glover
, et al. (16 additional authors not shown)
Abstract:
We study young massive clusters (YMCs) in their embedded "infant" phase with $\sim0.\!^{\prime\prime}1$ ALMA, HST, and JWST observations targeting the central starburst ring in NGC 3351, a nearby Milky Way analog galaxy. Our new ALMA data reveal 18 bright and compact (sub-)millimeter continuum sources, of which 8 have counterparts in JWST images and only 6 have counterparts in HST images. Based on…
▽ More
We study young massive clusters (YMCs) in their embedded "infant" phase with $\sim0.\!^{\prime\prime}1$ ALMA, HST, and JWST observations targeting the central starburst ring in NGC 3351, a nearby Milky Way analog galaxy. Our new ALMA data reveal 18 bright and compact (sub-)millimeter continuum sources, of which 8 have counterparts in JWST images and only 6 have counterparts in HST images. Based on the ALMA continuum and molecular line data, as well as ancillary measurements for the HST and JWST counterparts, we identify 14 sources as infant star clusters with high stellar and/or gas masses (${\sim}10^5\;\mathrm{M_\odot}$), small radii (${\lesssim}\,5\;\mathrm{pc}$), large escape velocities ($6{-}10\;\mathrm{km/s}$), and short free-fall times ($0.5{-}1\;\mathrm{Myr}$). Their multiwavelength properties motivate us to divide them into four categories, likely corresponding to four evolutionary stages from starless clumps to exposed HII region-cluster complexes. Leveraging age estimates for HST-identified clusters in the same region, we infer an evolutionary timeline going from $\sim$1-2 Myr before cluster formation as starless clumps, to $\sim$4-6 Myr after as exposed HII region-cluster complexes. Finally, we show that the YMCs make up a substantial fraction of recent star formation across the ring, exhibit an non-uniform azimuthal distribution without a very coherent evolutionary trend along the ring, and are capable of driving large-scale gas outflows.
△ Less
Submitted 10 April, 2024; v1 submitted 25 January, 2024;
originally announced January 2024.
-
The JWST Galactic Center Survey -- A White Paper
Authors:
Rainer Schoedel,
Steve Longmore,
Jonny Henshaw,
Adam Ginsburg,
John Bally,
Anja Feldmeier,
Matt Hosek,
Francisco Nogueras Lara,
Anna Ciurlo,
Mélanie Chevance,
J. M. Diederik Kruijssen,
Ralf Klessen,
Gabriele Ponti,
Pau Amaro-Seoane,
Konstantina Anastasopoulou,
Jay Anderson,
Maria Arias,
Ashley T. Barnes,
Cara Battersby,
Giuseppe Bono,
Lucía Bravo Ferres,
Aaron Bryant,
Miguel Cano Gonzáalez,
Santi Cassisi,
Leonardo Chaves-Velasquez
, et al. (85 additional authors not shown)
Abstract:
The inner hundred parsecs of the Milky Way hosts the nearest supermassive black hole, largest reservoir of dense gas, greatest stellar density, hundreds of massive main and post main sequence stars, and the highest volume density of supernovae in the Galaxy. As the nearest environment in which it is possible to simultaneously observe many of the extreme processes shaping the Universe, it is one of…
▽ More
The inner hundred parsecs of the Milky Way hosts the nearest supermassive black hole, largest reservoir of dense gas, greatest stellar density, hundreds of massive main and post main sequence stars, and the highest volume density of supernovae in the Galaxy. As the nearest environment in which it is possible to simultaneously observe many of the extreme processes shaping the Universe, it is one of the most well-studied regions in astrophysics. Due to its proximity, we can study the center of our Galaxy on scales down to a few hundred AU, a hundred times better than in similar Local Group galaxies and thousands of times better than in the nearest active galaxies. The Galactic Center (GC) is therefore of outstanding astrophysical interest. However, in spite of intense observational work over the past decades, there are still fundamental things unknown about the GC. JWST has the unique capability to provide us with the necessary, game-changing data. In this White Paper, we advocate for a JWST NIRCam survey that aims at solving central questions, that we have identified as a community: i) the 3D structure and kinematics of gas and stars; ii) ancient star formation and its relation with the overall history of the Milky Way, as well as recent star formation and its implications for the overall energetics of our galaxy's nucleus; and iii) the (non-)universality of star formation and the stellar initial mass function. We advocate for a large-area, multi-epoch, multi-wavelength NIRCam survey of the inner 100\,pc of the Galaxy in the form of a Treasury GO JWST Large Program that is open to the community. We describe how this survey will derive the physical and kinematic properties of ~10,000,000 stars, how this will solve the key unknowns and provide a valuable resource for the community with long-lasting legacy value.
△ Less
Submitted 14 March, 2024; v1 submitted 18 October, 2023;
originally announced October 2023.
-
Fuelling the nuclear ring of NGC 1097
Authors:
Mattia C. Sormani,
Ashley T. Barnes,
Jiayi Sun,
Sophia K. Stuber,
Eva Schinnerer,
Eric Emsellem,
Adam K. Leroy,
Simon C. O. Glover,
Jonathan D. Henshaw,
Sharon E. Meidt,
Justus Neumann,
Miguel Querejeta,
Thomas G. Williams,
Frank Bigiel,
Cosima Eibensteiner,
Francesca Fragkoudi,
Rebecca C. Levy,
Kathryn Grasha,
Ralf S. Klessen,
J. M. Diederik Kruijssen,
Nadine Neumayer,
Francesca Pinna,
Erik W. Rosolowsky,
Rowan J. Smith,
Yu-Hsuan Teng
, et al. (2 additional authors not shown)
Abstract:
Galactic bars can drive cold gas inflows towards the centres of galaxies. The gas transport happens primarily through the so-called bar ``dust lanes'', which connect the galactic disc at kpc scales to the nuclear rings at hundreds of pc scales much like two gigantic galactic rivers. Once in the ring, the gas can fuel star formation activity, galactic outflows, and central supermassive black holes.…
▽ More
Galactic bars can drive cold gas inflows towards the centres of galaxies. The gas transport happens primarily through the so-called bar ``dust lanes'', which connect the galactic disc at kpc scales to the nuclear rings at hundreds of pc scales much like two gigantic galactic rivers. Once in the ring, the gas can fuel star formation activity, galactic outflows, and central supermassive black holes. Measuring the mass inflow rates is therefore important to understanding the mass/energy budget and evolution of galactic nuclei. In this work, we use CO datacubes from the PHANGS-ALMA survey and a simple geometrical method to measure the bar-driven mass inflow rate onto the nuclear ring of the barred galaxy NGC~1097. The method assumes that the gas velocity in the bar lanes is parallel to the lanes in the frame co-rotating with the bar, and allows one to derive the inflow rates from sufficiently sensitive and resolved position-position-velocity diagrams if the bar pattern speed and galaxy orientations are known. We find an inflow rate of $\dot{M}=(3.0 \pm 2.1)\, \rm M_\odot\, yr^{-1}$ averaged over a time span of 40 Myr, which varies by a factor of a few over timescales of $\sim$10 Myr. Most of the inflow appears to be consumed by star formation in the ring which is currently occurring at a rate of ${\rm SFR}\simeq~1.8$-$2 \rm M_\odot\, yr^{-1}$, suggesting that the inflow is causally controlling the star formation rate in the ring as a function of time.
△ Less
Submitted 23 May, 2023;
originally announced May 2023.
-
Kinematics of Galactic Centre clouds shaped by shear-seeded solenoidal turbulence
Authors:
Maya A. Petkova,
J. M. Diederik Kruijssen,
Jonathan D. Henshaw,
Steven N. Longmore,
Simon C. O. Glover,
Mattia C. Sormani,
Lucia Armillotta,
Ashley T. Barnes,
Ralf S. Klessen,
Francisco Nogueras-Lara,
Robin G. Tress,
Jairo Armijos-Abendaño,
Laura Colzi,
Christoph Federrath,
Pablo García,
Adam Ginsburg,
Christian Henkel,
Sergio Martín,
Denise Riquelme,
Víctor M. Rivilla
Abstract:
The Central Molecular Zone (CMZ; the central ~ 500 pc of the Galaxy) is a kinematically unusual environment relative to the Galactic disc, with high velocity dispersions and a steep size-linewidth relation of the molecular clouds. In addition, the CMZ region has a significantly lower star formation rate (SFR) than expected by its large amount of dense gas. An important factor in explaining the low…
▽ More
The Central Molecular Zone (CMZ; the central ~ 500 pc of the Galaxy) is a kinematically unusual environment relative to the Galactic disc, with high velocity dispersions and a steep size-linewidth relation of the molecular clouds. In addition, the CMZ region has a significantly lower star formation rate (SFR) than expected by its large amount of dense gas. An important factor in explaining the low SFR is the turbulent state of the star-forming gas, which seems to be dominated by rotational modes. However, the turbulence driving mechanism remains unclear. In this work, we investigate how the Galactic gravitational potential affects the turbulence in CMZ clouds. We focus on the CMZ cloud G0.253+0.016 (`the Brick'), which is very quiescent and unlikely to be kinematically dominated by stellar feedback. We demonstrate that several kinematic properties of the Brick arise naturally in a cloud-scale hydrodynamics simulation that takes into account the Galactic gravitational potential. These properties include the line-of-sight velocity distribution, the steepened size-linewidth relation, and the predominantly solenoidal nature of the turbulence. Within the simulation, these properties result from the Galactic shear in combination with the cloud's gravitational collapse. This is a strong indication that the Galactic gravitational potential plays a crucial role in shaping the CMZ gas kinematics, and is a major contributor to suppressing the SFR by inducing predominantly solenoidal turbulent modes.
△ Less
Submitted 1 August, 2023; v1 submitted 21 April, 2023;
originally announced April 2023.
-
PHANGS-JWST First Results: Rapid Evolution of Star Formation in the Central Molecular Gas Ring of NGC1365
Authors:
Eva Schinnerer,
Eric Emsellem,
Jonathan D. Henshaw,
Daizhong Liu,
Sharon E. Meidt,
Miguel Querejeta,
Florent Renaud,
Mattia C. Sormani,
Jiayi Sun,
Oleg V. Egorov,
Kirsten L. Larson,
Adam K. Leroy,
Erik Rosolowsky,
Karin M. Sandstrom,
T. G. Williams,
Ashley T. Barnes,
F. Bigiel,
Melanie Chevance,
Yixian Cao,
Rupali Chandar,
Daniel A. Dale,
Cosima Eibensteiner,
Simon C. O. Glover,
Kathryn Grasha,
Stephen Hannon
, et al. (14 additional authors not shown)
Abstract:
Large-scale bars can fuel galaxy centers with molecular gas, often leading to the development of dense ring-like structures where intense star formation occurs, forming a very different environment compared to galactic disks. We pair ~0.3" (30pc) resolution new JWST/MIRI imaging with archival ALMA CO(2-1) mapping of the central ~5kpc of the nearby barred spiral galaxy NGC1365, to investigate the p…
▽ More
Large-scale bars can fuel galaxy centers with molecular gas, often leading to the development of dense ring-like structures where intense star formation occurs, forming a very different environment compared to galactic disks. We pair ~0.3" (30pc) resolution new JWST/MIRI imaging with archival ALMA CO(2-1) mapping of the central ~5kpc of the nearby barred spiral galaxy NGC1365, to investigate the physical mechanisms responsible for this extreme star formation. The molecular gas morphology is resolved into two well-known bright bar lanes that surround a smooth dynamically cold gas disk (R_gal ~ 475pc) reminiscent of non-star-forming disks in early type galaxies and likely fed by gas inflow triggered by stellar feedback in the lanes. The lanes host a large number of JWST-identified massive young star clusters. We find some evidence for temporal star formation evolution along the ring. The complex kinematics in the gas lanes reveal strong streaming motions and may be consistent with convergence of gas streamlines expected there. Indeed, the extreme line-widths are found to be the result of inter-`cloud' motion between gas peaks; ScousePy decomposition reveals multiple components with line widths of <sigma_CO,scouse> ~ 19km/s and surface densities of <Sigma_H2,scouse> ~ 800M_sun/pc^2, similar to the properties observed throughout the rest of the central molecular gas structure. Tailored hydro-dynamical simulations exhibit many of the observed properties and imply that the observed structures are transient and highly time-variable. From our study of NGC1365, we conclude that it is predominantly the high gas inflow triggered by the bar that is setting the star formation in its CMZ.
△ Less
Submitted 18 December, 2022;
originally announced December 2022.
-
The PHANGS-JWST Treasury Survey: Star Formation, Feedback, and Dust Physics at High Angular resolution in Nearby GalaxieS
Authors:
Janice C. Lee,
Karin M. Sandstrom,
Adam K. Leroy,
David A. Thilker,
Eva Schinnerer,
Erik Rosolowsky,
Kirsten L. Larson,
Oleg V. Egorov,
Thomas G. Williams,
Judy Schmidt,
Eric Emsellem,
Gagandeep S. Anand,
Ashley T. Barnes,
Francesco Belfiore,
Ivana Beslic,
Frank Bigiel,
Guillermo A. Blanc,
Alberto D. Bolatto,
Mederic Boquien,
Jakob den Brok,
Yixian Cao,
Rupali Chandar,
Jeremy Chastenet,
Melanie Chevance,
I-Da Chiang
, et al. (52 additional authors not shown)
Abstract:
The PHANGS collaboration has been building a reference dataset for the multi-scale, multi-phase study of star formation and the interstellar medium in nearby galaxies. With the successful launch and commissioning of JWST, we can now obtain high-resolution infrared imaging to probe the youngest stellar populations and dust emission on the scales of star clusters and molecular clouds ($\sim$5-50 pc)…
▽ More
The PHANGS collaboration has been building a reference dataset for the multi-scale, multi-phase study of star formation and the interstellar medium in nearby galaxies. With the successful launch and commissioning of JWST, we can now obtain high-resolution infrared imaging to probe the youngest stellar populations and dust emission on the scales of star clusters and molecular clouds ($\sim$5-50 pc). In Cycle 1, PHANGS is conducting an 8-band imaging survey from 2-21$μ$m of 19 nearby spiral galaxies. CO(2-1) mapping, optical integral field spectroscopy, and UV-optical imaging for all 19 galaxies have been obtained through large programs with ALMA, VLT/MUSE, and Hubble. PHANGS-JWST enables a full inventory of star formation, accurate measurement of the mass and age of star clusters, identification of the youngest embedded stellar populations, and characterization of the physical state of small dust grains. When combined with Hubble catalogs of $\sim$10,000 star clusters, MUSE spectroscopic mapping of $\sim$20,000 HII regions, and $\sim$12,000 ALMA-identified molecular clouds, it becomes possible to measure the timescales and efficiencies of the earliest phases of star formation and feedback, build an empirical model of the dependence of small dust grain properties on local ISM conditions, and test our understanding of how dust-reprocessed starlight traces star formation activity, all across a diversity of galactic environments. Here we describe the PHANGS-JWST Treasury survey, present the remarkable imaging obtained in the first few months of science operations, and provide context for the initial results presented in the first series of PHANGS-JWST publications.
△ Less
Submitted 5 December, 2022;
originally announced December 2022.
-
PHANGS-JWST First Results: Multi-wavelength view of feedback-driven bubbles (The Phantom Voids) across NGC 628
Authors:
Ashley T. Barnes,
Elizabeth J. Watkins,
Sharon E. Meidt,
Kathryn Kreckel,
Mattia C. Sormani,
Robin G. Tress,
Simon C. O. Glover,
Frank Bigiel,
Rupali Chandar,
Eric Emsellem,
Janice C. Lee,
Adam K. Leroy,
Karin M. Sandstrom,
Eva Schinnerer,
Erik W. Rosolowsky,
Francesco Belfiore,
Guillermo Blanc,
Mederic Boquien,
Jakob S. den Brok,
Yixian Cao,
Mélanie Chevance,
Daniel A. Dale,
Oleg Egorov,
Cosima Eibensteiner,
Kathryn Grasha
, et al. (29 additional authors not shown)
Abstract:
We present a high-resolution view of bubbles within The Phantom Galaxy (NGC 628); a nearby (~10Mpc), star-forming (~2Msun/yr), face-on (i~9deg) grand-design spiral galaxy. With new data obtained as part of the PHANGS-JWST treasury program, we perform a detailed case-study of two regions of interest, one of which contains the largest and most prominent bubble in the galaxy (The Phantom Void; over 1…
▽ More
We present a high-resolution view of bubbles within The Phantom Galaxy (NGC 628); a nearby (~10Mpc), star-forming (~2Msun/yr), face-on (i~9deg) grand-design spiral galaxy. With new data obtained as part of the PHANGS-JWST treasury program, we perform a detailed case-study of two regions of interest, one of which contains the largest and most prominent bubble in the galaxy (The Phantom Void; over 1kpc in diameter), and the other being a smaller region that may be the precursor to such a large bubble (The Precursor Phantom Void). When comparing to matched resolution Halpha observations from the Hubble Space Telescope (HST), we see that the ionized gas is brightest in the shells of both bubbles, and is coincident with the youngest (~1Myr) and most massive (~100,000Msun) stellar associations. We also find an older generation (~20Myr) of stellar associations is present within the bubble of The Phantom Void. From our kinematic analysis of the HI, H2 (CO) and HII gas across The Phantom Void, we infer a high expansion speed of around 15 to 50km/s. The large size and high expansion speed of The Phantom Void suggest that the driving mechanism is sustained stellar feedback due to multiple mechanisms, where early feedback first cleared a bubble (as we observe now in The Precursor Phantom Void), and since then SNe have been exploding within the cavity, and have accelerated the shell. Finally, comparison to simulations shows a striking resemblance to our JWST observations, and suggests that such large-scale stellar feedback-driven bubbles should be common within other galaxies.
△ Less
Submitted 1 December, 2022;
originally announced December 2022.
-
Metal Mixing in Minihalos: The Descendants of Pair-Instability Supernovae
Authors:
Mattis Magg,
Anna T. P. Schauer,
Ralf S. Klessen,
Simon C. O. Glover,
Robin G. Tress,
Ondrej Jaura
Abstract:
The lack of observations of abundance patterns originating in pair-instability supernovae has been a long-standing problem in relation to the first stars. This class of supernovae is expected to have an abundance pattern with a strong odd-even effect, making it substantially different from present-day supernovae. In this study, we use a cosmological radiation hydrodynamics simulation to model such…
▽ More
The lack of observations of abundance patterns originating in pair-instability supernovae has been a long-standing problem in relation to the first stars. This class of supernovae is expected to have an abundance pattern with a strong odd-even effect, making it substantially different from present-day supernovae. In this study, we use a cosmological radiation hydrodynamics simulation to model such supernovae and the subsequent formation of the second generation of stars. We incorporate streaming velocities for the first time. There are 14 star-forming minihalos in our $1\,\mathrm{Mpc}\,h^{-1}$ box, leading to 14 supernovae occurring before redshift $z=19.5$, where we start reducing the complexity of the simulation. Following the explosions, extremely metal-poor stars form in 10 halos via internal and external enrichment, which makes it the most common outcome. Only one halo does not recollapse during the simulations. This result is at tension with the current (lack of) observations of metal-poor stars with pair instability supernova abundance patterns, suggesting that these very massive stars might be rare even in the early Universe. The results from this simulation also give us insights into what drives different modes of recollapse and what determines the mixing behavior of metals after very energetic supernovae.
△ Less
Submitted 2 March, 2022; v1 submitted 28 October, 2021;
originally announced October 2021.
-
Dynamically Driven Inflow onto the Galactic Center and its Effect upon Molecular Clouds
Authors:
H Perry Hatchfield,
Mattia C. Sormani,
Robin G. Tress,
Cara Battersby,
Rowan J. Smith,
Simon C. O. Glover,
Ralf S. Klessen
Abstract:
The Galactic bar plays a critical role in the evolution of the Milky Way's Central Molecular Zone (CMZ), driving mass inward toward the Galactic Center via gas flows known as dust lanes. To explore the interaction between the CMZ and the dust lanes, we run hydrodynamic simulations in AREPO, modeling the potential of the Milky Way's bar in the absence of gas self-gravity and star formation physics,…
▽ More
The Galactic bar plays a critical role in the evolution of the Milky Way's Central Molecular Zone (CMZ), driving mass inward toward the Galactic Center via gas flows known as dust lanes. To explore the interaction between the CMZ and the dust lanes, we run hydrodynamic simulations in AREPO, modeling the potential of the Milky Way's bar in the absence of gas self-gravity and star formation physics, and we study the flows of mass using Monte Carlo tracer particles. We estimate the efficiency of the inflow via the dust lanes, finding that only about a third (30$\pm$12%) of the dust lanes' mass initially accretes onto the CMZ, while the rest overshoots and accretes later. Given observational estimates of the amount of gas within the Milky Way's dust lanes, this suggests that the true total inflow rate onto the CMZ is 0.8$\pm$0.6 M$_\odot$~yr$^{-1}$. Clouds in this simulated CMZ have sudden peaks in their average density near the apocenter, where they undergo violent collisions with inflowing material. While these clouds tend to counter-rotate due to shear, co-rotating clouds occasionally occur due to the injection of momentum from collisions with inflowing material (~52% are strongly counter-rotating, and ~7% are strongly co-rotating of the 44 cloud sample). We investigate the formation and evolution of these clouds, finding that they are fed by many discrete inflow events, providing a consistent source of gas to CMZ clouds even as they collapse and form stars
△ Less
Submitted 9 November, 2021; v1 submitted 15 June, 2021;
originally announced June 2021.
-
SIGAME v3: Gas Fragmentation in Post-processing of Cosmological Simulations for More Accurate Infrared Line Emission Modeling
Authors:
Karen Pardos Olsen,
Blakesley Burkhart,
Mordecai-Mark Mac Low,
Robin G. Treß,
Thomas R. Greve,
David Vizgan,
Jay Motka,
Josh Borrow,
Gergö Popping,
Romeel Davé,
Rowan J. Smith,
Desika Narayanan
Abstract:
We present an update to the framework called SImulator of GAlaxy Millimeter/submillimeter Emission (SÍGAME). SÍGAME derives line emission in the far-infrared (FIR) for galaxies in particle-based cosmological hydrodynamics simulations by applying radiative transfer and physics recipes via a post-processing step after completion of the simulation. In this version, a new technique is developed to mod…
▽ More
We present an update to the framework called SImulator of GAlaxy Millimeter/submillimeter Emission (SÍGAME). SÍGAME derives line emission in the far-infrared (FIR) for galaxies in particle-based cosmological hydrodynamics simulations by applying radiative transfer and physics recipes via a post-processing step after completion of the simulation. In this version, a new technique is developed to model higher gas densities by parametrizing the gas density probability distribution function (PDF) in higher resolution simulations for use as a look-up table, allowing for more adaptive PDFs than in previous work. SÍGAME v3 is tested on redshift z = 0 galaxies drawn from the SIMBA cosmological simulation for eight FIR emission lines tracing vastly different interstellar medium phases. Including dust radiative transfer with SKIRT and high resolution photo-ionization models with Cloudy, this new method is able to self-consistently reproduce observed relations between line luminosity and star formation rate in all cases, except for [NII]122, [NII]205 and [OI]63, the luminosities of which are overestimated by median factors of 1.6, 1.2 and 1.2 dex, respectively. We attribute the remaining disagreement with observations to the lack of precise attenuation of the interstellar light on subgrid scales (<200 pc).
△ Less
Submitted 3 December, 2021; v1 submitted 4 February, 2021;
originally announced February 2021.
-
Simulations of the star-forming molecular gas in an interacting M51-like galaxy: cloud population statistics
Authors:
Robin G. Tress,
Mattia C. Sormani,
Rowan J. Smith,
Simon C. O. Glover,
Ralf S. Klessen,
Mordecai-Mark Mac Low,
Paul Clark,
Ana Duarte-Cabral
Abstract:
To investigate how molecular clouds react to different environmental conditions at a galactic scale, we present a catalogue of giant molecular clouds resolved down to masses of $\sim 10$~M$_{\odot}$ from a simulation of the entire disc of an interacting M51-like galaxy and a comparable isolated galaxy. Our model includes time-dependent gas chemistry, sink particles for star formation and supernova…
▽ More
To investigate how molecular clouds react to different environmental conditions at a galactic scale, we present a catalogue of giant molecular clouds resolved down to masses of $\sim 10$~M$_{\odot}$ from a simulation of the entire disc of an interacting M51-like galaxy and a comparable isolated galaxy. Our model includes time-dependent gas chemistry, sink particles for star formation and supernova feedback, meaning we are not reliant on star formation recipes based on threshold densities and can follow the physics of the cold molecular phase. We extract giant molecular clouds at a given timestep of the simulations and analyse their properties. In the disc of our simulated galaxies, spiral arms seem to act merely as snowplows, gathering gas and clouds without dramatically affecting their properties. In the centre of the galaxy, on the other hand, environmental conditions lead to larger, more massive clouds. While the galaxy interaction has little effect on cloud masses and sizes, it does promote the formation of counter-rotating clouds. We find that the identified clouds seem to be largely gravitationally unbound at first glance, but a closer analysis of the hierarchical structure of the molecular interstellar medium shows that there is a large range of virial parameters with a smooth transition from unbound to mostly bound for the densest structures. The common observation that clouds appear to be virialised entities may therefore be due to CO bright emission highlighting a specific level in this hierarchical binding sequence. The small fraction of gravitationally bound structures found suggests that low galactic star formation efficiencies may be set by the process of cloud formation and initial collapse.
△ Less
Submitted 28 June, 2021; v1 submitted 10 December, 2020;
originally announced December 2020.
-
The Cloud Factory II: Gravoturbulent Kinematics of Resolved Molecular Clouds in a Galactic Potential
Authors:
Andres F. Izquierdo,
Rowan J. Smith,
Simon C. O. Glover,
Ralf S. Klessen,
Robin G. Tress,
Mattia C. Sormani,
Paul C. Clark,
Ana Duarte-Cabral,
Catherine Zucker
Abstract:
We present a statistical analysis of the gravoturbulent velocity fluctuations in molecular cloud complexes extracted from our "Cloud Factory" galactic-scale ISM simulation suite. For this purpose, we produce non-LTE $^{12}$CO J=1-0 synthetic observations and apply the Principal Component Analysis (PCA) reduction technique on a representative sample of cloud complexes. The velocity fluctuations are…
▽ More
We present a statistical analysis of the gravoturbulent velocity fluctuations in molecular cloud complexes extracted from our "Cloud Factory" galactic-scale ISM simulation suite. For this purpose, we produce non-LTE $^{12}$CO J=1-0 synthetic observations and apply the Principal Component Analysis (PCA) reduction technique on a representative sample of cloud complexes. The velocity fluctuations are self-consistently generated by different physical mechanisms at play in our simulations, which include galactic-scale forces, gas self-gravity, and supernova feedback. The statistical analysis suggests that, even though purely gravitational effects are necessary to reproduce standard observational laws, they are not sufficient in most cases. We show that the extra injection of energy from supernova explosions plays a key role in establishing the global turbulent field and the local dynamics and morphology of molecular clouds. Additionally, we characterise structure function scaling parameters as a result of cloud environmental conditions: some of the complexes are immersed in diffuse (inter-arm) or dense (spiral-arm) environments, and others are influenced by embedded or external supernovae. In quiescent regions, we obtain time-evolving trajectories of scaling parameters driven by gravitational collapse and supersonic turbulent flows. Our findings suggests that a PCA-based statistical study is a robust method to diagnose the physical mechanisms that drive the gravoturbulent properties of molecular clouds. Also, we present a new open source module, the PCAFACTORY, which smartly performs PCA to extract velocity structure functions from simulated or real data of the ISM in a user-friendly way. Software DOI: 10.5281/zenodo.3822718
△ Less
Submitted 4 November, 2020;
originally announced November 2020.
-
Synthetic observations of spiral arm tracers of a simulated Milky Way analog
Authors:
Stefan Reissl,
Jeroen M. Stil,
En Chen,
Robin G. Treß,
Mattia C. Sormani,
Rowan J. Smith,
Ralf S. Klessen,
Megan Buick,
Simon C. O. Glover,
Russell Shanahan,
Stephen J. Lemmer,
Juan D. Soler,
Henrik Beuther,
James S. Urquhart,
L. D. Anderson,
Karl M. Menten,
Andreas Brunthaler,
Sarah Ragan,
Michael R. Rugel
Abstract:
Context: The Faraday rotation measure (RM) is often used to study the magnetic field strength and orientation within the ionized medium of the Milky Way. Observations indicate a RM in the spiral arms that exceeds the commonly assumed range. This raises the question of under what conditions spiral arms create such strong RM. Aims: We investigate the effect of spiral arms on Galactic RMs through sho…
▽ More
Context: The Faraday rotation measure (RM) is often used to study the magnetic field strength and orientation within the ionized medium of the Milky Way. Observations indicate a RM in the spiral arms that exceeds the commonly assumed range. This raises the question of under what conditions spiral arms create such strong RM. Aims: We investigate the effect of spiral arms on Galactic RMs through shock compression of the interstellar medium (ISM). It has recently been suggested that the Sagittarius spiral arm creates a strong peak in RM where the line of sight (LOS) is tangent to the arm, and that enhanced RM follows along an intersecting LOS. We seek to understand the physical conditions that give rise to this effect and the role of viewing geometry. Methods: We apply a MHD simulation of the multi-phase ISM in a Milky Way type spiral galaxy disk in combination with radiative transfer to evaluate different tracers of spiral arm structures. For observers embedded in the disk, dust intensity, synchrotron emission and the kinematics of molecular gas observations are derived to identify spiral arm tangents. RMs are calculated through the disk and evaluated for different observer positions. The observer's perspective is related to the parameters of the local bubble surrounding the observer. Results: We reproduce a scattering of tangent points for different tracers of about $6^\circ$ per spiral arm similar to the Milky Way. As for the RM, the model shows that compression of the ISM and associated amplification of the magnetic field in spiral arms enhances RM by a few hundred rad $m^{-2}$ on top of the mean contribution of the disk. The arm-inter-arm contrast in RM along the LOS is approximately 10 in the inner Galaxy, fading to ~2 in the outer Galaxy. We identify a shark-fin like pattern in the RM Milky Way observations as well as the synthetic data that is characteristic for spiral arms.
△ Less
Submitted 21 July, 2020;
originally announced July 2020.
-
Simulations of the Milky Way's central molecular zone -- II. Star formation
Authors:
Mattia C. Sormani,
Robin G. Tress,
Simon C. O. Glover,
Ralf S. Klessen,
Cara D. Battersby,
Paul C. Clark,
H Perry Hatchfield,
Rowan J. Smith
Abstract:
The Milky Way's central molecular zone (CMZ) has emerged in recent years as a unique laboratory for the study of star formation. Here we use the simulations presented in Tress et al. 2020 to investigate star formation in the CMZ. These simulations resolve the structure of the interstellar medium at sub-parsec resolution while also including the large-scale flow in which the CMZ is embedded. Our ma…
▽ More
The Milky Way's central molecular zone (CMZ) has emerged in recent years as a unique laboratory for the study of star formation. Here we use the simulations presented in Tress et al. 2020 to investigate star formation in the CMZ. These simulations resolve the structure of the interstellar medium at sub-parsec resolution while also including the large-scale flow in which the CMZ is embedded. Our main findings are as follows. (1) While most of the star formation happens in the CMZ ring at $R\gtrsim100 {\, \rm pc}$, a significant amount also occurs closer to SgrA* at $R \lesssim 10{\, \rm pc}$. (2) Most of the star formation in the CMZ happens downstream of the apocentres, consistent with the "pearls-on-a-string" scenario, and in contrast to the notion that an absolute evolutionary timeline of star formation is triggered by pericentre passage. (3) Within the timescale of our simulations ($\sim100$ Myr), the depletion time of the CMZ is constant within a factor of $\sim2$. This suggests that variations in the star formation rate are primarily driven by variations in the mass of the CMZ, caused for example by AGN feedback or externally-induced changes in the bar-driven inflow rate, and not by variations in the depletion time. (4) We study the trajectories of newly born stars in our simulations. We find several examples that have age and 3D velocity compatible with those of the Arches and Quintuplet clusters. Our simulations suggest that these prominent clusters originated near the collision sites where the bar-driven inflow accretes onto the CMZ, at symmetrical locations with respect to the Galactic centre, and that they have already decoupled from the gas in which they were born.
△ Less
Submitted 7 October, 2020; v1 submitted 14 April, 2020;
originally announced April 2020.
-
Simulations of the Milky Way's central molecular zone -- I. Gas dynamics
Authors:
Robin G. Tress,
Mattia C. Sormani,
Simon C. O. Glover,
Ralf S. Klessen,
Cara D. Battersby,
Paul C. Clark,
H Perry Hatchfield,
Rowan J. Smith
Abstract:
We use hydrodynamical simulations to study the Milky Way's central molecular zone (CMZ). The simulations include a non-equilibrium chemical network, the gas self-gravity, star formation and supernova feedback. We resolve the structure of the interstellar medium at sub-parsec resolution while also capturing the interaction between the CMZ and the bar-driven large-scale flow out to $R\sim 5\kpc$. Ou…
▽ More
We use hydrodynamical simulations to study the Milky Way's central molecular zone (CMZ). The simulations include a non-equilibrium chemical network, the gas self-gravity, star formation and supernova feedback. We resolve the structure of the interstellar medium at sub-parsec resolution while also capturing the interaction between the CMZ and the bar-driven large-scale flow out to $R\sim 5\kpc$. Our main findings are as follows: (1) The distinction between inner ($R\lesssim120$~pc) and outer ($120\lesssim R\lesssim450$~pc) CMZ that is sometimes proposed in the literature is unnecessary. Instead, the CMZ is best described as single structure, namely a star-forming ring with outer radius $R\simeq 200$~pc which includes the 1.3$^\circ$ complex and which is directly interacting with the dust lanes that mediate the bar-driven inflow. (2) This accretion can induce a significant tilt of the CMZ out of the plane. A tilted CMZ might provide an alternative explanation to the $\infty$-shaped structure identified in Herschel data by Molinari et al. 2011. (3) The bar in our simulation efficiently drives an inflow from the Galactic disc ($R\simeq 3$~kpc) down to the CMZ ($R\simeq200$~pc) of the order of $1\rm\,M_\odot\,yr^{-1}$, consistent with observational determinations. (4) Supernova feedback can drive an inflow from the CMZ inwards towards the circumnuclear disc of the order of $\sim0.03\,\rm M_\odot\,yr^{-1}$. (5) We give a new interpretation for the 3D placement of the 20 and 50 km s$^{-1}$ clouds, according to which they are close ($R\lesssim30$~pc) to the Galactic centre, but are also connected to the larger-scale streams at $R\gtrsim100$~pc.
△ Less
Submitted 7 October, 2020; v1 submitted 14 April, 2020;
originally announced April 2020.
-
The Cloud Factory I: Generating resolved filamentary molecular clouds from galactic-scale forces
Authors:
Rowan J. Smith,
Robin G. Treß,
Mattia C. Sormani,
Simon C. O. Glover,
Ralf S. Klessen,
Paul C. Clark,
Andrés F. Izquierdo,
Ana Duarte Cabral,
Catherine Zucker
Abstract:
We introduce a new suite of simulations, "The Cloud Factory", which self-consistently forms molecular cloud complexes at high enough resolution to resolve internal substructure (up to 0.25 Msol in mass) all while including galactic-scale forces. We use a version of the Arepo code modified to include a detailed treatment of the physics of the cold molecular ISM, and an analytical galactic gravitati…
▽ More
We introduce a new suite of simulations, "The Cloud Factory", which self-consistently forms molecular cloud complexes at high enough resolution to resolve internal substructure (up to 0.25 Msol in mass) all while including galactic-scale forces. We use a version of the Arepo code modified to include a detailed treatment of the physics of the cold molecular ISM, and an analytical galactic gravitational potential for computational efficiency. The simulations have nested levels of resolution, with the lowest layer tied to tracer particles injected into individual cloud complexes. These tracer refinement regions are embedded in the larger simulation so continue to experience forces from outside the cloud. This allows the simulations to act as a laboratory for testing the effect of galactic environment on star formation. Here we introduce our method and investigate the effect of galactic environment on filamentary clouds. We find that cloud complexes formed after a clustered burst of feedback, have shorter lengths and are less likely to fragment compared to quiescent clouds (e.g. the Musca filament) or those dominated by the galactic potential (e.g. Nessie). Spiral arms and differential rotation preferentially align filaments, but strong feedback randomises them. Long filaments formed within the cloud complexes are necessarily coherent with low internal velocity gradients, which has implications for the formation of filamentary star-clusters. Cloud complexes formed in regions dominated by supernova feedback have fewer star-forming cores, and these are more widely distributed. These differences show galactic-scale forces can have a significant impact on star formation within molecular clouds.
△ Less
Submitted 18 November, 2019; v1 submitted 13 November, 2019;
originally announced November 2019.
-
Simulations of the star-forming molecular gas in an interacting M51-like galaxy
Authors:
Robin G. Tress,
Rowan J. Smith,
Mattia C. Sormani,
Simon C. O. Glover,
Ralf S. Klessen,
Mordecai-Mark Mac Low,
Paul C. Clark
Abstract:
We present here the first of a series of papers aimed at better understanding the evolution and properties of giant molecular clouds (GMCs) in a galactic context. We perform high resolution, three-dimensional {\sc arepo} simulations of an interacting galaxy inspired by the well-observed M51 galaxy. Our fiducial simulations include a non-equilibrium, time-dependent, chemical network that follows th…
▽ More
We present here the first of a series of papers aimed at better understanding the evolution and properties of giant molecular clouds (GMCs) in a galactic context. We perform high resolution, three-dimensional {\sc arepo} simulations of an interacting galaxy inspired by the well-observed M51 galaxy. Our fiducial simulations include a non-equilibrium, time-dependent, chemical network that follows the evolution of atomic and molecular hydrogen as well as carbon and oxygen self-consistently. Our calculations also treat gas self-gravity and subsequent star formation (described by sink particles), and coupled supernova feedback. In the densest parts of the simulated interstellar medium (ISM) we reach sub-parsec resolution, granting us the ability to resolve individual GMCs and their formation and destruction self-consistently throughout the galaxy. In this initial work we focus on the general properties of the ISM with a particular focus on the cold star-forming gas. We discuss the role of the interaction with the companion galaxy in generating cold molecular gas and controlling stellar birth. We find that while the interaction drives large-scale gas flows and induces spiral arms in the galaxy, it is of secondary importance in determining gas fractions in the different ISM phases and the overall star-formation rate. The behaviour of the gas on small GMC scales instead is mostly controlled by the self-regulating property of the ISM driven by coupled feedback.
△ Less
Submitted 23 September, 2019;
originally announced September 2019.
-
The geometry of the gas surrounding the Central Molecular Zone: on the origin of localised molecular clouds with extreme velocity dispersions
Authors:
Mattia C. Sormani,
Robin G. Tress,
Simon C. O. Glover,
Ralf S. Klessen,
Ashley T. Barnes,
Cara D. Battersby,
Paul C. Clark,
H Perry Hatchfield,
Rowan J. Smith
Abstract:
Observations of molecular gas near the Galactic centre ($| l | < 10^\circ$, $| b | < 1^\circ$) reveal the presence of a distinct population of enigmatic compact clouds which are characterised by extreme velocity dispersions ($Δv > 100\, \rm km/s$). These Extended Velocity Features (EVFs) are very prominent in the datacubes and dominate the kinematics of molecular gas just outside the Central Molec…
▽ More
Observations of molecular gas near the Galactic centre ($| l | < 10^\circ$, $| b | < 1^\circ$) reveal the presence of a distinct population of enigmatic compact clouds which are characterised by extreme velocity dispersions ($Δv > 100\, \rm km/s$). These Extended Velocity Features (EVFs) are very prominent in the datacubes and dominate the kinematics of molecular gas just outside the Central Molecular Zone (CMZ). The prototypical example of such a cloud is Bania Clump 2. We show that similar features are naturally produced in simulations of gas flow in a realistic barred potential. We analyse the structure of the features obtained in the simulations and use this to interpret the observations. We find that the features arise from collisions between material that has been infalling rapidly along the dust lanes of the Milky Way bar and material that belongs to one of the following two categories: (i) material that has `overshot' after falling down the dust lanes on the opposite side; (ii) material which is part of the CMZ. Both types of collisions involve gas with large differences in the line-of-sight velocities, which is what produces the observed extreme velocity dispersions. Examples of both categories can be identified in the observations. If our interpretation is correct, we are directly witnessing (a) collisions of clouds with relative speeds of $\sim 200\, \rm km/s$ and (b) the process of accretion of fresh gas onto the CMZ.
△ Less
Submitted 24 June, 2019;
originally announced June 2019.
-
A dynamical mechanism for the origin of nuclear rings
Authors:
Mattia C. Sormani,
Emanuele Sobacchi,
Francesca Fragkoudi,
Matthew Ridley,
Robin G. Tress,
Simon C. O. Glover,
Ralf S. Klessen
Abstract:
We develop a dynamical theory for the origin of nuclear rings in barred galaxies. In analogy with the standard theory of accretion discs, our theory is based on shear viscous forces among nested annuli of gas. However, the fact that gas follows non circular orbits in an external barred potential has profound consequences: it creates a region of reverse shear in which it is energetically favourable…
▽ More
We develop a dynamical theory for the origin of nuclear rings in barred galaxies. In analogy with the standard theory of accretion discs, our theory is based on shear viscous forces among nested annuli of gas. However, the fact that gas follows non circular orbits in an external barred potential has profound consequences: it creates a region of reverse shear in which it is energetically favourable to form a stable ring which does not spread despite dissipation. Our theory allows us to approximately predict the size of the ring given the underlying gravitational potential. The size of the ring is loosely related to the location of the Inner Lindblad Resonance in the epicyclic approximation, but the predicted location is more accurate and is also valid for strongly barred potentials. By comparing analytical predictions with the results of hydrodynamical simulations, we find that our theory provides a viable mechanism for ring formation if the effective sound speed of the gas is low ($\cs\lesssim1\kms$), but that nuclear spirals/shocks created by pressure destroy the ring when the sound speed is high ($\cs\simeq10\kms$). We conclude that whether this mechanism for ring formation is relevant for real galaxies ultimately depends on the effective equation of state of the ISM. Promising confirmation comes from simulations in which the ISM is modelled using state-of-the-art cooling functions coupled to live chemical networks, but more tests are needed regarding the role of turbulence driven by stellar feedback. If the mechanism is relevant in real galaxies, it could provide a powerful tool to constrain the gravitational potential, in particular the bar pattern speed.
△ Less
Submitted 14 August, 2018; v1 submitted 21 May, 2018;
originally announced May 2018.
-
A theoretical explanation for the Central Molecular Zone asymmetry
Authors:
Mattia C. Sormani,
Robin G. Tress,
Matthew Ridley,
Simon C. O. Glover,
Ralf S. Klessen,
James Binney,
John Magorrian,
Rowan Smith
Abstract:
It has been known for more than thirty years that the distribution of molecular gas in the innermost 300 parsecs of the Milky Way, the Central Molecular Zone, is strongly asymmetric. Indeed, approximately three quarters of molecular emission comes from positive longitudes, and only one quarter from negative longitudes. However, despite much theoretical effort, the origin of this asymmetry has rema…
▽ More
It has been known for more than thirty years that the distribution of molecular gas in the innermost 300 parsecs of the Milky Way, the Central Molecular Zone, is strongly asymmetric. Indeed, approximately three quarters of molecular emission comes from positive longitudes, and only one quarter from negative longitudes. However, despite much theoretical effort, the origin of this asymmetry has remained a mystery. Here we show that the asymmetry can be neatly explained by unsteady flow of gas in a barred potential. We use high-resolution 3D hydrodynamical simulations coupled to a state-of-the-art chemical network. Despite the initial conditions and the bar potential being point-symmetric with respect to the Galactic Centre, asymmetries develop spontaneously due to the combination of a hydrodynamical instability known as the "wiggle instability" and the thermal instability. The observed asymmetry must be transient: observations made tens of megayears in the past or in the future would often show an asymmetry in the opposite sense. Fluctuations of amplitude comparable to the observed asymmetry occur for a large fraction of the time in our simulations, and suggest that the present is not an exceptional moment in the life of our Galaxy.
△ Less
Submitted 15 December, 2017; v1 submitted 12 July, 2017;
originally announced July 2017.
-
Periodicity makes galactic shocks unstable - I. Linear analysis
Authors:
Mattia C. Sormani,
Emanuele Sobacchi,
Steven N. Shore,
Robin G. Tress,
Ralf S. Klessen
Abstract:
We study the dynamical stability of stationary galactic spiral shocks. The steady-state equilibrium flow contains a shock of the type derived by Roberts in the tightly wound approximation. We find that boundary conditions are critical in determining whether the solutions are stable or not. Shocks are unstable if periodic boundary conditions are imposed. For intermediate strengths of the spiral pot…
▽ More
We study the dynamical stability of stationary galactic spiral shocks. The steady-state equilibrium flow contains a shock of the type derived by Roberts in the tightly wound approximation. We find that boundary conditions are critical in determining whether the solutions are stable or not. Shocks are unstable if periodic boundary conditions are imposed. For intermediate strengths of the spiral potential, the instability disappears if boundary conditions are imposed such that the upstream flow is left unperturbed as in the classic analysis of D'yakov and Kontorovich. This reconciles apparently contradictory findings of previous authors regarding the stability of spiral shocks. This also shows that the instability is distinct from the Kelvin-Helmholtz instability, confirming the findings of Kim et al. We suggest that instability is a general characteristics of periodic shocks, regardless of the presence of shear, and provide a physical picture as to why this is the case. For strong spiral potentials, high post-shock shear makes the system unstable also to parasitic Kelvin-Helmholtz instability regardless of the boundary conditions. Our analysis is performed in the context of a simplified problem that, while preserving all the important characteristics of the original problem, strips it from unnecessary complications, and assumes that the gas is isothermal, non self-gravitating, non-magnetised.
△ Less
Submitted 2 July, 2017;
originally announced July 2017.
-
Nuclear Spirals in the inner Milky Way
Authors:
Matthew Ridley,
Mattia C. Sormani,
Robin G. Treß,
John Magorrian,
Ralf S. Klessen
Abstract:
We use hydrodynamical simulations to construct a new coherent picture for the gas flow in the Central Molecular Zone (CMZ), the region of our Galaxy within $R\leq 500\, \mathrm{pc}$. We relate connected structures observed in $(l,b,v)$ data cubes of molecular tracers to nuclear spiral arms. These arise naturally in hydrodynamical simulations of barred galaxies, and are similar to those that can be…
▽ More
We use hydrodynamical simulations to construct a new coherent picture for the gas flow in the Central Molecular Zone (CMZ), the region of our Galaxy within $R\leq 500\, \mathrm{pc}$. We relate connected structures observed in $(l,b,v)$ data cubes of molecular tracers to nuclear spiral arms. These arise naturally in hydrodynamical simulations of barred galaxies, and are similar to those that can be seen in external galaxies such as NGC4303 or NGC1097. We discuss a face-on view of the CMZ including the position of several prominent molecular clouds, such as Sgr B2, the $20\,{\rm km\, s^{-1}}$ and $50\,{\rm km\, s^{-1}}$ clouds, the polar arc, Bania Clump 2 and Sgr C. Our model is also consistent with the larger scale gas flow, up to $R\simeq 3\,\rm kpc$, thus providing a consistent picture of the entire Galactic bar region.
△ Less
Submitted 21 April, 2017; v1 submitted 12 April, 2017;
originally announced April 2017.
-
A simple method to convert sink particles into stars
Authors:
Mattia C. Sormani,
Robin G. Treß,
Ralf S. Klessen,
Simon C. O. Glover
Abstract:
Hydrodynamical simulations of star formation often do not possess the dynamic range needed to fully resolve the build-up of individual stars and star clusters, and thus have to resort to subgrid models. A popular way to do this is by introducing Lagrangian sink particles, which replace contracting high density regions at the point where the resolution limit is reached. A common problem then is how…
▽ More
Hydrodynamical simulations of star formation often do not possess the dynamic range needed to fully resolve the build-up of individual stars and star clusters, and thus have to resort to subgrid models. A popular way to do this is by introducing Lagrangian sink particles, which replace contracting high density regions at the point where the resolution limit is reached. A common problem then is how to assign fundamental stellar properties to sink particles, such as the distribution of stellar masses.
We present a new and simple statistical method to assign stellar contents to sink particles. Once the stellar content is specified, it can be used to determine a sink particle's radiative output, supernovae rate or other feedback parameters that may be required in the calculations. Advantages of our method are (i) it is simple to implement, (ii) it guarantees that the obtained stellar populations are good samples of the initial mass function, (iii) it can easily deal with infalling mass accreted at later times, and (iv) it does not put restrictions on the sink particles' masses in order to be used. The method works very well for sink particles that represent large star clusters and for which the stellar mass function is well sampled, but can also handle the transition to sink particles that represent a small number of stars.
△ Less
Submitted 7 December, 2016; v1 submitted 8 October, 2016;
originally announced October 2016.
