Avalanche photodiodes (APDs) are well-suited for single-photon detection on quantum communication satellites as they are a mature technology with high detection efficiency without requiring cryogenic cooling. They are, however, prone to significantly increased thermal noise caused by in-orbit radiation damage. Previous work demonstrated that a one-time application of thermal annealing reduces radiation-damage-induced APD thermal noise. Here we examine the effect of cyclical proton irradiation and thermal annealing, emulating the realistic operating profile of a satellite in low-Earth-orbit over a two-year life span. We show that repeated thermal annealing is effective in maintaining thermal noise of silicon APDs within a range suitable for quantum key distribution throughout the nominal mission life, and beyond. We examine two strategies---annealing at a fixed period of time, and annealing only when the thermal noise exceeds a pre-defined limit---and find that the latter exhibits lower thermal noise at end-of-life for most samples. We also observe that afterpulsing probability of the detector increases with cumulative proton irradiation. This knowledge helps guide design and tasking decisions for future space-borne quantum communication applications.

White dwarfs containing orbiting planetesimals or their debris represent crucial benchmarks by which theoretical investigations of post-main-sequence planetary systems may be calibrated. The photometric transit signatures of likely planetary debris in the ZTF J0139+5245 white dwarf system has an orbital period of about 110 days. An asteroid which breaks up to produce this debris may spin itself to destruction through repeated close encounters with the star without entering its Roche radius and without influence from the white dwarf's luminosity. Here, we place coupled constraints on the orbital pericentre ($q$) and the ratio ($\beta$) of the middle to longest semiaxes of a triaxial asteroid which disrupts outside of this white dwarf's Roche radius ($r_{\rm Roche}$) soon after attaining its 110-day orbit. We find that disruption within tens of years is likely when $\beta \lesssim 0.6$ and $q\approx 1.0-2.0r_{\rm Roche}$, and when $\beta \lesssim 0.2$ out to $q\approx 2.5r_{\rm Roche}$. Analysing the longer-timescale disruption of triaxial asteroids around ZTF J0139+5245 is desirable but may require either an analytical approach relying on ergodic theory or novel numerical techniques.

KIC 5773205 is the least luminous eclipsing M dwarf found in the Villanova catalog of eclipsing binaries detected by the {\it Kepler} mission. We processed and analyzed the three available quarters of mission data for this star and discovered a persistent periodic variation of the light curve with a period, which is in exact 4:5 commensurability to the orbital period. Three routes of interpretation are considered: 1) non-radial pulsations excited by the tidal interaction at a specific eigenfrequency; 2) a high-order spin-orbit resonance caused by the tides; 3) an ellipsoidal deformation caused by an outer orbiting companion in a mean motion resonance. All three explanations meet considerable difficulties, but the available facts seem to favor the tidally driven pulsation scenario. The star may represent a new type of heartbeat binary with tidally excited pulsations that are close to the orbital motion in frequency.

Observed planetary debris in white dwarf atmospheres predominately originate from the destruction of small bodies on highly eccentric ($>0.99$) orbits. Despite their importance, these minor planets have coupled physical and orbital evolution which has remained largely unexplored. Here, we present a novel approach for estimating the influence of fast chaotic rotation on the orbital evolution of high-eccentricity triaxial asteroids, and formally characterize the propagation of their angular rotation velocities and orbital elements as random time processes. By employing the impulse approximation, we demonstrate that the violent gravitational interactions during periastron passages transfer energy between the orbit and asteroid's rotation. If the distribution of spin impulses were symmetric around zero, then the net result would be a secular decrease of the semimajor axis and a further increase of the eccentricity. We find evidence, however, that the chaotic rotation may be self-regulated in such a manner that these effects are reduced or nullified. We discover that asteroids on highly eccentric orbits can break themselves apart --- in a type of YORP-less rotational fission --- without actually entering the Roche radius, with potentially significant consequences for the distribution of debris and energy requirements for gravitational scattering in metal-polluted white dwarf planetary systems. This mechanism provides a steady stream of material impacting a white dwarf without rapidly depleting the number of small bodies in the stellar system.

A new all-sky visible and Near-InfraRed (NIR) space astrometry mission with a wavelength cutoff in the K-band is not just focused on a single or small number of key science cases. Instead, it is extremely broad, answering key science questions in nearly every branch of astronomy while also providing a dense and accurate visible-NIR reference frame needed for future astronomy facilities. For almost 2 billion common stars the combination of Gaia and a new all-sky NIR astrometry mission would provide much improved proper motions, answering key science questions -- from the solar system and stellar systems, including exoplanet systems, to compact galaxies, quasars, neutron stars, binaries and dark matter substructures. The addition of NIR will result in up to 8 billion newly measured stars in some of the most obscured parts of our Galaxy, and crucially reveal the very heart of the Galactic bulge region. In this white paper we argue that rather than improving on the accuracy, a greater overall science return can be achieved by going deeper than Gaia and by expanding the wavelength range to the NIR.

Gaia is a revolutionary space mission developed by ESA and is delivering 5 parameter astrometry, photometry and radial velocities over the whole sky with astrometric accuracies down to a few tens of micro-arcseconds. A weakness of Gaia is that it only operates at optical wavelengths. However, much of the Galactic centre and the spiral arm regions, important for certain studies, are obscured by interstellar extinction and this makes it difficult for Gaia to deeply probe. This problem can be overcome by switching to the Near Infra-Red (NIR) but this is not possible with silicon CCDs. Additionally, to scan the entire sky and make global absolute parallax measurements the spacecraft must have a constant rotation and this requires the detectors operate in Time Delayed Integration (TDI) mode or similar.

Gaia is currently revolutionizing modern astronomy. However, much of the Galactic plane, center and the spiral arm regions are obscured by interstellar extinction, rendering them inaccessible because Gaia is an optical instrument. An all-sky near infrared (NIR) space observatory operating in the optical NIR, separated in time from the original Gaia would provide microarcsecond NIR astrometry and millimag photometry to penetrate obscured regions unraveling the internal dynamics of the Galaxy.

Observational evidence of white dwarf planetary systems is dominated by the remains of exo-asteroids through accreted metals, debris discs, and orbiting planetesimals. However, exo-planets in these systems play crucial roles as perturbing agents, and can themselves be perturbed close to the white dwarf Roche radius. Here, we illustrate a procedure for computing the tidal interaction between a white dwarf and a near-spherical solid planet. This method determines the planet's inward and/or outward drift, and whether the planet will reach the Roche radius and be destroyed. We avoid constant tidal lag formulations and instead employ the self-consistent secular Darwin-Kaula expansions from Bou\'{e} & Efroimsky (2019), which feature an arbitrary frequency dependence on the quality functions. We adopt wide ranges of dynamic viscosities and spin rates for the planet in order to straddle many possible outcomes, and provide a foundation for the future study of individual systems with known or assumed rheologies. We find that: (i) massive Super-Earths are destroyed more readily than minor planets (such as the ones orbiting WD 1145+017 and SDSS J1228+1040), (ii) low-viscosity planets are destroyed more easily than high-viscosity planets, and (iii) the boundary between survival and destruction is likely to be fractal and chaotic.

Motivated by the existence of binary systems where a stellar-mass black hole is bound to a normal star, we selected four red giants with large radial velocity (RV) variation from the survey of SIM grid stars and monitored their RVs for several months. None turned out to contain a massive companion above 2.5 solar masses. The red giant TYC 9299-1080-1 with a large RV and a large proper motion is a single-lined spectroscopic binary with a period of 81 days. It is an extreme halo object moving at 350 km s$^{-1}$ almost directly toward the Galactic center. HD 206092 is a double-lined binary with a short period of 4.37 days. It belongs to the rare class of active RS CVn-type binaries with evolved primary components, apparently undergoing mass transfer. The X-ray luminosity of HD 206092 is about twice as high as the most luminous coronal X-ray emitters observed by ROSAT, including II Peg and the prototype star RS CVn. HD~318347 has a variable double-peaked emission-line spectrum (not a giant), while HD~324668 has a constant RV. Despite the overall good quality of the SIM survey data confirmed by a comparison with Gaia DR2 mean radial velocities, the few large RV variations are explained, mostly, by erroneous data. We discuss the significance of the non-detection of massive companions in the SIM grid sample and the associated work.

We investigate a sample of 3412 {\it International Celestial Reference Frame} (ICRF3) extragalactic radio-loud sources with accurate positions determined by VLBI in the S/X band, mostly active galactic nuclei (AGN) and quasars, which are cross-matched with optical sources in the second Gaia data release (Gaia DR2). The main goal of this study is to determine a core sample of astrometric objects that define the mutual orientation of the two fundamental reference frames, the Gaia (optical) and the ICRF3 (radio) frames. The distribution of normalized offsets between the VLBI sources and their optical counterparts is non-Rayleigh, with a deficit around the modal value and a tail extending beyond the 3-sigma confidence level. A few filters are applied to the sample in order to discard double cross-matches, confusion sources, and Gaia astrometric solutions of doubtful quality. {\it Panoramic Survey Telescope and Rapid Response System} (Pan-STARRS) and {\it Dark Energy Survey} (DES) stacked multi-color images are used to further deselect objects that are less suitable for precision astrometry, such as extended galaxies, double and multiple sources, and obvious misidentifications. After this cleaning, 2642 quasars remain, of which 20\% still have normalized offset magnitudes exceeding 3, or 99\% confidence level. We publish a list of 2118 radio-loud quasars of prime astrometric quality. The observed dependence of binned median offset on redshift shows the expected decline at small redshifts, but also an unexpected rise at $z\sim 1.6$, which may be attributed to the emergence of the C IV emission line in the Gaia's $G$ band. The Gaia DR2 parallax zero-point is found to be color-dependent, suggesting an uncorrected instrumental calibration effect.

We perform numerical simulations of the TRAPPIST-1 system of seven exoplanets orbiting a nearby M dwarf, starting with a previously suggested stable configuration. The long-term stability of this configuration is confirmed, but the motion of planets is found to be chaotic. The eccentricity values are found to vary within finite ranges. The rates of tidal dissipation and tidal evolution of orbits are estimated, assuming an Earth-like rheology for the planets. We find that under this assumption the planets b, d, e were captured in the 3:2 or higher spin-orbit resonances during the initial spin-down but slipped further down into the 1:1 resonance. Dependent on its rheology, the innermost planet b may be captured in a stable pseudosynchronous rotation. Non-synchronous rotation ensures higher levels of tidal dissipation and internal heating. The positive feedback between the viscosity and the dissipation rate -- and the ensuing runaway heating -- are terminated by a few self-regulation processes. When the temperature is high and the viscosity is low enough, the planet spontaneously leaves the 3:2 resonance. Further heating is stopped either by passing the peak dissipation or by the emergence of partial melt in the mantle. In the post-solidus state, the tidal dissipation is limited to the levels supported by the heat transfer efficiency. The tides on the host star are unlikely to have had a significant dynamical impact. The tides on the synchronized inner planets tend to reduce these planets' orbital eccentricity, possibly contributing thereby to the system's stability.

The VLBI USNO 2016A (U16A) solution is part of a work-in-progress effort by USNO towards the preparation of the ICRF3. Most of the astrometric improvement with respect to the ICRF2 is due to the re-observation of the VCS sources. Our objective in this paper is to assess U16A's astrometry. A comparison with ICRF2 shows statistically significant offsets of size 0.1 mas between the two solutions. While Gaia DR1 positions are not precise enough to resolve these offsets, they are found to be significantly closer to U16A than ICRF2. In particular, the trend for typically larger errors for Southern sources in VLBI solutions are decreased in U16A. Overall, the VLBI-Gaia offsets are reduced by 21%. The U16A list includes 718 sources not previously included in ICRF2. Twenty of those new sources have statistically significant radio-optical offsets. In two-thirds of the cases, these offsets can be explained from PanSTARRS images.

Highly volcanic exoplanets, which can be variously characterized as 'lava worlds', 'magma ocean worlds', or 'super-Ios' are high priority targets for investigation. The term 'lava world' may refer to any planet with extensive surface lava lakes, while the term 'magma ocean world' refers to planets with global or hemispherical magma oceans at their surface. 'Highly volcanic planets', including super-Ios, may simply have large, or large numbers of, active explosive or extrusive volcanoes of any form. They are plausibly highly diverse, with magmatic processes across a wide range of compositions, temperatures, activity rates, volcanic eruption styles, and background gravitational force magnitudes. Worlds in all these classes are likely to be the most characterizable rocky exoplanets in the near future due to observational advantages that stem from their preferential occurrence in short orbital periods and their bright day-side flux in the infrared. Transit techniques should enable a level of characterization of these worlds analogous to hot Jupiters. Understanding processes on highly volcanic worlds is critical to interpret imminent observations. The physical states of these worlds are likely to inform not just geodynamic processes, but also planet formation, and phenomena crucial to habitability. Volcanic and magmatic activity uniquely allows chemical investigation of otherwise spectroscopically inaccessible interior compositions. These worlds will be vital to assess the degree to which planetary interior element abundances compare to their stellar hosts, and may also offer pathways to study both the very young Earth, and the very early form of many silicate planets where magma oceans and surface lava lakes are expected to be more prevalent. We suggest that highly volcanic worlds may become second only to habitable worlds in terms of both scientific and public long-term interest.

The geophysics of extrasolar planets is a scientific topic often regarded as standing largely beyond the reach of near-term observations. This reality in no way diminishes the central role of geophysical phenomena in shaping planetary outcomes, from formation, to thermal and chemical evolution, to numerous issues of surface and near-surface habitability. We emphasize that for a balanced understanding of extrasolar planets, it is important to look beyond the natural biases of current observing tools, and actively seek unique pathways to understand exoplanet interiors as best as possible during the long interim prior to a time when internal components are more directly accessible. Such pathways include but are not limited to: (a) enhanced theoretical and numerical modeling, (b) laboratory research on critical material properties, (c) measurement of geophysical properties by indirect inference from imprints left on atmospheric and orbital properties, and (d) the purpose-driven use of Solar System object exploration expressly for its value in comparative planetology toward exoplanet-analogs. Breaking down barriers that envision local Solar System exploration, including the study of Earth's own deep interior, as separate from and in financial competition with extrasolar planet research, may greatly improve the rate of needed scientific progress for exoplanet geophysics. As the number of known rocky and icy exoplanets grows in the years ahead, we expect demand for expertise in 'exogeoscience' will expand at a commensurately intense pace. We highlight key topics, including: how water oceans below ice shells may dominate the total habitability of our galaxy by volume, how free-floating nomad planets may often attain habitable subsurface oceans supported by radionuclide decay, and how deep interiors may critically interact with atmospheric mass loss via dynamo-driven magnetic fields.

Single-photon detectors in space must retain useful performance characteristics despite being bombarded with sub-atomic particles. Mitigating the effects of this space radiation is vital to enabling new space applications which require high-fidelity single-photon detection. To this end, we conducted proton radiation tests of various models of avalanche photodiodes (APDs) and one model of photomultiplier tube potentially suitable for satellite-based quantum communications. The samples were irradiated with 106 MeV protons at doses approximately equivalent to lifetimes of 0.6 , 6, 12 and 24 months in a low-Earth polar orbit. Although most detection properties were preserved, including efficiency, timing jitter and afterpulsing probability, all APD samples demonstrated significant increases in dark count rate (DCR) due to radiation-induced damage, many orders of magnitude higher than the 200 counts per second (cps) required for ground-to-satellite quantum communications. We then successfully demonstrated the mitigation of this DCR degradation through the use of deep cooling, to as low as -86 degrees C. This achieved DCR below the required 200 cps over the 24 months orbit duration. DCR was further reduced by thermal annealing at temperatures of +50 to +100 degrees C.

Accurate parallaxes from Gaia DR1 (TGAS) are combined with GALEX visual Nuv magnitudes to produce absolute Mnuv magnitudes and an ultraviolet HR diagram for a large sample of astrometric stars. A functional fit is derived of the lower envelope main sequence of the nearest 1403 stars (distance < 40 pc), which should be reddening-free. Using this empirical fit, 50 nearby stars are selected with significant Nuv excess. These are predominantly late K and early M dwarfs, often associated with X-ray sources, and showing other manifestations of magnetic activity. The sample may include systems with hidden white dwarfs, stars younger than the Pleiades, or, most likely, tight interacting binaries of the BY Dra-type. A separate collection of 40 stars with precise trigonometric parallaxes and Nuv - G colors bluer than 2 mag is presented. It includes several known novae, white dwarfs, and binaries with hot subdwarf (sdOB) components, but most remain unexplored.

KIC 7341653 is one of several late-type M dwarfs observed in the main mission of Kepler with peculiar infrared colors placing them in the domain of suspected young stellar objects (YSO). It is likely associated with a powerful X-ray emitter with X-ray flares. Kepler light curves reveal two distinct types of activity: frequent flares lasting from less than 30 min to a few hours and a periodic variability with a period of 0.5463441(7) d. The largest flare detected increased the flux in the Kepler passband by a factor of 2.8 and released an estimated 4$\times10^{34}$ erg of energy in the Kepler band. Segmented periodogram analysis reveals that the amplitude of the periodic variation was subject to secular changes, dropping from peak values around 20 ppt to below 5 ppt toward the end of the mission, while the phase varied periodically with an amplitude of 0.15 radians and period 362(3) d. Two possible interpretations of the phase periodicity are discussed: a migrating long-lived photospheric spot and a Doppler frequency shift generated by a solar-mass faint companion, such as a white dwarf.

In the context of the ESA M5 (medium mission) call we proposed a new satellite mission, Theia, based on relative astrometry and extreme precision to study the motion of very faint objects in the Universe. Theia is primarily designed to study the local dark matter properties, the existence of Earth-like exoplanets in our nearest star systems and the physics of compact objects. Furthermore, about 15 $\%$ of the mission time was dedicated to an open observatory for the wider community to propose complementary science cases. With its unique metrology system and "point and stare" strategy, Theia's precision would have reached the sub micro-arcsecond level. This is about 1000 times better than ESA/Gaia's accuracy for the brightest objects and represents a factor 10-30 improvement for the faintest stars (depending on the exact observational program). In the version submitted to ESA, we proposed an optical (350-1000nm) on-axis TMA telescope. Due to ESA Technology readiness level, the camera's focal plane would have been made of CCD detectors but we anticipated an upgrade with CMOS detectors. Photometric measurements would have been performed during slew time and stabilisation phases needed for reaching the required astrometric precision.

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.

A significant number of double stars with separations up to 2.5 arcsec are present in the Gaia Data Release 1 astrometric catalogs. Limiting our analysis to a well-studied sample of 1124 doubles resolved by Hipparcos, provided with individual Tycho component photometry, and cross-matched with the TGAS catalog, we estimate a rate of at least 3% for brighter double stars in Gaia DR1, which should be resolved in the future data releases. Gaia astrometric results are affected by unresolved duplicity. The variance-normalized quadratic differences of proper motion between Gaia and Hipparcos do not follow the expected chi^2 distribution and show signs of powerful degradation in the components aligned with the axes of the double systems. This concerns only pairs with separation below 1.2 -- 1.5 arcsec, which mostly remain unresolved in Gaia DR1. On the other hand, the orthogonal proper motion components and parallaxes do not have any detectable perturbation, as well as all astrometry for separations above 1.5 arcsec. Gaia parallaxes do not seem to be perturbed by duplicity, with Gaia - Hipparcos differences being systematically smaller than the expectation. The rate of incorrectly identified, or swapped, companions is estimated at 0.4%.

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.

We investigate a sample of 2293 ICRF2 extragalactic radio-loud sources with accurate positions determined by VLBI, mostly active galactic nuclei (AGN) and quasars, which are cross-matched with optical sources in the first Gaia release (Gaia DR1). The distribution of offsets between the VLBI sources and their optical counterparts is strongly non-Gaussian, with powerful wings extending beyond 1 arcsecond. Limiting our analysis to only high-confidence difference detections, we find (and publish) a list of 188 objects with normalized variances above 12 and offsets below 1 arcsecond. Pan-STARRS stacked and monochromatic images resolve some of these sources indicating the presence of double sources, confusion sources, or pronounced extended structures. Some 89 high-quality objects, however, do not show any perturbations and appear to be star-like single sources, yet displaced by multiples of the expected error from the radio-loud AGN. We conclude that a fraction of luminous AGN (more than 4%) can be physically dislodged from the optical centers of their parent galaxies.

We apply a PCA-based pre-whitening method to the entire collection of main Kepler mission long-cadence data for KIC 8462852 spanning four years. This technique removes the correlated variations of instrumental origin in both the detected light curves and astrometry, resolving intrinsic changes in flux and image position of less than 100 ppm and 1 mas, respectively. Beside the major dips in the light curve during mission quarters 8 and 16, when the flux dropped by up to 20%, we confirm multiple smaller dips across the time span of observation with amplitudes ranging from 0.1% to 7%. A variation of flux with a period of 0.88 d and a half-amplitude of approximately 90 ppm is confirmed in the PCA-cleaned data. We find that the phase of the wave is steady over the entire 15-month interval. We confidently detect a weak variability-induced motion (VIM) effect in the cleaned astrometric trajectories, when the moment-based centroids shift synchronously with the flux dips by up to 0.0008 pixels on the detector. The inconsistent magnitude and direction of VIM effects within the same quarter point at more than one source of photometric variability in the blended image. The 0.88 d periodicity comes from a different source, not from the target star KIC 8462852. We discuss a possible interpretation of the bizarre properties of the source as a swarm of interstellar junk (comets and planetoids) crossing the line of sight to the star and its optical companions at approximately 7 mas per year.

We describe development and application of a Global Astrometric Solution (GAS) to the problem of Pan-STARRS1 (PS1) astrometry. Current PS1 astrometry is based on differential astrometric measurements using 2MASS reference stars, thus PS1 astrometry inherits the errors of the 2MASS catalog. The GAS, based on a single, least squares adjustment to approximately 750k grid stars using over 3000 extragalactic objects as reference objects, avoids this catalog-to-catalog propagation of errors to a great extent. The GAS uses a relatively small number of Quasi-Stellar Objects (QSOs, or distant AGN) with very accurate (<1 mas) radio positions, referenced to the ICRF2. These QSOs provide a hard constraint in the global least squares adjustment. Solving such a system provides absolute astrometry for all the stars simultaneously. The concept is much cleaner than conventional astrometry but is not easy to perform for large catalogs. In this paper we describe our method and its application to Pan-STARRS1 data. We show that large-scale systematic errors are easily corrected but our solution residuals for position (~60 mas) are still larger than expected based on simulations (~10 mas). We provide a likely explanation for the reason the small-scale residual errors are not corrected in our solution as would be expected.

Variability induced motion (VIM) is an observable effect in simultaneous astrometric and photometric measurements caused by brightness variation in one of the components of a double source or blended image, which manifests itself as a strongly correlated shift of the optical photocenter. We have processed the entire collection of the Kepler long-cadence light curve data looking for correlated signals in astrometry and photometry on the time basis of a quarter year. Limiting the VIM correlation coefficient to 0.3, VIM events are detected for 129,525 Kepler stars at least in one quarter. Of 7305 Kepler objects of interest (KOI), 4440 are detected as VIM at least once. Known variable stars and resolved double stars have elevated rates of VIM detection. Confident VIM occurrences are found for stars with suggested superflare events, indicating possible signal contamination. We present a complete catalog of all quarterly VIM detections. This catalog should be checked for such astrophysically significant events as transits of exoplanets, new eclipsing stars, and superflares of solar-type stars.

A formula for the tidal dissipation rate in a spherical body is derived from first principles, to correct some mathematical inaccuracies found in the literature. The development is combined with the Darwin-Kaula formalism for tides. Our intermediate results are compared with those by Zschau (1978) and Platzman (1984). When restricted to the special case of an incompressible spherical planet spinning synchronously without libration, our final formula can be compared with the commonly used expression from Peale & Cassen (1978, Eqn. 31). The two turn out to differ. In our expression, the contributions from all Fourier modes are positive-definite, this not being the case of the formula from Ibid. (The presence of negative terms in their formula was noticed by Makarov 2013.) Examples of application of our expression for the tidal damping rate are provided in the work by Makarov and Efroimsky (2014).

Tidal dissipation of kinetic energy, when it is strong enough, tends to synchronize the rotation of planets and moons with the mean orbital motion, or drive it into long-term stable spin-orbit resonances. As the orbital motion undergoes periodic acceleration due to a finite orbital eccentricity, the spin rate oscillates around the equilibrium mean value too, giving rise to the forced, or eccentricity-driven, librations. Both the shape and amplitude of forced librations of synchronous viscoelastic planets and moons are defined by a combination of two different types of perturbative torque, the tidal torque and the triaxial torque. Consequently, forced librations can be tidally dominated (e.g., Io and possibly Titan) or deformation-dominated (e.g., the Moon) depending on a set of orbital, rheological, and other physical parameters. With small eccentricities, for the former kind, the largest term in the libration angle can be minus cosine of the mean anomaly, whereas for the latter kind, it is minus sine of the mean anomaly. The shape and the amplitude of tidal forced librations determine the rate of orbital evolution of synchronous planets and moons, i.e., the rate of dissipative damping of semimajor axis and eccentricity. The known super-Earth exoplanets can exhibit both kinds of libration, or a mixture thereof, depending on, for example, the effective Maxwell time of their rigid mantles. Our approach can be extended to estimate the amplitudes of other libration harmonics, as well as the forced libration in non-synchronous spin-orbit resonances.

We present an all-sky sample of ~ 1.4 million AGNs meeting a two color infrared photometric selection criteria for AGNs as applied to sources from the Wide-Field Infrared Survey Explorer final catalog release (AllWISE). We assess the spatial distribution and optical properties of our sample and find that the results are consistent with expectations for AGNs. These sources have a mean density of ~ 38 AGNs per square degree on the sky, and their apparent magnitude distribution peaks at g ~ 20, extending to objects as faint as g ~ 26. We test the AGN selection criteria against a large sample of optically-identified stars and determine the "leakage" (that is, the probability that a star detected in an optical survey will be misidentified as a QSO in our sample) rate to be < 4.0 x 10^-5. We conclude that our sample contains almost no optically-identified stars (< 0.041%), making this sample highly promising for future celestial reference frame work by significantly increasing the number of all-sky, compact extragalactic objects. We further compare our sample to catalogs of known AGNs/QSOs and find a completeness value of > 84% (that is, the probability of correctly identifying a known AGN/QSO is at least 84%) for AGNs brighter than a limiting magnitude of R < 19. Our sample includes approximately 1.1 million previously uncatalogued AGNs.

Proper motions are computed and collected in a catalog using the Hipparcos positions (epoch 1991.25) and URAT1 positions (epoch 2012.3 to 2014.6). The goal is to obtain a significant improvement on the proper motion accuracy of single stars in the northern hemisphere, and to identify new astrometric binaries perturbed by orbital motion. For binaries and multiple systems, the longer baseline of Tycho2 (~ 100 yr) makes it more reliable despite its larger formal uncertainties. The resulting proper motions obtained for 67,340 stars have a consequent gain in accuracy by a factor of ~ 3 compared to Hipparcos. Comparison between UrHip and Hipparcos shows that they are reasonably close, but also reveals stars with large discrepant proper motions, a fraction of which are potential binary candidates.

A wide range of exoplanet and exomoon models are characterized by a finite average rigidity and a viscosity much lower than the typical values for terrestrials. Such semiliquid bodies may or may not have rigid crusts with permanent figures. Unlike planets with solid mantles and Earth-like rheology, semiliquid bodies can be captured into stable pseudosynchronous spin resonance, where the average rate of rotation is higher than the synchronous 1:1 resonance. Two basic conditions are derived for capture of planets with a triaxial figure into pseudosynchronous rotation, one related to the characteristic tidal wave number (the product of the tidal frequency by the Maxwell time), and the other to the orbital eccentricity. If a semiliquid object does not satisfy either of the two conditions, it is captured into the synchronous resonance. For nearly axially symmetric bodies, only the first condition is in place, and the other is much relaxed, so they should predominantly be pseudosynchronous. It is also pointed out that the equilibrium pseudosychronous rotation rate can not reach the widely used asymptotic value from the constant time lag model but is in reality closer to the synchronous spin.

[Abridged] The Study Analysis Group 8 of the NASA Exoplanet Analysis Group was convened to assess the current capabilities and the future potential of the precise radial velocity (PRV) method to advance the NASA goal to "search for planetary bodies and Earth-like planets in orbit around other stars.: (U.S. National Space Policy, June 28, 2010). PRVs complement other exoplanet detection methods, for example offering a direct path to obtaining the bulk density and thus the structure and composition of transiting exoplanets. Our analysis builds upon previous community input, including the ExoPlanet Community Report chapter on radial velocities in 2008, the 2010 Decadal Survey of Astronomy, the Penn State Precise Radial Velocities Workshop response to the Decadal Survey in 2010, and the NSF Portfolio Review in 2012. The radial-velocity detection of exoplanets is strongly endorsed by both the Astro 2010 Decadal Survey "New Worlds, New Horizons" and the NSF Portfolio Review, and the community has recommended robust investment in PRVs. The demands on telescope time for the above mission support, especially for systems of small planets, will exceed the number of nights available using instruments now in operation by a factor of at least several for TESS alone. Pushing down towards true Earth twins will require more photons (i.e. larger telescopes), more stable spectrographs than are currently available, better calibration, and better correction for stellar jitter. We outline four hypothetical situations for PRV work necessary to meet NASA mission exoplanet science objectives.

We present analysis of precision radial velocities (RV) of 1134 mostly red giant stars in the southern sky, selected as candidate astrometric grid objects for the Space Interferometry Mission (SIM). Only a few (typically, 2 or 3) spectroscopic observations per star have been collected, with the main goal of screening binary systems. The estimated rate of spectroscopic binarity in this sample of red giants is 32% at the 0.95 confidence level, and 46% at the 0.75 confidence. The true binarity rate is likely to be higher, because our method is not quite sensitive to very wide binaries and low-mass companions. The estimated lower and upper bounds of stellar RV jitter for the entire sample are 24 and 51 m/s, respectively; the adopted mean value is 37 m/s. A few objects of interest are identified with large variations of radial velocities, implying abnormally high mass ratios.

In Efroimsky & Makarov (2014), we derived from the first principles a formula for the tidal heating rate in a tidally perturbed homogeneous sphere. We compared it with the formulae used in the literature, and pointed out the differences. Using this result, we now present three case studies - Mercury, Kepler-10b, and a triaxial Io. A very sharp frequency-dependence of k2/Q near spin-orbit resonances yields a similarly sharp dependence of k2/Q on the spin rate. This indicates that physical libration may play a major role in tidal heating of synchronously rotating bodies. The magnitude of libration in the spin rate being defined by the planet's triaxiality, the latter should be a factor determining the dissipation rate. Other parameters equal, a synchronously rotating body with a stronger triaxiality should generate more heat than a similar body of a more symmetrical shape. Further in the paper, we discuss scenarios where initially triaxial objects melt and lose their triaxiality. Thereafter, dissipation in them becomes less intensive; so the bodies freeze. The tidal bulge becomes a new permanent figure, with a new triaxiality lower than the original. In the paper, we also derive simplified, approximate expressions for dissipation rate in a rocky planet of the Maxwell rheology, with a not too small Maxwell time. The three expressions derived pertain to the cases of a synchronous spin, a 3:2 resonance, and a nonresonant rotation; so they can be applied to most close-in super-Earth exoplanets detected thus far. In such bodies, the rate of tidal heating outside of synchronous rotation is weakly dependent on the eccentricity and obliquity, provided both these parameters are small or moderate. According to our calculation, Kepler-10b could hardly survive the great amount of tidal heating without being synchronised, circularised and also reshaped through a complete or partial melt-down.

While it is accepted that the eccentricity of Mercury (0.206) favours entrapment into the 3:2 spin-orbit resonance, open is the question how and when the capture took place. A recent work by Makarov (2012) has demonstrated that trapping into this resonance is certain if the eccentricity is larger than 0.2, provided that we use a realistic tidal model, the one which is based on the Darwin-Kaula expansion of the tidal torque. The physics-based tidal model changes dramatically the statistics of the possible final spin states. First, we discover that after only one encounter with the spin-orbit 3:2 resonance this resonance becomes the most probable end-state. Second, if a capture into this (or any other) resonance takes place, the capture becomes final, several crossings of the same state being forbidden by our model. Third, within our model the trapping of Mercury happens much faster than previously believed: for most histories, 10 - 20 Myr are sufficient. Fourth, even a weak laminar friction between the solid mantle and a molten core would most likely result in a capture in the 2:1 or even higher resonance. So the principal novelty of our paper is that the 3:2 end-state is more ancient than the same end-state obtained when the constant time lag model is employed. The swift capture justifies our treatment of Mercury as a homogeneous, unstratified body whose liquid core had not yet formed by the time of trapping. We also provide a critical analysis of the hypothesis by Wieczorek et al. (2012) that the early Mercury might had been retrograde, whereafter it synchronised its spin and then accelerated it to the 3:2 resonance. Accurate processing of the available data on cratering does not support that hypothesis, while the employment of a realistic rheology invalidates a key element of the hypothesis, an intermediate pseudosynchronous state needed to spin-up to the 3:2 resonance.

We investigate the dynamical evolution of the potentially habitable super-earth GJ 667Cc in the multiple system of at least two exoplanets orbiting a nearby M dwarf, paying special attention to its spin-orbital state. The published radial velocities for this star are re-analyzed and evidence is found for additional periodic signals, which could be taken for two additional planets on eccentric orbits. Such systems are not dynamically viable and break up quickly in numerical integrations. Limiting the scope to the two originally detected planets, we assess the dynamical stability of the system and find no evidence for bounded chaos in the orbital motion. The orbital eccentricity of the planets b and c is found to change cyclicly in the range 0.06 - 0.28 and 0.05 - 0.25, respectively, with a period of approximately 0.46 yr. Taking the eccentricity variation into account, numerical integrations are performed of the differential equations modeling the spin-orbit interaction of the planet GJ 667Cc with its host star assuming a terrestrial composition of its mantle. Depending on the interior temperature of the planet, it is likely to be entrapped in the 3:2 (probability 0.51) or even higher spin-orbit resonance. It is less likely to reach the 1:1 resonance (probability 0.24). The estimated characteristic spin-down times are quite short for the two planets, i.e., within 1 Myr for planet c and even shorter for planet b. Both planet arrived at their current and, most likely, ultimate spin-orbit states a long time ago. The rate of tidal dissipation of energy is formidably high in the planets of GJ 667, estimated at $10^{23.7}$ and $10^{26.7}$ J yr$^{-1}$ for c and b, respectively. This raises a question of how such relatively massive, close super-Earths could survive overheating and destruction.

If the Moon's spin evolved from faster prograde rates, it could have been captured into a higher spin-orbit resonance than the current 1:1 resonance. At the current value of orbital eccentricity, the probability of capture into the 3:2 resonance is as high as 0.6, but it strongly depends on the temperature and average viscosity of the Moon's interior. A warmer, less viscous Moon on a higher-eccentricity orbit is even more easily captured into supersynchronous resonances. We discuss two likely scenarios for the present spin-orbit state: a cold Moon on a low-eccentricity orbit and a retrograde initial rotation.

GJ 581d is a potentially habitable super-Earth in the multiple system of exoplanets orbiting a nearby M dwarf. We investigate this planet's long-term dynamics, with an emphasis on its probable final rotation states acquired via tidal interaction with the host. The published radial velocities for the star are re-analysed with a benchmark planet detection algorithm, to confirm that there is no evidence for the recently proposed two additional planets (f and g). Limiting the scope to the four originally detected planets, we assess the dynamical stability of the system and find bounded chaos in the orbital motion. For the planet d, the characteristic Lyapunov time is 38 yr. Long-term numerical integration reveals that the system of four planets is stable, with the eccentricity of the planet d changing quasi-periodically in a tight range around 0.27, and with its semimajor axis varying only a little. The spin-orbit interaction of GJ 581d with its host star is dominated by the tides exerted by the star on the planet. We model this interaction, assuming a terrestrial composition of the mantle. Besides the customarily included secular parts of the triaxiality-caused and tidal torques, we also include these torques' oscillating components. It turns out that, dependent on the mantle temperature, the planet gets trapped into the 2:1 or an even higher spin-orbit resonance. It is very improbable that the planet could have reached the 1:1 resonance. This enhances the possibility of the planet being suitable for sustained life.

We reexamine the popular belief that a telluric planet or satellite on an eccentric orbit can, outside a spin-orbit resonance, be captured in a quasi-static tidal equilibrium called pseudosynchronous rotation. The existence of such configurations was deduced from oversimplified tidal models assuming either a constant tidal torque or a torque linear in the tidal frequency. A more accurate treatment requires that the torque be decomposed into the Darwin-Kaula series over the tidal modes, and that this decomposition be combined with a realistic choice of rheological properties of the mantle. This development demonstrates that there exist no stable equilibrium states for solid planets and moons, other than spin-orbit resonances.

Tidal torques play a key role in rotational dynamics of celestial bodies. They govern these bodies' tidal despinning, and also participate in the subtle process of entrapment of these bodies into spin-orbit resonances. This makes tidal torques directly relevant to the studies of habitability of planets and their moons. Our work begins with an explanation of how friction and lagging should be built into the theory of bodily tides. Although much of this material can be found in various publications, a short but self-consistent summary on the topic has been lacking in the hitherto literature, and we are filling the gap. After these preparations, we address a popular concise formula for the tidal torque, which is often used in the literature, for planets or stars.We explain why the derivation of this expression, offered in the paper by Goldreich (1966; AJ 71, 1 - 7) and in the books by Kaula (1968, eqn. 4.5.29), and Murray & Dermott (1999, eqn. 4.159), implicitly sets the time lag to be frequency independent. Accordingly, the ensuing expression for the torque can be applied only to bodies having a very special (and very hypothetical) rheology which makes the time lag frequency independent, i.e, the same for all Fourier modes in the spectrum of tide. This expression for the torque should not be used for bodies of other rheologies. Specifically, the expression cannot be combined with an extra assertion of the geometric lag (or the phase lag) being constant, because at finite eccentricities the said assumption is incompatible with the constant-time-lag condition.

We propose a computer-efficient and accurate method of estimation of spatially correlated errors in astrometric positions, parallaxes and proper motions obtained by space and ground-based astrometry missions. In our method, the simulated observational equations are set up and solved for the coefficients of scalar and vector spherical harmonics representing the output errors, rather than for individual objects in the output catalog. Both accidental and systematic correlated errors of astrometric parameters can be accurately estimated. The method is demonstrated on the example of the JMAPS mission, but can be used for other projects of space astrometry, such as SIM or JASMINE.

A subset of 51 Hipparcos astrometric binaries among FG dwarfs within 67pc has been surveyed with the NICI adaptive optics system at Gemini-S, directly resolving for the first time 17 sub-arcsecond companions and 7 wider ones. Using these data together with published speckle interferometry of 57 stars, we compare the statistics of resolved astrometric companions with those of a simulated binary population. The fraction of resolved companions is slightly lower than expected from binary statistics. About 10% of astrometric companions could be "dark" (white dwarfs and close pairs of late M-dwarfs). To our surprise, several binaries are found with companions too wide to explain the acceleration. Re-analysis of selected intermediate astrometric data shows that some acceleration solutions in the original Hipparcos catalog are spurious.

The dynamical evolution of terrestrial planets resembling Mercury in the vicinity of spin-orbit resonances is investigated using comprehensive harmonic expansions of the tidal torque taking into account the frequency-dependent quality factors and Love numbers. The torque equations are integrated numerically with a small step in time, includng the oscillating triaxial torque components but neglecting the layered structure of the planet and assuming a zero obliquity. We find that a Mercury-like planet with its current value of orbital eccentricity (0.2056) is always captured in the 3:2 resonance. The probability of capture in the higher 2:1 resonance is approximately 0.23. These results are confirmed by a semi-analytical estimation of capture probabilities as functions of eccentricity for both prograde and retrograde evolution of spin rate. As follows from analysis of equilibrium torques, entrapment in the 3:2 resonance is inevitable at eccentricities between 0.2 and 0.41. Considering the phase space parameters at the times of periastron, the range of spin rates and phase angles, for which an immediate resonance passage is triggered, is very narrow, and yet, a planet like Mercury rarely fails to align itself into this state of unstable equilibrium before it traverses the 2:1 resonance.

Quasometry is precision measurement of celestial positions and apparent motion of very distant extragalactic objects, such as quasars, galactic nuclei, and QSOs. We use this term to identify a specific area of research, the methodology of which differs from that of general astrometry. The main purpose of quasometry is to link the sub-milliarcsecond radio frame (ICRF) with the existing and emerging optical reference frames of similar accuracy, constructed by astrometric satellites. Some of the main difficulties in achieving this goal are discussed, e.g., the extended structures of quasar hosts, apparent motion on the sky, optical variability, galactic companions, faintness. Besides the strategic purpose, quasometry is undoubtedly useful for global astrometric surveys, as it helps to verify or even correct the resulting reference frames. There are two options of using measurements of distant quasars in a global astrometric solution: 1) hard constraints embedded in the fabric of observational equations; 2) {\it a posteriori} fitting of zonal errors. The relative benefits and shortcoming of the two options are reviewed. A relatively small set of about 200 carefully selected reference quasars can go a long way in improving the astrometric value of a space mission, if they are sufficiently bright, stable, fairly uniformly distributed on the sky, and are defining sources in the ICRF. We present an ongoing program at the USNO to construct a quality set of optical quasars with the required properties and to enhance the ICRF with new sources in the areas where known, well-observed quasars are scarce.

The relative motion of stars and other celestial objects in very wide pairs, separated by distances of the order of 1 pc, is strongly influenced by the tidal gravitational potential of the Galaxy. The Coriolis component of the horizontal tidal force in the rotating reference frame tends to disrupt such marginally bound pairs. However, even extremely wide pairs of bodies can be bound over intervals of time comparable to the Hubble time, under appropriate initial conditions. Here we show that for arbitrary chosen initial coordinates of a pair of stars, there exists a volume of the space of initial velocity components where the orbits remain bound in the planar tidal field for longer than 10 Gyr, even though the initial separation is well outside the Jacobi radius. The boundary of this phase space of stable orbits is fractal, and the motion at the boundary conditions is clearly chaotic. We found that the pairs may remain confined for several Gyr, and then suddenly disintegrate due to a particularly close rendezvous. By reversing such long-term stable orbits, we find that entrapment of unrelated stars into wide pairs is possible, but should be quite rare. Careful analysis of precision astrometry surveys revealed that extremely wide pairs of stars are present in significant numbers in the Galaxy. These results are expected to help discriminating the cases of genuine binarity and chance entrapment, and to make inroads in testing the limits of Newtonian gravitation.

The non-orthogonality of algebraic polynomials of field coordinates traditionally used to model field-dependent corrections to astrometric measurements, gives rise to subtle adverse effects. In particular, certain field dependent perturbations in the observational data propagate into the adjusted coefficients with considerable magnification. We explain how the worst perturbation, resulting in the largest solution error, can be computed for a given non-orthogonal distortion model. An algebraic distortion model of full rank can be converted into a fully orthonormal model based on the Zernike polynomials for a circular field of view, or a basis of functions constructed from the original model by a variant of the Gram-Schmidt orthogonalization process for a rectangular field of view. The relative significance of orthonormal distortion terms is assessed simply by the numerical values of the corresponding coefficients. Orthonormal distortion models are easily extendable when the distribution of residuals indicate the presence of higher order terms.

The daily variation of the solar photocenter over some 11 years is derived from the Mount Wilson data reprocessed by Ulrich et al. 2010 to closely match the surface distribution of solar irradiance. The standard deviations of astrometric jitter are 0.52 $\mu$AU and 0.39 $\mu$AU in the equatorial and the axial dimensions, respectively. The overall dispersion is strongly correlated with the solar cycle, reaching $0.91 \mu$AU at the maximum activity in 2000. The largest short-term deviations from the running average (up to 2.6 $\mu$AU) occur when a group of large spots happen to lie on one side with respect to the center of the disk. The amplitude spectrum of the photocenter variations never exceeds 0.033 $\mu$AU for the range of periods 0.6--1.4 yr, corresponding to the orbital periods of planets in the habitable zone. Astrometric detection of Earth-like planets around stars as quiet as the Sun is not affected by star spot noise, but the prospects for more active stars may be limited to giant planets.

The published Mount Wilson Doppler-shift measurements of the solar velocity field taken in 1967--1982 are revisited with a more accurate model, which includes two terms representing the meridional flow and three terms corresponding to the convective limb shift. Integration of the recomputed data over the visible hemisphere reveals significant variability of the net radial velocity at characteristic time scales of 0.1--10 years, with a standard deviation of 1.4 \ms. This result is supported by independent published observations. The implications for exoplanet detection include reduced sensitivity of the Doppler method to Earth-like planets in the habitable zone, and an elevated probability of false detections at periods of a few to several years.

The third US Naval Observatory (USNO) CCD Astrograph Catalog, UCAC3 was released at the IAU General Assembly on 2009 August 10. It is the first all-sky release in this series and contains just over 100 million objects, about 95 million of them with proper motions, covering about R = 8 to 16 magnitudes. Current epoch positions are obtained from the observations with the 20 cm aperture USNO Astrograph's "red lens", equipped with a 4k by 4k CCD. Proper motions are derived by combining these observations with over 140 ground- and space-based catalogs, including Hipparcos/Tycho and the AC2000.2, as well as unpublished measures of over 5000 plates from other astrographs. For most of the faint stars in the Southern Hemisphere the Yale/San Juan first epoch plates from the SPM program (YSJ1) form the basis for proper motions. These data are supplemented by all-sky Schmidt plate survey astrometry and photometry obtained from the SuperCOSMOS project, as well as 2MASS near-IR photometry. Major differences of UCAC3 data as compared to UCAC2 include a completely new raw data reduction with improved control over systematic errors in positions, significantly improved photometry, slightly deeper limiting magnitude, coverage of the north pole region, greater completeness by inclusion of double stars and weak detections. This of course leads to a catalog which is not as "clean" as UCAC2 and problem areas are outlined for the user in this paper. The positional accuracy of stars in UCAC3 is about 15 to 100 mas per coordinate, depending on magnitude, while the errors in proper motions range from 1 to 10 mas/yr depending on magnitude and observing history, with a significant improvement over UCAC2 achieved due to the re-reduced SPM data and inclusion of more astrograph plate data unavailable at the time of UCAC2.

Analytical relations are derived for the amplitude of astrometric, photometric and radial velocity perturbations caused by a single rotating spot. The relative power of the star spot jitter is estimated and compared with the available data for $\kappa^1$ Ceti and HD 166435, as well as with numerical simulations for $\kappa^1$ Ceti and the Sun. A Sun-like star inclined at $i=90\degr$ at 10 pc is predicted to have a RMS jitter of 0.087 \uas in its astrometric position along the equator, and 0.38 m s$^{-1}$ in radial velocities. If the presence of spots due to stellar activity is the ultimate limiting factor for planet detection, the sensitivity of SIM Lite to Earth-like planets in habitable zones is about an order of magnitude higher that the sensitivity of prospective ultra-precise radial velocity observations of nearby stars.

Luminous X-ray stars are very often found in visual double or multiple stars. Binaries with periods of a few days possess the highest degree of coronal X-ray activity among regular, non-relativistic stars because of their fast, tidally driven rotation. But the orbital periods in visual double stars are too large for any direct interaction between the companions to take place. We suggest that most of the strongest X-ray components in resolved binaries are yet-undiscovered short-period binaries, and that a few are merged remnants of such binaries. The omnipresence of short-period active stars, e.g. of BY-Dra-type binaries, in multiple systems is explained via the dynamical evolution of triple stars with large mutual inclinations. The dynamical perturbation on the inner pair pumps up the eccentricity in a cyclic manner, a phenomenon known as Kozai cycling. At times of close periapsis, tidal friction reduces the angular momentum of the binary, causing it to shrink. When the orbital period of the inner pair drops to a few days, fast surface rotation of the companions is driven by tidal forces, boosting activity by a few orders of magnitude. If the period drops still further, a merger may take place leaving a rapidly-rotating active dwarf with only a distant companion.

We ask if Earth-like planets (terrestrial mass and habitable-zone orbit) can be detected in multi-planet systems, using astrometric and radial velocity observations. We report here the preliminary results of double-blind calculations designed to answer this question.

Abridged: The golden age of astrophysics is upon us with both grand discoveries (extra-solar planets, dark matter, dark energy) and precision cosmology. Fundamental understanding of the working of stars and galaxies is within reach, thanks to newly available precision measurements. We highlight the importance of distances and model independent distances and masses. Distances are fundamental in astrophysics and their knowledge can change our perception of phenomena dramatically: e.g., in antiquity, the Heliocentric model was rejected because the predicted stellar parallaxes were not observed. Distance measurements are directly related to the history & fate of the uni- verse as they provide 2 of 3 methods available to date the universe. The 1st method is based on the ages of stars, which can be ascertained if their lumi- nosities (distances) are accurately known. The 2nd method relies on cosmolo- gical methods. To 1st order, the age of the universe is the inverse H_0. As stressed by the previous decadal report, "the fundamental goal of ... astrophysics is to understand how the universe ... galaxies [and] stars ... formed, how they evolved, and what their destiny will be." These questions can be answered partly by micro-arcsecond astrometry: 1) Galactic archeology: a detailed reconstruction of the formation history of the Milky Way and other Local Group galaxies, 2) the very oldest stars in the Milky Way and the age of the Universe, and 3) H_0 and concordance cosmology. These goals are achievable by combining muas-arcsecond astrometry from the proposed SIM-Lite mission supplemented with ground-based spectroscopy. The results of our proposed project will force the biggest reassessment of stellar astrophysics in 50 years, which will affect most branches of astrophysics.

A set of 41 nearby stars (closer than 25 pc) is investigated which have very wide binary and common proper motion (CPM) companions at projected separations between 1000 and $200 000$ AU. These companions are identified by astrometric positions and proper motions from the NOMAD catalog. Based mainly on measures of chromospheric and X-ray activity, age estimation is obtained for most of 85 identified companions. Color -- absolute magnitude diagrams are constructed to test if CPM companions are physically related to the primary nearby stars and have the same age. Our carefully selected sample includes three remote white dwarf companions to main sequence stars and two systems (55 Cnc and GJ 777A) of multiple planets and distant stellar companions. Ten new CPM companions, including three of extreme separations, are found. Multiple hierarchical systems are abundant; more than 25% of CPM components are spectroscopic or astrometric binaries or multiples themselves. Two new astrometric binaries are discovered among nearby CPM companions, GJ 264 and HIP 59000 and preliminary orbital solutions are presented. The Hyades kinematic group (or stream) is presented broadly in the sample, but we find few possible thick disk objects and none halo stars. It follows from our investigation that moderately young (age $\lesssim 1$ Gyr) thin disk dwarfs are the dominating species in the near CPM systems, in general agreement with the premises of the dynamical survival paradigm. Some of the multiple stellar systems with remote CPM companions probably undergo the dynamical evolution on non-coplanar orbits, known as the Kozai cycle.

Precision astrometry at microarcsecond accuracy has application to a wide range of astrophysical problems. This paper is a study of the science questions that can be addressed using an instrument that delivers parallaxes at about 4 microarcsec on targets as faint as V = 20, differential accuracy of 0.6 microarcsec on bright targets, and with flexible scheduling. The science topics are drawn primarily from the Team Key Projects, selected in 2000, for the Space Interferometry Mission PlanetQuest (SIM PlanetQuest). We use the capabilities of this mission to illustrate the importance of the next level of astrometric precision in modern astrophysics. SIM PlanetQuest is currently in the detailed design phase, having completed all of the enabling technologies needed for the flight instrument in 2005. It will be the first space-based long baseline Michelson interferometer designed for precision astrometry. SIM will contribute strongly to many astronomical fields including stellar and galactic astrophysics, planetary systems around nearby stars, and the study of quasar and AGN nuclei. SIM will search for planets with masses as small as an Earth orbiting in the `habitable zone' around the nearest stars using differential astrometry, and could discover many dozen if Earth-like planets are common. It will be the most capable instrument for detecting planets around young stars, thereby providing insights into how planetary systems are born and how they evolve with time. SIM will observe significant numbers of very high- and low-mass stars, providing stellar masses to 1%, the accuracy needed to challenge physical models. Using precision proper motion measurements, SIM will probe the galactic mass distribution and the formation and evolution of the Galactic halo. (abridged)

We report a late M-type, common proper motion companion to a nearby young visual binary HIP 115147 (V368 Cep), separated by 963 arcseconds from the primary K0 dwarf. This optically dim star has been identified as a candidate high proper motion, nearby dwarf LSPM J2322+7847 by L{\'e}pine in 2005. The wide companion is one of the latest post-T Tauri low mass stars found within 20 pc. We obtain a trigonometric parallax of $51.6\pm0.8$ mas, in good agreement with the Hipparcos parallax of the primary star ($50.7\pm0.6$ mas). Our $BVRI$ photometric data and near-infrared data from 2MASS are consistent with LSPM J2322+7847 being brighter by 1 magnitude in $K_s$ than field M dwarfs at $V-K_s=6.66$, which indicates its pre-main sequence status. We conclude that the most likely age of the primary HIP 115147 and its 11-arcsecond companion HIP 115147B is 20-50 Myr. The primary appears to be older than its close analog PZ Tel (age 12-20 Myr) and members of the TWA association (7 Myr).

We present further arguments that the Hipparcos parallaxes for some of the clusters and associations represented in the Hipparcos catalog should be used with caution in the study of the Galactic structure. It has been already shown that the discrepancy between the Hipparcos and ground based parallaxes for several clusters including the Pleiades, Coma Ber and NGC 6231 can be resolved by recomputing the Hipparcos astrometric solutions with an improved algorithm diminishing correlated errors in the attitude parameters. Here we present new parallaxes obtained with this algorithm for another group of stars with discrepant data - the galactic cluster Cr 121. The original Hipparcos parallaxes led de Zeeuw et al. to conclude that Cr 121 and the surrounding association of OB stars form a relatively compact and coherent moving group at a distance of 550 -- 600 pc. Our corrected parallaxes reveal a different spatial distribution of young stellar populace in this area. Both the cluster Cr 121 and the extended OB association are considerably more distant (750 -- 1000 pc), and the latter has a large depth probably extending beyond 1 kpc. Therefore, not only are the recalculated parallaxes in complete agreement with the photometric uvbybeta parallaxes, but the structure of the field they reveal is no longer in discrepancy with that found by the photometric method.

Taking advantage of an improved genetic optimization algorithm for fitting unconstrained Kepler orbits to the Hipparcos intermediate astrometric data, we obtain additional orbital solutions for 81 Hipparcos stars with previous stochastic (failed) solutions. The sample includes astrophysically interesting objects, including the old disk wide binary HIP 754, the nearby AGB star HIP 34922 (L2 Pup), and the nearby M2 dwarf HIP 5496 (GJ 54, at 8 pc from the Sun) which has a resolved M dwarf companion.

We analyze the local field of stellar tangential velocities for a sample of $42 339$ non-binary Hipparcos stars with accurate parallaxes, using a vector spherical harmonic formalism. We derive simple relations between the parameters of the classical linear model (Ogorodnikov-Milne) of the local systemic field and low-degree terms of the general vector harmonic decomposition. Taking advantage of these relationships we determine the solar velocity with respect to the local stars of $(V_X,V_Y,V_Z)=(10.5, 18.5, 7.3)\pm 0.1$ \kms. The Oort's parameters determined by a straightforward least-squares adjustment in vector spherical harmonics, are $A=14.0\pm 1.4$, $B=-13.1\pm 1.2$, $K=1.1\pm 1.8$, and $C=-2.9\pm 1.4$ \kmspc. We find a few statistically significant higher degree harmonic terms, which do not correspond to any parameters in the classical linear model. One of them, a third-degree electric harmonic, is tentatively explained as the response to a negative linear gradient of rotation velocity with distance from the Galactic plane, which we estimate at $\sim -20$ \kmspc. The most unexpected and unexplained term within the Ogorodnikov-Milne model is the first-degree magnetic harmonic representing a rigid rotation of the stellar field about the axis $-Y$ pointing opposite to the direction of rotation. This harmonic comes out with a statistically robust coefficient $6.2 \pm 0.9$ \kmspc, and is also present in the velocity field of more distant stars. The ensuing upward vertical motion of stars in the general direction of the Galactic center and the downward motion in the anticenter direction are opposite to the vector field expected from the stationary Galactic warp model.

General relativistic deflection of light by mass, dipole, and quadrupole moments of gravitational field of a moving massive planet in the Solar system is derived. All terms of order 1 microarcsecond are taken into account, parametrized, and classified in accordance with their physical origin. We calculate the instantaneous patterns of the light-ray deflections caused by the monopole, the dipole and the quadrupole moments, and derive equations describing apparent motion of the deflected position of the star in the sky plane as the impact parameter of the light ray with respect to the planet changes due to its orbital motion. The present paper gives the physical interpretation of the observed light-ray deflections and discusses the observational capabilities of the near-future optical (SIM) and radio (SKA) interferometers for detecting the Doppler modulation of the radial deflection, and the dipolar and quadrupolar light-ray bendings by the Jupiter and the Saturn.

A considerable number of astrometric binaries whose positions on the sky do not obey the standard model of mean position, parallax and linear proper motion, were observed by the Hipparcos satellite. Some of them remain non-discovered, and their observational data have not been properly processed with the more adequate astrometric model that includes nonlinear orbital motion. We develop an automated algorithm based on "genetic optimization", to solve the orbital fitting problem with no prior information about the orbital elements is available (from, e.g., spectroscopic data or radial velocity monitoring). We test this method on Hipparcos stars with known orbital solutions in the catalog, and further apply it to stars with stochastic solutions, which may be unresolved binaries. At a confidence level of 99%, orbital fits are obtained for 65 stars, most of which have not been known as binary. A few of the new probable binaries with A-type primaries with periods 444-2015 d are chemically peculiar stars, including Ap and \lambda Boo type. The anomalous spectra of these stars are explained as admixture of the light from the unresolved, sufficiently bright and massive companions. We estimate the apparent orbits of four stars which have been identified as members of the 300 Myr-old UMa kinematic group. Another four new nearby binaries may include low-mass M-type or brown dwarf companions. Similar astrometric models and algorithms can be used for binary stars and planet hosts observed by SIM PlanetQuest and Gaia.

One of the main endeavors of fundamental astrometry is to establish a practical realization of an inertial reference frame anchored to celestial objects whose positions are defined in the barycentric coordinates of the solar system. The development of astrometric facilities operating from space at a microarcsecond level of precision makes the non-uniformity of the galactic motion of the barycenter an observable effect that violates the inertiality of the barycentric frame. Most of the observable effect is caused by secular acceleration of the barycenter with respect to the center of the Galaxy. The acceleration results in a pattern of secular aberration which is observable astrometrically as a systematic vector field of the proper motions of distant quasars. We employ the vector spherical harmonics to describe the predicted field of the proper motions and evaluate its amplitude at each point on the celestial sphere. It is shown that the pattern of secular aberration is fully represented by three low-order electric-type vector harmonics, and hence, it is easily distinguishable from the residual rotations of the reference frame and other possible effects, such as the hypothetical long-period gravitational waves. Comprehensive numerical simulations of the grid astrometry with SIM PlanetQuest are conducted. The full covariance matrix of the simulated grid solution is used to evaluate the covariances of the three electric harmonic coefficients, representing the secular aberration pattern of proper motions. We conclude that the grid astrometry with SIM PlanetQuest will be sensitive to the main galactocentric component of secular acceleration, while the peculiar acceleration of the Sun with respect to LSR is expected to be too small to be detected with this astrometric space interferometer.

Hipparcos observations of some variable stars, and especially of long-period (e.g. Mira) variables, reveal a motion of the photocenter correlated with the brightness variation ({variability-induced mover -- VIM), suggesting the presence of a binary companion. A re-analysis of the Hipparcos photometric and astrometric data does not confirm the VIM solution for 62 among the 288 VIM objects (21%) in the Hipparcos catalogue. Most of these 288 VIMs are long-period (e.g. Mira) variables (LPV). The effect of a revised chromaticity correction, which accounts for the color variations along the light cycle, was then investigated. It is based on `instantaneous' $V-I$ color indices derived from Hipparcos and Tycho-2 epoch photometry. Among the 188 LPVs flagged as VIM in the Hipparcos catalogue, 89 (47%) are not confirmed as VIM after this improved chromaticity correction is applied. This dramatic decrease in the number of VIM solutions is not surprising, since the chromaticity correction applied by the Hipparcos reduction consortia was based on a fixed $V-I$ color. Astrophysical considerations lead us to adopt a more stringent criterion for accepting a VIM solution (first-kind risk of 0.27% instead of 10% as in the Hipparcos catalogue). With this more severe criterion, only 27 LPV stars remain VIM, thus rejecting 161 of the 188 (86%) of the LPVs defined as VIMs in the Hipparcos catalogue.

For Hipparcos M, S, and C spectral type stars, we provide calibrated instantaneous (epoch) Cousins $V-I$ color indices using newly derived $HpV_{T2}$ photometry. Three new sets of ground-based Cousins $VI$ data have been obtained for more than 170 carbon and red M giants. These datasets in combination with the published sources of $VI$ photometry served to obtain the calibration curves linking Hipparcos/Tycho $Hp-V_{T2}$ with the Cousins $V-I$ index. In total, 321 carbon stars and 4464 M- and S-type stars have new $V-I$ indices. The standard error of the mean $V-I$ is about 0.1 mag or better down to $Hp\approx9$ although it deteriorates rapidly at fainter magnitudes. These $V-I$ indices can be used to verify the published Hipparcos $V-I$ color indices. Thus, we have identified a handful of new cases where, instead of the real target, a random field star has been observed. A considerable fraction of the DMSA/C and DMSA/V solutions for red stars appear not to be warranted. Most likely such spurious solutions may originate from usage of a heavily biased color in the astrometric processing.

We consider the propagation of galactic cosmic rays in the fractal interstellar medium. Steady state solution of the fractional diffusion equation, describing cosmic ray propagation, is found. We show that the exponent of the steady state spectrum turns out to be equal to the exponent of the cosmic ray spectrum above the "knee", that is approximately equal to 3.1.

Accuracy in the absolute position in the sky is one of the limiting factors for pulsar timing, and timing parameters have a direct impact on the understanding of the physics of Isolated Neutron Stars (INS). We report here on a high-accuracy measurement of the optical position of Geminga (V=25.5), the only known radio-quiet INS. The procedure combines the Hipparcos and Tycho catalogues, ground-based astrometric data,and Hubble Space Telescope (HST) Wide Field Planetary Camera (WFPC2) images, to yield Geminga's absolute position to within ~40 mas (per coordinate). Such a positional accuracy, unprecedented for the optical position of a pulsar or an object this faint, is needed to combine in phase gamma-ray photons collected over more than 20 years, i.e. over 2.5 billions of star' revolutions. Although quite a difficult task, this is the only way to improve our knowledge of the timing parameters of this radio silent INS.