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.

We obtained thorough photometric observations of two binary near-Earth asteroids (66391) Moshup = 1999 KW4 and (88710) 2001 SL9 from June 2000 to June 2019 and derived physical and dynamical properties of the binary systems. We found that the data for 1999 KW4 are inconsistent with a constant orbital period and we obtained unique solution with a quadratic drift of the mean anomaly of the satellite of -0.65 +/- 0.16 deg/yr2 (all quoted uncertainties are 3sigma). This means that the semimajor axis of the mutual orbit of the components of this binary system increases in time with a mean rate of 1.2 +/- 0.3 cm/yr. The data for 2001 SL9 are also inconsistent with a constant orbital period and we obtained two solutions for the quadratic drift of the mean anomaly: 2.8 +/- 0.2 and 5.2 +/- 0.2 deg/yr2, implying that the semimajor axis of the mutual orbit of the components decreases in time with a mean rate of -2.8 +/- 0.2 or -5.1 +/- 0.2 cm/yr for the two solutions, respectively. The expanding orbit of 1999 KW4 may be explained by mutual tides interplaying with binary YORP (BYORP) effect (McMahon and Scheeres, 2010). However, a modeling of the BYORP drift using radar-derived shapes of the binary components predicted a much higher value of the orbital drift than the observed one. It suggests that either the radar-derived shape model of the secondary is inadequate for computing the BYORP effect, or the present theory of BYORP overestimates it. It is possible that the BYORP coefficient has instead an opposite sign than predicted; in that case, the system may be moving into an equilibrium between the BYORP and the tides. In the case of 2001 SL9, the BYORP effect is the only known physical mechanism that can cause the inward drift of its mutual orbit.

We studied 93 asteroid pairs. We estimated times elapsed since separation of pair members that are between 7*10^3 and a few 10^6 yr. We derived the rotation periods for all the primaries and a sample of secondaries. We derived the absolute magnitude differences of the asteroid pairs that provide their mass ratios. We refined their WISE geometric albedos and estimated their taxonomic classifications. For 17 pairs, we determined their pole positions. In 2 pairs where we obtained the spin poles for both components, we saw the same sense of rotation for both components and constrained the angles between their original spin vectors at the time of their separation. We found that the primaries of 13 pairs are actually binary or triple systems, i.e., they have one or two bound secondaries (satellites). As by-product, we found 3 new young asteroid clusters (each of them consisting of three known asteroids on highly similar orbits). We compared the obtained asteroid pair data with theoretical predictions and discussed their implications. We found that 86 of the 93 studied pairs follow the trend of primary rotation period vs mass ratio that was found by Pravec et al. (2010). Of the 7 outliers, 3 appear insignificant (may be due to our uncertain or incomplete knowledge), but 4 are high mass ratio pairs that were unpredicted by the theory of asteroid pair formation by rotational fission. We discuss a (remotely) possible way that they could be created by rotational fission of flattened parent bodies followed by re-shaping of the formed components. The 13 pairs with binary primaries are particularly interesting systems that place important constraints on formation and evolution of asteroid pairs. We present two hypotheses for their formation: The pairs having both bound and unbound secondaries could be `failed asteroid clusters', or they could be formed by a cascade primary spin fission process.

Parallaxes for 331 classical Cepheids, 31 Type II Cepheids and 364 RR Lyrae stars in common between Gaia and the Hipparcos and Tycho-2 catalogues are published in Gaia Data Release 1 (DR1) as part of the Tycho-Gaia Astrometric Solution (TGAS). In order to test these first parallax measurements of the primary standard candles of the cosmological distance ladder, that involve astrometry collected by Gaia during the initial 14 months of science operation, we compared them with literature estimates and derived new period-luminosity ($PL$), period-Wesenheit ($PW$) relations for classical and Type II Cepheids and infrared $PL$, $PL$-metallicity ($PLZ$) and optical luminosity-metallicity ($M_V$-[Fe/H]) relations for the RR Lyrae stars, with zero points based on TGAS. The new relations were computed using multi-band ($V,I,J,K_{\mathrm{s}},W_{1}$) photometry and spectroscopic metal abundances available in the literature, and applying three alternative approaches: (i) by linear least squares fitting the absolute magnitudes inferred from direct transformation of the TGAS parallaxes, (ii) by adopting astrometric-based luminosities, and (iii) using a Bayesian fitting approach. TGAS parallaxes bring a significant added value to the previous Hipparcos estimates. The relations presented in this paper represent first Gaia-calibrated relations and form a "work-in-progress" milestone report in the wait for Gaia-only parallaxes of which a first solution will become available with Gaia's Data Release 2 (DR2) in 2018.

Context. The first Gaia Data Release contains the Tycho-Gaia Astrometric Solution (TGAS). This is a subset of about 2 million stars for which, besides the position and photometry, the proper motion and parallax are calculated using Hipparcos and Tycho-2 positions in 1991.25 as prior information. Aims. We investigate the scientific potential and limitations of the TGAS component by means of the astrometric data for open clusters. Methods. Mean cluster parallax and proper motion values are derived taking into account the error correlations within the astrometric solutions for individual stars, an estimate of the internal velocity dispersion in the cluster, and, where relevant, the effects of the depth of the cluster along the line of sight. Internal consistency of the TGAS data is assessed. Results. Values given for standard uncertainties are still inaccurate and may lead to unrealistic unit-weight standard deviations of least squares solutions for cluster parameters. Reconstructed mean cluster parallax and proper motion values are generally in very good agreement with earlier Hipparcos-based determination, although the Gaia mean parallax for the Pleiades is a significant exception. We have no current explanation for that discrepancy. Most clusters are observed to extend to nearly 15 pc from the cluster centre, and it will be up to future Gaia releases to establish whether those potential cluster-member stars are still dynamically bound to the clusters. Conclusions. The Gaia DR1 provides the means to examine open clusters far beyond their more easily visible cores, and can provide membership assessments based on proper motions and parallaxes. A combined HR diagram shows the same features as observed before using the Hipparcos data, with clearly increased luminosities for older A and F dwarfs.

V471 Tauri, a white dwarf--red dwarf eclipsing binary in the Hyades, is well known for stimulating development of common envelope theory, whereby novae and other cataclysmic variables form from much wider binaries by catastrophic orbit shrinkage. Our evaluation of a recent imaging search that reported negative results for a much postulated third body shows that the object could have escaped detection or may have actually been seen. The balance of evidence continues to favor a brown dwarf companion about 12 AU from the eclipsing binary. A recently developed algorithm finds unified solutions from three datatypes. New radial velocities (RVs) of the red dwarf and BV RCIC light curves are solved simultaneously along with white dwarf and red dwarf RVs from the literature, uvby data, the MOST mission light curve, and 40 years of eclipse timings. Precision-based weighting is the key to proper information balance among the various datasets. Timewise variation of modeled starspots allows unified solution of multiple data eras. Light curve amplitudes strongly suggest decreasing spottedness from 1976 to about 1980, followed by approximately constant spot coverage from 1981 to 2005. An explanation is proposed for lack of noticeable variation in 1981 light curves, in terms of competition between spot and tidal variations. Photometric spectroscopic distance is estimated. The red dwarf mass comes out larger than normal for a K2V star, and even larger than adopted in several structure and evolution papers. An identified cause for this result is that much improved red dwarf RVs curves now exist.

We provide APASS photometry in the Landolt BV and Sloan g'r'i' bands for all the 425,743 stars included in the latest 4th RAVE Data Release. The internal accuracy of the APASS photometry of RAVE stars, expressed as error of the mean of data obtained and separately calibrated over a median of 4 distinct observing epochs and distributed between 2009 and 2013, is 0.013, 0.012, 0.012, 0.014 and 0.021 mag for B, V, g', r' and i' band, respectively. The equally high external accuracy of APASS photometry has been verified on secondary Landolt and Sloan photometric standard stars not involved in the APASS calibration process, and on a large body of literature data on field and cluster stars, confirming the absence of offsets and trends. Compared with the Carlsberg Meridian Catalog (CMC-15), APASS astrometry of RAVE stars is accurate to a median value of 0.098 arcsec. Brightness distribution functions for the RAVE stars have been derived in all bands. APASS photometry of RAVE stars, augmented by 2MASS JHK infrared data, has been chi2 fitted to a densely populated synthetic photometric library designed to widely explore in temperature, surface gravity, metallicity and reddening. Resulting Teff and E(B-V), computed over a range of options, are provided and discussed, and will be kept updated in response to future APASS and RAVE data releases. In the process it is found that the reddening caused by an homogeneous slab of dust, extending for 140 pc on either side of the Galactic plane and responsible for E(B-V,poles)=0.036 +/- 0.002 at the galactic poles, is a suitable approximation of the actual reddening encountered at Galactic latitudes |b|>=25 deg.

KOI-227, KOI-319 and KOI-884 are identified here as (at least) two planet systems. For KOI-319 and KOI-884, the observed Transit Timing Variations (TTVs) of the inner transiting planet are used to detect an outer non-transiting planet. The outer planet in KOI-884 is 2.6 Jupiter masses and has the orbital period just narrow of the 3:1 resonance with the inner planet (orbital period ratio 2.93). The distribution of parameters inferred from KOI-319.01's TTVs is bimodal with either a 1.6 Neptune-mass (M_N) planet wide of the 5:3 resonance (period 80.1 d) or a 1 Saturn-mass planet wide of the 7:3 resonance (period 109.2 d). The radial velocity measurements can be used in this case to determine which of these parameter modes is correct. KOI-227.01's TTVs with large 10 hour amplitude can be obtained for planetary-mass companions in various major resonances. Based on the Bayesian evidence, the current TTV data favor the outer 2:1 resonance with a companion mass 1.5 M_N, but this solution implies a very large density of KOI-227.01. The inner and outer 3:2 resonance solutions with sub-Neptune-mass companions are physically more plausible, but will need to be verified.

We present new radial velocities and photometry of the short-period Algol TT Herculis. Previous attempts to model the light curves of the system have met with limited success, primarily because of the lack of a reliable mass ratio. Our spectroscopic observations are the first to result in radial velocities for the secondary star, and thus provide a spectroscopic mass ratio. Simultaneous analysis of the radial velocities and new photometry shows that the system is a double contact binary, with a rapidly rotating primary that fills its limiting lobe.

The Transit Timing Variations (TTVs) can be used as a diagnostic of gravitational interactions between planets in a multi-planet system. Many Kepler Objects of Interest (KOIs) exhibit significant TTVs, but KOI-142.01 stands out among them with an unrivaled, 12-hour TTV amplitude. Here we report a thorough analysis of KOI-142.01's transits. We discover periodic Transit Duration Variations (TDVs) of KOI-142.01 that are nearly in phase with the observed TTVs. We show that KOI-142.01's TTVs and TDVs uniquely detect a non-transiting companion with a mass 0.7 that of Jupiter (KOI-142c). KOI-142.01's mass inferred from the transit variations is consistent with the measured transit depth, suggesting a Neptune class planet (KOI-142b). The orbital period ratio P_c/P_b=2.03 indicates that the two planets are just wide of the 2:1 resonance. The present dynamics of this system, characterized here in detail, can be used to test various formation theories that have been proposed to explain the near-resonant pairs of exoplanets.

We report the discovery and confirmation of a transiting circumbinary planet (PH1b) around KIC 4862625, an eclipsing binary in the Kepler field. The planet was discovered by volunteers searching the first six Quarters of publicly available Kepler data as part of the Planet Hunters citizen science project. Transits of the planet across the larger and brighter of the eclipsing stars are detectable by visual inspection every ~137 days, with seven transits identified in Quarters 1-11. The physical and orbital parameters of both the host stars and planet were obtained via a photometric-dynamical model, simultaneously fitting both the measured radial velocities and the Kepler light curve of KIC 4862625. The 6.18 +/- 0.17 Earth radii planet orbits outside the 20-day orbit of an eclipsing binary consisting of an F dwarf (1.734 +/- 0.044 Solar radii, 1.528 +/- 0.087 Solar masses) and M dwarf (0.378+/- 0.023 Solar radii, 0.408 +/- 0.024 Solar masses). For the planet, we find an upper mass limit of 169 Earth masses (0.531 Jupiter masses) at the 99.7% confidence level. With a radius and mass less than that of Jupiter, PH1b is well within the planetary regime. Outside the planet's orbit, at ~1000 AU,a previously unknown visual binary has been identified that is likely bound to the planetary system, making this the first known case of a quadruple star system with a transiting planet.

The spin state of small asteroids can change on a long timescale by the Yarkovsky-O'Keefe-Radzievskii-Paddack (YORP) effect, the net torque that arises from anisotropically scattered sunlight and proper thermal radiation from an irregularly-shaped asteroid. The secular change in the rotation period caused by the YORP effect can be detected by analysis of asteroid photometric lightcurves. We analyzed photometric lightcurves of near-Earth asteroids (1865) Cerberus, (2100) Ra-Shalom, and (3103) Eger with the aim to detect possible deviations from the constant rotation caused by the YORP effect. We carried out new photometric observations of the three asteroids, combined the new lightcurves with archived data, and used the lightcurve inversion method to model the asteroid shape, pole direction, and rotation rate. The YORP effect was modeled as a linear change in the rotation rate in time d\omega /dt. Values of d\omega/ dt derived from observations were compared with the values predicted by theory. We derived physical models for all three asteroids. We had to model Eger as a nonconvex body because the convex model failed to fit the lightcurves observed at high phase angles. We probably detected the acceleration of the rotation rate of Eger d\omega / dt = (1.4 +/- 0.6) x 10^{-8} rad/d (3\sigma error), which corresponds to a decrease in the rotation period by 4.2 ms/yr. The photometry of Cerberus and Ra-Shalom was consistent with a constant-period model, and no secular change in the spin rate was detected. We could only constrain maximum values of |d\omega / dt| < 8 x 10^{-9} rad/d for Cerberus, and |d\omega / dt| < 3 x 10^{-8} rad/d for Ra-Shalom.

We report on a BVRcIc survey of field W Ursae Majoris binary stars and present accurate colors for 606 systems that have been observed on at least three photometric nights from a robotic observatory in southern Arizona. Comparison with earlier photometry for a subset of the systems shows good agreement. We investigate two independent methods of determining the interstellar reddening, although both have limitations that can render them less effective than desired. A subset of 101 systems shows good agreement between the two reddening methods.

We present a preliminary spectroscopic analysis of the binary system V577Oph, observed during the summer of 2007 on the 2.6m NOT telescope on La Palma. We have obtained time series spectroscopic observations, which show clear binary motion as well as radial velocity variations due to pulsation in the primary star. By modelling the radial velocities we determine a full orbital solution of the system, which yields M_A sin^3 i = 1.562 +/- 0.012 M_solar and M_B sin^3 i = 1.461 +/- 0.020 M_solar. An estimate of inclination from photometry yields a primary mass of 1.6 M_solar. Using this derived mass, and the known pulsation frequency we can impose a lower limit of 1 Gyr on the age of the system, and constrain the parameters of the oscillation mode. We show that with further analysis of the spectra (extracting the atmospheric parameters), tighter constraints could be imposed on the age, metallicity and the mode parameters. This work emphasizes the power that a single pulsation frequency can have for constraining stellar parameters in an eclipsing binary system.

A distance is found for the W UMa type binary U Pegasi, with a newly modified version of the Wilson-Devinney program (W-D) that makes use of the direct distance estimation (DDE) algorithm. The reported distance of d = 123.6pc is an average based on solutions for B and V data and a primary star temperature of 5800K. Standardized light curves (not differential), radial velocities, and a spectroscopic primary star temperature are input to the pro- gram. Differential corrections were performed for each light curve band along with the velocities for two primary temperatures that span 100K. Log10d is a model parameter like many others that are adjustable in W-D. The eclipsing binary distance agrees with the Hipparcos parallax distance and is more precise.

New and previously published observations of the bright eclipsing binary VV Orionis are analyzed. We present new radial velocities and interstellar reddening measurements from high-resolution spectra of this detached, short-period (P=1.48 d) binary. We discuss the validity of prior claims for the existence of a third body and show that our new velocities and light curve solution cast doubt on them. The components of VV Ori are shown to be a B1 V primary with a mass $M{_1}=10.9 \pm 0.1 M_{\sun}$ and a radius $R_{1}=4.98 \pm 0.02 R_{\sun}$ and a B4.5 V secondary with a mass $M{_2}=4.09 \pm 0.05 M_{\sun}$ and a radius $R_{2}=2.41 \pm 0.01 R_{\sun}$.

Observations of Pluto and its solar-tidal stability zone were made using the Advanced Camera for Surveys' (ACS) Wide Field Channel (WFC) on the Hubble Space Telescope on UT 2005 May 15 and UT 2005 May 18. Two small satellites of Pluto, provisionally designated S/2005 P 1 and S/2005 P 2, were discovered, as discussed by Weaver et al. (2006) and Stern et al. (2006a). Confirming observations of the newly discovered moons were obtained using the ACS in the High Resolution Channel (HRC) mode on 2006 Feb 15 (Mutchler et al. 2006). Both sets of observations provide strong constraints on the existence of any additional satellites in the Pluto system. Based on the May 2005 observations using the ACS/WFC, we place a 90%-confidence lower limit of m_V = 26.8 (m_V = 27.4 for a 50%-confidence lower limit) on the magnitude of undiscovered satellites greater than 5" (1.1x10^5 km) from Pluto. Using the 2005 Feb 15 ACS/HRC observations we place 90%-confidence lower limits on the apparent magnitude of any additional satellites of m_V = 26.4 between 3"-5" (6.9x10^4-1.1x10^5 km) from Pluto, m_V = 25.7 between 1"-3" (2.3x10^4-6.9x10^4 km) from Pluto, and m_V = 24. between 0.3"-1" (6.9x10^3-2.3x10^4 km) from Pluto. The 90%-confidence magnitude limits translate into upper limits on the diameters of undiscovered satellites of 29 km outside of 5" from Pluto, 36 km between 3"-5" from Pluto, 49 km between 1"-3" from Pluto, and 115 km between 0.3"-1" for a comet-like albedo of p_V = 0.04. If potential satellites are assumed to have a Charon-like albedo of p_V = 0.38, the diameter limits are 9 km, 12 km, 16 km, and 37 km, respectively.

We present photometric and spectroscopic data on the early-type binary MP Centauri. The photometric data are analyzed simultaneously with radial velocities to derive preliminary absolute dimensions for the binary components. Analysis of the spectra shows that the stars rotate synchronously and that the line of sight to the system crosses two kinematically sharp and well-separated interstellar reddening sources. It is shown that MP Cen consists of a B3 primary with M_1 = 11.4 +/- 0.4 M_sun, R_1 = 7.7 +/- 0.1 R_sun and a lobe-filling B6-B7 secondary with M_2 = 4.4 +/- 0.2 M_sun, R_2 = 6.6 +/- 0.1 R_sun

We present new spectroscopic and photometric data on the early-type overcontact binary TU Muscae. The analysis of the spectroscopic data shows that the line of sight to the system crosses three kinematically sharp and well-separated interstellar reddening sources and that the stars rotate synchronously. We present new radial velocities that are in good agreement with earlier optical velocities and, thus, do not confirm the systematically smaller velocities obtained from IUE spectra. The optical velocities are analyzed simultaneously with the photometric data to derive accurate absolute dimesions for the binary components.The results show that TU Mus consists of an O7.5 primary with M_1=23.5 +/- 0.8 M_sun, R_1=7.48 +/- 0.08 R_sun and an O9.5 secondary with M_2=15.3 +/- 0.4 M_sun, R_2=6.15 +/- 0.07 R_sun in an overcontact configuration and that the orbital period has remained constant over the three decades covered by the observations. These results might imply that the mass transfer seen in late-type overcontact binaries does not occur in their early-type counterparts.