The goal of this study is to present the development of a machine learning based approach that utilizes phase space alone to separate the Gaia DR2 stars into two categories: those accreted onto the Milky Way from those that are in situ. Traditional selection methods that have been used to identify accreted stars typically rely on full 3D velocity, metallicity information, or both, which significantly reduces the number of classifiable stars. The approach advocated here is applicable to a much larger portion of Gaia DR2. A method known as "transfer learning" is shown to be effective through extensive testing on a set of mock Gaia catalogs that are based on the FIRE cosmological zoom-in hydrodynamic simulations of Milky Way-mass galaxies. The machine is first trained on simulated data using only 5D kinematics as inputs and is then further trained on a cross-matched Gaia/RAVE data set, which improves sensitivity to properties of the real Milky Way. The result is a catalog that identifies around 767,000 accreted stars within Gaia DR2. This catalog can yield empirical insights into the merger history of the Milky Way and could be used to infer properties of the dark matter distribution.

In Ostdiek et al. (2019), we developed a deep neural network classifier that only relies on phase-space information to obtain a catalog of accreted stars based on the second data release of Gaia (DR2). In this paper, we apply two clustering algorithms to identify velocity substructure within this catalog. We focus on the subset of stars with line-of-sight velocity measurements that fall in the range of Galactocentric radii $r \in [6.5, 9.5]$ kpc and vertical distances $|z| < 3$ kpc. Known structures such as Gaia Enceladus and the Helmi stream are identified. The largest previously-unknown structure, Nyx, first introduced in Necib et al. (2019a), is a vast stream consisting of at least 500 stars in the region of interest. This study displays the power of the machine learning approach by not only successfully identifying known features, but also discovering new kinematic structures that may shed light on the merger history of the Milky Way.

Massive dwarf galaxies that merge with the Milky Way on prograde orbits can be dragged into the disk plane before being completely disrupted. Such mergers can contribute to an accreted stellar disk and a dark matter disk. We present evidence for Nyx, a vast new stellar stream in the vicinity of the Sun, that may provide the first indication that such an event occurred in the Milky Way. We identify about 500 stars that have coherent radial and prograde motion in this stream using a catalog of accreted stars built by applying deep learning methods to the second Gaia data release. Nyx is concentrated within $\pm 2$ kpc of the Galactic midplane and spans the full radial range studied (6.5-9.5 kpc). The kinematics of Nyx stars are distinct from those of both the thin and thick disk. In particular, its rotational speed lags the disk by $\sim 80$ km/s and its stars follow more eccentric orbits. A small number of Nyx stars have chemical abundances or inferred ages; from these, we deduce that Nyx stars have a peak metallicity of [Fe/H] $\sim -0.5$ and ages $\sim $10-13 Gyr. Taken together with the kinematic observations, these results strongly favor the interpretation that Nyx is the remnant of a disrupted dwarf galaxy. To further justify this interpretation, we explicitly demonstrate that metal-rich, prograde streams like Nyx can be found in the disk plane of Milky Way-like galaxies using the FIRE hydrodynamic simulations. Future spectroscopic studies will be able to validate whether Nyx stars originate from a single progenitor.

We consider a twin WIMP scenario whose twin sector contains a full dark copy of the SM hadrons, where the lightest twin particles are twin pions. By analogy to the standard WIMP paradigm, the dark matter (DM) freezes out through twin electroweak interactions, and annihilates into a dark shower of light twin hadrons. These are either stable or decay predominantly to standard model (SM) photons. We show that this 'hadrosymmetric' scenario can be consistent with all applicable astrophysical, cosmological and collider constraints. In order to decay the twin hadrons before the big-bang nucleosynthesis epoch, an additional portal between the SM and twin sector is required. In most cases we find this additional mediator is within reach of either the LHC or future intensity frontier experiments. Furthermore, we conduct simulations of the dark shower and consequent photon spectra. We find that fits of these spectra to the claimed galactic center gamma-ray excess seen by Fermi-LAT non-trivially coincide with regions of parameter space that both successfully generate the observed DM abundance and exhibit minimal fine-tuning.

Magnetars are surrounded by diffuse plasma in magnetic field strengths well above the quantum electrodynamic critical value. We derive equations of "quantum force-free electrodynamics" for this plasma using effective field theory arguments. We argue that quantum effects do not modify the large scale structure of the magnetosphere, and in particular that the spin-down rate does not deviate significantly from the classical result. We provide definite evolution equations that can be used to explore potentially important small-scale corrections, such as shock formation, which has been proposed as a mechanism for both burst and quiescent emission from magnetars.