Source author record

David Vokrouhlicky

David Vokrouhlicky appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

23works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

23 published item(s)

preprint2022arXiv

Asteroid Families: properties, recent advances and future opportunities

Collisions are one of the key processes shaping planetary systems. Asteroid families are outcomes of such collisions still identifiable across our solar system. The families provide a unique view of catastrophic disruption phenomena and have been in the focus of planetary scientists for more than a century. Most of them are located in the main belt, a ring of asteroids between Mars and Jupiter. Here we review the basic properties of the families, discuss some recent advances, and anticipate future challenges. This review pays more attention to dynamic aspects such as family identification, age determination, and long-term evolution. The text, however, goes beyond that. Especially, we cover the details of young families that see the major advances in the last years, and we anticipate it will develop even faster in the future. We also discuss the relevance of asteroid families for water-ice content in the asteroid belt and our current knowledge on links between families and main-belt comets.

preprint2022arXiv

Dynamical Implantation of Blue Binaries in the Cold Classical Kuiper Belt

Colors and binarity provide important constraints on the Kuiper belt formation. The cold classical objects at radial distance r=42-47 au from the Sun are predominantly very red (spectral slope s>17%) and often exist as equal-size binaries (~30% observed binary fraction). This has been taken as evidence for the in-situ formation of cold classicals. Interestingly, a small fraction (~10%) of cold classicals is less red with s<17%, and these "blue" bodies are often found in wide binaries. Here we study the dynamical implantation of blue binaries from r<42 au. We find that they can be implanted into the cold classical belt from a wide range of initial radial distances, but the survival of the widest blue binaries -- 2001 QW322 and 2003 UN284 -- implies formation at r>30 au. This would be consistent with the hypothesized less-red to very-red transition at 30<r<40 au. For any reasonable choice of parameters (Neptune's migration history, initial disk profile, etc.), however, our model predicts a predominance of blue singles, rather than blue binaries, which contradicts existing observations. We suggest that wide blue binaries formed in situ at r=42-47 au and their color reflects early formation in a protoplanetary gas disk. The predominantly VR colors of cold classicals may be related to the production of methanol and other hydrocarbons during the late stages of the disk, when the temperature at 45 au dropped to 20 K and carbon monoxide was hydrogenated.

preprint2020arXiv

Clarissa Family Age from the Yarkovsky Effect Chronology

The Clarissa family is a small collisional family composed of primitive C-type asteroids. It is located in a dynamically stable zone of the inner asteroid belt. In this work we determine the formation age of the Clarissa family by modeling planetary perturbations as well as thermal drift of family members due to the Yarkovsky effect. Simulations were carried out using the Swift-rmvs4 integrator modified to account for the Yarkovsky and Yarkovsky-O'Keefe-Radzievskii-Paddack (YORP) effects. We ran multiple simulations starting with different ejection velocity fields of fragments, varying proportion of initially retrograde spins, and also tested different Yarkovsky/YORP models. Our goal was to match the observed orbital structure of the Clarissa family which is notably asymmetrical in the proper semimajor axis. The best fits were obtained with the initial ejection velocities < ~20 m/s of diameter D=2 km fragments, 4:1 preference for spin-up by YORP, and assuming that 80% of small family members initially had retrograde rotation. The age of the Clarissa family was found to be 56+/-6 Myr for the assumed asteroid density 1.5 g/cm3. Small variation of density to smaller or larger value would lead to slightly younger or older age estimates. This is the first case where the Yarkovsky effect chronology has been successfully applied to an asteroid family younger than 100 Myr.

preprint2020arXiv

OSSOS XX: The Meaning of Kuiper Belt Colors

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.

preprint2016arXiv

Neptune's Orbital Migration Was Grainy, Not Smooth

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

preprint2016arXiv

The orbital distribution of trans-Neptunian objects beyond 50 au

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.

preprint2015arXiv

Tilting Jupiter (a bit) and Saturn (a lot) During Planetary Migration

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.

preprint2014arXiv

Excitation of the Orbital Inclination of Iapetus during Planetary Encounters

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.

preprint2014arXiv

Orbit and Bulk Density of the OSIRIS-REx Target Asteroid (101955) Bennu

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σ$ 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 $ρ= 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.

preprint2014arXiv

The Effect of Conjunctions on the Transit Timing Variations of Exoplanets

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.

preprint2014arXiv

Transit timing variations for planets coorbiting in the horseshoe regime

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.

preprint2013arXiv

Capture of Trojans by Jumping Jupiter

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.

preprint2013arXiv

Stress field and spin axis relaxation for inelastic triaxial ellipsoids

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.

preprint2011arXiv

Dynamical Model for the Zodiacal Cloud and Sporadic Meteors

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.

preprint2011arXiv

Dynamics of Dust Particles Released from Oort Cloud Comets and Their Contribution to Radar Meteors

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.

preprint2011arXiv

Observed Binary Fraction Sets Limits on the Extent of Collisional Grinding in the Kuiper Belt

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.

preprint2011arXiv

Secular theory of the orbital evolution of the young stellar disc in the Galactic Centre

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.

preprint2011arXiv

Spin vector and shape of (6070) Rheinland and their implications

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.

preprint2010arXiv

Do Planetary Encounters Reset Surfaces of Near Earth Asteroids?

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.

preprint2010arXiv

Searching for Trojan Asteroids in the HD 209458 System: Space-based MOST Photometry and Dynamical Modeling

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.

preprint2010arXiv

YORP effect with anisotropic radiation

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.

preprint2009arXiv

Cometary Origin of the Zodiacal Cloud and Carbonaceous Micrometeorites

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 $μ$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 $μ$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-$μ$m excess observed around mature stars may be an indication of massive cometary populations existing in those systems.

preprint1996arXiv

Testing for gravitationally preferred directions using the lunar orbit

As gravity is a long-range force, it is {\it a priori} conceivable that the Universe's global matter distribution select a preferred rest frame for local gravitational physics. At the post-Newtonian approximation, the phenomenology of preferred-frame effects is described by two parameters, $α_1$ and $α_2$, the second of which is already very tightly constrained. Confirming previous suggestions, we show through a detailed Hill-Brown type calculation of a perturbed lunar orbit that lunar laser ranging data have the potential of constraining $α_1$ at the $10^{-4}$ level. It is found that certain retrograde planar orbits exhibit a resonant sensitivity to external perturbations linked to a fixed direction in space. The lunar orbit being quite far from such a resonance exhibits no significant enhancement due to solar tides. Our Hill-Brown analysis is extended to the perturbation linked to a possible differential acceleration toward the galactic center. It is, however, argued that there are strong {\it a priori} theoretical constraints on the conceivable magnitude of such an effect.