Search for: All records

Abstract Discoveries of gaps in data have been important in astrophysics. For example, there are kinematic gaps opened by resonances in dynamical systems, or exoplanets of a certain radius that are empirically rare. A gap in a data set is a kind of anomaly, but in an unusual sense: instead of being a single outlier data point, situated far from other data points, it is a region of the space, or a set of points, that is anomalous compared to its surroundings. Gaps are both interesting and hard to find and characterize, especially when they have nontrivial shapes. We present in this paper a statistic that can be used to estimate the (local) “gappiness” of a point in the data space. It uses the gradient and Hessian of the density estimate (and thus requires a twice-differentiable density estimator). This statistic can be computed at (almost) any point in the space and does not rely on optimization; it allows us to highlight underdense regions of any dimensionality and shape in a general and efficient way. We illustrate our method on the velocity distribution of nearby stars in the Milky Way disk plane, which exhibits gaps that could originate from different processes. Identifying and characterizing those gaps could help determine their origins. We provide in an appendix implementation notes and additional considerations for finding underdensities in data, using critical points and the properties of the Hessian of the density. 7 7 A Python implementation of t methods presented here is available at https://github.com/contardog/FindTheGap . 

Abstract The Galactic bulge is 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 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 ofR₀= 8108 ± 106_stat± 93_syspc, consistent with direct measurements 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 theV_Rquadrupole from the bar’s near side is reflected with respect to the Galactic center, indicating that the bar is bisymmetric and aligned with the inner disk. We also find that the vertical heightV_Zmap has no major structure in the region of the Galactic bulge, which is inconsistent with a current episode of bar buckling. Finally, we demonstrate withN-body simulations that distance uncertainty plays a factor in the alignment of the major and kinematic axes of the bar, necessitating caution when interpreting results for distant stars. 

Abstract In the coming decade, thousands of stellar streams will be observed in the halos of external galaxies. What fundamental discoveries will we make about dark matter from these streams? As a first attempt to look at these questions, we model Magellan/Megacam imaging of the Centaurus A (Cen A) disrupting dwarf companion Dwarf 3 (Dw3) and its associated stellar stream, to find out what can be learned about the Cen A dark matter halo. We develop a novel external galaxy stream-fitting technique and generate model stellar streams that reproduce the stream morphology visible in the imaging. We find that there are many viable stream models that fit the data well, with reasonable parameters, provided that Cen A has a halo mass larger than M 200 > 4.70 × 10 12 M ⊙ . There is a second stream in Cen A’s halo that is also reproduced within the context of this same dynamical model. However, stream morphology in the imaging alone does not uniquely determine the mass or mass distribution for the Cen A halo. In particular, the stream models with high likelihood show covariances between the inferred Cen A mass distribution, the inferred Dw3 progenitor mass, the Dw3 velocity, and the Dw3 line-of-sight position. We show that these degeneracies can be broken with radial-velocity measurements along the stream, and that a single radial velocity measurement puts a substantial lower limit on the halo mass. These results suggest that targeted radial-velocity measurements will be critical if we want to learn about dark matter from extragalactic stellar streams. 

ABSTRACT Gaia Data Release 2 revealed that the Milky Way contains significant indications of departures from equilibrium in the form of asymmetric features in the phase space density of stars in the Solar neighbourhood. One such feature is the z–vz phase spiral, interpreted as the response of the disc to the influence of a perturbation perpendicular to the disc plane, which could be external (e.g. a satellite) or internal (e.g. the bar or spiral arms). In this work, we use Gaia Data Release 3 to dissect the phase spiral by dividing the local data set into groups with similar azimuthal actions, Jϕ, and conjugate angles, θϕ, which selects stars on similar orbits and at similar orbital phases, thus having experienced similar perturbations in the past. These divisions allow us to explore areas of the Galactic disc larger than the surveyed region. The separation improves the clarity of the z–vz phase spiral and exposes changes to its morphology across the different action-angle groups. In particular, we discover a transition to two armed ‘breathing spirals’ in the inner Milky Way. We conclude that the local data contain signatures of not one, but multiple perturbations with the prospect to use their distinct properties to infer the properties of the interactions that caused them. 

ABSTRACT We present a novel method for constraining the length of the Galactic bar using 6D phase-space information to directly integrate orbits. We define a pseudo-length for the Galactic bar, named RFreq, based on the maximal extent of trapped bar orbits. We find the RFreq measured from orbits is consistent with the RFreq of the assumed potential only when the length of the bar and pattern speed of said potential is similar to the model from which the initial phase-space coordinates of the orbits are derived. Therefore, one can measure the model’s or the Milky Way’s bar length from 6D phase-space coordinates by determining which assumed potential leads to a self-consistent measured RFreq. When we apply this method to ≈210 000 stars in APOGEE DR17 and Gaia eDR3 data, we find a consistent result only for potential models with a dynamical bar length of ≈3.5 kpc. We find the Milky Way’s trapped bar orbits extend out to only ≈3.5 kpc, but there is also an overdensity of stars at the end of the bar out to 4.8 kpc which could be related to an attached spiral arm. We also find that the measured orbital structure of the bar is strongly dependent on the properties of the assumed potential. 

ABSTRACT In this work, we present two new ∼109 particle self-consistent simulations of the merger of a Sagittarius-like dwarf galaxy with a Milky Way (MW)-like disc galaxy. One model is a violent merger creating a thick disc, and a Gaia–Enceladus/Sausage-like remnant. The other is a highly stable disc which we use to illustrate how the improved phase space resolution allows us to better examine the formation and evolution of structures that have been observed in small, local volumes in the MW, such as the z−vz phase spiral and clustering in the vR−vϕ plane when compared to previous works. The local z−vz phase spirals are clearly linked to the global asymmetry across the disc: we find both 2-armed and 1-armed phase spirals, which are related to breathing and bending behaviours, respectively. Hercules-like moving groups are common, clustered in vR−vϕ in local data samples in the simulation. These groups migrate outwards from the inner galaxy, matching observed metallicity trends even in the absence of a galactic bar. We currently release the best-fitting ‘present-day’ merger snapshots along with the unperturbed galaxies for comparison. 

Abstract Signatures of vertical disequilibrium have been observed across the Milky Way’s (MW’s) disk. These signatures manifest locally as unmixed phase spirals inz–v_zspace (“snails-in-phase”), and globally as nonzero meanzandv_z, wrapping around the disk into physical spirals in thex–yplane (“snails-in-space”). We explore the connection between these local and global spirals through the example of a satellite perturbing a test-particle MW-like disk. We anticipate our results to broadly apply to any vertical perturbation. Using az–v_zasymmetry metric, we demonstrate that in test-particle simulations: (a) multiple local phase-spiral morphologies appear when stars are binned by azimuthal actionJ_ϕ, excited by a single event (in our case, a satellite disk crossing); (b) these distinct phase spirals are traced back to distinct disk locations; and (c) they are excited at distinct times. Thus, local phase spirals offer a global view of the MW’s perturbation history from multiple perspectives. Using a toy model for a Sagittarius (Sgr)–like satellite crossing the disk, we show that the full interaction takes place on timescales comparable to orbital periods of disk stars withinR≲ 10 kpc. Hence such perturbations have widespread influence, which peaks in distinct regions of the disk at different times. This leads us to examine the ongoing MW–Sgr interaction. While Sgr has not yet crossed the disk (currently,z_Sgr≈ −6 kpc,v_z,Sgr≈ 210 km s⁻¹), we demonstrate that the peak of the impact has already passed. Sgr’s pull over the past 150 Myr creates a globalv_zsignature with amplitude ∝M_Sgr, which might be detectable in future spectroscopic surveys. 

Abstract Gaia DR2 has provided an unprecedented wealth of information about the positions and motions of stars in our Galaxy, and has highlighted the degree of disequilibria in the disc. As we collect data over a wider area of the disc it becomes increasingly appealing to start analysing stellar actions and angles, which specifically label orbit space, instead of their current phase space location. Conceptually, while $$\bar{x}$$ and $$\bar{v}$$ tell us about the potential and local interactions, grouping in action puts together stars that have similar frequencies and hence similar responses to dynamical effects occurring over several orbits. Grouping in actions and angles refines this further to isolate stars which are travelling together through space and hence have shared histories. Mixing these coordinate systems can confuse the interpretation. For example, it has been suggested that by moving stars to their guiding radius, the Milky Way spiral structure is visible as ridge-like overdensities in the Gaia data (Khoperskov et al. 2020). However, in this work, we show that these features are in fact the known kinematic moving groups, both in the Lz − φ and the vR − vφ planes. Using simulations we show how this distinction will become even more important as we move to a global view of the Milky Way. As an example, we show that the radial velocity wave seen in the Galactic disc in Gaia and APOGEE should become stronger in the action-angle frame, and that it can be reproduced by transient spiral structure.