Search

The improved astrometry precision of Gaia-eDR3 allows us to perform a detailed study of the Upper Scorpius OB association and revisit its spatial, kinematic, and age substructure. We achieve this by combining clustering techniques and complementing with age estimations based on Gaia photometry. Our census retrieves 3661 candidate members for Upper Scorpius with contamination $\sim$9\%. We also extract an astrometrically clean sample of 3004 sources with contamination $\sim$6\%. We show that Upper Scorpius can be divided into at least 3 main kinematic groups. We systematically investigate and characterize the Upper Scorpius' internal structure, revealing that at least $\sim 34\%$ of its stellar populations are contained in 7 spatial substructures, with well defined boundaries, kinematics and relative ages, with suggested names: $\pi$ Scorpii (20 $^{\pm2}_{\pm1}$ Myr), $\alpha$ Scorpii (14$^{\pm2}_{\pm1}$ Myr), $\delta$ Scorpii (9$^{\pm2}_{\pm1}$ Myr), $\beta$ Scorpii (8$^{\pm1}_{\pm1}$ Myr), $\omega$ Scorpii (8$^{\pm1}_{\pm1}$ Myr), $\nu$ Scorpii (7$^{\pm1}_{\pm1}$ Myr), after their brightest member, and the well known $\rho$ Ophiuchi (4$^{\pm1}_{\pm1}$ Myr). We find a clear correlation in (1) density-age, providing an empirical expansion law to be tested in other associations, and (2) tangential velocity-age, providing constrains on the dynamics of these substructures and the position of potential past triggering events. We estimate the time at which 4 potential supernovae events occurred in Upper Scorpius. Based on these findings, we tie together previous work on the region and suggest a star formation history with unprecedented temporal resolution., Comment: Accepted for publication in MNRAS, 2023 February 15

Although not predicted by standard stellar evolution, the surface abundance of light elements, such as lithium (Li), carbon, and nitrogen, changes during the red giant branch (RGB) as a result of extra-mixing. This is associated usually with thermohaline mixing acting after the RGB bump. Peculiar Li-enriched RGB stars might also be related to either enhanced extra-mixing or pollution from external sources. We measure the Li abundance and carbon isotopic ratio 12C/13C in a sample of 166 field red giants with -0.3, 15 pages, 11 figures, accepted for publication in A&A

While most moving groups are young and nearby, a small number have been identified in the Galactic halo. Understanding the origin and evolution of these groups is an important piece of reconstructing the formation history of the halo. Here we report on our analysis of three putative halo moving groups: G03-37, G18-39, and G21-22. Based on Gaia EDR3 data, the stars associated with each group show some scatter in velocity (e.g., Toomre diagram) and integrals of motion (energy, angular momentum) spaces, counter to expectations of moving-group stars. We choose the best candidate of the three groups, G21-22, for follow-up chemical analysis based on high-resolution spectroscopy of six presumptive members. Using a new Python code that uses a Bayesian method to self-consistently propagate uncertainties from stellar atmosphere solutions in calculating individual abundances and spectral synthesis, we derive the abundances of $\alpha$- (Mg, Si, Ca, Ti), Fe-peak (Cr, Sc, Mn, Fe, Ni), odd-$Z$ (Na, Al, V), and neutron-capture (Ba, Eu) elements for each star. We find that the G21-22 stars are not chemically homogeneous. Based on the kinematic analysis for all three groups and the chemical analysis for G21-22, we conclude the three are not genuine moving groups. The case for G21-22 demonstrates the benefit of combining kinematic and chemical information in identifying conatal populations when either alone may be insufficient. Comparing the integrals of motion and velocities of the six G21-22 stars with those of known structures in the halo, we tentatively associate them with the Gaia-Enceladus accretion event., Comment: 21 pages, 12 figures, seven tables; accepted to AJ

We present observational constraints for the initial-to-final mass relation (IFMR) derived from 11 white dwarfs (WDs) in wide binaries (WBs) that contain a turnoff/subgiant primary. Because the components of WBs are coeval to a good approximation, the age of the WD progenitor can be determined from the study of its wide companion. However, previous works that used WBs to constrain the IFMR suffered from large uncertainties in the initial masses because their MS primaries are difficult to age-date with good precision. Our selection of WBs with slightly evolved primaries avoids this problem by restricting to a region of parameter space where isochrone ages are significantly easier to determine with precision. The WDs of two of our originally selected binaries were found to be close double degenerates, and are not used in the IFMR analysis. We obtained more precise constraints than existing ones in the mass range 1-2 M$_{\odot}$, corresponding to a previously poorly constrained region of the IFMR. Having introduced the use of turnoff/subgiant-WD binaries, the study of the IFMR is not limited anymore by the precision in initial mass, but now the pressure is on final mass, i.e., the mass of the WD today. Looking at the full dataset, our results would suggest a relatively large dispersion in the IFMR at low initial masses. More precise determinations of the mass of the WD components of our targets are necessary for settling this question., Accepted for publication on ApJ. 20 pages, 12 figures

Our view of the variety of stellar structures pervading the local Milky Way has been transformed by the application of clustering algorithms to the Gaia catalog. In particular, several stellar streams have been recently discovered that are comprised of hundreds to thousands of stars and span several hundred parsecs. We analyze one such structure, Theia 456, a low-density stellar stream extending nearly 200 pc and 20$^{\circ}$ across the sky. By supplementing Gaia astrometric data with spectroscopic metallicities from LAMOST and photometric rotation periods from the Zwicky Transient Facility (ZTF) and the Transiting Exoplanet Survey Satellite (TESS), we establish Theia 456's radial velocity coherence, and we find strong evidence that members of Theia 456 have a common age ($\simeq$175 Myr), common dynamical origin, and formed from chemically homogeneous pre-stellar material ([Fe/H] = $-$0.07 dex). Unlike well-known stellar streams in the Milky Way, which are in its halo, Theia 456 is firmly part of the thin disk. If our conclusions about Theia 456 can be applied to even a small fraction of the remaining $\simeq$8300 independent structures in the Theia catalog, such low-density stellar streams may be ubiquitous. We comment on the implications this has for the nature of star-formation throughout the Galaxy., Comment: 12 pages, 10 figures, 1 table; accepted for publication in AJ

High lithium-7 ($\mathrm{^7Li}$) abundances in giants are indicative of non-standard physical processes affecting the star. Mechanisms that could produce this signature include contamination from an external source, such as planets, or internal production and subsequent mixing to the stellar surface. However, distinguishing between different families of solutions has proven challenging, and there is no current consensus model that explains all the data. The lithium-6 ($\mathrm{^6Li}$) abundance may be a potentially important discriminant, as the relative $\mathrm{^6Li}$ and $\mathrm{^7Li}$ abundances are expected to be different if the enrichment were to come from internal production or from engulfment. In this work, we model the $\mathrm{^6Li}$ and $\mathrm{^7Li}$ abundances of different giants after the engulfment of a substellar mass companion. Given that $\mathrm{^6Li}$ is more strongly affected by Galactic chemical evolution than $\mathrm{^7Li}$, $\mathrm{^6Li}$ is not a good discriminant at low metallicities, where it is expected to be low in both star and planet. For modeled metallicities ([Fe/H]$>-0.5$), we use a "best case" initial $\mathrm{^6Li/^7Li}$ ratio equal to the solar value. $\mathrm{^6Li}$ increases significantly after the engulfment of a companion. However, at metallicities close to solar and higher, the $\mathrm{^6Li}$ signal does not last long in the stellar surface. As such, detection of surface $\mathrm{^6Li}$ in metal-rich red giants would most likely indicate the action of a mechanism for $\mathrm{^6Li}$-enrichment other than planet engulfment. At the same time, $\mathrm{^6Li}$ should not be used to reject the hypothesis of engulfment in a $\mathrm{^7Li}$-enriched giant or to support a particular $\mathrm{^7Li}$-enhancement mechanism.

The S2 stream is a kinematically cold stream that is plunging downwards through the Galactic disc. It may be part of a hotter and more diffuse structure called the Helmi stream. We present a multi-instrument chemical analysis of the stars in the metal-poor S2 stream using both high- and low-resolution spectroscopy, complemented with a re-analysis of the archival data to give a total sample of 62 S2 members. Our high-resolution program provides alpha-elements (C, Mg, Si, Ca and Ti), iron-peak elements (V, Cr, Mn, Fe, Ni), n-capture process elements (Sr, Ba) and other elements such as Li, Na, Al, and Sc for a subsample of S2 objects. We report coherent abundance patterns over a large metallicity spread (~1 dex) confirming that the S2 stream was produced by a disrupted dwarf galaxy. The combination of S2's $\alpha$-elements displays a mildly decreasing trend with increasing metallicity which can be tentatively interpreted as a ``knee'' at [Fe/H], Comment: Accepted for publitacion in MNRAS. 20 pages and 12 figures

The abundance ratios of some chemical species have been found to correlate with stellar age, leading to the possibility of using stellar atmospheric abundances as stellar age indicators. These chemical clocks have been calibrated with solar-twins, open clusters and red giants, but it remains to be seen whether they can be effective at identifying coeval stars in a field population that spans a broad parameter space (i.e., the promise of chemical tagging). Since the components of wide binaries are known to be stars of common origins, they constitute ideal laboratories for testing the usefulness of chemical clocks for the age dating of field stars. We determined the abundances of a new sample of 5 binaries and collected data for other 31 systems from the literature in order to test the applicability of chemical clocks. We recover the well known result that the components of wide binaries have more consistent chemistry than that of random pairs. However, we also show for the first time that abundance ratios designed as chemical clocks are even more consistent among the components of wide binaries than their [X/Fe] ratios. Not only that, but the special case of the pair HIP 34426/HIP 34407 may indicate that chemical clocks are consistent for coeval stars even when the individual abundances are not. If the assumption that chemical clocks are reliable age indicators is correct, this would constitute first quantitative, statistically significant evidence that the components of wide binaries in the Galactic field are indeed coeval, validating a large body of published work that relies on that to be the case. Moreover, our results provide strong evidence that chemical clocks indeed carry important information about stellar birthplaces and chemical evolution, and thus we propose that including them in chemical tagging efforts may facilitate the identification of nowadays dissolved stellar groups., 16 pages, 9 figures. Accepted for publication in ApJ

From our recent catalog based on the first Gaia data release (TGAS), we select wide binaries in which both stars have been observed by the Radial Velocity Experiment (RAVE) or the Large Sky Area Multi-Object Fibre Spectroscopic Telescope (LAMOST). Using RAVE and LAMOST metallicities and RAVE Mg, Al, Si, Ti, and Fe abundances, we find that the differences in the metallicities and elemental abundances of components of wide binaries are consistent with being due to observational uncertainties, in agreement with previous results for smaller and more restricted samples. The metallicity and elemental abundance consistency between wide binary components presented in this work confirms their common origin and bolsters the status of wide binaries as "mini-open clusters." Furthermore, this is evidence that wide binaries are effectively co-eval and co-chemical, supporting their use for e.g., constraining age-activity-rotation relations, the initial-final mass relation for white dwarfs, and M-dwarf metallicity indicators. Additionally, we demonstrate that the common proper motion, common parallax pairs in TGAS with the most extreme separations (s $\gtrsim$ 0.1 pc) typically have inconsistent metallicities, radial velocities or both and are therefore likely to be predominantly comprised of random alignments of unassociated stars with similar astrometry, in agreement with our previous results. Finally, we propose that wide binaries form an ideal data set with which to test chemical tagging as a method to identify stars of common origin, particularly because the stars in wide binaries span a wide range of metallicities, much wider than that spanned by nearby open clusters., Comment: 15 pages, 11 figures, 3 tables, accepted for publication in MNRAS

Half of all stars reside in binary systems, many of which have orbital separations in excess of 1000 AU. Such binaries are typically identified in astrometric catalogs by matching the proper motions vectors of close stellar pairs. We present a fully Bayesian method that properly takes into account positions, proper motions, parallaxes, and their correlated uncertainties to identify widely separated stellar binaries. After applying our method to the >2 × 106stars in the Tycho-Gaia astrometric solution from Gaia DR1, we identify over 6000 candidate wide binaries. For those pairs with separations less than 40,000 AU, we determine the contamination rate to be ≈5%. This sample has an orbital separation (a) distribution that is roughly flat in log space for separations less than ~5000 AU and follows a power law ofa−1.6at larger separations.

The presence of Jupiter is crucial to the architecture of the Solar System and models underline this to be a generic feature of planetary systems. We find the detection of the difference between the position and motion recorded by the contemporary astrometric satellite Gaia and its precursor Hipparcos can be used to discover Jupiter-like planets. We illustrate how observations of the nearby star $\varepsilon$ Indi A giving astrometric and radial velocity data can be used to independently find the orbit of its suspected companion. The radial velocity and astrometric data provide complementary detections which allow for a much stronger solution than either technique would provide individually. We quantify $\varepsilon$ Indi A b as the closest Jupiter-like exoplanet with a mass of 3 $M_{Jup}$ on a slightly eccentric orbit with an orbital period of 45 yr. While other long-period exoplanets have been discovered, $\varepsilon$ Indi A b provides a well constrained mass and along with the well-studied brown dwarf binary in orbit around $\varepsilon$ Indi A means that the system provides a benchmark case for our understanding of the formation of gas giant planets and brown dwarfs., 18 pages, 9 figures, 2 tables, accepted to be published in MNRAS. arXiv admin note: text overlap with arXiv:1803.08163

We conducted a high-precision elemental abundance analysis of the twin-star comoving pair HIP34407/HIP34426. With mean error of 0.013 dex in the differential abundances (D[X/H]), a significant difference was found: HIP34407 is more metal-rich than HIP34426. The elemental abundance differences correlate strongly with condensation temperature, with the lowest for the volatile elements like carbon around 0.05+/-0.02 dex, and the highest up to about 0.22+/-0.01 dex for the most refractory elements like aluminum. Dissimilar chemical composition for stars in twin-star comoving pairs are not uncommon, thus we compile previously-published results like ours and look for correlations between abundance differences and stellar parameters, finding no significant trends with average effective temperature, surface gravity, iron abundance, or their differences. Instead, we found a weak correlation between the absolute value of abundance difference and the projected distance between the stars in each pair that appears to be more important for elements which have a low absolute abundance. If confirmed, this correlation could be an important observational constraint for binary star system formation scenarios., MNRAS, in press

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., 10 pages, 5 figures, Accepted in ApJ Letters

Stars of a common origin are thought to have similar, if not nearly identical, chemistry. Chemical tagging seeks to exploit this fact to identify Milky Way subpopulations through their unique chemical fingerprints. In this work, we compare the chemical abundances of dwarf stars in wide binaries to test the abundance consistency of stars of a common origin. Our sample of 31 wide binaries is identified from a catalog produced by cross-matching APOGEE stars with UCAC5 astrometry, and we confirm the fidelity of this sample with precision parallaxes from Gaia DR2. For as many as 14 separate elements, we compare the abundances between components of our wide binaries, finding they have very similar chemistry (typically within 0.1 dex). This level of consistency is more similar than can be expected from stars with different origins (which show typical abundance differences of 0.3-0.4 dex within our sample). For the best measured elements, Fe, Si, K, Ca, Mn, and Ni, these differences are reduced to 0.05-0.08 dex when selecting pairs of dwarf stars with similar temperatures. Our results suggest that APOGEE dwarf stars may currently be used for chemical tagging at the level of $\sim$0.1 dex or at the level of $\sim$0.05 dex when restricting for the best-measured elements in stars of similar temperatures. Larger wide binary catalogs may provide calibration sets, in complement to open cluster samples, for on-going spectroscopic surveys., 21 pages, 14 figures, accepted for publication in ApJ

We perform a search for wide binaries in the Kepler field with the prospect of providing new constraints for gyrochronology. First, we construct our base catalog by compiling astrometry for the stars observed by Kepler, and supplement it with parallaxes, radial velocities (RVs), and metallicities. We then mine our base catalog for wide binary candidates by matching the stars' proper motions, as well as parallaxes, RVs, and metallicities, if available. We mitigate the presence of chance alignments among our candidates by performing a comprehensive data-based contamination analysis in the proper motion versus angular separation phase space. Our final sample contains 55 binary candidates. A crossmatch of our pairs with the Second Data Release (DR2) from Gaia validates our candidates and confirms the reliability of our search method, particularly for $\varpi \gtrsim 2$ mas. Due to the implicit Kepler selection function and image scale per pixel, our binary search is incomplete for angular separations of less than 20 arcsec. We crossmatch our candidates with rotation period and asteroseismic ages catalogs, and find that our binary candidates do not follow a simple period-color relation, in agreement with previous studies. Two pairs have an age estimate for one component star and rotation period for its companion, positioning them as potentially new gyrochronology constraints at old ages. This is the first study that uses RVs and metallicities as criteria, rather than as a confirmation, in a binary search., 32 pages, 17 figures; accepted for publication in MNRAS

Samples of reliably identified halo wide binaries are scarce. If reasonably free from selection effects and with a small degree of contamination by chance alignments, these wide binaries become a powerful dynamical tool, having provided one of the very few experiments capable of constraining the nature of dark matter in the Milky Way halo. Currently, however, the best available sample of halo wide binaries is confined to the solar neighborhood, and is plagued by small number statistics at the widest separations. We present the results of a program aimed to probe the wide binary population of the Galactic halo at significantly larger distances, and which informs future searches that could improve the statistics by orders of magnitude. Halo stars were taken from the Sloan Digital Sky Survey after analysing the Galactic orbits of stars in the reduced proper motion diagram. We then select candidate binaries by searching for pairs with small differences in proper motion and small projected separation on the sky. Using medium-resolution spectroscopy, a subsample of candidates is validated via radial velocities, finding that about 68\% of candidate pairs up to 20$^{\prime\prime}$ separation can be considered genuine halo wide binaries, to the limits of the available data. Precise parallaxes from \textit{Gaia} confirm that most of our selected pairs have their components at the same distances, independently confirming the robustness of our selection method. These results should prove valuable to guide the optimal assembly of larger catalogs of halo wide binaries from upcoming large databases, e.g., \textit{Gaia} and LSST., Comment: 13 pages, 18 figures. Accepted for publication in MNRAS

We address the existence and origin of the lithium (Li) desert, a region in the Li - T_eff plane sparsely populated by stars. Here we analyze some of the explanations that have been suggested for this region, including mixing in the late main sequence, a Li dip origin for stars with low Li abundances in the region, and a possible relation with the presence of planets. To study the Li desert, we measured the atmospheric parameters and Li abundance of 227 late-F dwarfs and subgiants, chosen to be in the T_eff range of the desert and without previous Li abundance measurements. Subsequently, we complemented those with literature data to obtain a homogeneous catalog of 2318 stars, for which we compute masses and ages. We characterize stars surrounding the region of the Li desert. We conclude that stars with low Li abundances below the desert are more massive and more evolved than stars above the desert. Given the unexpected presence of low Li abundance stars in this effective temperature range, we concentrate on finding their origin. We conclude that these stars with low Li abundance do not evolve from stars above the desert: at a given mass, stars with low Li (i.e., below the desert) are more metal-poor. Instead, we suggest that stars below the Li desert are consistent with having evolved from the Li dip, discarding the need to invoke additional mixing to explain this feature. Thus, stars below the Li desert are not peculiar and are only distinguished from other subgiants evolved from the Li dip in that their combination of atmospheric parameters locates them in a range of effective temperatures where otherwise only high Li abundance stars would be found (i.e., stars above the desert)., Comment: Accepted for publication in A&A

We report the peculiar chemical abundance patterns of eleven atypical Milky Way (MW) field red giant stars observed by the Apache Point Observatory Galactic Evolution Experiment (APOGEE). These atypical giants exhibit strong Al and N enhancements accompanied by C and Mg depletions, strikingly similar to those observed in the so-called second-generation (SG) stars of globular clusters (GCs). Remarkably, we find low-Mg abundances ([Mg/Fe]$, Comment: Accepted for publication in ApJ Letters, 11 pages, 2 Figures, 2 Tables

We mine the Tycho-{\it Gaia} astrometric solution (TGAS) catalog for wide stellar binaries by matching positions, proper motions, and astrometric parallaxes. We separate genuine binaries from unassociated stellar pairs through a Bayesian formulation that includes correlated uncertainties in the proper motions and parallaxes. Rather than relying on assumptions about the structure of the Galaxy, we calculate Bayesian priors and likelihoods based on the nature of Keplerian orbits and the TGAS catalog itself. We calibrate our method using radial velocity measurements and obtain 6196 high-confidence candidate wide binaries with projected separations $s\lesssim1$ pc. The normalization of this distribution suggests that at least 0.6\% of TGAS stars have an associated, distant TGAS companion in a wide binary. We demonstrate that {\it Gaia}'s astrometry is precise enough that it can detect projected orbital velocities in wide binaries with orbital periods as large as 10$^6$ yr. For pairs with $s\ \lesssim\ 4\times10^4$~AU, characterization of random alignments indicate our contamination to be $\approx$5\%. For $s \lesssim 5\times10^3$~AU, our distribution is consistent with \"{O}pik's Law. At larger separations, the distribution is steeper and consistent with a power-law $P(s)\propto s^{-1.6}$; there is no evidence in our data of any bimodality in this distribution for $s \lesssim$ 1 pc. Using radial velocities, we demonstrate that at large separations, i.e., of order $s \sim$ 1 pc and beyond, any potential sample of genuine wide binaries in TGAS cannot be easily distinguished from ionized former wide binaries, moving groups, or contamination from randomly aligned stars., Comment: 27 pages, 17 figures, 2 tables, submitted to MNRAS

In this work we use structural properties of Milky Way's outer halo ($R_G > 25\,\mathrm{kpc}$) satellites (dwarf spheroidal galaxies, ultra-faint dwarf galaxies and globular clusters) derived from deep, wide-field and homogeneous data, to present evidence of a correlation in the S\'ersic index v/s effective radius plane followed by a large fraction of outer halo globular clusters and satellite dwarf galaxies. We show that this correlation can be entirely reproduced by fitting empirical relations in the central surface brightness v/s absolute magnitude and S\'ersic index v/s absolute magnitude parameter spaces, and by assuming the existence of two types of outer halo globular clusters: one of high surface brightness (HSB group), with properties similar to inner halo clusters; and another of low surface brightness (LSB group), which share characteristics with dwarf spheroidal and ultra-faint dwarf galaxies. Given the similarities of LSB clusters with dwarf spheroidal and ultra-faint dwarf galaxies, we discuss the possibility that outer halo clusters also originated inside dark matter halos and that tidal forces from different galaxy host's potentials are responsible for the different properties between HSB and LSB clusters., Comment: 20 pages, 9 figures, 3 tables

The LMC hosts a rich variety of star clusters seen in close projected proximity. Ages have been derived for few of them showing differences up to few million years, hinting at being binary star clusters. However, final confirmation needs to be done through spectroscopic analysis. Here we focus on the LMC cluster pair NGC2006-SL538 and aim to determine whether the star cluster pair is a bound entity (binary star cluster) or a chance alignment. Using the MIKE echelle spectrograph at LCO we have acquired integrated-light spectra for each cluster. We have measured radial velocities by two methods: a) direct line profile measurement yields v$_r=300.3\pm5\pm6$ km/s for NGC2006 and $v_r=310.2\pm4\pm6$ km/s for SL538. b) By comparing observed spectra with synthetic bootstrapped spectra yielding $v_r=311.0\pm0.6$ km/s for NGC2006 and $v_r=309.4\pm0.5$ km/s for SL538. Finally when spectra are directly compared, we find a ${\Delta}v=1.08\pm0.47$ km/s. Full-spectrum SED fits reveal that the stellar population ages lie in the range 13-21 Myr with a metallicity of Z=0.008. We find indications for differences in the chemical abundance patterns as revealed by the helium absorption lines between the two clusters. The dynamical analysis shows that the two clusters are likely to merge within the next $\sim$150 Myr. The NGC2006-SL538 cluster pair shows radial velocities, stellar population and dynamical parameters consistent with a gravitational bound entity. We conclude that this is a genuine binary cluster pair, and we propose that their differences in ages and stellar population chemistry is most likely due to variances in their chemical enrichment history within their environment. We suggest that their formation may have taken place in a loosely bound star-formation complex which saw initial fragmentation but then had its clusters become a gravitationally bound pair by tidal capture., Comment: Accepted for publication in Astronomy & Astrophysics. 15 pages, 10 figures in low resolution

We present results from an observing campaign to identify low-metallicity stars in the Best & Brightest Survey. From medium-resolution (R ~ 1, 200 - 2, 000) spectroscopy of 857 candidates, we estimate the stellar atmospheric parameters (Teff, log g, and [Fe/H]), as well as carbon and alpha-element abundances. We find that 69% of the observed stars have [Fe/H], 57 pages, accepted for publication in ApJ

The Gaia-ESO survey recently reported on a large sample of lithium (Li) abundance determinations for evolved stars in the rich open cluster Trumpler 20. They argue for a scenario where virtually all stars experience post main sequence mixing and Li is preserved in only two objects. We present an alternate explanation, where Li is normal in the vast majority of cluster stars and anomalously high in these two cases. We demonstrate that the Li upper limits in the red giants can be explained with a combination of main sequence depletion and standard dredge-up, and that they are close to the detected levels in other systems of similar age. In our framework, the two detected giants are anomalously Li-rich, and we propose that both could have been produced by the engulfment of a substellar mass companion of 16(+6-10) M_J. This would imply that ~5% of 1.8 solar mass stars in this system, and by extension elsewhere, should have substellar mass companions of high mass that could be engulfed at some point in their lifetimes. We discuss future tests that could confirm or refute this scenario., 6 pages, 3 figures, accepted

We carry out controlled $N$-body simulations that follow the dynamical evolution of binary stars in the dark matter (DM) haloes of ultra-faint dwarf spheroidals (dSphs). We find that wide binaries with semi-major axes $a\gtrsim a_t$ tend to be quickly disrupted by the tidal field of the halo. In smooth potentials the truncation scale, $a_t$, is mainly governed by (i) the mass enclosed within the dwarf half-light radius ($R_h$) and (ii) the slope of the DM halo profile at $R\approx R_h$, and is largely independent of the initial eccentricity distribution of the binary systems and the anisotropy of the stellar orbits about the galactic potential. For the reported velocity dispersion and half-light radius of Segue I, the closest ultra-faint, our models predict $a_t$ values that are a factor 2--3 smaller in cuspy haloes than in any of the cored models considered here. Using mock observations of Segue I we show that measuring the projected two-point correlation function of stellar pairs with sub-arcsecond resolution may provide a useful tool to constrain the amount and distribution of DM in the smallest and most DM-dominated galaxies., Accepted to MNRAS Letters

A small fraction of red giants are known to be lithium (Li) rich, in contradiction with expectations from stellar evolutionary theory. A possible explanation for these atypical giants is the engulfment of a Li-rich planet or brown dwarf by the star. In this work, we model the evolution of Li abundance in canonical red giants including the accretion of a sub-stellar mass companion. We consider a wide range of stellar and companion masses, Li abundances, stellar metallicities, and planetary orbital periods. Based on our calculations, companions with masses lower than 15 M_J dissolve in the convective envelope and can induce Li enrichment in regimes where extra mixing does not operate. Our models indicate that the accretion of a substellar companion can explain abundances up to A(Li)~2.2, setting an upper limit for Li-rich giants formed by this mechanism. Giants with higher abundances need another mechanism to be explained. For reasonable planetary distributions, we predict the Li abundance distribution of low-mass giants undergoing planet engulfment, finding that between 1% to 3% of them should have A(Li)>1.5. We show that depending on the stellar mass range, this traditional definition of Li-rich giants is misleading, as isolated massive stars would be considered anomalous while giants engulfing a companion would be set aside, flagged as normal. We explore the detectability of companion engulfment, finding that planets with masses higher than ~7 M_J produce a distinct signature, and that descendants of stars originating in the Li-dip and low luminosity red giants are ideal tests of this channel., Accepted for publication, ApJ. 25 pages, 13 figures

We present results from a medium-resolution (R ~ 2, 000) spectroscopic follow-up campaign of 1,694 bright (V < 13.5), very metal-poor star candidates from the RAdial Velocity Experiment (RAVE). Initial selection of the low-metallicity targets was based on the stellar parameters published in RAVE Data Releases 4 and 5. Follow-up was accomplished with the Gemini-N and Gemini-S, the ESO/NTT, the KPNO/Mayall, and the SOAR telescopes. The wavelength coverage for most of the observed spectra allows for the determination of carbon and ��-element abundances, which are crucial for con- sidering the nature and frequency of the carbon-enhanced metal-poor (CEMP) stars in this sample. We find that 88% of the observed stars have [Fe/H], 134 pages, accepted for publication in AJ

Featuring greatly improved astrometry and precise radial velocities (RVs), the second data release (DR2) from Gaia affords us a unique opportunity to test the validity of samples of stellar wide binaries. In a previous work, we presented a set of wide binaries identified in the joint Tycho-Gaia Astrometric Solution catalog. Here, we use DR2 to confirm the low contamination rate of this sample, thereby verifying the effectiveness of our algorithm for application to DR2 data., Accepted for publication in Research Notes of the AAS

Establishing or ruling out, either through solid mass measurements or upper limits, the presence of intermediate-mass black holes (IMBHs) at the centers of star clusters would profoundly impact our understanding of problems ranging from the formation and long-term dynamical evolution of stellar systems, to the nature of the seeds and the growth mechanisms of supermassive black holes. While there are sound theoretical arguments both for and against their presence in today's clusters, observational studies have so far not yielded truly conclusive IMBH detections nor upper limits. We argue that the most promising approach to solving this issue is provided by the combination of measurements of the proper motions of stars at the centers of Galactic globular clusters and dynamical models able to take full advantage of this type of data set. We present a program based on HST observations and recently developed tools for dynamical analysis designed to do just that., 7 pages; Invited talk at IAU Symposium 266 "Star clusters: basic galactic building blocks" (Rio de Janeiro, 10-14 August 2009), R. de Grijs and R.D. Lepine, eds; v2 reflects a better citation of some references and corrected typos

Wide binaries, particularly in large numbers and as free from selection biases as possible, constitute a largely overlooked tool for studying the Galaxy. The goal of this review is to highlight the potential inherent to large samples of field wide binaries for research on problems as varied as star formation in the early Galaxy, the nature of halo dark matter, the evolution of the stellar halo, new geometric distances, metallicities, masses, and ages of field stars and white dwarfs, and much more. Using the Revised NLTT as an illustrative example, I review the main steps in the assembly of a large catalog of wide binaries useful for multiple applications. The capability of cleanly separating between the Galactic disk and halo populations using good colors and proper motions is emphasized. The critical role of large surveys for research on wide binaries as well as for the better understanding of the Galaxy in general is stressed throughout. Finally, I point out the potential for assembling new samples of wide binaries from available proper-motion surveys, and report on current efforts of using the SDSS towards this goal.

(ABRIDGED) We use a fully self-consistent evolutionary code to follow the rotational evolution of red giants, making a comprehensive attempt to assess the role of rotationally induced mixing in the development of abundance anomalies in giants with a range of masses and metallicities in stellar clusters and the field. Unlike most previous work, we do not focus on the determination of combinations of mixing rate and depth that reproduce the data on a particular stellar type. Instead, we concentrate on the more fundamental problem of the simultaneous reproduction of the CNO surface patterns in both Population I and Population II giants using the same physics and models. A general result of all our models is that rotational mixing, although present in small amounts, is inefficient on the lower RGB independently of any inhibiting effect of composition barriers. Models with differentially rotating envelopes are able to reproduce the carbon isotope data on M67 giants with initial rotation rates adequate to their progenitors, but fail to do so for open clusters of larger turnoff mass as well as for metal-poor giants in the field and globular clusters. Possible solutions are discussed. Our favored scenario is one in which the overall strength of canonical extra mixing has been underestimated by existent derivations, but which additionally needs to be coupled with a much lower efficiency for rotational mixing among the rapidly rotating open cluster giants than in the slowly rotating ones in the field and globular clusters. We hypothesize that this last requirement is provided by the interaction between convection and rotation in the envelopes of giants, in the sense that rapidly rotating stars would develop much shallower angular velocity profiles in their envelopes than do slowly rotating stars., Comment: 83 pages, 21 figures; replaced with version accepted for publication in the ApJ

We propose a simple test for the existence of a cluster of black hole remnants around Sgr A* that is based on a small sample of any type of Galactic Center objects, provided they are substantially less massive than the black holes and constitute part of an old (> 1 Gyr) population. The test relies on the fact that, under the presence of such a cluster of heavy remnants and because of energy equipartition, lower mass objects would be expelled from the central regions and settle into a distribution very different than the cusp expected to be induced by the supermassive black hole alone. We show that with a sample of just 50 objects and using only their angular positions on the sky relative to Sgr A* it is possible to clearly differentiate between a distribution consistent with the presence of the cluster of black holes and a power-law cusp distribution. We argue that millisecond pulsars might currently be the best candidate to perform this test, because of the large uncertainties involved in the age determination of less exotic objects. In addition, by measuring their first and second period derivatives, millisecond pulsars offer the rare opportunity of determining the complete phase space information of the objects. We show that this extra information improves the detection of mass segregation by about 30%., 16 pages, 1 Postscript figure; version accepted for publication in ApJ

Gas-rich galaxies in dense environments such as galaxy clusters and massive groups are affected by a number of possible types of interactions with the cluster environment, which make their evolution radically different than that of field galaxies. The dIrr galaxy NGC 1427A, presently infalling towards the core of the Fornax galaxy cluster, offers a unique opportunity to study those processes in a level of detail not possible to achieve for galaxies at higher redshifts. Using HST/ACS and auxiliary VLT/FORS ground-based observations, we study the properties of the most recent episodes of star formation in this gas-rich galaxy, the only one of its type near the core of the Fornax cluster. We study the structural and photometric properties of young star cluster complexes in NGC 1427A, identifying 12 bright such complexes with exceptionally blue colors. The comparison of our broadband near-UV/optical photometry with simple stellar population models yields ages below ~4x10^6 yr and stellar masses from a few thousand up to ~3x10^4 Msun, slightly dependent on the assumption of cluster metallicity and IMF. Their grouping is consistent with hierarchical and fractal star cluster formation. We use deep Ha imaging data to determine the current Star Formation Rate (SFR) in NGC 1427A and estimate the ratio, Gamma, of star formation occurring in these star cluster complexes to that in the entire galaxy. We find Gamma to be among the largest such values available in the literature, consistent with starburst galaxies. Thus a large fraction of the current star formation in NGC 1427A is occurring in star clusters, with the peculiar spatial arrangement of such complexes strongly hinting at the possibility that the starburst is being triggered by the passage of the galaxy through the cluster environment., 14 pages, 15 figures, accepted for publication in The Astronomical Journal. Updated document

We present the first study of high-precision internal proper motions (PMs) in a large sample of globular clusters, based on Hubble Space Telescope (HST) data obtained over the past decade with the ACS/WFC, ACS/HRC, and WFC3/UVIS instruments. We determine PMs for over 1.3 million stars in the central regions of 22 clusters, with a median number of ~60,000 stars per cluster. These PMs have the potential to significantly advance our understanding of the internal kinematics of globular clusters by extending past line-of-sight (LOS) velocity measurements to two- or three-dimensional velocities, lower stellar masses, and larger sample sizes. We describe the reduction pipeline that we developed to derive homogeneous PMs from the very heterogeneous archival data. We demonstrate the quality of the measurements through extensive Monte-Carlo simulations. We also discuss the PM errors introduced by various systematic effects, and the techniques that we have developed to correct or remove them to the extent possible. We provide in electronic form the catalog for NGC 7078 (M 15), which consists of 77,837 stars in the central 2.4 arcmin. We validate the catalog by comparison with existing PM measurements and LOS velocities, and use it to study the dependence of the velocity dispersion on radius, stellar magnitude (or mass) along the main sequence, and direction in the plane of the sky (radial/tangential). Subsequent papers in this series will explore a range of applications in globular-cluster science, and will also present the PM catalogs for the other sample clusters., 34 pages, 22 figures (3 in low res), 30 tables, accepted for publication in ApJ on October 20, 2014

We show that at any given time, the Galactocentric black hole Sgr A* is expected to be microlensing N(lens)=1.7 bulge stars if the threshold of detectability of the fainter image is K(thr)=21, and about N(lens)=8 sources if K(thr)=23. The lensed images then provide a unique way to detect stellar-mass black holes predicted to cluster around Sgr A*. If a black hole passes close to a microlensed image, it will give rise to a short (weeks long) microlensing event. We show that the mass and projected velocity of the stellar-mass black hole can both be measured unambiguously from such an event, provided that either a caustic crossing is observed or the astrometric displacement is measured. For K(thr)=23 and moderate magnifications by Sgr A*, the microlensing event rate from a cluster of 20000 black holes within a radius of 0.7 pc is only 0.06 events per year; however, if highly magnified images of a star were found, the rate of events by the stellar black holes would be much higher. In addition, the N(lens) sources lensed by Sgr A* provide a unique probe of extinction behind the Galactic center along 2N(lens) lines of sight.

We report the discovery of BD+20594b, a high density sub-Neptune exoplanet, made using photometry from Campaign 4 of the two-wheeled Kepler (K2) mission, ground-based radial velocity follow-up from HARPS and high resolution lucky and adaptive optics imaging obtained using AstraLux and MagAO, respectively. The host star is a bright ($V=11.04$, $K_s = 9.37$), slightly metal poor ([Fe/H]$=-0.15\pm 0.05$ dex) solar analogue located at $152.1^{+9.7}_{-7.4}$ pc from Earth, for which we find a radius of $R_*=0.928^{+0.055}_{-0.040}R_\odot$ and a mass of $M_* = 0.961^{+0.032}_{-0.029}M_\odot$. A joint analysis of the K2 photometry and HARPS radial velocities reveal that the planet is in a $\approx 42$ day orbit around its host star, has a radius of $2.23^{+0.14}_{-0.11}R_\oplus$, and a mass of $16.3^{+6.0}_{-6.1}M_\oplus$. Although the data at hand puts the planet in the region of the mass-radius diagram where we could expect planets with a pure rock (i.e. magnesium silicate) composition using two-layer models (i.e., between rock/iron and rock/ice compositions), we discuss more realistic three-layer composition models which can explain the high density of the discovered exoplanet. The fact that the planet lies in the boundary between "possibly rocky" and "non-rocky" exoplanets, makes it an interesting planet for future RV follow-up., Comment: 12 pages, 11 figures. Accepted for publication in ApJ

Schwarzschild's orbit-superposition technique is the most developed and well-tested method available for constraining the detailed mass distributions of equilibrium stellar systems. Here I provide a very short overview of the method and its existing implementations, and briefly discuss their viability as a tool for modeling the Galaxy using Gaia data., Comment: 2 pages; invited review at "Modelling the Galaxy in the era of Gaia", proceedings of Joint Discussion 5 at the XXVIIth IAU General Assembly (Rio de Janeiro, August 2009); James Binney, Ian F. Corbett, eds.

(ABRIDGED) We present a new Schwarzschild orbit-superposition code designed to model discrete datasets composed of velocities of individual kinematic tracers in a dynamical system. This constitutes an extension of previous implementations that can only address continuous data in the form of (the moments of) velocity distributions, thus avoiding potentially important losses of information due to data binning. Furthermore, the code can handle any combination of available velocity components, i.e., only line-of-sight velocities, only proper motions, or a combination of both. It can also handle a combination of discrete and continuous data. The code finds the distribution function (DF, a function of the three integrals of motion E, Lz, and I3) that best reproduces the available kinematic and photometric observations in a given axisymmetric gravitational potential. The fully numerical approach ensures considerable freedom on the form of the DF f(E,Lz,I3). This allows a very general modeling of the orbital structure, thus avoiding restrictive assumptions about the degree of (an)isotropy of the orbits. We describe the implementation of the discrete code and present a series of tests of its performance based on the modeling of simulated datasets generated from a known DF. We find that the discrete Schwarzschild code recovers the original orbital structure, M/L ratios, and inclination of the input datasets to satisfactory accuracy, as quantified by various statistics. The code will be valuable, e.g., for modeling stellar motions in Galactic globular clusters, and those of individual stars, planetary nebulae, or globular clusters in nearby galaxies. This can shed new light on the total mass distributions of these systems, with central black holes and dark matter halos being of particular interest., ApJ, in press; 51 pages, 11 figures; manuscript revised following comments by referee

The Kapteyn moving group has been postulated as tidal debris from $\omega$ Centauri. If true, members of the group should show some of the chemical abundance patterns known for stars in the cluster. We present an optical and near-infrared high-resolution, high-S/N spectroscopic study of 14 stars of the Kapteyn group, plus 10 additional stars (the $\omega$ Cen-group) that, while not listed as members of the Kapteyn group as originally defined, have been nevertheless associated dynamically with $\omega$ Centauri. Abundances for Na, O, Mg, Al, Ca and Ba were derived from the optical spectra, while the strength of the chromospheric He I 10830 {\AA} line is studied as a possible helium abundance indicator. The resulting Na-O and Mg-Al patterns for stars of the combined Kapteyn and $\omega$ Cen-group samples do not resemble those of $\omega$ Centauri, and are not different from those of field stars of the Galactic halo. The distribution of equivalent widths of the He I 10830 {\AA} line is consistent with that found among non-active field stars. Therefore, no evidence is found for second-generation stars within our samples, which most likely rules out a globular-cluster origin. Moreover, no hint of the unique Ba-overabundance at the metal-rich end, well-established for $\omega$ Centauri stars, is seen among stars of the combined samples. Because this specific Ba pattern is present in $\omega$ Centauri irrespective of stellar generation, this would rule out the possibility that our entire sample might be composed of only first generation stars from the cluster. Finally, for the stars of the Kapteyn group, the possibility of an origin in the hypothetical $\omega$ Centauri's parent galaxy is disfavored by the different run of $\alpha$-elements with metallicity between our targets and stars from present-day dwarf galaxies., Comment: 23 pages, 11 figures, ApJ, 808, 103

We present a study of the old globular cluster (GC) population of the dwarf irregular galaxy NGC 1427A using multi-wavelength VLT observations in U, B, V, I, H_alpha, J, H, and Ks bands under excellent observing conditions. We applied color and size selection criteria to select old GC candidates and made use of archival ACS images taken with the Hubble Space Telescope to reject contaminating background sources and blended objects from the GC candidates' list. The H_alpha observations were used to check for contamination due to compact, highly reddened young star clusters whose colors and sizes could mimic those of old GCs. After accounting for contamination we obtain a total number of 38+/-8 GC candidates with colors consistent with an old (~10 Gyr) and metal-poor (Z < 0.4xZ_solar) population as judged by simple stellar population models. Our contamination analysis indicates that the density distribution of GCs in the outskirts of the Fornax central cD galaxy NGC1399 may not be spherically symmetric. We derive a present-day specific frequency S_N of 1.6+/-0.23 for NGC 1427A, a value significantly larger than what is observed in the Local Group dwarf irregular galaxies and comparable with the values found for the same galaxy types in the Virgo and Fornax clusters. Assuming a universal globular cluster luminosity function turnover magnitude, we derive a distance modulus to NGC 1427A of 31.01+/-0.21 mag which places it 3.2+/-2.5(statistic)+/-1.6(systematic) Mpc in front of the Fornax central cD galaxy NGC 1399. The implications of this result for the relationship between NGC 1427A and the cluster environment are briefly discussed., 16 pages, 2 tables, 14 figures; accepted for publication in A&A

We present a catalog of 1147 candidate common proper motion binaries selected from the revised New Luyten Two-Tenths Catalog. Among these, we identify 999 genuine physical pairs using the measured proper-motion difference and the relative positions of each binary's components on a reduced proper-motion (RPM) diagram. The RPM positions also serve to classify them as either disk main-sequence (801), halo subdwarf (116), or pairs containing at least one white dwarf (82). The disk and halo samples are complete to separations of \theta=500" and \theta=900", which correspond to ~0.1 pc and ~1 pc, respectively. At wide separations, both distributions are well described by single power laws, dN/d\theta ~ \theta^{-\alpha}: \alpha=1.67+-0.07 for the disk and \alpha=1.55+-0.10 for the halo. The fact that these distributions have similar slopes (and similar normalizations as well) argues for similarity of the star-formation conditions of these two populations. The fact that the halo binaries obey a single power law out to ~1 pc permits strong constraints on halo dark-matter candidates. At somewhat closer separations (10", Comment: 56 pages, 16 figures; replaced with version accepted for publication in ApJ

We present a catalog of 424 common proper motion companions to Hipparcos stars with good (>3 sigma) parallaxes, thereby effectively providing new parallaxes for these companions. Compared to stars in the Hipparcos catalog, these stars are substantially dimmer. The catalog includes 20 WDs and an additional 29 stars with M_V>14, the great majority of the latter being M dwarfs., Submitted to ApJS, 20 pages

NGC 1427A is a LMC-like irregular galaxy in the Fornax cluster with an extended pattern of strong star formation around one of its edges, which is probably due to some kind of interaction with the cluster environment. We present H-alpha velocities within NGC 1427A, obtained through long-slit spectroscopy at seven different positions, chosen to fall on the brightest HII regions of the galaxy. Due to its location very near the center of the cluster this object is an excellent candidate to study the effects that the cluster environment has on gas-rich galaxies embedded in it. The rotation of NGC 1427A is modeled in two different ways. The global ionized gas kinematics is reasonably well described by solid-body rotation, although on small scales it shows a chaotic behaviour. In this simple model, the collision with a smaller member of the cluster as being responsible for the peculiar morphology of NGC 1427A is very unlikely, since the only candidate intruder falls smoothly into the general velocity pattern of the main galaxy. In a more elaborate model, for which we obtain a better solution, this object does not lie in the same plane of NGC 1427A, in which case we identify it as a satellite bound to the galaxy. These results are discussed in the context of a normal irregular versus one interacting with some external agent. Based on several arguments and quantitative estimates, we argue that the passage through the hot intracluster gas of the Fornax cluster is a very likely scenario to explain the morphological properties of NGC 1427A., Comment: 31 pages, LaTeX2e, uses aas2pp4.sty and psfig.sty, including 7 Postscript figures; accepted for publication in ApJ, Vol. 530, February 2000

We present some first results of stellar models of low mass, solar metallicity giants that include rotation and the mixing associated with it. We base our choice of initial conditions and input physics on observational constraints obtained for main sequence (MS) and horizontal branch (HB) stars, and obtain what can be regarded as a “best case scenario” for rotational mixing. These observations suggest local conservation of angular momentum in the cores and constant specific angular momentum in convective regions on the red giant branch (RGB). We show that under these circumstances and with moderate rotation rates at the turnoff it is possible to obtain enough mixing in order to reproduce the anomalous abundance ratios of the CNO species seen in giants. It is also shown that solid body rotation in the convective envelope fails to account for the observations. It is found that rotational mixing is inefficient on the lower RGB, so there is no need to invoke any inhibiting role of μ–gradients.

Oxygen abundances of 67 dwarf stars in the metallicity range -1.6, Comment: To appear in The Astrophysical Journal

We present a program designed to obtain age-rotation measurements of solar-type dwarfs to be used in the calibration of gyrochronology relations at ages of several Gyr. This is a region of parameter space crucial for the large-scale study of the Milky Way, and where the only constraint available today is that provided by the Sun. Our program takes advantage of a set of wide binaries selected so that one component is an evolved star and the other is a main-sequence star of FGK type. In this way, we obtain the age of the system from the evolved star, while the rotational properties of the main sequence component provides the information relevant for gyrochronology regarding the spin-down of solar-type stars. By mining currently available catalogs of wide binaries, we assemble a sample of 37 pairs well positioned for our purposes: 19 with turnoff or subgiant primaries, and 18 with white dwarf components. Using high-resolution optical spectroscopy, we measure precise stellar parameters for a subset of 15 of the pairs with turnoff/subgiant components, and use these to derive isochronal ages for the corresponding systems. Ages for 16 of the 18 pairs with white dwarf components are taken from the literature. The ages of this initial sample of 31 wide binaries range from 1 to 9 Gyr, with precisions better than 20% for almost half of these systems. When combined with measurements of the rotation period of their main sequence components, these wide binary systems would potentially provide a similar number of points useful for the calibration of gyrochronology relations at very old ages., Comment: 58 pages, 7 figures, 7 tables; replaced with version accepted for publication in ApJ

We have developed a code for modeling 3-dimensional kinematical data (i.e., both radial velocity and proper motions) of individual tracers in stellar systems. The code relies on the Schwarzschild method of orbit superposition to find the distribution function (a function of energy, angular momentum, and third integral of motion) that best reproduces all the observations (the positions and velocities of the tracers as well as the overall light distribution) in a given potential. Previous implementations of this method have generally focused on fitting integrated light measurements, instead of data for individual tracers. However, the latter is the only data that are generally available in a large class of problems (e.g., kinematics of tracers in nearby globular clusters, kinematics of planetary nebulae or globular clusters in galaxy halos, etc.). We have performed extensive tests of the fitting process by making use of data generated by a well-tested N-body simulation of an axisymmetric galaxy, quantifying how well the code can recover the known input distribution function and potential.