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.

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.

[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.

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.

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.

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.

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 in the approximation of the linearized Einstein equations. All terms of order 1 microarcsecond are taken into account, parametrized, and classified in accordance with their physical origin. We discuss the observational capabilities of the near-future optical and radio interferometers for detecting the Doppler modulation of the radial deflection, and the dipolar and quadrupolar light-ray bendings by Jupiter and the Saturn.

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.