Observations show that 100-km-class Kuiper belt objects (KBOs) can be divided in (at least) two color groups, hereafter red (R, g-i<1.2) and very red (VR, g-i>1.2), reflecting a difference in their surface composition. This is thought to imply that KBOs formed over a relatively wide range of radial distance, r. The cold classicals at 42<r<47 au are predominantly VR and known Neptune Trojans at r=30 au are mostly R. Intriguingly, however, the dynamically hot KBOs show a mix of R and VR colors and no correlation of color with r. Here we perform migration/instability simulations where the Kuiper belt is populated from an extended planetesimal disk. We find that the color observations can be best understood if R objects formed at r<r* and VR objects at r>r*, with 30<r*<40 au. The proposed transition at 30<r*<40 au would explain why the VR objects in the dynamically hot population have smaller orbital inclinations than the R objects, because the orbital excitation from Neptune weakens for orbits starting beyond 30 au. Possible causes of the R-VR color bimodality are discussed.

The Colombo top is a basic model in the rotation dynamics of a celestial body moving on a precessing orbit and perturbed by a gravitational torque. The paper presents a detailed study of analytical solution to this problem. By solving algebraic equations of degree 4, we provide the expressions for the extreme points of trajectories as functions of their energy. The location of stationary points (known as the Cassini states) is found as the function of the two parameters of the problem. Analytical solution in terms the Weierstrass and the Jacobi elliptic functions is given for regular trajectories. Some trajectories are expressible through elementary functions: not only the homoclinic orbits, as expected, but also a special periodic solution whose energy is equal to that of the first Cassini state (unnoticed in previous studies).

A recent examination of K2 lightcurves indicates that ~15% of Jupiter Trojans have very slow rotation (spin periods Ps>100 h). Here we consider the possibility that these bodies formed as equal-size binaries in the massive outer disk at ~20-30 au. Prior to their implantation as Jupiter Trojans, tight binaries tidally evolved toward a synchronous state with Ps~Pb, where Pb is the binary orbit period. They may have been subsequently dissociated by impacts and planetary encounters with at least one binary component retaining its slow rotation. Surviving binaries on Trojan orbits would continue to evolve by tides and spin-changing impacts over 4.5 Gyr. To explain the observed fraction of slow rotators, we find that at least ~15-20% of outer disk bodies with diameters 15<D<50 km would have to form as equal-size binaries with 12<ab/R<30, where ab is the binary semimajor axis and R=D/2. The mechanism proposed here could also explain very slow rotators found in other small body populations.

We studied asteroid clusters suggesting a possibility of at least two disruption events in their recent history (< 5 Myr). We searched for new members of known asteroid pairs and clusters and we verified their membership. We found four asteroid clusters, namely the clusters of (11842) Kap'bos, (14627) Emilkowalski, (63440) 2001 MD30 and (157123) 2004 NW5 that show at least two secondary separation events that occurred at significantly different times. We considered a possible formation mechanism for these clusters: The parent of an asteroid cluster was spun up to its critical rotation frequency, underwent a rotation fission and was slowed down by escape of the newly formed secondary/ies. Then the YORP effect spun up the primary again and it reached its critical rotation frequency and underwent another fission. We created a simple model to test whether this scenario is possible for the four clusters. We obtained a good agreement between the model and the cluster properties for the clusters of Kap'bos and (63440). For the cluster of Emilkowalski, the time needed for the primary to reach its critical frequency after the first fission event was predicted to be too long by a factor of several. We suspect, considering also its D type taxonomic classification and the existence of a dust band associated with the cluster, that the asteroid Emilkowalski may actually be a cometary nucleus. For the cluster of (157123), the final rotational frequency of the primary after the last fission event predicted by our model is in a good agreement with the observed rotation frequency of (157123). However, a separation of the older secondary is not possible in our model due to the deficiency of free energy needed for an escape of the large secondary. This could be due to an error in the H value of the secondary or the possibility that we did not find the real primary of this cluster.

Aims: Our goal was to increase number of known doubly eclipsing systems such that the resulting dataset would allow to study them via statistical means, as well as prove that they constitute gravitationally bound 2+2 quadruple system. Methods: We analysed photometric data for eclipsing binaries provided by the OGLE survey in the LMC fields. We found a large number of new doubly eclipsing systems (our discoveries are 3x more numerous than previous studies). With a typical orbital period of days for the binaries, we sought eclipse time variations (ETVs) on the timescale of years. In the cases where we were able to detect the ETV period, the difference between the inner and outer periods in the quadruple system is large enough. This allows us to interpret ETVs primarily as the light-time effect, thus providing an interesting constraint on masses of the binaries. Results: In addition to significantly enlarging the database of known doubly eclipsing systems, we performed a thorough analysis of 72 cases. ETVs for 28 of them (39% of the studied cases) showed evidence of relative motion. We note OGLE BLG-ECL-145467 as the most interesting case; it is bright (I=12.6 mag), consists of two detached binaries with periods of about 3.3 d and 4.9 d (making it a candidate for a 3:2 resonant system), mutual period about 1538 d. Distribution of the orbital period ratio P_A/P_B of binaries in 2+2 quadruples shows statistically significant excess at 1 and 1.5. The former is likely a natural statistical preference in weakly interacting systems with periods within the same range. The latter is thought to be evidence of a capture in the 3:2 mean motion resonance of the two binaries. This sets important constraints on evolutionary channels in these systems. The total number of doubly eclipsing systems increased to 146, more than 90% of which are at low declinations on the southern sky.

We use published models of the early Solar System evolution with a slow, long-range and grainy migration of Neptune to predict the orbital element distributions and the number of modern-day Centaurs. The model distributions are biased by the OSSOS survey simulator and compared with the OSSOS Centaur detections. We find an excellent match to the observed orbital distribution, including the wide range of orbital inclinations which was the most troublesome characteristic to fit in previous models. A dynamical model, in which the original population of outer disk planetesimals was calibrated from Jupiter Trojans, is used to predict that OSSOS should detect 11+/-4 Centaurs with semimajor axis a<30 au, perihelion distance q>7.5 au and diameter D>10 km (absolute magnitude H_r<13.7 for a 6% albedo). This is consistent with 15 actual OSSOS Centaur detections with H_r<13.7. The population of Centaurs is estimated to be 21,000+/-8,000 for D>10 km. The inner scattered disk at 50<a<200 au should contain (2.0+/-0.8)x10^7 D>10 km bodies and the Oort cloud should contain (5.0+/-1.9)x10^8 D>10 km comets. Population estimates for different diameter cutoffs can be obtained from the size distribution of Jupiter Trojans (N(>D) proportional to D^(-2.1) for 5<D<100 km). We discuss model predictions for the Large Synoptic Survey Telescope observations of Centaurs.

As indicated by their special characteristics, the cold classical Kuiper belt objects (KBOs) formed and survived at 42-47 au. Notably, they show a large fraction of equal-size binaries whose formation is probably related to the accretion of KBOs themselves. These binaries are uncommon in other --hot, resonant, scattered-- populations, which are thought to have been implanted from the massive disk below 30 au to >30 au during Neptune's migration. Here we highlight the possibility that equal-size binaries formed in the disk but were subsequently removed by impacts and/or dynamical effects (e.g., scattering encounters with Neptune). We determine the dependence of these processes on the size and separation of binary components. Our results indicate that tighter binaries, if they formed in the massive disk, have relatively good chances of survival (unless the disk was long-lived). In contrast, the widest binaries in the hot population, such as 2002 VF130, have a very low survival probability (<1%) even if the massive disk was short-lived. They may represent a trace of lucky survivors of a much larger population of the original disk binaries, or they formed at ~30-40~au and dodged the impact- and encounter-related perturbations that we studied here. We find that all known satellites of the largest KBOs would survive during the dynamical implantation of these bodies in the Kuiper belt. The low orbital eccentricities of Pluto's small moons may have been excited by impacts and/or encounters of the Pluto system to Neptune.

The recent close approach of the NEA (3200) Phaethon offered a rare opportunity to obtain high-quality observational data. We used the newly obtained optical light curves to improve the spin and shape model of Phaethon and to determine its surface physical properties derived by thermophysical modeling. We also used the available astrometric observations of Phaethon, including those obtained by the Arecibo radar and the Gaia spacecraft, to constrain the secular drift of the orbital semimajor axis. This constraint allowed us to estimate the bulk density by assuming that the drift is dominated by the Yarkovsky effect. We used the convex inversion model to derive the 3D shape model of Phaethon, and a detailed numerical approach for an accurate analysis of the Yarkovsky effect. We obtained a unique solution for Phaethon's pole orientation at $(318,-47)^{\circ}$ ecliptic longitude and latitude (uncertainty of $5^{\circ}$), and confirm the previously reported thermophysical properties ($D=5.1\pm0.2$ km, $\Gamma=600\pm200$ SI). Phaethon has a top-like shape with possible north-south asymmetry. The characteristic size of the regolith grains is 1-2 cm. The orbit analysis reveals a secular drift of the semimajor axis of $-(6.9\pm1.9)\times 10^{-4}$ au Myr$^{-1}$. With the derived volume-equivalent size of 5.1~km, the bulk density $\rho$ is $1.67\pm0.47$ g cm$^{-3}$. If the size is slightly larger $\sim5.7$ km, as suggested by radar data, $\rho$ would decrease to $1.48\pm0.42$ g cm$^{-3}$. We further investigated the suggestion that Phaethon may be in a cluster with asteroids (155140) 2005 UD and (225416) 1999 YC that was formed by rotational fission of a critically spinning parent body. Phaethon's $\rho$ is consistent with typical values for large ($>100$ km) C-complex asteroids and supports its association with asteroid (2) Pallas. These findings render a cometary origin unlikely for Phaethon.

The orbital distribution of trans-Neptunian objects provides strong evidence for the radial migration of Neptune. The outer planets' orbits are thought to have become unstable during the early stages with Jupiter having scattering encounters with a Neptune-class planet. As a consequence, Jupiter jumped inward by a fraction of an au, as required from inner solar system constraints, and obtained its current orbital eccentricity. The timing of these events is often linked to the lunar Late Heavy Bombardment that ended ~700 Myr after the dispersal of the protosolar nebula (t_0). Here we show instead that planetary migration started shortly after t_0. Such early migration is inferred from the survival of the Patroclus-Menoetius binary Jupiter Trojan. The binary formed at t <~ t_0 within a massive planetesimal disk once located beyond Neptune. The longer the binary stayed in the disk, the greater the likelihood that collisions would strip its components from one another. The simulations of its survival indicate that the disk had to have been dispersed by migrating planets within <100 Myr of t_0. This constraint implies that the planetary migration is unrelated to the formation of the youngest lunar basins.

The debiased absolute-magnitude and orbit distributions as well as source regions for near-Earth objects (NEOs) provide a fundamental frame of reference for studies of individual NEOs and more complex population-level questions. We present a new four-dimensional model of the NEO population that describes debiased steady-state distributions of semimajor axis, eccentricity, inclination, and absolute magnitude $H$ in the range $17<H<25$. The modeling approach improves upon the methodology originally developed by Bottke et al. (2000; Science 288, 2190-2194) in that it is, for example, based on more realistic orbit distributions and uses source-specific absolute-magnitude distributions that allow for a power-law slope that varies with $H$. We divide the main asteroid belt into six different entrance routes or regions (ER) to the NEO region: the $\nu_6$, 3:1J, 5:2J and 2:1J resonance complexes as well as Hungarias and Phocaeas. In addition we include the Jupiter-family comets as the primary cometary source of NEOs. We calibrate the model against NEO detections by Catalina Sky Surveys' stations 703 and G96 during 2005-2012, and utilize the complementary nature of these two systems to quantify the systematic uncertainties associated to the resulting model. The most important ERs are the $\nu_6$ and 3:1J resonance complexes with JFCs contributing a few percent of NEOs on average. We predict that there are $962^{+52}_{-56}$ ($802^{+48}_{-42}\times10^3$) NEOs with $H<17.75$ ($H<25$) and these numbers are in agreement with the most recent estimates found in the literature (the uncertainty estimates only account for the random component).

The Rosetta spacecraft observations revealed that the nucleus of comet 67P/Churyumov-Gerasimenko consists of two similarly sized lobes connected by a narrow neck. Here we evaluate the possibility that 67P is a collapsed binary. We assume that the progenitor of 67P was a binary and consider various physical mechanisms that could have brought the binary components together, including small-scale impacts and gravitational encounters with planets. We find that 67P could be a primordial body (i.e., not a collisional fragment) if the outer planetesimal disk lasted <10 Myr before it was dispersed by migrating Neptune. The probability of binary collapse by impact is 30% for tightly bound binaries. Most km-class binaries become collisionally dissolved. Roughly 10% of the surviving binaries later evolve to become contact binaries during the disk dispersal, when bodies suffer gravitational encounters to Neptune. Overall, the processes described in this work do not seem to be efficient enough to explain the large fraction (~67%) of bi-lobed cometary nuclei inferred from spacecraft imaging.

On February 6th, 2018 SpaceX launched a Tesla Roadster on a Mars-crossing orbit. We perform N-body simulations to determine the fate of the object over the next 15 Myr. The orbital evolution is initially dominated by close encounters with the Earth. While a precise orbit can not be predicted beyond the next several centuries due to these repeated chaotic scatterings, one can reliably predict the long-term outcomes by statistically analyzing a large suite of possible trajectories with slightly perturbed initial conditions. Repeated gravitational scatterings with Earth lead to a random walk. Collisions with the Earth, Venus and the Sun represent primary sinks for the Roadster's orbital evolution. Collisions with Mercury and Mars, or ejections from the Solar System by Jupiter, are highly unlikely. We calculate a dynamical half-life of the Tesla of approximately 15 Myr, with some 22%, 12% and 12% of Roadster orbit realizations impacting the Earth, Venus, and the Sun within one half-life, respectively. Because the eccentricities and inclinations in our ensemble increase over time due to mean-motion and secular resonances, the impact rates with the terrestrial planets decrease beyond a few million years, whereas the impact rate on the Sun remains roughly constant.

The rotation states of small asteroids are affected by a net torque arising from an anisotropic sunlight reflection and thermal radiation from the asteroids' surfaces. On long timescales, this so-called YORP effect can change asteroid spin directions and their rotation periods. We analyzed lightcurves of four selected near-Earth asteroids with the aim of detecting secular changes in their rotation rates that are caused by YORP. We use the lightcurve inversion method to model the observed lightcurves and include the change in the rotation rate $\mathrm{d} \omega / \mathrm{d} t$ as a free parameter of optimization. We collected more than 70 new lightcurves. For asteroids Toro and Cacus, we used thermal infrared data from the WISE spacecraft and estimated their size and thermal inertia. We also used the currently available optical and radar astrometry of Toro, Ra-Shalom, and Cacus to infer the Yarkovsky effect. We detected a YORP acceleration of $\mathrm{d}\omega / \mathrm{d} t = (1.9 \pm 0.3) \times 10^{-8}\,\mathrm{rad}\,\mathrm{d}^{-2}$ for asteroid Cacus. For Toro, we have a tentative ($2\sigma$) detection of YORP from a significant improvement of the lightcurve fit for a nonzero value of $\mathrm{d}\omega / \mathrm{d} t = 3.0 \times 10^{-9}\,\mathrm{rad}\,\mathrm{d}^{-2}$. For asteroid Eger, we confirmed the previously published YORP detection with more data and updated the YORP value to $(1.1 \pm 0.5) \times 10^{-8}\,\mathrm{rad}\,\mathrm{d}^{-2}$. We also updated the shape model of asteroid Ra-Shalom and put an upper limit for the change of the rotation rate to $|\mathrm{d}\omega / \mathrm{d} t| \lesssim 1.5 \times 10^{-8}\,\mathrm{rad}\,\mathrm{d}^{-2}$. Ra-Shalom has a greater than $3\sigma$ Yarkovsky detection with a theoretical value consistent with observations assuming its size and/or density is slightly larger than the nominally expected values.

Comets are icy objects that orbitally evolve from the trans-Neptunian region (the Kuiper belt and beyond) into the inner Solar System, where they are heated by solar radiation and become active due to sublimation of water ice. Here we perform end-to-end simulations in which cometary reservoirs are formed in the early Solar System and evolved over 4.5 Gyr. The gravitational effects of Planet 9 (P9), hypothesized to circle the Sun on a wide orbit, are included in some of our simulations. Different models are considered for comets to be active, including a simple assumption that comets remain active for Np(q) perihelion passages with perihelion distance q<2.5 au. The orbital distribution and number of active comets produced in our model is compared to observations. The orbital distribution of ecliptic comets (ECs) is well reproduced in models with Np(2.5)=500 and without P9. With P9, the inclination distribution of model ECs is wider than the observed one. We find that the known Halley-type comets (HTCs) have a nearly isotropic inclination distribution (with only a slight preference for prograde orbits). In our model, the HTCs appear to be an extension of the population of returning Oort-cloud comets (OCCs) to shorter orbital periods. The inclination distribution of model HTCs becomes broader with increasing Np, but the existing observational data are not good enough to constrain Np from orbital fits. Np(2.5)>1000 is required to obtain a steady-state population of large active HTCs that is consistent with observations. To fit the ratio of the returning-to-new OCCs, by contrast, our model implies that Np(2.5)<10, possibly because the detected long-period comets are smaller and much easier to disrupt than observed HTCs.

The dynamical structure of the Kuiper belt beyond 50 au is not well understood. Here we report results of a numerical model with long-range, slow and grainy migration of Neptune. The model implies that bodies scattered outward by Neptune to semimajor axes a>50 au often evolve into resonances which subsequently act to raise the perihelion distances of orbits to q>40 au. The implication of the model is that the orbits with 50<a<100 au and q>40 au should cluster near (but not in) the resonances with Neptune (3:1 at a=62.6 au, 4:1 at a=75.9 au, 5:1 at a=88.0 au, etc.). The recent detection of several distant Kuiper Belt Objects (KBOs) near resonances is consistent with this prediction, but it is not yet clear whether the orbits are really non-resonant as our model predicts. We estimate from the model that there should presently be ~1600-2400 bodies at the 3:1 resonance and ~1000-1400 bodies at the 4:1 resonance (for q>40 au and diameters D>100 km). These results favorably compare with the population census of distant KBOs inferred from existing observations.

The near-Earth asteroid (3200) Phaethon is an intriguing object: its perihelion is at only 0.14 au and is associated with the Geminid meteor stream. We aim to use all available disk-integrated optical data to derive a reliable convex shape model of Phaethon. By interpreting the available space- and ground-based thermal infrared data and Spitzer spectra using a thermophysical model, we also aim to further constrain its size, thermal inertia, and visible geometric albedo. We applied the convex inversion method to the new optical data obtained by six instruments and to previous observations. The convex shape model was then used as input for the thermophysical modeling. We also studied the long-term stability of Phaethon's orbit and spin axis with a numerical orbital and rotation-state integrator. We present a new convex shape model and rotational state of Phaethon: a sidereal rotation period of 3.603958(2) h and ecliptic coordinates of the preferred pole orientation of (319$^{\circ}$, $-$39$^{\circ}$) with a 5$^{\circ}$ uncertainty. Moreover, we derive its size ($D$=5.1$\pm$0.2 km), thermal inertia ($\Gamma$=600$\pm$200 J m$^{-2}$ s$^{-1/2}$ K$^{-1}$), geometric visible albedo ($p_{\mathrm{V}}$=0.122$\pm$0.008), and estimate the macroscopic surface roughness. We also find that the Sun illumination at the perihelion passage during the past several thousand years is not connected to a specific area on the surface, which implies non-preferential heating.

The Transit Timing Variations (TTVs) are deviations of the measured mid-transit times from the exact periodicity. One of the most interesting causes of TTVs is the gravitational interaction between planets. Here we consider a case of two planets in a mean motion resonance (orbital periods in a ratio of small integers). This case is important because the resonant interaction can amplify the TTV effect and allow planets to be detected more easily. We develop an analytic model of the resonant dynamics valid for small orbital eccentricities and use it to derive the principal TTV terms. We find that a resonant system should show TTV terms with two basic periods (and their harmonics). The resonant TTV period is proportional (m/M_*)^(-2/3), where m and M_* are the planetary and stellar masses. For m=10^(-4) M_*, for example, the TTV period exceeds the orbital period by ~2 orders of magnitude. The amplitude of the resonant TTV terms scales linearly with the libration amplitude. The ratio of the TTV amplitudes of two resonant planets is inversely proportional to the ratio of their masses. These and other relationships discussed in the main text can be used to aid the interpretation of TTV observations.

The Kuiper belt is a population of icy bodies beyond the orbit of Neptune. The complex orbital structure of the Kuiper belt, including several categories of objects inside and outside of resonances with Neptune, emerged as a result of Neptune's migration into an outer planetesimal disk. An outstanding problem with the existing migration models is that they invariably predict excessively large resonant populations, while observations show that the non-resonant orbits are in fact common (e.g., the main belt population is 2-4 times larger than Plutinos in the 3:2 resonance). Here we show that this problem can resolved if it is assumed that Neptune's migration was grainy, as expected from scattering encounters of Neptune with massive planetesimals. The grainy migration acts to destabilize resonant bodies with large libration amplitudes, a fraction of which ends up on stable non-resonant orbits. Thus, the non-resonant--to--resonant ratio obtained with the grainy migration is higher, up to ~10 times higher for the range of parameters investigated here, than in a model with smooth migration. The grainy migration leads to a narrower distribution of the libration amplitudes in the 3:2 resonance. The best fit to observations is obtained when it is assumed that the outer planetesimal disk below 30 AU contained 1000-4000 Plutos, or ~1000 bodies twice as massive as Pluto. We find that the probability for an outer disk object to end up on a stable orbit in the Kuiper belt is ~10^{-3}. Together, these results are consistent with having only two (known) Pluto-mass bodies in the Kuiper belt today (Pluto and Eris). The combined mass of Pluto-class objects in the original disk was ~2-8 Earth masses (M_Earth), which represents 10-40% of the estimated disk mass (M_disk ~ 20 M_Earth).

We study the effects of planetary late migration on the gas giants obliquities. We consider the planetary instability models from Nesvorny & Morbidelli (2012), in which the obliquities of Jupiter and Saturn can be excited when the spin-orbit resonances occur. The most notable resonances occur when the $s_7$ and $s_8$ frequencies, changing as a result of planetary migration, become commensurate with the precession frequencies of Jupiter's and Saturn's spin vectors. We show that Jupiter may have obtained its present obliquity by crossing of the $s_8$ resonance. This would set strict constrains on the character of migration during the early stage. Additional effects on Jupiter's obliquity are expected during the last gasp of migration when the $s_7$ resonance was approached. The magnitude of these effects depends on the precise value of the Jupiter's precession constant. Saturn's large obliquity was likely excited by capture into the $s_8$ resonance. This probably happened during the late stage of planetary migration when the evolution of the $s_8$ frequency was very slow, and the conditions for capture into the spin-orbit resonance with $s_8$ were satisfied. However, whether or not Saturn is in the spin-orbit resonance with $s_8$ at the present time is not clear, because the existing observations of Saturn's spin precession and internal structure models have significant uncertainties.

The Yarkovsky effect describes a small but significant force that affects the orbital motion of meteoroids and asteroids smaller than $30-40$ kilometers in diameter. It is caused by sunlight; when these bodies heat up in the Sun, they eventually re-radiate the energy away in the thermal waveband, which in turn creates a tiny thrust. This recoil acceleration is much weaker than solar and planetary gravitational forces, but it can produce measurable orbital changes over decades and substantial orbital effects over millions to billions of years. The same physical phenomenon also creates a thermal torque that, complemented by a torque produced by scattered sunlight, can modify the rotation rates and obliquities of small bodies as well. This rotational variant has been coined the Yarkovsky-O'Keefe-Radzievskii-Paddack (YORP) effect. During the past decade or so, the Yarkovsky and YORP effects have been used to explore and potentially resolve a number of unsolved mysteries in planetary science dealing with small bodies. Here we review the main results to date, and preview the goals for future work.

Saturn's moon Iapetus has an orbit in a transition region where the Laplace surface is bending from the equator to the orbital plane of Saturn. The orbital inclination of Iapetus to the local Laplace plane is ~8 deg, which is unexpected, because the inclination should be ~0 if Iapetus formed from a circumplanetary disk on the Laplace surface. It thus appears that some process has pumped up Iapetus's inclination while leaving its eccentricity near zero (e=0.03 at present). Here we examined the possibility that Iapetus's inclination was excited during the early solar system instability when encounters between Saturn and ice giants occurred. We found that the dynamical effects of planetary encounters on Iapetus's orbit sensitively depend on the distance of the few closest encounters. In four out of ten instability cases studied here, the orbital perturbations were too large to be plausible. In one case, Iapetus's orbit was practically unneffected. In the remaining five cases, the perturbations of Iapetus's inclination were adequate to explain its present value. In three of these cases, however, Iapetus's eccentricity was excited to 0.1-0.25, and it is not clear whether it could have been damped to its present value (0.03) by some subsequent process (e.g., tides and dynamical friction from captured irregular satellites do not seem to be strong enough). Our results therefore imply that only 2 out of 10 instability cases (~20%) can excite Iapetus's inclination to its present value (~30% of trials lead to >5 deg) while leaving its orbital eccentricity low.

While not detected yet, pairs of exoplanets in the 1:1 mean motion resonance probably exist. Low eccentricity, near-planar orbits, which in the comoving frame follow the horseshoe trajectories, are one of the possible stable configurations. Here we study transit timing variations produced by mutual gravitational interaction of planets in this orbital architecture, with the goal to develop methods that can be used to recognize this case in observational data. In particular, we use a semi-analytic model to derive parametric constraints that should facilitate data analysis. We show that characteristic traits of the transit timing variations can directly constrain the (i) ratio of planetary masses, and (ii) their total mass (divided by that of the central star) as a function of the minimum angular separation as seen from the star. In an ideal case, when transits of both planets are observed and well characterized, the minimum angular separation can also be inferred from the data. As a result, parameters derived from the observed transit timing series alone can directly provide both planetary masses scaled to the central star mass.

We develop an analytic model for transit timing variations produced by orbital conjunctions between gravitationally interacting planets. If the planetary orbits have tight orbital spacing, which is a common case among the Kepler planets, the effect of a single conjunction can be best described as: (1) a step-like change of the transit timing ephemeris with subsequent transits of the inner planet being delayed and those of the outer planet being sped up, and (2) a discrete change in sampling of the underlying oscillations from eccentricity-related interaction terms. In the limit of small orbital eccentricities, our analytic model gives explicit equations for these effects as a function of the mass and orbital separation of planets. We point out that a detection of the conjunction effect in real data is of crucial importance for the physical characterization of planetary systems from transit timing variations.

The Nice model of the dynamical instability and migration of the giant planets can explain many properties of the present Solar System, and can be used to constrain its early architecture. In the jumping-Jupiter version of the Nice model, required from the terrestrial planet constraint and dynamical structure of the asteroid belt, Jupiter has encounters with an ice giant. Here we study the survival of the Galilean satellites in the jumping-Jupiter model. This is an important concern because the ice-giant encounters, if deep enough, could dynamically perturb the orbits of the Galilean satellites, and lead to implausible results. We performed numerical integrations where we tracked the effect of planetary encounters on the Galilean moons. We considered three instability cases from Nesvorny & Morbidelli (2012) that differed in the number and distribution of encounters. We found that in one case, where the number of close encounters was relatively small, the Galilean satellite orbits were not significantly affected. In the other two, the orbital eccentricities of all moons were excited by encounters, Callisto's semimajor axis changed, and, in a large fraction of trials, the Laplace resonance of the inner three moons was disrupted. The subsequent evolution by tides damps eccentricities and can recapture the moons in the Laplace resonance. A more important constraint is represented by the orbital inclinations of the moons, which can be excited during the encounters and not appreciably damped by tides. We find that one instability case taken from Nesvorny & Morbidelli (2012) clearly fails this constraint. This shows how the regular satellites of Jupiter can be used to set limits on the properties of encounters in the jumping-Jupiter model, and help us to better understand how the early Solar System evolved.

The target asteroid of the OSIRIS-REx asteroid sample return mission, (101955) Bennu (formerly 1999 RQ$_{36}$), is a half-kilometer near-Earth asteroid with an extraordinarily well constrained orbit. An extensive data set of optical astrometry from 1999--2013 and high-quality radar delay measurements to Bennu in 1999, 2005, and 2011 reveal the action of the Yarkovsky effect, with a mean semimajor axis drift rate $da/dt = (-19.0 \pm 0.1)\times 10^{-4}$ au/Myr or $284\pm 1.5\;\rm{m/yr}$. The accuracy of this result depends critically on the fidelity of the observational and dynamical model. As an example, neglecting the relativistic perturbations of the Earth during close approaches affects the orbit with $3\sigma$ significance in $da/dt$. The orbital deviations from purely gravitational dynamics allow us to deduce the acceleration of the Yarkovsky effect, while the known physical characterization of Bennu allows us to independently model the force due to thermal emissions. The combination of these two analyses yields a bulk density of $\rho = 1260\pm70\,\rm{kg/m^3}$, which indicates a macroporosity in the range $40\pm10$% for the bulk densities of likely analog meteorites, suggesting a rubble-pile internal structure. The associated mass estimate is $(7.8\pm0.9)\times 10^{10}\, \rm{kg}$ and $GM = 5.2\pm0.6\,\rm{m^3/s^2}$. Bennu's Earth close approaches are deterministic over the interval 1654--2135, beyond which the predictions are statistical in nature. In particular, the 2135 close approach is likely within the lunar distance and leads to strong scattering and therefore numerous potential impacts in subsequent years, from 2175--2196. The highest individual impact probability is $9.5\times 10^{-5}$ in 2196, and the cumulative impact probability is $3.7\times 10^{-4}$, leading to a cumulative Palermo Scale of -1.70.

The irregular satellites of outer planets are thought to have been captured from heliocentric orbits. The exact nature of the capture process, however, remains uncertain. We examine the possibility that irregular satellites were captured from the planetesimal disk during the early Solar System instability when encounters between the outer planets occurred (Nesvorny, Vokrouhlicky & Morbidelli 2007, AJ 133; hereafter NVM07). NVM07 already showed that the irregular satellites of Saturn, Uranus and Neptune were plausibly captured during planetary encounters. Here we find that the current instability models present favorable conditions for capture of irregular satellites at Jupiter as well, mainly because Jupiter undergoes a phase of close encounters with an ice giant. We show that the orbital distribution of bodies captured during planetary encounters provides a good match to the observed distribution of irregular satellites at Jupiter. The capture efficiency for each particle in the original transplanetary disk is found to be (1.3-3.6)x10^-8. This is roughly enough to explain the observed population of jovian irregular moons. We also confirm NVM07's results for the irregular satellites of Saturn, Uranus and Neptune.

A compact formula for the stress tensor inside a self-gravitating, triaxial ellipsoid in an arbitrary rotation state is given. It contains no singularity in the incompressible medium limit. The stress tensor and the quality factor model are used to derive a solution for the energy dissipation resulting in the damping (short axis mode) or excitation (long axis) of wobbling. In the limit of an ellipsoid of revolution, we compare our solution with earlier ones and show that, with appropriate corrections, the differences in damping times estimates are much smaller than it has been claimed. This version implements corrections of misprints found in the MNRAS published text.

Jupiter Trojans are thought to be survivors of a much larger population of planetesimals that existed in the planetary region when planets formed. They can provide important constraints on the mass and properties of the planetesimal disk, and its dispersal during planet migration. Here we tested a possibility that the Trojans were captured during the early dynamical instability among the outer planets (aka the Nice model), when the semimajor axis of Jupiter was changing as a result of scattering encounters with an ice giant. The capture occurs in this model when Jupiter's orbit and its Lagrange points become radially displaced in a scattering event and fall into a region populated by planetesimals (that previously evolved from their natal transplanetary disk to ~5 AU during the instability). Our numerical simulations of the new capture model, hereafter jump capture, satisfactorily reproduce the orbital distribution of the Trojans and their total mass. The jump capture is potentially capable of explaining the observed asymmetry in the number of leading and trailing Trojans. We find that the capture probability is (6-8)x10^-7 for each particle in the original transplanetary disk, implying that the disk contained (3-4)x10^7 planetesimals with absolute magnitude H<9 (corresponding to diameter D=80 km for a 7% albedo). The disk mass inferred from this work, M_disk=14-28 M_Earth, is consistent with the mass deduced from recent dynamical simulations of the planetary instability.

We seek evidence of the Yarkovsky effect among Near Earth Asteroids (NEAs) by measuring the Yarkovsky-related orbital drift from the orbital fit. To prevent the occurrence of unreliable detections we employ a high precision dynamical model, including the Newtonian attraction of 16 massive asteroids and the planetary relativistic terms, and a suitable astrometric data treatment. We find 21 NEAs whose orbital fits show a measurable orbital drift with a signal to noise ratio (SNR) greater than 3. The best determination is for asteroid (101955) 1999 RQ36, resulting in the recovery of one radar apparition and an orbit improvement by two orders of magnitude. In addition, we find 16 cases with a lower SNR that, despite being less reliable, are good candidates for becoming stronger detections in the future. In some cases it is possible to constrain physical quantities otherwise unknown by means of the detected orbital drift. Furthermore, the distribution of the detected orbital drifts shows an excess of retrograde rotators that can be connected to the delivery mechanism from the most important NEA feeding resonances and allows us to infer the distribution for NEAs obliquity. We discuss the implications of the Yarkovsky effect for impact predictions. In particular, for asteroid (29075) 1950 DA our results favor a retrograde rotation that would rule out an impact in 2880.

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

Main belt asteroids (6070) Rheinland and (54827) 2001NQ8 belong to a small population of couples of bodies which reside on very similar heliocentric orbits. Vokrouhlicky & Nesvorny (2008, AJ 136, 280) promoted a term "asteroid pairs", pointing out their common origin within the past tens to hundreds of ky. Previous attempts to reconstruct the initial configuration of Rheinland and 2001NQ8 at the time of their separation have led to the prediction that Rheinland's rotation should be retrograde. Here we report extensive photometric observations of this asteroid and use the lightcurve inversion technique to directly determine its rotation state and shape. We confirm the retrograde sense of rotation of Rheinland, with obliquity value constrained to be >= 140 deg. The ecliptic longitude of the pole position is not well constrained as yet. The asymmetric behavior of Rheinland's lightcurve reflects a sharp, near-planar edge in our convex shape representation of this asteroid. Our calibrated observations in the red filter also allow us to determine $H_R = 13.68\pm 0.05$ and $G = 0.31\pm 0.05$ values of the H-G system. With the characteristic color index $V-R = 0.49\pm 0.05$ for the S-type asteroids, we thus obtain $H = 14.17\pm 0.07$ for the absolute magnitude of (6070) Rheinland. This a significantly larger value than previously obtained from analysis of the astrometric survey observations. We next use the obliquity constraint for Rheinland to eliminate some degree of uncertainty in the past propagation of its orbit. This is because the sign of the past secular change of its semimajor axis due to the Yarkovsky effect is now constrained. Determination of the rotation state of the secondary component, asteroid (54827) 2001NQ8, is the key element in further constraining the age of the pair and its formation process.

The solar system is dusty, and would become dustier over time as asteroids collide and comets disintegrate, except that small debris particles in interplanetary space do not last long. They can be ejected from the solar system by Jupiter, thermally destroyed near the Sun, or physically disrupted by collisions. Also, some are swept by the Earth (and other planets), producing meteors. Here we develop a dynamical model for the solar system meteoroids and use it to explain meteor radar observations. We find that the Jupiter Family Comets (JFCs) are the main source of the prominent concentrations of meteors arriving to the Earth from the helion and antihelion directions. To match the radiant and orbit distributions, as measured by the Canadian Meteor Orbit Radar (CMOR) and Advanced Meteor Orbit Radar (AMOR), our model implies that comets, and JFCs in particular, must frequently disintegrate when reaching orbits with low perihelion distance. Also, the collisional lifetimes of millimeter particles may be longer (>10^5 yr at 1 AU) than postulated in the standard collisional models (10^4 yr at 1 AU), perhaps because these chondrule-sized meteoroids are stronger than thought before. Using observations of the Infrared Astronomical Satellite (IRAS) to calibrate the model, we find that the total cross section and mass of small meteoroids in the inner solar system are (1.7-3.5)x10^11 km^2 and 4x10^19 g, respectively, in a good agreement with previous studies. The mass input required to keep the Zodiacal Cloud (ZC) in a steady state is estimated to be 10^4-10^5 kg/s. The input is up to 10 times larger than found previously, mainly because particles released closer to the Sun have shorter collisional lifetimes, and need to be supplied at a faster rate.

The Oort Cloud Comets (OCCs), exemplified by the Great Comet of 1997 (Hale-Bopp), are occasional visitors from the heatless periphery of the solar system. Previous works hypothesized that a great majority of OCCs must physically disrupt after one or two passages through the inner solar system, where strong thermal gradients can cause phase transitions or volatile pressure buildup. Here we study the fate of small debris particles produced by OCC disruptions to determine whether the imprints of a hypothetical population of OCC meteoroids can be found in the existing meteor radar data. We find that OCC particles with diameters D<10 um are blown out from the solar system by radiation pressure, while those with D>1 um have a very low Earth-impact probability. The intermediate particle sizes, D=100 um, represent a sweet spot. About 1% of these particles orbitally evolve by Poynting-Robertson drag to reach orbits with semimajor axis a=1 AU. They are expected to produce meteors with radiants near the apex of the Earth's orbital motion. We find that the model distributions of their impact speeds and orbits provide a good match to radar observations of apex meteors, except for the eccentricity distribution, which is more skewed toward e=1 in our model. Finally, we propose an explanation for the long-standing problem in meteor science related to the relative strength of apex and helion/antihelion sources. As we show in detail, the observed trend, with the apex meteors being more prominent in observations of highly sensitive radars, can be related to orbital dynamics of particles released on the long-period orbits.

We investigate the orbital evolution of a system of N mutually interacting stars on initially circular orbits around the dominating central mass. We include perturbative influence of a distant axisymmetric source and an extended spherical potential. In particular, we focus on the case when the secular evolution of orbital eccentricities is suppressed by the spherical perturbation. By means of standard perturbation methods, we derive semi-analytic formulae for the evolution of normal vectors of the individual orbits. We find its two qualitatively different modes. Either the orbits interact strongly and, under such circumstances, they become dynamically coupled, precessing synchronously in the potential of the axisymmetric perturbation. Or, if their mutual interaction is weaker, the orbits precess independently, interchanging periodically their angular momentum, which leads to oscillations of inclinations. We argue that these processes may have been fundamental for the evolution of the disc of young stars orbiting the supermassive black hole in the centre of the Milky Way.

The size distribution in the cold classical Kuiper belt can be approximated by two idealized power laws: one with steep slope for radii R>R* and one with shallow slope for R<R*, where R*~25-50 km. Previous works suggested that the SFD roll-over at R* can be the result of extensive collisional grinding in the Kuiper belt that led to the catastrophic disruption of most bodies with R<R*. Here we use a new code to test the effect of collisions in the Kuiper belt. We find that the observed roll-over could indeed be explained by collisional grinding provided that the initial mass in large bodies was much larger than the one in the present Kuiper belt, and was dynamically depleted. In addition to the size distribution changes, our code also tracks the effects of collisions on binary systems. We find that it is generally easier to dissolve wide binary systems, such as the ones existing in the cold Kuiper belt today, than to catastrophically disrupt objects with R~R*. Thus, the binary survival sets important limits on the extent of collisional grinding in the Kuiper belt. We find that the extensive collisional grinding required to produce the SFD roll-over at R* would imply a strong gradient of the binary fraction with R and separation, because it is generally easier to dissolve binaries with small components and/or those with wide orbits. The expected binary fraction for R<R* is <0.1. The present observational data do not show such a gradient. Instead, they suggest a large binary fraction of ~0.4 for R=30-40 km. This may indicate that the roll-over was not produced by disruptive collisions, but is instead a fossil remnant of the KBO formation process.

Asteroid pairs sharing similar heliocentric orbits were found recently. Backward integrations of their orbits indicated that they separated gently with low relative velocities, but did not provide additional insight into their formation mechanism. A previously hypothesized rotational fission process4 may explain their formation - critical predictions are that the mass ratios are less than about 0.2 and, as the mass ratio approaches this upper limit, the spin period of the larger body becomes long. Here we report photometric observations of a sample of asteroid pairs revealing that primaries of pairs with mass ratios much less than 0.2 rotate rapidly, near their critical fission frequency. As the mass ratio approaches 0.2, the primary period grows long. This occurs as the total energy of the system approaches zero requiring the asteroid pair to extract an increasing fraction of energy from the primary's spin in order to escape. We do not find asteroid pairs with mass ratios larger than 0.2. Rotationally fissioned systems beyond this limit have insufficient energy to disrupt. We conclude that asteroid pairs are formed by the rotational fission of a parent asteroid into a proto-binary system which subsequently disrupts under its own internal system dynamics soon after formation.

The influence of optical scattering and thermal radiation models an the Yarkovsky-O'Keefe-Radzievskii-Paddack (YORP) effect is studied. Lambertian formulation is compared with Hapke scattering and emission laws and Lommel-Seeliger reflection. Although the form of reflectivity function strongly influences mean torques due to scattering or thermal radiation alone, their combined contribution to the rotation period YORP is not much different from the standard Lambertian values. For higher albedo values the differences between the Hapke and Lambert models become significant for the YORP in attitude.

Processes such as the solar wind sputtering and micrometeorite impacts can modify optical properties of surfaces of airless bodies. This explains why spectra of the main belt asteroids, exposed to these `space weathering' processes over eons, do not match the laboratory spectra of ordinary chondrite (OC) meteorites. In contrast, an important fraction of Near Earth Asteroids (NEAs), defined as Q-types in the asteroid taxonomy, display spectral attributes that are a good match to OCs. Here we study the possibility that the Q-type NEAs underwent recent encounters with the terrestrial planets and that the tidal gravity (or other effects) during these encounters exposed fresh OC material on the surface (thus giving it the Q-type spectral properties). We used numerical integrations to determine the statistics of encounters of NEAs to planets. The results were used to calculate the fraction and orbital distribution of Q-type asteroids expected in the model as a function of the space weathering timescale, t_sw (see main text for definition), and maximum distance, r*, at which planetary encounters can reset the surface. We found that t_sw ~ 1e6 yr (at 1 AU) and r* ~ 5 R_pl, where R_pl is the planetary radius, best fit the data. Values t_sw < 1e5 yr would require that r* > 20 R_pl, which is probably implausible because these very distant encounters should be irrelevant. Also, the fraction of Q-type NEAs would be probably much larger than the one observed if t_sw > 1e7 yr. We found that t_sw ~ q^2, where q is the perihelion distance, expected if the solar wind sputtering controls t_sw, provides a better match to the orbital distribution of Q-type NEAs than models with fixed t_sw. We also discuss how the Earth magnetosphere and radiation effects such as YORP can influence the spectral properties of NEAs.

We have searched Microvariability and Oscillations of STars (MOST) satellite photometry obtained in 2004, 2005, and 2007 of the solar-type star HD 209458 for Trojan asteroid swarms dynamically coupled with the system's transiting "hot Jupiter" HD 209458b. Observations of the presence and nature of asteroids around other stars would provide unique constraints on migration models of exoplanetary systems. Our results set an upper limit on the optical depth of Trojans in the HD 209458 system that can be used to guide current and future searches of similar systems by upcoming missions. Using cross-correlation methods with artificial signals implanted in the data, we find that our detection limit corresponds to a relative Trojan transit depth of 1\times10-4, equivalent to ~1 lunar mass of asteroids, assuming power-law Trojan size distributions similar to Jupiter's Trojans in our solar system. We confirm with dynamical interpretations that some asteroids could have migrated inward with the planet to its current orbit at 0.045 AU, and that the Yarkovsky effect is ineffective at eliminating objects of > 1 m in size. However, using numerical models of collisional evolution we find that, due to high relative speeds in this confined Trojan environment, collisions destroy the vast majority of the asteroids in <10 Myr. Our modeling indicates that the best candidates to search for exoTrojan swarms in 1:1 mean resonance orbits with "hot Jupiters" are young systems (ages of about 1 Myr or less). Years of Kepler satellite monitoring of such a system could detect an asteroid swarm with a predicted transit depth of 3\times10-7.

The multiple system V505 Sagittarii is composed of at least three stars: a compact eclipsing pair and a distant component, which orbit is measured directly using speckle interferometry. In order to explain the observed orbit of the third body in V505 Sagittarii and also other observable quantities, namely the minima timings of the eclipsing binary and two different radial velocities in the spectrum, we thoroughly test a fourth-body hypothesis - a perturbation by a dim, yet-unobserved object. We use an N-body numerical integrator to simulate future and past orbital evolution of 3 or 4 components in this system. We construct a suitable chi^2 metric from all available speckle-interferometry, minima-timings and radial-velocity data and we scan a part of a parameter space to get at least some of allowed solutions. In principle, we are able to explain all observable quantities by a presence of a fourth body, but the resulting likelihood of this hypothesis is very low. We also discuss other theoretical explanations of the minima timings variations. Further observations of the minima timings during the next decade or high-resolution spectroscopic data can significantly constrain the model.

The radioscience experiment is one of the on board experiment of the Mercury ESA mission BepiColombo that will be launched in 2014. The goals of the experiment are to determine the gravity field of Mercury and its rotation state, to determine the orbit of Mercury, to constrain the possible theories of gravitation (for example by determining the post-Newtonian (PN) parameters), to provide the spacecraft position for geodesy experiments and to contribute to planetary ephemerides improvement. This is possible thanks to a new technology which allows to reach great accuracies in the observables range and range rate; it is well known that a similar level of accuracy requires studying a suitable model taking into account numerous relativistic effects. In this paper we deal with the modelling of the space-time coordinate transformations needed for the light-time computations and the numerical methods adopted to avoid rounding-off errors in such computations.

The zodiacal cloud is a thick circumsolar disk of small debris particles produced by asteroid collisions and comets. Here, we present a zodiacal cloud model based on the orbital properties and lifetimes of comets and asteroids, and on the dynamical evolution of dust after ejection. The model is quantitatively constrained by IRAS observations of thermal emission, but also qualitatively consistent with other zodiacal cloud observations. We find that 85-95% of the observed mid-infrared emission is produced by particles from the Jupiter-family comets (JFCs) and $<$10% by dust from long period comets. Asteroidal dust is found to be present at $<$10%. We suggest that spontaneous disruptions of JFCs, rather than the usual cometary activity driven by sublimating volatiles, is the main mechanism that librates cometary particles into the zodiacal cloud. Our results imply that JFC particles represent $\sim$85% of the total mass influx at Earth. Since their atmospheric entry speeds are typically low ($\approx$14.5 km s$^{-1}$ mean for D=100-200 $\mu$m with $\approx$12 km s$^{-1}$ being the most common case), many JFC grains should survive frictional heating and land on the Earth's surface. This explains why most micrometeorites collected in antarctic ice have primitive carbonaceous composition. The present mass of the inner zodiacal cloud at $<$5 AU is estimated to be 1-$2\times10^{19}$ g, mainly in D=100-200 $\mu$m particles. The inner zodiacal cloud should have been $>10^4$ times brighter during the Late Heavy Bombardment (LHB) epoch $\approx$3.8 Gyr ago, when the outer planets scattered numerous comets into the inner solar system. The bright debris disks with a large 24-$\mu$m excess observed around mature stars may be an indication of massive cometary populations existing in those systems.

The recent discovery of the 2003EL61 collisional family in the Kuiper belt (Brown et al., 2007) is surprising because the formation of such a family is a highly improbable event in today's belt. Assuming Brown et al.'s estimate of the size of progenitors, we find that the probability that a Kuiper belt object was involved in such a collision since primordial times is less than roughly 0.001. In addition, it is not possible for the collision to have occurred in a massive primordial Kuiper belt because the dynamical coherence of the family would not have survived whatever event produced the currently observed orbital excitation. Here we suggest that the family is the result of a collision between two scattered disk objects. We show that the probability that a collision occurred between two such objects with sizes similar to those advocated in Brown et al. (2007) and that the center of mass of the resulting family is on an orbit typical of the Kuiper belt can be as large as 47%. Given the large uncertainties involved in this estimate, this result is consistent with the existence of one such family. If true, this result has implications far beyond the origin of a single collisional family, because it shows that collisions played an important role in shaping the dynamical structure of the small body populations that we see today.

This paper discusses the IAU Resolutions B1.3, B1.4, B1.5 and B1.9 (2000) that were adopted during the 24th General Assembly in Manchester, 2000 and provides details and explanations for these Resolutions. It is explained why they present significant progress over the corresponding IAU 1991 Resolutions and why they are necessary in the light of present accuracies in astrometry, celestial mechanics and metrology. In fact most of these Resolutions are consistent with astronomical models and software already in use. The metric tensors and gravitational potentials of both the Barycentric Celestial Reference System and Geocentric Celestial Reference System are defined and discussed. The necessity and relevance of the two celestial reference systems are explained. The transformations of coordinates and gravitational potentials are discussed. Potential coefficients parameterizing the post-Newtonian gravitational potentials are expounded. Simplified versions of the time transformations suitable for modern clock accuracies are elucidated. Various approximations used in the Resolutions are explicated and justified. Some models (e.g. for higher spin moments) that serve the purpose for estimating orders of magnitude have actually never been published before.

The long-term evolution of stellar orbits bound to a massive centre is studied in order to understand the cores of star clusters in central regions of galaxies. Stellar trajectories undergo tiny perturbation, the origin of which is twofold: (i) gravitational field of a thin gaseous disc surrounding the galactic centre, and (ii) cumulative drag due to successive interactions of the stars with material of the disc. Both effects are closely related because they depend on the total mass of the disc, assumed to be a small fraction of the central mass. It is shown that, in contrast to previous works, most of the retrograde (with respect to the disc) orbits are captured by the central object, presumably a massive black hole. Initially prograde orbits are also affected, so that statistical properties of the central star cluster in quasi-equilibrium may differ significantly from those deduced in previous analyses.

Long-term evolution of a stellar orbit captured by a massive galactic center via successive interactions with an accretion disc has been examined. An analytical solution describing evolution of the stellar orbital parameters during the initial stage of the capture was found. Our results are applicable to thin Keplerian discs with an arbitrary radial distribution of density and rather general prescription for the star-disc interaction. Temporal evolution is given in the form of quadrature which can be carried out numerically.

We have constructed general relativistic models of a stationary, axially symmnetric, Keplerian thin disk around a rotating black hole. We computed profiles of a spectral line, emitted in the inner region of the disk. In our models we have taken into account also the self-gravity of the disk. The aim of this work is to study gravitational effects on the line profiles in connection with the X-ray features observed in spectra of active galactic nuclei. In some cases, the calculated profiles are clearly affected by the disk gravity but relativistic dragging effects are found to be negligible.

This preprint is based on two articles accepted by MNRAS and ApJ. Our aim is to study the evolution of the orbit of a star under the influence of interactions with an accretion disc in an AGN. The model considered consists of a low-mass compact object orbiting a supermassive black hole and colliding periodically with the accretion disc. Gravitational field of the nucleus is described by the Kerr metric. The star is assumed to move along a geodesic arc between successive interactions with an equatorial accretion disc. We study the gravitomagnetic (Lense-Thirring) precession of the trajectory. We do not assume any particular value for the eccentricity or inclination of the orbit or the angular momentum of the black hole. We also discuss the periodicity related to the relativistic shift of the pericenter. The model of star-disc interactions has been suggested to explain the X-ray variability observed in the Seyfert galaxy NGC~6814. We briefly discuss on this subject.