We have used the galaxy modeling algorithm \texttt{GravPot16}, to explore the more probable orbital elements of a sample of 64 selected N-rich stars across the Milky Way. Using the newly measured proper motions from \texttt{Gaia} DR2 with existing line-of-sight velocities from APOGEE-2 survey and spectrophotometric distance estimations from the \texttt{StarHorse}. We adopted a set of high-resolution particle simulations evolved in the same steady-state Galactic potential model with a bar, in order to identify the groups of N-rich stars that have a high probability of belonging to the bulge/bar, disk, and stellar halo component. We find that the vast majority of the N-rich stars show typically maximum height from the Galactic plane below 3 kpc, and develop rather eccentric orbits (\textit{e}$>$0.5), which means these stars appear to have bulge/bar-like and/or halo-like orbits. We also show that $\sim66$\% of the selected N-rich stars currently lives in the inner Galaxy inside the corotation radius (C.R.), whilst $\sim14$\% of the N-rich star resides in halo-like orbits. Among the N-rich in the inner Galaxy, $\sim27\%$ of them share orbital properties in the boundary between bulge/bar and disk, depending on the bar pattern speeds. Our dynamical analysis also indicates that some of the N-rich are likely halo interlopers and therefore suggest that halo contamination is not insignificant within the bulge area.

We revisit the subject of Be star candidates towards the Magellanic Clouds, previously studied by the authors using SPM4 proper motions. We obtain GAIA DR2 parallaxes and proper motions for 2357 and 994 LMC and SMC Be candidates, respectively. Parallaxes and proper motions vs. color V-I easily reveal the presence of the redder galactic contaminant foreground, as concluded in our previous work, but this time we do find a few red Be stars candidates consistent with being true Magellanic objects. A membership assessment to each Magellanic Cloud is done for each Be candidate, based on the distribution of their parallaxes and proper motions. From a compilation of published catalogues of spectroscopically confirmed Be stars, we found that 40 (LMC) and 64 (SMC) of these Be candidates, are in fact Be stars. Near infrared IRSF JHKs magnitudes were obtained for about 70% the Be stars candidates with GAIA DR2 astrometric data. Mid-infrared SAGE IRAC 3.6, 4.5, 5.8 and 8.0 $\mu$m magnitudes were obtained for about 85% as well. 6 LMC and 7 SMC confirmed Be stars show optical, near- or mid-infrared colours redder than what has been typically measured for Classical Be stars. Several of the Be candidates follow those redder-than-expected colours distributions suggesting the existence of more red Magellanic Be stars.

We explore the origins of the young B-type stars found by Casetti-Dinescu et al.(2014) at the outskirts of the Milky-Way disk in the sky region of Leading Arm of the Magellanic Stream. High-resolution spectroscopic observations made with the MIKE instrument on the Magellan Clay 6.5m telescope for nine stars are added to the previous sample analyzed by Zhang et al. (2017). We compile a sample of fifteen young stars with well-determined stellar types, ages, abundances and kinematics. With proper motions from Gaia DR2 we also derive orbits in a realistic Milky-Way potential. We find that our previous radial-velocity selected LA candidates have substantial orbital angular momentum. The substantial amount of rotational component for these stars is in contrast with the near-polar Magellanic orbit, thus rendering these stars unlikely members of the LA. There are four large orbital-energy stars in our sample. The highest orbital-energy one has an age shorter than the time to disk crossing, with a birthplace $z=2.5$~kpc and $R_{\rm GC}\sim 28$~kpc. Therefore, the origin of this star is uncertain. The remaining three stars have disk runaway origin with birthplaces between 12 and 25 kpc from the Galactic center. Also, the most energetic stars are more metal poor ([Mg/H] =$-0.50\pm0.07$) and with larger He scatter ($\sigma_{\rm [He/H]} = 0.72$) than the inner disk ones ([Mg/H] $=0.12\pm0.36$, $\sigma_{\rm [He/H]} = 0.15$). While the former group's abundance is compatible with that of the Large Magellanic Cloud, it could also reflect the metallicity gradient of the MW disk and their runaway status via different runaway mechanisms.

We explore the Gaia DR2 proper motions of six young, main-sequence stars, members of the Large Magellanic Cloud (LMC) reported by Moni Bidin et al. (2017). These stars are located in the outskirts of the disk, between 7 and 13 degrees from the LMC's center where there is very low H I content. Gaia DR2 proper motions confirm that four stars formed locally, in situ, while two are consistent with being expelled via dynamical interactions from inner, more gas-rich regions of the LMC. This finding establishes that recent star formation occurred in the periphery of the LMC, where thus far only old populations are known.

IC 166 is an intermediate-age open cluster ($\sim 1$ Gyr) which lies in the transition zone of the metallicity gradient in the outer disc. Its location, combined with our very limited knowledge of its salient features, make it an interesting object of study. We present the first high-resolution spectroscopic and precise kinematical analysis of IC 166, which lies in the outer disc with $R_{GC} \sim 12.7$ kpc. High resolution \textit{H}-band spectra were analyzed using observations from the SDSS-IV Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey. We made use of the Brussels Automatic Stellar Parameter (BACCHUS) code to provide chemical abundances based on a line-by-line approach for up to eight chemical elements (Mg, Si, Ca, Ti, Al, K, Mn and Fe). The $\alpha-$element (Mg, Si, Ca and whenever available Ti) abundances, and their trends with Fe abundances have been analysed for a total of 13 high-likelihood cluster members. No significant abundance scatter was found in any of the chemical species studied. Combining the positional, heliocentric distance, and kinematic information we derive, for the first time, the probable orbit of IC 166 within a Galactic model including a rotating boxy bar, and found that it is likely that IC 166 formed in the Galactic disc, supporting its nature as an unremarkable Galactic open cluster with an orbit bound to the Galactic plane.

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]$<$0.0) together with strong Al and N overabundances in the majority (5/7) of the metal-rich ([Fe/H]$\gtrsim - 1.0$) sample stars, which is at odds with actual observations of SG stars in Galactic CGs of similar metallicities. This chemical pattern is unique and unprecedented among MW stars, posing urgent questions about its origin. These atypical stars could be former SG stars of dissolved GCs formed with intrinsically lower abundances of Mg and enriched Al (subsequently self-polluted by massive AGB stars) or the result of exotic binary systems. We speculate that the stars Mg-deficiency as well as the orbital properties suggest that they could have an extragalactic origin. This discovery should guide future dedicated spectroscopic searches of atypical stellar chemical patterns in our Galaxy; a fundamental step forward to understand the Galactic formation and evolution.

We present a proper motion investigation of a sample of Be star candidates towards the Magellanic Clouds, which has resulted in the identification of separate populations, in the Galactic foreground and in the Magellanic background. Be stars are broadly speaking B-type stars that have shown emission lines in their spectra. In this work, we studied a sample of 2446 and 1019 Be star candidates towards the LMC and SMC respectively, taken from the literature and proposed as possible Be stars due to their variability behaviour in the OGLE-II I band. JHKs magnitudes from the IRSF catalog and proper motions from the SPM4 catalog, were obtained for 1188 and 619 LMC and SMC Be stars candidates, respectively. Color-color and vector-point diagrams were used to identify different populations among the Be star candidates. In the LMC sample, two populations with distinctive infrared colours and kinematics were found, the bluer sample is consistent with being in the LMC and the redder one with belonging to the Milky Way disk. This settles the nature of the redder sample which had been described in previous publications as a possible unknown subclass of stars among the Be candidates in the LMC. In the SMC sample, a similar but less evident result was obtained, since this apparent unknown subclass was not seen in this galaxy. We confirm that in the selection of Be stars by their variability, although generally successful, there is a higher risk of contamination by Milky Way objects towards redder B$-$V and V$-$I colors.

Chemical abundances of eight O- and B-type stars are determined from high-resolution spectra obtained with the MIKE instrument on the Magellan 6.5m Clay telescope. The sample is selected from 42 candidates of membership in the Leading Arm of the Magellanic System. Stellar parameters are measured by two independent grids of model atmospheres and analysis procedures, confirming the consistency of the stellar parameter results. Abundances of seven elements (He, C, N, O, Mg, Si, and S) are determined for the stars, as are their radial velocities and estimates of distances and ages. Among the seven B-type stars analyzed, the five that have radial velocities compatible with membership to the LA have an average [Mg/H] of $-0.42\pm0.16$, significantly lower than the average of the remaining two [Mg/H] = $-0.07\pm0.06$ that are kinematical members of the Galactic disk. Among the five LA members, four have individual [Mg/H] abundance compatible with that in the LMC. Within errors, we can not exclude the possibility that one of these stars has a [Mg/H] consistent with the more metal-poor, SMC-like material. The remaining fifth star has a [Mg/H] close to MW values. Distances to the LA members indicate that they are at the edge of the Galactic disk, while ages are of the order of $\sim 50-70$ Myr, lower than the dynamical age of the LA, suggesting a single star-forming episode in the LA. V$_{\rm LSR}$ the LA members decreases with decreasing Magellanic longitude, confirming the results of previous LA gas studies.

Despite their close proximity, the complex interplay between the two Magellanic Clouds, the Milky Way, and the resulting tidal features, is still poorly understood. Recent studies have shown that the Large Magellanic Cloud (LMC) has a very extended disk strikingly perturbed in its outskirts. We search for recent star formation in the far outskirts of the LMC, out to ~30 degrees from its center. We have collected intermediate-resolution spectra of thirty-one young star candidates in the periphery of the LMC and measured their radial velocity, stellar parameters, distance and age. Our measurements confirm membership to the LMC of six targets, for which the radial velocity and distance values match well those of the Cloud. These objects are all young (10-50 Myr), main-sequence stars projected between 7 and 13 degrees from the center of the parent galaxy. We compare the velocities of our stars with those of a disk model, and find that our stars have low to moderate velocity differences with the disk model predictions, indicating that they were formed in situ. Our study demonstrates that recent star formation occurred in the far periphery of the LMC, where thus far only old objects were known. The spatial configuration of these newly-formed stars appears ring-like with a radius of 12 kpc, and a displacement of 2.6 kpc from the LMC's center. This structure, if real, would be suggestive of a star-formation episode triggered by an off-center collision between the Small Magellanic Cloud and the LMC's disk.

We report on detection, from observations obtained with the APOGEE spectroscopic survey, of a metal-poor ([Fe/H] $= -1.3$ dex) field giant star with an extreme Mg-Al abundance ratio ([Mg/Fe] $= -0.31$ dex; [Al/Fe] $= 1.49$ dex). Such low Mg/Al ratios are seen only among the second-generation population of globular clusters, and are not present among Galactic disk field stars. The light element abundances of this star, 2M16011638-1201525, suggest that it could have been born in a globular cluster. We explore several origin scenarios, in particular studying the orbit of the star to check the probability of it being kinematically related to known globular clusters. We performed simple orbital integrations assuming the estimated distance of 2M16011638-1201525 and the available six-dimensional phase-space coordinates of 63 globular clusters, looking for close encounters in the past with a minimum distance approach within the tidal radius of each cluster. We found a very low probability that 2M16011638-1201525 was ejected from most globular clusters; however, we note that the best progenitor candidate to host this star is globular cluster $\omega$ Centauri (NGC 5139). Our dynamical investigation demonstrates that 2M16011638-1201525 reaches a distance $|Z_{max}| < 3 $ kpc from the Galactic plane and a minimum and maximum approach to the Galactic center of $R_{min}<0.62$ kpc and $R_{max}<7.26$ kpc in an eccentric ($e\sim0.53$) and retrograde orbit. Since the extreme chemical anomaly of 2M16011638-1201525 has also been observed in halo field stars, this object could also be considered a halo contaminant, likely been ejected into the Milky Way disk from the halo. We conclude that, 2M16011638-20152 is also kinematically consistent with the disk but chemically consistent with halo field stars.

The most accurate 6D phase-space information from the Radial Velocity Experiment (RAVE) was used to integrate the orbits of 105 stars around the galactic globular cluster 47 Tucanae, to look for close encounters between them in the past, with a minimum distance approach less than the cluster tidal radius. The stars are currently over the distance range 3.0 kpc $<$ d $<$ 5.5 kpc. Using the uncertainties in the current position and velocity vector for both, star and cluster, 105 pairs of star-cluster orbits were generated in a Monte Carlo numerical scheme, integrated over 2 Gyr and considering an axisymmetric and non-axisymmetric Milky-Way-like Galactic potential, respectively. In this scheme, we identified 20 potential cluster members that had close encounters with the globular cluster 47 Tucanae, all of which have a relative velocity distribution (V$_{rel}$) less than 200 km s$^{-1}$ at the minimum distance approach. Among these potential members, 9 had close encounters with the cluster with velocities less than the escape velocity of 47 Tucanae, therefore a scenario of tidal stripping seems likely. These stars have been classified with a 93\% confidence level, leading to the identification of extratidal cluster stars. For the other 11 stars, V$_{rel}$ exceeds the escape velocity of the cluster, therefore they were likely ejected or are unassociated interlopers.

We present preliminary results of the CIDA Equatorial Variability Survey (CEVS), looking for quasar (hereafter QSO) candidates near the Galactic plane. The CEVS contains photometric data from extended and adjacent regions of the Milky Way disk ($\sim$ 500 sq. deg.). In this work 2.5 square degrees with moderately high temporal sampling in the CEVS were analyzed. The selection of QSO candidates was based on the study of intrinsic optical photometric variability of 14,719 light curves. We studied samples defined by cuts in the variability index (Vindex $>$ 66.5), periodicity index (Q $>$ 2), and the distribution of these sources in the plane (AT , ${\gamma}$), using a slight modification of the first-order of the structure function for the temporal sampling of the survey. Finally, 288 sources were selected as QSO candidates. The results shown in this work are a first attempt to develop a robust method to detect QSO towards the Galactic plane in the era of massive surveys such as VISTA and Gaia.

Using six-dimesional phase-space information from the Fourth Data release of the Radial Velocity Experiment (RAVE) over the range of Galactic longitude 240$^{\circ}< l <$ 360$^{\circ}$ and $V_{LSR} < -239$ kms$^{-1}$, we have computed orbits for 329 RAVE stars that were originally selected as chemically and kinematically related to $\omega$ Centauri. The orbits were integrated in a Milky-Way-like axisymmetric Galactic potential, ignoring the effects of the dynamical evolution of $\omega$ Centauri due to the tidal effects of the Galaxy disk on the cluster along time. We also ignored secular changes in the Milky Way potential over time. In a Monte Carlo scheme, and under the assumption that the stars may have been ejected with velocities greater than the escape velocity ($V_{rel}>V_{esc,0}$) from the cluster, we identified 15 stars as having close encounters with $\omega$ Centauri: (\textit{i}) 8 stars with relative velocities $V_{rel}< 200 $ kms$^{-1}$ may have been ejected $\sim$ 200 Myr ago from $\omega$ Centauri; (\textit{ii}) other group of 7 stars were identified with high relative velocity $V_{rel}> 200 $ kms$^{-1}$ during close encounters, and seems unlikely that they have been ejected from $\omega$ Centauri. We also confirm the link between J131340.4-484714 as potential member of $\omega$ Centauri, and probably ejected $\sim$ 2.0 Myr ago, with a relative velocity $V_{rel}\sim80$ kms$^{-1}$.

Strongly interacting galaxies undergo a short-lived but dramatic phase of evolution characterized by enhanced star formation, tidal tails, bridges and other morphological peculiarities. The nearest example of a pair of interacting galaxies is the Magellanic Clouds, whose dynamical interaction produced the gaseous features known as the Magellanic Stream trailing the pair's orbit about the Galaxy, the Bridge between the Clouds, and the Leading Arm, a wide and irregular feature leading the orbit. Young, newly formed stars in the Bridge are known to exist, giving witness to the recent interaction between the Clouds. However, the interaction of the Clouds with the Milky Way is less well understood. In particular, the Leading Arm must have a tidal origin, however no purely gravitational model is able to reproduce its morphology and kinematics. A hydrodynamical interaction with the gaseous hot halo and disk of the Galaxy is plausible as suggested by some models and supporting neutral hydrogen observations. Here we show for the first time that young, recently formed stars exist in the Leading Arm, indicating that the interaction between the Clouds and our Galaxy is strong enough to trigger star formation in certain regions of the Leading Arm --- regions in the outskirts of the Milky Way disk (R ~ 18 kpc), far away from the Clouds and the Bridge.

We present our spectroscopic program aimed to study some new interesting features recently discovered in the Magellanic Cloud System. These were revealed by the spatial distribution of OB-type candidate stars, selected based on UV, optical and IR photometry and proper motions from existing large-area catalogs. As a pilot study of our project, we are studying OB-star candidates in the Leading Arm (LA) of the Magellanic Stream, a gaseous tidal structure with no stellar counterpart known so far. Our targets group in three clumps near regions of high HI density in the LA. If confirmed, these young stars would evidence recent star formation in the LA, and they would help better understand and constrain the formation of the LA and its interactions with the Milky Way.

Young, OB-type candidates are identified in a ~7900 sq-deg. region encompassing the Large and Small Magellanic Clouds (LMC/SMC) periphery, the Bridge, part of the Magellanic Stream (MS) and Leading Arm (LA). Selection is based on UV, optical and IR photometry from existing large-area surveys and proper motions from the Southern Proper Motion 4 catalog (SPM4). The spatial distribution of these young star candidates shows: 1) a well-populated SMC wing which continues westward with two branches partially surrounding the SMC, 2) a rather narrow path from the SMC wing eastward toward the LMC which is offset by 1 to 2 deg. from the high-density H I ridge in the Bridge, 3) a well-populated periphery of the LMC dominated by clumps of stars at the ends of the LMC bar and 4) a few scattered candidates in the MS and two overdensities in the LA regions above and below the Galactic plane. Additionally, a proper-motion analysis is made of a radial-velocity selected sample of red giants and supergiants in the LMC, previously shown to be a kinematically and chemically distinct subgroup, most likely captured from the SMC. SPM4 proper motions of these stars also indicate they are distinct from the LMC population. The observational results presented here, combined with the known orbits of the Clouds, and other aspects of the LMC morphology, suggest an off-center, moderate to highly-inclined collision between the SMC and the LMC's disk that took place between 100 and 200 Myr ago.

We present the fourth installment of the Yale/San Juan Southern Proper Motion Catalog, SPM4. The SPM4 contains absolute proper motions, celestial coordinates, and (B,V) photometry for over 103 million stars and galaxies between the south celestial pole and -20 deg declination. The catalog is roughly complete to V=17.5 and is based on photographic and CCD observations taken with the Yale Southern Observatory's double-astrograph at Cesco Observatory in El Leoncito, Argentina. The proper-motion precision, for well-measured stars, is estimated to be 2 to 3 mas/yr, depending on the type of second-epoch material. At the bright end, proper motions are on the International Celestial Reference System by way of Hipparcos Catalog stars, while the faint end is anchored to the inertial system using external galaxies. Systematic uncertainties in the absolute proper motions are on the order of 1 mas/yr.

Absolute proper motions are determined for stars and galaxies to V=17.5 over a 450 square-degree area that encloses both Magellanic Clouds. The proper motions are based on photographic and CCD observations of the Yale/San Juan Southern Proper Motion program, which span over a baseline of 40 years. Multiple, local relative proper motion measures are combined in an overlap solution using photometrically selected Galactic Disk stars to define a global relative system that is then transformed to absolute using external galaxies and Hipparcos stars to tie into the ICRS. The resulting catalog of 1.4 million objects is used to derive the mean absolute proper motions of the Large Magellanic Cloud and the Small Magellanic Cloud; $(\mu_\alpha\cos\delta,\mu_\delta)_{LMC}=(1.89,+0.39)\pm (0.27,0.27)\;\;\{mas yr}^{-1}$ and $(\mu_\alpha\cos\delta,\mu_\delta)_{SMC}=(0.98,-1.01)\pm (0.30,0.29)\;\;\{mas yr}^{-1}$. These mean motions are based on best-measured samples of 3822 LMC stars and 964 SMC stars. A dominant portion (0.25 mas yr$^{-1}$) of the formal errors is due to the estimated uncertainty in the inertial system of the Hipparcos Catalog stars used to anchor the bright end of our proper motion measures. A more precise determination can be made for the proper motion of the SMC {\it relative} to the LMC; $(\mu_{\alpha\cos\delta},\mu_\delta)_{SMC-LMC} = (-0.91,-1.49) \pm (0.16,0.15)\;\;\{mas yr}^{-1}$. This differential value is combined with measurements of the proper motion of the LMC taken from the literature to produce new absolute proper-motion determinations for the SMC, as well as an estimate of the total velocity difference of the two clouds to within $\pm$54 kms$^{-1}$.

A proper motion study of a field of 20' x 20' inside Plaut's low extinction window (l,b)=(0 deg,-8 deg), has been completed. Relative proper motions and photographic BV photometry have been derived for ~21,000 stars reaching to V~20.5 mag, based on the astrometric reduction of 43 photographic plates, spanning over 21 years of epoch difference. Proper motion errors are typically 1 mas/yr and field dependent systematics are below 0.2 mas/yr. Cross-referencing with the 2MASS catalog yielded a sample of ~8,700 stars, from which predominantly disk and bulge subsamples were selected photometrically from the JH color-magnitude diagram. The two samples exhibited different proper-motion distributions, with the disk displaying the expected reflex solar motion as a function of magnitude. Galactic rotation was also detected for stars between ~2 and ~3 kpc from us. The bulge sample, represented by red giants, has an intrinsic proper motion dispersion of (sigma_l,sigma_b)=(3.39, 2.91)+/-(0.11,0.09) mas/yr, which is in good agreement with previous results, and indicates a velocity anisotropy consistent with either rotational broadening or tri-axiality. A mean distance of 6.37^{+0.87}_{-0.77} kpc has been estimated for the bulge sample, based on the observed K magnitude of the horizontal branch red clump. The metallicity [M/H] distribution was also obtained for a subsample of 60 bulge giants stars, based on calibrated photometric indices. The observed [M/H] shows a peak value at [M/H]~-0.1 with an extended metal poor tail and around 30% of the stars with supersolar metallicity. No change in proper motion dispersion was observed as a function of [M/H]. We are currently in the process of obtaining CCD UBVRI photometry for the entire proper-motion sample of ~21,000 stars.

Two motions in the galactic plane have been detected and characterized, based on the determination of a common systematic component in Hipparcos catalogue proper motions. The procedure is based only on positions, proper motions and parallaxes, plus a special algorithm which is able to reveal systematic trends. Our results come from two stellar samples. Sample 1 has 4566 stars and defines a motion of apex (l,b)=(177.8,3.7)+/-(1.5,1.0) and space velocity V=27+/-1 km/s. Sample 2 has 4083 stars and defines a motion of apex (l,b)=(5.4,-0.6)+/-(1.9,1.1) and space velocity V=32+/-2 km/s. Both groups are distributed all over the sky and cover a large variety of spectral types, which means that they do not belong to a specific stellar population. Herschel's method is used to define the initial samples of stars and later to compute the common space velocity. The intermediate process is based on the use of a special algorithm to determine systematic components in the proper motions. As an important contribution, this paper sets out a new way to study the kinematics of the solar neighborhood, in the search for streams, associations, clusters and any other space motion shared by a large number of stars, without being restricted by the availability of radial velocities.

We have measured the periods and light curves of 148 RR Lyrae variables from V=13.5 to 19.7 from the first 100 sq. degrees of the QUEST RR Lyrae survey. Approximately 55% of these stars belong to the clump of stars detected earlier by the Sloan Digital Sky Survey. According to our measurements, this feature has ~10 times the background density of halo stars, spans at least 37.5 deg by 3.5 deg in right ascension and declination (>=30 by >=3 kpc), lies ~50 kpc from the Sun, and has a depth along the line of sight of ~5 kpc (1 sigma). These properties are consistent with the recent models that suggest it is a tidal stream from the Sgr dSph galaxy. The mean period of the type ab variables, 0.58 d, is also consistent. In addition, we have found two smaller over-densities in the halo, one of which may be related to the globular cluster Pal 5.

We describe the discovery circumstances and photometric properties of 2000 EB173, now one of the brightest trans-Neptunian objects (TNOs) with opposition magnitude m_R=18.9 and also one of the largest Plutinos, found with the drift-scanning camera of the QUEST Collaboration, attached to the 1-m Schmidt telescope of the National Observatory of Venezuela. We measure B-V = 0.99 +/- 0.14 and V-R = 0.57 +/- 0.05, a red color observed for many fainter TNOs. At our magnitude limit m_R = 20.1 +/- 0.20, our single detection reveals a sky density of 0.015 (+0.034, -0.012) TNOs per deg^2 (the error bars are 68% confidence limits), consistent with fainter surveys showing a cumulative number proportional to 10^0.5m_R. Assuming an inclination distribution of TNOs with FWHM exceeding 30 deg, it is likely that one hundred to several hundred objects brighter than m_R=20.1 remain to be discovered.