Showing 1–36 of 36 results for author: Tepper-García, T

Nexus: A framework for simulations of idealised galaxies

Authors: Thor Tepper-Garcia, Joss Bland-Hawthorn, Eugene Vasiliev, Oscar Agertz, Romain Teyssier, Christoph Federrath

Abstract: Motivated by the need for realistic, dynamically self-consistent, evolving galaxy models that avoid the inherent complexity of full, and zoom-in, cosmological simulations, we have developed Nexus, an integral, flexible framework to create synthetic galaxies made of both collisionless and gaseous components. Nexus leverages the power of publicly available, tried-and-tested packages: i) the stellar-… ▽ More Motivated by the need for realistic, dynamically self-consistent, evolving galaxy models that avoid the inherent complexity of full, and zoom-in, cosmological simulations, we have developed Nexus, an integral, flexible framework to create synthetic galaxies made of both collisionless and gaseous components. Nexus leverages the power of publicly available, tried-and-tested packages: i) the stellar-dynamics, action-based library AGAMA; and ii) the Adaptive Mesh Refinement, N-body/hydrodynamical code Ramses, modified to meet our needs. In addition, we make use of a proprietary module to account for realistic galaxy formation (sub-grid) physics, including star formation, stellar feedback, and chemical enrichment. As a framework to perform controlled experiments of idealised galaxies, Nexus' basic functionality consists in the generation of bespoke initial conditions (ICs) for any desired galaxy model, which are advanced in time to simulate the system's evolution. The fully self-consistent ICs are generated with a distribution-function based approach, as implemented in the galaxy modelling module of AGAMA -- up to now restricted to collisionless components, extended in this work to treat two types of gaseous configurations: (i) hot halos; and (ii) gas disks. For the first time, we are able to construct equilibrium models with disc gas fractions in the range $0 < f_{\rm gas} < 1$, needed to model high-redshift galaxies. The framework is ideally suited to the study of galactic ecology, specifically how stars and gas work together over billions of years. As a validation of our framework, we reproduce - and improve upon - several isolated galaxy model setups reported in earlier studies. Finally, we showcase Nexus by presenting an interesting type of `nested bar' galaxy class. Future upgrades of Nexus will include magneto-hydrodynamics and highly energetic particle (`cosmic ray') heating. △ Less

Submitted 1 June, 2024; originally announced June 2024.

Comments: To be submitted to MNRAS; we welcome general comments and pointers to relevant, missing references

Discovery of a dormant 33 solar-mass black hole in pre-release Gaia astrometry

Authors: Gaia Collaboration, P. Panuzzo, T. Mazeh, F. Arenou, B. Holl, E. Caffau, A. Jorissen, C. Babusiaux, P. Gavras, J. Sahlmann, U. Bastian, Ł. Wyrzykowski, L. Eyer, N. Leclerc, N. Bauchet, A. Bombrun, N. Mowlavi, G. M. Seabroke, D. Teyssier, E. Balbinot, A. Helmi, A. G. A. Brown, A. Vallenari, T. Prusti, J. H. J. de Bruijne , et al. (390 additional authors not shown)

Abstract: Gravitational waves from black-hole merging events have revealed a population of extra-galactic BHs residing in short-period binaries with masses that are higher than expected based on most stellar evolution models - and also higher than known stellar-origin black holes in our Galaxy. It has been proposed that those high-mass BHs are the remnants of massive metal-poor stars. Gaia astrometry is exp… ▽ More Gravitational waves from black-hole merging events have revealed a population of extra-galactic BHs residing in short-period binaries with masses that are higher than expected based on most stellar evolution models - and also higher than known stellar-origin black holes in our Galaxy. It has been proposed that those high-mass BHs are the remnants of massive metal-poor stars. Gaia astrometry is expected to uncover many Galactic wide-binary systems containing dormant BHs, which may not have been detected before. The study of this population will provide new information on the BH-mass distribution in binaries and shed light on their formation mechanisms and progenitors. As part of the validation efforts in preparation for the fourth Gaia data release (DR4), we analysed the preliminary astrometric binary solutions, obtained by the Gaia Non-Single Star pipeline, to verify their significance and to minimise false-detection rates in high-mass-function orbital solutions. The astrometric binary solution of one source, Gaia BH3, implies the presence of a 32.70 \pm 0.82 M\odot BH in a binary system with a period of 11.6 yr. Gaia radial velocities independently validate the astrometric orbit. Broad-band photometric and spectroscopic data show that the visible component is an old, very metal-poor giant of the Galactic halo, at a distance of 590 pc. The BH in the Gaia BH3 system is more massive than any other Galactic stellar-origin BH known thus far. The low metallicity of the star companion supports the scenario that metal-poor massive stars are progenitors of the high-mass BHs detected by gravitational-wave telescopes. The Galactic orbit of the system and its metallicity indicate that it might belong to the Sequoia halo substructure. Alternatively, and more plausibly, it could belong to the ED-2 stream, which likely originated from a globular cluster that had been disrupted by the Milky Way. △ Less

Submitted 19 April, 2024; v1 submitted 16 April, 2024; originally announced April 2024.

Comments: 23 pages, accepted fro publication in A&A Letters. New version with small fixes

Turbulent gas-rich disks at high redshift: bars & bulges in a radial shear flow

Authors: Joss Bland-Hawthorn, Thor Tepper-Garcia, Oscar Agertz, Christoph Federrath

Abstract: Recent observations of high-redshift galaxies ($z \lesssim 7$) reveal that a substantial fraction have turbulent, gas-rich disks with well-ordered rotation and elevated levels of star formation. In some instances, disks show evidence of spiral arms, with bar-like structures. These remarkable observations have encouraged us to explore a new class of dynamically self-consistent models using our hydr… ▽ More Recent observations of high-redshift galaxies ($z \lesssim 7$) reveal that a substantial fraction have turbulent, gas-rich disks with well-ordered rotation and elevated levels of star formation. In some instances, disks show evidence of spiral arms, with bar-like structures. These remarkable observations have encouraged us to explore a new class of dynamically self-consistent models using our hydrodynamic N-body simulation framework that mimic a plausible progenitor of the Milky Way at high redshift. We explore disk gas fractions of $f_{\rm gas} = 0, 20, 40, 60, 80, 100\%$ and track the creation of stars and metals. The high gas surface densities encourage vigorous star formation, which in turn couples with the gas to drive turbulence. We explore three distinct histories: (i) there is no ongoing accretion and the gas is used up by the star formation; (ii) the star-forming gas is replenished by cooling in the hot halo gas; (iii) in a companion paper, we revisit the models in the presence of a strong perturbing force. At low $f_{\rm disk}$ ($<0.3$), where $f_{\rm disk}$ is the mass fraction of stars relative to dark matter within 2.2 $R_{\rm disk}$, a bar does not form in a stellar disk; this remains true even when gas dominates the inner disk potential. For a dominant baryon disk ($f_{\rm disk} \gtrsim 0.5$) at all gas fractions, the turbulent gas forms a strong "radial shear flow" that leads to an intermittent star-forming bar within about 500 Myr; turbulent gas speeds up the formation of bars compared to gas-free models. For $f_{\rm gas} \lesssim 60\%$, all bars survive, but for higher gas fractions, the bar devolves into a central bulge after 1 Gyr. The star-forming bars are reminiscent of recent discoveries in high-redshift ALMA observations of gaseous disks. △ Less

Submitted 1 May, 2024; v1 submitted 8 February, 2024; originally announced February 2024.

Comments: 40 pages, 21 figures, ApJ, in press (April 2024) - simulation movies at http://www.physics.usyd.edu.au/turbo_disks

Galactic-Seismology Substructures and Streams Hunter with LAMOST and Gaia. I. Methodology and Local Halo Results

Authors: Guan-Yu Wang, Hai-Feng Wang, Yang-Ping Luo, Yuan-Sen Ting, Thor Tepper-García, Joss Bland-Hawthorn, Jeffrey Carlin

Abstract: We present a novel, deep-learning based method -- dubbed Galactic-Seismology Substructures and Streams Hunter, or GS$^{3}$ Hunter for short, to search for substructures and streams in stellar kinematics data. GS$^{3}$ Hunter relies on a combined application of Siamese Neural Networks to transform the phase space information and the K-means algorithm for the clustering. As a validation test, we app… ▽ More We present a novel, deep-learning based method -- dubbed Galactic-Seismology Substructures and Streams Hunter, or GS$^{3}$ Hunter for short, to search for substructures and streams in stellar kinematics data. GS$^{3}$ Hunter relies on a combined application of Siamese Neural Networks to transform the phase space information and the K-means algorithm for the clustering. As a validation test, we apply GS$^{3}$ Hunter to a subset of the Feedback in Realistic Environments (FIRE) cosmological simulations. The stellar streams and substructures thus identified are in good agreement with corresponding results reported earlier by the FIRE team. In the same vein, we apply our method to a subset of local halo stars from the Gaia Early Data Release 3 and GALAH DR3 datasets, and recover several, previously known dynamical groups, such as Thamnos 1+2, Hot Thick Disk, ED-1, L-RL3, Helmi 1+2, and Gaia-Sausage-Enceladus, Sequoia, VRM, Cronus, Nereus. Finally, we apply our method without fine-tuning to a subset of K-giant stars located in the inner halo region, obtained from the LAMOST Data Release 5 (DR5) dataset. We recover three, previously known structures (Sagittarius, Hercules-Aquila Cloud, and the Virgo Overdensity), but we also discover a number of new substructures. We anticipate that GS$^{3}$ Hunter will become a useful tool for the community dedicated to the search of stellar streams and structures in the Milky Way (MW) and the Local group, thus helping advance our understanding of the stellar inner and outer halos, and of the assembly and tidal stripping history in and around the MW. △ Less

Submitted 10 January, 2024; originally announced January 2024.

Comments: 40 pages, 33 figures, 3 tables

The far side of the Galactic bar/bulge revealed through semi-regular variables

Authors: Daniel R. Hey, Daniel Huber, Benjamin J. Shappee, Joss Bland-Hawthorn, Thor Tepper-García, Robyn Sanderson, Sukanya Chakrabarti, Nicholas Saunders, Jason A. S. Hunt, Timothy R. Bedding, John Tonry

Abstract: The Galactic bulge and bar are critical to our understanding of the Milky Way. However, due to the lack of reliable stellar distances, the structure and kinematics of the bulge/bar beyond the Galactic center have remained largely unexplored. Here, we present a method to measure distances of luminous red giants using a period-amplitude-luminosity relation anchored to the Large Magellanic Cloud, wit… ▽ More The Galactic bulge and bar are critical to our understanding of the Milky Way. However, due to the lack of reliable stellar distances, the structure and kinematics of the bulge/bar beyond the Galactic center have remained largely unexplored. Here, we present a method to measure distances of luminous red giants using a period-amplitude-luminosity relation anchored to the Large Magellanic Cloud, with random uncertainties of 10-15% and systematic errors below 1-2%. We apply this method to data from the Optical Gravitational Lensing Experiment (OGLE) to measure distances to $190,302$ stars in the Galactic bulge and beyond out to 20 kpc. Using this sample we measure a distance to the Galactic center of $R_0$ = $8108\pm106_{\rm stat}\pm93_{\rm sys}$ pc, consistent with astrometric monitoring of stars orbiting Sgr A*. We cross-match our distance catalog with Gaia DR3 and use the subset of $39,566$ overlapping stars to provide the first constraints on the Milky Way's velocity field ($V_R,V_φ,V_z$) beyond the Galactic center. We show that the $V_R$ quadrupole from the bar's near side is reflected with respect to the Galactic center, indicating that the bar is both bi-symmetric and aligned with the inner disk, and therefore dynamically settled along its full extent. We also find that the vertical height $V_Z$ map has no major structure in the region of the Galactic bulge, which is inconsistent with a current episode of bar buckling. Finally, we demonstrate with N-body simulations that distance uncertainty plays a major factor in the alignment of the major and kinematic axes of the bar and distribution of velocities, necessitating caution when interpreting results for distant stars. △ Less

Submitted 21 November, 2023; v1 submitted 30 May, 2023; originally announced May 2023.

Comments: Accepted to the Astronomical Journal

Prospective Dark Matter Annihilation Signals From the Sagittarius Dwarf Spheroidal

Authors: Thomas A. A. Venville, Alan R. Duffy, Roland M. Crocker, Oscar Macias, Thor Tepper-García

Abstract: The Sagittarius Dwarf Spheroidal galaxy (Sgr) is investigated as a target for DM annihilation searches utilising J-factor distributions calculated directly from a high-resolution hydrodynamic simulation of the infall and tidal disruption of Sgr around the Milky Way. In contrast to past studies, the simulation incorporates DM, stellar and gaseous components for both the Milky Way and the Sgr progen… ▽ More The Sagittarius Dwarf Spheroidal galaxy (Sgr) is investigated as a target for DM annihilation searches utilising J-factor distributions calculated directly from a high-resolution hydrodynamic simulation of the infall and tidal disruption of Sgr around the Milky Way. In contrast to past studies, the simulation incorporates DM, stellar and gaseous components for both the Milky Way and the Sgr progenitor galaxy. The simulated distributions account for significant tidal disruption affecting the DM density profile. Our estimate of the J-factor value for Sgr, $J_{\text{Sgr}}=1.48\times 10^{10}$ M$_\odot^2$ kpc$^{-5}$ ($6.46\times10^{16}\ \text{GeV}\ \text{cm}^{-5}$), is significantly lower than found in prior studies. This value, while formally a lower limit, is likely close to the true J-factor value for Sgr. It implies a DM cross-section incompatibly large in comparison with existing constraints would be required to attribute recently observed $γ$-ray emission from Sgr to DM annihilation. We also calculate a J-factor value using a NFW profile fitted to the simulated DM density distribution to facilitate comparison with past studies. This NFW J-factor value supports the conclusion that most past studies have overestimated the dark matter density of Sgr on small scales. This, together with the fact that the Sgr has recently been shown to emit $γ$-rays of astrophysical origin, complicate the use of Sgr in indirect DM detection searches. △ Less

Submitted 13 November, 2023; v1 submitted 23 March, 2023; originally announced March 2023.

Comments: 15 pages, 15 figures. Accepted for publication in MNRAS. Accepted manuscript uploaded

The rapid onset of stellar bars in the baryon-dominated centers of disk galaxies

Authors: J. Bland-Hawthorn, T. Tepper-Garcia, O. Agertz, K. Freeman

Abstract: Recent observations of high-redshift galactic disks ($z\approx 1-3$) show a strong negative trend in the dark matter fraction $f_{DM}$ with increasing baryonic surface density. For this to be true, the inner baryons must dominate over dark matter in early massive galaxies, as observed in the Milky Way today. If disks are dominant at early times, we show that stellar bars form promptly within these… ▽ More Recent observations of high-redshift galactic disks ($z\approx 1-3$) show a strong negative trend in the dark matter fraction $f_{DM}$ with increasing baryonic surface density. For this to be true, the inner baryons must dominate over dark matter in early massive galaxies, as observed in the Milky Way today. If disks are dominant at early times, we show that stellar bars form promptly within these disks, leading to a high bar fraction at early times. New JWST observations provide the best evidence to date for mature stellar bars in this redshift range. The disk mass fraction $f_{disk}$ within $R_s=2.2 R_{disk}$ is the dominant factor in determining how rapidly a bar forms. Using 3D hydro simulations of halo-disk-bulge galaxies, we confirm the "Fujii relation" for the exponential dependence of the bar formation time $τ_{bar}$ as a function of $f_{disk}$. For $f_{disk} > 0.3$, the bar formation time declines exponentially fast with increasing $f_{disk}$. This relation is a challenge to simulators - barred models with inadequate resolution fall off this curve. Instead of Fujii's arbitrary threshold for when a bar forms, for the first time, we exploit the exponential growth timescale associated with a positive feedback cycle as the bar emerges from the underlying disk. A modified, mass-dependent trend is observed for halos relevant to systems at cosmic noon ($10.5 < \log M_{halo} < 12$), where the bar onset is slower for higher mass halos at a fixed $f_{disk}$. If baryons dominate over dark matter within $R \approx R_s$, we predict that a high fraction of bars will be found in high-redshift disks long before $z = 1$. Due to its widespread use in simulations, we investigate the Efstathiou-Lake-Negroponte criterion for bar instability: this sub-optimal parameter is inversely related to $f_{disk}$, with a secondary dependence on $M_{halo}$. △ Less

Submitted 13 March, 2023; v1 submitted 9 March, 2023; originally announced March 2023.

Comments: 27 pages, 8 figures, 1 table - Astrophysical Journal, accepted (9 March 2023)

A new resonance-like feature in the outer disc of the Milky Way

Authors: Ronald Drimmel, Shourya Khanna, Elena D'Onghia, Thorsten Tepper-García, Joss Bland-Hawthorn, Laurent Chemin, Vincenzo Ripepi, Mercé Romero-Gómez, Pau Ramos, Eloisa Poggio, Rene Andrae, Ronny Blomme, Tristan Cantat-Gaudin, Alfred Castro-Ginard, Gisella Clementini, Francesca Fiqueras, Yves Frémat, Morgan Fouesneau, Alex Lobel, Douglas Marshall, Tatiana Muraveva

Abstract: Modern astrometric and spectroscopic surveys have revealed a wealth of structure in the phase space of stars in the Milky Way, with evidence of resonance features and non-equilibrium processes. Using Gaia's third data release, we present evidence of a new resonance-like feature in the outer disc of the Milky Way. The feature is most evident in the angular momentum distribution of the young Classic… ▽ More Modern astrometric and spectroscopic surveys have revealed a wealth of structure in the phase space of stars in the Milky Way, with evidence of resonance features and non-equilibrium processes. Using Gaia's third data release, we present evidence of a new resonance-like feature in the outer disc of the Milky Way. The feature is most evident in the angular momentum distribution of the young Classical Cepheids, a population for which we can derive accurate distances over much of the Galactic disc. We then search for similar features in the outer disc using a much larger sample of red giant stars, as well as a compiled list of over 31 million stars with spectroscopic line-of-sight velocity measurements. While much less evident in these two older samples, the distribution of stars in action-configuration space suggests that resonance features are present here as well. The position of the feature in action-configuration space suggests that the new feature may be related to the Galactic bar, but other possibilities are discussed. △ Less

Submitted 20 January, 2023; v1 submitted 26 July, 2022; originally announced July 2022.

Comments: 13 pages, 11 figures, in press for A&A

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

The disturbed outer Milky Way disc

Authors: Paul J. McMillan, Jonathan Petersson, Thor Tepper-Garcia, Joss Bland-Hawthorn, Teresa Antoja, Laurent Chemin, Francesca Figueras, Shourya Khanna, Georges Kordopatis, Pau Ramos, Merce Romero-Gómez, George Seabroke

Abstract: The outer parts of the Milky Way's disc are significantly out of equilibrium. Using only distances and proper motions of stars from Gaia's Early Data Release 3, in the range |b|<10°, 130°<l<230°, we show that for stars in the disc between around 10 and 14 kpc from the Galactic centre, vertical velocity is strongly dependent on the angular momentum, azimuth, and position above or below the Galactic… ▽ More The outer parts of the Milky Way's disc are significantly out of equilibrium. Using only distances and proper motions of stars from Gaia's Early Data Release 3, in the range |b|<10°, 130°<l<230°, we show that for stars in the disc between around 10 and 14 kpc from the Galactic centre, vertical velocity is strongly dependent on the angular momentum, azimuth, and position above or below the Galactic plane. We further show how this behaviour translates into a bimodality in the velocity distribution of stars in the outer Milky Way disc. We use an N-body model of an impulse-like interaction of the Milky Way disc with a perturber similar to the Sagittarius dwarf to demonstrate that this mechanism can generate a similar disturbance. It has already been shown that this interaction can produce a phase spiral similar to that seen in the Solar neighbourhood. We argue that the details of this substructure in the outer galaxy will be highly sensitive to the timing of the perturbation or the gravitational potential of the Galaxy, and therefore may be key to disentangling the history and structure of the Milky Way. △ Less

Submitted 22 September, 2022; v1 submitted 8 June, 2022; originally announced June 2022.

Comments: 15 pages, MNRAS accepted version. 7 minute video summary available at https://youtu.be/6zQvMivZONk

Kinematic-Chemical analysis and Time tagging for the Diagonal Ridge Structure of the Galactic Outer Disk with LAMOST Red Giant Branch Stars

Authors: Peng Yang, Hai-Feng Wang, Zhi-Quan Luo, Thor Tepper-García, Yang-Ping Luo

Abstract: We investigate the kinematic-chemical distribution of Red Giant Branch (RGB) stars from the LAMOST survey crossed matched with Gaia DR2 proper motions, and present time tagging for the well-known ridge structures (diagonal distributions for $V_R$ in the $R$, $V_φ$ plane) in the range of Galactocentric distance $R$ = 8 to 15 kpc. We detect six ridge structures, including five ridges apparent in the… ▽ More We investigate the kinematic-chemical distribution of Red Giant Branch (RGB) stars from the LAMOST survey crossed matched with Gaia DR2 proper motions, and present time tagging for the well-known ridge structures (diagonal distributions for $V_R$ in the $R$, $V_φ$ plane) in the range of Galactocentric distance $R$ = 8 to 15 kpc. We detect six ridge structures, including five ridges apparent in the radial velocity distribution and three ridges apparent in the vertical velocity, the sensitive time of which to the perturbations are from young population (0$-$3 Gyr) to old population (9$-$14 Gyr). Based on an analysis of the evolution of angular momentum distribution, we find that four ridges are relatively stationary, while another is evolving with time, which is confirmed by the difference analysis at different populations and supporting that there might be two kinds of dynamical origins. Furthermore, ridge features are also vividly present in the chemical properties ([Fe/H], [$α$/Fe]). The comparison between the north and south hemispheres of the Galaxy does show some differences and the ridge features are asymmetrical. Moreover, we find that diagonal ridge structures may affect the shape of the rotation curve, which is manifested as fluctuations and undulations on top of a smooth profile. Finally we speculate that the bar dynamics should be not enough to explain all ridge properties including the break feature in the $V_Z$-$L_Z$ plane. △ Less

Submitted 21 December, 2022; v1 submitted 18 May, 2022; originally announced May 2022.

Comments: 17 pages, 14 figures. Accepted publication in The Astronomical Journal

Galactic seismology: joint evolution of impact-triggered stellar and gaseous disc corrugations

Authors: Thor Tepper-Garcia, Joss Bland-Hawthorn, Ken Freeman

Abstract: Evidence for wave-like corrugations are well established in the Milky Way and in nearby disc galaxies. These were originally detected as a displacement of the interstellar medium about the midplane, either in terms of vertical distance or vertical velocity. Over the past decade, similar patterns have emerged in the Milky Way's stellar disc. We investigate how these vertical waves are triggered by… ▽ More Evidence for wave-like corrugations are well established in the Milky Way and in nearby disc galaxies. These were originally detected as a displacement of the interstellar medium about the midplane, either in terms of vertical distance or vertical velocity. Over the past decade, similar patterns have emerged in the Milky Way's stellar disc. We investigate how these vertical waves are triggered by a passing satellite. Using high-resolution N-body/hydrodynamical simulations, we systematically study how the corrugations set up and evolve jointly in the stellar and gaseous discs. We find that the gas corrugations follow the stellar corrugations, i.e. they are initially in phase although, after a few rotation periods (500-700 Myr), the distinct waves separate and thereafter evolve in different ways. The spatial and kinematic amplitudes (and thus the energy) of the corrugations dampen with time, with the gaseous corrugation settling at a faster rate (~800 Myr versus ~1 Gyr). In contrast, the vertical energy of individual disc stars is fairly constant throughout the galaxy's evolution. This difference arises because corrugations are an emergent phenomenon supported by the collective, ordered motions of co-spatial ensembles of stars. We show that the damping of the stellar corrugations can be understood as a consequence of incomplete phase mixing, while the damping of the gaseous corrugations is a natural consequence of the dissipative nature of the gas. We suggest that - in the absence of further, strong perturbations - the degree of correlation between the stellar and gaseous waves may help to age-date the phenomenon. △ Less

Submitted 26 July, 2022; v1 submitted 26 April, 2022; originally announced April 2022.

Comments: MNRAS, in press. Title slightly shortened. References added. Matches accepted version

The role of the halo magnetic field on accretion through High-Velocity Clouds

Authors: Asger Grønnow, Thor Tepper-García, Joss Bland-Hawthorn, Filippo Fraternali

Abstract: High-Velocity Clouds (HVCs) are believed to be an important source of gas accretion for star formation in the Milky Way. Earlier numerical studies have found that the Galactic magnetic field and radiative cooling strongly affects accretion. However, these effects have not previously been included together in the context of clouds falling through the Milky Way's gravitational potential. We explore… ▽ More High-Velocity Clouds (HVCs) are believed to be an important source of gas accretion for star formation in the Milky Way. Earlier numerical studies have found that the Galactic magnetic field and radiative cooling strongly affects accretion. However, these effects have not previously been included together in the context of clouds falling through the Milky Way's gravitational potential. We explore this by simulating an initially stationary cloud falling through the hot hydrostatic corona towards the disc. This represents a HVC that has condensed out of the corona. We include the magnetic field in the corona to examine its effect on accretion of the HVC and its associated cold gas. Remnants of the original cloud survive in all cases, although a strong magnetic field causes it to split into several fragments. We find that mixing of cold and hot gas leads to cooling of coronal gas and an overall growth with time in cold gas mass, despite the low metallicity of the cloud and corona. The role of the magnetic field is to (moderately to severely) suppress the mixing and subsequent cooling, which in turn leads to less accretion compared to when the field is absent. A stronger field leads to less suppression of condensation because it enhances Rayleigh-Taylor instability. However, magnetic tension in a stronger field substantially decelerates condensed cloudlets. These have velocities typically a factor 3-8 below the velocity of the main cloud remnants by the end of the simulation. Some of these cloudlets likely disperse before reaching the disc. △ Less

Submitted 24 November, 2021; originally announced November 2021.

Comments: 15 pages, 14 figures. Accepted for publication in MNRAS

A barred Milky Way surrogate from an N-body simulation

Authors: T. Tepper-Garcia, J. Bland-Hawthorn, E. Vasiliev, E. Athanassoula, O. Gerhard, A. Quillen, P. McMillan, K. Freeman, G. F. Lewis, R. Teyssier, S. Sharma, M. R. Hayden, S. Buder

Abstract: We present an N-body model for the barred Milky Way (MW) galaxy that reproduces many of its properties, including the overall mass distribution, the disc kinematics, and the properties of the central bar. Our high-resolution (N ~ 10^8 particles) simulation, performed with the Ramses code, starts from an axisymmetric non-equilibrium configuration constructed within the AGAMA framework. This is a se… ▽ More We present an N-body model for the barred Milky Way (MW) galaxy that reproduces many of its properties, including the overall mass distribution, the disc kinematics, and the properties of the central bar. Our high-resolution (N ~ 10^8 particles) simulation, performed with the Ramses code, starts from an axisymmetric non-equilibrium configuration constructed within the AGAMA framework. This is a self-consistent dynamical model of the MW defined by the best available parameters for the dark matter halo, the stellar disc and the bulge. For the known (stellar and gas) disc mass (4.5 x 10^10 Msun) and disc mass fraction at R ~ 2.2 R_d (f_d ~ 0.3 - 0.6), the low mass limit does not yield a bar in a Hubble time. The high mass limit adopted here produces a box/peanut bar within about 2 Gyr with the correct mass (~10^10 Msun), size (~5 kpc) and peak pattern speed (~ 40-45 km/s/kpc). In agreement with earlier work, the bar formation timescale scales inversely with f_d (i.e. log [T/Gyr] ~ 0.60/f_d - 0.83 for 1 < f_d < 0.3). The disc radial heating is strong, but, in contrast to earlier claims, we find that disc vertical heating outside of the box/peanut bulge structure is negligible. The synthetic barred MW exhibits long-term stability, except for the slow decline (roughly -2 km/s/kpc/Gyr) of the bar pattern speed, consistent with recent estimates. If our model is indicative of the Milky Way, we estimate that the bar first emerged 3-4 Gyr ago. △ Less

Submitted 9 November, 2021; originally announced November 2021.

Comments: We welcome comments and missing references

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

The GALAH Survey: Chemical Clocks

Authors: Michael R. Hayden, Sanjib Sharma, Joss Bland-Hawthorn, Lorenzo Spina, Sven Buder, Martin Asplund, Andrew R. Casey, Gayandhi M. De Silva, Valentina D'Orazi, Ken C. Freeman, Janez Kos, Geraint F. Lewis, Jane Lin, Karin Lind, Sarah L. Martell, Katharine J. Schlesinger, Jeffrey D. Simpson, Daniel B. Zucker, Tomaz Zwitter, Boquan Chen, Klemen Cotar, Diane Feuillet, Jonti Horner, Meridith Joyce, Thomas Nordlander , et al. (5 additional authors not shown)

Abstract: Previous studies have found that the elemental abundances of a star correlate directly with its age and metallicity. Using this knowledge, we derive ages for a sample of 250,000 stars taken from GALAH DR3 using only their overall metallicity and chemical abundances. Stellar ages are estimated via the machine learning algorithm $XGBoost$, using main sequence turnoff stars with precise ages as our i… ▽ More Previous studies have found that the elemental abundances of a star correlate directly with its age and metallicity. Using this knowledge, we derive ages for a sample of 250,000 stars taken from GALAH DR3 using only their overall metallicity and chemical abundances. Stellar ages are estimated via the machine learning algorithm $XGBoost$, using main sequence turnoff stars with precise ages as our input training set. We find that the stellar ages for the bulk of the GALAH DR3 sample are accurate to 1-2 Gyr using this method. With these ages, we replicate many recent results on the age-kinematic trends of the nearby disk, including the age-velocity dispersion relationship of the solar neighborhood and the larger global velocity dispersion relations of the disk found using $Gaia$ and GALAH. The fact that chemical abundances alone can be used to determine a reliable age for a star have profound implications for the future study of the Galaxy as well as upcoming spectroscopic surveys. These results show that the chemical abundance variation at a given birth radius is quite small, and imply that strong chemical tagging of stars directly to birth clusters may prove difficult with our current elemental abundance precision. Our results highlight the need of spectroscopic surveys to deliver precision abundances for as many nucleosynthetic production sites as possible in order to estimate reliable ages for stars directly from their chemical abundances. Applying the methods outlined in this paper opens a new door into studies of the kinematic structure and evolution of the disk, as ages may potentially be estimated for a large fraction of stars in existing spectroscopic surveys. This would yield a sample of millions of stars with reliable age determinations, and allow precise constraints to be put on various kinematic processes in the disk, such as the efficiency and timescales of radial migration. △ Less

Submitted 27 November, 2020; originally announced November 2020.

Comments: 13 pages, 15 figures, submitted to MNRAS

Combined APOGEE-GALAH stellar catalogues using the Cannon

Authors: Govind Nandakumar, Michael R. Hayden, Sanjib Sharma, Sven Buder, Martin Asplund, Joss Bland-Hawthorn, Gayandhi M. De Silva, Valentina D'Orazi, Ken C. Freeman, Janez Kos, Geraint F. Lewis, Sarah L. Martell, Katharine J. Schlesinger, Jane Lin, Jeffrey D. Simpson, Daniel B. Zucker, Tomaz Zwitter, Thomas Nordlander, Luca Casagrande, Karin Lind, Klemen Cotar, Dennis Stello, Robert A. Wittenmyer, Thor Tepper-Garcia

Abstract: APOGEE and GALAH are two high resolution multi-object spectroscopic surveys that provide fundamental stellar parameters and multiple elemental abundance estimates for about half a million stars in the Milky Way. Both surveys observe in different wavelength regimes and use different data reduction pipelines leading to significant offsets and trends in stellar parameters and abundances for the commo… ▽ More APOGEE and GALAH are two high resolution multi-object spectroscopic surveys that provide fundamental stellar parameters and multiple elemental abundance estimates for about half a million stars in the Milky Way. Both surveys observe in different wavelength regimes and use different data reduction pipelines leading to significant offsets and trends in stellar parameters and abundances for the common stars observed in both surveys. Such systematic differences/offsets in stellar parameters and abundances make it difficult to effectively utilise them to investigate Galactic abundance trends in spite of the unique advantage provided by their complementary sky coverage and different Milky Way components they observe. Hence, we use the \textit{Cannon} data-driven method selecting a training set of 4418 common stars observed by both surveys. This enables the construction of two catalogues, one with the APOGEE scaled and the other with the GALAH scaled stellar parameters. Using repeat observations in APOGEE and GALAH, we find high precision in metallicity (~ 0.02-0.4 dex) and alpha abundances (~ 0.02-0.03 dex) for spectra with good signal-to-noise ratio (SNR > 80 for APOGEE, SNR > 40 for GALAH). We use open and globular clusters to validate our parameter estimates and find small scatter in metallicity (0.06 dex) and alpha abundances (0.03 dex) in APOGEE scaled case. The final catalogues have been cross matched with the Gaia EDR3 catalogue to enable their use to carry out detailed chemo-dynamic studies of the Milky Way from perspectives of APOGEE and GALAH. △ Less

Submitted 29 March, 2022; v1 submitted 5 November, 2020; originally announced November 2020.

Comments: 24 pages, 20 figures, Accepted for publication in MNRAS after major revision of the previous version of the manuscript titled "The GALAH survey: Milky Way disc metallicity and alpha-abundance trends in combined APOGEE-GALAH catalogues"

The GALAH Survey: Accreted stars also inhabit the Spite Plateau

Authors: Jeffrey D. Simpson, Sarah L. Martell, Sven Buder, Sanjib Sharma, Martin Asplund, Joss Bland-Hawthorn, Andrew R. Casey, Gayandhi M. De Silva, Valentina D'Orazi, Ken C. Freeman, Michael Hayden, Janez Kos, Geraint F. Lewis, Karin Lind, Dennis Stello, Daniel B. Zucker, Tomaž Zwitter, Katharine J. Schlesinger, Yuan-Sen Ting, Thomas Nordlander, Gary Da Costa, Klemen Čotar, Jonathan Horner, Thor Tepper-García, The GALAH Collaboration

Abstract: The ESA Gaia astrometric mission has enabled the remarkable discovery that a large fraction of the stars near the Solar neighbourhood appear to be debris from a single in-falling system, the so-called Gaia-Enceladus-Sausage (GSE). One exciting feature of this result is that it gives astronomers for the first time a large sample of easily observable unevolved stars that formed in an extra-Galactic… ▽ More The ESA Gaia astrometric mission has enabled the remarkable discovery that a large fraction of the stars near the Solar neighbourhood appear to be debris from a single in-falling system, the so-called Gaia-Enceladus-Sausage (GSE). One exciting feature of this result is that it gives astronomers for the first time a large sample of easily observable unevolved stars that formed in an extra-Galactic environment, which can be compared to stars that formed within our Milky Way. Here we use these stars to investigate the "Spite Plateau" -- the near-constant lithium abundance observed in metal-poor dwarf stars across a wide range of metallicities (-3<[Fe/H]<-1). In particular our aim is to test whether the stars that formed in the GSE show a different Spite Plateau to other Milky Way stars that inhabit the disk and halo. Individual galaxies could have different Spite Plateaus -- e.g., the ISM could be more depleted in lithium in a lower galactic mass system due to it having a smaller reservoir of gas. We identified 76 GSE dwarf stars observed and analyzed by the GALactic Archeology with HERMES (GALAH) survey as part of its Third Data Release. Orbital actions were used to select samples of Gaia-Enceladus stars, and comparison samples of halo and disk stars. We find that the Gaia-Enceladus stars show the same lithium abundance as other likely accreted stars and in situ Milky Way stars, strongly suggesting that the "lithium problem" is not a consequence of the formation environment. This result fits within the growing consensus that the Spite Plateau, and more generally the "cosmological lithium problem" -- the observed discrepancy between the amount of lithium in warm, metal-poor dwarf stars in our Galaxy, and the amount of lithium predicted to have been produced by Big Bang Nucleosynthesis -- is the result of lithium depletion processes within stars. △ Less

Submitted 4 July, 2021; v1 submitted 4 November, 2020; originally announced November 2020.

Comments: 12 pages, 7 figures, updated version resubmitted to MNRAS after responding to positive referee report. Also now has attached galah_spite_stars.csv file which is full version of Table 1

The GALAH survey: tracing the Galactic disk with Open Clusters

Authors: Lorenzo Spina, Yuan-Sen Ting, Gayandhi M. De Silva, Neige Frankel, Sanjib Sharma, Tristan Cantat-Gaudin, Meridith Joyce, Dennis Stello, Amanda I. Karakas, Martin B. Asplund, Thomas Nordlander, Luca Casagrande, Valentina D'Orazi, Andrew R. Casey, Peter Cottrell, Thor Tepper-García, Martina Baratella, Janez Kos, Klemen Čotar, Joss Bland-Hawthorn, Sven Buder, Ken C. Freeman, Michael R. Hayden, Geraint F. Lewis, Jane Lin , et al. (6 additional authors not shown)

Abstract: Open clusters are unique tracers of the history of our own Galaxy's disk. According to our membership analysis based on \textit{Gaia} astrometry, out of the 226 potential clusters falling in the footprint of GALAH or APOGEE, we find that 205 have secure members that were observed by at least one of the survey. Furthermore, members of 134 clusters have high-quality spectroscopic data that we use to… ▽ More Open clusters are unique tracers of the history of our own Galaxy's disk. According to our membership analysis based on \textit{Gaia} astrometry, out of the 226 potential clusters falling in the footprint of GALAH or APOGEE, we find that 205 have secure members that were observed by at least one of the survey. Furthermore, members of 134 clusters have high-quality spectroscopic data that we use to determine their chemical composition. We leverage this information to study the chemical distribution throughout the Galactic disk of 21 elements, from C to Eu. The radial metallicity gradient obtained from our analysis is $-$0.076$\pm$0.009 dex kpc$^{-1}$, which is in agreement with previous works based on smaller samples. Furthermore, the gradient in the [Fe/H] - guiding radius (r$_{\rm guid}$) plane is $-$0.073$\pm$0.008 dex kpc$^{-1}$. We show consistently that open clusters trace the distribution of chemical elements throughout the Galactic disk differently than field stars. In particular, at given radius, open clusters show an age-metallicity relation that has less scatter than field stars. As such scatter is often interpreted as an effect of radial migration, we suggest that these differences are due to the physical selection effect imposed by our Galaxy: clusters that would have migrated significantly also had higher chances to get destroyed. Finally, our results reveal trends in the [X/Fe]$-$r$_{\rm guid}$$-$age space, which are important to understand production rates of different elements as a function of space and time. △ Less

Submitted 16 February, 2021; v1 submitted 4 November, 2020; originally announced November 2020.

Comments: 18 pages, 10 figures, accepted for publication MNRAS

MSC Class: 85-02

The GALAH+ Survey: Third Data Release

Authors: Sven Buder, Sanjib Sharma, Janez Kos, Anish M. Amarsi, Thomas Nordlander, Karin Lind, Sarah L. Martell, Martin Asplund, Joss Bland-Hawthorn, Andrew R. Casey, Gayandhi M. De Silva, Valentina D'Orazi, Ken C. Freeman, Michael R. Hayden, Geraint F. Lewis, Jane Lin, Katharine J. Schlesinger, Jeffrey D. Simpson, Dennis Stello, Daniel B. Zucker, Tomaz Zwitter, Kevin L. Beeson, Tobias Buck, Luca Casagrande, Jake T. Clark , et al. (22 additional authors not shown)

Abstract: The ensemble of chemical element abundance measurements for stars, along with precision distances and orbit properties, provides high-dimensional data to study the evolution of the Milky Way. With this third data release of the Galactic Archaeology with HERMES (GALAH) survey, we publish 678 423 spectra for 588 571 mostly nearby stars (81.2% of stars are within <2 kpc), observed with the HERMES spe… ▽ More The ensemble of chemical element abundance measurements for stars, along with precision distances and orbit properties, provides high-dimensional data to study the evolution of the Milky Way. With this third data release of the Galactic Archaeology with HERMES (GALAH) survey, we publish 678 423 spectra for 588 571 mostly nearby stars (81.2% of stars are within <2 kpc), observed with the HERMES spectrograph at the Anglo-Australian Telescope. This release (hereafter GALAH+ DR3) includes all observations from GALAH Phase 1 (bright, main, and faint survey, 70%), K2-HERMES (17%), TESS-HERMES (5%), and a subset of ancillary observations (8%) including the bulge and >75 stellar clusters. We derive stellar parameters $T_\text{eff}$, $\log g$, [Fe/H], $v_\text{mic}$, $v_\text{broad}$ & $v_\text{rad}$ using our modified version of the spectrum synthesis code Spectroscopy Made Easy (SME) and 1D MARCS model atmospheres. We break spectroscopic degeneracies in our spectrum analysis with astrometry from $Gaia$ DR2 and photometry from 2MASS. We report abundance ratios [X/Fe] for 30 different elements (11 of which are based on non-LTE computations) covering five nucleosynthetic pathways. We describe validations for accuracy and precision, flagging of peculiar stars/measurements and recommendations for using our results. Our catalogue comprises 65% dwarfs, 34% giants, and 1% other/unclassified stars. Based on unflagged chemical composition and age, we find 62% young low-$α$, 9% young high-$α$, 27% old high-$α$, and 2% stars with $\mathrm{[Fe/H]} \leq -1$. Based on kinematics, 4% are halo stars. Several Value-Added-Catalogues, including stellar ages and dynamics, updated after $Gaia$ eDR3, accompany this release and allow chrono-chemodynamic analyses, as we showcase. △ Less

Submitted 28 April, 2021; v1 submitted 4 November, 2020; originally announced November 2020.

Comments: 48 pages, 27 figures, 5 pages appendix, accepted for publication in MNRAS, comments welcome, GALAH+ DR3 can be found at https://docs.datacentral.org.au/galah/dr3/overview/, Catalogs can be downloaded from https://cloud.datacentral.org.au/teamdata/GALAH/public/GALAH_DR3/

Galactic seismology: the evolving "phase spiral" after the Sagittarius dwarf impact

Authors: Joss Bland-Hawthorn, Thor Tepper-Garcia

Abstract: In 2018, the ESA \Gaia\ satellite discovered a remarkable spiral pattern ("phase spiral") in the $z-V_z$ phase plane throughout the solar neighbourhood, where $z$ and $V_z$ are the displacement and velocity of a star perpendicular to the Galactic disc. In response to Binney \& Schönrich's analytic model of a disc-crossing satellite to explain the \Gaia\ data, we carry out a high-resolution, N-body… ▽ More In 2018, the ESA \Gaia\ satellite discovered a remarkable spiral pattern ("phase spiral") in the $z-V_z$ phase plane throughout the solar neighbourhood, where $z$ and $V_z$ are the displacement and velocity of a star perpendicular to the Galactic disc. In response to Binney \& Schönrich's analytic model of a disc-crossing satellite to explain the \Gaia\ data, we carry out a high-resolution, N-body simulation (N$\:\approx 10^8$ particles) of an impulsive mass ($2\times 10^{10}$ \Msun) that interacts with a cold stellar disc at a single transit point. The disc response is complex since the impulse triggers a superposition of two distinct bisymmetric ($m=2$) modes $-$ a density wave and a corrugated bending wave $-$ that wrap up at different rates. Stars in the {\it faster} density wave wrap up with time $T$ according to $φ_D(R,T)=(Ω_D(R) + Ω_{\rm o})\:T$ where $φ_D$ describes the spiral pattern and $Ω_D =Ω(R) -κ(R)/2$, where $κ$ is the epicyclic frequency. While the pattern speed $Ω_{\rm o}$ is small, it is non-zero. The {\it slower} bending wave wraps up according to $Ω_B\approxΩ_D/2$ producing a corrugated wave. The bunching effect of the density wave triggers the phase spiral as it rolls up and down on the bending wave ("rollercoaster" model). The phase spiral emerges slowly about $ΔT \approx 400$ Myr after impact. It appears to be a long-lived, disc-wide phenomenon that continues to evolve over most of the 2~Gyr simulation. Thus, given Sagittarius' (Sgr) low total mass today ($M_{\rm tot}\sim 3\times 10^8$ \Msun\ within 10 kpc diameter), we believe the phase spiral was excited by the disc-crossing dwarf some $1-2$ Gyr {\it before} the recent transit. For this to be true, Sgr must be losing mass at 0.5-1 dex per orbit loop. △ Less

Submitted 8 March, 2021; v1 submitted 4 September, 2020; originally announced September 2020.

Comments: 21 pages, 16 figures, MNRAS accepted (Feb 2021). We have left modifications in bold print. A brief history of Galactic seismology is given in the Appendix

The M31/M33 tidal interaction: A hydrodynamic simulation of the extended gas distribution

Authors: Thor Tepper-García, Joss Bland-Hawthorn, Di Li

Abstract: We revisit the orbital history of the Triangulum galaxy (M33) around the Andromeda galaxy (M31) in view of the recent Gaia Data Release 2 proper motion measurements for both Local Group galaxies. Earlier studies consider highly idealised dynamical friction, but neglect the effects of dynamical mass loss. We show the latter process to be important using mutually consistent orbit integration and N-b… ▽ More We revisit the orbital history of the Triangulum galaxy (M33) around the Andromeda galaxy (M31) in view of the recent Gaia Data Release 2 proper motion measurements for both Local Group galaxies. Earlier studies consider highly idealised dynamical friction, but neglect the effects of dynamical mass loss. We show the latter process to be important using mutually consistent orbit integration and N-body simulations. Following this approach we find an orbital solution that brings these galaxies to within ~50 kpc of each other in the past, ~6.5 Gyr ago. We explore the implications of their interaction using an N-body/hydrodynamical simulation with a focus on the origin of two prominent features: 1) M31's Giant Stellar Stream; and 2) the M31-M33 HI filament. We find that the tidal interaction does not produce a structure reminiscent of the stellar stream that survives up to the present day. In contrast, the M31-M33 HI filament is likely a fossil structure dating back to the time of the ancient encounter between these galaxies. Similarly, the observed outer disc warp in M33 may well be a relic of this past event. Our model suggests the presence of a tidally induced gas envelope around these galaxies, and the existence of a diffuse gas stream, the 'Triangulum stream', stretching for tens of kpc from M33 away from M31. We anticipate upcoming observations with the recently commissioned, Five-hundred-meter Aperture Spherical radio Telescope (FAST) that will target the putative stream in its first years of operation. △ Less

Submitted 3 February, 2020; v1 submitted 26 November, 2019; originally announced November 2019.

Comments: MNRAS, in press. Radical change of title requested by referee; results and conclusions unchanged. Matches accepted version

The Magellanic Stream as a Probe of Astrophysics

Authors: Andrew J. Fox, Kathleen A. Barger, Joss Bland-Hawthorn, Dana Casetti-Dinescu, Elena D'Onghia, Felix J. Lockman, Naomi McClure-Griffiths, David Nidever, Mary Putman, Philipp Richter, Snezana Stanimirovic, Thorsten Tepper-Garcia

Abstract: Extending for over 200 degrees across the sky, the Magellanic Stream together with its Leading Arm is the most spectacular example of a gaseous stream in the local Universe. The Stream is an interwoven tail of filaments trailing the Magellanic Clouds as they orbit the Milky Way. Thought to be created by tidal forces, ram pressure, and halo interactions, the Stream is a benchmark for dynamical mode… ▽ More Extending for over 200 degrees across the sky, the Magellanic Stream together with its Leading Arm is the most spectacular example of a gaseous stream in the local Universe. The Stream is an interwoven tail of filaments trailing the Magellanic Clouds as they orbit the Milky Way. Thought to be created by tidal forces, ram pressure, and halo interactions, the Stream is a benchmark for dynamical models of the Magellanic System, a case study for gas accretion and dwarf-galaxy accretion onto galaxies, a probe of the outer halo, and the bearer of more gas mass than all other Galactic high velocity clouds combined. If it survives to reach the Galactic disk, it may maintain or even elevate the Galactic star-formation rate. In this white paper, we emphasize the Stream's importance for many areas of Galactic astronomy, summarize key unanswered questions, and identify future observations and simulations needed to resolve them. We stress the importance of ultraviolet, optical, and radio spectroscopy, and the need for computational models that capture full particle and radiation treatments within an MHD environment. △ Less

Submitted 12 March, 2019; originally announced March 2019.

Comments: Science White Paper submitted to Astro2020 Decadal Survey. 7 pages, 2 figures

The GALAH survey and Gaia DR2: Linking ridges, arches and vertical waves in the kinematics of the Milky Way

Authors: Shourya Khanna, Sanjib Sharma, Thor Tepper-Garcia, Joss Bland-Hawthorn, Michael Hayden, Martin Asplund, Sven Buder, Boquan Chen, Gayandhi M. De Silva, Ken C. Freeman, Janez Kos, Jane Lin, Sarah L. Martell, Jeffrey D. Simpson, Dennis Stello, Yuan-Sen Ting, Daniel B. Zucker, Tomaz Zwitter

Abstract: Gaia DR2 has revealed new small-scale and large-scale patterns in the phase-space distribution of stars in the Milky Way. In cylindrical Galactic coordinates $(R,φ,z)$, ridge-like structures can be seen in the \vphiR{} plane and asymmetric arch-like structures in the \vphivR{} plane. We show that the ridges are also clearly present when the third dimension of the \vphiR{} plane is represented by… ▽ More Gaia DR2 has revealed new small-scale and large-scale patterns in the phase-space distribution of stars in the Milky Way. In cylindrical Galactic coordinates $(R,φ,z)$, ridge-like structures can be seen in the \vphiR{} plane and asymmetric arch-like structures in the \vphivR{} plane. We show that the ridges are also clearly present when the third dimension of the \vphiR{} plane is represented by $\langle z \rangle$, $\langle V_z \rangle$, $\langle V_R \rangle$, $\langle$[Fe/H]$\rangle$ and $\langle[α/{\rm Fe}]\rangle$. The maps suggest that stars along the ridges lie preferentially close to the Galactic midplane ($|z|<0.2$ kpc), and have metallicity and $α$ elemental abundance similar to that of the Sun. We show that phase mixing of disrupting spiral arms can generate both the ridges and the arches. It also generates discrete groupings in orbital energy $-$ the ridges and arches are simply surfaces of constant energy. We identify 8 distinct ridges in the \gaia{} data: six of them have constant energy while two have constant angular momentum. Given that the signature is strongest for stars close to the plane, the presence of ridges in $\langle z \rangle$ and $\langle V_z \rangle$ suggests a coupling between planar and vertical directions. We demonstrate, using N-body simulations that such coupling can be generated both in isolated discs and in discs perturbed by an orbiting satellite like the Sagittarius dwarf galaxy. △ Less

Submitted 1 September, 2019; v1 submitted 26 February, 2019; originally announced February 2019.

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

The Magellanic System: the puzzle of the leading gas stream

Authors: Thor Tepper-García, Joss Bland-Hawthorn, Marcel S. Pawlowski, Tobias K. Fritz

Abstract: The Magellanic Clouds (MCs) are the most massive gas-bearing systems falling into the Galaxy at the present epoch. They show clear signs of interaction, manifested in particular by the Magellanic Stream, a spectacular gaseous wake that trails from the MCs extending more than 150 degree across the sky. Ahead of the MCs is the "Leading Arm" usually interpreted as the tidal counterpart of the Magella… ▽ More The Magellanic Clouds (MCs) are the most massive gas-bearing systems falling into the Galaxy at the present epoch. They show clear signs of interaction, manifested in particular by the Magellanic Stream, a spectacular gaseous wake that trails from the MCs extending more than 150 degree across the sky. Ahead of the MCs is the "Leading Arm" usually interpreted as the tidal counterpart of the Magellanic Stream, an assumption we now call into question. We revisit the formation of these gaseous features in a first-infall scenario, including for the first time a Galactic model with a weakly magnetised, spinning hot corona. In agreement with previous studies, we broadly recover the location and the extension of the Stream on the sky. In contrast, we find that the formation of the Leading Arm -- that is otherwise present in models without a corona -- is inhibited by the hydrodynamic interaction with the hot component. These results hold with or without coronal rotation or a weak, ambient magnetic field. Since the existence of the hot corona is well established, we are led to two possible interpretations: (i) the Leading Arm survives because the coronal density beyond 20 kpc is a factor of 10 lower than required by conventional spheroidal coronal x-ray models, consistent with recent claims of rapid coronal rotation; or (ii) the `Leading Arm' is cool gas trailing from a frontrunner, a satellite moving ahead of the MCs, consistent with its higher metallicity compared to the trailing stream. Both scenarios raise issues that we discuss. △ Less

Submitted 13 June, 2019; v1 submitted 17 January, 2019; originally announced January 2019.

Comments: MNRAS accepted. Matches accepted version

The GALAH survey and Gaia DR2: dissecting the stellar disc's phase space by age, action, chemistry and location

Authors: Joss Bland-Hawthorn, Sanjib Sharma, Thor Tepper-Garcia, James Binney, Ken Freeman, Michael Hayden, Janez Kos, Gayandhi De Silva, Simon Ellis, Shourya Khanna, Geraint Lewis, Martin Asplund, Sven Buder, Andrew Casey, Valentina D'Orazi, Ly Duong, Jane Lin, Karin Lind, Sarah Martell, Melissa Ness, Jeff Simpson, Daniel Zucker, Tomaz Zwitter, Prajwal Kafle, Alice Quillen , et al. (2 additional authors not shown)

Abstract: We use the second data releases of the ESA Gaia astrometric survey and the high-resolution GALAH spectroscopic survey to analyse the structure of our Galaxy's disc components. With GALAH, we separate the alpha-rich and alpha-poor discs (with respect to Fe), which are superposed in both position and velocity space, and examine their distributions in action space. We study the distribution of stars… ▽ More We use the second data releases of the ESA Gaia astrometric survey and the high-resolution GALAH spectroscopic survey to analyse the structure of our Galaxy's disc components. With GALAH, we separate the alpha-rich and alpha-poor discs (with respect to Fe), which are superposed in both position and velocity space, and examine their distributions in action space. We study the distribution of stars in the zV_z phase plane, for both V_phi and V_R, and recover the remarkable "phase spiral" discovered by Gaia. We identify the anticipated quadrupole signature in zV_z of a tilted velocity ellipsoid for stars above and below the Galactic plane. By connecting our work with earlier studies, we show that the phase spiral is likely to extend well beyond the narrow solar neighbourhood cylinder in which it was found. The phase spiral is a signature of corrugated waves that propagate through the disc, and the associated non-equilibrium phase mixing. The radially asymmetric distribution of stars involved in the phase spiral reveals that the corrugation, which is mostly confined to the alpha-poor disc, grows in z-amplitude with increasing radius. We present new simulations of tidal disturbance of the Galactic disc by the Sagittarius (Sgr) dwarf. The effect on the zV_z phase plane lasts >2 Gyr but a subsequent disc crossing wipes out the coherent structure. We find that the phase spiral was excited < 0.5 Gyr ago by an object like Sgr with total mass 3 x 10^10 Msun (stripped down from 5 x 10^10 Msun when it first entered the halo) passing through the plane. △ Less

Submitted 17 January, 2019; v1 submitted 7 September, 2018; originally announced September 2018.

Comments: 27 pages, 24 figures, MNRAS (accepted, Jan. 2019)

Magnetic fields in the Galactic halo restrict fountain-driven recycling and accretion

Authors: Asger Grønnow, Thor Tepper-García, Joss Bland-Hawthorn

Abstract: The Galactic halo contains a complex ecosystem of multiphase intermediate-velocity and high-velocity gas clouds whose origin has defied clear explanation. They are generally believed to be involved in a Galaxy-wide recycling process, either through an accretion flow or a large-scale fountain flow, or both. We examine the evolution of these clouds in light of recent claims that they may trigger con… ▽ More The Galactic halo contains a complex ecosystem of multiphase intermediate-velocity and high-velocity gas clouds whose origin has defied clear explanation. They are generally believed to be involved in a Galaxy-wide recycling process, either through an accretion flow or a large-scale fountain flow, or both. We examine the evolution of these clouds in light of recent claims that they may trigger condensation of gas from the Galactic corona as they move through it. We measure condensation along a cloud's wake, with and without the presence of an ambient magnetic field, using two- (2D) and three-dimensional (3D), high-resolution simulations. We find that 3D simulations are essential to correctly capture the condensation in all cases. Magnetic fields significantly inhibit condensation in the wake of clouds at $t \gtrsim 25$ Myr, preventing the sharp upturn in cold gas mass seen in previous non-magnetic studies. The magnetic field suppresses the Kelvin-Helmholtz instability responsible for the ablation and consequent mixing of a cloud with halo gas which drives the condensation. This effect is universal across different cloud properties (density, metallicity, velocity) and magnetic field properties (strength and orientation). Simple convergence tests demonstrate that resolving the gas on progressively smaller scales leads to even less condensation. While condensation still occurs in all cases, our results show that an ambient magnetic field drastically lowers the efficiency of fountain-driven accretion and likely also accretion from condensation around high-velocity clouds. These lower specific accretion rates are in better agreement with observational constraints compared to 3D, non-magnetic simulations. △ Less

Submitted 31 August, 2018; v1 submitted 10 May, 2018; originally announced May 2018.

Comments: 22 pages, 12 figures. Accepted for publication in ApJ

The Sagittarius dwarf galaxy: Where did all the gas go?

Authors: Thor Tepper-García, Joss Bland-Hawthorn

Abstract: The remarkable 1994 discovery of the Sagittarius dwarf galaxy (Sgr) revealed that, together with the Magellanic Clouds (MCs), there are at least three major dwarf galaxies, each with a total mass of order 10^10 - 10^11 solar masses, falling onto the Galaxy in the present epoch. Beyond a Galactic radius of 300 kpc, dwarfs tend to retain their gas. At roughly 50 kpc, the MCs have experienced substan… ▽ More The remarkable 1994 discovery of the Sagittarius dwarf galaxy (Sgr) revealed that, together with the Magellanic Clouds (MCs), there are at least three major dwarf galaxies, each with a total mass of order 10^10 - 10^11 solar masses, falling onto the Galaxy in the present epoch. Beyond a Galactic radius of 300 kpc, dwarfs tend to retain their gas. At roughly 50 kpc, the MCs have experienced substantial gas stripping as evidenced by the Magellanic Stream which extends from them. Since Sgr experienced star formation long after it fell into the Galaxy, it is interesting to explore just how and when this dwarf lost its gas. To date, there has been no definitive detection of an associated gas component. We revisit recent simulations of the stellar and dark matter components of Sgr but, for the first time, include gas that is initially bound to the infalling galaxy. We find that the gas stripping was 30 - 50 % complete at its first disc crossing ~2.7 Gyr ago, then entirely stripped at its last disc crossing ~1 Gyr ago. Our timeline is consistent with the last substantial burst of star formation in Sgr which occurred about the time of the last disc crossing. We discuss the consequences of gas stripping and conclude that the vast majority of the stripped gas was fully settled onto the Galaxy by ~300 Myr ago. It is highly unlikely that any of the high- or intermediate-velocity clouds have a direct association with the Sgr dwarf. △ Less

Submitted 22 May, 2018; v1 submitted 26 February, 2018; originally announced February 2018.

Comments: MNRAS accepted, virtually without changes after referee report

Revealing the Ionization Properties of the Magellanic Stream using Optical Emission

Authors: K. A. Barger, G. J. Madsen, A. J. Fox, B. P. Wakker, J. Bland-Hawthorn, D. Nidever, L. M. Haffner, J. Antwi-Danso, M. Hernandez, N. Lehner, A. S. Hill, A. Curzons, T. Tepper-Garcia

Abstract: The Magellanic Stream, a gaseous tail that trails behind the Magellanic Clouds, could replenish the Milky Way with a tremendous amount of gas if it reaches the Galactic disk before it evaporates into the halo. To determine how the Magellanic Stream's properties change along its length, we have conducted an observational study of the H-alpha emission, along with other optical warm ionized gas trace… ▽ More The Magellanic Stream, a gaseous tail that trails behind the Magellanic Clouds, could replenish the Milky Way with a tremendous amount of gas if it reaches the Galactic disk before it evaporates into the halo. To determine how the Magellanic Stream's properties change along its length, we have conducted an observational study of the H-alpha emission, along with other optical warm ionized gas tracers, toward 39 sight lines. Using the Wisconsin H-alpha Mapper telescope, we detect H-alpha emission brighter than 30 - 50 mR in 26 of our 39 sight lines. This H-alpha emission extends more than 2-degree away from the HI emission. By comparing H-alpha and [OI] intensities, we find that regions with log NHI = 19.5 - 20.0 are 16 - 67% ionized. Most of the H-alpha intensities along the Magellanic Stream are much higher than expected if the primary ionization source is photoionization from Magellanic Clouds, the Milky Way, and the extragalactic background. We find that the additional contribution from self ionization through a "shock cascade" that results as the Stream plows through the halo might be sufficient to reproduce the underlying level of H-alpha emission along the Stream. In the sparsely sampled region below the South Galactic Pole, there exists a subset of sight lines with uncharacteristically bright emission, which suggest that gas is being ionized further by an additional source that could be a linked to energetic processes associated with the Galactic center. △ Less

Submitted 12 November, 2017; originally announced November 2017.

Comments: Accepted for publication in the Astrophysical Journal on Nov 8, 2017

The Smith Cloud: surviving a high speed transit of the Galactic disc

Authors: Thor Tepper-Garcia, Joss Bland-Hawthorn

Abstract: The origin and survival of the Smith high-velocity HI cloud has so far defied explanation. This object has several remarkable properties: (i) its prograde orbit is ~100 km/s faster than the underlying Galactic rotation; (ii) its total gas mass ($\gtrsim 4 \times 10^6 ~{\rm M}_{\odot}$) exceeds the mass of all other high-velocity clouds (HVC) outside of the Magellanic Stream; (iii) its head-tail mo… ▽ More The origin and survival of the Smith high-velocity HI cloud has so far defied explanation. This object has several remarkable properties: (i) its prograde orbit is ~100 km/s faster than the underlying Galactic rotation; (ii) its total gas mass ($\gtrsim 4 \times 10^6 ~{\rm M}_{\odot}$) exceeds the mass of all other high-velocity clouds (HVC) outside of the Magellanic Stream; (iii) its head-tail morphology extends to the Galactic HI disc, indicating some sort of interaction. The Smith Cloud's kinetic energy rules out models based on ejection from the disc. We construct a dynamically self-consistent, multi-phase model of the Galaxy with a view to exploring whether the Smith Cloud can be understood in terms of an infalling, compact HVC that has transited the Galactic disc. We show that while a dark-matter (DM) free HVC of sufficient mass and density can reach the disc, it does not survive the transit. The most important ingredient to survival during a transit is a confining DM subhalo around the cloud; radiative gas cooling and high spatial resolution ($\lesssim$ 10 pc) are also essential. In our model, the cloud develops a head-tail morphology within ~10 Myr before and after its first disc crossing; after the event, the tail is left behind and accretes onto the disc within ~400 Myr. In our interpretation, the Smith Cloud corresponds to a gas 'streamer' that detaches, falls back and fades after the DM subhalo, distorted by the disc passage, has moved on. We conclude that subhalos with ${\rm M}_{\rm DM} \lesssim 10^9 ~{\rm M}_{\odot}$ have accreted $\sim10^9 ~{\rm M}_{\odot}$ of gas into the Galaxy over cosmic time - a small fraction of the total baryon budget. △ Less

Submitted 12 October, 2017; originally announced October 2017.

Comments: MNRAS accepted

Magnetized High Velocity Clouds in the Galactic Halo: A New Distance Constraint

Authors: Asger Grønnow, Thor Tepper-García, Joss Bland-Hawthorn, Naomi McClure-Griffiths

Abstract: High velocity gas that does not conform to Galactic rotation is observed throughout the Galaxy's halo. One component of this gas, HI high velocity clouds (HVCs), have attracted attention since their discovery in the 1960s and remain controversial in terms of their origins, largely due to the lack of reliable distance estimates. The recent discovery of enhanced magnetic fields towards HVCs has enco… ▽ More High velocity gas that does not conform to Galactic rotation is observed throughout the Galaxy's halo. One component of this gas, HI high velocity clouds (HVCs), have attracted attention since their discovery in the 1960s and remain controversial in terms of their origins, largely due to the lack of reliable distance estimates. The recent discovery of enhanced magnetic fields towards HVCs has encouraged us to explore their connection to cloud evolution, kinematics, and survival as they fall through the magnetized Galactic halo. For a reasonable model for the halo magnetic field, most infalling clouds see transverse rather than radial field lines. We find that significant compression (and thereby amplification) of the ambient magnetic field occurs in front of the cloud and in the tail of material stripped from the cloud. The compressed transverse field attenuates hydrodynamical instabilities. This delays cloud destruction, though not indefinitely. The observed B-field compression is related to the cloud's distance from the Galactic plane. As a result, the observed rotation measure provides useful distance information on a cloud's location. △ Less

Submitted 23 August, 2017; v1 submitted 10 July, 2017; originally announced July 2017.

Comments: 18 pages, 14 figures. Published in ApJ

Journal ref: ApJ.845.69(2017)

An HST/COS legacy survey of intervening SiIII absorption in the extended gaseous halos of low-redshift galaxies

Authors: P. Richter, B. P. Wakker, C. Fechner, P. Herenz, T. Tepper-Garcia, A. J. Fox

Abstract: Doubly ionized silicon (SiIII) is a powerful tracer of diffuse ionized gas inside and outside of galaxies. It can be observed in the local Universe in ultraviolet (UV) absorption against bright extragalactic background sources. We here present an extensive study of intervening SiIII-selected absorbers and their relation to the circumgalactic medium (CGM) of galaxies at low redshift (z<=0.1), based… ▽ More Doubly ionized silicon (SiIII) is a powerful tracer of diffuse ionized gas inside and outside of galaxies. It can be observed in the local Universe in ultraviolet (UV) absorption against bright extragalactic background sources. We here present an extensive study of intervening SiIII-selected absorbers and their relation to the circumgalactic medium (CGM) of galaxies at low redshift (z<=0.1), based on the analysis of UV absorption spectra along 303 extragalactic lines of sight obtained with the Cosmic Origins Spectrograph (COS) on board the Hubble Space Telescope (HST). Along a total redshift path of Dz=24 we identify 69 intervening SiIII systems that all show associated absorption from other low and high ions. We derive a bias-corrected number density of dN/dz(SiIII)=2.5 for absorbers with column densities log N(SiIII)>12.2. We develop a geometrical model for the absorption-cross section of the CGM around the local galaxy population and find excellent agreement between the model predictions and the observations. We further compare redshifts and positions of the absorbers with that of ~64,000 galaxies using archival galaxy-survey data. For the majority of the absorbers we identify possible host galaxies within 300 km/s of the absorbers and derive impact parameters rho<200 kpc, demonstrating that the spatial distributions of SiIII absorbers and galaxies are highly correlated. Our study indicates that the majority of SiIII-selected absorbers in our sample trace the CGM of nearby galaxies within their virial radii at a typical covering fraction of ~70 per cent. From a detailed ionization model we estimate that diffuse gas in the CGM around galaxies, as traced by SiIII, contains substantially more baryonic mass than their neutral interstellar medium. △ Less

Submitted 17 March, 2016; v1 submitted 22 July, 2015; originally announced July 2015.

Comments: 32 pages, 17 figures; final version accepted for publication in A&A

Journal ref: A&A 590, A68 (2016)

The Magellanic Stream: break up and accretion onto the hot Galactic corona

Authors: Thor Tepper-Garcia, Joss Bland-Hawthorn, Ralph S. Sutherland

Abstract: The Magellanic HI Stream (~2x10^9 Msun [d/55 kpc]^2) encircling the Galaxy at a distance 'd' is arguably the most important tracer of what happens to gas accreting onto a disk galaxy. Recent observations reveal that the Stream's mass is in fact dominated (3:1) by its ionised component. Here we revisit the origin of the mysterious H-alpha recombination emission observed along much of its length tha… ▽ More The Magellanic HI Stream (~2x10^9 Msun [d/55 kpc]^2) encircling the Galaxy at a distance 'd' is arguably the most important tracer of what happens to gas accreting onto a disk galaxy. Recent observations reveal that the Stream's mass is in fact dominated (3:1) by its ionised component. Here we revisit the origin of the mysterious H-alpha recombination emission observed along much of its length that is overly bright (~150-200 milli-Rayleigh) for the known Galactic ultraviolet background (~20-40 mR / [d/55 kpc]^2). In an earlier model, we proposed that a slow shock cascade was operating along the Stream due to its interaction with the extended Galactic hot corona. We find that, for a smooth coronal density profile, this model can explain the bright H-alpha emission if the coronal density satisfies 2 < (n / 10^{-4} cm^{-3}) < 4 at d = 55 kpc. But in view of updated parameters for the Galactic halo and mounting evidence that most of the Stream must lie far beyond the Magellanic Clouds (d>55 kpc), we revisit the shock cascade model in detail. At lower densities, the HI gas is broken down by the shock cascade but mostly mixes with the hot corona without significant recombination. At higher densities, the hot coronal mass (including the other baryonic components) exceeds the baryon budget of the Galaxy. If the H-alpha emission arises from the shock cascade, the upper limit on the smooth coronal density constrains the Stream's mean distance to < 75 kpc. If, as some models indicate, the Stream is even further out, either the shock cascade is operating in a regime where the corona is substantially mass-loaded with recent gas debris, or an entirely different ionization mechanism is responsible. △ Less

Submitted 23 September, 2015; v1 submitted 7 May, 2015; originally announced May 2015.

Comments: Significant expansion of the parameter space explored in response to referee's comments. ApJ accepted

Absorption signatures of warm-hot gas at low redshift: NeVIII

Authors: Thor Tepper-Garcia, Philipp Richter, Joop Schaye

Abstract: At z<1 a large fraction of the baryons is thought to reside in diffuse gas that has been shock-heated to high temperatures (log T=5-6). Absorption by the 770.41, 780.32 Å doublet of NeVIII in quasar spectra represents a unique tool to study this elusive warm-hot phase. We have developed an analytic model for the properties of NeVIII absorbers that allows for an inhomogeneous metal distribution. Ou… ▽ More At z<1 a large fraction of the baryons is thought to reside in diffuse gas that has been shock-heated to high temperatures (log T=5-6). Absorption by the 770.41, 780.32 Å doublet of NeVIII in quasar spectra represents a unique tool to study this elusive warm-hot phase. We have developed an analytic model for the properties of NeVIII absorbers that allows for an inhomogeneous metal distribution. Our model agrees with the predictions of a simulation from the OWLS project indicating that the average line-of-sight metal-filling fraction within the absorbing gas is low (c_L~0.1). Most of the NeVIII in our model is produced in low-density, collisionally ionized gas (log n=-6 to -4 cm^-3, log T=5-6). Strong NeVIII absorbers (log N(NeVIII)>14), like those recently detected by HST/COS, are found to arise in higher density gas (log n>-4, log T=5.75). NeVIII cloudlets harbour only 1 per cent of the cosmic baryon budget. The baryon content of the surrounding gas (which has similar densities and temperatures as the NeVIII cloudlets) is a factor 1/c_L higher. We conclude that NeVIII absorbers are robust probes of shock-heated diffuse gas, but that spectra with signal-to-noise ratios S/N>100 would be required to detect the bulk of the baryons in warm-hot gas. △ Less

Submitted 7 October, 2013; v1 submitted 25 March, 2013; originally announced March 2013.

Comments: 19 pages, 10 figures; MNRAS in press; minor changes with respect to v1

Absorption signatures of warm-hot gas at low redshift: Broad HI Lyman-Alpha Absorbers

Authors: Thorsten Tepper-García, Philipp Richter, Joop Schaye, Craig M. Booth, Claudio Dalla Vecchia, Tom Theuns

Abstract: We investigate the physical state of HI absorbing gas at low redshift (z = 0.25) using a subset of cosmological, hydrodynamic simulations from the OWLS project, focusing in particular on broad (b_HI > 40 km/s) HI Lyman-Alpha absorbers (BLAs), which are believed to originate in shock-heated gas in the warm-hot intergalactic medium (WHIM). Our fiducial model, which includes radiative cooling by heav… ▽ More We investigate the physical state of HI absorbing gas at low redshift (z = 0.25) using a subset of cosmological, hydrodynamic simulations from the OWLS project, focusing in particular on broad (b_HI > 40 km/s) HI Lyman-Alpha absorbers (BLAs), which are believed to originate in shock-heated gas in the warm-hot intergalactic medium (WHIM). Our fiducial model, which includes radiative cooling by heavy elements and feedback by supernovae and active galactic nuclei, predicts that by z = 0.25 nearly 60 per cent of the gas mass ends up at densities and temperatures characteristic of the WHIM and we find that half of this fraction is due to outflows. The standard HI observables (distribution of HI column densities N_HI, distribution of Doppler parameters b_HI, b_HI - N_HI correlation) and the BLA line number density predicted by our simulations are in remarkably good agreement with observations. BLAs arise in gas that is hotter, more highly ionised and more enriched than the gas giving rise to typical Lyman-Alpha forest absorbers. The majority of the BLAs arise in warm-hot (log (T/K) ~ 5) gas at low (log Delta < 1.5) overdensities. On average, thermal broadening accounts for at least 60 per cent of the BLA line width, which in turn can be used as a rough indicator of the thermal state of the gas. Detectable BLAs account for only a small fraction of the true baryon content of the WHIM at low redshift. In order to detect the bulk of the mass in this gas phase, a sensitivity at least one order of magnitude better than achieved by current ultraviolet spectrographs is required. We argue that BLAs mostly trace gas that has been shock-heated and enriched by outflows and that they therefore provide an important window on a poorly understood feedback process. △ Less

Submitted 2 July, 2012; v1 submitted 26 January, 2012; originally announced January 2012.

Comments: 25 pages (20 pages main text), 19 figures; accepted for publication in MNRAS

Absorption signatures of warm-hot gas at low redshift: OVI

Authors: Thorsten Tepper-Garcia, Philipp Richter, Joop Schaye, C. M. Booth, Claudio Dalla Vecchia, Tom Theuns, Robert P. C. Wiersma

Abstract: We investigate the origin and physical properties of OVI absorbers at low redshift (z = 0.25) using a subset of cosmological, hydrodynamical simulations from the OverWhelmingly Large Simulations (OWLS) project. Intervening OVI absorbers are believed to trace shock-heated gas in the Warm-Hot Intergalactic Medium (WHIM) and may thus play a key role in the search for the missing baryons in the presen… ▽ More We investigate the origin and physical properties of OVI absorbers at low redshift (z = 0.25) using a subset of cosmological, hydrodynamical simulations from the OverWhelmingly Large Simulations (OWLS) project. Intervening OVI absorbers are believed to trace shock-heated gas in the Warm-Hot Intergalactic Medium (WHIM) and may thus play a key role in the search for the missing baryons in the present-day Universe. When compared to observations, the predicted distributions of the different OVI line parameters (column density, Doppler parameter, rest equivalent width) from our simulations exhibit a lack of strong OVI absorbers. This suggests that physical processes on sub-grid scales (e.g. turbulence) may strongly influence the observed properties of OVI systems. We find that the intervening OVI absorption arises mainly in highly metal-enriched (0.1 << Z/Z_sun < 1) gas at typical overdensities of 1 << rho/<rho> < 100. One third of the OVI absorbers in our simulation are found to trace gas at temperatures T < 10^5 K, while the rest arises in gas at higher temperatures around T =10^5.3 K. The OVI resides in a similar region of (rho,T)-space as much of the shock-heated baryonic matter, but the vast majority of this gas has a lower metal content and does not give rise to detectable OVI absorption As a consequence of the patchy metal distribution, OVI absorbers in our simulations trace only a very small fraction of the cosmic baryons (<2 percent) and the cosmic metals. Instead, these systems presumably trace previously shock-heated, metal-rich material from galactic winds that is now cooling. The common approach of comparing OVI and HI column densities to estimate the physical conditions in intervening absorbers from QSO observations may be misleading, as most of the HI (and most of the gas mass) is not physically connected with the high-metallicity patches that give rise to the OVI absorption. △ Less

Submitted 2 December, 2010; v1 submitted 16 July, 2010; originally announced July 2010.

Comments: Accepted for publication in MNRAS

Journal ref: 2011MNRAS.413..190T

Stochastic Absorption of the Light of Background Sources due to Intergalactic Neutral Hydrogen I. Testing different line-number evolution models via the cosmic flux decrement

Authors: Thorsten Tepper-Garcia, Uta Fritze

Abstract: [Abridged] We test the accuracy of different models of the attenuation of light due to resonant scattering by intergalactic neutral hydrogen by comparing their predictions of the evolution of the mean cosmic flux decrement, D_A, to measurements of this quantity based on observations. To this end, we use data available in the literature and our own measurements of the cosmic flux decrement for 25… ▽ More [Abridged] We test the accuracy of different models of the attenuation of light due to resonant scattering by intergalactic neutral hydrogen by comparing their predictions of the evolution of the mean cosmic flux decrement, D_A, to measurements of this quantity based on observations. To this end, we use data available in the literature and our own measurements of the cosmic flux decrement for 25 quasars in the redshift range 2.71 < z < 5.41 taken from the SDSS Data Release 5. In order to perform the measurements of D_A, we fit a power-law to the continuum redward of the Lya emission line, and extrapolate this fit to region blueward of it, where the flux is severely affected by absorption due to intervening HI absorbers. We compute, using numerical simulations, the redshift evolution of D_A accounting for the presence of Lya Forest absorbers and Lyman limit systems randomly distributed along the line-of-sight, and compute its intrinsic scatter at the 1-, 2-, and 3-sigma level due to fluctuations in the absorber properties (column density, Doppler parameter, redshift) along different lines-of-sight. The numerical simulations consist of Monte Carlo realizations of distributions of the absorber properties constrained from observations. The results from the models considered here confirm our theoretical expectation that the distribution of D_A at any given redshift be well described by a lognormal distribution function. This implies that the effective optical depth, usually defined as the negative logarithm of the average flux, 1 - D_A, is very accurately Gaussian distributed, in contrast to previous studies. This result is independent to the form of the input distribution functions, and rather insensitive to the presence of high-column density absorbers, such as the Lyman limit systems. △ Less

Submitted 9 November, 2007; v1 submitted 9 May, 2007; originally announced May 2007.

Comments: Accepted by MNRAS. Matches accepted version after changes following referee's report