BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

David W. Hogg

David W. Hogg contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.GAastro-ph.IMastro-ph.SRastro-ph.COApplicationsastro-ph.EPComputationhep-phMachine Learningphysics.data-anphysics.pop-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
21works
0followers
11topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

21 published item(s)

preprint2022arXiv

An empirical model of the Gaia DR3 selection function

Interpreting and modelling astronomical catalogues requires an understanding of the catalogues' completeness or selection function: objects of what properties had a chance to end up in the catalogue. Here we set out to empirically quantify the completeness of the overall Gaia DR3 catalogue. This task is not straightforward because Gaia is the all-sky optical survey with the highest angular resolution to date and no consistent ``ground truth'' exists to allow direct comparisons. However, well-characterised deeper imaging enables an empirical assessment of Gaia's $G$-band completeness across parts of the sky. On this basis, we devised a simple analytical completeness model of Gaia as a function of the observed $G$ magnitude and position over the sky, which accounts for both the effects of crowding and the complex Gaia scanning law. Our model only depends on a single quantity: the median magnitude $M_{10}$ in a patch of the sky of catalogued sources with $\texttt{astrometric_matched_transits}$ $\leq 10$. $M_{10}$ reflects elementary completeness decisions in the Gaia pipeline and is computable from the Gaia DR3 catalogue itself and therefore applicable across the whole sky. We calibrate our model using the Dark Energy Camera Plane Survey (DECaPS) and test its predictions against Hubble Space Telescope observations of globular clusters. We find that our model predicts Gaia's completeness values to a few per cent across the sky. We make the model available as a part of the $\texttt{gaiasf}$ Python package built and maintained by the GaiaUnlimited project: $\texttt{https://github.com/gaia-unlimited/gaiaunlimited}$

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

GREAT3 results I: systematic errors in shear estimation and the impact of real galaxy morphology

We present first results from the third GRavitational lEnsing Accuracy Testing (GREAT3) challenge, the third in a sequence of challenges for testing methods of inferring weak gravitational lensing shear distortions from simulated galaxy images. GREAT3 was divided into experiments to test three specific questions, and included simulated space- and ground-based data with constant or cosmologically-varying shear fields. The simplest (control) experiment included parametric galaxies with a realistic distribution of signal-to-noise, size, and ellipticity, and a complex point spread function (PSF). The other experiments tested the additional impact of realistic galaxy morphology, multiple exposure imaging, and the uncertainty about a spatially-varying PSF; the last two questions will be explored in Paper II. The 24 participating teams competed to estimate lensing shears to within systematic error tolerances for upcoming Stage-IV dark energy surveys, making 1525 submissions overall. GREAT3 saw considerable variety and innovation in the types of methods applied. Several teams now meet or exceed the targets in many of the tests conducted (to within the statistical errors). We conclude that the presence of realistic galaxy morphology in simulations changes shear calibration biases by $\sim 1$ per cent for a wide range of methods. Other effects such as truncation biases due to finite galaxy postage stamps, and the impact of galaxy type as measured by the Sérsic index, are quantified for the first time. Our results generalize previous studies regarding sensitivities to galaxy size and signal-to-noise, and to PSF properties such as seeing and defocus. Almost all methods' results support the simple model in which additive shear biases depend linearly on PSF ellipticity.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Magnitudes, distance moduli, bolometric corrections, and so much more

This pedagogical document about stellar photometry - aimed at those for whom astronomical arcana seem arcane - endeavours to explain the concepts of magnitudes, color indices, absolute magnitudes, distance moduli, extinctions, attenuations, color excesses, K corrections, and bolometric corrections. I include some discussion of observational technique, and some discussion of epistemology, but the primary focus here is on the theoretical or interpretive connections between the observational astronomical quantities and the physical properties of the observational targets.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Mapping Interstellar Dust with Gaussian Processes

Interstellar dust corrupts nearly every stellar observation, and accounting for it is crucial to measuring physical properties of stars. We model the dust distribution as a spatially varying latent field with a Gaussian process (GP) and develop a likelihood model and inference method that scales to millions of astronomical observations. Modeling interstellar dust is complicated by two factors. The first is integrated observations. The data come from a vantage point on Earth and each observation is an integral of the unobserved function along our line of sight, resulting in a complex likelihood and a more difficult inference problem than in classical GP inference. The second complication is scale; stellar catalogs have millions of observations. To address these challenges we develop ziggy, a scalable approach to GP inference with integrated observations based on stochastic variational inference. We study ziggy on synthetic data and the Ananke dataset, a high-fidelity mechanistic model of the Milky Way with millions of stars. ziggy reliably infers the spatial dust map with well-calibrated posterior uncertainties.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Stellar Abundance Maps of the Milky Way Disk

To understand the formation of the Milky Way&#39;s prominent bar it is important to know whether stars in the bar differ in the chemical element composition of their birth material as compared to disk stars. This requires stellar abundance measurements for large samples across the Milky Way&#39;s body. Such samples, e.g. luminous red giant stars observed by SDSS&#39;s Apogee survey, will inevitably span a range of stellar parameters; as a consequence, both modelling imperfections and stellar evolution may preclude consistent and precise estimates of their chemical composition at a level of purported bar signatures, which has left current analyses of a chemically distinct bar inconclusive. Here, we develop a new self-calibration approach to eliminate both modelling and astrophysical abundance systematics among red giant branch (RGB) stars of different luminosities (and hence surface gravity $\log g$). We apply our method to $48,853$ luminous Apogee DR16 RGB stars to construct spatial abundance maps of $20$ chemical elements near the Milky Way&#39;s mid-plane, covering Galactocentric radii of $0\,{\rm kpc}<R_{\rm GC}<20\,\rm kpc$. Our results indicate that there are no abundance variations whose geometry matches that of the bar, and that the mean abundance gradients vary smoothly and monotonically with Galactocentric radius. We confirm that the high-$α$ disk is chemically homogeneous, without spatial gradients. Furthermore, we present the most precise [Fe/H] vs. $R_{\rm GC}$ gradient to date with a slope of $-0.057\pm0.001\rm~dex\,kpc^{-1}$ out to approximately $15$ kpc.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The EXPRES Stellar Signals Project II. State of the Field in Disentangling Photospheric Velocities

Measured spectral shifts due to intrinsic stellar variability (e.g., pulsations, granulation) and activity (e.g., spots, plages) are the largest source of error for extreme precision radial velocity (EPRV) exoplanet detection. Several methods are designed to disentangle stellar signals from true center-of-mass shifts due to planets. The EXPRES Stellar Signals Project (ESSP) presents a self-consistent comparison of 22 different methods tested on the same extreme-precision spectroscopic data from EXPRES. Methods derived new activity indicators, constructed models for mapping an indicator to the needed RV correction, or separated out shape- and shift-driven RV components. Since no ground truth is known when using real data, relative method performance is assessed using the total and nightly scatter of returned RVs and agreement between the results of different methods. Nearly all submitted methods return a lower RV RMS than classic linear decorrelation, but no method is yet consistently reducing the RV RMS to sub-meter-per-second levels. There is a concerning lack of agreement between the RVs returned by different methods. These results suggest that continued progress in this field necessitates increased interpretability of methods, high-cadence data to capture stellar signals at all timescales, and continued tests like the ESSP using consistent data sets with more advanced metrics for method performance. Future comparisons should make use of various well-characterized data sets -- such as solar data or data with known injected planetary and/or stellar signals -- to better understand method performance and whether planetary signals are preserved.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The Thresher: Lucky Imaging without the Waste

In traditional lucky imaging (TLI), many consecutive images of the same scene are taken with a high frame-rate camera, and all but the sharpest images are discarded before constructing the final shift-and-add image. Here we present an alternative image analysis pipeline -- The Thresher -- for these kinds of data, based on online multi-frame blind deconvolution. It makes use of all available data to obtain a best estimate of the astronomical scene in the context of reasonable computational limits; it does not require prior estimates of the point-spread functions in the images, or knowledge of point sources in the scene that could provide such estimates. Most importantly, the scene it aims to return is the optimum of a justified scalar objective based on the likelihood function. Because it uses the full set of images in the stack, The Thresher outperforms TLI in signal-to-noise; as it accounts for the individual-frame PSFs, it does this without loss of angular resolution. We demonstrate the effectiveness of our algorithm on both simulated data and real Electron-Multiplying CCD images obtained at the Danish 1.54m telescope (hosted by ESO, La Silla). We also explore the current limitations of the algorithm, and find that for the choice of image model presented here, non-linearities in flux are introduced into the returned scene. Ongoing development of the software can be viewed at https://github.com/jah1994/TheThresher.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The unpopular Package: a Data-driven Approach to De-trend TESS Full Frame Image Light Curves

The majority of observed pixels on the Transiting Exoplanet Survey Satellite (TESS) are delivered in the form of full frame images (FFI). However, the FFIs contain systematic effects such as pointing jitter and scattered light from the Earth and Moon that must be removed before downstream analysis. We present unpopular, an open-source Python package to de-trend TESS FFI light curves based on the causal pixel model method. Under the assumption that shared flux variations across multiple distant pixels are likely to be systematics, unpopular removes these common (i.e., popular) trends by modeling the systematics in a given pixel&#39;s light curve as a linear combination of light curves from many other distant pixels. To prevent overfitting we employ ridge regression and a train-and-test framework where the data points being de-trended are separated from those used to obtain the model coefficients. We also allow for simultaneous fitting with a polynomial model to capture any long-term astrophysical trends. We validate our method by de-trending different sources (e.g., supernova, tidal disruption event, exoplanet-hosting star, fast rotating star) and comparing our light curves to those obtained by other pipelines when appropriate. We also show that unpopular is able to preserve sector-length astrophysical signals, allowing for the extraction of multi-sector light curves from the FFI data. The unpopular source code and tutorials are freely available online.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

An unsupervised method for identifying $X$-enriched stars directly from spectra: Li in LAMOST

Stars with peculiar element abundances are important markers of chemical enrichment mechanisms. We present a simple method, tangent space projection (TSP), for the detection of $X$-enriched stars, for arbitrary elements $X$, even from blended lines. Our method does not require stellar labels, but instead directly estimates the counterfactual unrenriched spectrum from other unlabelled spectra. As a case study, we apply this method to the $6708~$Å Li doublet in LAMOST DR5, identifying 8,428 Li-enriched stars seamlessly across evolutionary state. We comment on the explanation for Li-enrichement for different subpopulations, including planet accretion, nonstandard mixing, and youth.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Dimensionality reduction, regularization, and generalization in overparameterized regressions

Overparameterization in deep learning is powerful: Very large models fit the training data perfectly and yet often generalize well. This realization brought back the study of linear models for regression, including ordinary least squares (OLS), which, like deep learning, shows a &#34;double-descent&#34; behavior: (1) The risk (expected out-of-sample prediction error) can grow arbitrarily when the number of parameters $p$ approaches the number of samples $n$, and (2) the risk decreases with $p$ for $p>n$, sometimes achieving a lower value than the lowest risk for $p<n$. The divergence of the risk for OLS can be avoided with regularization. In this work, we show that for some data models it can also be avoided with a PCA-based dimensionality reduction (PCA-OLS, also known as principal component regression). We provide non-asymptotic bounds for the risk of PCA-OLS by considering the alignments of the population and empirical principal components. We show that dimensionality reduction improves robustness while OLS is arbitrarily susceptible to adversarial attacks, particularly in the overparameterized regime. We compare PCA-OLS theoretically and empirically with a wide range of projection-based methods, including random projections, partial least squares (PLS), and certain classes of linear two-layer neural networks. These comparisons are made for different data generation models to assess the sensitivity to signal-to-noise and the alignment of regression coefficients with the features. We find that methods in which the projection depends on the training data can outperform methods where the projections are chosen independently of the training data, even those with oracle knowledge of population quantities, another seemingly paradoxical phenomenon that has been identified previously. This suggests that overparameterization may not be necessary for good generalization.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Snails Across Scales: Local and Global Phase-Mixing Structures as Probes of the Past and Future Milky Way

Signatures of vertical disequilibrium have been observed across the Milky Way&#39;s disk. These signatures manifest locally as unmixed phase-spirals in $z$--$v_z$ space (&#34;snails-in-phase&#34;) and globally as nonzero mean $z$ and $v_z$ which wraps around as a physical spiral across the $x$--$y$ plane (&#34;snails-in-space&#34;). We explore the connection between these local and global spirals through the example of a satellite perturbing a test-particle Milky Way (MW)-like disk. We anticipate our results to broadly apply to any vertical perturbation. Using a $z$--$v_z$ asymmetry metric we demonstrate that in test-particle simulations: (a) multiple local phase-spiral morphologies appear when stars are binned by azimuthal action $J_ϕ$, 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&#39;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 within $R \lesssim 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} \approx -6$ kpc, $v_{z,Sgr} \approx 210$ km/s), we demonstrate that the peak of the impact has already passed. Sgr&#39;s pull over the past 150 Myr creates a global $v_z$ signature with amplitude $\propto M_{Sgr}$, which might be detectable in future spectroscopic surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Close Binary Companions to APOGEE DR16 Stars: 20,000 Binary-star Systems Across the Color-Magnitude Diagram

Many problems in contemporary astrophysics---from understanding the formation of black holes to untangling the chemical evolution of galaxies---rely on knowledge about binary stars. This, in turn, depends on discovery and characterization of binary companions for large numbers of different kinds of stars in different chemical and dynamical environments. Current stellar spectroscopic surveys observe hundreds of thousands to millions of stars with (typically) few observational epochs, which allows binary discovery but makes orbital characterization challenging. We use a custom Monte Carlo sampler (The Joker) to perform discovery and characterization of binary systems through radial-velocities, in the regime of sparse, noisy, and poorly sampled multi-epoch data. We use it to generate posterior samplings in Keplerian parameters for 232,531 sources released in APOGEE Data Release 16. Our final catalog contains 19,635 high-confidence close-binary (P < few years, a < few AU) systems that show interesting relationships between binary occurrence rate and location in the color-magnitude diagram. We find notable faint companions at high masses (black-hole candidates), at low masses (substellar candidates), and at very close separations (mass-transfer candidates). We also use the posterior samplings in a (toy) hierarchical inference to measure the long-period binary-star eccentricity distribution. We release the full set of posterior samplings for the entire parent sample of 232,531 stars. This set of samplings involves no heuristic &#34;discovery&#34; threshold and therefore can be used for myriad statistical purposes, including hierarchical inferences about binary-star populations and sub-threshold searches.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Data Analysis Recipes: Products of multivariate Gaussians in Bayesian inferences

A product of two Gaussians (or normal distributions) is another Gaussian. That&#39;s a valuable and useful fact! Here we use it to derive a refactoring of a common product of multivariate Gaussians: The product of a Gaussian likelihood times a Gaussian prior, where some or all of those parameters enter the likelihood only in the mean and only linearly. That is, a linear, Gaussian, Bayesian model. This product of a likelihood times a prior pdf can be refactored into a product of a marginalized likelihood (or a Bayesian evidence) times a posterior pdf, where (in this case) both of these are also Gaussian. The means and variance tensors of the refactored Gaussians are straightforward to obtain as closed-form expressions; here we deliver these expressions, with discussion. The closed-form expressions can be used to speed up and improve the precision of inferences that contain linear parameters with Gaussian priors. We connect these methods to inferences that arise frequently in physics and astronomy. If all you want is the answer, the question is posed and answered at the beginning of Section 3. We show two toy examples, in the form of worked exercises, in Section 4. The solutions, discussion, and exercises in this Note are aimed at someone who is already familiar with the basic ideas of Bayesian inference and probability.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Forward modeling the orbits of companions to pulsating stars from their light travel time variations

Mutual gravitation between a pulsating star and an orbital companion leads to a time-dependent variation in path length for starlight traveling to Earth. These variations can be used for coherently pulsating stars, such as the δ Scuti variables, to constrain the masses and orbits of their companions. Observing these variations for δ Scuti stars has previously relied on subdividing the light curve and measuring the average pulsation phase in equally sized subdivisions, which leads to under-sampling near periapsis. We introduce a new approach that simultaneously forward-models each sample in the light curve and show that this method improves upon current sensitivity limits - especially in the case of highly eccentric and short-period binaries. We find that this approach is sensitive enough to observe Jupiter mass planets around δ Scuti stars under ideal conditions, and use gravity-mode pulsations in the subdwarf B star KIC 7668647 to detect its companion without radial velocity data. We further provide robust detection limits as a function of the SNR of the pulsation mode and determine that the minimum detectable light travel time amplitude for a typical Kepler δ Scuti is around 2 s. This new method significantly enhances the application of light travel time variations to detecting short period binaries with pulsating components, and pulsating A-type exoplanet host stars, especially as a tool for eliminating false positives.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

High-resolution spectroscopy of the GD-1 stellar stream localizes the perturber near the orbital plane of Sagittarius

The $100^\circ$-long thin stellar stream in the Milky Way halo, GD-1, has an ensemble of features that may be due to dynamical interactions. Using high-resolution MMT/Hectochelle spectroscopy we show that a spur of GD-1-like stars outside of the main stream are kinematically and chemically consistent with the main stream. In the spur, as in the main stream, GD-1 has a low intrinsic radial velocity dispersion, $σ_{V_r}\lesssim1\,\rm km\,s^{-1}$, is metal-poor, $\rm [Fe/H]\approx-2.3$, with little $\rm [Fe/H]$ spread and some variation in $\rm [α/Fe]$ abundances, which point to a common globular cluster progenitor. At a fixed location along the stream, the median radial velocity offset between the spur and the main stream is smaller than $0.5\,\rm km\,s^{-1}$, comparable to the measurement uncertainty. A flyby of a massive, compact object can change orbits of stars in a stellar stream and produce features like the spur observed in GD-1. In this scenario, the radial velocity of the GD-1 spur relative to the stream constrains the orbit of the perturber and its current on-sky position to $\approx5,000\,\rm deg^2$. The family of acceptable perturber orbits overlaps the stellar and dark-matter debris of the Sagittarius dwarf galaxy in present-day position and velocity. This suggests that GD-1 may have been perturbed by a globular cluster or an extremely compact dark-matter subhalo formerly associated with Sagittarius.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

How to obtain the redshift distribution from probabilistic redshift estimates

A trustworthy estimate of the redshift distribution $n(z)$ is crucial for using weak gravitational lensing and large-scale structure of galaxy catalogs to study cosmology. Spectroscopic redshifts for the dim and numerous galaxies of next-generation weak-lensing surveys are expected to be unavailable, making photometric redshift (photo-$z$) probability density functions (PDFs) the next-best alternative for comprehensively encapsulating the nontrivial systematics affecting photo-$z$ point estimation. The established stacked estimator of $n(z)$ avoids reducing photo-$z$ PDFs to point estimates but yields a systematically biased estimate of $n(z)$ that worsens with decreasing signal-to-noise, the very regime where photo-$z$ PDFs are most necessary. We introduce Cosmological Hierarchical Inference with Probabilistic Photometric Redshifts (CHIPPR), a statistically rigorous probabilistic graphical model of redshift-dependent photometry, which correctly propagates the redshift uncertainty information beyond the best-fit estimator of $n(z)$ produced by traditional procedures and is provably the only self-consistent way to recover $n(z)$ from photo-$z$ PDFs. We present the $\texttt{chippr}$ prototype code, noting that the mathematically justifiable approach incurs computational expense. The CHIPPR approach is applicable to any one-point statistic of any random variable, provided the prior probability density used to produce the posteriors is explicitly known; if the prior is implicit, as may be the case for popular photo-$z$ techniques, then the resulting posterior PDFs cannot be used for scientific inference. We therefore recommend that the photo-$z$ community focus on developing methodologies that enable the recovery of photo-$z$ likelihoods with support over all redshifts, either directly or via a known prior probability density.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Principled point-source detection in collections of astronomical images

We review the well-known matched filter method for the detection of point sources in astronomical images. This is shown to be optimal (that is, to saturate the Cramer--Rao bound) under stated conditions that are very strong: an isolated source in background-dominated imaging with perfectly known background level, point-spread function, and noise models. We show that the matched filter produces a maximum-likelihood estimate of the brightness of a purported point source, and this leads to a simple way to combine multiple images---taken through the same bandpass filter but with different noise levels and point-spread functions---to produce an optimal point source detection map. We then extend the approach to images taken through different bandpass filters, introducing the SED-matched filter, which allows us to combine images taken through different filters, but requires us to specify the colors of the objects we wish to detect. We show that this approach is superior to some methods traditionally employed, and that other traditional methods can be seen as instances of SED-matched filtering with implied (and often unreasonable) priors. We present a Bayesian formulation, including a flux prior that leads to a closed-form expression with low computational cost.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Temperatures and Metallicities of M dwarfs in the APOGEE Survey

M dwarfs have enormous potential for our understanding of structure and formation on both Galactic and exoplanetary scales through their properties and compositions. However, current atmosphere models have limited ability to reproduce spectral features in stars at the coolest temperatures ($T_{\rm eff} < 4200\,$K) and to fully exploit the information content of current and upcoming large-scale spectroscopic surveys. Here we present a catalog of spectroscopic temperatures, metallicities, and spectral types for 5875 M dwarfs in the Apache Point Observatory Galactic Evolution Experiment (APOGEE) and Gaia DR2 surveys using The Cannon: a flexible, data-driven spectral-modeling and parameter-inference framework demonstrated to estimate stellar-parameter labels ($T_{\rm eff}$, log$g$, [Fe/H], and detailed abundances) to high precision. Using a training sample of 87 M dwarfs with optically derived labels spanning $2860 < T_{\rm eff} < 4130\,$K calibrated with bolometric temperatures, and $-0.5 < $[Fe/H]$ < 0.5\,$dex calibrated with FGK binary metallicities, we train a two-parameter model with predictive accuracy (in cross-validation) to $77\,$K and $0.09\,$dex respectively. We also train a one-dimensional spectral classification model using 51 M dwarfs with Sloan Digital Sky Survey optical spectral types ranging from M0 to M6, to predictive accuracy of 0.7 types. We find Cannon temperatures to be in agreement to within $60\,$K compared to a subsample of 1702 sources with color-derived temperatures, and Cannon metallicities to be in agreement to within $0.08\,$dex metallicity compared to a subsample of 15 FGK+M or M+M binaries. Finally, our comparison between Cannon and APOGEE pipeline (ASPCAP DR14) labels finds that ASPCAP is systematically biased toward reporting higher temperatures and lower metallicities for M dwarfs.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The power of co-ordinate transformations in dynamical interpretations of Galactic structure

$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 \citep{Khoperskov+19b}. However, in this work, we show that these features are in fact the known kinematic moving groups, both in the $L_z-ϕ$ and the $v_{\mathrm{R}}-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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Sixteenth Data Release of the Sloan Digital Sky Surveys: First Release from the APOGEE-2 Southern Survey and Full Release of eBOSS Spectra

This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the main cosmological program of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and all raw and reduced spectra from that project are released here. DR16 also includes all the data from the Time Domain Spectroscopic Survey (TDSS) and new data from the SPectroscopic IDentification of ERosita Survey (SPIDERS) programs, both of which were co-observed on eBOSS plates. DR16 has no new data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (or the MaNGA Stellar Library &#34;MaStar&#34;). We also preview future SDSS-V operations (due to start in 2020), and summarize plans for the final SDSS-IV data release (DR17).

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Strength of the Dynamical Spiral Perturbation in the Galactic Disk

The mean Galactocentric radial velocities $\langle v_{R}\rangle(R,φ)$ of luminous red giant stars within the mid-plane of the Milky Way reveal a spiral signature, which could plausibly reflect the response to a non-axisymmetric perturbation of the gravitational potential in the Galactic disk. We apply a simple steady-state toy model of a logarithmic spiral to interpret these observations, and find a good qualitative and quantitative match. Presuming that the amplitude of the gravitational potential perturbation is proportionate to that in the disk&#39;s surface mass density, we estimate the surface mass density amplitude to be $Σ_{\rm max} (R_{\odot})\approx 5.5\,\rm M_{\odot}\,pc^{-2}$ at the solar radius when choosing a fixed pattern speed of $Ω_{\mathrm p}=12\,\rm km\,s^{-1}\,kpc^{-1}$. Combined with the local disk density, this implies a surface mass density contrast between the arm and inter-arm regions of approximately $\pm 10\%$ at the solar radius, with an increases towards larger radii. Our model constrains the pitch angle of the dynamical spiral arms to be approximately $12^{\circ}$.

0 citations0 reviewsNo code linkedOpen access