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.

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.

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.

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.

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.

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.

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.

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.

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.

[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.