The Cold Classical Kuiper Belt, a class of small bodies in undisturbed orbits beyond Neptune, are primitive objects preserving information about Solar System formation. The New Horizons spacecraft flew past one of these objects, the 36 km long contact binary (486958) Arrokoth (2014 MU69), in January 2019. Images from the flyby show that Arrokoth has no detectable rings, and no satellites (larger than 180 meters diameter) within a radius of 8000 km, and has a lightly-cratered smooth surface with complex geological features, unlike those on previously visited Solar System bodies. The density of impact craters indicates the surface dates from the formation of the Solar System. The two lobes of the contact binary have closely aligned poles and equators, constraining their accretion mechanism.

The problem of denoising a one-dimensional signal possessing varying degrees of smoothness is ubiquitous in time-domain astronomy and astronomical spectroscopy. For example, in the time domain, an astronomical object may exhibit a smoothly varying intensity that is occasionally interrupted by abrupt dips or spikes. Likewise, in the spectroscopic setting, a noiseless spectrum typically contains intervals of relative smoothness mixed with localized higher frequency components such as emission peaks and absorption lines. In this work, we present trend filtering, a modern nonparametric statistical tool that yields significant improvements in this broad problem space of denoising $spatially$ $heterogeneous$ signals. When the underlying signal is spatially heterogeneous, trend filtering is superior to any statistical estimator that is a linear combination of the observed data---including kernel smoothers, LOESS, smoothing splines, Gaussian process regression, and many other popular methods. Furthermore, the trend filtering estimate can be computed with practical and scalable efficiency via a specialized convex optimization algorithm, e.g. handling sample sizes of $n\gtrsim10^7$ within a few minutes. In a companion paper, we explicitly demonstrate the broad utility of trend filtering to observational astronomy by carrying out a diverse set of spectroscopic and time-domain analyses.

Trend filtering---first introduced into the astronomical literature in Paper I of this series---is a state-of-the-art statistical tool for denoising one-dimensional signals that possess varying degrees of smoothness. In this work, we demonstrate the broad utility of trend filtering to observational astronomy by discussing how it can contribute to a variety of spectroscopic and time-domain studies. The observations we discuss are (1) the Lyman-$\alpha$ forest of quasar spectra; (2) more general spectroscopy of quasars, galaxies, and stars; (3) stellar light curves with planetary transits; (4) eclipsing binary light curves; and (5) supernova light curves. We study the Lyman-$\alpha$ forest in the greatest detail---using trend filtering to map the large-scale structure of the intergalactic medium along quasar-observer lines of sight. The remaining studies share broad themes of: (1) estimating observable parameters of light curves and spectra; and (2) constructing observational spectral/light-curve templates. We also briefly discuss the utility of trend filtering as a tool for one-dimensional data reduction and compression.

We present the results from four stellar occultations by (486958) Arrokoth, the flyby target of the New Horizons extended mission. Three of the four efforts led to positive detections of the body, and all constrained the presence of rings and other debris, finding none. Twenty-five mobile stations were deployed for 2017 June 3 and augmented by fixed telescopes. There were no positive detections from this effort. The event on 2017 July 10 was observed by SOFIA with one very short chord. Twenty-four deployed stations on 2017 July 17 resulted in five chords that clearly showed a complicated shape consistent with a contact binary with rough dimensions of 20 by 30 km for the overall outline. A visible albedo of 10% was derived from these data. Twenty-two systems were deployed for the fourth event on 2018 Aug 4 and resulted in two chords. The combination of the occultation data and the flyby results provides a significant refinement of the rotation period, now estimated to be 15.9380 $\pm$ 0.0005 hours. The occultation data also provided high-precision astrometric constraints on the position of the object that were crucial for supporting the navigation for the New Horizons flyby. This work demonstrates an effective method for obtaining detailed size and shape information and probing for rings and dust on distant Kuiper Belt objects as well as being an important source of positional data that can aid in spacecraft navigation that is particularly useful for small and distant bodies.

Though pairs of dynamically associated asteroids in the Main Belt have been identified and studied for over a decade, very few pair systems have been identified in the near-Earth asteroid population. We present data and analysis that supports the existence of two genetically related pairs in near-Earth space. The members of the individual systems, 2015 EE7 -- 2015 FP124 and 2017 SN16 -- 2018 RY7, are found to be of the same spectral taxonomic class, and both pairs are interpreted to have volatile-poor compositions. In conjunction with dynamical arguments, this suggests that these two systems formed via YORP spin-up and/or dissociation of a binary precursor. Backwards orbital integrations suggest a separation age of <10 kyr for the pair 2017 SN16 -- 2018 RY7, making these objects amongst the youngest multiple asteroid systems known to date. A unique separation age was not realized for 2015 EE7 -- 2015 FP124 due to large uncertainties associated with these objects' orbits. Determining the ages of such young pairs is of great value for testing models of space weathering and asteroid spin-state evolution. As the NEO catalog continues to grow with current and future discovery surveys, it is expected that more NEO pairs will be found, thus providing an ideal laboratory for studying time dependent evolutionary processes that are relevant to asteroids throughout the Solar System.

Context. Pluto's tenuous nitrogen (N2) atmosphere undergoes strong seasonal effects due to high obliquity and orbital eccentricity, and has been recently (July 2015) observed by the New Horizons spacecraft. Goals are (i) construct a well calibrated record of the seasonal evolution of surface pressure on Pluto and (ii) constrain the structure of the lower atmosphere using a central flash observed in 2015. Method: eleven stellar occultations by Pluto observed between 2002 and 2016 are used to retrieve atmospheric profiles (density, pressure, temperature) between $\sim$5 km and $\sim$380 km altitude levels (i.e. pressures from about 10 microbar to 10 nanobar). Results: (i) Pressure has suffered a monotonic increase from 1988 to 2016, that is compared to a seasonal volatile transport model, from which tight constraints on a combination of albedo and emissivity of N2 ice are derived; (ii) A central flash observed on 2015 June 29 is consistent with New Horizons REX profiles, provided that (a) large diurnal temperature variations (not expected by current models) occur over Sputnik Planitia and/or (b) hazes with tangential optical depth of about 0.3 are present at 4-7 km altitude levels and/or (c) the nominal REX density values are overestimated by an implausibly large factor of about 20% and/or (d) higher terrains block part of the flash in the Charon facing hemisphere.

We observed the July 17, 2017 stellar occultation of HD 168233 by the Kuiper Belt Object (486958) 2014 MU$_{69}$, the close flyby target of the extended New Horizons mission. Rather than capture a solid body occultation by the KBO itself, our program aimed to constrain the opacity of rings, moons, or other debris in the nearby environment. We used the Hubble Space Telescope Fine Guidance Sensors (HST FGS) instrument in TRANS F583W mode to collect 40 Hz time resolution photometry of the stellar occultation star for two HST orbits during this observation. We present the results of reduction and calibration of the HST FGS photometry, and set upper limits on rings or other dust opacity within the Hill sphere of (486958) 2014 MU$_{69}$ at distances ranging from $\sim$20,000 km to $\sim$75,000 km from the main body.

We present the first double-lined orbital solution for the close binary in the GW Ori triple system. Using 12 epochs of infrared spectroscopy, we detected the lines of both stars in the inner pair, previously known as single-lined only. Our preliminary infrared orbital solution has an eccentricity of e = 0.21 +/- 0.10, a period of P = 241.15 +/- 0.72 days, and a mass ratio of q = 0.66 +/- 0.13. We find a larger semi-amplitude for the primary star, K1 = 6.57 +/- 1.00 km s^-1, with an infrared-only solution compared to K1 = 4.41 +/- 0.33 km s^-1 with optical data from the literature, likely the result of line blending and veiling in the optical. The component spectral types correspond to G3 and K0 stars, with vsini values of 43 km s^-1 and 50 km s^-1, respectively. We obtained a flux ratio of alpha = 0.58 +/- 0.14 in the H-band, allowing us to estimate individual masses of 3.2 and 2.7 M_sun for the primary and secondary, respectively, using evolutionary tracks. The tracks also yield a coeval age of 1 Myr for both components to within 1 sigma. GW Ori is surrounded by a circumbinary/circumtriple disk. A tertiary component has been detected in previous studies; however, we did not detect this component in our near-infrared spectra, probably the result of its relative faintness and blending in the absorption lines of these rapidly rotating stars. With these results, GW Ori joins the small number of classical T Tauri, double-lined spectroscopic binaries.

We have characterized the spectroscopic orbit of the TWA 3A binary and provide preliminary families of probable solutions for the TWA 3A visual orbit as well as for the wide TWA 3A--B orbit. TWA 3 is a hierarchical triple located at 34 pc in the $\sim$10 Myr old TW Hya association. The wide component separation is 1."55; the close pair was first identified as a possible binary almost 20 years ago. We initially identified the 35-day period orbital solution using high-resolution infrared spectroscopy which angularly resolved the A and B components. We then refined the preliminary orbit by combining the infrared data with a re-analysis of our high-resolution optical spectroscopy. The orbital period from the combined spectroscopic solution is $\sim$35 days, the eccentricity is $\sim$0.63, and the mass ratio is $\sim$0.84; although this high mass ratio would suggest that optical spectroscopy alone should be sufficient to identify the orbital solution, the presence of the tertiary B component likely introduced confusion in the blended optical spectra. Using millimeter imaging from the literature, we also estimate the inclinations of the stellar orbital planes with respect to the TWA 3A circumbinary disk inclination and find that all three planes are likely misaligned by at least $\sim$30 degrees. The TWA 3A spectroscopic binary components have spectral types of M4.0 and M4.5; TWA 3B is an M3. We speculate that the system formed as a triple, is bound, and that its properties were shaped by dynamical interactions between the inclined orbits and disk.

We study the effects of filaments on galaxy properties in the Sloan Digital Sky Survey (SDSS) Data Release 12 using filaments from the `Cosmic Web Reconstruction' catalogue (Chen et al. 2016), a publicly available filament catalogue for SDSS. Since filaments are tracers of medium-to-high density regions, we expect that galaxy properties associated with the environment are dependent on the distance to the nearest filament. Our analysis demonstrates that a red galaxy or a high-mass galaxy tend to reside closer to filaments than a blue or low-mass galaxy. After adjusting the effect from stellar mass, on average, early-forming galaxies or large galaxies have a shorter distance to filaments than late-forming galaxies or small galaxies. For the Main galaxy sample (MGS), all signals are very significant ($>6\sigma$). For the LOWZ and CMASS sample, the stellar mass and size are significant ($>2 \sigma$). The filament effects we observe persist until $z = 0.7$ (the edge of the CMASS sample). Comparing our results to those using the galaxy distances from redMaPPer galaxy clusters as a reference, we find a similar result between filaments and clusters. Moreover, we find that the effect of clusters on the stellar mass of nearby galaxies depends on the galaxy's filamentary environment. Our findings illustrate the strong correlation of galaxy properties with proximity to density ridges, strongly supporting the claim that density ridges are good tracers of filaments.

We present a Mass-Luminosity Relation (MLR) for red dwarfs spanning a range of masses from 0.62 Msun to the end of the stellar main sequence at 0.08 Msun. The relation is based on 47 stars for which dynamical masses have been determined, primarily using astrometric data from Fine Guidance Sensors (FGS) 3 and 1r, white-light interferometers on the Hubble Space Telescope (HST), and radial velocity data from McDonald Observatory. For our HST/FGS sample of 15 binaries component mass errors range from 0.4% to 4.0% with a median error of 1.8%. With these and masses from other sources, we construct a V-band MLR for the lower main sequence with 47 stars, and a K-band MLR with 45 stars with fit residuals half of those of the V-band. We use GJ 831 AB as an analysis example, obtaining an absolute trigonometric parallax, pi_abs = 125.3 +/- 0.3 milliseconds of arc, with orbital elements yielding MA = 0.270 +/- 0.004 Msun and MB = 0.145 +/- 0.002 Msun. The mass precision rivals that derived for eclipsing binaries. A remaining major task is the interpretation of the intrinsic cosmic scatter in the observed MLR for low mass stars in terms of physical effects. In the meantime, useful mass values can be estimated from the MLR for the ubiquitous red dwarfs that account for 75% of all stars, with applications ranging from the characterization of exoplanet host stars to the contribution of red dwarfs to the mass of the Universe.

We present results derived from the first multi-chord stellar occultation by the trans-Neptunian object (229762) 2007 UK$_{126}$, observed on 2014 November 15. The event was observed by the Research and Education Collaborative Occultation Network (RECON) project and International Occultation Timing Association (IOTA) collaborators throughout the United States. Use of two different data analysis methods obtain a satisfactory fit to seven chords, yelding an elliptical fit to the chords with an equatorial radius of $R=338_{-10} ^{+15}$ km and equivalent radius of $R_{eq}=319_{-7} ^{+14}$ km. A circular fit also gives a radius of $R=324_{-23} ^{+30}$ km. Assuming that the object is a Maclaurin spheroid with indeterminate aspect angle, and using two published absolute magnitudes for the body, we derive possible ranges for geometric albedo between $p_{V}=0.159_{-0.013} ^{+0.007}$ and $p_{R}=0.189_{-0.015}^{+0.009}$, and for the body oblateness between $\epsilon=0.105_{-0.040} ^{+0.050}$ and $\epsilon=0.118_{-0.048} ^{+0.055}$. For a nominal rotational period of 11.05 h, an upper limit for density of $\rho=1740$ kg~m$^{-3}$ is estimated for the body.

We report dynamical mass measurements for the components of the previously known double-lined spectroscopic subdwarfs G 006-026 B and C using the Fine Guidance Sensors (FGS) on the Hubble Space Telescope. To build the empirical mass-luminosity relation for low metallicity subdwarfs, we collect four other subdwarf systems with dynamical masses that we compare to theoretical models for various metallicities on the mass-luminosity relation. For most stars, they fall in the regions where the models predict to be low metallicity. This effort highlights the scarcity of dynamical masses for subdwarfs and that much work remains to be done to improve the mass errors and metallicity measurements of low mass subdwarfs in our Galaxy.

We construct a catalogue for filaments using a novel approach called SCMS (subspace constrained mean shift; Ozertem & Erdogmus 2011; Chen et al. 2015). SCMS is a gradient-based method that detects filaments through density ridges (smooth curves tracing high-density regions). A great advantage of SCMS is its uncertainty measure, which allows an evaluation of the errors for the detected filaments. To detect filaments, we use data from the Sloan Digital Sky Survey, which consist of three galaxy samples: the NYU main galaxy sample (MGS), the LOWZ sample and the CMASS sample. Each of the three dataset covers different redshift regions so that the combined sample allows detection of filaments up to z = 0.7. Our filament catalogue consists of a sequence of two-dimensional filament maps at different redshifts that provide several useful statistics on the evolution cosmic web. To construct the maps, we select spectroscopically confirmed galaxies within 0.050 < z < 0.700 and partition them into 130 bins. For each bin, we ignore the redshift, treating the galaxy observations as a 2-D data and detect filaments using SCMS. The filament catalogue consists of 130 individual 2-D filament maps, and each map comprises points on the detected filaments that describe the filamentary structures at a particular redshift. We also apply our filament catalogue to investigate galaxy luminosity and its relation with distance to filament. Using a volume-limited sample, we find strong evidence (6.1$\sigma$ - 12.3$\sigma$) that galaxies close to filaments are generally brighter than those at significant distance from filaments.

The detection and characterization of filamentary structures in the cosmic web allows cosmologists to constrain parameters that dictates the evolution of the Universe. While many filament estimators have been proposed, they generally lack estimates of uncertainty, reducing their inferential power. In this paper, we demonstrate how one may apply the Subspace Constrained Mean Shift (SCMS) algorithm (Ozertem and Erdogmus (2011); Genovese et al. (2012)) to uncover filamentary structure in galaxy data. The SCMS algorithm is a gradient ascent method that models filaments as density ridges, one-dimensional smooth curves that trace high-density regions within the point cloud. We also demonstrate how augmenting the SCMS algorithm with bootstrap-based methods of uncertainty estimation allows one to place uncertainty bands around putative filaments. We apply the SCMS method to datasets sampled from the P3M N-body simulation, with galaxy number densities consistent with SDSS and WFIRST-AFTA and to LOWZ and CMASS data from the Baryon Oscillation Spectroscopic Survey (BOSS). To further assess the efficacy of SCMS, we compare the relative locations of BOSS filaments with galaxy clusters in the redMaPPer catalog, and find that redMaPPer clusters are significantly closer (with p-values $< 10^{-9}$) to SCMS-detected filaments than to randomly selected galaxies.

In this paper, we study the filamentary structures and the galaxy alignment along filaments at redshift $z=0.06$ in the MassiveBlack-II simulation, a state-of-the-art, high-resolution hydrodynamical cosmological simulation which includes stellar and AGN feedback in a volume of (100 Mpc$/h$)$^3$. The filaments are constructed using the subspace constrained mean shift (SCMS; Ozertem & Erdogmus (2011) and Chen et al. (2015a)). First, we show that reconstructed filaments using galaxies and reconstructed filaments using dark matter particles are similar to each other; over $50\%$ of the points on the galaxy filaments have a corresponding point on the dark matter filaments within distance $0.13$ Mpc$/h$ (and vice versa) and this distance is even smaller at high-density regions. Second, we observe the alignment of the major principal axis of a galaxy with respect to the orientation of its nearest filament and detect a $2.5$ Mpc$/h$ critical radius for filament's influence on the alignment when the subhalo mass of this galaxy is between $10^9M_\odot/h$ and $10^{12}M_\odot/h$. Moreover, we find the alignment signal to increase significantly with the subhalo mass. Third, when a galaxy is close to filaments (less than $0.25$ Mpc$/h$), the galaxy alignment toward the nearest galaxy group depends on the galaxy subhalo mass. Finally, we find that galaxies close to filaments or groups tend to be rounder than those away from filaments or groups.

We present an observational and dynamical study of newly discovered main-belt comet 313P/Gibbs. We find that the object is clearly active both in observations obtained in 2014 and in precovery observations obtained in 2003 by the Sloan Digital Sky Survey, strongly suggesting that its activity is sublimation-driven. This conclusion is supported by a photometric analysis showing an increase in the total brightness of the comet over the 2014 observing period, and dust modeling results showing that the dust emission persists over at least three months during both active periods, where we find start dates for emission no later than 2003 July 24+/-10 for the 2003 active period and 2014 July 28+/-10 for the 2014 active period. From serendipitous observations by the Subaru Telescope in 2004 when the object was apparently inactive, we estimate that the nucleus has an absolute R-band magnitude of H_R=17.1+/-0.3, corresponding to an effective nucleus radius of r_e~1.00+/-0.15 km. The object's faintness at that time means we cannot rule out the presence of activity, and so this computed radius should be considered an upper limit. We find that 313P's orbit is intrinsically chaotic, having a Lyapunov time of T_l=12000 yr and being located near two 3-body mean-motion resonances with Jupiter and Saturn, 11J-1S-5A and 10J+12S-7A, yet appears stable over >50 Myr in an apparent example of stable chaos. We furthermore find that 313P is the second main-belt comet, after P/2012 T1 (PANSTARRS), to belong to the ~155 Myr old Lixiaohua asteroid family.

The incredible variety of galaxy shapes cannot be summarized by human defined discrete classes of shapes without causing a possibly large loss of information. Dictionary learning and sparse coding allow us to reduce the high dimensional space of shapes into a manageable low dimensional continuous vector space. Statistical inference can be done in the reduced space via probability distribution estimation and manifold estimation.

The Deep Ecliptic Survey (DES) discovered hundreds of Kuiper Belt objects from 1998-2005. Follow-up observations yielded 304 objects with good dynamical classifications (Classical, Scattered, Centaur, or 16 mean-motion resonances with Neptune). The DES search fields are well documented, enabling us to calculate the probability of detecting objects with particular orbital parameters and absolute magnitudes at a randomized point in each orbit. Grouping objects together by dynamical class leads, we estimate the orbital element distributions (a, e, i) for the largest three classes (Classical, 3:2, and Scattered) using maximum likelihood. Using H-magnitude as a proxy for the object size, we fit a power law to the number of objects for 8 classes with at least 5 detected members (246 objects). The best Classical slope is alpha=1.02+/-0.01 (observed from 5<=H<=7.2). Six dynamical classes (Scattered plus 5 resonances) are consistent in slope with the Classicals, though the absolute number of objects is scaled. The exception to the power law relation are the Centaurs (non-resonant with perihelia closer than Neptune, and thus detectable at smaller sizes), with alpha=0.42+/-0.02 (7.5<H<11). This is consistent with a knee in the H-distribution around H=7.2 as reported elsewhere (Bernstein et al. 2004, Fraser et al. 2014). Based on the Classical-derived magnitude distribution, the total number of objects (H<=7) in each class are: Classical (2100+/-300 objects), Scattered (2800+/-400), 3:2 (570+/-80), 2:1 (400+/-50), 5:2 (270+/-40), 7:4 (69+/-9), 5:3 (60+/-8). The independent estimate for the number of Centaurs in the same H range is 13+/-5. If instead all objects are divided by inclination into "Hot" and "Cold" populations, following Fraser et al. (2014), we find that alphaHot=0.90+/-0.02, while alphaCold=1.32+/-0.02, in good agreement with that work.

The Lyman-alpha forest is a portion of the observed light spectrum of distant galactic nuclei which allows us to probe remote regions of the Universe that are otherwise inaccessible. The observed Lyman-alpha forest of a quasar light spectrum can be modeled as a noisy realization of a smooth curve that is affected by a `damping effect' which occurs whenever the light emitted by the quasar travels through regions of the Universe with higher matter concentration. To decode the information conveyed by the Lyman-alpha forest about the matter distribution, we must be able to separate the smooth `continuum' from the noise and the contribution of the damping effect in the quasar light spectra. To predict the continuum in the Lyman-alpha forest, we use a nonparametric functional regression model in which both the response and the predictor variable (the smooth part of the damping-free portion of the spectrum) are function-valued random variables. We demonstrate that the proposed method accurately predicts the unobservable continuum in the Lyman-alpha forest both on simulated spectra and real spectra. Also, we introduce distribution-free prediction bands for the nonparametric functional regression model that have finite sample guarantees. These prediction bands, together with bootstrap-based confidence bands for the projection of the mean continuum on a fixed number of principal components, allow us to assess the degree of uncertainty in the model predictions.

Visualizing the high-redshift Universe is difficult due to the dearth of available data; however, the Lyman-alpha forest provides a means to map the intergalactic medium at redshifts not accessible to large galaxy surveys. Large-scale structure surveys, such as the Baryon Oscillation Spectroscopic Survey (BOSS), have collected quasar (QSO) spectra that enable the reconstruction of HI density fluctuations. The data fall on a collection of lines defined by the lines-of-sight (LOS) of the QSO, and a major issue with producing a 3D reconstruction is determining how to model the regions between the LOS. We present a method that produces a 3D map of this relatively uncharted portion of the Universe by employing local polynomial smoothing, a nonparametric methodology. The performance of the method is analyzed on simulated data that mimics the varying number of LOS expected in real data, and then is applied to a sample region selected from BOSS. Evaluation of the reconstruction is assessed by considering various features of the predicted 3D maps including visual comparison of slices, PDFs, counts of local minima and maxima, and standardized correlation functions. This 3D reconstruction allows for an initial investigation of the topology of this portion of the Universe using persistent homology.

By analyzing brightness variation with ecliptic longitude and using the Lowell Observatory photometric database, we estimate spin-axis longitudes for more than 350 000 asteroids. Hitherto, spin-axis longitude estimates have been made for fewer than 200 asteroids. We investigate longitude distributions in different dynamical groups and asteroid families. We show that asteroid spin-axis longitudes are not isotropically distributed as previously considered. We find that the spin-axis longitude distribution for main-belt asteroids is clearly non-random, with an excess of longitudes from the interval 30{\deg}-110{\deg} and a paucity between 120{\deg}-180{\deg}. The explanation of the non-isotropic distribution is unknown at this point. Further studies have to be conducted to determine if the shape of the distribution can be explained by observational bias, selection effects, a real physical process or other mechanism.

We report orbital parameters for two low-mass, pre-main sequence, double-lined spectroscopic binaries VSB 111 and VSB 126. These systems were originally identified as single-lined on the basis of visible-light observations. We obtained high-resolution, infrared spectra with the 10-m Keck II telescope, detected absorption lines of the secondary stars, and measured radial velocities of both components in the systems. The visible light spectra were obtained on the 1.5-m Wyeth reflector at the Oak Ridge Observatory, the 1.5-m Tillinghast reflector at the F. L. Whipple Observatory, and the 4.5-m equivalent Multiple Mirror Telescope. The combination of our visible and infrared observations of VSB 111 leads to a period of 902.1+/-0.9 days, an eccentricity of 0.788+/-0.008, and a mass ratio of 0.52+/-0.05. VSB 126 has a period of 12.9244+/-0.0002 days, an eccentricity of 0.18+/-0.02, and a mass ratio of 0.29+/-0.02. Visible-light photometry, using the 0.8-m telescope at Lowell Observatory, provided rotation periods for the primary stars in both systems, 3.74+/-0.02 days for VSB 111 and 5.71+/-0.07 days for VSB 126. Both binaries are located in the young, active star-forming cluster NGC 2264 at a distance of ~800 pc. The difference in the center-of-mass velocities of the two systems is consistent with the radial velocity gradient seen across NGC 2264. To test the evolutionary models for accuracy and consistency, we compare the stellar properties derived from several sets of theoretical calculations for pre-main sequence evolution with our dynamical results.

Combining stellar occultation observations probing Pluto's atmosphere from 1988 to 2013 and models of energy balance between Pluto's surface and atmosphere, we conclude that Pluto's atmosphere does not collapse at any point in its 248-year orbit. The occultation results show an increasing atmospheric pressure with time in the current epoch, a trend present only in models with a high thermal inertia and a permanent N2 ice cap at Pluto's north rotational pole.

We report the discovery that the pre-main sequence object LkCa3 in the Taurus-Auriga star-forming region is a hierarchical quadruple system of M stars. It was previously known to be a close (~0.5 arc sec) visual pair, with one component being a moderately eccentric 12.94-day single-lined spectroscopic binary. A re-analysis of archival optical spectra complemented with new near-infrared spectroscopy shows both visual components to be double-lined, the second one having a period of 4.06 days and a circular orbit. In addition to the orbital elements, we determine optical and near-infrared flux ratios, effective temperatures, and projected rotational velocities for all four stars. Using existing photometric monitoring observations of the system that had previously revealed the rotational period of the primary in the longer-period binary, we detect also the rotational signal of the primary in the 4.06-day binary, which is synchronized with the orbital motion. With only the assumption of coevality, a comparison of all of these constraints with current stellar evolution models from the Dartmouth series points to an age of 1.4 Myr and a distance of 133 pc, consistent with previous estimates for the region and suggesting the system is on the near side of the Taurus complex. Similar comparisons of the properties of LkCa3 and of the well-known quadruple pre-main sequence system GG Tau with the widely used models from the Lyon series for a mixing length parameter of alpha_ML = 1.0 strongly favor the Dartmouth models.

As part of our continuing campaign to measure the masses of pre-main sequence stars dynamically and thus to assess the reliability of the discrepant theoretical calculations of contraction to the main sequence, we present new results for NTTS 045251+3016, a visual and double-lined spectroscopic binary in the Taurus Star Forming Region. We obtained new high angular resolution astrometry and high spectral resolution spectroscopy at Keck Observatory. The new data lead to a significant revision of previously published orbital parameters. In particular, we find that the masses of the primary and secondary are 0.86+/-0.11 and 0.55+/-0.05 M_sun, respectively, smaller than previously reported, and that the system lies 158.7+/-3.9 pc from the sun, further than previously reported.

We present the discovery of a long-term stable L5 (trailing) Neptune Trojan in data acquired to search for candidate Trans-Neptunian objects for the New Horizons spacecraft to fly by during an extended post-Pluto mission. This Neptune Trojan, 2011 HM102, has the highest inclination (29.4 degrees) of any known member of this population. It is intrinsically brighter than any single L5 Jupiter Trojan at H~8.18. We have determined its gri colors (a first for any L5 Neptune Trojan), which we find to be similar to the moderately red colors of the L4 Neptune Trojans, suggesting similar surface properties for members of both Trojan clouds. We also present colors derived from archival data for two L4 Neptune Trojans (2006 RJ103 and 2007 VL305), better refining the overall color distribution of the population. In this document we describe the discovery circumstances, our physical characterization of 2011 HM102, and this object's implications for the Neptune Trojan population overall. Finally, we discuss the prospects for detecting 2011 HM102 from the New Horizons spacecraft during their close approach in mid- to late-2013.

RX J0513.1+0851 and RX J0539.9+0956 were previously identified as young, low-mass, single-lined spectroscopic binary systems and classified as weak-lined T Tauri stars at visible wavelengths. Here we present radial velocities, spectral types, vsini values, and flux ratios for the components in these systems resulting from two-dimensional cross-correlation analysis. These results are based on high-resolution, near-infrared spectroscopy taken with the Keck II telescope to provide a first characterization of these systems as double-lined rather than single-lined. It applies the power of infrared spectroscopy to the detection of cool secondaries; the flux scales as a less steep function of mass in the infrared than in the visible, thus enabling an identification of low-mass secondaries. We found that the RX J0513.1+0851 and RX J0539.9+0956 primary stars are fast rotators, 60 km/s and 80 km/s respectively; this introduces extra difficulty in the detection of the secondary component as a result of the quite broad absorption lines. To date, these are the highest rotational velocities measured for a pre-main sequence spectroscopic binary. The orbital parameters and mass ratios were determined by combining new visible light spectroscopy with our infrared data for both systems. For RX J0513.1+0851, we derived a period of ~4 days and a mass ratio of q = 0.46 +/- 0.01 and for RX J0539.9+0956, a period of ~1117 days and a mass ratio of q = 0.66 +/- 0.01. Based on our derived properties for the stellar components, we estimate the luminosities and hence distances to these binaries at 220 pc and 90 pc. They appear to be significantly closer than previously estimated.

The young, low-mass, triple system NTTS 155808-2219 (ScoPMS 20) was previously identified as a ~17-day period single-lined spectroscopic binary with a tertiary component at 0.21 arcseconds. Using high-resolution infrared spectra, acquired with NIRSPEC on Keck II, both with and without adaptive optics, we measured radial velocities of all three components. Reanalysis of the single-lined visible light observations, made from 1987 to 1993, also yielded radial velocity detections of the three stars. Combining visible light and infrared data to compute the orbital solution produces orbital parameters consistent with the single-lined solution and a mass ratio of q = 0.78 +/- 0.01 for the SB. We discuss the consistency between our results and previously published data on this system, our radial-velocity analysis with both observed and synthetic templates, and the possibility that this system is eclipsing, providing a potential method for the determination of the stars' absolute masses. Over the ~20 year baseline of our observations, we have measured the acceleration of the SB's center-of-mass in its orbit with the tertiary. Long-term, adaptive optics imaging of the tertiary will eventually yield dynamical data useful for component mass estimates.

Hubble Space Telescope observations between 2001 and 2010 resolved the binary components of the Cold Classical transneptunian object (79360) Sila-Nunam (provisionally designated 1997 CS29). From these observations we have determined the circular, retrograde mutual orbit of Nunam relative to Sila with a period of 12.50995 \pm 0.00036 days and a semimajor axis of 2777 \pm 19 km. A multi-year season of mutual events, in which the two near-equal brightness bodies alternate in passing in front of one another as seen from Earth, is in progress right now, and on 2011 Feb. 1 UT, one such event was observed from two different telescopes. The mutual event season offers a rich opportunity to learn much more about this barely-resolvable binary system, potentially including component sizes, colors, shapes, and albedo patterns. The low eccentricity of the orbit and a photometric lightcurve that appears to coincide with the orbital period are consistent with a system that is tidally locked and synchronized, like the Pluto-Charon system. The orbital period and semimajor axis imply a system mass of (10.84 \pm 0.22) \times 10^18 kg, which can be combined with a size estimate based on Spitzer and Herschel thermal infrared observations to infer an average bulk density of 0.72 +0.37 -0.23 g cm^-3, comparable to the very low bulk densities estimated for small transneptunian binaries of other dynamical classes.

We present observations of comet-like main-belt object P/2010 R2 (La Sagra) obtained by Pan-STARRS 1 and the Faulkes Telescope-North on Haleakala in Hawaii, the University of Hawaii 2.2 m, Gemini-North, and Keck I telescopes on Mauna Kea, the Danish 1.54 m telescope at La Silla, and the Isaac Newton Telescope on La Palma. An antisolar dust tail is observed from August 2010 through February 2011, while a dust trail aligned with the object's orbit plane is also observed from December 2010 through August 2011. Assuming typical phase darkening behavior, P/La Sagra is seen to increase in brightness by >1 mag between August 2010 and December 2010, suggesting that dust production is ongoing over this period. These results strongly suggest that the observed activity is cometary in nature (i.e., driven by the sublimation of volatile material), and that P/La Sagra is therefore the most recent main-belt comet to be discovered. We find an approximate absolute magnitude for the nucleus of H_R=17.9+/-0.2 mag, corresponding to a nucleus radius of ~0.7 km, assuming an albedo of p=0.05. Using optical spectroscopy, we find no evidence of sublimation products (i.e., gas emission), finding an upper limit CN production rate of Q_CN<6x10^23 mol/s, from which we infer an H2O production rate of Q_H2O<10^26 mol/s. Numerical simulations indicate that P/La Sagra is dynamically stable for >100 Myr, suggesting that it is likely native to its current location and that its composition is likely representative of other objects in the same region of the main belt, though the relatively close proximity of the 13:6 mean-motion resonance with Jupiter and the (3,-2,-1) three-body mean-motion resonance with Jupiter and Saturn mean that dynamical instability on larger timescales cannot be ruled out.

We explore the correlation between an asteroid's taxonomy and photometric phase curve using the H, G12 photometric phase function, with the shape of the phase function described by the single parameter G12. We explore the usability of G12 in taxonomic classification for individual objects, asteroid families, and dynamical groups. We conclude that the mean values of G12 for the considered taxonomic complexes are statistically different, and also discuss the overall shape of the G12 distribution for each taxonomic complex. Based on the values of G12 for about half a million asteroids, we compute the probabilities of C, S, and X complex membership for each asteroid. For an individual asteroid, these probabilities are rather evenly distributed over all of the complexes, thus preventing meaningful classification. We then present and discuss the G12 distributions for asteroid families, and predict the taxonomic complex preponderance for asteroid families given the distribution of G12 in each family. For certain asteroid families, the probabilistic prediction of taxonomic complex preponderance can clearly be made. The Nysa-Polana family shows two distinct regions in the proper element space with different G12 values dominating in each region. We conclude that the G12-based probabilistic distribution of taxonomic complexes through the main belt agrees with the general view of C complex asteroid proportion increasing towards the outer belt. We conclude that the G12 photometric parameter cannot be used in determining taxonomic complex for individual asteroids, but it can be utilized in the statistical treatment of asteroid families and different regions of the main asteroid belt.

With the NEOWISE portion of the \emph{Wide-field Infrared Survey Explorer} (WISE) project, we have carried out a highly uniform survey of the near-Earth object (NEO) population at thermal infrared wavelengths ranging from 3 to 22 $\mu$m, allowing us to refine estimates of their numbers, sizes, and albedos. The NEOWISE survey detected NEOs the same way whether they were previously known or not, subject to the availability of ground-based follow-up observations, resulting in the discovery of more than 130 new NEOs. The survey's uniformity in sensitivity, observing cadence, and image quality have permitted extrapolation of the 428 near-Earth asteroids (NEAs) detected by NEOWISE during the fully cryogenic portion of the WISE mission to the larger population. We find that there are 981$\pm$19 NEAs larger than 1 km and 20,500$\pm$3000 NEAs larger than 100 m. We show that the Spaceguard goal of detecting 90% of all 1 km NEAs has been met, and that the cumulative size distribution is best represented by a broken power law with a slope of 1.32$\pm$0.14 below 1.5 km. This power law slope produces $\sim13,200\pm$1,900 NEAs with $D>$140 m. Although previous studies predict another break in the cumulative size distribution below $D\sim$50-100 m, resulting in an increase in the number of NEOs in this size range and smaller, we did not detect enough objects to comment on this increase. The overall number for the NEA population between 100-1000 m are lower than previous estimates. The numbers of near-Earth comets will be the subject of future work.

We present initial results from time series imaging at infrared wavelengths of 0.9 sq. degrees in the Orion Nebula Cluster (ONC). During Fall 2009 we obtained 81 epochs of Spitzer 3.6 and 4.5 micron data over 40 consecutive days. We extracted light curves with ~3% photometric accuracy for ~2000 ONC members ranging from several solar masses down to well below the hydrogen burning mass limit. For many of the stars, we also have time-series photometry obtained at optical (Ic) and/or near-infrared (JKs) wavelengths. Our data set can be mined to determine stellar rotation periods, identify new pre-main-sequence (PMS) eclipsing binaries, search for new substellar Orion members, and help better determine the frequency of circumstellar disks as a function of stellar mass in the ONC. Our primary focus is the unique ability of 3.6 & 4.5 micron variability information to improve our understanding of inner disk processes and structure in the Class I and II young stellar objects (YSOs). In this paper, we provide a brief overview of the YSOVAR Orion data obtained in Fall 2009, and we highlight our light curves for AA-Tau analogs - YSOs with narrow dips in flux, most probably due to disk density structures passing through our line of sight. Detailed follow-up observations are needed in order to better quantify the nature of the obscuring bodies and what this implies for the structure of the inner disks of YSOs.

We apply statistically rigorous methods of nonparametric risk estimation to the problem of inferring the local peculiar velocity field from nearby supernovae (SNIa). We use two nonparametric methods - Weighted Least Squares (WLS) and Coefficient Unbiased (CU) - both of which employ spherical harmonics to model the field and use the estimated risk to determine at which multipole to truncate the series. We show that if the data are not drawn from a uniform distribution or if there is power beyond the maximum multipole in the regression, a bias is introduced on the coefficients using WLS. CU estimates the coefficients without this bias by including the sampling density making the coefficients more accurate but not necessarily modeling the velocity field more accurately. After applying nonparametric risk estimation to SNIa data, we find that there are not enough data at this time to measure power beyond the dipole. The WLS Local Group bulk flow is moving at 538 +- 86 km/s towards (l,b) = (258 +- 10 deg, 36 +- 11 deg) and the CU bulk flow is moving at 446 +- 101 km/s towards (l,b) = (273 +- 11 deg, 46 +- 8 deg). We find that the magnitude and direction of these measurements are in agreement with each other and previous results in the literature.

We consider the problem of estimating filamentary structure from planar point process data. We make some connections with computational geometry and we develop nonparametric methods for estimating the filaments. We show that, under weak conditions, the filaments have a simple geometric representation as the medial axis of the data distribution's support. Our methods convert an estimator of the support's boundary into an estimator of the filaments. We also find the rates of convergence of our estimators.

We report the orbital parameters for ROXR1 14 and RX J1622.7-2325Nw, two young, low-mass, and double-lined spectroscopic binaries recently discovered in the Ophiuchus star forming region. Accurate orbital solutions were determined from over a dozen high-resolution spectra taken with the Keck II and Gemini South telescopes. These objects are T Tauri stars with mass ratios close to unity and periods of ~5 and ~3 days, respectively. In particular, RX J1622.7-2325Nw shows a non-circularized orbit with an eccentricity of 0.30, higher than any other short-period pre-main sequence spectroscopic binary known to date. We speculate that orbit of RX J1622.7-2325Nw has not yet circularized because of the perturbing action of a ~1" companion, itself a close visual pair. A comparison of known young spectroscopic binaries and main sequence (MS) spectroscopic binaries in the eccentricity-period plane shows an indistinguishable distribution of the two populations, implying that orbital circularization occurs in the first 1 Myr of a star's lifetime. With the results presented in this paper we are increasing by ~4% the small sample of PMS spectroscopic binary stars with known orbital elements.

Using data from the Deep Ecliptic Survey (DES), we investigate the inclination distributions of objects in the Kuiper Belt. We present a derivation for observational bias removal and use this procedure to generate unbiased inclination distributions for Kuiper Belt objects (KBOs) of different DES dynamical classes, with respect to the Kuiper Belt Plane. Consistent with previous results, we find that the inclination distribution for all DES KBOs is well fit by the sum of two Gaussians, or a Gaussian plus a generalized Lorentzian, multiplied by sin i. Approximately 80% of KBOs are in the high-inclination grouping. We find that Classical object inclinations are well fit by sin i multiplied by the sum of two Gaussians, with roughly even distribution between Gaussians of widths 2.0 -0.5/+0.6 degrees and 8.1 -2.1/+2.6 degrees. Objects in different resonances exhibit different inclination distributions. The inclinations of Scattered objects are best matched by sin i multiplied by a single Gaussian that is centered at 19.1 -3.6/+3.9 degrees with a width of 6.9 -2.7/+4.1 degrees. Centaur inclinations peak just below 20 degrees, with one exceptionally high-inclination object near 80 degrees. The currently observed inclination distribution of the Centaurs is not dissimilar to that of the Scattered Extended KBOs and Jupiter-family comets, but is significantly different from the Classical and Resonant KBOs. While the sample sizes of some dynamical classes are still small, these results should begin to serve as a critical diagnostic for models of Solar System evolution.

The surprising discovery of an accelerating universe led cosmologists to posit the existence of "dark energy"--a mysterious energy field that permeates the universe. Understanding dark energy has become the central problem of modern cosmology. After describing the scientific background in depth, we formulate the task as a nonlinear inverse problem that expresses the comoving distance function in terms of the dark-energy equation of state. We present two classes of methods for making sharp statistical inferences about the equation of state from observations of Type Ia Supernovae (SNe). First, we derive a technique for testing hypotheses about the equation of state that requires no assumptions about its form and can distinguish among competing theories. Second, we present a framework for computing parametric and nonparametric estimators of the equation of state, with an associated assessment of uncertainty. Using our approach, we evaluate the strength of statistical evidence for various competing models of dark energy. Consistent with current studies, we find that with the available Type Ia SNe data, it is not possible to distinguish statistically among popular dark-energy models, and that, in particular, there is no support in the data for rejecting a cosmological constant. With much more supernova data likely to be available in coming years (e.g., from the DOE/NASA Joint Dark Energy Mission), we address the more interesting question of whether future data sets will have sufficient resolution to distinguish among competing theories.

The young system RX J0529.3+1210 was initially identified as a single-lined spectroscopic binary. Using high-resolution infrared spectra, acquired with NIRSPEC on Keck II, we measured radial velocities for the secondary. The method of using the infrared regime to convert single-lined spectra into double-lined spectra, and derive the mass ratio for the binary system, has been successfully used for a number of young, low-mass binaries. For RX J0529.3+1210, a long- period(462 days) and highly eccentric(0.88) binary system, we determine the mass ratio to be 0.78+/-0.05 using the infrared double-lined velocity data alone, and 0.73+/-0.23 combining visible light and infrared data in a full orbital solution. The large uncertainty in the latter is the result of the sparse sampling in the infrared and the high eccentricity: the stars do not have a large velocity separation during most of their ~1.3 year orbit. A mass ratio close to unity, consistent with the high end of the one sigma uncertainty for this mass ratio value, is inconsistent with the lack of a visible light detection of the secondary component. We outline several scenarios for a color difference in the two stars, such as one heavily spotted component, higher order multiplicity, or a unique evolutionary stage, favoring detection of only the primary star in visible light, even in a mass ratio ~1 system. However, the evidence points to a lower ratio. Although RX J0529.3+1210 exhibits no excess at near-infrared wavelengths, a small 24 micron excess is detected, consistent with circumbinary dust. The properties of this binary and its membership in Lambda Ori versus a new nearby stellar moving group at ~90 pc are discussed. We speculate on the origin of this unusual system and on the impact of such high eccentricity on the potential for planet formation.

The distributions of galaxy properties vary with environment, and are often multimodal, suggesting that the galaxy population may be a combination of multiple components. The behaviour of these components versus environment holds details about the processes of galaxy development. To release this information we apply a novel, nonparametric statistical technique, identifying four components present in the distribution of galaxy H$\alpha$ emission-line equivalent-widths. We interpret these components as passive, star-forming, and two varieties of active galactic nuclei. Independent of this interpretation, the properties of each component are remarkably constant as a function of environment. Only their relative proportions display substantial variation. The galaxy population thus appears to comprise distinct components which are individually independent of environment, with galaxies rapidly transitioning between components as they move into denser environments.

With the Navy Prototype Optical Interferometer (NPOI), the binary system Theta 1 Orionis C, the most massive member of the Trapezium, was spatially resolved over a time period extending from February 2006 to March 2007. The data show significant orbital motion over the 14 months, and, after combining the NPOI data with previous measurements of the system from the literature, the observations span 10 years of the orbit. Our results indicate that the secondary did not experience an unusually close periastron passage this year, in contradiction to the prediction of a recently published, highly eccentric ~11 year orbit. Future observations of this source will be required to improve the orbital solution. Possible implications of the results in terms of system distance are discussed, although a main conclusion of this work is that a definitive orbit solution will require more time to obtain sufficient phase coverage, and that the interaction effects expected at periastron did not occur in 2007.

We present a new technique to compute simultaneously valid confidence intervals for a set of model parameters. We apply our method to the Wilkinson Microwave Anisotropy Probe's (WMAP) Cosmic Microwave Background (CMB) data, exploring a seven dimensional space (tau, Omega_DE, Omega_M, omega_DM, omega_B, f_nu, n_s). We find two distinct regions-of-interest: the standard Concordance Model, and a region with large values of omega_DM, omega_B and H_0. This second peak in parameter space can be rejected by applying a constraint (or a prior) on the allowable values of the Hubble constant. Our new technique uses a non-parametric fit to the data, along with a frequentist approach and a smart search algorithm to map out a statistical confidence surface. The result is a confidence ``ball'': a set of parameter values that contains the true value with probability at least 1-alpha. Our algorithm performs a role similar to the often used Markov Chain Monte Carlo (MCMC), which samples from the posterior probability function in order to provide Bayesian credible intervals on the parameters. While the MCMC approach samples densely around a peak in the posterior, our new technique allows cosmologists to perform efficient analyses around any regions of interest: e.g., the peak itself, or, possibly more importantly, the 1-alpha confidence surface.

We outline our first steps towards marrying two new and emerging technologies; the Virtual Observatory (e.g, AstroGrid) and the computational grid. We discuss the construction of VOTechBroker, which is a modular software tool designed to abstract the tasks of submission and management of a large number of computational jobs to a distributed computer system. The broker will also interact with the AstroGrid workflow and MySpace environments. We present our planned usage of the VOTechBroker in computing a huge number of n-point correlation functions from the SDSS, as well as fitting over a million CMBfast models to the WMAP data.

There is a growing need for massive computational resources for the analysis of new astronomical datasets. To tackle this problem, we present here our first steps towards marrying two new and emerging technologies; the Virtual Observatory (e.g, AstroGrid) and the computational grid (e.g. TeraGrid, COSMOS etc.). We discuss the construction of VOTechBroker, which is a modular software tool designed to abstract the tasks of submission and management of a large number of computational jobs to a distributed computer system. The broker will also interact with the AstroGrid workflow and MySpace environments. We discuss our planned usages of the VOTechBroker in computing a huge number of n-point correlation functions from the SDSS data and massive model-fitting of millions of CMBfast models to WMAP data. We also discuss other applications including the determination of the XMM Cluster Survey selection function and the construction of new WMAP maps.

We analyze first-year data of WMAP to determine the significance of asymmetry in summed power between arbitrarily defined opposite hemispheres, using maps that we create ourselves with software developed independently of the WMAP team. We find that over the multipole range l=[2,64], the significance of asymmetry is ~ 10^-4, a value insensitive to both frequency and power spectrum. We determine the smallest multipole ranges exhibiting significant asymmetry, and find twelve, including l=[2,3] and [6,7], for which the significance -> 0. In these ranges there is an improbable association between the direction of maximum significance and the ecliptic plane (p ~ 0.01). Also, contours of least significance follow great circles inclined relative to the ecliptic at the largest scales. The great circle for l=[2,3] passes over previously reported preferred axes and is insensitive to frequency, while the great circle for l=[6,7] is aligned with the ecliptic poles. We examine how changing map-making parameters affects asymmetry, and find that at large scales, it is rendered insignificant if the magnitude of the WMAP dipole vector is increased by approximately 1-3 sigma (or 2-6 km/s). While confirmation of this result would require data recalibration, such a systematic change would be consistent with observations of frequency-independent asymmetry. We conclude that the use of an incorrect dipole vector, in combination with a systematic or foreground process associated with the ecliptic, may help to explain the observed asymmetry.

We report on Keck Interferometer observations of the double-lined binary (B) component of the quadruple pre-main sequence (PMS) system HD 98800. With these interferometric observations combined with astrometric measurements made by the Hubble Space Telescope Fine Guidance Sensors (FGS), and published radial velocity observations we have estimated preliminary visual and physical orbits of the HD 98800 B subsystem. Our orbit model calls for an inclination of 66.8 $\pm$ 3.2 deg, and allows us to infer the masses and luminosities of the individual components. In particular we find component masses of 0.699 $\pm$ 0.064 and 0.582 $\pm$ 0.051 M$_{\sun}$ for the Ba (primary) and Bb (secondary) components respectively. Modeling of the component SEDs finds temperatures and luminosities in agreement with previous studies, and coupled with the component mass estimates allows for comparison with PMS models in the low-mass regime with few empirical constraints. Solar abundance models seem to under-predict the inferred component temperatures and luminosities, while assuming slightly sub-solar abundances bring the models and observations into better agreement. The present preliminary orbit does not yet place significant constraints on existing pre-main sequence stellar models, but prospects for additional observations improving the orbit model and component parameters are very good.

The Cosmic Microwave Background (CMB), which permeates the entire Universe, is the radiation left over from just 380,000 years after the Big Bang. On very large scales, the CMB radiation field is smooth and isotropic, but the existence of structure in the Universe - stars, galaxies, clusters of galaxies - suggests that the field should fluctuate on smaller scales. Recent observations, from the Cosmic Microwave Background Explorer to the Wilkinson Microwave Anisotropy Project, have strikingly confirmed this prediction. CMB fluctuations provide clues to the Universe's structure and composition shortly after the Big Bang that are critical for testing cosmological models. For example, CMB data can be used to determine what portion of the Universe is composed of ordinary matter versus the mysterious dark matter and dark energy. To this end, cosmologists usually summarize the fluctuations by the power spectrum, which gives the variance as a function of angular frequency. The spectrum's shape, and in particular the location and height of its peaks, relates directly to the parameters in the cosmological models. Thus, a critical statistical question is how accurately can these peaks be estimated. We use recently developed techniques to construct a nonparametric confidence set for the unknown CMB spectrum. Our estimated spectrum, based on minimal assumptions, closely matches the model-based estimates used by cosmologists, but we can make a wide range of additional inferences. We apply these techniques to test various models and to extract confidence intervals on cosmological parameters of interest. Our analysis shows that, even without parametric assumptions, the first peak is resolved accurately with current data but that the second and third peaks are not.

We report the detections of low-mass companions to five M and L dwarfs within 10 pc of the Sun using the Hubble Space Telescope (HST) Near-Infrared Camera and Multi-Object Spectrometer (NICMOS). Follow-up observations using the HST Advanced Camera for Surveys (ACS) and Fine Guidance Sensor 1r (FGS1r) confirm our NICMOS discoveries of companions to the L4.5 dwarf GJ 1001B (LHS 102B) and the M5 dwarf LHS 224, respectively. Images obtained with the Astrophysical Research Consortium 3.5 m telescope at Apache Point Observatory verify our discovery of a companion to the M3 dwarf G 239-25. Our NICMOS images confirm the previously suspected duplicity of the M3 dwarfs GJ 54 and GJ 84. The components of GJ 1001BC and LHS 224AB have nearly equal luminosities in all the ACS and/or NICMOS bandpasses in which they were observed. The magnitudes of GJ 54A and B in one FGS1r bandpass and four NICMOS bandpasses differ by <~ 1. GJ 84B and G 239-25B are ~4 magnitudes fainter than their M3 companions in the NICMOS bandpasses. We estimate from our NICMOS photometry that GJ 84B and G 239-25B have spectral types M7 and M8, respectively, and masses near the lower limit for sustained hydrogen burning. The apparent separations of GJ 1001BC, GJ 54AB, and LHS 224AB suggest orbital periods less than 5 yr. These binary dwarfs are ideal candidates for continued astrometric monitoring that will quickly yield accurate dynamic masses needed to constrain the mass-luminosity relation for low-mass stars and brown dwarfs.

Many fundamental statistical methods have become critical tools for scientific data analysis yet do not scale tractably to modern large datasets. This paper will describe very recent algorithms based on computational geometry which have dramatically reduced the computational complexity of 1) kernel density estimation (which also extends to nonparametric regression, classification, and clustering), and 2) the n-point correlation function for arbitrary n. These new multi-tree methods typically yield orders of magnitude in speedup over the previous state of the art for similar accuracy, making millions of data points tractable on desktop workstations for the first time.

The Deep Ecliptic Survey is a project whose goal is to survey a large area of the near-ecliptic region to a faint limiting magnitude (R ~ 24) in search of objects in the outer solar system. We are collecting a large homogeneous data sample from the Kitt Peak Mayall 4-m and Cerro Tololo Blanco 4-m telescopes with the Mosaic prime-focus CCD cameras. Our goal is to collect a sample of 500 objects with good orbits to further our understanding of the dynamical structure of the outer solar system. This survey has been in progress since 1998 and is responsible for 272 designated discoveries as of March 2003. We summarize our techniques, highlight recent results, and describe publically available resources.

We review ongoing efforts to identify occupants of mean-motion resonances (MMRs) and collisional families in the Edgeworth-Kuiper belt. Direct integrations of trajectories of Kuiper belt objects (KBOs) reveal the 1:1 (Trojan), 5:4, 4:3, 3:2 (Plutino), 5:3, 7:4, 9:5, 2:1 (Twotino), and 5:2 MMRs to be inhabited. Apart from the Trojan, resonant KBOs typically have large orbital eccentricities and inclinations. The observed pattern of resonance occupation is consistent with resonant capture and adiabatic excitation by a migratory Neptune; however, the dynamically cold initial conditions prior to resonance sweeping that are typically assumed by migration simulations are probably inadequate. Given the dynamically hot residents of the 5:2 MMR and the substantial inclinations observed in all exterior MMRs, a fraction of the primordial belt was likely dynamically pre-heated prior to resonance sweeping. A pre-heated population may have arisen as Neptune gravitationally scattered objects into trans-Neptunian space. The spatial distribution of Twotinos offers a unique diagnostic of Neptune's migration history. The Neptunian Trojan population may rival the Jovian Trojan population, and the former's existence is argued to rule out violent orbital histories for Neptune. Finally, lowest-order secular theory is applied to several hundred non-resonant KBOs with well-measured orbits to update proposals of collisional families. No convincing family is detected.

We present an absolute parallax and relative proper motion for the central star of the planetary nebula NGC 6853 (The Dumbell). We obtain these with astrometric data from FGS 3, a white-light interferometer on {\it HST}. Spectral classifications and VRIJHKT$_2$M and DDO51 photometry of the stars comprising the astrometric reference frame provide spectrophotometric estimates of their absolute parallaxes. Introducing these into our model as observations with error, we find $\pi_{abs} = 2.10 \pm 0.48$ mas for the DAO central star of NGC 6853. A weighted average with a previous ground-based USNO determination yields $\pi_{abs} = 2.40 \pm 0.32$. We assume that the extinction suffered by the reference stars nearest (in angular separation and distance) to the central star is the same as for the central star. Correcting for color differences, we find $<$A$_V>$ = 0.30 $ \pm $ 0.06 for the central star, hence, an absolute magnitude M$_V = 5.48^{-0.16}_{+0.15}$. A recent determination of the central star effective temperature aided in estimating the central star radius, $R_{*}=0.055 \pm 0.02R_{\sun}$, a star that may be descending to the white dwarf cooling track.

As part of our ongoing Deep Ecliptic Survey (DES) of the Kuiper belt, we report on the occupation of the 1:1 (Trojan), 4:3, 3:2, 7:4, 2:1, and 5:2 Neptunian mean-motion resonances (MMRs). The occupation of the 1:1 and 5:2 MMRs is not easily understood within the standard model of resonance sweeping by a migratory Neptune over an initially dynamically cold belt. Our dynamically hot, 5:2 resonant objects can librate with modest amplitudes of 90 deg within the resonance for at least 1 Gyr. Their trajectories cannot be explained by close encounters with Neptune alone, given the latter's current orbit. The hot orbits of such 5:2 resonant KBOs, unlike hot orbits of previously known resonant KBOs, may imply that these objects were pre-heated to large inclination and large eccentricity prior to resonance capture by a migratory Neptune. Our first discovered Neptunian Trojan, 2001QR322, may not owe its existence to Neptune's migration at all. The trajectory of 2001QR322 is remarkably stable; the object can undergo tadpole-type libration about Neptune's leading Lagrange (L4) point for at least 1 Gyr with a libration amplitude of 24 deg. Trojan capture probably occurred while Neptune accreted the bulk of its mass. For an assumed albedo of 12--4%, our Trojan is 130--230 km in diameter. Model-dependent estimates place the total number of Neptune Trojans resembling 2001QR322 at 20--60. Their existence might rule out violent orbital histories for Neptune.

We present an absolute parallax and relative proper motion for the fundamental distance scale calibrator, delta Cep. We obtain these with astrometric data from FGS 3, a white-light interferometer on HST. Utilizing spectrophotometric estimates of the absolute parallaxes of our astrometric reference stars and constraining delta Cep and reference star HD 213307 to belong to the same association (Cep OB6, de Zeeuw et al. 1999), we find pi_{abs} = 3.66 +/- 0.15 mas. The larger than typical astrometric residuals for the nearby astrometric reference star HD 213307 are found to satisfy Keplerian motion with P = 1.07 +/- 0.02 years, a perturbation and period that could be due to a F0V companion ~7 mas distant from and ~4 magnitudes fainter than the primary. Spectral classifications and VRIJHKT$_2$M and DDO51 photometry of the astrometric reference frame surrounding delta Cep indicate that field extinction is high and variable along this line of sight. However the extinction suffered by the reference star nearest (in angular separation and distance) to delta Cep, HD 213307, is lower and nearly the same as for delta Cep. Correcting for color differences, we find <A_V> = 0.23 +/- 0.03 for delta Cep, hence, an absolute magnitude M_V = -3.47 +/- 0.10. Adopting an average V magnitude, <V> = 15.03 +/- 0.03, for Cepheids with log P = 0.73 in the LMC from Udalski et al. (1999), we find a V-band distance modulus for the LMC, m-M = 18.50 +/- 0.13 or, 18.58 +/- 0.15, where the latter value results from a highly uncertain metallicity correction (Freedman et al. 2001). These agree with our previous RR Lyr HST parallax-based determination of the distance modulus of the LMC.

We discuss non-parametric density estimation and regression for astrophysics problems. In particular, we show how to compute non-parametric confidence intervals for the location and size of peaks of a function. We illustrate these ideas with recent data on the Cosmic Microwave Background. We also briefly discuss non-parametric Bayesian inference.

We examine Cosmic Microwave Background (CMB) temperature power spectra from the BOOMERANG, MAXIMA, and DASI experiments. We non-parametrically estimate the true power spectrum with no model assumptions. This is a significant departure from previous research which used either cosmological models or some other parameterized form (e.g. parabolic fits). Our non-parametric estimate is practically indistinguishable from the best fit cosmological model, thus lending independent support to the underlying physics that governs these models. We also generate a confidence set for the non-parametric fit and extract confidence intervals for the numbers, locations, and heights of peaks and the successive peak-to-peak height ratios. At the 95%, 68%, and 40% confidence levels, we find functions that fit the data with one, two, and three peaks respectively (0 <= l <= 1100). Therefore, the current data prefer two peaks at the 1 sigma level. However, we also rule out a constant temperature function at the >8 sigma level. If we assume that there are three peaks in the data, we find their locations to be within l_1 = (118,300), l_2 = (377,650), and l_3 = (597,900). We find the ratio of the first peak-height to the second (Delta T_1)/(Delta T_2)^2= (1.06, 4.27) and the second to the third (Delta T_2)/(Delta T_3)^2= (0.41, 2.5). All measurements are for 95% confidence. If the standard errors on the temperature measurements were reduced to a third of what they are currently, as we expect to be achieved by the MAP and Planck CMB experiments, we could eliminate two-peak models at the 95% confidence limit. The non-parametric methodology discussed in this paper has many astrophysical applications.

The False Discovery Rate (FDR) method has recently been described by Miller et al (2001), along with several examples of astrophysical applications. FDR is a new statistical procedure due to Benjamini and Hochberg (1995) for controlling the fraction of false positives when performing multiple hypothesis testing. The importance of this method to source detection algorithms is immediately clear. To explore the possibilities offered we have developed a new task for performing source detection in radio-telescope images, Sfind 2.0, which implements FDR. We compare Sfind 2.0 with two other source detection and measurement tasks, Imsad and SExtractor, and comment on several issues arising from the nature of the correlation between nearby pixels and the necessary assumption of the null hypothesis. The strong suggestion is made that implementing FDR as a threshold defining method in other existing source-detection tasks is easy and worthwhile. We show that the constraint on the fraction of false detections as specified by FDR holds true even for highly correlated and realistic images. For the detection of true sources, which are complex combinations of source-pixels, this constraint appears to be somewhat less strict. It is still reliable enough, however, for a priori estimates of the fraction of false source detections to be robust and realistic.

We present an absolute parallax and relative proper motion for the fundamental distance scale calibrator, RR Lyr. We obtain these with astrometric data from FGS 3, a white-light interferometer on HST. We find $\pi_{abs} = 3.82 \pm 0.2$ mas. Spectral classifications and VRIJHKT$_2$M and DDO51 photometry of the astrometric reference frame surrounding RR Lyr indicate that field extinction is low along this line of sight. We estimate <A_V>=0.07\pm0.03 for these reference stars. The extinction suffered by RR Lyr becomes one of the dominant contributors to the uncertainty in its absolute magnitude. Adopting the average field absorption, <A_V>=0.07 \pm 0.03, we obtain M_V^{RR} = 0.61 ^{-0.11}_{+0.10}. This provides a distance modulus for the LMC, m-M = 18.38 - 18.53^{-0.11}_{+0.10} with the average extinction-corrected magnitude of RR Lyr variables in the LMC, <V(RR)>, remaining a significant uncertainty. We compare this result to more than 80 other determinations of the distance modulus of the LMC.

I present here a review of past and present multi-disciplinary research of the Pittsburgh Computational AstroStatistics (PiCA) group. This group is dedicated to developing fast and efficient statistical algorithms for analysing huge astronomical data sources. I begin with a short review of multi-resolutional kd-trees which are the building blocks for many of our algorithms. For example, quick range queries and fast n-point correlation functions. I will present new results from the use of Mixture Models (Connolly et al. 2000) in density estimation of multi-color data from the Sloan Digital Sky Survey (SDSS). Specifically, the selection of quasars and the automated identification of X-ray sources. I will also present a brief overview of the False Discovery Rate (FDR) procedure (Miller et al. 2001a) and show how it has been used in the detection of ``Baryon Wiggles'' in the local galaxy power spectrum and source identification in radio data. Finally, I will look forward to new research on an automated Bayes Network anomaly detector and the possible use of the Locally Linear Embedding algorithm (LLE; Roweis & Saul 2000) for spectral classification of SDSS spectra.

TV Columbae (TV Col) is a 13th magnitude Intermediate Polar (IP) Cataclysmic Variable (CV), with multiple periods found in the light curves. Past estimates predicted a distance of 400 parsec to greater than 500 parsec. Recently completed Hubble Space Telescope (HST) Fine Guidance Sensor (FGS) interferometric observations allow us to determine the first trigonometric parallax to TV Col. This determination puts the distance of TV Col at 368 -15+17 parsecs. CD-32 2376, a 10th magnitude Tycho Catalog star, is a reference star in the TV Col frame. We find a distance of $127.7 -1+1 parsecs.

The False Discovery Rate (FDR) is a new statistical procedure to control the number of mistakes made when performing multiple hypothesis tests, i.e. when comparing many data against a given model hypothesis. The key advantage of FDR is that it allows one to a priori control the average fraction of false rejections made (when comparing to the null hypothesis) over the total number of rejections performed. We compare FDR to the standard procedure of rejecting all tests that do not match the null hypothesis above some arbitrarily chosen confidence limit, e.g. 2 sigma, or at the 95% confidence level. When using FDR, we find a similar rate of correct detections, but with significantly fewer false detections. Moreover, the FDR procedure is quick and easy to compute and can be trivially adapted to work with correlated data. The purpose of this paper is to introduce the FDR procedure to the astrophysics community. We illustrate the power of FDR through several astronomical examples, including the detection of features against a smooth one-dimensional function, e.g. seeing the ``baryon wiggles'' in a power spectrum of matter fluctuations, and source pixel detection in imaging data. In this era of large datasets and high precision measurements, FDR provides the means to adaptively control a scientifically meaningful quantity -- the number of false discoveries made when conducting multiple hypothesis tests.

We present here a new algorithm for the fast computation of N-point correlation functions in large astronomical data sets. The algorithm is based on kdtrees which are decorated with cached sufficient statistics thus allowing for orders of magnitude speed-ups over the naive non-tree-based implementation of correlation functions. We further discuss the use of controlled approximations within the computation which allows for further acceleration. In summary, our algorithm now makes it possible to compute exact, all-pairs, measurements of the 2, 3 and 4-point correlation functions for cosmological data sets like the Sloan Digital Sky Survey (SDSS; York et al. 2000) and the next generation of Cosmic Microwave Background experiments (see Szapudi et al. 2000).

We present an analysis of astrometric data from FGS 3, a white-light interferometer on {\it HST}, and of radial velocity data from two ground-based campaigns. We model the astrometric and radial velocity measurements simultaneously to obtain parallax, proper motion and component masses for Wolf 1062 = Gl 748 AB (M3.5V). To derive the mass fraction, we relate FGS 3 fringe scanning observations of the science target to a reference frame provided by fringe tracking observations of a surrounding star field. We obtain an absolute parallax $\pi_{abs} = 98.0 \pm 0.4$ milliseconds of arc, yielding ${\cal M}_A = 0.379 \pm 0.005{\cal M}_{\sun}$ and ${\cal M}_B= 0.192 \pm 0.003 {\cal M}_{\sun}$, high quality component masses with errors of only 1.5%.

We outline here the next generation of cluster-finding algorithms. We show how advances in Computer Science and Statistics have helped develop robust, fast algorithms for finding clusters of galaxies in large multi-dimensional astronomical databases like the Sloan Digital Sky Survey (SDSS). Specifically, this paper presents four new advances: (1) A new semi-parametric algorithm - nicknamed ``C4'' - for jointly finding clusters of galaxies in the SDSS and ROSAT All-Sky Survey databases; (2) The introduction of the False Discovery Rate into Astronomy; (3) The role of kernel shape in optimizing cluster detection; (4) A new determination of the X-ray Cluster Luminosity Function which has bearing on the existence of a ``deficit'' of high redshift, high luminosity clusters. This research is part of our ``Computational AstroStatistics'' collaboration (see Nichol et al. 2000) and the algorithms and techniques discussed herein will form part of the ``Virtual Observatory'' analysis toolkit.

In this paper, we outline the use of Mixture Models in density estimation of large astronomical databases. This method of density estimation has been known in Statistics for some time but has not been implemented because of the large computational cost. Herein, we detail an implementation of the Mixture Model density estimation based on multi-resolutional KD-trees which makes this statistical technique into a computationally tractable problem. We provide the theoretical and experimental background for using a mixture model of Gaussians based on the Expectation Maximization (EM) Algorithm. Applying these analyses to simulated data sets we show that the EM algorithm - using the AIC penalized likelihood to score the fit - out-performs the best kernel density estimate of the distribution while requiring no ``fine--tuning'' of the input algorithm parameters. We find that EM can accurately recover the underlying density distribution from point processes thus providing an efficient adaptive smoothing method for astronomical source catalogs. To demonstrate the general application of this statistic to astrophysical problems we consider two cases of density estimation: the clustering of galaxies in redshift space and the clustering of stars in color space. From these data we show that EM provides an adaptive smoothing of the distribution of galaxies in redshift space (describing accurately both the small and large-scale features within the data) and a means of identifying outliers in multi-dimensional color-color space (e.g. for the identification of high redshift QSOs). Automated tools such as those based on the EM algorithm will be needed in the analysis of the next generation of astronomical catalogs (2MASS, FIRST, PLANCK, SDSS) and ultimately in in the development of the National Virtual Observatory.

We present initial results on the use of Mixture Models for density estimation in large astronomical databases. We provide herein both the theoretical and experimental background for using a mixture model of Gaussians based on the Expectation Maximization (EM) Algorithm. Applying these analyses to simulated data sets we show that the EM algorithm - using the both the AIC & BIC penalized likelihood to score the fit - can out-perform the best kernel density estimate of the distribution while requiring no ``fine-tuning'' of the input algorithm parameters. We find that EM can accurately recover the underlying density distribution from point processes thus providing an efficient adaptive smoothing method for astronomical source catalogs. To demonstrate the general application of this statistic to astrophysical problems we consider two cases of density estimation; the clustering of galaxies in redshift space and the clustering of stars in color space. From these data we show that EM provides an adaptive smoothing of the distribution of galaxies in redshift space (describing accurately both the small and large-scale features within the data) and a means of identifying outliers in multi-dimensional color-color space (e.g. for the identification of high redshift QSOs). Automated tools such as those based on the EM algorithm will be needed in the analysis of the next generation of astronomical catalogs (2MASS, FIRST, PLANCK, SDSS) and ultimately in the development of the National Virtual Observatory.

With fourteen epochs of fringe tracking data spanning 1.7y from Fine Guidance Sensor 3 we have obtained a parallax (pi_abs=113.1 +- 0.3 mas) and perturbation orbit for Gl 791.2A. Contemporaneous fringe scanning observations yield only three clear detections of the secondary on both interferometer axes. They provide a mean component magnitude difference, Delta V = 3.27 +- 0.10. The period (P = 1.4731 yr) from the perturbation orbit and the semi-major axis (a = 0.963 +- 0.007 AU) from the measured component separations with our parallax provide a total system mass M_A + M_B = 0.412 +- 0.009 M_sun. Component masses are M_A=0.286 +- 0.006 M_sun and M_B = 0.126 +- 0.003 M_sun. Gl 791.2A and B are placed in a sparsely populated region of the lower main sequence mass-luminosity relation where they help define the relation because the masses have been determined to high accuracy, with errors of only 2%.

With HST FGS 3 we have determined a parallax for the white dwarf - M dwarf interacting binary, Feige 24. The white dwarf (DA) component has an effective temperature, T_eff of approximately 56,000 K degrees. A weighted average with past parallax determinations (pi_abs = 14.6 +- 0.4 milliseconds of arc) narrows the range of possible radius values, compared to past estimates. We obtain R_DA = 0.0185 +- 0.0008 R_solar with uncertainty in the temperature and bolometric correction the dominant contributors to the error. FGS photometry provides a light curve entirely consistent with reflection effects. A recently refined model Mass-Luminosity Relation (Baraffe et al. 1998) for low mass stars provides a mass estimate for the M dwarf companion, M_dM = 0.37 +- 0.20 M_solar, where the mass range is due to metallicity and age uncertainties. Radial velocities from Vennes and Thorstensen (1994) provide a mass ratio from which we obtain M_DA = 0.49+0.19-0.05 ~ M_solar. Independently, our radius and recent log_g determinations yield 0.44 < M_DA < 0.47 M_solar. In each case the minimum DA mass is that derived by Vennes & Thorstensen from their radial velocities and Keplerian circular orbits with inclination <= 90 degrees. Locating Feige 24 on an M - R plane suggests a carbon core. Our radius and these mass estimates yield a gamma_grav inconsistent with that derived by Vennes & Thorstensen. We speculate on the nature of a third component whose existence would resolve the discrepancy.

We present an improved orbit determination for the visual pair in the M-dwarf triple system Gliese 866 that is based on new speckle-interferometric and HST observations. The system mass is M = 0.34 +- 0.03 solar masses. The masses of the components derived using the mass-luminosity relation are consistent with this mass sum. All three components of Gliese 866 seem to have masses not far from the hydrogen burning mass limit.

We report on a sub-stellar companion search utilizing interferometric fringe-tracking astrometry acquired with Fine Guidance Sensor 3 (FGS 3) on the Hubble Space Telescope. Our targets were Proxima Centauri and Barnard's Star. We obtain absolute parallax values for Proxima Cen pi_{abs} = 0.7687 arcsecond and for Barnard's Star pi_{abs} = 0.5454 arcsecond. Once low-amplitude instrumental systematic errors are identified and removed, our companion detection sensitivity is less than or equal to one Jupiter mass for periods longer than 60 days for Proxima Cen. Between the astrometry and the radial velocity results we exclude all companions with M > 0.8M_{Jup} for the range of periods 1 < P < 1000 days. For Barnard's Star our companion detection sensitivity is less than or equal to one Jupiter mass for periods long er than 150 days. Our null results for Barnard's Star are consistent with those of Gatewood (1995).

RW Triangulum (RW Tri) is a 13th magnitude Nova-like Cataclysmic Variable star with an orbital period of 0.2319 days (5.56 hours). Infrared observations of RW Tri indicate that its secondary is most likely a late K-dwarf. Past analyses predicted a distance of 270 parsec, derived from a black-body fit to spectrum of the central part of the disk. Recently completed Hubble Space Telescope Fine Guidance Sensor interferometric observations allow us to determine the first trigonometric parallax to RW Tri. This determination puts the distance of RW Tri at 341, one of the most distant objects with a direct parallax measurement.

We have observed Proxima Centauri and Barnard's Star with Hubble Space Telescope Fine Guidance Sensor 3. Proxima Centauri exhibits small-amplitude, periodic photometric variations. Once several sources of systematic photometric error are corrected, we obtain 2 milli-magnitude internal photometric precision. We identify two distinct behavior modes over the past four years: higher amplitude, longer period; smaller amplitude, shorter period. Within the errors one period (P ~ 83d) is twice the other. Barnard's Star shows very weak evidence for periodicity on a timescale of approximately 130 days. If we interpret these periodic phenomena as rotational modulation of star spots, we identify three discrete spots on Proxima Cen and possibly one spot on Barnard's Star. We find that the disturbances change significantly on time scales as short as one rotation period.