Source author record

Sergey E. Koposov

Sergey E. Koposov appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

45works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

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

45 published item(s)

preprint2026arXiv

DESI Data Release 1: Stellar Catalogue

In this paper we present the stellar Value-Added Catalogue (VAC) based on the DESI Data Release 1. This VAC contains stellar parameter, abundance and radial velocity measurements for more than 4 million stars. It also contains, for the first time, measurements from individual epochs for more than a million stars with at least two observations. The main contribution to the catalogue comes from the bright program of the main survey, which includes $\sim $2.5 million stars, and the backup program, which includes $\sim $ 1 million stars. The combined magnitude range for the stars in the catalogue extends from ${\it Gaia}$ G $\sim 12$ to G $\sim 21$. For the magnitude range $17.5<G<21$ this catalogue represents a factor of 10 increase in the number of stars with radial velocity and abundance measurements compared to existing surveys. Despite DESI's resolution (R $\sim 2500-5000$), the median radial velocity uncertainty for stars in the catalogue is better than 1 km s$^{-1}$. The stellar parameters and abundances of stars in DESI are measured by two independent pipelines, and after applying a temperature-dependent calibration, [Fe/H] abundances of high signal-to-noise stars are accurate to better than $\sim$ 0.1 dex when compared to high-resolution surveys. The catalogue probes different Galactic components including a particularly large number of distant stars: tens of thousands of stars further than 10 kpc, and thousands further than 50 kpc. The catalogue also contains several thousand extremely metal-poor stars with ${\rm [Fe/H]}<-3$. The released sample of stars includes measurements for thousands of stars that are members of dwarf galaxies, open and globular clusters as well as members of several dozen stellar streams. The next public DESI data release is expected in less than two years and will contain three times as many stars as DR1.

preprint2026arXiv

Multiple Populations and a CH Star Found in the 300S Globular Cluster Stellar Stream

Milky Way globular clusters (GCs) display chemical enrichment in a phenomenon called multiple stellar populations (MSPs). While the enrichment mechanism is not fully understood, there is a correlation between a cluster's mass and the fraction of enriched stars found therein. However, present-day GC masses are often smaller than their masses at the time of formation due to dynamical mass loss. In this work, we explore the relationship between mass and MSPs using the stellar stream 300S. We present the chemical abundances of eight red giant branch member stars in 300S with high-resolution spectroscopy from Magellan/MIKE. We identify one enriched star characteristic of MSPs and no detectable metallicity dispersion, confirming that the progenitor of 300S was a globular cluster. The fraction of enriched stars (12.5\%) observed in our 300S stars is less than the 50\% of stars found enriched in Milky Way GCs of comparable present-day mass ($\sim10^{4.5}$\msun). We calculate the mass of 300S's progenitor and compare it to the initial masses of intact GCs, finding that 300S aligns well with the trend between the system mass at formation and enrichment. 300S's progenitor may straddle the critical mass threshold for the formation of MSPs and can therefore serve as a benchmark for the stellar enrichment process. Additionally, we identify a CH star, with high abundances of \textit{s}-process elements, probably accreted from a binary companion. The rarity of such binaries in intact GCs may imply stellar streams permit the survival of binaries that would otherwise be disrupted.

preprint2023arXiv

The Target-selection Pipeline for the Dark Energy Spectroscopic Instrument

In 2021 May, the Dark Energy Spectroscopic Instrument (DESI) began a 5 yr survey of approximately 50 million total extragalactic and Galactic targets. The primary DESI dark-time targets are emission line galaxies (ELGs), luminous red galaxies (LRGs) and quasars (QSOs). In bright time, DESI will focus on two surveys known as the Bright Galaxy Survey (BGS) and the Milky Way Survey (MWS). DESI also observes a selection of "secondary" targets for bespoke science goals. This paper gives an overview of the publicly available pipeline (desitarget) used to process targets for DESI observations. Highlights include details of the different DESI survey targeting phases, the targeting ID (TARGETID) used to define unique targets, the bitmasks used to indicate a particular type of target, the data model and structure of DESI targeting files, and examples of how to access and use the desitarget code base. This paper will also describe "supporting" DESI target classes, such as standard stars, sky locations, and random catalogs that mimic the angular selection function of DESI targets. The DESI target selection pipeline is complex and sizable; this paper attempts to summarize the most salient information required to understand and work with DESI targeting data.

preprint2022arXiv

$S^5$: The Orbital and Chemical Properties of One Dozen Stellar Streams

We report the kinematic, orbital, and chemical properties of 12 stellar streams with no evident progenitors, using line-of-sight velocities and metallicities from the Southern Stellar Stream Spectroscopic Survey ($S^5$), proper motions from $Gaia$ EDR3, and distances derived from distance tracers or the literature. This data set provides the largest homogeneously analyzed set of streams with full 6D kinematics and metallicities. All streams have heliocentric distances between ${\sim}10-50$ kpc. The velocity and metallicity dispersions show that half of the stream progenitors were disrupted dwarf galaxies (DGs), while the other half originated from disrupted globular clusters (GCs), hereafter referred to as DG and GC streams. Based on the mean metallicities of the streams and the mass-metallicity relation, the luminosities of the progenitors of the DG streams range between Carina and Ursa Major I ($-9.5\lesssim M_V\lesssim-5.5$). Four of the six GC streams have mean metallicities of [Fe/H]$< -2$, more metal-poor than typical Milky Way (MW) GCs at similar distances. Interestingly, the 300S and Jet GC streams are the only streams on retrograde orbits in our dozen stream sample. Finally, we compare the orbital properties of the streams with known DGs and GCs in the MW, finding several possible associations. Some streams appear to have been accreted with the recently discovered Gaia-Enceladus-Sausage system, and others suggest that GCs were formed in and accreted together with the progenitors of DG streams whose stellar masses are similar to Draco to Carina ($\sim10^5-10^6M_\odot$).

preprint2022arXiv

Catalogues of Active Galactic Nuclei From Gaia and unWISE Data

We present two catalogues of active galactic nucleus (AGN) candidates selected from the latest data of two all-sky surveys -- Data Release 2 (DR2) of the \emph{Gaia} mission and the unWISE catalogue of the \emph{Wide-field Infrared Survey Explorer} (\emph{WISE}). We train a random forest classifier to predict the probability of each source in the \emph{Gaia}-unWISE joint sample being an AGN, $P_{\rm RF}$, based on \emph{Gaia} astrometric and photometric measurements and unWISE photometry. The two catalogues, which we designate C75 and R85, are constructed by applying different $P_{\rm RF}$ threshold cuts to achieve an overall completeness of 75\% ($\approx$90\% at \emph{Gaia} $G\leq20$ mag) and reliability of 85\% respectively. The C75 (R85) catalogue contains 2,734,464 (2,182,193) AGN candidates across the effective 36,000 deg$^2$ sky, of which $\approx$0.91 (0.52) million are new discoveries. Photometric redshifts of the AGN candidates are derived by a random forest regressor using \emph{Gaia} and \emph{WISE} magnitudes and colours. The estimated overall photometric redshift accuracy is 0.11. Cross-matching the AGN candidates with a sample of known bright cluster galaxies, we identify a high-probability strongly-lensed AGN candidate system, SDSS\,J1326$+$4806, with a large image separation of 21\farcs06. All the AGN candidates in our catalogues will have $\sim$5-year long light curves from \emph{Gaia} by the end of the mission, and thus will be a great resource for AGN variability studies. Our AGN catalogues will also be helpful in AGN target selections for future spectroscopic surveys, especially ones in the southern hemisphere. The C75 catalogue can be downloaded at https://zenodo.org/record/6837642#.YtEgjy8RrzI.

preprint2022arXiv

Constraining the shape of Milky Way satellites with distance gradients

We combine the Dark Energy Camera Legacy Survey (DECaLS) DR8 photometry with Gaia photometry to study the 3-D structure of Bootes I, Draco, Ursa Minor, Sextans and Sculptor dwarf galaxies using blue horizontal branch (BHB) stars as distance indicators. We construct a new colour-absolute magnitude of BHB stars that we use to measure the distance gradients within the body of the dwarf galaxies. We detect a statistically significant non-zero gradient only in Sextans and Sculptor. Through modeling of the gradient and 2-D density of the systems by triaxial Plummer models we find that the distance gradients in both dwarf galaxies are inconsistent with prolate shape, but compatible with oblate or triaxial shapes. In order to explain the observed gradients, oblate models of Sextans and Sculptor need to have a significant intrinsic ellipticity larger than $0.47$ for Sextans and $0.46$ for Sculptor. The flattened oblate shape may imply a significant anisotropy in velocity distribution in order to be consistent with the lack of significant velocity gradients in these systems.

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.

preprint2022arXiv

Snowmass2021 Cosmic Frontier White Paper: Prospects for obtaining Dark Matter Constraints with DESI

Despite efforts over several decades, direct-detection experiments have not yet led to the discovery of the dark matter (DM) particle. This has led to increasing interest in alternatives to the Lambda CDM (LCDM) paradigm and alternative DM scenarios (including fuzzy DM, warm DM, self-interacting DM, etc.). In many of these scenarios, DM particles cannot be detected directly and constraints on their properties can ONLY be arrived at using astrophysical observations. The Dark Energy Spectroscopic Instrument (DESI) is currently one of the most powerful instruments for wide-field surveys. The synergy of DESI with ESA's Gaia satellite and future observing facilities will yield datasets of unprecedented size and coverage that will enable constraints on DM over a wide range of physical and mass scales and across redshifts. DESI will obtain spectra of the Lyman-alpha forest out to z~5 by detecting about 1 million QSO spectra that will put constraints on clustering of the low-density intergalactic gas and DM halos at high redshift. DESI will obtain radial velocities of 10 million stars in the Milky Way (MW) and Local Group satellites enabling us to constrain their global DM distributions, as well as the DM distribution on smaller scales. The paradigm of cosmological structure formation has been extensively tested with simulations. However, the majority of simulations to date have focused on collisionless CDM. Simulations with alternatives to CDM have recently been gaining ground but are still in their infancy. While there are numerous publicly available large-box and zoom-in simulations in the LCDM framework, there are no comparable publicly available WDM, SIDM, FDM simulations. DOE support for a public simulation suite will enable a more cohesive community effort to compare observations from DESI (and other surveys) with numerical predictions and will greatly impact DM science.

preprint2022arXiv

Stellar kinematics of dwarf galaxies from multi-epoch spectroscopy: application to Triangulum II

We present new MMT/Hectochelle spectroscopic measurements for 257 stars observed along the line of sight to the ultra-faint dwarf galaxy Triangulum II. Combining with results from previous Keck/DEIMOS spectroscopy, we obtain a sample that includes 16 likely members of Triangulum II, with up to 10 independent redshift measurements per star. To this multi-epoch kinematic data set we apply methodology that we develop in order to infer binary orbital parameters from sparsely sampled radial velocity curves with as few as two epochs. For a previously-identified (spatially unresolved) binary system in Tri~II, we infer an orbital solution with period $296.0_{-3.3}^{+3.8} \rm~ days$ , semi-major axis $1.12^{+0.41}_{-0.24}\rm~AU$, and a systemic velocity $ -380.0 \pm 1.7 \rm~km ~s^{-1}$ that we then use in the analysis of Tri~II's internal kinematics. Despite this improvement in the modeling of binary star systems, the current data remain insufficient to resolve the velocity dispersion of Triangulum II. We instead find a 95% confidence upper limit of $σ_{v} \lesssim 3.4 \rm ~km~s^{-1}$.

preprint2022arXiv

Stellar multiplicity and stellar rotation: Insights from APOGEE

We measure rotational broadening in spectra taken by the Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey to characterise the relationship between stellar multiplicity and rotation. We create a sample of 2786 giants and 24 496 dwarfs with stellar parameters and multiple radial velocities from the APOGEE pipeline, projected rotation speeds \vsini\ determined from our own pipeline, and distances, masses, and ages measured by Sanders \& Das. We use the statistical distribution of the maximum shift in the radial velocities, \drvm, as a proxy for the close binary fraction to explore the interplay between stellar evolution, rotation, and multiplicity. Assuming that the minimum orbital period allowed is the critical period for Roche Lobe overflow and rotational synchronization, we calculate theoretical upper limits on expected \vsini\ and \drvm\ values. These expectations agree with the positive correlation between the maximum \drvm\ and \vsini\ values observed in our sample as a function of \logg. We find that the fast rotators in our sample have a high occurrence of short-period ($\log(P/\text{d})\lesssim 4$) companions. We also find that old, rapidly-rotating main sequence stars have larger completeness-corrected close binary fractions than their younger peers. Furthermore, rapidly-rotating stars with large \drvm\ consistently show differences of 1-10 Gyr between the predicted gyrochronological and measured isochronal ages. These results point towards a link between rapid rotation and close binarity through tidal interactions. We conclude that stellar rotation is strongly correlated with stellar multiplicity in the field, and caution should be taken in the application of gyrochronology relations to cool stars.

preprint2022arXiv

Uniform Modelling of the Stellar Density of Thirteen Tidal Streams within the Galactic Halo

We present the results of fitting a flexible stellar stream density model to a collection of thirteen streams around the Milky Way, using photometric data from DES, DECaLS, and Pan-STARRS. We construct density maps for each stream and characterise their tracks on the sky, width, and distance modulus curves along the length of each stream. We use these measurements to compute lengths and total luminosities of streams and identify substructures. Several streams show prominent substructures, such as stream broadening, gaps, large deviations of stream tracks and sharp changes in stream densities. Examining the group of streams as a population, as expected we find that streams with globular cluster progenitors are typically narrower than those with dwarf galaxy progenitors, with streams around 100 pc wide showing overlap between the two populations. We also note the average luminosity of globular cluster streams is significantly lower than the typical luminosity of intact globular clusters. The likely explanation is that observed globular cluster streams preferentially come from lower luminosity and lower density clusters. The stream measurements done in a uniform manner presented here will be helpful for more detailed stream studies such as identifying candidate stream members for spectroscopic follow up and stellar stream dynamical modeling.

preprint2021arXiv

Elevated r-process enrichment in Gaia Sausage and Sequoia

The Gaia Sausage (GS) and the Sequoia represent the major accretion events that formed the stellar halo of the Milky Way. A detailed chemical study of these main building blocks provides a pristine view of the early steps of the Galaxy's assembly. We present the results of the analysis of the UVES high-resolution spectroscopic observations at the 8.2m VLT of 9 Sausage/Sequoia members selected kinematically using Gaia DR2. We season this set of measurements with archival data from Nissen & Schuster (2011) and GALAH DR3 (2020). Here, we focus on the neutron-capture process by analysing Sr, Y, Ba and Eu behavior. We detect clear enhancement in Eu abundance ([Eu/Fe]~0.6-0.7) indicative of large prevalence of r-process in the stellar n-capture makeup. We are also able to trace the evolution of the heavy element production across a wide range of metallicity. The barium to europium changes from a tight, flat sequence with [Ba/Eu]=-0.7 reflecting dominant contribution from exploding massive stars, to a clear upturn at higher iron abundances, betraying the onset of contamination from asymptotic giant branch (AGB) ejecta. Additionally, we discover two clear sequences in [Fe/H]-[Ba/Fe] plane likely caused by distinct levels of s-process pollution and mixing within the GS progenitor.

preprint2021arXiv

Measuring the Mass of the Large Magellanic Cloud with Stellar Streams Observed by ${S}^5$

Stellar streams are excellent probes of the underlying gravitational potential in which they evolve. In this work, we fit dynamical models to five streams in the Southern Galactic hemisphere, combining observations from the Southern Stellar Stream Spectroscopic Survey (${S}^5$), Gaia EDR3, and the Dark Energy Survey (DES), to measure the mass of the Large Magellanic Cloud (LMC). With an ensemble of streams, we find a mass of the LMC ranging from 14 to $19 \times 10^{10}\ \mathrm{M}_{\odot}$, probed over a range of closest approach times and distances. With the most constraining stream (Orphan-Chenab), we measure an LMC mass of $18.8^{+ 3.5}_{- 4.0} \times 10^{10}\ \mathrm{M}_{\odot}$, probed at a closest approach time of 310 Myr and a closest approach distance of 25.4 kpc. This mass is compatible with previous measurements, showing that a consistent picture is emerging of the LMC's influence on structures in the Milky Way. Using this sample of streams, we find that the LMC's effect depends on the relative orientation of the stream and LMC at their point of closest approach. To better understand this, we present a simple model based on the impulse approximation and we show that the LMC's effect depends both on the magnitude of the velocity kick imparted to the stream and the direction of this kick.

preprint2021arXiv

Two-Point Separation Functions for Modeling Wide Binary Systems in Nearby Dwarf Galaxies

We use a geometric method to derive (two-dimensional) separation functions amongst pairs of objects within populations of specified position function $\mathrm{d} N/\mathrm{d} \vec{R}$. We present analytic solutions for separation functions corresponding to a uniform surface density within a circular field, a Plummer sphere (viewed in projection), and the mixture thereof -- including contributions from binary objects within both sub-populations. These results enable inferences about binary object populations via direct modeling of object position and pair separation data, without resorting to standard estimators of the two-point correlation function. Analyzing mock data sets designed to mimic known dwarf spheroidal galaxies, we demonstrate the ability to recover input properties including the number of wide binary star systems and, in cases where the number of resolved binary pairs is assumed to be $\gtrsim$ a few hundred, characteristic features (e.g., steepening and/or truncation) of their separation function. Combined with forthcoming observational capabilities, this methodology opens a window onto the formation and/or survival of wide binary populations in dwarf galaxies, and offers a novel probe of dark matter substructure on the smallest galactic scales.

preprint2020arXiv

A Mystery in Chamaeleon: Serendipitous Discovery of a Galactic Symbiotic Nova

We present the serendipitous discovery of a low luminosity nova occurring in a symbiotic binary star system in the Milky Way. We lay out the extensive archival data alongside new follow-up observations related to the stellar object V$^*$ CN Cha in the constellation of Chamaeleon. The object had long period ($\sim\! 250\,$day), high amplitude ($\sim\! 3\,$mag) optical variability in its recent past, preceding an increase in optical brightness by $\sim\! 8\,$magnitudes and a persistence at this luminosity for about 3 years, followed by a period of $\sim\! 1.4\,{\rm mag}\,{\rm yr}^{-1}$ dimming. The object's current optical luminosity seems to be dominated by H$α$ emission, which also exhibits blue-shifted absorption (a P-Cygni-like profile). After consideration of a number of theories to explain these myriad observations, we determine that V$^*$ CN Cha is most likely a symbiotic (an evolved star-white dwarf binary) system which has undergone a long-duration, low luminosity, nova. Interpreted in this way, the outburst in V$^*$ CN Cha is among the lowest luminosity novae ever observed.

preprint2020arXiv

Discovery of a disrupting open cluster far into the Milky Way halo: a recent star formation event in the leading arm of the Magellanic stream?

We report the discovery of a young (${\rm age} \sim 130~{\rm Myr}$), low-mass ($M \sim 1200~{\rm M}_\odot$), metal-poor ($[{\rm Fe}/{\rm H}] \sim -1.1$) stellar association at a heliocentric distance $D \approx 29~{\rm kpc}$, placing it far into the Milky Way halo. At its present Galactocentric position $(R, z) \sim (23, 15)~{\rm kpc}$, the association is (on the sky) near the leading arm of the gas stream emanating from the Magellanic cloud system, but is located $\approx 60^\circ$ from the Large Magellanic Cloud (LMC) center on the other side of the Milky Way disk. If we assume that the cluster is co-located with HI gas in the stream, we directly measure the distance to the leading arm of the Magellanic stream. The measured distance is inconsistent with stream predictions from models of the LMC/SMC interaction and infall into the Milky Way that do not account for ram pressure and gas interaction with Milky Way disk. The estimated age of the cluster is consistent with the time of last passage of the leading arm gas through the Galactic midplane, and we therefore speculate that this star-formation event was triggered by its the last disk midplane passage. Most details of this idea remain a puzzle: the Magellanic stream has low column density, the Milky Way disk at this large radius has low gas density, and the relative velocity of the leading arm gas and Milky Way gas is large. However it formed, the discovery of a young stellar cluster in the Milky Way halo presents an interesting opportunity for study. This cluster was discovered with Gaia astrometry and photometry alone, but folow-up DECam photometry was crucial for measuring its properties.

preprint2020arXiv

Exploring Chemical Homogeneity in Dwarf Galaxies: A VLT-MUSE study of JKB 18

Deciphering the distribution of metals throughout galaxies is fundamental in our understanding of galaxy evolution. Nearby, low-metallicity, star-forming dwarf galaxies in particular can offer detailed insight into the metal-dependent processes that may have occurred within galaxies in the early Universe. Here we present VLT/MUSE observations of one such system, JKB 18, a blue diffuse dwarf galaxy with a metallicity of only 12+log(O/H)=7.6$\pm$0.2 (~0.08 Zsol). Using high spatial-resolution integral field spectroscopy of the entire system, we calculate chemical abundances for individual HII regions using the direct method and derive oxygen abundance maps using strong-line metallicity diagnostics. With large-scale dispersions in O/H, N/H and N/O of ~0.5-0.6 dex and regions harbouring chemical abundances outside this 1$σ$ distribution, we deem JKB 18 to be chemically inhomogeneous. We explore this finding in the context of other chemically inhomogeneous dwarf galaxies and conclude that neither the accretion of metal-poor gas, short mixing timescales, or self-enrichment from Wolf-Rayet stars are accountable. Using a galaxy-scale, multi-phase, hydrodynamical simulation of a low-mass dwarf galaxy, we find that chemical inhomogeneities of this level may be attributable to the removal of gas via supernovae and the specific timing of the observations with respect to star-formation activity. This study not only draws attention to the fact that dwarf galaxies can be chemically inhomogeneous, but also that the methods used in the assessment of this characteristic can be subject to bias.

preprint2020arXiv

Mapping the tilt of the Milky Way bulge velocity ellipsoids with ARGOS and $Gaia$ DR2

Until the recent advent of $Gaia$ Data Release 2 (DR2) and deep multi-object spectroscopy, it has been difficult to obtain 6-D phase space information for large numbers of stars beyond 4 kpc, in particular towards the Galactic centre, where dust and crowding effects are significant. In this study we combine line-of-sight velocities from the Abundances and Radial velocity Galactic Origins Survey (ARGOS) spectroscopic survey with proper motions from $Gaia$ DR2, to obtain a sample of $\sim$ 7,000 red clump stars with 3-D velocities. We perform a large scale stellar kinematics study of the Milky Way (MW) bulge to characterize the bulge velocity ellipsoids. We measure the tilt $l_{v}$ of the major-axis of the velocity ellipsoid in the radial-longitudinal velocity plane in 20 fields across the bulge. The tilt or vertex deviation, is characteristic of non-axisymmetric systems and a significant tilt is a robust indicator of non-axisymmetry or bar presence. We compare the observations to the predicted kinematics of an N-body boxy-bulge model formed from dynamical instabilities. In the model, the $l_{v}$ values are strongly correlated with the angle ($α$) between the bulge major-axis and the Sun-Galactic centre line-of-sight. We use a maximum likelihood method to obtain an independent measurement of $α$, from bulge stellar kinematics alone. The most likely value of $α$ given our model is $α= (29 \pm 3)^{\circ}$. In the Baade's window, the metal-rich stars display a larger vertex deviation ($l_{v} = -40^{\circ}$) than the metal-poor stars ($l_{v} = 10^{\circ}$) but we do not detect significant $l_{v}-$metallicity trends in the other fields.

preprint2020arXiv

The Halo's Ancient Metal-Rich Progenitor Revealed with BHB Stars

Using the data from the Sloan Digital Sky Survey and the Gaia satellite, we assemble a pure sample of $\sim$3000 Blue Horizontal Branch (BHB) stars with 7-D information, including positions, velocities and metallicities. We demonstrate that, as traced with BHBs, the Milky Way's stellar halo is largely unmixed and can not be well represented with a conventional Gaussian velocity distribution. A single-component model fails because the inner portions of the halo are swamped with metal-rich tidal debris from an ancient, head-on collision, known as the "Gaia Sausage". Motivated by the data, we build a flexible mixture model which allows us to track the evolution of the halo make-up across a wide range of radii. It is built from two components, one representing the radially anisotropic Sausage stars with their lobed velocity distribution, the other representing a more metal-poor and more isotropic component built up from minor mergers. We show that inside 25 kpc the "Sausage" contributes at least 50 % of the Galactic halo. The fraction of "Sausage" stars diminishes sharply beyond 30 kpc, which is the long-established break radius of the classical stellar halo.

preprint2020arXiv

The tidal remnant of an unusually metal-poor globular cluster

Globular clusters are some of the oldest bound stellar structures observed in the Universe. They are ubiquitous in large galaxies and are believed to trace intense star formation events and the hierarchical build-up of structure. Observations of globular clusters in the Milky Way, and a wide variety of other galaxies, have found evidence for a `metallicity floor', whereby no globular clusters are found with chemical (`metal') abundances below approximately 0.3 to 0.4 per cent of that of the Sun. The existence of this metallicity floor may reflect a minimum mass and a maximum redshift for surviving globular clusters to form, both critical components for understanding the build-up of mass in the universe. Here we report measurements from the Southern Stellar Streams Spectroscopic Survey of the spatially thin, dynamically cold Phoenix stellar stream in the halo of the Milky Way. The properties of the Phoenix stream are consistent with it being the tidally disrupted remains of a globular cluster. However, its metal abundance ([Fe/H] = -2.7) is substantially below that of the empirical metallicity floor. The Phoenix stream thus represents the debris of the most metal-poor globular cluster discovered so far, and its progenitor is distinct from the present-day globular cluster population in the local Universe. Its existence implies that globular clusters below the metallicity floor have probably existed, but were destroyed during Galactic evolution.

preprint2020arXiv

Unresolved stellar companions with Gaia DR2 astrometry

For stars with unresolved companions, motions of the centre of light and that of mass decouple, causing a single-source astrometric model to perform poorly. We show that such stars can be easily detected with the reduced chi2 statistic, or RUWE, provided as part of Gaia DR2. We convert RUWE into the amplitude of the image centroid wobble, which, if scaled by the source distance, is proportional to the physical separation between companions (for periods up to several years). We test this idea on a sample of known spectroscopic binaries and demonstrate that the amplitude of the centroid perturbation scales with the binary period and the mass ratio as expected. We apply this technique to the Gaia DR2 data and show how the binary fraction evolves across the Hertzsprung--Russell diagram. The observed incidence of unresolved companions is high for massive young stars and drops steadily with stellar mass, reaching its lowest levels for white dwarfs. We highlight the elevated binary fraction for the nearby Blue Stragglers and Blue Horizontal Branch stars. We also illustrate how unresolved hierarchical triples inflate the relative velocity signal in wide binaries. Finally, we point out a hint of evidence for the existence of additional companions to the hosts of extrasolar hot jupiters.

preprint2019arXiv

Chemo-dynamical properties of the Anticenter Stream: a surviving disc fossil from a past satellite interaction

Using Gaia DR2, we trace the Anticenter Stream (ACS) in various stellar populations across the sky and find that it is kinematically and spatially decoupled from the Monoceros Ring. Using stars from {\sc lamost} and {\sc segue}, we show that the ACS is systematically more metal-poor than Monoceros by $0.1$ dex with indications of a narrower metallicity spread. Furthermore, the ACS is predominantly populated of old stars ($\sim 10\,\rm{Gyr}$), whereas Monoceros has a pronounced tail of younger stars ($6-10\, \rm{Gyr}$) as revealed by their cumulative age distributions. Put togehter, all of this evidence support predictions from simulations of the interaction of the Sagittarius dwarf with the Milky Way, which argue that the Anticenter Stream (ACS) is the remains of a tidal tail of the Galaxy excited during Sgr's first pericentric passage after it crossed the virial radius, whereas Monoceros consists of the composite stellar populations excited during the more extended phases of the interaction. We suggest that the ACS can be used to constrain the Galactic potential, particularly its flattening, setting strong limits on the existence of a dark disc. Importantly, the ACS can be viewed as a stand-alone fossil of the chemical enrichment history of the Galactic disc.

preprint2019arXiv

Detailed Abundances in the Ultra-faint Magellanic Satellites Carina II and III

We present the first detailed elemental abundances in the ultra-faint Magellanic satellite galaxies Carina II (Car II) and Carina III (Car III). With high-resolution Magellan/MIKE spectroscopy, we determined abundances of nine stars in Car II including the first abundances of an RR Lyrae star in an ultra-faint dwarf galaxy; and two stars in Car III. The chemical abundances demonstrate that both systems are clearly galaxies and not globular clusters. The stars in these galaxies mostly display abundance trends matching those of other similarly faint dwarf galaxies: enhanced but declining [alpha/Fe] ratios, iron-peak elements matching the stellar halo, and unusually low neutron-capture element abundances. One star displays a low outlying [Sc/Fe] = -1.0. We detect a large Ba scatter in Car II, likely due to inhomogeneous enrichment by low-mass AGB star winds. The most striking abundance trend is for [Mg/Ca] in Car II, which decreases from +0.4 to -0.4 and indicates clear variation in the initial progenitor masses of enriching core-collapse supernovae. So far, the only ultra-faint dwarf galaxies displaying a similar [Mg/Ca] trend are likely satellites of the Large Magellanic Cloud. We find two stars with [Fe/H] < -3.5, whose abundances likely trace the first generation of metal-free Population III stars and are well-fit by Population III core-collapse supernova yields. An appendix describes our new abundance uncertainty analysis that propagates line-by-line stellar parameter uncertainties.

preprint2016arXiv

A 10 kpc stellar substructure at the edge of the Large Magellanic Cloud: perturbed outer disk or evidence for tidal stripping?

We report the discovery of a substantial stellar overdensity in the periphery of the Large Magellanic Cloud (LMC), found using public imaging from the first year of the Dark Energy Survey. The structure appears to emanate from the edge of the outer LMC disk at a radius $\approx 13.5$ degrees due north of its centre, and stretches more than $10$ kpc towards the east. It is roughly $1.5$ kpc wide and has an integrated $V$-band luminosity of at least $M_V = -7.4$. The stellar populations in the feature are indistinguishable from those in the outer LMC disk. We attempt to quantify the geometry of the outer disk using simple planar models, and find that only a disk with mild intrinsic ellipticity can simultaneously explain the observed stellar density on the sky and the azimuthal line-of-sight distance profile. We also see possible non-planar behaviour in the outer disk that may reflect a warp and/or flare, as well as deviations that resemble a ring-like structure between $\sim9-12$ degrees from the LMC centre. Based on all these observations, we conclude that our remote, stream-like feature is likely comprised of material that has been stripped from the outskirts of the LMC disk, although we cannot rule out that it represents a transient overdensity in the disk itself. We conduct a simple $N$-body simulation to show that either type of structure could plausibly arise due to the tidal force of the Milky Way; however we also recognize that a recent close interaction between the LMC and the SMC may be the source of the stripping or perturbation. Finally, we observe evidence for extremely diffuse LMC populations extending to radii of $\sim 18.5$ kpc in the disk plane ($\approx 20$ degrees on the sky), corroborating previous spectroscopic detections at comparable distances.

preprint2015arXiv

A Hubble Space Telescope Study of the Enigmatic Milky Way Halo Globular Cluster Crater

We analyze the resolved stellar populations of the faint stellar system, Crater, based on deep optical imaging taken with the Hubble Space Telescope. The HST/ACS-based color-magnitude diagram (CMD) of Crater extends $\sim$4 magnitudes below the oldest main sequence turnoff, providing excellent leverage on Crater's physical properties. Structurally, Crater has a half-light radius of $\sim$20 pc and shows no evidence for tidal distortions. Crater is well-described by a simple stellar population with an age of $\sim$7.5 Gyr, [M/H]$\sim-1.65$, a M$_{\star}\sim10^4$ M$_{\odot}$, M$_{\rm V}\sim -5.3$, located at a distance of d$\sim$ 145 kpc, with modest uncertainties in these properties due to differences in the underlying stellar evolution models. The sparse sampling of stars above the turnoff and sub-giant branch are likely to be 1.0-1.4 M$_{\odot}$ binary star systems (blue stragglers) and their evolved descendants, as opposed to intermediate age main sequence stars. Confusion of these populations highlights a substantial challenge in accurately characterizing sparsely populated stellar systems. Our analysis shows that Crater is not a dwarf galaxy, but instead is an unusually young cluster given its location in the Milky Way's very outer stellar halo. Crater is similar to SMC cluster Lindsay 38, and its position and velocity are in good agreement with observations and models of the Magellanic stream debris, suggesting it may have accreted from the Magellanic Clouds. However, its age and metallicity are also in agreement with the age-metallicity relationships of lower mass dwarf galaxies such as Leo I or Carina. Despite uncertainty over its progenitor system, Crater appears to have been incorporated into the Galaxy more recently than $z\sim1$ (8 Gyr ago), providing an important new constraint on the accretion history of the Milky Way. [abridged]

preprint2015arXiv

Beasts of the Southern Wild: Discovery of nine Ultra Faint satellites in the vicinity of the Magellanic Clouds

We have used the publicly released Dark Energy Survey data to hunt for new satellites of the Milky Way in the Southern hemisphere. Our search yielded a large number of promising candidates. In this paper, we announce the discovery of 9 new unambiguous ultra-faint objects, whose authenticity can be established with the DES data alone. Based on the morphological properties, three of the new satellites are dwarf galaxies, one of which is located at the very outskirts of the Milky Way, at a distance of 380 kpc. The remaining 6 objects have sizes and luminosities comparable to the Segue~1 satellite and can not be classified straightforwardly without follow-up spectroscopic observations. The satellites we have discovered cluster around the LMC and the SMC. We show that such spatial distribution is unlikely under the assumption of isotropy, and, therefore, conclude that at least some of the new satellites must have been associated with the Magellanic Clouds in the past.

preprint2015arXiv

Indication of Gamma-ray Emission from the Newly Discovered Dwarf Galaxy Reticulum II

We present a search for gamma-ray emission from the direction of the newly discovered dwarf galaxy Reticulum II. Using Fermi-LAT data, we detect a signal that exceeds expected backgrounds between ~2-10 GeV and is consistent with annihilation of dark matter for particle masses less than a few x 10^2 GeV. Modeling the background as a Poisson process based on Fermi-LAT diffuse models, and taking into account trials factors, we detect emission with p-value less than 9.8 x 10^-5 (>3.7 sigma). An alternative, model-independent treatment of background reduces the significance, raising the p-value to 9.7 x 10^-3 (2.3 sigma). Even in this case, however, Reticulum II has the most significant gamma-ray signal of any known dwarf galaxy. If Reticulum II has a dark matter halo that is similar to those inferred for other nearby dwarfs, the signal is consistent with the s-wave relic abundance cross section for annihilation.

preprint2015arXiv

Kinematics and chemistry of recently discovered Reticulum 2 and Horologium 1 dwarf galaxies

Photometry alone is not sufficient to unambiguously distinguish between ultra-faint star clusters and dwarf galaxies because of their overlap in morphological properties. Here we report on VLT/GIRAFFE spectra of candidate member stars in two recently discovered ultra-faint satellites Reticulum 2 and Horologium 1, obtained as part of the ongoing Gaia-ESO Survey. We identify 18 members in Reticulum 2 and 5 in Horologium 1. We find Reticulum 2 to have a velocity dispersion of ~3.22 km/s, implying a M/L ratio of ~ 500. We have inferred stellar parameters for all candidates and we find Reticulum 2 to have a mean metallicity of [Fe/H] = -2.46+/-0.1, with an intrinsic dispersion of ~ 0.29, and is alpha-enhanced to the level of [alpha/Fe]~0.4. We conclude that Reticulum 2 is a dwarf galaxy. We also report on the serendipitous discovery of four stars in a previously unknown stellar substructure near Reticulum 2 with [Fe/H] ~ -2 and V_hel ~ 220 km/s, far from the systemic velocity of Reticulum 2. For Horologium 1 we infer a velocity dispersion of 4.9^{+2.8}_{-0.9} km/s and a consequent M/L ratio of ~ 600, leading us to conclude that Horologium 1 is also a dwarf galaxy. Horologium 1 is slightly more metal-poor than Reticulum 2 [Fe/H] = -2.76 +/- 0.1 and is similarly alpha-enhanced: [alpha/Fe] ~ 0.3. Despite a large error-bar, we also measure a significant spread of metallicities of 0.17 dex which strengthen the evidence that Horologium 1 is indeed a dwarf galaxy. The line-of-sight velocity of Reticulum 2 is offset by some 100 km/s from the prediction of the orbital velocity of the LMC, thus making its association with the Cloud uncertain. However, at the location of Horologium 1, both the backward integrated LMC's orbit and the LMC's halo are predicted to have radial velocities similar to that of the dwarf. Therefore, it is likely that Horologium 1 is or once was a member of the Magellanic Family.

preprint2015arXiv

Magellan/M2FS Spectroscopy of the Reticulum 2 Dwarf Spheroidal Galaxy

We present results from spectroscopic observations with the Michigan/Magellan Fiber System (M2FS) of 182 stellar targets along the line of sight to the newly-discovered `ultrafaint' object Reticulum 2 (Ret 2). For 38 of these targets, the spectra are sufficient to provide simultaneous estimates of line-of-sight velocity ($v_{\rm los}$, median random error $δ_{v_{\rm los}}=1.3$ km s$^{-1}$), effective temperature ($T_{\rm eff}$, $δ_{T_{\rm eff}}=464$ K), surface gravity ($\log g$, $δ_{\rm logg}=0.54$ dex) and iron abundance ([Fe/H], $δ_{\mathrm{[Fe/H]}}=0.45$ dex). We use these results to confirm 18 stars as members of Ret 2. From the member sample we estimate a velocity dispersion of $σ_{v_{\rm los}}=3.6_{-0.6}^{+0.9}$ km s$^{-1}$ about a mean of $\langle v_{\rm los}\rangle =64.8_{-1.0}^{+1.1}$ km s$^{-1}$ in the solar rest frame ($\sim -90.9$ km s$^{-1}$ in the Galactic rest frame), and a metallicity dispersion of $σ_{\rm [Fe/H]}=0.50_{-0.13}^{+0.17}$ dex about a mean of $\langle \mathrm{[Fe/H]} \rangle =-2.67_{-0.34}^{+0.34}$. These estimates marginalize over possible velocity and metallicity gradients, which are consistent with zero. Our results place Ret 2 on chemodynamical scaling relations followed by the Milky Way's dwarf-galactic satellites. Under assumptions of dynamic equilibrium and negligible contamination from binary stars---both of which must be checked with deeper imaging and repeat spectroscopic observations---the estimated velocity dispersion suggests a dynamical mass of $M(R_{\rm h})\approx 5R_{\rm h}σ_{v_{\rm los}}^2/(2G)=2.4_{-0.8}^{+1.3}\times 10^5$ $M_{\odot}$ enclosed within projected halflight radius $R_{\rm h}\sim 32$ pc, with mass-to-light ratio $\approx 2M(R_{\rm h})/L_{\rm V}=462_{-157}^{+264}$ in solar units.

preprint2015arXiv

Magellan/M2FS Spectroscopy of Tucana 2 and Grus 1

We present results from spectroscopic observations with the Michigan/Magellan Fiber System (M2FS) of $147$ stellar targets along the line of sight to the newly-discovered `ultrafaint' stellar systems Tucana 2 (Tuc 2) and Grus 1 (Gru 1). Based on simultaneous estimates of line-of-sight velocity and stellar-atmospheric parameters, we identify 8 and 7 stars as probable members of Tuc 2 and and Gru 1, respectively. Our sample for Tuc 2 is sufficient to resolve an internal velocity dispersion of $8.6_{-2.7}^{+4.4}$ km s$^{-1}$ about a mean of $-129.1_{-3.5}^{+3.5}$ km s$^{-1}$ (solar rest frame), and to estimate a mean metallicity of [Fe/H]= $-2.23_{-0.12}^{+0.18}$. These results place Tuc 2 on chemodynamical scaling relations followed by dwarf galaxies, suggesting a dominant dark matter component with dynamical mass $2.7_{-1.3}^{+3.1}\times 10^6$ $\mathrm{M}_{\odot}$ enclosed within the central $\sim 160$ pc, and dynamical mass-to-light ratio $1900_{-900}^{+2200}$ $\mathrm{M}_{\odot}/L_{V,\odot}$. For Gru 1 we estimate a mean velocity of $-140.5_{-1.6}^{+2.4}$ km s$^{-1}$ and a mean metallicity of [Fe/H]=$-1.42_{-0.42}^{+0.55}$, but our sample does not resolve Gru 1's velocity dispersion. The radial coordinates of Tuc 2 and Gru 1 in Galactic phase space suggest that their orbits are among the most energetic within distance $\leq 300$ kpc. Moreover, their proximity to each other in this space arises naturally if both objects are trailing the Large Magellanic Cloud.

preprint2014arXiv

Balancing mass and momentum in the Local Group

In the rest frame of the Local Group (LG), the total momentum of the Milky Way (MW) and Andromeda (M31) should balance to zero. We use this fact to constrain new solutions for the solar motion with respect to the LG centre-of-mass, the total mass of the LG, and the individual masses of M31 and the MW. Using the set of remote LG galaxies at $>350$ kpc from the MW and M31, we find that the solar motion has amplitude $V_{\odot}=299\pm 15 {\rm ~km~s^{-1}}$ in a direction pointing toward galactic longitude $l_{\odot}=98.4^{\circ}\pm 3.6^{\circ}$ and galactic latitude $b_{\odot}=-5.9^{\circ}\pm 3.0^{\circ}$. The velocities of M31 and the MW in this rest frame give a direct measurement of their mass ratio, for which we find $\log_{10} (M_{\rm M31}/M_{\rm MW})=0.36 \pm 0.29$. We combine these measurements with the virial theorem to estimate the total mass within the LG as $M_{\rm LG}=(2.5\pm 0.4)\times 10^{12}~{\rm M}_{\odot}$. Our value for $M_{\rm LG}$ is consistent with the sum of literature values for $M_{\rm MW}$ and $M_{\rm M31}$. This suggests that the mass of the LG is almost entirely located within the two largest galaxies rather than being dispersed on larger scales or in a background medium. The outskirts of the LG are seemingly rather empty. Combining our measurement for $M_{\rm LG}$ and the mass ratio, we estimate the individual masses of the MW and M31 to be $M_{\rm MW}=(0.8\pm 0.5)\times 10^{12}~{\rm M}_{\odot}$ and $M_{\rm M31}=(1.7\pm 0.3)\times 10^{12}~{\rm M}_{\odot}$, respectively. Our analysis favours M31 being more massive than the MW by a factor of $\sim$2.3, and the uncertainties allow only a small probability (9.8%) that the MW is more massive. This is consistent with other properties such as the maximum rotational velocities, total stellar content, and numbers of globular clusters and dwarf satellites, which all suggest that $M_{\rm M31}/M_{\rm MW}>1$.

preprint2014arXiv

Strong RR Lyrae excess in the Hercules-Aquila Cloud

We map the large-scale sub-structure in the Galactic stellar halo using accurate 3D positions of ~14,000 RR Lyrae reported by the Catalina Sky Survey. In the heliocentric distance range of 10-25 kpc, in the region of the sky approximately bounded by 30° < l < 55° and -45° < b < -25°, there appears to be a strong excess of RRab stars. This overdensity, peaking at 18 kpc, is most likely associated with the so-called Hercules-Aquila Cloud, previously detected using Main Sequence tracers at similar distances in the Sloan Digital Sky Survey data. Our analysis of the period-amplitude distribution of RR Lyrae in this region indicates that the HAC is dominated by the Oosterhoff I type population. By comparing the measured RR Lyrae number density to models of a smooth stellar halo, we estimate the significance of the observed excess and provide an updated estimate of the total luminosity of the Cloud's progenitor.

preprint2014arXiv

Touching The Void: A Striking Drop in Stellar Halo Density Beyond 50 kpc

We use A-type stars selected from Sloan Digital Sky Survey data release 9 photometry to measure the outer slope of the Milky Way stellar halo density profile beyond 50 kpc. A likelihood-based analysis is employed that models the ugr photometry distribution of blue horizontal branch (BHB) and blue straggler (BS) stars. In the magnitude range, 18.5 < g < 20.5, these stellar populations span a heliocentric distance range of: 10 kpc < D_BS < 75 kpc, 40 kpc < D_BHB < 100 kpc. Contributions from contaminants, such as QSOs, and the effect of photometric uncertainties, are also included in our modeling procedure. We find evidence for a very steep outer halo profile, with power-law index alpha ~ 6 beyond Galactocentric radii r=50 kpc, and even steeper slopes favored (alpha ~ 6-10) at larger radii. This result holds true when stars belonging to known overdensities, such as the Sagittarius stream, are included or excluded. We show that, by comparison to numerical simulations, stellar halos with shallower slopes at large distances tend to have more recent accretion activity. Thus, it is likely that the Milky Way has undergone a relatively quiet accretion history over the past several Gyr. Our measurement of the outer stellar halo profile may have important implications for dynamical mass models of the Milky Way, where the tracer density profile is strongly degenerate with total mass-estimates.

preprint2013arXiv

Sagittarius stream 3-d kinematics from SDSS Stripe 82

Using multi-epoch observations of the Stripe 82 region done by Sloan Digital Sky Survey, we measure precise statistical proper motions of the stars in the Sagittarius stellar stream. The multi-band photometry and SDSS radial velocities allow us to efficiently select Sgr members and thus enhance the proper motion precision to ~0.1 mas/yr. We measure separately the proper motion of a photometrically selected sample of the main sequence turn-off stars, as well as of a spectroscopically selected Sgr giants. The data allow us to determine the proper motion separately for the two Sgr streams in the South found in Koposov et al.(2012). Together with the precise velocities from SDSS, our proper motion provide exquisite constraints of the 3-D motions of the stars in the Sgr streams.

preprint2012arXiv

A statistical method for measuring the Galactic potential and testing gravity with cold tidal streams

We introduce the Minimum Entropy Method, a simple statistical technique for constraining the Milky Way gravitational potential and simultaneously testing different gravity theories directly from 6D phase-space surveys and without adopting dynamical models. We demonstrate that orbital energy distributions that are separable (i.e. independent of position) have an associated entropy that increases under wrong assumptions about the gravitational potential and/or gravity theory. Of known objects, `cold' tidal streams from low-mass progenitors follow orbital distributions that most nearly satisfy the condition of separability. Although the orbits of tidally stripped stars are perturbed by the progenitor's self-gravity, systematic variations of the energy distribution can be quantified in terms of the cross-entropy of individual tails, giving further sensitivity to theoretical biases in the host potential. The feasibility of using the Minimum Entropy Method to test a wide range of gravity theories is illustrated by evolving restricted N-body models in a Newtonian potential and examining the changes in entropy introduced by Dirac, MONDian and f(R) gravity modifications.

preprint2012arXiv

Hyper-luminous Reddened Broad-Line Quasars at z~2 from the VISTA Hemisphere Survey and WISE All Sky Survey

We present the first sample of spectroscopically confirmed heavily reddened broad-line quasars selected using the new near infra-red VISTA Hemisphere Survey and WISE All-Sky Survey. Observations of four candidates with (J-K)>2.5 and K<=16.5 over ~180 sq deg, leads to confirmation that two are highly dust-reddened broad-line Type 1 quasars at z~2. The typical dust extinctions are AV ~ 2-2.5 mags. We measure black-hole masses of ~10^9 M0 and extinction corrected bolometric luminosities of ~10^47 erg/s, making these among the brightest Type 1 quasars currently known. Despite this, these quasars lie well below the detection limits of wide-field optical surveys like the SDSS with iAB > 22. We also present WISE photometry at 3--22um, for our full sample of spectroscopically confirmed reddened quasars including those selected from the UKIDSS Large Area Survey (Banerji et al. 2012a). We demonstrate that the rest-frame infrared SEDs of these reddened quasars are similar to UV-luminous Type 1 quasars with significant hot dust emission and starburst quasar hosts like Mrk231. The average 12um flux density of our reddened quasars is similar to that of the recently discovered HyLIRG WISE1814+3412 (z=2.452) at similar redshifts, with two of our reddened quasars also having comparable 22um flux densities to this extreme HyLIRG. These optically faint, heavily reddened broad-line quasars are therefore among the most mid infrared luminous galaxies at z~2, now being discovered using WISE.

preprint2012arXiv

Photometric study of open clusters Berkeley 96, Berkeley 97, King 12, NGC 7261, NGC 7296, and NGC 7788

The BVRcIc CCD photometry in the fields of six open clusters toward the Perseus spiral arm is presented. These data along with JHKs magnitudes taken from 2MASS catalog, have been used to determine cluster's ages, distances and color excesses. In addition, the gaps in mass function of Be 97, King 12 and NGC 7788 clusters have been revealed in mass intervals of 1.3-1.5, 1.4-1.6 and 1.5-1.7 solar masses, respectively.

preprint2012arXiv

The coupling between the core/cusp and missing satellite problems

We calculate the energy that baryons must inject in cold dark matter (CDM) haloes in order to remove centrally-divergent DM cusps on scales relevant to observations of dwarf spheroidal galaxies (dSphs). We estimate that the CDM haloes often associated with the Milky Way's dSphs (M_vir/M_\odot \sim 10^{9-10}) require ΔE/erg \sim 10^{53-55} in order to form cores on scales comparable to the luminous size of these galaxies. While supernova type II (SNeII) explosions can in principle generate this energy, the actual contribution is limited by the low star formation efficiency implied by the abundance of luminous satellites. Considering that CDM's well-known `core/cusp' and `missing satellite' problems place opposing demands on star formation efficiencies, existing observational evidences for large cores in the most luminous dSphs require that CDM models invoke some combination of the following: (i) efficient (of order unity) coupling of SNeII energy into dark matter particles, (ii) star formation histories peaking at unexpectedly high redshifts (z>6), (iii) a top-heavy stellar IMF, and/or (iv) substantial satellite disruption or other stochastic effects to ease the substructure abundance constraints. Our models show that the tension between CDM problems on small scales would increase if cored DM profiles were to be found in fainter dwarves.

preprint2012arXiv

The Sagittarius Streams in the Southern Galactic Hemisphere

The structure of the Sagittarius stream in the Southern Galactic hemisphere is analysed with the Sloan Digital Sky Survey Data Release 8. Parallel to the Sagittarius tidal track, but ~ 10deg away, there is another fainter and more metal-poor stream. We provide evidence that the two streams follow similar distance gradients but have distinct morphological properties and stellar populations. The brighter stream is broader, contains more metal-rich stars and has a richer colour-magnitude diagram with multiple turn-offs and a prominent red clump as compared to the fainter stream. Based on the structural properties and the stellar population mix, the stream configuration is similar to the Northern "bifurcation". In the region of the South Galactic Cap, there is overlapping tidal debris from the Cetus Stream, which crosses the Sagittarius stream. Using both photometric and spectroscopic data, we show that the blue straggler population belongs mainly to Sagittarius and the blue horizontal branch stars belong mainly to the Cetus stream in this confused location in the halo.

preprint2011arXiv

Accurate Stellar Kinematics at Faint Magnitudes: application to the Bootes~I dwarf spheroidal galaxy

We develop, implement and characterise an enhanced data reduction approach which delivers precise, accurate, radial velocities from moderate resolution spectroscopy with the fibre-fed VLT/FLAMES+GIRAFFE facility. This facility, with appropriate care, delivers radial velocities adequate to resolve the intrinsic velocity dispersions of the very faint dSph dwarf galaxies. Importantly, repeated measurements let us reliably calibrate our individual velocity errors ($0.2 \leq δ_V\leq 5$ km s$^{-1}$) and directly detect stars with variable radial velocities. We show, by application to the Bootes-1 dwarf spheroidal, that the intrinsic velocity dispersion of this system is significantly below 6.5\,km/s reported by previous studies. Our data favor a two-population model of Bootes-1, consisting of a majority `cold' stellar component, with velocity dispersion $2.4^{+0.9}_{-0.5}$\,km/s, and a minority `hot' stellar component, with velocity dispersion $\sim 9$\,km/s, although we can not completely rule out a single component distribution with velocity dispersion $4.6^{0.8}_{-0.6}$\,km/s. We speculate this complex velocity distribution actually reflects the distribution of velocity anisotropy in Bootes-1, which is a measure of its formation processes.

preprint2011arXiv

Quantifying Kinematic Substructure in the Milky Way's Stellar Halo

We present and analyze the positions, distances, and radial velocities for over 4000 blue horizontal-branch (BHB) stars in the Milky Way's halo, drawn from SDSS DR8. We search for position-velocity substructure in these data, a signature of the hierarchical assembly of the stellar halo. Using a cumulative "close pair distribution" (CPD) as a statistic in the 4-dimensional space of sky position, distance, and velocity, we quantify the presence of position-velocity substructure at high statistical significance among the BHB stars: pairs of BHB stars that are close in position on the sky tend to have more similar distances and radial velocities compared to a random sampling of these overall distributions. We make analogous mock-observations of 11 numerical halo formation simulations, in which the stellar halo is entirely composed of disrupted satellite debris, and find a level of substructure comparable to that seen in the actually observed BHB star sample. This result quantitatively confirms the hierarchical build-up of the stellar halo through a signature in phase (position-velocity) space. In detail, the structure present in the BHB stars is somewhat less prominent than that seen in most simulated halos, quite possibly because BHB stars represent an older sub-population. BHB stars located beyond 20 kpc from the Galactic center exhibit stronger substructure than at $\rm r_{gc} < 20$ kpc.

preprint2010arXiv

Constraining the Milky Way potential with a 6-D phase-space map of the GD-1 stellar stream

The narrow GD-1 stream of stars, spanning 60 deg on the sky at a distance of ~10 kpc from the Sun and ~15 kpc from the Galactic center, is presumed to be debris from a tidally disrupted star cluster that traces out a test-particle orbit in the Milky Way halo. We combine SDSS photometry, USNO-B astrometry, and SDSS and Calar Alto spectroscopy to construct a complete, empirical 6-dimensional phase-space map of the stream. We find that an eccentric orbit in a flattened isothermal potential describes this phase-space map well. Even after marginalizing over the stream orbital parameters and the distance from the Sun to the Galactic center, the orbital fit to GD-1 places strong constraints on the circular velocity at the Sun's radius V_c=224 \pm 13 km/s and total potential flattening q_Φ=0.87^{+0.07}_{-0.04}. When we drop any informative priors on V_c the GD-1 constraint becomes V_c=221 \pm 18 km/s. Our 6-D map of GD-1 therefore yields the best current constraint on V_c and the only strong constraint on q_Φat Galactocentric radii near R~15 kpc. Much, if not all, of the total potential flattening may be attributed to the mass in the stellar disk, so the GD-1 constraints on the flattening of the halo itself are weak: q_{Φ,halo}>0.89 at 90% confidence. The greatest uncertainty in the 6-D map and the orbital analysis stems from the photometric distances, which will be obviated by Gaia.

preprint2009arXiv

A quantitative explanation of the observed population of Milky Way satellite galaxies

We revisit the well known discrepancy between the observed number of Milky Way (MW) dwarf satellite companions and the predicted population of cold dark matter (CDM) sub-halos, in light of the dozen new low luminosity satellites found in SDSS imaging data and our recent calibration of the SDSS satellite detection efficiency, which implies a total population far larger than these dozen discoveries. We combine a dynamical model for the CDM sub-halo population with simple, physically motivated prescriptions for assigning stellar content to each sub-halo, then apply observational selection effects and compare to the current observational census. As expected, models in which the stellar mass is a constant fraction F(Omega_b/Omega_m) of the sub-halo mass M_sat at the time it becomes a satellite fail for any choice of F. However, previously advocated models that invoke suppression of gas accretion after reionization in halos with circular velocity v_c <~ 35 km/s can reproduce the observed satellite counts for -15 < M_V < 0, with F ~ 10^{-3}. Successful models also require strong suppression of star formation BEFORE reionization in halos with v_c <~ 10 km/s; models without pre-reionization suppression predict far too many satellites with -5 < M_V < 0. Our models also reproduce the observed stellar velocity dispersions ~ 5-10 km/s of the SDSS dwarfs given the observed sizes of their stellar distributions, and model satellites have M(<300 pc) ~ 10^7 M_sun as observed even though their present day total halo masses span more than two orders of magnitude. Our modeling shows that natural physical mechanisms acting within the CDM framework can quantitatively explain the properties of the MW satellite population as it is presently known, thus providing a convincing solution to the `missing satellite' problem.

preprint2009arXiv

Luminosity function and radial distribution of Milky Way Satellites in a LCDM Universe

We study the luminosity function and the radial distribution of satellite galaxies within Milky Way sized haloes as predicted in Cold Dark Matter based models of galaxy formation, making use of numerical N-body techniques as well as three different semi-analytic model (SAMs) galaxy formation codes. We extract merger trees from very high-resolution dissipationless simulations of four Galaxy-sized DM haloes, and use these as common input for the semi-analytic models. We present a detailed comparison of our predictions with the observational data recently obtained on the Milky Way satellite luminosity function (LF). We find that semi-analytic models with rather standard astrophysical ingredients are able to reproduce the observed luminosity function over six orders of magnitude in luminosity, down to magnitudes as faint as M_V=-2. We also perform a comparison with the actual observed number of satellites as a function of luminosity, by applying the selection criteria of the SDSS survey to our simulations instead of correcting the observations for incompleteness. Using this approach we again find good agreement for both the luminosity and radial distributions of MW satellites. We investigate which physical processes in our models are responsible for shaping the predicted satellite LF, and find that tidal destruction, suppression of gas infall by a photo-ionizing background, and supernova feedback all make important contributions. We conclude that the number and luminosity of Milky Way satellites can be naturally accounted for within the (Lambda) Cold Dark Matter paradigm, and this should no longer be considered a problem.

preprint2007arXiv

GEMS: Galaxy fitting catalogues and testing parametric galaxy fitting codes

In the context of measuring structure and morphology of intermediate redshift galaxies with recent HST/ACS surveys, we tune, test, and compare two widely used fitting codes (GALFIT and GIM2D) for fitting single-component Sersic models to the light profiles of both simulated and real galaxy data. We find that fitting accuracy depends sensitively on galaxy profile shape. Exponential disks are well fit with Sersic models and have small measurement errors, whereas fits to de Vaucouleurs profiles show larger uncertainties owing to the large amount of light at large radii. We find that both codes provide reliable fits and little systematic error, when the effective surface brightness is above that of the sky. Moreover, both codes return errors that significantly underestimate the true fitting uncertainties, which are best estimated with simulations. We find that GIM2D suffers significant systematic errors for spheroids with close companions owing to the difficulty of effectively masking out neighboring galaxy light; there appears to be no work around to this important systematic in GIM2D's current implementation. While this crowding error affects only a small fraction of galaxies in GEMS, it must be accounted for in the analysis of deeper cosmological images or of more crowded fields with GIM2D. In contrast, GALFIT results are robust to the presence of neighbors because it can simultaneously fit the profiles of multiple companions thereby deblending their effect on the fit to the galaxy of interest. We find GALFIT's robustness to nearby companions and factor of >~20 faster runtime speed are important advantages over GIM2D for analyzing large HST/ACS datasets. Finally we include our final catalog of fit results for all 41,495 objects detected in GEMS.