Showing 1–46 of 46 results for author: di Teodoro, E

Turbulence statistics of HI clouds entrained in the Milky Way's nuclear wind

Authors: Isabella A. Gerrard, Karlie A. Noon, Christoph Federrath, Enrico M. Di Teodoro, Antoine Marchal, N. M. McClure-Griffiths

Abstract: The interstellar medium (ISM) is ubiquitously turbulent across many physically distinct environments within the Galaxy. Turbulence is key in controlling the structure and dynamics of the ISM, regulating star formation, and transporting metals within the Galaxy. We present the first observational measurements of turbulence in neutral hydrogen entrained in the hot nuclear wind of the Milky Way. Usin… ▽ More The interstellar medium (ISM) is ubiquitously turbulent across many physically distinct environments within the Galaxy. Turbulence is key in controlling the structure and dynamics of the ISM, regulating star formation, and transporting metals within the Galaxy. We present the first observational measurements of turbulence in neutral hydrogen entrained in the hot nuclear wind of the Milky Way. Using recent MeerKAT observations of two extra-planar HI clouds above (gal. lat.$\,\sim7.0^{\circ}$) and below (gal. lat.$\,\sim-3.9^{\circ}$) the Galactic disc, we analyse centroid velocity and column density maps to estimate the velocity dispersion ($σ_{v,\mathrm{3D}}$), the turbulent sonic Mach number ($\mathcal{M}$), the volume density dispersion ($σ_{ρ/ρ_0}$), and the turbulence driving parameter ($b$). We also present a new prescription for estimating the spatial temperature variations of HI in the presence of related molecular gas. We measure these turbulence quantities on the global scale of each cloud, but also spatially map their variation across the plane-of-sky extent of each cloud by using a roving kernel method. We find that the two clouds share very similar characteristics of their internal turbulence, despite their varying latitudes. Both clouds are in the sub-to-trans-sonic Mach regime, and have primarily compressively-driven ($b\sim1$) turbulence. Given that there is no known active star-formation present in these clouds, this may be indicative of the way the cloud-wind interaction injects energy into the entrained atomic material on parsec scales. △ Less

Submitted 28 April, 2024; originally announced April 2024.

Comments: Accepted to MNRAS, 15 pages, 13 figures

Feedback and ionized gas outflows in four low-radio power AGN at z $\sim$0.15

Authors: L. Ulivi, G. Venturi, G. Cresci, A. Marconi, C. Marconcini, A. Amiri, F. Belfiore, E. Bertola, S. Carniani, Q. D Amato, E. Di Teodoro, M. Ginolfi, A. Girdhar, C. Harrison, R. Maiolino, F. Mannucci, M. Mingozzi, M. Perna, M. Scialpi, N. Tomicic, G. Tozzi, E. Treister

Abstract: An increasing number of observations and simulations suggests that low-power (<10$^{44}$ erg s$^{-1}$) jets may be a significant channel of feedback produced by active galactic nuclei (AGN), but little is known about their actual effect on their host galaxies from the observational point of view. We targeted four luminous type 2 AGN hosting moderately powerful radio emission ($\sim$10$^{44}$ erg s… ▽ More An increasing number of observations and simulations suggests that low-power (<10$^{44}$ erg s$^{-1}$) jets may be a significant channel of feedback produced by active galactic nuclei (AGN), but little is known about their actual effect on their host galaxies from the observational point of view. We targeted four luminous type 2 AGN hosting moderately powerful radio emission ($\sim$10$^{44}$ erg s$^{-1}$), two of which and possibly a third are associated with jets, with optical integral field spectroscopy observations from the Multi Unit Spectroscopic Explorer (MUSE) at the Very Large Telescope (VLT) to analyze the properties of their ionized gas as well as the properties and effects of ionized outflows. We combined these observations with Very Large Array (VLA) and e-MERLIN data to investigate the relations and interactions between the radio jets and host galaxies. We detected ionized outflows as traced by the fast bulk motion of the gas. The outflows extended over kiloparsec scales in the direction of the jet, when present. In the two sources with resolved radio jets, we detected a strong enhancement in the emission-line velocity dispersion (up to 1000 km s$^{-1}$) perpendicular to the direction of the radio jets. We also found a correlation between the mass and the energetics of this high-velocity dispersion gas and the radio power, which supports the idea that the radio emission may cause the enhanced turbulence. This phenomenon, which is now being observed in an increasing number of objects, might represent an important channel for AGN feedback on galaxies. △ Less

Submitted 2 March, 2024; originally announced March 2024.

Comments: 18 pages, 17 figures

Journal ref: A&A 685, A122 (2024)

Braving the Storm: Quantifying Disk-wide Ionized Outflows in the Large Magellanic Cloud with ULLYSES

Authors: Yong Zheng, Kirill Tchernyshyov, Knut Olsen, Yumi Choi, Chad Bustard, Julia Roman-Duval, Robert Zhu, Enrico M. Di Teodoro, Jessica Werk, Mary Putman, Anna F. McLeod, Yakov Faerman, Raymond C. Simons, Joshua Peek

Abstract: The Large Magellanic Cloud (LMC) is home to many HII regions, which may lead to significant outflows. We examine the LMC's multiphase gas ($T\sim10^{4-5}$ K) in HI, SII, SiIV, and CIV using 110 stellar sight lines from the HST's Ultraviolet Legacy Library of Young Stars as Essential Standards (ULLYSES) program. We develop a continuum fitting algorithm based on the concept of Gaussian Process regre… ▽ More The Large Magellanic Cloud (LMC) is home to many HII regions, which may lead to significant outflows. We examine the LMC's multiphase gas ($T\sim10^{4-5}$ K) in HI, SII, SiIV, and CIV using 110 stellar sight lines from the HST's Ultraviolet Legacy Library of Young Stars as Essential Standards (ULLYSES) program. We develop a continuum fitting algorithm based on the concept of Gaussian Process regression and identify reliable LMC interstellar absorption over $v_{\rm helio}=175-375$ km s$^{-1}$. Our analyses show disk-wide ionized outflows in SiIV and CIV across the LMC with bulk velocities of $|v_{\rm out, bulk}|\sim20-60$ km s$^{-1}$, which indicates that most of the outflowing mass is gravitationally bound. The outflows' column densities correlate with the LMC's star formation rate surface densities ($Σ_{\rm SFR}$), and the outflows with higher $Σ_{\rm SFR}$ tend to be more ionized. Considering outflows from both sides of the LMC as traced by CIV, we conservatively estimate a total outflow rate of $\dot{M}_{\rm out}\gtrsim 0.03~M_\odot {\rm yr}^{-1}$ and a mass loading factor of $η\gtrsim 0.15$. We compare the LMC's outflows with those detected in starburst galaxies and simulation predictions, and find a universal scaling relation of $|v_{\rm out, bulk}|\propto Σ_{\rm SFR}^{0.23}$ over a wide range of star-forming conditions ($Σ_{\rm SFR}\sim10^{-4.5}-10^{2}~M_\odot {\rm yr}^{-1} {\rm kpc}^{-2}$). Lastly, we find that the outflows are co-rotating with the LMC's young stellar disk and the velocity field does not seem to be significantly impacted by external forces; we thus speculate on the existence of a bow shock leading the LMC, which may have shielded the outflows from ram pressure as the LMC orbits the Milky Way. △ Less

Submitted 4 August, 2024; v1 submitted 6 February, 2024; originally announced February 2024.

Comments: Accepted for publication in ApJ. Key findings can be found in Figures 9-12. Main updates from the previous version include new estimates on star formation rate surface densities. Normalized SII, SiIV, and CIV line spectra derived for this work are published as a High Level Science Product called LMC-FLOWS (doi: 10.17909/hz0m-np43), available on website: https://archive.stsci.edu/hlsp/lmc-flows

HOMERUN a new approach to photoionization modelling. I -- reproducing observed emission lines with percent accuracy and obtaining accurate physical properties of the ionized gas

Authors: A. Marconi, A. Amiri, A. Feltre, F. Belfiore, G. Cresci, M. Curti, F. Mannucci, E. Bertola, M. Brazzini, S. Carniani, E. Cataldi, Q. D'Amato, G. de Rosa, E. Di Teodoro, M. Ginolfi, N. Kumari, C. Marconcini, R. Maiolino, L. Magrini, A. Marasco, M. Mingozzi, B. Moreschini, T. Nagao, E. Oliva, M. Scialpi , et al. (4 additional authors not shown)

Abstract: We present HOMERUN (Highly Optimized Multi-cloud Emission-line Ratios Using photo-ionizatioN), a new approach to modelling emission lines from photoionized gas that can simultaneously reproduce all observed line intensities from a wide range of ionization levels and with high accuracy. Our approach is based on the weighted combination of multiple single-cloud photoionization models and, contrary t… ▽ More We present HOMERUN (Highly Optimized Multi-cloud Emission-line Ratios Using photo-ionizatioN), a new approach to modelling emission lines from photoionized gas that can simultaneously reproduce all observed line intensities from a wide range of ionization levels and with high accuracy. Our approach is based on the weighted combination of multiple single-cloud photoionization models and, contrary to previous works, the novelty of our approach consists in using the weights as free parameters of the fit and constraining them with the observed data. One of the main applications of HOMERUN is the accurate determination of gas-phase metallicities and we show that a critical point is to allow for a variation of the N/O and S/O abundance ratios which can significantly improve the quality of the fit and the accuracy of the results. Moreover, our approach provides a major improvement compared to the single-cloud, constant-pressure models commonly used in the literature. By using high-quality literature spectra of H ii regions where 10 to 20 emission lines (including several auroral lines) are detected with high signal-to-noise ratio, we show that all lines are reproduced by the model with an accuracy better than 10%. In particular, the model is able to simultaneously reproduce [O i]6300, 6363, [O ii]3726, 3729, [O iii]4959, 5007, [S ii]6717, 6731, and [S iii]9069, 9532 emission lines which, to our knowledge, is an unprecedented result. Finally, we show that the gas metallicities estimated with our models for HII regions in the Milky Way are in agreement with the stellar metallicities than the estimates based on the Te-method. Overall, our method provides a new accurate tool to estimate the metallicity and the physical conditions of the ionized gas. It can be applied to many different science cases from HII regions to AGN and wherever there are emission lines from photoionized gas. △ Less

Submitted 26 June, 2024; v1 submitted 23 January, 2024; originally announced January 2024.

Comments: A&A in press, new version following referee report

Journal ref: A&A 689, A78 (2024)

A Galactic Eclipse: The Small Magellanic Cloud is Forming Stars in Two, Superimposed Systems

Authors: Claire E. Murray, Sten Hasselquist, Joshua E. G. Peek, Christina Willecke Lindberg, Andres Almeida, Yumi Choi, Jessica E. M. Craig, Helga Denes, John M. Dickey, Enrico M. Di Teodoro, Christoph Federrath, Isabella A. Gerrard, Steven J. Gibson, Denis Leahy, Min-Young Lee, Callum Lynn, Yik Ki Ma, Antoine Marchal, N. M. McClure-Griffiths, David Nidever, Hiep Nguyen, Nickolas M. Pingel, Elizabeth Tarantino, Lucero Uscanga, Jacco Th. van Loon

Abstract: The structure and dynamics of the star-forming disk of the Small Magellanic Cloud (SMC) have long confounded us. The SMC is widely used as a prototype for galactic physics at low metallicity, and yet we fundamentally lack an understanding of the structure of its interstellar medium (ISM). In this work, we present a new model for the SMC by comparing the kinematics of young, massive stars with the… ▽ More The structure and dynamics of the star-forming disk of the Small Magellanic Cloud (SMC) have long confounded us. The SMC is widely used as a prototype for galactic physics at low metallicity, and yet we fundamentally lack an understanding of the structure of its interstellar medium (ISM). In this work, we present a new model for the SMC by comparing the kinematics of young, massive stars with the structure of the ISM traced by high-resolution observations of neutral atomic hydrogen (HI) from the Galactic Australian Square Kilometer Array Pathfinder survey (GASKAP-HI). Specifically, we identify thousands of young, massive stars with precise radial velocity constraints from the Gaia and APOGEE surveys and match these stars to the ISM structures in which they likely formed. By comparing the average dust extinction towards these stars, we find evidence that the SMC is composed of two structures with distinct stellar and gaseous chemical compositions. We construct a simple model that successfully reproduces the observations and shows that the ISM of the SMC is arranged into two, superimposed, star-forming systems with similar gas mass separated by ~5 kpc along the line of sight. △ Less

Submitted 12 December, 2023; originally announced December 2023.

Comments: ApJ accepted. 20 pages, 18 figures

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.

Comments: This White Paper will be updated when required (e.g. new authors joining, editing of content). Most recent update: 24 Oct 2023

MOKA3D: An innovative approach to 3D gas kinematic modelling. I. Application to AGN ionized outflows

Authors: C. Marconcini, A. Marconi, G. Cresci, G. Venturi, L. Ulivi, F. Mannucci, F. Belfiore, G. Tozzi, M. Ginolfi, A. Marasco, S. Carniani, A. Amiri, E. Di Teodoro, M. Scialpi, N. Tomicic, M. Mingozzi, M. Brazzini, B. Moreschini

Abstract: Studying the feedback process of Active Galactic Nuclei (AGN) requires characterising multiple kinematical components, such as rotating gas and stellar disks, outflows, inflows, and jets. To compare the observed properties with theoretical predictions of galaxy evolution and feedback models and to assess the mutual interaction and energy injection rate into the interstellar medium (ISM), one usual… ▽ More Studying the feedback process of Active Galactic Nuclei (AGN) requires characterising multiple kinematical components, such as rotating gas and stellar disks, outflows, inflows, and jets. To compare the observed properties with theoretical predictions of galaxy evolution and feedback models and to assess the mutual interaction and energy injection rate into the interstellar medium (ISM), one usually relies on simplified kinematic models. These models have several limitations, as they often do not take into account projection effects, beam smearing and the surface brightness distribution of the emitting medium. Here, we present MOKA3D, an innovative approach to model the 3D gas kinematics from integral field spectroscopy observations. In this first paper, we discuss its application to the case of AGN ionised outflows, whose observed clumpy emission and apparently irregular kinematics are only marginally accounted for by existing kinematical models. Unlike previous works, our model does not assume the surface brightness distribution of the gas, but exploits a novel procedure to derive it from the observations by reconstructing the 3D distribution of emitting clouds and providing accurate estimates of the spatially resolved outflow physical properties (e.g. mass rate, kinetic energy). As an example, we demonstrate the capabilities of our method by applying it to three nearby Seyfert-II galaxies observed with MUSE at the VLT and selected from the MAGNUM survey, showing that the complex kinematic features observed can be described by a conical outflow with a constant radial velocity field and a clumpy distribution of clouds. △ Less

Submitted 4 July, 2023; originally announced July 2023.

Comments: 17 pages, 14 figures

Journal ref: A&A 677, A58 (2023)

Direct observations of the atomic-molecular phase transition in the Milky Way's nuclear wind

Authors: Karlie A. Noon, Mark R. Krumholz, Enrico M. Di Teodoro, Naomi M. McClure-Griffiths, Felix J. Lockman, Lucia Armillotta

Abstract: Hundreds of high-velocity atomic gas clouds exist above and below the Galactic Centre, with some containing a molecular component. However, the origin of these clouds in the Milky Way's wind is unclear. This paper presents new high-resolution MeerKAT observations of three atomic gas clouds and studies the relationship between the atomic and molecular phases at $\sim 1$ pc scales. The clouds' atomi… ▽ More Hundreds of high-velocity atomic gas clouds exist above and below the Galactic Centre, with some containing a molecular component. However, the origin of these clouds in the Milky Way's wind is unclear. This paper presents new high-resolution MeerKAT observations of three atomic gas clouds and studies the relationship between the atomic and molecular phases at $\sim 1$ pc scales. The clouds' atomic hydrogen column densities, $N_{\mathrm{HI}}$, are less than a $\mbox{few}\times 10^{20}$ cm$^{-2}$, but the two clouds closest to the Galactic Centre nonetheless have detectable CO emission. This implies the presence of H$_{2}$ at levels of $N_{\mathrm{HI}}$ at least a factor of ten lower than in the typical Galactic interstellar medium. For the cloud closest to the Galactic Centre, there is little correlation between the $N_{\mathrm{HI}}$ and the probability that it will harbour detectable CO emissions. In contrast, for the intermediate cloud, detectable CO is heavily biased toward the highest values of $N_{\mathrm{HI}}$. The cloud most distant from the Galactic Centre has no detectable CO at similar $N_{\mathrm{HI}}$ values. Moreover, we find that the two clouds with detectable CO are too molecule-rich to be in chemical equilibrium, given the depths of their atomic shielding layers, which suggests a scenario whereby these clouds consist of pre-existing molecular gas from the disc that the Galactic wind has swept up, and that is dissociating into atomic hydrogen as it flows away from the Galaxy. We estimate that entrained molecular material of this type has a $\sim \mathrm{few}-10$ Myr lifetime before photodissociating. △ Less

Submitted 13 July, 2023; v1 submitted 13 April, 2023; originally announced April 2023.

Comments: 11 pages, 6 figures, 2 tables. Accepted for publication in MNRAS

Journal ref: Volume 524, Issue 1, 2023, Pages 1258 - 1268

Cold gas disks in main-sequence galaxies at cosmic noon: Low turbulence, flat rotation curves, and disk-halo degeneracy

Authors: Federico Lelli, Zhi-Yu Zhang, Thomas G. Bisbas, Lingrui Lin, Padelis Papadopoulos, James M. Schombert, Enrico Di Teodoro, Antonino Marasco, Stacy S. McGaugh

Abstract: We study the dynamics of cold molecular gas in two main-sequence galaxies at cosmic noon (zC-488879 at $z\simeq1.47$ and zC-400569 at $z\simeq2.24$) using new high-resolution ALMA observations of multiple $^{12}$CO transitions. For zC-400569 we also re-analyze high-quality H$α$ data from the SINS/zC-SINF survey. We find that (1) Both galaxies have regularly rotating CO disks and their rotation cur… ▽ More We study the dynamics of cold molecular gas in two main-sequence galaxies at cosmic noon (zC-488879 at $z\simeq1.47$ and zC-400569 at $z\simeq2.24$) using new high-resolution ALMA observations of multiple $^{12}$CO transitions. For zC-400569 we also re-analyze high-quality H$α$ data from the SINS/zC-SINF survey. We find that (1) Both galaxies have regularly rotating CO disks and their rotation curves are flat out to $\sim$8 kpc contrary to previous results pointing to outer declines in the rotation speed $V_{\rm rot}$; (2) The intrinsic velocity dispersions are low ($σ_{\rm CO}\lesssim15$ km/s for CO and $σ_{\rm Hα}\lesssim37$ km/s for H$α$) and imply $V_{\rm rot}/σ_{\rm CO}\gtrsim17-22$ yielding no significant pressure support; (3) Mass models using HST images display a severe disk-halo degeneracy: models with inner baryon dominance and models with "cuspy" dark matter halos can fit the rotation curves equally well due to the uncertainties on stellar and gas masses; (4) Milgromian dynamics (MOND) can successfully fit the rotation curves with the same acceleration scale $a_0$ measured at $z\simeq0$. The question of the amount and distribution of dark matter in high-$z$ galaxies remains unsettled due to the limited spatial extent of the available kinematic data; we discuss the suitability of various emission lines to trace extended rotation curves at high $z$. Nevertheless, the properties of these two high-$z$ galaxies (high $V_{\rm rot}/σ_{\rm V}$ ratios, inner rotation curve shapes, bulge-to-total mass ratios) are remarkably similar to those of massive spirals at $z\simeq0$, suggesting weak dynamical evolution over more than 10 Gyr of the Universe's lifetime. △ Less

Submitted 25 May, 2023; v1 submitted 31 January, 2023; originally announced February 2023.

Comments: 18 pages, 11 figures, 4 tables, 2 appendices. Minor revisions to match the published version in Astronomy and Astrophysics

Journal ref: A&A 672, A106 (2023)

Identification of galaxy shreds in large photometric catalogs using Convolutional Neural Networks

Authors: Enrico M. Di Teodoro, Josh E. G. Peek, John F. Wu

Abstract: Contamination from galaxy fragments, identified as sources, is a major issue in large photometric galaxy catalogs. In this paper, we prove that this problem can be easily addressed with computer vision techniques. We use image cutouts to train a convolutional neural network (CNN) to identify catalogued sources that are in reality just star formation regions and/or shreds of larger galaxies. The CN… ▽ More Contamination from galaxy fragments, identified as sources, is a major issue in large photometric galaxy catalogs. In this paper, we prove that this problem can be easily addressed with computer vision techniques. We use image cutouts to train a convolutional neural network (CNN) to identify catalogued sources that are in reality just star formation regions and/or shreds of larger galaxies. The CNN reaches an accuracy ~98% on our testing datasets. We apply this CNN to galaxy catalogs from three amongst the largest surveys available today: the Sloan Digital Sky Survey (SDSS), the DESI Legacy Imaging Surveys and the Panoramic Survey Telescope and Rapid Response System Survey (Pan-STARSS). We find that, even when strict selection criteria are used, all catalogs still show a ~5% level of contamination from galaxy shreds. Our CNN gives a simple yet effective solution to clean galaxy catalogs from these contaminants. △ Less

Submitted 23 January, 2023; originally announced January 2023.

Comments: Accepted for publication in AJ

Study of the excess Fe XXV line emission in the central degrees of the Galactic centre using XMM-Newton data

Authors: K. Anastasopoulou, G. Ponti, M. C. Sormani, N. Locatelli, F. Haberl, M. R. Morris, E. M. Churazov, R. Sch�del, C. Maitra, S. Campana, E. M. Di Teodoro, C. Jin, I. Khabibullin, S. Mondal, M. Sasaki, Y. Zhang, X. Zheng

Abstract: The diffuse Fe XXV (6.7 keV) line emission observed in the Galactic ridge is widely accepted to be produced by a superposition of a large number of unresolved X-ray point sources. In the very central degrees of our Galaxy, however, the existence of an extremely hot ($\sim$7 keV) diffuse plasma is still under debate. In this work we measure the Fe XXV line emission using all available XMM-Newton ob… ▽ More The diffuse Fe XXV (6.7 keV) line emission observed in the Galactic ridge is widely accepted to be produced by a superposition of a large number of unresolved X-ray point sources. In the very central degrees of our Galaxy, however, the existence of an extremely hot ($\sim$7 keV) diffuse plasma is still under debate. In this work we measure the Fe XXV line emission using all available XMM-Newton observations of the Galactic centre (GC) and inner disc ($-10^{\circ}$$<\ell<10^{\circ}$, $-2^{\circ}<b<2^{\circ}$). We use recent stellar mass distribution models to estimate the amount of X-ray emission originating from unresolved point sources, and find that within a region of $\ell=\pm1^{\circ}$ and $b=\pm0.25^\circ$ the 6.7 keV emission is 1.3 to 1.5 times in excess of what is expected from unresolved point sources. The excess emission is enhanced towards regions where known supernova remnants are located, suggesting that at least a part of this emission is due to genuine diffuse very hot plasma. If the entire excess is due to very hot plasma, an energy injection rate of at least $\sim6\times10^{40}$ erg s$^{-1}$ is required, which cannot be provided by the measured supernova explosion rate or past Sgr A$^{*}$ activity alone. However, we find that almost the entire excess we observe can be explained by assuming GC stellar populations with iron abundances $\sim$1.9 times higher than those in the bar/bulge, a value that can be reproduced by fitting diffuse X-ray spectra from the corresponding regions. Even in this case, a leftover X-ray excess is concentrated within $\ell=\pm0.3^{\circ}$ and $b=\pm0.15^\circ$, corresponding to a thermal energy of $\sim2\times10^{52}$ erg, which can be reproduced by the estimated supernova explosion rate in the GC. Finally we discuss a possible connection to the observed GC Fermi-LAT excess. △ Less

Submitted 10 January, 2023; originally announced January 2023.

Comments: 21 pages, 11 figures, 7 tables, accepted for publication in A&A

Journal ref: A&A 671, A55 (2023)

WALLABY Pilot Survey: Public release of HI kinematic models for more than 100 galaxies from phase 1 of ASKAP pilot observations

Authors: N. Deg, K. Spekkens, T. Westmeier, T. N. Reynolds, P. Venkataraman, S. Goliath, A. X. Shen, R. Halloran, A. Bosma, B. Catinella, W. J. G. de Blok, H. D�nes, E. M. Di Teodoro, A. Elagali, B. -Q. For, C. Howlett, G. I. G. J�zsa, P. Kamphuis, D. Kleiner, B. Koribalski, K. Lee-Waddell, F. Lelli, X. Lin, C. Murugeshan, S. Oh , et al. (7 additional authors not shown)

Abstract: We present the Widefield ASKAP L-band Legacy All-sky Blind surveY (WALLABY) Pilot Phase I HI kinematic models. This first data release consists of HI observations of three fields in the direction of the Hydra and Norma clusters, and the NGC 4636 galaxy group. In this paper, we describe how we generate and publicly release flat-disk tilted-ring kinematic models for 109/592 unique HI detections in t… ▽ More We present the Widefield ASKAP L-band Legacy All-sky Blind surveY (WALLABY) Pilot Phase I HI kinematic models. This first data release consists of HI observations of three fields in the direction of the Hydra and Norma clusters, and the NGC 4636 galaxy group. In this paper, we describe how we generate and publicly release flat-disk tilted-ring kinematic models for 109/592 unique HI detections in these fields. The modelling method adopted here - which we call the WALLABY Kinematic Analysis Proto-Pipeline (WKAPP) and for which the corresponding scripts are also publicly available - consists of combining results from the homogeneous application of the FAT and 3DBAROLO algorithms to the subset of 209 detections with sufficient resolution and S/N in order to generate optimized model parameters and uncertainties. The 109 models presented here tend to be gas rich detections resolved by at least 3-4 synthesized beams across their major axes, but there is no obvious environmental bias in the modelling. The data release described here is the first step towards the derivation of similar products for thousands of spatially-resolved WALLABY detections via a dedicated kinematic pipeline. Such a large publicly available and homogeneously analyzed dataset will be a powerful legacy product that that will enable a wide range of scientific studies. △ Less

Submitted 14 November, 2022; originally announced November 2022.

Comments: Accepted to PASA

On the patchy appearance of the circum-Galactic medium and the influence of foreground absorption

Authors: G. Ponti, J. S. Sanders, N. Locatelli, X. Zheng, Y. Zhang, M. Yeung, M. Freyberg, K. Dennerl, J. Comparat, A. Merloni, E. Di Teodoro, M. Sasaki, T. H. Reiprich

Abstract: Recent studies have demonstrated that the emission from the circumgalactic medium of the Milky Way displays a relatively high degree of patchiness on angular scales of $\sim10^\circ$. Taking advantage of the Spectrum Roentgen Gamma eROSITA Final Equatorial Depth Survey (eFEDS), we aim to constrain any variation in the soft X-ray surface brightness on scales going from sub-degree to a hundred squar… ▽ More Recent studies have demonstrated that the emission from the circumgalactic medium of the Milky Way displays a relatively high degree of patchiness on angular scales of $\sim10^\circ$. Taking advantage of the Spectrum Roentgen Gamma eROSITA Final Equatorial Depth Survey (eFEDS), we aim to constrain any variation in the soft X-ray surface brightness on scales going from sub-degree to a hundred square degrees. We observe modulations of about $60$% on scales of several degrees and decreasing for higher energies. The observed patchiness is stable over a period of two years, therefore excluding that it is induced by Solar wind charge exchange. We also observe no correlation between such excess and the density of galaxies in the Local Universe, suggesting no strong contribution from the hot baryons in the filaments of the Cosmic web. Instead, the soft X-ray emission is anti-correlated with the column density of absorbing material. Indeed, we can reproduce the spectrum of the bright and dark regions by simply varying the column density of the matter absorbing the emission components located beyond the Local Hot Bubble, while no modulation of the intrinsic emission is required. At high Galactic latitudes, the eROSITA all sky map shows patchiness of the soft X-ray diffuse emission similar to the one observed in the eFEDS field, it is therefore likely that the same "absorption-modulation" is present over the entire sky. These results highlight the importance of an accurate treatment of the absorption effects, to determine the patchiness of the circumgalactic medium. △ Less

Submitted 3 November, 2022; originally announced November 2022.

Comments: Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 670, A99 (2023)

Dark matter halos and scaling relations of extremely massive spiral galaxies from extended HI rotation curves

Authors: Enrico M. Di Teodoro, Lorenzo Posti, S. Michael Fall, Patrick M. Ogle, Thomas Jarrett, Philip N. Appleton, Michelle E. Cluver, Martha P. Haynes, Ute Lisenfeld

Abstract: We present new and archival atomic hydrogen (\hi) observations of \galnum\ of the most massive spiral galaxies in the local Universe ($M_\star>10^{11} \, \mathrm{M}_\odot$). From 3D kinematic modeling of the datacubes, we derive extended \hi\ rotation curves, and from these, we estimate masses of the dark matter halos and specific angular momenta of the discs. We confirm that massive spiral galaxi… ▽ More We present new and archival atomic hydrogen (\hi) observations of \galnum\ of the most massive spiral galaxies in the local Universe ($M_\star>10^{11} \, \mathrm{M}_\odot$). From 3D kinematic modeling of the datacubes, we derive extended \hi\ rotation curves, and from these, we estimate masses of the dark matter halos and specific angular momenta of the discs. We confirm that massive spiral galaxies lie at the upper ends of the Tully-Fisher relation (mass vs velocity, $M \propto V^{4}$) and Fall relation (specific angular momentum vs mass, $j \propto M^{0.6}$), in both stellar and baryonic forms, with no significant deviations from single power laws. We study the connections between baryons and dark matter through the stellar (and baryon)-to-halo ratios of mass $f_\mathrm{M} \equiv M_\star/M_\mathrm{h}$ and specific angular momentum $f_\mathrm{j} \equiv j_\star/j_\mathrm{h}$ and $f_\mathrm{j,bar} \equiv j_\mathrm{bar}/j_\mathrm{h}$. Combining our sample with others from the literature for less massive disc-dominated galaxies, we find that $f_\mathrm{M}$ rises monotonically with $M_\star$ and $M_\mathrm{h}$ (instead of the inverted-U shaped $f_\mathrm{M}$ for spheroid-dominated galaxies), while $f_\mathrm{j}$ and $f_\mathrm{j,bar}$ are essentially constant near unity over four decades in mass. Our results indicate that disc galaxies constitute a self-similar population of objects closely linked to the self-similarity of their dark halos. This picture is reminiscent of early analytical models of galaxy formation wherein discs grow by relatively smooth and gradual inflow, isolated from disruptive events such as major mergers and strong AGN feedback, in contrast to the more chaotic growth of spheroids. △ Less

Submitted 21 November, 2022; v1 submitted 6 July, 2022; originally announced July 2022.

Comments: Accepted for publication on MNRAS

The impact of gas disc flaring on rotation curve decomposition and revisiting baryonic and dark-matter relations for nearby galaxies

Authors: Pavel E. Mancera Pi�a, Filippo Fraternali, Tom Oosterloo, Elizabeth A. K. Adams, Enrico di Teodoro, Cecilia Bacchini, Giuliano Iorio

Abstract: Gas discs of late-type galaxies are flared, with scale heights increasing with the distance from the galaxy centres and often reaching kpc scales. We study the effects of gas disc flaring on the recovered dark matter halo parameters from rotation curve decomposition. For this, we carefully select a sample of 32 dwarf and spiral galaxies with high-quality neutral gas, molecular gas, and stellar mas… ▽ More Gas discs of late-type galaxies are flared, with scale heights increasing with the distance from the galaxy centres and often reaching kpc scales. We study the effects of gas disc flaring on the recovered dark matter halo parameters from rotation curve decomposition. For this, we carefully select a sample of 32 dwarf and spiral galaxies with high-quality neutral gas, molecular gas, and stellar mass profiles, robust H\,{\sc i} rotation curves obtained via 3D kinematic modelling, and reliable bulge-disc decomposition. By assuming vertical hydrostatic equilibrium, we derive the scale heights of the atomic and molecular gas discs and fit dark matter haloes to the rotation curves self-consistently. We find that the effect of the gas flaring in the rotation curve decomposition can play an important role only for the smallest, gas-dominated dwarfs, while for most of the galaxies the effect is minor and can be ignored. We revisit the stellar- and baryon-to-halo mass relations ($M_\ast-M_{200}$ and $M_{\rm bar}-M_{200}$). Both relations increase smoothly up to $M_{200} \approx 10^{12}~\rm{ M_\odot}$, with galaxies at this end having high $M_\ast/M_{200}$ and $M_{\rm bar}/M_{200}$ ratios approaching the cosmological baryon fraction. At higher $M_{200}$ the relations show a larger scatter. Most haloes of our galaxy sample closely follow the concentration-mass ($c_{200}-M_{\rm 200}$) relation resulting from N-body cosmological simulations. Interestingly, the galaxies deviating above and below the relation have the highest and lowest stellar and baryon factions, respectively, which suggests that the departures from the $c_{200}-M_{\rm 200}$ law are regulated by adiabatic contraction and an increasing importance of feedback. △ Less

Submitted 25 May, 2022; originally announced May 2022.

Comments: accepted for publication in MNRAS. Supplementary material will be available from MNRAS and meantime at this link: https://www.dropbox.com/sh/768phag6ttgym0j/AABZwerBgFofkkZ5npHW17YBa

The PUMA project. III. Incidence and properties of ionised gas disks in ULIRGs, associated velocity dispersion and its dependence on starburstiness

Authors: M. Perna, S. Arribas, L. Colina, M. Pereira Santaella, I. Lamperti, E. Di Teodoro, H. �bler, L. Costantin, R. Maiolino, G. Cresci, E. Bellocchi, C. Catal�n-Torrecilla, S. Cazzoli, J. Piqueraz L�pez

Abstract: A classical scenario suggests that ULIRGs transform colliding spiral galaxies into a spheroid dominated early-type galaxy. Recent high-resolution simulations have instead shown that, under some circumstances, rotation disks can be preserved during the merging process or rapidly regrown after coalescence. Our goal is to analyze in detail the ionised gas kinematics in a sample of ULIRGs to infer the… ▽ More A classical scenario suggests that ULIRGs transform colliding spiral galaxies into a spheroid dominated early-type galaxy. Recent high-resolution simulations have instead shown that, under some circumstances, rotation disks can be preserved during the merging process or rapidly regrown after coalescence. Our goal is to analyze in detail the ionised gas kinematics in a sample of ULIRGs to infer the incidence of gas rotational dynamics in late-stage interacting galaxies and merger remnants. We analysed MUSE data of a sample of 20 nearby (z<0.165) ULIRGs, as part of the "Physics of ULIRGs with MUSE and ALMA" (PUMA) project. We found that 27% individual nuclei are associated with kpc-scale disk-like gas motions. The rest of the sample displays a plethora of gas kinematics, dominated by winds and merger-induced flows, which make the detection of rotation signatures difficult. On the other hand, the incidence of stellar disk-like motions is ~2 times larger than gaseous disks, as the former are probably less affected by winds and streams. The eight galaxies with a gaseous disk present relatively high intrinsic gas velocity dispersion (sigma = 30-85 km/s), rotationally-supported motions (with gas rotation velocity over velocity dispersion vrot/sigma > 1-8), and dynamical masses in the range (2-7)x1e10 Msun. By combining our results with those of local and high-z disk galaxies from the literature, we found a significant correlation between sigma and the offset from the main sequence (MS), after correcting for their evolutionary trends. Our results confirm the presence of kpc-scale rotating disks in interacting galaxies and merger remnants, with an incidence going from 27% (gas) to ~50% (stars). The ULIRGs gas velocity dispersion is up to a factor of ~4 higher than in local normal MS galaxies, similar to high-z starbursts as presented in the literature. △ Less

Submitted 4 February, 2022; originally announced February 2022.

Comments: Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 662, A94 (2022)

GASKAP-HI Pilot Survey Science I: ASKAP Zoom Observations of HI Emission in the Small Magellanic Cloud

Authors: N. M. Pingel, J. Dempsey, N. M. McClure-Griffiths, J. M. Dickey, K. E. Jameson, H. Arce, G. Anglada, J. Bland-Hawthorn, S. L. Breen, F. Buckland-Willis, S. E. Clark, J. R. Dawson, H. Dénes, E. M. Di Teodoro, B. -Q. For, Tyler J. Foster, J. F. Gómez, H. Imai, G. Joncas, C. -G. Kim, M. -Y. Lee, C. Lynn, D. Leahy, Y. K. Ma, A. Marchal , et al. (31 additional authors not shown)

Abstract: We present the most sensitive and detailed view of the neutral hydrogen (HI) emission associated with the Small Magellanic Cloud (SMC), through the combination of data from the Australian Square Kilometre Array Pathfinder (ASKAP) and Parkes (Murriyang), as part of the Galactic Australian Square Kilometre Array Pathfinder (GASKAP) pilot survey. These GASKAP-HI pilot observations, for the first time… ▽ More We present the most sensitive and detailed view of the neutral hydrogen (HI) emission associated with the Small Magellanic Cloud (SMC), through the combination of data from the Australian Square Kilometre Array Pathfinder (ASKAP) and Parkes (Murriyang), as part of the Galactic Australian Square Kilometre Array Pathfinder (GASKAP) pilot survey. These GASKAP-HI pilot observations, for the first time, reveal HI in the SMC on similar physical scales as other important tracers of the interstellar medium, such as molecular gas and dust. The resultant image cube possesses an rms noise level of 1.1 K (1.6 mJy/beam) per 0.98 km s$^{-1}$ spectral channel with an angular resolution of 30$''$ ($\sim$10 pc). We discuss the calibration scheme and the custom imaging pipeline that utilizes a joint deconvolution approach, efficiently distributed across a computing cluster, to accurately recover the emission extending across the entire $\sim$25 deg$^2$ field-of-view. We provide an overview of the data products and characterize several aspects including the noise properties as a function of angular resolution and the represented spatial scales by deriving the global transfer function over the full spectral range. A preliminary spatial power spectrum analysis on individual spectral channels reveals that the power-law nature of the density distribution extends down to scales of 10 pc. We highlight the scientific potential of these data by comparing the properties of an outflowing high velocity cloud with previous ASKAP+Parkes HI test observations. △ Less

Submitted 10 December, 2021; v1 submitted 9 November, 2021; originally announced November 2021.

Comments: Accepted for publication in PASA, 34 pages, 18 figures, 5 tables

Radial motions and radial gas flows in local spiral galaxies

Authors: Enrico M. Di Teodoro, Josh E. G. Peek

Abstract: We determine radial velocities and mass flow rates in a sample of 54 local spiral galaxies by modelling high-resolution and high-sensitivity data of the atomic hydrogen emission line. We found that, although radial inflow motions seem to be slightly preferred over outflow motions, their magnitude is generally small. Most galaxies show radial flows of only a few km/s throughout their HI disks, eith… ▽ More We determine radial velocities and mass flow rates in a sample of 54 local spiral galaxies by modelling high-resolution and high-sensitivity data of the atomic hydrogen emission line. We found that, although radial inflow motions seem to be slightly preferred over outflow motions, their magnitude is generally small. Most galaxies show radial flows of only a few km/s throughout their HI disks, either inwards or outwards, without any clear increase in magnitude in the outermost regions, as we would expect for continuous radial accretion. Gas mass flow rates for most galaxies are less than 1 M$_\odot$/yr. Over the entire sample, we estimated an average inflow rate of 0.3 M$_\odot$/yr outside the optical disk and of 0.1 M$_\odot$/yr in the outskirts of the HI disks. These inflow rates are about 5-10 times smaller than the average star formation rate of 1.4 M$_\odot$/yr. Our study suggests that there is no clear evidence for systematic radial accretion inflows that alone could feed and sustain the star formation process in the inner regions of local spiral galaxies at its current rate. △ Less

Submitted 4 October, 2021; originally announced October 2021.

Comments: Accepted for publications in ApJ. 12 pages, 6 figures + appendixes

MIGHTEE-HI: The baryonic Tully-Fisher relation over the last billion years

Authors: Anastasia A. Ponomareva, Wanga Mulaudzi, Natasha Maddox, Bradley S. Frank, Matt J. Jarvis, Enrico M. Di Teodoro, Marcin Glowacki, Renée C. Kraan-Korteweg, Tom A. Oosterloo, Elizabeth A. K. Adams, Hengxing Pan, Isabella Prandoni, Sambatriniaina H. A. Rajohnson, Francesco Sinigaglia, Nathan J. Adams, Ian Heywood, Rebecca A. A. Bowler, Peter W. Hatfield, Jordan D. Collier, Srikrishna Sekhar

Abstract: Using a sample of 67 galaxies from the MIGHTEE Survey Early Science data we study the HI-based baryonic Tully-Fisher relation (bTFr), covering a period of $\sim$one billion years ($0 \leq z \leq 0.081 $). We consider the bTFr based on two different rotational velocity measures: the width of the global HI profile and $\rm V_{out}$, measured as the outermost rotational velocity from the resolved HI… ▽ More Using a sample of 67 galaxies from the MIGHTEE Survey Early Science data we study the HI-based baryonic Tully-Fisher relation (bTFr), covering a period of $\sim$one billion years ($0 \leq z \leq 0.081 $). We consider the bTFr based on two different rotational velocity measures: the width of the global HI profile and $\rm V_{out}$, measured as the outermost rotational velocity from the resolved HI rotation curves. Both relations exhibit very low intrinsic scatter orthogonal to the best-fit relation ($σ_{\perp}=0.07\pm0.01$), comparable to the SPARC sample at $z \simeq 0$. The slopes of the relations are similar and consistent with the $ z \simeq 0$ studies ($3.66^{+0.35}_{-0.29}$ for $\rm W_{50}$ and $3.47^{+0.37}_{-0.30}$ for $\rm V_{out}$). We find no evidence that the bTFr has evolved over the last billion years, and all galaxies in our sample are consistent with the same relation independent of redshift and the rotational velocity measure. Our results set up a reference for all future studies of the HI-based bTFr as a function of redshift that will be conducted with the ongoing deep SKA pathfinders surveys. △ Less

Submitted 10 September, 2021; originally announced September 2021.

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

Rotation curves and scaling relations of extremely massive spiral galaxies

Authors: Enrico M. Di Teodoro, Lorenzo Posti, Patrick M. Ogle, S. Michael Fall, Thomas Jarrett

Abstract: We study the kinematics and scaling relations of a sample of 43 giant spiral galaxies that have stellar masses exceeding $10^{11}$ $M_\odot$ and optical discs up to 80 kpc in radius. We use a hybrid 3D-1D approach to fit 3D kinematic models to long-slit observations of the H$α$-[NII] emission lines and we obtain robust rotation curves of these massive systems. We find that all galaxies in our samp… ▽ More We study the kinematics and scaling relations of a sample of 43 giant spiral galaxies that have stellar masses exceeding $10^{11}$ $M_\odot$ and optical discs up to 80 kpc in radius. We use a hybrid 3D-1D approach to fit 3D kinematic models to long-slit observations of the H$α$-[NII] emission lines and we obtain robust rotation curves of these massive systems. We find that all galaxies in our sample seem to reach a flat part of the rotation curve within the outermost optical radius. We use the derived kinematics to study the high-mass end of the two most important scaling relations for spiral galaxies: the stellar/baryonic mass Tully-Fisher relation and the Fall (mass-angular momentum) relation. All galaxies in our sample, with the possible exception of the two fastest rotators, lie comfortably on both these scaling relations determined at lower masses, without any evident break or bend at the high-mass regime. When we combine our high-mass sample with lower-mass data from the Spitzer Photometry & Accurate Rotation Curves catalog, we find a slope of $α=4.25\pm0.19$ for the stellar Tully-Fisher relation and a slope of $γ=0.64\pm0.11$ for the Fall relation. Our results indicate that most, if not all, of these rare, giant spiral galaxies are scaled up versions of less massive discs and that spiral galaxies are a self-similar population of objects up to the very high-mass end. △ Less

Submitted 8 September, 2021; originally announced September 2021.

Comments: 12 pages, 5 figures + Appendix. Accepted for publication in MNRAS

A massive stellar bulge in a regularly rotating galaxy 1.2 billion years after the Big Bang

Authors: Federico Lelli, Enrico M. Di Teodoro, Filippo Fraternali, Allison W. S. Man, Zhi-Yu Zhang, Carlos De Breuck, Timothy A. Davis, Roberto Maiolino

Abstract: Cosmological models predict that galaxies forming in the early Universe experience a chaotic phase of gas accretion and star formation, followed by gas ejection due to feedback processes. Galaxy bulges may assemble later via mergers or internal evolution. Here we present submillimeter observations (with spatial resolution of 700 parsecs) of ALESS 073.1, a starburst galaxy at redshift z~5, when the… ▽ More Cosmological models predict that galaxies forming in the early Universe experience a chaotic phase of gas accretion and star formation, followed by gas ejection due to feedback processes. Galaxy bulges may assemble later via mergers or internal evolution. Here we present submillimeter observations (with spatial resolution of 700 parsecs) of ALESS 073.1, a starburst galaxy at redshift z~5, when the Universe was 1.2 billion years old. This galaxy's cold gas forms a regularly rotating disk with negligible noncircular motions. The galaxy rotation curve requires the presence of a central bulge in addition to a star-forming disk. We conclude that massive bulges and regularly rotating disks can form more rapidly in the early Universe than predicted by models of galaxy formation. △ Less

Submitted 11 February, 2021; originally announced February 2021.

Comments: Published in Science. This preprint corresponds to the accepted and language edited version of the manuscript. 36 pages, 9 figures, 2 tables

Cold gas in the Milky Way's nuclear wind

Authors: Enrico M. Di Teodoro, Naomi M. McClure-Griffiths, Felix J. Lockman, Lucia Armillotta

Abstract: The centre of the Milky Way is the site of several high-energy processes that have strongly impacted the inner regions of our Galaxy. Activity from the super-massive black hole, Sgr A*, and/or stellar feedback from the inner molecular ring expel matter and energy from the disc in the form of a galactic wind. Multiphase gas has been observed within this outflow, from hot highly-ionized, to warm ion… ▽ More The centre of the Milky Way is the site of several high-energy processes that have strongly impacted the inner regions of our Galaxy. Activity from the super-massive black hole, Sgr A*, and/or stellar feedback from the inner molecular ring expel matter and energy from the disc in the form of a galactic wind. Multiphase gas has been observed within this outflow, from hot highly-ionized, to warm ionized and cool atomic gas. To date, however, there has been no evidence of the cold and dense molecular phase. Here we report the first detection of molecular gas outflowing from the centre of our Galaxy. This cold material is associated with atomic hydrogen clouds travelling in the nuclear wind. The morphology and the kinematics of the molecular gas, resolved on ~1 pc scale, indicate that these clouds are mixing with the warmer medium and are possibly being disrupted. The data also suggest that the mass of molecular gas driven out is not negligible and could impact the rate of star formation in the central regions. The presence of this cold, dense, high-velocity gas is puzzling, as neither Sgr A* at its current level of activity, nor star formation in the inner Galaxy seem viable sources for this material. △ Less

Submitted 20 August, 2020; originally announced August 2020.

Comments: Published in the August 19 issue of Nature. This is the authors' version before final edits. Published version is available at http://www.nature.com/articles/s41586-020-2595-z

Journal ref: Nature, 2020, vol. 584, pp. 364-367

Robust HI kinematics of gas-rich ultra-diffuse galaxies: hints of a weak-feedback formation scenario

Authors: Pavel E. Mancera Piña, Filippo Fraternali, Kyle A. Oman, Elizabeth A. K. Adams, Cecilia Bacchini, Antonino Marasco, Tom Oosterloo, Gabriele Pezzulli, Lorenzo Posti, Lukas Leisman, John M. Cannon, Enrico M. di Teodoro, Lexi Gault, Martha P. Haynes, Kameron Reiter, Katherine L. Rhode, John J. Salzer, Nicholas J. Smith

Abstract: We study the gas kinematics of a sample of six isolated gas-rich low surface brightness galaxies, of the class called ultra-diffuse galaxies (UDGs). These galaxies have recently been shown to be outliers from the baryonic Tully-Fisher relation (BTFR), as they rotate much slower than expected given their baryonic mass, and to have baryon fractions similar to the cosmological mean. By means of a 3D… ▽ More We study the gas kinematics of a sample of six isolated gas-rich low surface brightness galaxies, of the class called ultra-diffuse galaxies (UDGs). These galaxies have recently been shown to be outliers from the baryonic Tully-Fisher relation (BTFR), as they rotate much slower than expected given their baryonic mass, and to have baryon fractions similar to the cosmological mean. By means of a 3D kinematic modelling fitting technique, we show that the HI in our UDGs is distributed in "thin" regularly rotating discs and we determine their rotation velocity and gas velocity dispersion. We revisit the BTFR adding galaxies from other studies. We find a previously unknown trend between the deviation from the BTFR and the disc scale length valid for dwarf galaxies with circular speeds < 45 km/s, with our UDGs being at the extreme end. Based on our findings, we suggest that the high baryon fractions of our UDGs may originate due to the fact that they have experienced weak stellar feedback, likely due to their low star formation rate surface densities, and as a result they did not eject significant amounts of gas out of their discs. At the same time, we find indications that our UDGs may have higher-than-average stellar specific angular momentum, which can explain their large optical scale lengths. △ Less

Submitted 3 May, 2020; v1 submitted 29 April, 2020; originally announced April 2020.

Comments: Accepted for publication in MNRAS. v2: a few typos have been corrected and a couple of references added

The Life Cycle of the Central Molecular Zone. II: Distribution of atomic and molecular gas tracers

Authors: Lucia Armillotta, Mark R. Krumholz, Enrico M. Di Teodoro

Abstract: We use the hydrodynamical simulation of our inner Galaxy presented in Armillotta et al. (2019) to study the gas distribution and kinematics within the CMZ. We use a resolution high enough to capture the gas emitting in dense molecular tracers such as NH3 and HCN, and simulate a time window of 50 Myr, long enough to capture phases during which the CMZ experiences both quiescent and intense star for… ▽ More We use the hydrodynamical simulation of our inner Galaxy presented in Armillotta et al. (2019) to study the gas distribution and kinematics within the CMZ. We use a resolution high enough to capture the gas emitting in dense molecular tracers such as NH3 and HCN, and simulate a time window of 50 Myr, long enough to capture phases during which the CMZ experiences both quiescent and intense star formation. We then post-process the simulated CMZ to calculate its spatially-dependent chemical and thermal state, producing synthetic emission data cubes and maps of both HI and the molecular gas tracers CO, NH3 and HCN. We show that, as viewed from Earth, gas in the CMZ is distributed mainly in two parallel and elongated features extending from positive longitudes and velocities to negative longitudes and velocities. The molecular gas emission within these two streams is not uniform, and it is mostly associated to the region where gas flowing towards the Galactic Center through the dust lanes collides with gas orbiting within the ring. Our simulated data cubes reproduce a number of features found in the observed CMZ. However, some discrepancies emerge when we use our results to interpret the position of individual molecular clouds. Finally, we show that, when the CMZ is near a period of intense star formation, the ring is mostly fragmented as a consequence of supernova feedback, and the bulk of the emission comes from star-forming molecular clouds. This correlation between morphology and star formation rate should be detectable in observations of extragalactic CMZs. △ Less

Submitted 14 February, 2020; originally announced February 2020.

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

Observation of Acceleration of HI Clouds Within the Fermi Bubbles

Authors: Felix J. Lockman, Enrico M. Di Teodoro, N. M. McClure-Griffiths

Abstract: The ~200 HI clouds observed to be entrained in the Fermi Bubble wind show a trend of increasing maximum |VLSR| with Galactic latitude. We analyze previous observations and present new data from the Green Bank Telescope that rule out systematic effects as the source of this phenomenon. Instead, it is likely evidence for acceleration of the clouds. The data suggest that clouds in the lower 2 kpc of… ▽ More The ~200 HI clouds observed to be entrained in the Fermi Bubble wind show a trend of increasing maximum |VLSR| with Galactic latitude. We analyze previous observations and present new data from the Green Bank Telescope that rule out systematic effects as the source of this phenomenon. Instead, it is likely evidence for acceleration of the clouds. The data suggest that clouds in the lower 2 kpc of the Fermi Bubbles, within the Bubble boundaries established from X-ray studies, have an outflow velocity that rises from ~150 - 200 km/s close to the Galactic Center and reaches ~330 km/s at a distance of 2.5 - 3.5 kpc. These parameters are also consistent with the kinematics of UV absorption lines from highly ionized species observed against two targets behind the Fermi Bubbles at $b = -6.6^{\circ}$, and $b = +11.2^{\circ}$. The implied neutral cloud lifetime is 4 - 10 Myr. △ Less

Submitted 15 November, 2019; originally announced November 2019.

Comments: Accepted for publication in the Astrophysical Journal

The 3D Kinematics of Gas in the Small Magellanic Cloud

Authors: Claire E. Murray, J. E. G. Peek, Enrico M. Di Teodoro, N. M. McClure-Griffiths, John M. Dickey, Helga Denes

Abstract: We investigate the kinematics of neutral gas in the Small Magellanic Cloud (SMC) and test the hypothesis that it is rotating in a disk. To trace the 3D motions of the neutral gas distribution, we identify a sample of young, massive stars embedded within it. These are stars with radial velocity measurements from spectroscopic surveys and proper motion measurements from Gaia, whose radial velocities… ▽ More We investigate the kinematics of neutral gas in the Small Magellanic Cloud (SMC) and test the hypothesis that it is rotating in a disk. To trace the 3D motions of the neutral gas distribution, we identify a sample of young, massive stars embedded within it. These are stars with radial velocity measurements from spectroscopic surveys and proper motion measurements from Gaia, whose radial velocities match with dominant HI components. We compare the observed radial and tangential velocities of these stars with predictions from the state-of-the-art rotating disk model based on high-resolution 21 cm observations of the SMC from the Australian Square Kilometer Array Pathfinder telescope. We find that the observed kinematics of gas-tracing stars are inconsistent with disk rotation. We conclude that the kinematics of gas in the SMC are more complex than can be inferred from the integrated radial velocity field. As a result of violent tidal interactions with the LMC, non-rotational motions are prevalent throughout the SMC, and it is likely composed of distinct sub-structures overlapping along the line of sight. △ Less

Submitted 24 October, 2019; originally announced October 2019.

Comments: 9 pages, 5 figures, 1 Appendix; ApJ accepted

Molecular gas in the outflow of the Small Magellanic Cloud

Authors: Enrico M. Di Teodoro, N. M. McClure-Griffiths, C. De Breuck, L. Armillotta, N. M. Pingel, K. E. Jameson, J. M. Dickey, M. Rubio, S. Stanimirovic, L. Staveley-Smith

Abstract: We report the first evidence of molecular gas in two atomic hydrogen (HI) clouds associated with gas outflowing from the Small Magellanic Cloud (SMC). We used the Atacama Pathfinder Experiment (APEX) to detect and spatially resolve individual clumps of CO(2-1) emission in both clouds. CO clumps are compact (~ 10 pc) and dynamically cold (linewidths < 1 km/s). Most CO emission appears to be offset… ▽ More We report the first evidence of molecular gas in two atomic hydrogen (HI) clouds associated with gas outflowing from the Small Magellanic Cloud (SMC). We used the Atacama Pathfinder Experiment (APEX) to detect and spatially resolve individual clumps of CO(2-1) emission in both clouds. CO clumps are compact (~ 10 pc) and dynamically cold (linewidths < 1 km/s). Most CO emission appears to be offset from the peaks of the HI emission, some molecular gas lies in regions without a clear HI counterpart. We estimate a total molecular gas mass of 10^3-10^4 Msun in each cloud and molecular gas fractions up to 30% of the total cold gas mass (molecular + neutral). Under the assumption that this gas is escaping the galaxy, we calculated a cold gas outflow rate of 0.3-1.8 Msun/yr and mass loading factors of 3 -12 at a distance larger than 1 kpc. These results show that relatively weak star-formation-driven winds in dwarf galaxies like the SMC are able to accelerate significant amounts of cold and dense matter and inject it into the surrounding environment. △ Less

Submitted 23 October, 2019; originally announced October 2019.

Comments: Accepted for publication in ApJ Letters

Off the baryonic Tully-Fisher relation: a population of baryon-dominated ultra-diffuse galaxies

Authors: Pavel E. Mancera Pi�a, Filippo Fraternali, Elizabeth A. K. Adams, Antonino Marasco, Tom Oosterloo, Kyle A. Oman, Lukas Leisman, Enrico M. di Teodoro, Lorenzo Posti, Michael Battipaglia, John M. Cannon, Lexi Gault, Martha P. Haynes, Steven Janowiecki, Elizabeth McAllan, Hannah J. Pagel, Kameron Reiter, Katherine L. Rhode, John J. Salzer, Nicholas J. Smith

Abstract: We study the gas kinematics traced by the 21-cm emission of a sample of six HI$-$rich low surface brightness galaxies classified as ultra-diffuse galaxies (UDGs). Using the 3D kinematic modelling code $\mathrm{^{3D}}$Barolo we derive robust circular velocities, revealing a startling feature: HI$-$rich UDGs are clear outliers from the baryonic Tully-Fisher relation, with circular velocities much lo… ▽ More We study the gas kinematics traced by the 21-cm emission of a sample of six HI$-$rich low surface brightness galaxies classified as ultra-diffuse galaxies (UDGs). Using the 3D kinematic modelling code $\mathrm{^{3D}}$Barolo we derive robust circular velocities, revealing a startling feature: HI$-$rich UDGs are clear outliers from the baryonic Tully-Fisher relation, with circular velocities much lower than galaxies with similar baryonic mass. Notably, the baryon fraction of our UDG sample is consistent with the cosmological value: these UDGs are compatible with having no "missing baryons" within their virial radii. Moreover, the gravitational potential provided by the baryons is sufficient to account for the amplitude of the rotation curve out to the outermost measured point, contrary to other galaxies with similar circular velocities. We speculate that any formation scenario for these objects will require very inefficient feedback and a broad diversity in their inner dark matter content. △ Less

Submitted 7 September, 2019; v1 submitted 3 September, 2019; originally announced September 2019.

Comments: Accepted for publication by The Astrophysical Journal Letters (ApJL). V2: Acknowledgments have been updated, and a typo has been corrected

A multi-wavelength study of a massive, active galaxy at $z\sim 2$: coupling the kinematics of the ionized and molecular gas

Authors: Federica Loiacono, Margherita Talia, Filippo Fraternali, Andrea Cimatti, Enrico Maria Di Teodoro, Gabriel Bartosch Caminha

Abstract: We report a multi-wavelength study of the massive ($M_{\star} \gtrsim 10^{11} \rm{M}_{\odot}$), $z\sim 2$ star-forming galaxy GMASS 0953, which hosts an obscured AGN. We combined near-infrared observations of the GNIRS, SINFONI and KMOS spectrographs to study the kinematics of the [O III]$λ5007$ and H$α$ emission lines. Our analysis shows that GMASS 0953 may host an ionized disc extending up to 13… ▽ More We report a multi-wavelength study of the massive ($M_{\star} \gtrsim 10^{11} \rm{M}_{\odot}$), $z\sim 2$ star-forming galaxy GMASS 0953, which hosts an obscured AGN. We combined near-infrared observations of the GNIRS, SINFONI and KMOS spectrographs to study the kinematics of the [O III]$λ5007$ and H$α$ emission lines. Our analysis shows that GMASS 0953 may host an ionized disc extending up to 13 kpc, which rotates at a velocity of $V_{\rm{ion}} = 203^{+17}_{-20} {\ \rm kms}^{-1}$ at the outermost radius. Evidence of rotation on a smaller scale ($R \sim 1$ kpc) arises from the CO(J=6-5) line. The central velocity $V_{\rm{CO}} = 320^{+ 92}_{-53} {\ \rm kms}^{-1}$ traced by the molecular gas is higher than $V_{\rm{ion}}$, suggesting that the galaxy harbors a multi-phase disc with a rotation curve that peaks in the very central regions. The galaxy appears well located on the $z = 0$ baryonic Tully-Fisher relation. We also discuss the possibility that the [O III]$λ5007$ and H$α$ velocity gradients are due to a galactic-scale wind. Besides, we found evidence of an AGN-driven outflow traced by a broad blueshifted wing affecting the [O III]$λ5007$ line, which presents a velocity offset $Δv = -535 \pm 152 {\ \rm kms}^{-1}$ from the systemic velocity. Because of the short depletion timescale ($τ_{\rm{dep}}\sim 10^8$ yr) due to gas ejection and gas consumption by star formation activity, GMASS 0953 may likely evolve into a passive galaxy. However, the role of the AGN in depleting the gas reservoir of the galaxy is quite unclear because of the uncertainties affecting the outflow rate. △ Less

Submitted 4 August, 2019; originally announced August 2019.

Comments: 18 pages, 12 figures

Separating Line Emission from Star Formation, Shocks, and AGN Ionisation in NGC 1068

Authors: Joshua J. D'Agostino, Lisa J. Kewley, Brent A. Groves, Anne M. Medling, Enrico Di Teodoro, Michael A. Dopita, Adam D. Thomas, Ralph S. Sutherland, Santiago Garcia-Burillo

Abstract: In the optical spectra of galaxies, the separation of line emission from gas ionised by star formation and an AGN, or by star formation and shocks, are very well-understood problems. However, separating line emission between AGN and shocks has proven difficult. With the aid of a new three-dimensional diagnostic diagram, we show the simultaneous separation of line emission from star formation, shoc… ▽ More In the optical spectra of galaxies, the separation of line emission from gas ionised by star formation and an AGN, or by star formation and shocks, are very well-understood problems. However, separating line emission between AGN and shocks has proven difficult. With the aid of a new three-dimensional diagnostic diagram, we show the simultaneous separation of line emission from star formation, shocks, and AGN in NGC 1068, and quantify the ratio of star formation, shocks, and AGN in each spaxel. The AGN, shock, and star formation luminosity distributions across the galaxy accurately align with X-ray, radio, and CO(3-2) observations, respectively. Comparisons with previous separation methods show that the shocked emission heavily mixes with the AGN emission. We also show that if the H$α$ flux is to be used as a star formation rate indicator, separating line emission from as many sources as possible should be attempted to ensure accurate results. △ Less

Submitted 19 June, 2019; originally announced June 2019.

Comments: 17 pages, 16 figures, accepted for publication in MNRAS

Testing Star Formation Laws on Spatially Resolved Regions in a $z \approx 4.3$ Starburst Galaxy

Authors: Piyush Sharda, Elisabete da Cunha, Christoph Federrath, Emily Wisnioski, Enrico di Teodoro, Ken-ichi Tadaki, Min Yun, Itziar Aretxaga, Ryohei Kawabe

Abstract: We probe the star formation properties of the gas in AzTEC-1 in the COSMOS field, one of the best resolved and brightest starburst galaxies at $z \approx 4.3$, forming stars at a rate > 1000 $\mathrm{M_{\odot}}\,\mathrm{yr^{-1}}$. Using recent ALMA observations, we study star formation in the galaxy nucleus and an off-center star-forming clump and measure a median star formation rate (SFR) surface… ▽ More We probe the star formation properties of the gas in AzTEC-1 in the COSMOS field, one of the best resolved and brightest starburst galaxies at $z \approx 4.3$, forming stars at a rate > 1000 $\mathrm{M_{\odot}}\,\mathrm{yr^{-1}}$. Using recent ALMA observations, we study star formation in the galaxy nucleus and an off-center star-forming clump and measure a median star formation rate (SFR) surface density of $Σ^{\mathrm{nucleus}}_{\mathrm{SFR}} = 270\pm54$ and $Σ^{\mathrm{sfclump}}_{\mathrm{SFR}} = 170\pm38\,\mathrm{M_{\odot}}\,\mathrm{yr}^{-1}\,\mathrm{kpc}^{-2}$, respectively. Following the analysis by Sharda et al. (2018), we estimate the molecular gas mass, freefall time and turbulent Mach number in these regions to predict $Σ_{\mathrm{SFR}}$ from three star formation relations in the literature. The Kennicutt-Schmidt (Kennicutt 1998, KS) relation, which is based on the gas surface density, underestimates the $Σ_{\mathrm{SFR}}$ in these regions by a factor 2-3. The $Σ_{\mathrm{SFR}}$ we calculate from the single-freefall model of Krumholz et al. 2012 (KDM) is consistent with the measured $Σ_{\mathrm{SFR}}$ in the nucleus and the star-forming clump within the uncertainties. The turbulence-regulated star formation relation by Salim et al. 2015 (SFK) agrees slightly better with the observations than the KDM relation. Our analysis reveals that an interplay between turbulence and gravity can help sustain high SFRs in high-redshift starbursts. It can also be extended to other high- and low-redshift galaxies thanks to the high angular resolution and sensitivity of ALMA observations. △ Less

Submitted 3 June, 2019; originally announced June 2019.

Comments: 9 pages, 5 figures, accepted for publication in MNRAS

The Life Cycle of the Central Molecular Zone. I: Inflow, Star Formation, and Winds

Authors: L. Armillotta, M. R. Krumholz, E. M. Di Teodoro, N. M. McClure-Griffiths

Abstract: We present a study of the gas cycle and star formation history in the central 500 pc of the Milky Way, known as Central Molecular Zone (CMZ). Through hydrodynamical simulations of the inner 4.5 kpc of our Galaxy, we follow the gas cycle in a completely self-consistent way, starting from gas radial inflow due to the Galactic bar, the channelling of this gas into a dense, star-forming ring/stream at… ▽ More We present a study of the gas cycle and star formation history in the central 500 pc of the Milky Way, known as Central Molecular Zone (CMZ). Through hydrodynamical simulations of the inner 4.5 kpc of our Galaxy, we follow the gas cycle in a completely self-consistent way, starting from gas radial inflow due to the Galactic bar, the channelling of this gas into a dense, star-forming ring/stream at ~ 200 - 300 pc from the Galactic centre, and the launching of galactic outflows powered by stellar feedback. We find that star formation activity in the CMZ goes through oscillatory burst/quench cycles, with a period of tens to hundreds of Myr, characterised by roughly constant gas mass but order-of-magnitude level variations in the star formation rate. Comparison with the observed present-day star formation rate of the CMZ suggests that we are currently near a minimum of this cycle. Stellar feedback drives a mainly two-phase wind off the Galactic disc. The warm phase dominates the mass flux, and carries 100 - 200 % of the gas mass converted into stars. However, most of this gas goes into a fountain and falls back onto the disc rather than escaping the Galaxy. The hot phase carries most of the energy, with a time-averaged energy outflow rate of 10 - 20 % of the supernova energy budget. △ Less

Submitted 10 October, 2019; v1 submitted 3 May, 2019; originally announced May 2019.

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

The SAMI Galaxy Survey: Bayesian Inference for Gas Disk Kinematics using a Hierarchical Gaussian Mixture Model

Authors: Mathew R. Varidel, Scott M. Croom, Geraint F. Lewis, Brendon J. Brewer, Enrico M. Di Teodoro, Joss Bland-Hawthorn, Julia J. Bryant, Christoph Federrath, Caroline Foster, Karl Glazebrook, Michael Goodwin, Brent Groves, Andrew M. Hopkins, Jon S. Lawrence, Ángel R. López-Sánchez, Anne M. Medling, Matt S. Owers, Samuel N. Richards, Richard Scalzo, Nicholas Scott, Sarah M. Sweet, Dan S. Taranu, Jesse van de Sande

Abstract: We present a novel Bayesian method, referred to as Blobby3D, to infer gas kinematics that mitigates the effects of beam smearing for observations using Integral Field Spectroscopy (IFS). The method is robust for regularly rotating galaxies despite substructure in the gas distribution. Modelling the gas substructure within the disk is achieved by using a hierarchical Gaussian mixture model. To acco… ▽ More We present a novel Bayesian method, referred to as Blobby3D, to infer gas kinematics that mitigates the effects of beam smearing for observations using Integral Field Spectroscopy (IFS). The method is robust for regularly rotating galaxies despite substructure in the gas distribution. Modelling the gas substructure within the disk is achieved by using a hierarchical Gaussian mixture model. To account for beam smearing effects, we construct a modelled cube that is then convolved per wavelength slice by the seeing, before calculating the likelihood function. We show that our method can model complex gas substructure including clumps and spiral arms. We also show that kinematic asymmetries can be observed after beam smearing for regularly rotating galaxies with asymmetries only introduced in the spatial distribution of the gas. We present findings for our method applied to a sample of 20 star-forming galaxies from the SAMI Galaxy Survey. We estimate the global H$α$ gas velocity dispersion for our sample to be in the range $\barσ_v \sim $[7, 30] km s$^{-1}$. The relative difference between our approach and estimates using the single Gaussian component fits per spaxel is $Δ\barσ_v / \barσ_v = - 0.29 \pm 0.18$ for the H$α$ flux-weighted mean velocity dispersion. △ Less

Submitted 10 March, 2019; v1 submitted 7 March, 2019; originally announced March 2019.

Comments: 23 pages, 12 figures, accepted for MNRAS

The angular momentum of disc galaxies at z=1

Authors: A. Marasco, F. Fraternali, L. Posti, M. Ijtsma, E. M. Di Teodoro, T. Oosterloo

Abstract: We investigate the relation between stellar mass and specific stellar angular momentum, or `Fall relation', for a sample of 17 isolated, regularly rotating disc galaxies at z=1. All galaxies have a) rotation curves determined from Halpha emission-line data; b) HST imaging in optical and infrared filters; c) robust determinations of their stellar masses. We use HST images in f814w and f160w filters… ▽ More We investigate the relation between stellar mass and specific stellar angular momentum, or `Fall relation', for a sample of 17 isolated, regularly rotating disc galaxies at z=1. All galaxies have a) rotation curves determined from Halpha emission-line data; b) HST imaging in optical and infrared filters; c) robust determinations of their stellar masses. We use HST images in f814w and f160w filters, roughly corresponding to rest-frames B and I bands, to extract surface brightness profiles for our systems. We robustly bracket the specific angular momentum by assuming that rotation curves beyond the outermost Halpha rotation point stay either flat or follow a Keplerian fall-off. By comparing our measurements with those determined for disc galaxies in the local Universe, we find no evolution in the Fall relation in the redshift range 0<z<1, regardless of the band used and despite the uncertainties in the stellar rotation curves at large radii. This result holds unless stellar masses at z=1 are systematically underestimated by more than 50%. Our findings are compatible with expectations based on a LCDM cosmological framework and support a scenario where both the stellar Tully-Fisher and mass-size relations for spirals do not evolve significantly in this redshift range. △ Less

Submitted 20 December, 2018; originally announced December 2018.

Comments: 11 pages, 4 figures, 1 table. Accepted by A&A

Journal ref: A&A 621, L6 (2019)

On the dynamics of the Small Magellanic Cloud through high-resolution ASKAP HI observations

Authors: E. M. Di Teodoro, N. M. McClure-Griffiths, K. E. Jameson, H. Denes, John M. Dickey, S. Stanimirovic, L. Staveley-Smith, C. Anderson, J. D. Bunton, A. Chippendale, K. Lee-Waddell, A. MacLeod, M. A Voronkov

Abstract: We use new high-resolution HI data from the Australian Square Kilometre Array Pathfinder (ASKAP) to investigate the dynamics of the Small Magellanic Cloud (SMC). We model the HI gas component as a rotating disc of non-negligible angular size, moving into the plane of the sky and undergoing nutation/precession motions. We derive a high-resolution (~ 10 pc) rotation curve of the SMC out to R ~ 4 kpc… ▽ More We use new high-resolution HI data from the Australian Square Kilometre Array Pathfinder (ASKAP) to investigate the dynamics of the Small Magellanic Cloud (SMC). We model the HI gas component as a rotating disc of non-negligible angular size, moving into the plane of the sky and undergoing nutation/precession motions. We derive a high-resolution (~ 10 pc) rotation curve of the SMC out to R ~ 4 kpc. After correcting for asymmetric drift, the circular velocity slowly rises to a maximum value of Vc ~ 55 km/s at R ~ 2.8 kpc and possibly flattens outwards. In spite of the SMC undergoing strong gravitational interactions with its neighbours, its HI rotation curve is akin to that of many isolated gas-rich dwarf galaxies. We decompose the rotation curve and explore different dynamical models to deal with the unknown three-dimensional shape of the mass components (gas, stars and dark matter). We find that, for reasonable mass-to-light ratios, a dominant dark matter halo with mass M(R<4 kpc) = 1-1.5 x 10^9 solar masses is always required to successfully reproduce the observed rotation curve, implying a large baryon fraction of 30%-40%. We discuss the impact of our assumptions and the limitations of deriving the SMC kinematics and dynamics from HI observations. △ Less

Submitted 23 November, 2018; originally announced November 2018.

Comments: Accepted for publication in MNRAS

Cold gas outflows from the Small Magellanic Cloud traced with ASKAP

Authors: N. M. McClure-Griffiths, H. Dénes, J. M. Dickey, S. Stanimirović, L. Staveley-Smith, Katherine Jameson, Enrico Di Teodoro, James R. Allison, J. D. Collier, A. P. Chippendale, T. Franzen, Gülay Gürkan, G. Heald, A. Hotan, D. Kleiner, K. Lee-Waddell, D. McConnell, A. Popping, Jonghwan Rhee, C. J. Riseley, M. A. Voronkov, M. Whiting

Abstract: Feedback from massive stars plays a critical role in the evolution of the Universe by driving powerful outflows from galaxies that enrich the intergalactic medium and regulate star formation. An important source of outflows may be the most numerous galaxies in the Universe: dwarf galaxies. With small gravitational potential wells, these galaxies easily lose their star-forming material in the prese… ▽ More Feedback from massive stars plays a critical role in the evolution of the Universe by driving powerful outflows from galaxies that enrich the intergalactic medium and regulate star formation. An important source of outflows may be the most numerous galaxies in the Universe: dwarf galaxies. With small gravitational potential wells, these galaxies easily lose their star-forming material in the presence of intense stellar feedback. Here, we show that the nearby dwarf galaxy, the Small Magellanic Cloud (SMC), has atomic hydrogen outflows extending at least 2 kiloparsecs (kpc) from the star-forming bar of the galaxy. The outflows are cold, $T<400~{\rm K}$, and may have formed during a period of active star formation $25 - 60$ million years (Myr) ago. The total mass of atomic gas in the outflow is $\sim 10^7$ solar masses, ${\rm M_{\odot}}$, or $\sim 3$% of the total atomic gas of the galaxy. The inferred mass flux in atomic gas alone, $\dot{M}_{HI}\sim 0.2 - 1.0~{\rm M_{\odot}~yr^{-1}}$, is up to an order of magnitude greater than the star formation rate. We suggest that most of the observed outflow will be stripped from the SMC through its interaction with its companion, the Large Magellanic Cloud (LMC), and the Milky Way, feeding the Magellanic Stream of hydrogen encircling the Milky Way. △ Less

Submitted 5 November, 2018; originally announced November 2018.

Comments: Published in Nature Astronomy, 29 October 2018, http://dx.doi.org/10.1038/s41550-018-0608-8

Journal ref: Nature Astronomy 2, pg 901 - 906 (2018)

A novel 3D technique to study the kinematics of lensed galaxies

Authors: Francesca Rizzo, Simona Vegetti, Filippo Fraternali, Enrico Di Teodoro

Abstract: We present a 3D Bayesian method to model the kinematics of strongly lensed galaxies from spatially-resolved emission-line observations. This technique enables us to simultaneously recover the lens-mass distribution and the source kinematics directly from the 3D data cube. We have tested this new method with simulated OSIRIS observations for nine star-forming lensed galaxies with different kinemati… ▽ More We present a 3D Bayesian method to model the kinematics of strongly lensed galaxies from spatially-resolved emission-line observations. This technique enables us to simultaneously recover the lens-mass distribution and the source kinematics directly from the 3D data cube. We have tested this new method with simulated OSIRIS observations for nine star-forming lensed galaxies with different kinematic properties. The simulated rotation curves span a range of shapes which are prototypes of different morphological galaxy types, from dwarf to massive spiral galaxies. We have found that the median relative accuracy on the inferred lens and kinematic parameters are at the level of 1 and 2 per cent, respectively. We have also tested the robustness of the technique against different inclination angles, signal-to-noise ratios, the presence of warps or non-circular motions and we have found that the accuracy stays within a few per cent in most cases. This technique represents a significant step forward with respect to the methods used until now, as the lens parameters and the kinematics of the source are derived from the same 3D data. This enables us to study the possible degeneracies between the two and estimate the uncertainties on all model parameters consistently. △ Less

Submitted 19 September, 2018; originally announced September 2018.

Comments: Accepted for publication in MNRAS

The angular momentum-mass relation: a fundamental law from dwarf irregulars to massive spirals

Authors: Lorenzo Posti, Filippo Fraternali, Enrico di Teodoro, Gabriele Pezzulli

Abstract: In a $Λ$CDM Universe, the specific stellar angular momentum ($j_\ast$) and stellar mass ($M_\ast$) of a galaxy are correlated as a consequence of the scaling existing for dark matter haloes ($j_{\rm h}\propto M_{\rm h}^{2/3}$). The shape of this law is crucial to test galaxy formation models, which are currently discrepant especially at the lowest masses, allowing to constrain fundamental paramete… ▽ More In a $Λ$CDM Universe, the specific stellar angular momentum ($j_\ast$) and stellar mass ($M_\ast$) of a galaxy are correlated as a consequence of the scaling existing for dark matter haloes ($j_{\rm h}\propto M_{\rm h}^{2/3}$). The shape of this law is crucial to test galaxy formation models, which are currently discrepant especially at the lowest masses, allowing to constrain fundamental parameters, e.g. the retained fraction of angular momentum. In this study, we accurately determine the empirical $j_\ast-M_\ast$ relation (Fall relation) for 92 nearby spiral galaxies (from S0 to Irr) selected from the Spitzer Photometry and Accurate Rotation Curves (SPARC) sample in the unprecedented mass range $7 \lesssim \log M_\ast/M_\odot \lesssim 11.5$. We significantly improve all previous estimates of the Fall relation by determining $j_\ast$ profiles homogeneously for all galaxies, using extended HI rotation curves, and selecting only galaxies for which a robust $j_\ast$ could be measured (converged $j_\ast(<R)$ radial profile). We find the relation to be well described by a single, unbroken power-law $j_\ast\propto M_\ast^α$ over the entire mass range, with $α=0.55\pm 0.02$ and orthogonal intrinsic scatter of $0.17\pm 0.01$ dex. We finally discuss some implications for galaxy formation models of this fundamental scaling law and, in particular, the fact that it excludes models in which discs of all masses retain the same fraction of the halo angular momentum. △ Less

Submitted 12 April, 2018; originally announced April 2018.

Comments: A&A Letters, accepted

Journal ref: A&A 612, L6 (2018)

Blowing in the Milky Way wind: neutral hydrogen clouds tracing the Galactic nuclear outflow

Authors: Enrico M. Di Teodoro, Naomi M. McClure-Griffiths, Felix J. Lockman, Sara R. Denbo, Ryan Endsley, H. Alyson Ford, Kevin Harrington

Abstract: We present the results of a new sensitive survey of neutral hydrogen above and below the Galactic Center with the Green Bank Telescope. The observations extend up to Galactic latitude | b | < 10 deg with an effective angular resolution of 9.5' and an average rms brightness temperature noise of 40 mK in a 1 km/s channel. The survey reveals the existence of a population of anomalous high-velocity cl… ▽ More We present the results of a new sensitive survey of neutral hydrogen above and below the Galactic Center with the Green Bank Telescope. The observations extend up to Galactic latitude | b | < 10 deg with an effective angular resolution of 9.5' and an average rms brightness temperature noise of 40 mK in a 1 km/s channel. The survey reveals the existence of a population of anomalous high-velocity clouds extending up to heights of about 1.5 kpc from the Galactic Plane and showing no signature of Galactic rotation. These clouds have local standard of rest velocities | Vlsr | < 360 km/s and, assuming a Galactic Center origin, they have sizes of a few tens of parsecs and neutral hydrogen masses spanning $10-10^5 \, M_\odot$. Accounting for selection effects, the cloud population is symmetric in longitude, latitude, and Vlsr. We model the cloud kinematics in terms of an outflow expanding from the Galactic Center and find the population consistent with being material moving with radial velocity Vw ~ 330 km/s distributed throughout a bi-cone with opening angle $α>140$ deg. This simple model implies an outflow luminosity $Lw > 3 \times 10^{40}$ erg/s over the past 10 Myr, consistent with star formation feedback in the inner region of the Milky Way, with a cold gas mass-loss rate $\lesssim 0.1 \, M_\odot$/yr. These clouds may represent the cold gas component accelerated in the nuclear wind driven by our Galaxy, although some of the derived properties challenge current theoretical models of the entrainment process. △ Less

Submitted 6 February, 2018; originally announced February 2018.

Comments: 16 pages, 9 figures. Accepted for publication in ApJ

Kinematics of the SN Refsdal host revealed by MUSE: a regularly rotating spiral galaxy at z~1.5

Authors: E. M. Di Teodoro, C. Grillo, F. Fraternali, R. Gobat, W. Karman, A. Mercurio, P. Rosati, I. Balestra, G. B. Caminha, K. I. Caputi, M. Lombardi, S. H. Suyu, T. Treu, E. Vanzella

Abstract: We use Multi Unit Spectroscopic Explorer (MUSE) observations of the galaxy cluster MACS J1149.5+2223 to explore the kinematics of the grand-design spiral galaxy Sp1149 hosting the SN Refsdal. Sp1149 lies at $z\simeq1.49$, has a stellar mass $M_*\simeq5\times10^9 \, \mathrm{M_\odot}$, a star-formation rate $\mathrm{SFR} \simeq1-6 \, \mathrm{M_\odot/yr}$ and represents a likely progenitor of a Milky… ▽ More We use Multi Unit Spectroscopic Explorer (MUSE) observations of the galaxy cluster MACS J1149.5+2223 to explore the kinematics of the grand-design spiral galaxy Sp1149 hosting the SN Refsdal. Sp1149 lies at $z\simeq1.49$, has a stellar mass $M_*\simeq5\times10^9 \, \mathrm{M_\odot}$, a star-formation rate $\mathrm{SFR} \simeq1-6 \, \mathrm{M_\odot/yr}$ and represents a likely progenitor of a Milky-Way-like galaxy. All the four multiple images of Sp1149 in our data show strong OII-line emissions pointing to a clear rotation pattern. We take advantage of the gravitational lensing magnification effect ($\simeq 4 \times$) on the OII emission of the least distorted image to fit 3D kinematic models to the MUSE data-cube and derive the rotation curve and the velocity dispersion profile of Sp1149. We find that the rotation curve steeply rises, peaks at $R\simeq1$ kpc and then (initially) declines and flattens to an average $V_\mathrm{flat} = 128^{+29}_{-19}$ km/s. The shape of the rotation curve is well determined but the actual value of $V_\mathrm{flat}$ is quite uncertain because of the nearly face-on configuration of the galaxy. The intrinsic velocity dispersion due to gas turbulence is almost constant across the entire disc with an average of $27\pm5$ km/s. This value is consistent with $z=0$ measurements in the ionized gas component and a factor 2-4 lower than other estimates in different galaxies at similar redshifts. The average stellar-to-total mass fraction is of the order of one fifth. Our kinematic analysis returns the picture of a regular star-forming, mildly turbulent, rotation-dominated ($V / σ\simeq5$) spiral galaxy in a 4 Gyr old Universe. △ Less

Submitted 5 February, 2018; v1 submitted 19 January, 2018; originally announced January 2018.

Comments: 10 pages, 5 figures. Accepted for publication in MNRAS

Galaxy spin as a formation probe: the stellar-to-halo specific angular momentum relation

Authors: Lorenzo Posti, Gabriele Pezzulli, Filippo Fraternali, Enrico di Teodoro

Abstract: We derive the stellar-to-halo specific angular momentum relation (SHSAMR) of galaxies at $z=0$ by combining i) the standard $Λ$CDM tidal torque theory ii) the observed relation between stellar mass and specific angular momentum (Fall relation) and iii) various determinations of the stellar-to-halo mass relation (SHMR). We find that the ratio $f_j = j_\ast/j_{\rm h}$ of the specific angular momentu… ▽ More We derive the stellar-to-halo specific angular momentum relation (SHSAMR) of galaxies at $z=0$ by combining i) the standard $Λ$CDM tidal torque theory ii) the observed relation between stellar mass and specific angular momentum (Fall relation) and iii) various determinations of the stellar-to-halo mass relation (SHMR). We find that the ratio $f_j = j_\ast/j_{\rm h}$ of the specific angular momentum of stars to that of the dark matter i) varies with mass as a double power-law, ii) it always has a peak in the mass range explored and iii) it is $3-5$ times larger for spirals than for ellipticals. The results have some dependence on the adopted SHMR and we provide fitting formulae in each case. For any choice of the SHMR, the peak of $f_j$ occurs at the same mass where the stellar-to-halo mass ratio $f_\ast = M_\ast/M_{\rm h}$ has a maximum. This is mostly driven by the straightness and tightness of the Fall relation, which requires $f_j$ and $f_\ast$ to be correlated with each other roughly as $f_j\propto f_\ast^{2/3}$, as expected if the outer and more angular momentum rich parts of a halo failed to accrete onto the central galaxy and form stars (biased collapse). We also confirm that the difference in the angular momentum of spirals and ellipticals at a given mass is too large to be ascribed only to different spins of the parent dark-matter haloes (spin bias). △ Less

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

Comments: matches MNRAS published version

Journal ref: Monthly Notices of the Royal Astronomical Society, 2018, 475, 232

LITTLE THINGS in 3D: robust determination of the circular velocity of dwarf irregular galaxies

Authors: Giuliano Iorio, Filippo Fraternali, Carlo Nipoti, Enrico Di Teodoro, Justin I. Read, Giuseppina Battaglia

Abstract: Dwarf Irregular galaxies (dIrrs) are the smallest stellar systems with extended HI discs. The study of the kinematics of such discs is a powerful tool to estimate the total matter distribution at these very small scales. In this work, we study the HI kinematics of 17 galaxies extracted from the `Local Irregulars That Trace Luminosity Extremes, The HI Nearby Galaxy Survey' (LITTLE THINGS). Our appr… ▽ More Dwarf Irregular galaxies (dIrrs) are the smallest stellar systems with extended HI discs. The study of the kinematics of such discs is a powerful tool to estimate the total matter distribution at these very small scales. In this work, we study the HI kinematics of 17 galaxies extracted from the `Local Irregulars That Trace Luminosity Extremes, The HI Nearby Galaxy Survey' (LITTLE THINGS). Our approach differs significantly from previous studies in that we directly fit 3D models (two spatial dimensions plus one spectral dimension) using the software $^\text{3D}$BAROLO, fully exploiting the information in the HI datacubes. For each galaxy we derive the geometric parameters of the HI disc (inclination and position angle), the radial distribution of the surface density, the velocity-dispersion ($σ_v$) profile and the rotation curve. The circular velocity (V$_{\text{c}}$), which traces directly the galactic potential, is then obtained by correcting the rotation curve for the asymmetric drift. As an initial application, we show that these dIrrs lie on a baryonic Tully-Fisher relation in excellent agreement with that seen on larger scales. The final products of this work are high-quality, ready-to-use kinematic data ($\textrm{V}_\textrm{c}$ and $σ_v$) that we make publicly available. These can be used to perform dynamical studies and improve our understanding of these low-mass galaxies. △ Less

Submitted 11 November, 2016; originally announced November 2016.

Comments: 36 pages, 28 figures, 2 tables. Submitted to MNRAS (revised version after the referee report). The final rotation curves can be downloaded from http://www.filippofraternali.com/styled-9/index.html

Journal ref: Mon Not R Astron Soc (2017) 466 (4): 4159-4192

Flat rotation curves and low velocity dispersions in KMOS star-forming galaxies at z~1

Authors: E. M. Di Teodoro, F. Fraternali, S. H. Miller

Abstract: The study of the evolution of star-forming galaxies requires the determination of accurate kinematics and scaling relations out to high redshift. In this paper we select a sample of 18 galaxies at z~1, observed in the H-alpha emission-line with KMOS, to derive accurate kinematics using a novel 3D analysis technique. We use the new code 3D-Barolo, that models the galaxy emission directly in the 3D… ▽ More The study of the evolution of star-forming galaxies requires the determination of accurate kinematics and scaling relations out to high redshift. In this paper we select a sample of 18 galaxies at z~1, observed in the H-alpha emission-line with KMOS, to derive accurate kinematics using a novel 3D analysis technique. We use the new code 3D-Barolo, that models the galaxy emission directly in the 3D observational space, without the need to extract kinematic maps. This technique's major advantage is that it is not affected by beam smearing and thus it enables the determination of rotation velocity and intrinsic velocity dispersion, even at low spatial resolution. We find that: 1) the rotation curves of these z~1 galaxies rise very steeply within few kiloparsecs and remain flat out to the outermost radius and 2) the H-alpha velocity dispersions are low, ranging from 15 to 40 km/s, which leads to V/sigma = 3-10. These characteristics are similar to those of disc galaxies in the local Universe. Finally, we also report no significant evolution of the stellar-mass Tully-Fisher relation. Our results show that disc galaxies are kinematically mature and rotation-dominated already at z~1. △ Less

Submitted 18 July, 2016; v1 submitted 16 February, 2016; originally announced February 2016.

Comments: Accepted for publication in A&A, 11 pages, 6 figures, 1 table

Journal ref: A&A 594, A77 (2016)

A highly-ionized region surrounding SN Refsdal revealed by MUSE

Authors: W. Karman, C. Grillo, I. Balestra, P. Rosati, K. I. Caputi, E. Di Teodoro, F. Fraternali, R. Gavazzi, A. Mercurio, J. X. Prochaska, S. Rodney, T. Treu

Abstract: Supernova (SN) Refsdal is the first multiply-imaged, highly-magnified, and spatially-resolved SN ever observed. The SN exploded in a highly-magnified spiral galaxy at z=1.49 behind the Frontier Fields Cluster MACS1149, and provides a unique opportunity to study the environment of SNe at high z. We exploit the time delay between multiple images to determine the properties of the SN and its environm… ▽ More Supernova (SN) Refsdal is the first multiply-imaged, highly-magnified, and spatially-resolved SN ever observed. The SN exploded in a highly-magnified spiral galaxy at z=1.49 behind the Frontier Fields Cluster MACS1149, and provides a unique opportunity to study the environment of SNe at high z. We exploit the time delay between multiple images to determine the properties of the SN and its environment, before, during, and after the SN exploded. We use the integral-field spectrograph MUSE on the VLT to simultaneously target all observed and model-predicted positions of SN Refsdal. We find MgII emission at all positions of SN Refsdal, accompanied by weak FeII* emission at two positions. The measured ratios of [OII] to MgII emission of 10-20 indicate a high degree of ionization with low metallicity. Because the same high degree of ionization is found in all images, and our spatial resolution is too coarse to resolve the region of influence of SN Refsdal, we conclude that this high degree of ionization has been produced by previous SNe or a young and hot stellar population. We find no variability of the [OII] line over a period of 57 days. This suggests that there is no variation in the [OII] luminosity of the SN over this period, or that the SN has a small contribution to the integrated [OII] emission over the scale resolved by our observations. △ Less

Submitted 10 November, 2015; v1 submitted 24 September, 2015; originally announced September 2015.

Comments: 5 pages, 4 figures, accepted for publication in A&A

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

3D-Barolo: a new 3D algorithm to derive rotation curves of galaxies

Authors: Enrico Di Teodoro, Filippo Fraternali

Abstract: We present 3D-Barolo, a new code that derives rotation curves of galaxies from emission-line observations. This software fits 3D tilted-ring models to spectroscopic data-cubes and can be used with a variety of observations: from HI and molecular lines to optical/IR recombination lines. We describe the structure of the main algorithm and show that it performs much better than the standard 2D approa… ▽ More We present 3D-Barolo, a new code that derives rotation curves of galaxies from emission-line observations. This software fits 3D tilted-ring models to spectroscopic data-cubes and can be used with a variety of observations: from HI and molecular lines to optical/IR recombination lines. We describe the structure of the main algorithm and show that it performs much better than the standard 2D approach on velocity fields. A number of successful applications, from high to very low spatial resolution data are presented and discussed. 3D-Barolo can recover the true rotation curve and estimate the intrinsic velocity dispersion even in barely resolved galaxies (about 2 resolution elements) provided that the signal to noise of the data is larger that 2-3. It can also be run automatically thanks to its source-detection and first-estimate modules, which make it suitable for the analysis of large 3D datasets. These features make 3D-Barolo a uniquely useful tool to derive reliable kinematics for both local and high-redshift galaxies from a variety of different instruments including the new-generation IFUs, ALMA and the SKA pathfinders. △ Less

Submitted 28 May, 2015; originally announced May 2015.

Comments: 13 pages + 5 (appendix). Accepted for publication in MNRAS. Code available at http://editeodoro.github.io/Bbarolo

Gas accretion from minor mergers in local spiral galaxies

Authors: Enrico Di Teodoro, Filippo Fraternali

Abstract: In this paper we quantify the gas accretion rate from minor mergers onto star-forming galaxies in the Local Universe using HI observations of 148 nearby spiral galaxies (WHISP sample). We developed a dedicated code that iteratively analyses HI data-cubes, finds dwarf gas-rich satellites around larger galaxies and estimates an upper limit to the gas accretion rate. We found that 22% of the galaxies… ▽ More In this paper we quantify the gas accretion rate from minor mergers onto star-forming galaxies in the Local Universe using HI observations of 148 nearby spiral galaxies (WHISP sample). We developed a dedicated code that iteratively analyses HI data-cubes, finds dwarf gas-rich satellites around larger galaxies and estimates an upper limit to the gas accretion rate. We found that 22% of the galaxies have at least one detected dwarf companion. We made the very stringent assumption that all satellites are going to merge in the shortest possible time transferring all their gas to the main galaxies. This leads to an estimate of the maximum gas accretion rate of 0.28 solar masses per year, about five times lower than the average SFR of the sample. Given the assumptions, our accretion rate is clearly an overestimate. Our result strongly suggests that minor mergers do not play a significant role in the total gas accretion budget in local galaxies. △ Less

Submitted 3 June, 2014; originally announced June 2014.

Comments: 15 pages, 5 figures, 2 tables

Journal ref: A&A 567, A68 (2014)