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.

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.

Perturbations of all sorts destabilise the rotation of a small body and leave it in a non-principal spin state. In such a state, the body experiences alternating stresses generated by the inertial forces. This yields nutation relaxation, i.e., evolution of the spin towards the principal rotation about the maximal-inertia axis. Knowledge of the timescales needed to damp the nutation is crucial in studies of small bodies' dynamics. In the literature hitherto, nutation relaxation has always been described with aid of an empirical quality factor $\,Q\,$ introduced to parameterise the energy dissipation rate. Among the drawbacks of this approach was its inability to describe the dependence of the relaxation rate upon the current nutation angle. This inability stemmed from our lack of knowledge of the quality factor's dependence on the forcing frequency. In this article, we derive our description of nutation damping directly from the rheological law obeyed by the material. This renders us the nutation damping rate as a function of the current nutation angle, as well as of the shape and the rheological parameters of the body. In contradistinction from the approach based on an empirical $\,Q\,$-factor, our development gives a zero damping rate in the spherical-shape limit. Our method is generic and applicable to any shape and to any linear rheological law. However, to simplify the developments, here we consider a dynamically oblate rotator with a Maxwell rheology.

The Darwin-Kaula theory of bodily tides is intended for celestial bodies rotating without libration. We demonstrate that this theory, in its customary form, is inapplicable to a librating body. Specifically, in the presence of libration in longitude, the actual spectrum of Fourier tidal modes differs from the conventional spectrum rendered by the Darwin-Kaula theory for a non-librating celestial object. This necessitates derivation of formulae for the tidal torque and the tidal heating rate, that are applicable under libration. We derive the tidal spectrum for longitudinal forced libration with one and two main frequencies, generalisation to more main frequencies being straightforward. (By main frequencies we understand those emerging due to the triaxiality of the librating body.) Separately, we consider a case of free libration at one frequency (once again, generalisation to more frequencies being straightforward). We also calculate the tidal torque. This torque provides correction to the triaxiality-caused physical libration. Our theory is not self-consistent: we assume that the tidal torque is much smaller than the permanent-triaxiality-caused torque; so the additional libration due to tides is much weaker than the main libration due to the permanent triaxiality. Finally, we calculate the tidal dissipation rate in a body experiencing forced libration at the main mode, or free libration at one frequency, or superimposed forced and free librations.

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

New astrometric reductions of the US Naval Observatory CCD Astrograph Catalog (UCAC) all-sky observations were performed from first principles using the TGAS stars in the 8 to 11 magnitude range as reference star catalog. Significant improvements in the astrometric solutions were obtained and the UCAC5 catalog of mean positions at a mean epoch near 2001 was generated. By combining UCAC5 with Gaia DR1 data new proper motions on the Gaia coordinate system for over 107 million stars were obtained with typical accuracies of 1 to 2 mas/yr (R = 11 to 15 mag), and about 5 mas/yr at 16th mag. Proper motions of most TGAS stars are improved over their Gaia data and the precision level of TGAS proper motions is extended to many millions more, fainter stars. External comparisons were made using stellar cluster fields and extragalactic sources. The TGAS data allow us to derive the limiting precision of the UCAC x,y data, which is significantly better than 1/100 pixel.

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 use numerical simulations to measure the sensitivity of the tidal spin down rate of a homogeneous triaxial ellipsoid to its axis ratios by comparing the drift rate in orbital semi-major axis to that of a spherical body with the same mass, volume and simulated rheology. We use a mass-spring model approximating a viscoelastic body spinning around its shortest body axis, with spin aligned with orbital spin axis, and in circular orbit about a point mass. The torque or drift rate can be estimated from that predicted for a sphere with equivalent volume if multiplied by $0.5 (1 + b^4/a^4)(b/a)^{-4/3} (c/a)^{-\alpha_c}$ where $b/a$ and $c/a$ are the body axis ratios and index $\alpha_c \approx 1.05$ is consistent with the random lattice mass spring model simulations but $\alpha_c = 4/3$ suggested by scaling estimates. A homogeneous body with axis ratios 0.5 and and 0.8, like Haumea, has orbital semi-major axis drift rate about twice as fast as a spherical body with the same mass, volume and material properties. A simulation approximating a mostly rocky body but with 20\% of its mass as ice concentrated at its ends has a drift rate 10 times faster than the equivalent homogeneous rocky sphere. However, this increase in drift rate is not enough to allow Haumea's satellite, Hi'iaka, to have tidally drifted away from Haumea to its current orbital semi-major axis.

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.

We use a damped mass-spring model within an N-body code to simulate the tidal evolution of the spin and orbit of a self-gravitating viscoelastic spherical body moving around a point-mass perturber. The damped mass-spring model represents a Kelvin-Voigt viscoelastic solid. We measure the tidal quality function (the dynamical Love number $\,k_2\,$ divided by the tidal quality factor $\,Q\,$) from the numerically computed tidal drift of the semimajor axis of the binary. The shape of $\,k_2/Q\,$, as a function of the principal tidal frequency, reproduces the kink shape predicted by Efroimsky (2012a; CeMDA 112$\,:\,$283) for the tidal response of near-spherical homogeneous viscoelastic rotators. We demonstrate that we can directly simulate the tidal evolution of spinning viscoelastic objects. In future, the mass-spring N-body model can be generalised to inhomogeneous and/or non-spherical bodies.

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.

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 study the stability of the (87) Sylvia system and of the neighborhood of its two satellites. We use numerical integrations considering the non-sphericity of Sylvia, as well as the mutual perturbation of the satellites and the solar perturbation. Two numerical models have been used, which describe respectively the short and long-term evolution of the system. We show that the actual system is in a deeply stable zone, but surrounded by both fast and secular chaotic regions due to resonances. We then investigate how tidal and BYORP effects modify the location of the system over time with respect to the instability zones. Finally, we briefly generalize this study to other known triple systems and to satellites of asteroids in general, and discuss about their distance from mean-motion and evection resonances.