Search

Based on the rate of change of its orbital period, PSR J2043+1711 has a substantial peculiar acceleration of 3.5 $\pm$ 0.8 mm/s/yr, which deviates from the acceleration predicted by equilibrium Milky Way models at a $4\sigma$ level. The magnitude of the peculiar acceleration is too large to be explained by disequilibrium effects of the Milky Way interacting with orbiting dwarf galaxies ($\sim$1 mm/s/yr), and too small to be caused by period variations due to the pulsar being a redback. We identify and examine two plausible causes for the anomalous acceleration: a stellar flyby, and a long-period orbital companion. We identify a main-sequence star in \textit{Gaia} DR3 and Pan-STARRS DR2 with the correct mass, distance, and on-sky position to potentially explain the observed peculiar acceleration. However, the star and the pulsar system have substantially different proper motions, indicating that they are not gravitationally bound. However, it is possible that this is an unrelated star that just happens to be located near J2043+1711 along our line of sight (chance probability of 1.6\%). Therefore, we also constrain possible orbital parameters for a circumbinary companion in a hierarchical triple system with J2043+1711; the changes in the spindown rate of the pulsar are consistent with an outer object that has an orbital period of 80 kyr, a companion mass of 0.3 $M_\odot$ (indicative of a white dwarf or low-mass star), and a semi-major axis of 2000 AU. Continued timing and/or future faint optical observations of J2043+1711 may eventually allow us to differentiate between these scenarios.

Measurements of the accelerations of stars enabled by time-series extreme-precision spectroscopic observations, from pulsar timing, and from eclipsing binary stars in the Solar Neighborhood offer insights into the mass distribution of the Milky Way that do not rely on traditional equilibrium modeling. Given the measured accelerations, we can determine a total mass density, and from this, by accounting for the mass in stars, gas, and dust, we can infer the amount of dark matter. Leveraging the FIRE-2 simulations of Milky Way-mass galaxies, we compare vertical acceleration profiles between cold dark matter (CDM) and self-interacting dark matter (SIDM) with constant cross-section of 1 cm$^2$ g$^{-1}$ across three halos with diverse assembly histories. Notably, significant asymmetries in vertical acceleration profiles near the midplane at fixed radii are observed in both CDM and SIDM, particularly in halos recently affected by mergers with satellites of Sagittarius/SMC-like masses or greater. These asymmetries offer a unique window into exploring the merger history of a galaxy. We show that SIDM halos consistently exhibit higher local stellar and dark matter densities and steeper vertical acceleration gradients, up to 30% steeper near the Solar Neighborhood. SIDM halos also manifest a more oblate halo shape in the Solar Neighborhood. Furthermore, enhanced precision in acceleration measurements and larger datasets promise to provide better constraints on the local dark matter density, complementing our understanding from kinematic analysis of their distribution within galaxies., Comment: 14 pages, 5 figures, and 2 tables. Submitted to APJ

We measure the line-of-sight accelerations of 26 binary pulsars due to the Milky Way's gravitational potential, and produce a 3-dimensional map of the acceleration field of the Galaxy. Acceleration measurements directly give us the change in the line-of-sight velocity at present day, without requiring any assumptions inherent to kinematic modeling. We measure the Oort limit ($\rho_0=0.062\pm0.017$ \msun/pc$^3$) and the dark matter density in the midplane ($\rho_{0,\textrm{DM}}=-0.010\pm0.018$ \msun/pc$^3$); these values are similar to, but have smaller uncertainties than previous pulsar timing measurements of these quantities. Here, we provide for the first time, values for the Oort constants and the slope of the rotation curve from direct acceleration measurements. We find that $A=15.4\pm2.6$ km/s/kpc and $B=-13.1\pm2.6$ km/s/kpc (consistent with results from \textit{Gaia}), and the slope of the rotation curve near the Sun is $-2\pm5$ km/s/kpc. We show that the Galactic acceleration field is clearly asymmetric, but due to data limitations it is not yet clear which physical processes drive this asymmetry. We provide updated models of the Galactic potential that account for various sources of disequilibrium; these models are incompatible with commonly used kinematic potentials. This indicates that use of kinematically derived Galactic potentials in precision tests (e.g., in tests of general relativity with pulsar timing) may be subject to larger uncertainties than reported. The acceleration data indicates that the mass of the Galaxy within the Solar circle is $2.3 \times 10^{11}$ M$_\odot$, roughly twice as large as currently accepted models. Additionally, the residuals of the acceleration data compared to existing Galactic models have a dependence on radial position; this trend can be explained if the Sun has an additional acceleration away from the Galactic center.

The Milky Way's (MW) inner stellar halo contains an [Fe/H]-rich component with highly eccentric orbits, often referred to as the "last major merger." Hypotheses for the origin of this component include Gaia-Sausage/Enceladus (GSE), where the progenitor collided with the MW proto-disk 8-11 Gyr ago, and the Virgo Radial Merger (VRM), where the progenitor collided with the MW disk within the last 3 Gyr. These two scenarios make different predictions about observable structure in local phase space, because the morphology of debris depends on how long it has had to phase mix. The recently-identified phase-space folds in Gaia DR3 have positive caustic velocities, making them fundamentally different than the phase-mixed chevrons found in simulations at late times. Roughly 20\% of the stars in the prograde local stellar halo are associated with the observed caustics. Based on a simple phase-mixing model, the observed number of caustics are consistent with a merger that occurred 1--2 Gyr ago. We also compare the observed phase-space distribution to FIRE-2 Latte simulations of GSE-like mergers, using a quantitative measurement of phase mixing (2D causticality). The observed local phase-space distribution best matches the simulated data 1--2 Gyr after collision, and certainly not later than 3 Gyr. This is further evidence that the progenitor of the "last major merger" did not collide with the MW proto-disk at early times, as is thought for the GSE, but instead collided with the MW disk within the last few Gyr, consistent with the body of work surrounding the VRM.

We find that the chemical abundances and dynamics of APOGEE and GALAH stars in the local stellar halo are inconsistent with a scenario in which the inner halo is primarily composed of debris from a single, massive, ancient merger event, as has been proposed to explain the Gaia-Enceladus/Gaia Sausage (GSE) structure. The data contains trends of chemical composition with energy which are opposite to expectations for a single massive, ancient merger event, and multiple chemical evolution paths with distinct dynamics are present. We use a Bayesian Gaussian mixture model regression algorithm to characterize the local stellar halo, and find that the data is best fit by a model with four components. We interpret these components as the VRM, Cronus, Nereus, and Thamnos; however, Nereus and Thamnos likely represent more than one accretion event because the chemical abundance distributions of their member stars contain many peaks. Although the Cronus and Thamnos components have different dynamics, their chemical abundances suggest they may be related. We show that the distinct low and high alpha halo populations from Nissen & Schuster (2010) are explained by VRM and Cronus stars, as well as some in-situ stars. Because the local stellar halo contains multiple substructures, different popular methods of selecting GSE stars will actually select different mixtures of these substructures, which may change the apparent chemodynamic properties of the selected stars. We also find that the Splash stars in the solar region are shifted to higher v_phi and slightly lower [Fe/H] than previously reported.

We use halo dwarf stars with photometrically determined metallicities that are located within 2 kpc of the Sun to identify local halo substructure. The kinematic properties of these stars do not indicate a single, dominant radial merger event (RME). The retrograde Virgo Radial Merger (VRM) component has [Fe/H] = -1.7. A second, non-rotating RME component we name Nereus is identified with [Fe/H] = -2.1 and has similar energy as the VRM. A possible third RME we name Cronus is identified that is co-rotating with the disk, has lower energy than the VRM, and has [Fe/H] = -1.2. We identify the Nyx Stream in the data. In addition to these substructures, we observe metal-poor halo stars ([Fe/H] ~ -2.0 and $\sigma_v$ ~ 180 km s$^{-1}$) and a disk/Splash component with lower rotational velocity than the disk and lower metallicity than typically associated with the Splash. An additional excess of halo stars with low velocity and metallicity of [Fe/H] = -1.5 could be associated with the shell of a lower energy RME or indicate that lower energy halo stars have higher metallicity. Stars which comprise the "Gaia Sausage" velocity structure are a combination of the components identified in this work.

Although Blue Horizontal Branch (BHB) stars are commonly used to trace halo substructure, the stars bluer than (g-r)<-0.3 are ignored due to the difficulty in determining their absolute magnitudes. The blue extention of the horizontal branch (HBX) includes BHB tail stars and Extreme Horizontal Branch (EHB) stars. We present a method for identifying HBX stars in the field, using spectra and photometry from the Sloan Digital Sky Survey Data Release 14 (SDSS DR14). We derive an estimate for the absolute magnitudes of BHB tail and EHB stars as a function of color, and use this relationship to calculate distances. We identify an overdensity of HBX stars that appears to trace the northern end of the Hercules-Aquila Cloud (HAC). We identify three stars that are likely part of a tidal stream, but this is not enough stars to explain the observed overdensity. Combining SDSS data with Gaia DR2 proper motions allows us to show that the orbits of the majority of the HBX stars in the overdensity are on high eccentricity orbits similar to those in the Virgo Radial Merger/Gaia-Enceladus/Gaia Sausage structure, and that the overdensity of high eccentricity orbits extends all the way to the Virgo Overdensity. We use stellar kinematics to separate the HBX stars into disk stars andhalo stars. The halo stars are primarily on highly eccentric (radial) orbits. The fraction of HBX stars that are EHBs is highest in the disk population and lowest in the low eccentricity halo stars., Comment: 22 pages, 12 figures, accepted for publication in The Astrophysical Journal

We identify shell structures in the Milky Way for the first time. We find 2 shells in the Virgo Overdensity (VOD) region and 2 shells in the Hercules Aquila Cloud (HAC) region using Sloan Digital Sky Survey, Gaia, and LAMOST data. These shell stars are a subset of the substructure previously identified as the Virgo Radial Merger (VRM). Timing arguments for these shells indicate that their progenitor dwarf galaxy passed through the Galactic center 2.7 +/- 0.2 Gyr ago. Based on the time of collision, it is also possible that the VRM is related to the phenomenon that created phase-space spirals in the vertical motion of the disk and/or the Splash, and could have caused a burst of star formation in the inner disk. We analyze phase mixing in a collection of radial merger N-body simulations, and find that shell structure similar to that observed in Milky Way data disappears by 5 Gyr after collision with the Galactic center. The method used to calculate the merger time of the VRM was able to reliably recover the correct merger times for these simulations. Previous work supports the idea that the VRM and the Gaia Sausage/Gaia-Enceladus Merger are the same. However, the Gaia Sausage is widely believed to be 8--11 Gyr old. The disparate ages could be reconciled if the larger age is associated with an infall time when the progenitor crossed the virial radius; we do not constrain the time at which the progenitor became bound to the Milky Way. Alternatively, the Gaia Sausage could be younger than previously thought.

We suggest that the Virgo Overdensity (VOD) of stars in the stellar halo is the result of a radial dwarf galaxy merger that we call the Virgo Radial Merger. Because the dwarf galaxy passed very near to the Galactic center, the debris has a large range of energies but nearly zero $L_z$ angular momentum. The debris appears to extend from 5 to 50 kpc from the Sun in the Virgo region. We connect different parts of this merger debris to the Perpendicular and Parallel Streams (the Virgo Stellar Stream is associated with either or both of these streams), the Hercules-Aquila Cloud (HAC), and possibly the Eridanus Phoenix Overdensity (EriPhe). This radial merger can explain the majority of the observed moving groups of RR Lyrae and blue horizontal branch stars that have previously been identified in Virgo. This merger also produces debris in the Solar neighborhood similar to that identified as the Gaia-Enceladus or it Gaia-sausage merger. Orbits are provided for components of the Virgo Radial Merger progenitor and for debris that appears to be related to the Cocytos Stream, which was also recovered in the Virgo region., Comment: 17 pages, 5 figures, submitted to AAS Interstellar Matter and the Local Universe corridor