We present astroplan-an open source, open development, Astropy affiliated package for ground-based observation planning and scheduling in Python. astroplan is designed to provide efficient access to ...common observational quantities such as celestial rise, set, and meridian transit times and simple transformations from sky coordinates to altitude-azimuth coordinates without requiring a detailed understanding of astropy's implementation of coordinate systems. astroplan provides convenience functions to generate common observational plots such as airmass and parallactic angle as a function of time, along with basic sky (finder) charts. Users can determine whether or not a target is observable given a variety of observing constraints, such as airmass limits, time ranges, Moon illumination/separation ranges, and more. A selection of observation schedulers are included that divide observing time among a list of targets, given observing constraints on those targets. Contributions to the source code from the community are welcome.
ABSTRACT The Ophiuchus stellar stream is peculiar: (1) its length is short given the age of its constituent stars, and (2) several probable member stars have dispersions in sky position and velocity ...that far exceed those seen within the stream. The stream's proximity to the Galactic center suggests that its dynamical history is significantly influenced by the Galactic bar. We explore this hypothesis with models of stream formation along orbits consistent with Ophiuchus' properties in a Milky Way potential model that includes a rotating bar. In all choices for the rotation parameters of the bar, orbits fit to the stream are strongly chaotic. Mock streams generated along these orbits qualitatively match the observed properties of the stream: because of chaos, stars stripped early generally form low-density, high-dispersion "fans" leaving only the most recently disrupted material detectable as a strong over-density. Our models predict that there should be a significant amount of low-surface-brightness tidal debris around the stream with a complex phase-space morphology. The existence of or lack of these features could provide interesting constraints on the Milky Way bar and would rule out formation scenarios for the stream. This is the first time that chaos has been used to explain the properties of a stellar stream and is the first demonstration of the dynamical importance of chaos in the Galactic halo. The existence of long, thin streams around the Milky Way, presumably formed along non- or weakly chaotic orbits, may represent only a subset of the total population of disrupted satellites.
In simple models of the Milky Way, tidally disrupting satellites produce long and thin-nearly one-dimensional-stellar streams. Using astrometric data from the Gaia second data release and photometry ...from the Dark Energy Survey, we demonstrate that the Jhelum stream, a stellar stream in the inner halo, is a two-dimensional structure. The spatial distribution of highly probable Jhelum members reveals a dense thin component and an associated diffuse, spatially offset component. These two spatial components have indistinguishable proper motions (at ∼ 1 mas yr−1 level) and a similar ratio of blue straggler to blue horizontal branch stars, which indicates a common origin for the two components. The best-fit orbit to the narrow component (pericenter 8 kpc, apocenter 24 kpc), however, does not explain the wide component of the Jhelum stream. On the other hand, an older orbital wrap of Jhelum's orbit traces the Indus stream, indicating a possible connection between these two structures and additional complexity in Jhelum's formation. Substructure in the Jhelum progenitor or precession of its tidal debris in the Milky Way potential may explain the observed structure of Jhelum. Future spectroscopic data will enable discrimination between these "nature" and "nurture" formation scenarios. Jhelum adds to the growing list of cold stellar streams that display complex morphologies and promise to reveal the dynamical history of the Milky Way.
Given sparse or low-quality radial velocity measurements of a star, there are often many qualitatively different stellar or exoplanet companion orbit models that are consistent with the data. The ...consequent multimodality of the likelihood function leads to extremely challenging search, optimization, and Markov chain Monte Carlo (MCMC) posterior sampling over the orbital parameters. Here we create a custom Monte Carlo sampler for sparse or noisy radial velocity measurements of two-body systems that can produce posterior samples for orbital parameters even when the likelihood function is poorly behaved. The six standard orbital parameters for a binary system can be split into four nonlinear parameters (period, eccentricity, argument of pericenter, phase) and two linear parameters (velocity amplitude, barycenter velocity). We capitalize on this by building a sampling method in which we densely sample the prior probability density function (pdf) in the nonlinear parameters and perform rejection sampling using a likelihood function marginalized over the linear parameters. With sparse or uninformative data, the sampling obtained by this rejection sampling is generally multimodal and dense. With informative data, the sampling becomes effectively unimodal but too sparse: in these cases we follow the rejection sampling with standard MCMC. The method produces correct samplings in orbital parameters for data that include as few as three epochs. The Joker can therefore be used to produce proper samplings of multimodal pdfs, which are still informative and can be used in hierarchical (population) modeling. We give some examples that show how the posterior pdf depends sensitively on the number and time coverage of the observations and their uncertainties.
Abstract
The dark matter halos that surround Milky Way–like galaxies in cosmological simulations are, to first order, triaxial. Nearly 30 yr ago it was predicted that such triaxial dark matter halos ...should exhibit steady figure rotation or tumbling motions for durations of several gigayears. The angular frequency of figure rotation predicted by cosmological simulations is described by a log-normal distribution of pattern speed Ω
p
with a median value 0.15
h
km s
−1
kpc
−1
(∼0.15
h
rad Gyr
−1
∼ 9°
h
Gyr
−1
) and a width of 0.83
h
km s
−1
kpc
−1
. These pattern speeds are so small that they have generally been considered both unimportant and undetectable. In this work we show that even extremely slow figure rotation can significantly alter the structure of extended stellar streams produced by the tidal disruption of satellites in the Milky Way halo. We simulate the behavior of a Sagittarius-like polar tidal stream in triaxial dark matter halos with different shapes, when the halos are rotated about the three principal axes. For pattern speeds typical of cosmological halos, we demonstrate, for the first time, that a Sagittarius-like tidal stream would be altered to a degree that is detectable even with current observations. This discovery will potentially allow for a future measurement of figure rotation of the Milky Way’s dark matter halo, perhaps enabling the first evidence of this relatively unexplored prediction of cold dark matter.
Abstract
Stellar streams in the Galactic halo are useful probes of the assembly of galaxies like the Milky Way. Many tidal stellar streams that have been found in recent years are accompanied by a ...known progenitor globular cluster or dwarf galaxy. However, the Orphan–Chenab (OC) stream is one case where a relatively narrow stream of stars has been found without a known progenitor. In an effort to find the parent of the OC stream, we use astrometry from the early third data release of ESA’s Gaia mission (Gaia EDR3) and radial velocity information from the Sloan Digital Sky Survey (SDSS)-IV Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey to find up to 13 stars that are likely members of the OC stream. We use the APOGEE survey to study the chemical nature (for up to 10 stars) of the OC stream in the
α
(O, Mg, Ca, Si, Ti, and S), odd-
Z
(Al, K, and V), Fe-peak (Fe, Ni, Mn, Co, and Cr), and neutron-capture (Ce) elemental groups. We find that the stars that make up the OC stream are not consistent with a monometallic population and have a median metallicity of −1.92 dex with a dispersion of 0.28 dex. Our results also indicate that the α elements are depleted compared to the known Milky Way populations and that its Mg/Al abundance ratio is not consistent with second-generation stars from globular clusters. The detailed chemical pattern of these stars, namely the
α
/Fe–Fe/H plane and the metallicity distribution, indicates that the OC stream progenitor is very likely to be a dwarf spheroidal galaxy with a mass of ∼10
6
M
⊙
.
We use photometry from the DECam Legacy Survey to detect candidate tidal tails extending ∼5° on either side of the Palomar 13 globular cluster. The tails are aligned with the proper motion of Palomar ...13 and are consistent with its old, metal-poor stellar population. We identify three RR Lyrae stars (RRLs) that are plausibly associated with the tails, in addition to four previously known in the cluster. From these RRLs, we find that the mean distance to the cluster and tails is 23.6 0.2 kpc and estimate the total (initial) luminosity of the cluster to be , consistent with previous claims that its initial luminosity was higher than its current luminosity. Combined with previously determined proper motion and radial velocity measurements of the cluster, we find that Palomar 13 is on a highly eccentric orbit (e ∼ 0.8) with a pericenter of ∼9 kpc and an apocenter of ∼69 kpc, and a recent pericentric passage of the cluster ∼75 Myr ago. We note a prominent linear structure in the interstellar dust map that runs parallel to the candidate tidal features, but conclude that reddening due to dust is unlikely to account for the structure that we observe. If confirmed, the Palomar 13 stellar stream would be one of very few streams with a known progenitor system, making it uniquely powerful for studying the disruption of globular clusters, the formation of the stellar halo, and the distribution of matter within our Galaxy.
We employ the earlier published proper motions of the newly discovered Antlia 2 dwarf galaxy derived from Gaia data to calculate its orbital distribution in the cosmologically recent past. Using ...these observationally motivated orbits, we calculate the effect of the Antlia 2 dwarf galaxy on the outer H i disk of the Milky Way, using both test particle and smoothed particle hydrodynamics simulations. We find that orbits with low pericenters, ∼10 kpc, produce disturbances that match the observed outer H i disk perturbations. We have independently recalculated the proper motion of the Antlia 2 dwarf from Gaia data and found a proper motion of ( cosδ, δ) = (−0.068, 0.032) (0.023, −0.031) mas yr−1, which agrees with results from Torrealba et al. within the errors, but gives lower mean pericenters, e.g., ∼15 kpc for our fiducial model of the Milky Way. We also show that the Sagittarius dwarf galaxy interaction does not match the observed perturbations in the outer gas disk. Thus, Antlia 2 may be the driver of the observed large perturbations in the outer gas disk of the Galaxy. The current location of the Antlia 2 dwarf galaxy closely matches that predicted by an earlier dynamical analysis of the dwarf galaxy that drove ripples in the outer Galaxy, and, in particular, its orbit is nearly coplanar to the Galactic disk. If the Antlia 2 dwarf galaxy is responsible for the perturbations in the outer Galactic disk, it would have a specific range of proper motions that we predict here; this can be tested soon with Gaia DR-3 and Gaia DR-4 data.
ABSTRACT
We analyse two binary systems containing giant stars, V723 Mon (‘the Unicorn’) and 2M04123153+6738486 (‘the Giraffe’). Both giants orbit more massive but less luminous companions, previously ...proposed to be mass-gap black holes. Spectral disentangling reveals luminous companions with star-like spectra in both systems. Joint modelling of the spectra, light curves, and spectral energy distributions robustly constrains the masses, temperatures, and radii of both components: the primaries are luminous, cool giants ($T_{\rm eff,\, giant} = 3800$ and $4000\, \rm K$, $R_{\rm giant}= 22.5$ and $25\, {\rm R}_{\odot }$) with exceptionally low masses ($M_{\rm giant} \approx 0.4\, {\rm M}_{\odot }$) that likely fill their Roche lobes. The secondaries are only slightly warmer subgiants ($T_{\rm eff,\, 2} = 5800$ and $5150\, \rm K$, $R_2= 8.3$ and $9\, {\rm R}_{\odot }$) and thus are consistent with observed UV limits that would rule out main-sequence stars with similar masses ($M_2 \approx 2.8$ and ${\approx}1.8\, {\rm M}_{\odot }$). In the Unicorn, rapid rotation blurs the spectral lines of the subgiant, making it challenging to detect even at wavelengths where it dominates the total light. Both giants have surface abundances indicative of CNO processing and subsequent envelope stripping. The properties of both systems can be reproduced by binary evolution models in which a $1{-}2\, {\rm M}_{\odot }$ primary is stripped by a companion as it ascends the giant branch. The fact that the companions are also evolved implies either that the initial mass ratio was very near unity, or that the companions are temporarily inflated due to rapid accretion. The Unicorn and Giraffe offer a window into into a rarely observed phase of binary evolution preceding the formation of wide-orbit helium white dwarfs, and eventually, compact binaries containing two helium white dwarfs.
The Monoceros Ring (also known as the Galactic Anticenter Stellar Structure) and A13 are stellar overdensities at estimated heliocentric distances of d ∼ 11 kpc and 15 kpc observed at low Galactic ...latitudes toward the anticenter of our Galaxy. While these overdensities were initially thought to be remnants of a tidally disrupted satellite galaxy, an alternate scenario is that they are composed of stars from the Milky Way (MW) disk kicked out to their current location due to interactions between a satellite galaxy and the disk. To test this scenario, we study the stellar populations of the Monoceros Ring and A13 by measuring the number of RR Lyrae and M giant stars associated with these overdensities. We obtain low-resolution spectroscopy for RR Lyrae stars in the two structures and measure radial velocities to compare with previously measured velocities for M giant stars in the regions of the Monoceros Ring and A13, to assess the fraction of RR Lyrae to M giant stars (fRR:MG) in A13 and Mon/GASS. We perform velocity modeling on 153 RR Lyrae stars (116 in the Monoceros Ring and 37 in A13) and find that both structures have very low fRR:MG. The results support a scenario in which stars in A13 and Mon/GASS formed in the MW disk. We discuss a possible association between Mon/GASS, A13, and the Triangulum-Andromeda overdensity based on their similar velocity distributions and fRR:MG.