Search

Motivated by the need for realistic, dynamically self-consistent, evolving galaxy models that avoid the complexity of full, and zoom-in, cosmological simulations, we have developed NEXUS, an integral framework to create and evolve synthetic galaxies made of collisionless and gaseous components. NEXUS leverages the power of publicly available, tried-and-tested packages: the stellar-dynamics, action-based library AGAMA; and 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 galaxy formation physics, including gas cooling and heating, star formation, stellar feedback, and chemical enrichment. NEXUS' basic functionality consists in the generation of bespoke initial conditions (ICs) for a diversity of galaxy models, which are advanced in time to simulate the galaxy'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: hot halos and gas discs. NEXUS allows constructing equilibrium models with disc gas fractions $0~\leq~f_{\rm gas}~\leq~1$, appropriate to model both low- and high-redshift galaxies. Similarly, the framework is ideally suited to the study of galactic ecology, i.e. the dynamical interplay between stars and gas over billions of years. As a validation and illustration of our framework, we reproduce several isolated galaxy model setups reported in earlier studies, and present a new, `nested bar' galaxy simulation. Future upgrades of NEXUS will include magneto-hydrodynamics and highly energetic particle (`cosmic ray') heating., Comment: MNRAS accepted

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., Comment: 40 pages, 21 figures, ApJ, in press (April 2024) - simulation movies at http://www.physics.usyd.edu.au/turbo_disks

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., Comment: 41 pages, 34 figures, 3 tables. Accepted for Publication in the Astrophysical Journal

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_\phi,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., Comment: Accepted to the Astronomical Journal

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 $\gamma$-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 $\gamma$-rays of astrophysical origin, complicate the use of Sgr in indirect DM detection searches., Comment: 15 pages, 15 figures. Accepted for publication in MNRAS. Accepted manuscript uploaded

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{\deg}, 130{\deg}

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_\phi$ 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], [$\alpha$/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., Comment: 17 pages, 14 figures. Accepted publication in The Astronomical Journal

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., Comment: MNRAS, in press. Title slightly shortened. References added. Matches accepted version

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., Comment: 15 pages, 14 figures. Accepted for publication in MNRAS

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., Comment: 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 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., Comment: 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

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., Comment: 18 pages, 10 figures, accepted for publication MNRAS

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-$\alpha$, 9% young high-$\alpha$, 27% old high-$\alpha$, 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., Comment: 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/

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\"onrich'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 $\phi_D(R,T)=(\Omega_D(R) + \Omega_{\rm o})\:T$ where $\phi_D$ describes the spiral pattern and $\Omega_D =\Omega(R) -\kappa(R)/2$, where $\kappa$ is the epicyclic frequency. While the pattern speed $\Omega_{\rm o}$ is small, it is non-zero. The {\it slower} bending wave wraps up according to $\Omega_B\approx\Omega_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 $\Delta 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., Comment: 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

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., Comment: MNRAS, in press. Radical change of title requested by referee; results and conclusions unchanged. Matches accepted version

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,\phi,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[\alpha/{\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 $\alpha$ 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., Comment: 18 pages, 16 figures, accepted for publication in MNRAS

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., Comment: MNRAS accepted. Matches accepted version

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., Comment: 27 pages, 24 figures, MNRAS (accepted, Jan. 2019)

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., Comment: 22 pages, 12 figures. Accepted for publication in ApJ

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., Comment: MNRAS accepted, virtually without changes after referee report

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., Comment: MNRAS accepted

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., Comment: 18 pages, 14 figures. Published in ApJ

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., Comment: Significant expansion of the parameter space explored in response to referee's comments. ApJ accepted

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 {\AA} 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., Comment: 19 pages, 10 figures; MNRAS in press; minor changes with respect to v1