Several dozen optical echelle spectra demonstrate that HR 6819 is a hierarchical triple. A classical Be star is in a wide orbit with an unconstrained period around an inner 40 d binary consisting of a B3 III star and an unseen companion in a circular orbit. The radial-velocity semi-amplitude of 61.3 km/s of the inner star and its minimum (probable) mass of 5.0 Msun (6.3 +- 0.7 Msun) imply a mass of the unseen object of >= 4.2 Msun (>= 5.0 +- 0.4 Msun), that is, a black hole (BH). The spectroscopic time series is stunningly similar to observations of LB-1. A similar triple-star architecture of LB-1 would reduce the mass of the BH in LB-1 from ~70 Msun to a level more typical of Galactic stellar remnant BHs. The BH in HR 6819 probably is the closest known BH to the Sun, and together with LB-1, suggests a population of quiet BHs. Its embedment in a hierarchical triple structure may be of interest for models of merging double BHs or BH + neutron star binaries. Other triple stars with an outer Be star but without BH are identified; through stripping, such systems may become a source of single Be stars.

Rapid rotation is a fundamental characteristic of classical Be stars and a crucial property allowing for the formation of their circumstellar disks. Past evolution in a mass and angular momentum transferring binary system offers a plausible solution to how Be stars attained their fast rotation. Although the subdwarf remnants of mass donors in such systems should exist in abundance, only a few have been confirmed due to tight observational constraints. An indirect method of detecting otherwise hidden companions is offered by their effect on the outer parts of Be star disks, which are expected to be disrupted or truncated. In the context of the IR and radio continuum excess radiation originating in the disk, the disk truncation can be revealed by a turndown in the spectral energy distribution due to reduced radio flux levels. In this work we search for signs of spectral turndown in a sample of 57 classical Be stars with radio data, which include new data for 23 stars and the longest wavelength detections so far ($\lambda \approx$ 10\,cm) for 2 stars. We confidently detect the turndown for all 26 stars with sufficient data coverage (20 of which are not known to have close binary companions). For the remaining 31 stars, data are inconclusive as to whether the turndown is present or not. The analysis suggests that many if not all Be stars have close companions influencing their outer disks. If confirmed to be subdwarf companions, the mass transfer spin-up scenario might explain the existence of the vast majority of classical Be stars.

We revisit the viscous decretion disk (VDD) model of the classical Be star $\beta$ CMi as presented by \citet{klement} using an updated version of the radiative transfer code {\ttfamily HDUST}. A software bug was causing the mean intensities to be slightly underestimated in the equatorial region of the disk, with small but detectable effects on the disk temperature and emergent spectrum. The new model fixes an unexplained feature of the original model, which was able to reproduce the observations only when considering a dual density behavior: a steep density fall-off in the very inner parts of the disk followed by a shallower density profile. The new model is able to reproduce all the observables reasonably well using a single power-law for the density profile throughout the whole disk, as predicted by the VDD model. All the other original conclusions, most importantly the reported truncation of the disk at the distance of 35 stellar equatorial radii from the central star, remain unchanged.

Context. The structure of the inner parts of Be star disks (20 stellar radii) is well explained by the viscous decretion disk (VDD) model, which is able to reproduce the observable properties of most of the objects studied so far. The outer parts, on the ther hand, are not observationally well-explored, as they are observable only at radio wavelengths. A steepening of the spectral slope somewhere between infrared and radio wavelengths was reported for several Be stars that were previously detected in the radio, but a convincing physical explanation for this trend has not yet been provided. Aims. We test the VDD model predictions for the extended parts of a sample of six Be disks that have been observed in the radio to address the question of whether the observed turndown in the spectral energy distribution (SED) can be explained in the framework of the VDD model, including recent theoretical development for truncated Be disks in binary systems. Methods. We combine new multi-wavelength radio observations from the Karl. G. Jansky Very Large Array (JVLA) and Atacama Pathfinder Experiment (APEX) with previously published radio data and archival SED measurements at ultraviolet, visual, and infrared wavelengths. The density structure of the disks, including their outer parts, is constrained by radiative transfer modeling of the observed spectrum using VDD model predictions. In the VDD model we include the presumed effects of possible tidal influence from faint binary companions. Results. For 5 out of 6 studied stars, the observed SED shows strong signs of SED turndown between far-IR and radio wavelengths. A VDD model that extends to large distances closely reproduces the observed SEDs up to far IR wavelengths, but fails to reproduce the radio SED. ... (abstract continues but did not fit here)

The late-type Be star $\beta$ CMi is remarkably stable compared to other Be stars that have been studied. This has led to a realistic model of the outflowing Be disk by Klement et al. These results showed that the disk is likely truncated at a finite radius from the star, which Klement et al.~suggest is evidence for an unseen binary companion in orbit. Here we report on an analysis of the Ritter Observatory spectroscopic archive of $\beta$ CMi to search for evidence of the elusive companion. We detect periodic Doppler shifts in the wings of the H$\alpha$ line with a period of 170 d and an amplitude of 2.25 km s$^{-1}$, consistent with a low-mass binary companion ($M\approx 0.42 M_\odot$). We then compared the small changes in the violet-to-red peak height changes ($V/R$) with the orbital motion. We find weak evidence that it does follow the orbital motion, as suggested by recent Be binary models by Panoglou et al. Our results, which are similar to those for several other Be stars, suggest that $\beta$ CMi may be a product of binary evolution where Roche lobe overflow has spun up the current Be star, likely leaving a hot subdwarf or white dwarf in orbit around the star. Unfortunately, no direct sign of this companion star is found in the very limited archive of {\it International Ultraviolet Explorer} spectra.

A global disk oscillation implemented in the viscous decretion disk (VDD) model has been used to reproduce most of the observed properties of the well known Be star $\zeta$ Tau. 48 Librae shares several similarities with $\zeta$ Tau -- they are both early-type Be stars, they display shell characteristics in their spectra, and they exhibit cyclic $V/R$ variations -- but has some marked differences as well, such as a much denser and more extended disk, a much longer $V/R$ cycle, and the absence of the so-called triple-peak features. We aim to reproduce the photometric, polarimetric, and spectroscopic observables of 48 Librae with a self-consistent model, and to test the global oscillation scenario for this target. Our calculations are carried out with the three-dimensional NLTE radiative transfer code HDUST. We employ a rotationally deformed, gravity-darkened central star, surrounded by a disk whose unperturbed state is given by the VDD model. A two-dimensional global oscillation code is then used to calculate the disk perturbation, and superimpose it on the unperturbed disk. A very good, self-consistent fit to the time-averaged properties of the disk is obtained with the VDD. The calculated perturbation has a period $P = 12$ yr, which agrees with the observed period, and the behaviour of the $V/R$ cycle is well reproduced by the perturbed model. The perturbed model improves the fit to the photometric data and reproduces some features of the observed spectroscopic data. Some suggestions to improve the synthesized spectroscopy in a future work are given.

The viscous decretion disk (VDD) model is able to explain most of the currently observable properties of the circumstellar disks of Be stars. However, more stringent tests, focusing on reproducing multitechnique observations of individual targets via physical modeling, are needed to study the predictions of the VDD model under specific circumstances. In the case of nearby, bright Be star $\beta$ CMi, these circumstances are a very stable low-density disk and a late-type (B8Ve) central star. The aim is to test the VDD model thoroughly, exploiting the full diagnostic potential of individual types of observations, in particular, to constrain the poorly known structure of the outer disk if possible, and to test truncation effects caused by a possible binary companion using radio observations. We use the Monte Carlo radiative transfer code HDUST to produce model observables, which we compare with a very large set of multitechnique and multiwavelength observations that include ultraviolet and optical spectra, photometry covering the interval between optical and radio wavelengths, optical polarimetry, and optical and near-IR (spectro)interferometry. Due to the absence of large scale variability, data from different epochs can be combined into a single dataset. A parametric VDD model with radial density exponent of $n$ = 3.5, which is the canonical value for isothermal flaring disks, is found to explain observables typically formed in the inner disk, while observables originating in the more extended parts favor a shallower, $n$ = 3.0, density falloff. Modeling of radio observations allowed for the first determination of the physical extent of a Be disk (35$^{+10}_{-5}$ stellar radii), which might be caused by a binary companion. Finally, polarization data allowed for an indirect measurement of the rotation rate of the star, which was found to be $W \gtrsim 0.98$, i.e., very close to critical.

Aims. We reanalyze FEROS observations of the star HIP 11952 to reassess the existence of the proposed planetary system. Methods. The radial velocity of the spectra were measured by cross-correlating the observed spectrum with a synthetic template. We also analyzed a large dataset of FEROS and HARPS archival data of the calibrator HD 10700 spanning over more than five years. We compared the barycentric velocities computed by the FEROS and HARPS pipelines. Results. The barycentric correction of the FEROS-DRS pipeline was found to be inaccurate and to introduce an artificial one-year period with a semi-amplitude of 62 m/s. Thus the reanalysis of the FEROS data does not support the existence of planets around HIP 11952.

Aims. We carried out a radial-velocity survey to search for planets around metal-poor stars. In this paper we report the discovery of two planets around HIP 11952, a metal-poor star with [Fe/H]= -1.9 that belongs to our target sample. Methods. Radial velocity variations of HIP 11952 were monitored systematically with FEROS at the 2.2 m telescope located at the ESO La Silla observatory from August 2009 until January 2011. We used a cross-correlation technique to measure the stellar radial velocities (RV). Results. We detected a long-period RV variation of 290 d and a short-period one of 6.95 d. The spectroscopic analysis of the stellar activity reveals a stellar rotation period of 4.8 d. The Hipparcos photometry data shows intra-day variabilities, which give evidence for stellar pulsations. Based on our analysis, the observed RV variations are most likely caused by the presence of unseen planetary companions. Assuming a primary mass of 0.83 M\odot, we computed minimum planetary masses of 0.78 MJup for the inner and 2.93 MJup for the outer planet. The semi-major axes are a1 = 0.07 AU and a2 = 0.81 AU, respectively. Conclusions. HIP 11952 is one of very few stars with [Fe/H]< -1.0 which have planetary companions. This discovery is important to understand planet formation around metal-poor stars

The Pan-STARRS1 survey is obtaining multi-epoch imaging in 5 bands (gps rps ips zps yps) over the entire sky North of declination -30deg. We describe here the implementation of the Photometric Classification Server (PCS) for Pan-STARRS1. PCS will allow the automatic classification of objects into star/galaxy/quasar classes based on colors, the measurement of photometric redshifts for extragalactic objects, and constrain stellar parameters for stellar objects, working at the catalog level. We present tests of the system based on high signal-to-noise photometry derived from the Medium Deep Fields of Pan-STARRS1, using available spectroscopic surveys as training and/or verification sets. We show that the Pan-STARRS1 photometry delivers classifications and photometric redshifts as good as the Sloan Digital Sky Survey (SDSS) photometry to the same magnitude limits. In particular, our preliminary results, based on this relatively limited dataset down to the SDSS spectroscopic limits and therefore potentially improvable, show that stars are correctly classified as such in 85% of cases, galaxies in 97% and QSOs in 84%. False positives are less than 1% for galaxies, ~19% for stars and ~28% QSOs. Moreover, photometric redshifts for 1000 luminous red galaxies up to redshift 0.5 are determined to 2.4% precision with just 0.4% catastrophic outliers and small (-0.5%) residual bias. PCS will create a value added catalog with classifications and photometric redshifts for eventually many millions sources.

We study the kinematics of the Galactic thin and thick disk populations using stars from the RAVE survey's second data release together with distance estimates from Breddels et al. (2009). The velocity distribution exhibits the expected moving groups present in the solar neighborhood. We separate thick and thin disk stars by applying the X (stellar-population) criterion of Schuster et al. (1993), which takes into account both kinematic and metallicity information. For 1906 thin disk and 110 thick disk stars classified in this way, we find a vertical velocity dispersion, mean rotational velocity and mean orbital eccentricity of (sigma_W, Vphi, e)_thin = (18\pm0.3 km/s, 223\pm0.4 km/s, 0.07\pm0.07) and (sigma_W, Vphi, e)_thick = (35\pm2 km/s, 163\pm2 km/s, 0.31\pm0.16), respectively. From the radial Jeans equation, we derive a thick disk scale length in the range 1.5-2.2 kpc, whose greatest uncertainty lies in the adopted form of the underlying potential. The shape of the orbital eccentricity distribution indicates that the thick disk stars in our sample most likely formed in situ with minor gas-rich mergers and/or radial migration being the most likely cause for their orbits. We further obtain mean metal abundances of <[M/H]>_thin = +0.03 \pm 0.17, and <[M/H]>_thick = -0.51\pm0.23, in good agreement with previous estimates. We estimate a radial metallicity gradient in the thin disk of -0.07 dex/kpc, which is larger than predicted by chemical evolution models where the disk grows insideout from infalling gas. It is, however, consistent with models where significant migration of stars shapes the chemical signature of the disk, implying that radial migration might play at least part of a role in the thick disk's formation.

We report the detection of a planetary companion around HIP 13044, a metal-poor star on the red Horizontal Branch. The detection is based on radial velocity observations with FEROS, a high-resolution spectrograph at the 2.2-m MPG/ESO telescope, located at ESO La Silla observatory in Chile. The periodic radial velocity variation of P = 16.2 days can be distinguished from the periods of the stellar activity indicators. We computed a minimum planetary mass of 1.25 MJup and an orbital semi-major axis of 0.116 AU for the planet. This discovery is unique in three aspects: First, it is the first planet detection around a star with a metallicity much lower than few percent of the solar value; second, the planet host star resides in a stellar evolutionary stage that is still unexplored in the exoplanet surveys; third, the star HIP 13044 belongs to one of the most significant stellar halo streams in the solar neighborhood, implying an extragalactic origin of the planetary system HIP 13044 in a disrupted former satellite of the Milky Way.

Stars in their late stage of evolution, such as Horizontal Branch stars, are still largely unexplored for planets. We report the detection of a planetary companion around HIP 13044, a very metal-poor star on the red Horizontal Branch, based on radial velocity observations with a high-resolution spectrograph at the 2.2-m MPG/ESO telescope. The star's periodic radial velocity variation of P=16.2 days caused by the planet can be distinguished from the periods of the stellar activity indicators. The minimum mass of the planet is 1.25 Jupiter masses and its orbital semi-major axis 0.116 AU. Because HIP 13044 belongs to a group of stars that have been accreted from a disrupted satellite galaxy of the Milky Way, the planet most likely has an extragalactic origin.

We report the detection of a planetary companion around HIP 13044, a metal-poor red horizontal branch star belonging to a stellar halo stream that results from the disruption of an ancient Milky Way satellite galaxy. The detection is based on radial velocity observations with FEROS at the 2.2-m MPG/ESO telescope. The periodic radial velocity variation of P=16.2 days can be distinguished from the periods of the stellar activity indicators. We computed a minimum planetary mass of 1.25 Jupiter masses and an orbital semimajor axis of 0.116 AU for the planet. This discovery is unique in three aspects: First, it is the first planet detection around a star with a metallicity much lower than few percent of the solar value; second, the planet host star resides in a stellar evolutionary stage that is still unexplored in the exoplanet surveys; third, the planetary system HIP 13044 most likely has an extragalactic origin in a disrupted former satellite of the Milky Way.

An efficient separation between dwarfs and giants in surveys of bright stars is important, especially for studies in which distances are estimated through photometric parallax relations. We use the available spectroscopic log g estimates from the second RAVE data release (DR2) to assign each star a probability for being a dwarf or subgiant/giant based on mixture model fits to the log g distribution in different color bins. We further attempt to use these stars as a labeled training set in order to classify stars which lack log g estimates into dwarfs and giants with a SVM algorithm. We assess the performance of this classification against different choices of the input feature vector. In particular, we use different combinations of reduced proper motions, 2MASS JHK, DENIS IJK and USNO-B B2R2 apparent magnitudes. Our study shows that -- for our color ranges -- the infrared bands alone provide no relevant information to separate dwarfs and giants. Even when optical bands and reduced proper motions are added, the fraction of true giants classified as dwarfs (the contamination) remains above 20%. Using only the dwarfs with available spectroscopic log g and distance estimates (the latter from Breddels et al. 2010), we then repeat the stream search by Klement, Fuchs & Rix (2008, KFR08), which assumed all stars were dwarfs and claimed the discovery of a new stellar stream at V = -160 km/s in a sample of 7015 stars from RAVE DR1. Our re-analysis of the pure DR2 dwarf sample exhibits an overdensity of 5 stars at the phase-space position of the KFR08 stream, with a metallicity distribution that appears inconsistent with that of stars at comparably low rotational velocities. Compared to several smooth Milky Way models, the mean standardized deviation of the KFR08 stream is only marginal at 1.6$\pm$0.4... (abbreviated)

We test competing models that aim at explaining the nature of stars in the Milky Way that are well away (|z|$\gtrsim$ 1kpc) from the midplane, the so-called thick disk: the stars may have gotten there through orbital migration, through satellite mergers and accretion, or through heating of pre-existing thin disk stars. Sales et al. (2009) proposed the eccentricity distribution of thick disk stars as a diagnostic to differentiate between these mechanisms. Drawing on SDSS DR7, we have assembled a sample of 34,223 G-dwarfs with 6-D phase-space information and metallicities, and have derived orbital eccentricities for them. Comparing the resulting eccentricity distributions, p(e|z), with the models, we find that: a) the observed p(e|z) is inconsistent with that predicted by orbital migration only, as there are more observed stars of high and of very low eccentricity; b) scenarios where the thick disk is made predominantly through abrupt heating of a pre-existing thin disk are also inconsistent, as they predict more high-eccentricity stars than observed; c) the observed p(e|z) fits well with a "gas-rich merger" scenario, where most thick disk stars were born from unsettled gas in situ.

We investigate the performance of some common machine learning techniques in identifying BHB stars from photometric data. To train the machine learning algorithms, we use previously published spectroscopic identifications of BHB stars from SDSS data. We investigate the performance of three different techniques, namely k nearest neighbour classification, kernel density estimation and a support vector machine (SVM). We discuss the performance of the methods in terms of both completeness and contamination. We discuss the prospect of trading off these values, achieving lower contamination at the expense of lower completeness, by adjusting probability thresholds for the classification. We also discuss the role of prior probabilities in the classification performance, and we assess via simulations the reliability of the dataset used for training. Overall it seems that no-prior gives the best completeness, but adopting a prior lowers the contamination. We find that the SVM generally delivers the lowest contamination for a given level of completeness, and so is our method of choice. Finally, we classify a large sample of SDSS DR7 photometry using the SVM trained on the spectroscopic sample. We identify 27,074 probable BHB stars out of a sample of 294,652 stars. We derive photometric parallaxes and demonstrate that our results are reasonable by comparing to known distances for a selection of globular clusters. We attach our classifications, including probabilities, as an electronic table, so that they can be used either directly as a BHB star catalogue, or as priors to a spectroscopic or other classification method. We also provide our final models so that they can be directly applied to new data.

The phase-space structure of our Galaxy holds the key to understand and reconstruct its formation. The Lambda-CDM model predicts a richly structured phase-space distribution of dark matter and (halo) stars, consisting of streams of particles torn from their progenitors during the process of hierarchical merging. While such streams quickly loose their spatial coherence in the process of phase mixing, the individual stars keep their common origin imprinted into their kinematic and chemical properties, allowing the recovery of the Galaxy's individual "building blocks". The field of Galactic Archeology has witnessed a dramatic boost over the last decade, thanks to the increasing quality and size of available data sets. This is especially true for the solar neighborhood, a volume of 1-2 kpc around the sun, where large scale surveys like SDSS/SEGUE continue to reveal the full 6D phase-space information of thousands of halo stars. In this review, I summarize the discoveries of stellar halo streams made so far and give a theoretical overview over the search strategies imployed. This paper is intended as an introduction to researchers new to field, but also as a reference illustrating the achievements made so far. I conclude that disentangling the individual fragments from which the Milky Way was built requires more precise data that will ultimately be delivered by the Gaia mission.

We have detected stellar halo streams in the solar neighborhood using data from the 7th public data release of the Sloan Digital Sky Survey (SDSS), which includes the directed stellar program SEGUE: Sloan Extension For Galactic Understanding and Exploration. In order to derive distances to each star, we used the metallicity-dependent photometric parallax relation from Ivezic et al. (2008) for which we examine and quantify the accuracy. Our final sample consists of 22,321 nearby (d < 2 kpc), metal-poor ([Fe/H] < -0.5) main-sequence stars with 6D estimates of position and space velocity. We characterize the orbits of these stars through suitable kinematic proxies for their "effective" integrals of motion, angular momentum, eccentricity, and orbital polar angle and compare the observed distribution to expectations from a smooth distribution in four [Fe/H] bins. On this basis we identify at least five significant "phase-space overdensities" of stars on very similar orbits in the solar neighborhood to which we can assign unambiguously peaked [Fe/H] distributions. Three of them have been identified previously, including the halo stream discovered by Helmi et al. (1999) at a significance level of 12.0. In addition, we find at least two new genuine halo streams, judged by their kinematics and [Fe/H], at significance levels of 2.9 and 4.8, respectively. For one stream the stars even show coherence in configuration space, matching a spatial overdensity of stars found by Juric et al. (2008) at (R,z) \approx (9.5,0.8) kpc. Our results demonstrate the practical power of our search method to detect substructure in the phase-space distribution of nearby stars without making a-priori assumptions about the detailed form of the gravitational potential.

We study the velocity distribution of Milky Way disk stars in a kiloparsec-sized region around the Sun, based on ~ 2 million M-type stars from DR7 of SDSS, which have newly re-calibrated absolute proper motions from combining SDSS positions with the USNO-B catalogue. We estimate photometric distances to all stars, accurate to ~ 20 %, and combine them with the proper motions to derive tangential velocities for this kinematically unbiased sample of stars. Based on a statistical de-projection method we then derive the vertical profiles (to heights of Z = 800 pc above the disk plane) for the first and second moments of the three dimensional stellar velocity distribution. We find that <W> = -7 +/- 1 km/s and <U> = -9 +/- 1 km/s, independent of height above the mid-plane, reflecting the Sun's motion with respect to the local standard of rest. In contrast, <V> changes distinctly from -20 +/- 2 km/s in the mid-plane to <V> = -32 km/s at Z = 800 pc, reflecting an asymmetric drift of the stellar mean velocity that increases with height. All three components of the M-star velocity dispersion show a strong linear rise away from the mid-plane, most notably \sigma_{ZZ}, which grows from 18 km/s (Z = 0) to 40 km/s (at Z = 800 pc). We determine the orientation of the velocity ellipsoid, and find a significant vertex deviation of 20 to 25 degrees, which decreases only slightly to heights of Z = 800 pc. Away from the mid-plane, our sample exhibits a remarkably large tilt of the velocity ellipsoid towards the Galactic plane, which reaches 20 deg. at Z = 800 pc and which is not easily explained. Finally, we determine the ratio \sigma^2_{\phi\phi}/\sigma^2_{RR} near the mid-plane, which in the epicyclic approximation implies an almost perfectly flat rotation curve at the Solar radius.

Context: Deep photometric surveys for substellar objects in young clusters and for high-redshift quasars are affected by contaminant sources at different heliocentric distances.If not correctly taken into account, the contamination may have a strong effect on the Initial Mass Function determination and on the identification of quasars. Aims: We calculate in detail the back- and foreground contamination by field dwarfs of very late spectral types (intermediate and late M, L and T) in deep surveys and provide the data and tools for the computation. Methods: Up-to-date models and data from the literature have been used: (i) a model of the Galactic thin disc by an exponential law; (ii) the length and height scales for late-type dwarfs; (iii) the local spatial densities, absolute magnitudes and colours of dwarfs for each spectral type. Results: We derive a simplified expression for the spatial density in the thin disc that depends on the heliocentric distance and the galactic coordinates (l, b) and integrate it in the truncated cone screened in the survey. As a practical application, we compute the numbers of L- and T-type field dwarfs in very deep (I = 21-29 mag) surveys in the direction of the young sigma Orionis cluster. The increasing number of contaminants at the faintest magnitudes could inhibit the study of the opacity mass limit at M <~ 0.003 M_sol in the cluster.

We searched for and detected stellar streams or moving groups in the solar neighbourhood, using the data provided by the 1st RAVE public data release. This analysis is based on distances to RAVE stars estimated from a color-magnitude relation that was calibrated on Hipparcos stars. Our final sample consists of 7015 stars selected to be within 500 pc of the Sun and to have distance errors better than 25%. Together with radial velocities from RAVE and proper motions from various data bases, there are estimates for all 6 phase-space coordinates of the stars in the sample. We characterize the orbits of these stars through suitable proxies for their angular momentum and eccentricity, and compare the observed distribution to the expectations from a smooth distribution. On this basis we identify at least four "phase space overdensities" of stars on very similar orbits in the Solar neighbourhood. We estimate the statistical significance of these overdensities by Monte Carlo simulations. Three of them have been identified previously: the Sirius and Hercules moving group and a stream found independently in 2006 by Arifyanto and Fuchs and Helmi et al. In addition, we have found a new stream candidate on a quite radial orbit, suggesting an origin external to the Milky Way's disk. Also, there is evidence for the Arcturus stream and the Hyades-Pleiades moving group in the sample. This analysis, using only a minute fraction of the final RAVE data set, shows the power of this experiment to probe the phase-space substructure of stars around the Sun.