BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Jason L. Sanders

Jason L. Sanders appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.GAastro-ph.SRastro-ph.COastro-ph.HE

Catalog footprint

What is connected

20works
4topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

20 published item(s)

preprint2026arXiv

Towards Galactic Archaeology with Inferred Ages of Giant Stars From Gaia Spectra

In the era of Gaia, the accurate determination of stellar ages is transforming Galactic archaeology. We demonstrate the feasibility of inferring stellar ages from Gaia's RVS spectra and the BP/RP (XP) spectrophotometric data, specifically for red giant branch and high-mass red clump stars. We successfully train two machine learning models, dubbed SIDRA: Stellar age Inference Derived from Gaia spectRA to predict the age. The SIDRA-RVS model uses the RVS spectra and SIDRA-XP the stellar parameters obtained from the XP spectra. Both models use BINGO, an APOGEE-derived stellar age as the training data. SIDRA-RVS estimates ages of stars whose age is around $τ_\mathrm{BINGO}=10$~Gyr with a standard deviation of residuals of $\sim$ 0.12 dex in the unseen test dataset, while SIDRA-XP achieves higher precision with residuals $\sim$ 0.064 dex for stars around $τ_\mathrm{BINGO}=10$~ Gyr. Since SIDRA-XP outperforms SIDRA-RVS, we apply SIDRA-XP to analyse the ages for 2,218,154 stars. This allowed us to map the chronological and chemical properties of Galactic disc stars, reproducing the known distinct features such as the Gaia-Sausage-Enceladus merger and a potential gas-rich interaction event linked to the first infall of the Sagittarius dwarf galaxy. This study demonstrates that machine learning techniques applied to Gaia's spectra can provide valuable individual age information, particularly for giant stars, thereby enhancing our understanding of the Milky Way's formation and evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2024arXiv

Nuclear rings are the inner edge of a gap around the Lindblad Resonance

Gaseous nuclear rings are large-scale coherent structures commonly found at the centres of barred galaxies. We propose that they are an accumulation of gas at the inner edge of an extensive gap that forms around the Inner Lindblad Resonance (ILR). The gap initially opens because the bar potential excites strong trailing waves near the ILR, which remove angular momentum from the gas disc and transport the gas inwards. The gap then widens because the bar potential continuously excites trailing waves at the inner edge of the gap, which remove further angular momentum, moving the edge further inwards until it stops at a distance of several wavelengths from the ILR. The gas accumulating at the inner edge of the gap forms the nuclear ring. The speed at which the gap edge moves and its final distance from the ILR strongly depend on the sound speed, explaining the puzzling dependence of the nuclear ring radius on the sound speed in simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Information content of BP/RP spectra in Gaia DR3

Gaia Data Release 3 has provided the astronomical community with the largest stellar spectroscopic survey to date ($>$ 220 million sources). The low resolution (R$\sim$50) blue photometer (BP) and red photometer (RP) spectra will allow for the estimation of stellar atmospheric parameters such as effective temperature, surface gravity and metallicity. We create mock Gaia BP/RP spectra and use Fisher information matrices to probe the resolution limit of stellar parameter measurements using BP/RP spectra. The best-case scenario uncertainties that this analysis provides are then used to produce a mock-observed stellar population in order to evaluate the false positive rate (FPR) of identifying extremely metal-poor (EMP) stars. We conclude that the community will be able to confidently identify metal-poor stars at magnitudes brighter than $G = 16$ using BP/RP spectra. At fainter magnitudes true detections will start to be overwhelmed by false positives. When adopting the commonly-used $G < 14$ limit for metal-poor star searches, we find a FPR for the low-metallicity regimes [Fe/H] $<$ -2, -2.5 and -3 of just 14$\%$, 33$\%$ and 56$\%$ respectively, offering the potential for significant improvements on previous targeting campaigns. Additionally, we explore the chemical sensitivity obtainable directly from BP/RP spectra for Carbon and $α$-elements. We find an absolute Carbon abundance uncertainty of $σ_{A(C)} < 1$ dex for Carbon-enriched metal-poor (CEMP) stars, indicating the potential to identify a CEMP stellar population for follow-up confirmation with higher resolution spectroscopy. Finally, we find that large uncertainties in $α$-element abundance measurements using BP/RP spectra means that efficiently obtaining these abundances will be challenging.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Mira variables in the Milky Way's nuclear stellar disc: discovery and classification

The properties of the Milky Way's nuclear stellar disc give crucial information on the epoch of bar formation. Mira variables are promising bright candidates to study the nuclear stellar disc, and through their period-age relation dissect its star formation history. We report on a sample of $1782$ Mira variable candidates across the central $3\times3\,\mathrm{deg}^2$ of the Galaxy using the multi-epoch infrared VISTA Variables in Via Lactea (VVV) survey. We describe the algorithms employed to select candidate variable stars and then model their light curves using periodogram and Gaussian process methods. By combining with WISE, 2MASS and other archival photometry, we model the multi-band light curves to refine the periods and inspect the amplitude variation between different photometric bands. The infrared brightness of the Mira variables means many are too bright and missed by VVV. However, our sample follows a well-defined selection function as expected from artificial star tests. The multi-band photometry is modelled using stellar models with circumstellar dust that characterise the mass loss rates. We demonstrate how $\gtrsim90$ per cent of our sample is consistent with O-rich chemistry. Comparison to period-luminosity relations demonstrates that the bulk of the short period stars are situated at the Galactic Centre distance. Many of the longer period variables are very dusty, falling significantly under the O-rich Magellanic Cloud and solar neighbourhood period-luminosity relations and exhibit high mass-loss rates of $\sim2.5\times10^{-5}M_\odot\,\mathrm{yr}^{-1}$. The period distribution appears consistent with the nuclear stellar disc forming $\gtrsim8\,\mathrm{Gyr}$ ago although it is not possible to disentangle the relative contributions of the nuclear stellar disc and the contaminating bulge.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The extinction law in the inner $3\times3$ deg$^2$ of the Milky Way and the red clump absolute magnitude in the inner bar-bulge

The extinction law from $0.9$ to $8$ microns in the inner $3\times3$ deg$^2$ of the Milky Way is measured using data from VISTA Variables in the Via Lactea, GLIMPSE and WISE. Absolute extinction ratios are found by requiring that the observed red clump density peaks at the GRAVITY collaboration distance to the Galactic centre. When combined with selective extinction ratios measured from the bulge giant colour-colour diagrams, we find an extinction law of $A_Z:A_Y:A_J:A_H:A_{K_s}:A_{W1}:A_{[3.6]}:A_{[4.5]}:A_{W2}:A_{[5.8]}:A_{[8.0]} =7.19(0.30):5.11(0.20):3.23(0.11):1.77(0.04):1:0.54(0.02):0.46(0.03):0.34(0.03):0.32(0.03):0.24(0.04):0.28(0.03)$ valid for low extinctions where non-linearities are unimportant. These results imply an extinction law from the Rayleigh Jeans colour excess method (RJCE) of $A_{K_s}=0.677(H-[4.5]-0.188)$. We find little evidence for significant selective extinction ratio variation over the inspected region (around $5\%$). Assuming the absolute extinction ratios do not vary across the inspected region gives an independent measurement of the absolute $K_s$ magnitude of the red clump at the Galactic Centre of $(-1.61\pm0.07)\,\mathrm{mag}$. This is very similar to the value measured for solar neighbourhood red clump stars giving confidence in the use of red clump stars as standard candles across the Galaxy. As part of our analysis, we inspect the completeness of PSF photometry from the VVV survey using artificial star tests, finding $90\%$ completeness at $K_s\approx16 \,(17)$ in high (low) density regions and good agreement with the number counts with respect to the GALACTICNUCLEUS and DECAPS catalogues over small regions of the survey.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Variable star classification across the Galactic bulge and disc with the VISTA Variables in the Vía Láctea survey

We present VIVACE, the VIrac VAriable Classification Ensemble, a catalogue of variable stars extracted from an automated classification pipeline for the Vista Variables in the Vía Láctea (VVV) infrared survey of the Galactic bar/bulge and southern disc. Our procedure utilises a two-stage hierarchical classifier to first isolate likely variable sources using simple variability summary statistics and training sets of non-variable sources from the Gaia early third data release, and then classify candidate variables using more detailed light curve statistics and training labels primarily from OGLE and VSX. The methodology is applied to point-spread-function photometry for $\sim490$ million light curves from the VIRAC v2 astrometric and photometric catalogue resulting in a catalogue of $\sim1.4$ million likely variable stars, of which $\sim39,000$ are high-confidence (classification probability $>0.9$) RR Lyrae ab stars, $\sim8000$ RR Lyrae c/d stars, $\sim187,000$ detached/semi-detached eclipsing binaries, $\sim18,000$ contact eclipsing binaries, $\sim1400$ classical Cepheid variables and $\sim2200$ Type II Cepheid variables. Comparison with OGLE-4 suggests a completeness of around $90\,\%$ for RRab and $\lesssim60\%$ for RRc/d, and a misclassification rate for known RR Lyrae stars of around $1\%$ for the high confidence sample. We close with two science demonstrations of our new VIVACE catalogue: first, a brief investigation of the spatial and kinematic properties of the RR Lyrae stars within the disc/bulge, demonstrating the spatial elongation of bar-bulge RR Lyrae stars is in the same sense as the more metal-rich red giant population whilst having a slower rotation rate of $\sim40\,\mathrm{km\,s}^{-1}\mathrm{kpc}^{-1}$; and secondly, an investigation of the Gaia EDR3 parallax zeropoint using contact eclipsing binaries across the Galactic disc plane and bulge.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The Biggest Splash

Using a large sample of bright nearby stars with accurate Gaia Data Release 2 astrometry and auxiliary spectroscopy we map out the properties of the principle Galactic components such as the "thin" and "thick" discs and the halo. We show that in the Solar neighborhood, there exists a large population of metal-rich ([Fe/H]>-0.7) stars on highly eccentric orbits. By studying the evolution of elemental abundances, kinematics and stellar ages in the plane of azimuthal velocity v_phi and metallicity [Fe/H], we demonstrate that this metal-rich halo-like component, which we dub the Splash, is linked to the alpha-rich (or "thick") disc. Splash stars have little to no angular momentum and many are on retrograde orbits. They are predominantly old, but not as old as the stars deposited into the Milky Way in the last major merger. We argue, in agreement with several recent studies, that the Splash stars may have been born in the Milky Way's proto-disc prior to the massive ancient accretion event which drastically altered their orbits. We can not, however, rule out other (alternative) formation channels. Taking advantage of the causal connection between the merger and the Splash, we put constraints of the epoch of the last massive accretion event to have finished 9.5 Gyr ago. The link between the local metal-rich and metal-poor retrograde stars is confirmed using a large suite of cutting-edge numerical simulations of the Milky Way's formation.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A review of action estimation methods for galactic dynamics

We review the available methods for estimating actions, angles and frequencies of orbits in both axisymmetric and triaxial potentials. The methods are separated into two classes. Unless an orbit has been trapped by a resonance, convergent, or iterative, methods are able to recover the actions to arbitrarily high accuracy given sufficient computing time. Faster non-convergent methods rely on the potential being sufficiently close to a separable potential and the accuracy of the action estimate cannot be improved through further computation. We critically compare the accuracy of the methods and the required computation time for a range of orbits in an axisymmetric multi-component Galactic potential. We introduce a new method for estimating actions that builds on the adiabatic approximation of Schönich & Binney (2012) and discuss the accuracy required for the actions, angles and frequencies using suitable distribution functions for the thin and thick discs, the stellar halo and a star stream. We conclude that for studies of the disc and smooth halo component of the Milky Way the most suitable compromise between speed and accuracy is the Stäckel Fudge, whilst when studying streams the non-convergent methods do not offer sufficient accuracy and the most suitable method is computing the actions from an orbit integration via a generating function. All the software used in this study can be downloaded from https://github.com/jls713/tact.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Dynamics of stream-subhalo interactions

We develop a formalism for modelling the impact of dark matter subhaloes on cold thin streams. Our formalism models the formation of a gap in a stream in angle-frequency space and is able to handle general stream and impact geometry. We analyse an $N$-body simulation of a cold stream formed from a progenitor on an eccentric orbit in an axisymmetric potential, which is perturbed by a direct impact from a $10^8 M_\odot$ subhalo, and produce a complete generative model of the perturbed stream that matches the simulation well at a range of times. We show how the results in angle-frequency space can be related to physical properties of the gaps and that previous results for more constrained simulations are recovered. We demonstrate how our results are dependent upon the mass of the subhalo and the location of the impact along the stream. We find that gaps formed far downstream grow more rapidly than those closer to the progenitor due to the more ordered nature of the stream members far from the progenitor. Additionally, we show that the minimum gap density plateaus in time at a value that decreases with increasing subhalo mass.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Indirect Dark Matter Detection for Flattened Dwarf Galaxies

Gamma-ray experiments seeking to detect evidence of dark matter annihilation in dwarf spheroidal galaxies require knowledge of the distribution of dark matter within these systems. We analyze the effects of flattening on the annihilation (J) and decay (D) factors of dwarf spheroidal galaxies with both analytic and numerical methods. Flattening has two consequences: first, there is a geometric effect as the squeezing (or stretching) of the dark matter distribution enhances (or diminishes) the J-factor; second, the line of sight velocity dispersion of stars must hold up the flattened baryonic component in the flattened dark matter halo. We provide analytic formulae and a simple numerical approach to estimate the correction to the J- and D-factors required over simple spherical modeling. The formulae are validated with a series of equilibrium models of flattened stellar distributions embedded in flattened dark-matter distributions. We compute corrections to the J- and D-factors for the Milky Way dwarf spheroidal galaxies under the assumption that they are all prolate or all oblate and find that the hierarchy of J-factors for the dwarf spheroidals is slightly altered (typical correction factors for an ellipticity of $0.4$ are $0.75$ for the oblate case and $1.6$ for the prolate case). We demonstrate that spherical estimates of the D-factors are very insensitive to the flattening and introduce uncertainties significantly less than the uncertainties in the D-factors from the other observables for all the dwarf spheroidals (for example, ${}^{+10\mathrm{\,percent}}_{-3\mathrm{\,percent}}$ for a typical ellipticity of $0.4$). We conclude by investigating the spread in correction factors produced by triaxial figures and provide uncertainties in the J-factors for the dwarf spheroidals using different physically-motivated assumptions for their intrinsic shape and axis alignments. (abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Mass estimators for flattened dispersion-supported galaxies

We investigate the reliability of mass estimators based on the observable velocity dispersion and half-light radius $R_\mathrm{h}$ for dispersion-supported galaxies. We show how to extend them to flattened systems and provide simple formulae for the mass within an ellipsoid under the assumption the dark matter density and the stellar density are stratified on the same self-similar ellipsoids. We demonstrate explicitly that the spherical mass estimators (Walker et al. 2009, Wolf et al. 2010) give accurate values for the mass within the half-light ellipsoid, provided $R_\mathrm{h}$ is replaced by its `circularized' analogue $R_\mathrm{h}\sqrt{1-ε}$. We provide a mathematical justification for this surprisingly simple and effective workaround. It means, for example, that the mass-to-light ratios are valid not just when the light and dark matter are spherically distributed, but also when they are flattened on ellipsoids of the same constant shape.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Stray, swing and scatter: angular momentum evolution of orbits and streams in aspherical potentials

In aspherical potentials orbital planes continuously evolve. The gravitational torques impel the angular momentum vector to precess, that is to slowly stray around the symmetry axis, and nutate, i.e. swing up and down periodically in the perpendicular direction. This familiar orbital pole motion - if detected and measured - can reveal the shape of the underlying gravitational potential, the quantity only crudely gauged in the Galaxy so far. Here we demonstrate that the debris poles of stellar tidal streams show a very similar straying and swinging behavior, and give analytic expressions to link the amplitude and the frequency of the pole evolution to the flattening of the dark matter distribution. While these results are derived for near-circular orbits, we show they are also valid for eccentric orbits. Most importantly, we explain how the differential orbital plane precession leads to the broadening of the stream and show that streams on polar orbits ought to scatter faster. We provide expressions for the stream width evolution as a function of the axisymmetric potential flattening and the angle from the symmetry plane and prove that our models are in good agreement with streams produced in N-body simulations. Interestingly, the same intuition applies to streams whose progenitors are on short or long-axis loops in a triaxial potential. Finally, we present a compilation of the Galactic cold stream data, and discuss how the simple picture developed here, along with stream modelling, can be used to constrain the symmetry axes and flattening of the Milky Way.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The number and size of subhalo-induced gaps in stellar streams

Ample observational capabilities exist today to detect the small density perturbations that low-mass dark matter subhaloes impart on stellar streams from disrupting Galactic satellites. In anticipation of these observations, we investigate the expected number and size of gaps by combining an analytic prescription for gap evolution on circular orbits with the flux of subhaloes near the stream. We explore the distribution of gap sizes and depths for a typical cold stream around the Milky Way and find that for a given stream age and gap depth, each subhalo mass produces a characteristic gap size. For a stream with an age of a few Gyr, orbiting at a distance of 10-20 kpc from the Galactic center, even modest subhaloes with a mass of $10^6-10^7 M_\odot$ produce gaps with sizes that are on the order of several degrees. We consider the number and distribution of gap sizes created by subhaloes with masses $10^5-10^9 M_\odot$, accounting for the expected depletion of subhaloes by the Milky Way disk, and present predictions for six cold streams around the Milky Way. For Pal 5, we forecast 0.7 gaps with a density depletion of at least 25\% and a typical gap size of $8^\circ$. Thus, there appears to be {\it no tension} between the recent non-detection of density depletions in the Pal 5 tidal tails and $Λ$CDM expectations. These predictions can be used to guide the scale of future gap searches.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Extended distribution functions for our Galaxy

We extend models of our Galaxy based on distribution functions (DFs) that are analytic functions of the action integrals to extended distribution functions (EDFs), which have an analytic dependence on metallicity as well. We use a simple, but physically-motivated, functional forms for the metallicity of the interstellar medium as a function of radius and time and for the star-formation rate, and a model for the diffusion of stars through phase space to suggest the required functional form of an EDF. We introduce a simple prescription for radial migration that preserves the overall profile of the disc while allowing individual stars to migrate throughout the disc. Our models explicitly consider the thin and thick discs as two distinct components separated in age. We show how an EDF can be used to incorporate realistic selection functions in models, and to construct mock catalogues of observed samples. We show that the selection function of the Geneva-Copenhagen Survey (GCS) biases in favour of young stars, which have atypically small random velocities. With the selection function taken into account our models produce good fits of the GCS data in chemo-dynamical space and the Gilmore and Reid (1983) density data. From our EDF, we predict the structure of the SEGUE G-dwarf sample. The kinematics are successfully predicted. The predicted metallicity distribution has too few stars with [Fe/H]$\simeq-0.5$ dex and too many metal-rich stars. A significant problem may be the lack of any chemical-kinematic correlations in our thick disc. We argue that EDFs will prove essential tools for the analysis of both observational data and sophisticated models of Galaxy formation and evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Self-consistent triaxial models

We present self-consistent triaxial stellar systems that have analytic distribution functions (DFs) expressed in terms of the actions. These provide triaxial density profiles with cores or cusps at the centre. They are the first self-consistent triaxial models with analytic DFs suitable for modelling giant ellipticals and dark haloes. Specifically, we study triaxial models that reproduce the Hernquist profile from Williams & Evans (2015), as well as flattened isochrones of the form proposed by Binney (2014). We explore the kinematics and orbital structure of these models in some detail. The models typically become more radially anisotropic on moving outwards, have velocity ellipsoids aligned in Cartesian coordinates in the centre and aligned in spherical polar coordinates in the outer parts. In projection, the ellipticity of the isophotes and the position angle of the major axis of our models generally changes with radius. So, a natural application is to elliptical galaxies that exhibit isophote twisting. As triaxial Stäckel models do not show isophote twists, our DFs are the first to generate mass density distributions that do exhibit this phenomenon, typically with a gradient of $\approx 10^\circ$/effective radius, which is comparable to the data. Triaxiality is a natural consequence of models that are susceptible to the radial orbit instability. We show how a family of spherical models with anisotropy profiles that transition from isotropic at the centre to radially anisotropic becomes unstable when the outer anisotropy is made sufficiently radial. Models with a larger outer anisotropy can be constructed but are found to be triaxial. We argue that the onset of the radial orbit instability can be identified with the transition point when adiabatic relaxation yields strongly triaxial rather than weakly spherical endpoints.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A fast algorithm for estimating actions in triaxial potentials

We present an approach to approximating rapidly the actions in a general triaxial potential. The method is an extension of the axisymmetric approach presented by Binney (2012), and operates by assuming that the true potential is locally sufficiently close to some Stäckel potential. The choice of Stäckel potential and associated ellipsoidal coordinates is tailored to each individual input phase-space point. We investigate the accuracy of the method when computing actions in a triaxial Navarro-Frenk-White potential. The speed of the algorithm comes at the expense of large errors in the actions, particularly for the box orbits. However, we show that the method can be used to recover the observables of triaxial systems from given distribution functions to sufficient accuracy for the Jeans equations to be satisfied. Consequently, such models could be used to build models of external galaxies as well as triaxial components of our own Galaxy. When more accurate actions are required, this procedure can be combined with torus mapping to produce a fast convergent scheme for action estimation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Actions, angles and frequencies for numerically integrated orbits

We present a method for extracting actions, angles and frequencies from an orbit's time series. The method recovers the generating function that maps an analytic phase-space torus to the torus to which the orbit is confined by simultaneously solving the constraints provided by each time step. We test the method by recovering the actions and frequencies of a triaxial Stäckel potential, and use it to investigate the structure of orbits in a triaxial potential that has been fitted to our Galaxy's Sagittarius stream. The method promises to be useful for analysing N-body simulations. It also takes a step towards constructing distribution functions for the triaxial components of our Galaxy, such as the bar and dark halo.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Probabilistic model for constraining the Galactic potential using tidal streams

We present a generative probabilistic model for a tidal stream and demonstrate how this model is used to constrain the Galactic potential. The model takes advantage of the simple structure of a stream in angle and frequency space for the correct potential. We investigate how the method performs on full 6D mock stream data, and mock data with outliers included. As currently formulated the technique is computationally costly when applied to data with large observational errors, but we describe several modifications that promise to make the technique computationally tractable.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Stream-orbit misalignment I: The dangers of orbit-fitting

Tidal streams don't, in general, delineate orbits. A stream-orbit misalignment is expected to lead to biases when using orbit-fitting to constrain models for the Galactic potential. In this first of two papers we discuss the expected magnitude of the misalignment and the resulting dangers of using orbit-fitting algorithms to constrain the potential. We summarize data for known streams which should prove useful for constraining the Galactic potential, and compute their actions in a realistic Galactic potential. We go on to discuss the formation of tidal streams in angle-action space, and explain why, in general, streams do not delineate orbits. The magnitude of the stream-orbit misalignment is quantified for a logarithmic potential and a multi-component Galactic potential. Specifically, we focus on the expected misalignment for the known streams. By introducing a two-parameter family of realistic Galactic potentials we demonstrate that assuming these streams delineate orbits can lead to order one errors in the halo flattening and halo-to-disc force ratio at the Sun. We present a discussion of the dependence of these results on the progenitor mass, and demonstrate that the misalignment is mass-independent for the range of masses of observed streams. Hence, orbit-fitting does not yield better constraints on the potential if one uses narrower, lower-mass streams.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Stream-orbit misalignment II: A new algorithm to constrain the Galactic potential

In the first of these two papers we demonstrated that assuming streams delineate orbits can lead to order one errors in potential parameters for realistic Galactic potentials. Motivated by the need for an improvement on orbit-fitting, we now present an algorithm for constraining the Galactic potential using tidal streams without assuming that streams delineate orbits. This approach is independent of the progenitor mass so is valid for all observed tidal streams. The method makes heavy use of angle-action variables and seeks the potential which recovers the expected correlations in angle space. We demonstrate that the method can correctly recover the parameters of a simple two-parameter logarithmic potential by analysing an N-body simulation of a stream. We investigate the magnitude of the errors in observational data for which the method can still recover the correct potential and compare this to current and future errors in data. The errors in the observables of individual stars for current and near future data are shown to be too large for the direct use of this method, but when the data are averaged in bins on the sky, the resulting averaged data are accurate enough to constrain correctly the potential parameters for achievable observational errors. From pseudo-data with errors comparable to those that will be furnished in the era of Gaia (20 per cent distance errors, 1.2 mas/yr proper motion errors, and 10 km/s line-of-sight velocity errors) we recover the circular velocity, V_c=220 km/s, and the flattening of the potential, q=0.9, to be V_c=223+/-10km/s and q=0.91+/-0.09.

0 citations0 reviewsNo code linkedOpen access