BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Renu Malhotra

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

ResearcherUnclaimed source record
astro-ph.EPastro-phastro-ph.SRastro-ph.IMastro-ph.GAmath.DS

Catalog footprint

What is connected

35works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

35 published item(s)

preprint2022arXiv

Pluto near the edge of chaos

Many of the unusual properties of Pluto's orbit are widely accepted as evidence for the orbital migration of the giant planets in early solar system history. However, some properties remain an enigma. Pluto's long term orbital stability is supported by two special properties of its orbit that limit the location of its perihelion in azimuth and in latitude. We revisit Pluto's orbital dynamics with a view to elucidating the individual and collective gravitational effects of the giant planets on constraining its perihelion location. While the resonant perturbations from Neptune account for the azimuthal constraint on Pluto's perihelion location, we demonstrate that the long term and steady persistence of the latitudinal constraint is possible only in a narrow range of additional secular forcing which arises fortuitously from the particular orbital architecture of the other giant planets. Our investigations also find that Jupiter has a largely stabilizing influence whereas Uranus has a largely destabilizing influence on Pluto's orbit. Overall, Pluto's orbit is rather surprisingly close to a zone of strong chaos.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

An Integrative Analysis of the HD 219134 Planetary System and the Inner Solar System: Extending DYNAMITE with Enhanced Orbital Dynamical Stability Criteria

Planetary architectures remain unexplored for the vast majority of exoplanetary systems, even among the closest ones, with potentially hundreds of planets still ``hidden" from our knowledge. DYNAMITE is a powerful software package that can predict the presence and properties of these yet undiscovered planets. We have significantly expanded the integrative capabilities of DYNAMITE, which now allows for (i) planets of unknown inclinations alongside planets of known inclinations, (ii) population statistics and model distributions for the eccentricity of planetary orbits, and (iii) three different dynamical stability criteria. We demonstrate the new capabilities with a study of the HD 219134 exoplanet system consisting of four confirmed planets and two likely candidates, where five of the likely planets are Neptune-size or below with orbital periods less than 100 days. By integrating the known data for the HD 219134 planetary system with contextual and statistical exoplanet population information, we tested different system architecture hypotheses to determine their likely dynamical stability. Our results provide support for the planet candidates, and we predict at least two additional planets in this system. We also deploy DYNAMITE on analogs of the inner Solar System by excluding Venus or Earth from the input parameters to test DYNAMITE's predictive power. Our analysis finds the system remains stable while also recovering the excluded planets, demonstrating the increasing capability of DYNAMITE to accurately and precisely model the parameters of additional planets in multi-planet systems.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

New results on orbital resonances

Perturbative analyses of planetary resonances commonly predict singularities and/or divergences of resonance widths at very low and very high eccentricities. We have recently re-examined the nature of these divergences using non-perturbative numerical analyses, making use of Poincaré sections but from a different perspective relative to previous implementations of this method. This perspective reveals fine structure of resonances which otherwise remains hidden in conventional approaches, including analytical, semi-analytical and numerical-averaging approaches based on the critical resonant angle. At low eccentricity, first order resonances do not have diverging widths but have two asymmetric branches leading away from the nominal resonance location. A sequence of structures called ``low-eccentricity resonant bridges" connecting neighboring resonances is revealed. At planet-grazing eccentricity, the true resonance width is non-divergent. At higher eccentricities, the new results reveal hitherto unknown resonant structures and show that these parameter regions have a loss of some -- though not necessarily entire -- resonance libration zones to chaos. The chaos at high eccentricities was previously attributed to the overlap of neighboring resonances. The new results reveal the additional role of bifurcations and co-existence of phase-shifted resonance zones at higher eccentricities. By employing a geometric point of view, we relate the high eccentricity phase space structures and their transitions to the shapes of resonant orbits in the rotating frame. We outline some directions for future research to advance understanding of the dynamics of mean motion resonances.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Refining the transit timing and photometric analysis of TRAPPIST-1: Masses, radii, densities, dynamics, and ephemerides

We have collected transit times for the TRAPPIST-1 system with the Spitzer Space Telescope over four years. We add to these ground-based, HST and K2 transit time measurements, and revisit an N-body dynamical analysis of the seven-planet system using our complete set of times from which we refine the mass ratios of the planets to the star. We next carry out a photodynamical analysis of the Spitzer light curves to derive the density of the host star and the planet densities. We find that all seven planets' densities may be described with a single rocky mass-radius relation which is depleted in iron relative to Earth, with Fe 21 wt% versus 32 wt% for Earth, and otherwise Earth-like in composition. Alternatively, the planets may have an Earth-like composition, but enhanced in light elements, such as a surface water layer or a core-free structure with oxidized iron in the mantle. We measure planet masses to a precision of 3-5%, equivalent to a radial-velocity (RV) precision of 2.5 cm/sec, or two orders of magnitude more precise than current RV capabilities. We find the eccentricities of the planets are very small; the orbits are extremely coplanar; and the system is stable on 10 Myr timescales. We find evidence of infrequent timing outliers which we cannot explain with an eighth planet; we instead account for the outliers using a robust likelihood function. We forecast JWST timing observations, and speculate on possible implications of the planet densities for the formation, migration and evolution of the planet system.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Dynamical instabilities in systems of multiple short-period planets are likely driven by secular chaos: a case study of Kepler-102

We investigated the dynamical stability of high-multiplicity Kepler and K2 planetary systems. Our numerical simulations find instabilities in $\sim20\%$ of the cases on a wide range of timescales (up to $5\times10^9$ orbits) and over an unexpectedly wide range of initial dynamical spacings. To identify the triggers of long-term instability in multi-planet systems, we investigated in detail the five-planet Kepler-102 system. Despite having several near-resonant period ratios, we find that mean motion resonances are unlikely to directly cause instability for plausible planet masses in this system. Instead, we find strong evidence that slow inward transfer of angular momentum deficit (AMD) via secular chaos excites the eccentricity of the innermost planet, Kepler-102 b, eventually leading to planet-planet collisions in $\sim80\%$ of Kepler-102 simulations. Kepler-102 b likely has a mass $>\sim0.1M_{\oplus}$, hence a bulk density exceeding about half Earth's, in order to avoid dynamical instability. To investigate the role of secular chaos in our wider set of simulations, we characterize each planetary system's AMD evolution with a "spectral fraction" calculated from the power spectrum of short integrations ($\sim5\times10^6$ orbits). We find that small spectral fractions ($\lesssim0.01$) are strongly associated with dynamical stability on long timescales ($5\times10^9$ orbits) and that the median time to instability decreases with increasing spectral fraction. Our results support the hypothesis that secular chaos is the driver of instabilities in many non-resonant multi-planet systems, and also demonstrate that the spectral analysis method is an efficient numerical tool to diagnose long term (in)stability of multi-planet systems from short simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Lunar close encounters compete with the circumterrestrial Lidov-Kozai effect

Luna 3 (or Lunik 3 in Russian sources) was the first spacecraft to perform a flyby of the Moon. Launched in October 1959 on a translunar trajectory with large semi-major axis and eccentricity, it collided with the Earth in late March 1960. The short, 6-month dynamical lifetime has often been explained through an increase in eccentricity due to the Lidov-Kozai effect. However, the classical Lidov-Kozai solution is only valid in the limit of small semi-major axis ratio, a condition that is satisfied only for solar (but not for lunar) perturbations. We undertook a study of the dynamics of Luna 3 with the aim of assessing the principal mechanisms affecting its evolution. We analyze the Luna 3 trajectory by generating accurate osculating solutions, and by comparing them to integrations of singly- and doubly-averaged equations of motion in vectorial form. Lunar close encounters, which cannot be reproduced in an averaging approach, decisively affect the trajectory and break the doubly-averaged dynamics. Solar perturbations induce oscillations of intermediate period that affect the geometry of the close encounters and cause the singly-averaged and osculating inclinations to change quadrants (the orbital plane ``flips''). We find that the peculiar evolution of Luna 3 can only be explained by taking into account lunar close encounters and intermediate-period terms; such terms are averaged out in the Lidov-Kozai solution, which is not adequate to describe translunar or cislunar trajectories. Understanding the limits of the Lidov-Kozai solution is of particular significance for the motion of objects in the Earth-Moon environment and of exoplanetary systems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On the Divergence of First Order Resonance Widths at Low Eccentricities

Orbital resonances play an important role in the dynamics of planetary systems. Classical theoretical analyses found in textbooks report that libration widths of first order mean motion resonances diverge for nearly circular orbits. Here we examine the nature of this divergence with a non-perturbative analysis of a few first order resonances interior to a Jupiter-mass planet. We show that a first order resonance has two branches, the pericentric and the apocentric resonance zone. As the eccentricity approaches zero, the centers of these zones diverge away from the nominal resonance location but their widths shrink. We also report a novel finding of "bridges" between adjacent first order resonances: at low eccentricities, the apocentric libration zone of a first order resonance smoothly connects with the pericentric libration zone of the neighboring first order resonance. These bridges may facilitate resonant migration across large radial distances in planetary systems, entirely in the low eccentricity regime.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Search for L5 Earth Trojans with DECam

Most of the major planets in the Solar System support populations of co-orbiting bodies, known as Trojans, at their L4 and L5 Lagrange points. In contrast, Earth has only one known co-orbiting companion. This paper presents the results from a search for Earth Trojans using the DECam instrument on the Blanco Telescope at CTIO. This search found no additional Trojans in spite of greater coverage compared to previous surveys of the L5 point. Therefore, the main result of this work is to place the most stringent constraints to date on the population of Earth Trojans. These constraints depend on assumptions regarding the underlying population properties, especially the slope of the magnitude distribution (which in turn depends on the size and albedo distributions of the objects). For standard assumptions, we calculate upper limits to a 90% confidence limit on the L5 population of $N_{ET}<1$ for magnitude $H<15.5$, $N_{ET}=60-85$ for $H<19.7$, and $N_{ET}\ $= 97 for $H=20.4$. This latter magnitude limit corresponds to Trojans $\sim$300 m in size for albedo $0.15$. At H=19.7, these upper limits are consistent with previous L4 Earth Trojan constraints and significantly improve L5 constraints.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The evolution of dust-disk sizes from a homogeneous analysis of 1-10 Myr-old stars

We utilize ALMA archival data to estimate the dust disk size of 152 protoplanetary disks in Lupus (1-3 Myr), Chamaeleon I (2-3 Myr), and Upper-Sco (5-11 Myr). We combine our sample with 47 disks from Tau/Aur and Oph whose dust disk radii were estimated, as here, through fitting radial profile models to visibility data. We use these 199 homogeneously derived disk sizes to identify empirical disk-disk and disk-host property relations as well as to search for evolutionary trends. In agreement with previous studies, we find that dust disk sizes and millimeter luminosities are correlated, but show for the first time that the relationship is not universal between regions. We find that disks in the 2-3 Myr-old Cha I are not smaller than disks in other regions of similar age, and confirm the Barenfeld et al. (2017) finding that the 5-10 Myr USco disks are smaller than disks belonging to younger regions. Finally, we find that the outer edge of the Solar System, as defined by the Kuiper Belt, is consistent with a population of dust disk sizes which have not experienced significant truncation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Scientific Impact of the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) for Solar System Science

Vera C. Rubin Observatory will be a key facility for small body science in planetary astronomy over the next decade. It will carry out the Legacy Survey of Space and Time (LSST), observing the sky repeatedly in u, g, r, i, z, and y over the course of ten years using a 6.5 m effective diameter telescope with a 9.6 square degree field of view, reaching approximately r = 24.5 mag (5-σ depth) per visit. The resulting dataset will provide extraordinary opportunities for both discovery and characterization of large numbers (10--100 times more than currently known) of small solar system bodies, furthering studies of planetary formation and evolution. This white paper summarizes some of the expected science from the ten years of LSST, and emphasizes that the planetary astronomy community should remain invested in the path of Rubin Observatory once the LSST is complete.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Corralling a distant planet with extreme resonant Kuiper belt objects

The four longest period Kuiper belt objects have orbital periods close to integer ratios with each other. A hypothetical planet with orbital period $\sim$17,117 years, semimajor axis $\sim$665 AU, would have N/1 and N/2 period ratios with these four objects. The orbital geometries and dynamics of resonant orbits constrain the orbital plane, the orbital eccentricity and the mass of such a planet, as well as its current location in its orbital path.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Lightcurves of the Karin family asteroids

The Karin family is a very young asteroid family created by an asteroid breakup 5.8 Myr ago. Since the members of this family probably have not experienced significant orbital or collisional evolution yet, it is possible that they still preserve properties of the original family-forming event in terms of their spin state. As we carried out a series of photometric observations of the Karin family asteroids, here we report an analysis result of lightcurves including the rotation period of eleven members. The mean rotation rate of the Karin family members turned out to be much lower than those of NEAs or smaller MBAs (diameter D<12 km), and even lower than that of larger MBAs (D>130 km). We investigated a correlation between the peak-to-peak variation magnitude reduced to zero solar phase angle and the rotation period of the eleven Karin family asteroids, and found a possible trend that elongated members have lower spin rate, and less elongated members have higher spin rate. However, this has to be confirmed by another series of future observations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The current impact flux on Mars and its seasonal variation

We calculate the present-day impact flux on Mars and its variation over the Martian year, using the current data on the orbital distribution of known Mars-crossing minor planets. We adapt the {Ö}pik-Wetherill formulation for calculating collision probabilities, paying careful attention to the non-uniform distribution of the perihelion longitude and the argument of perihelion owed to secular planetary perturbations. We find that these previously neglected non-uniformities have a significant effect on the mean annual impact flux as well as its seasonal variation. The impact flux peaks when Mars is at aphelion, but the near-alignment of Mars' eccentricity vector with the mean direction of the eccentricity vectors of Mars-crossers causes the mean annual impact flux as well as the amplitude of the seasonal variation to be significantly lower than the estimate based on a uniform random distribution of perihelion longitudes of Mars-crossers. We estimate that the flux of large impactors (of absolute magnitude $H<16$) within $\pm30^\circ$ of Mars' aphelion is about four times larger than when the planet is near perihelion. Extrapolation of our results to a model population of meter-size Mars-crossers shows that if these small impactors have a uniform distribution of their angular elements, then their aphelion-to-perihelion impact flux ratio would be as large as 25. These theoretical predictions can be tested with observational data of contemporary impacts that is becoming available from spacecraft currently in orbit about Mars.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Mass Distribution Function of Planets

The distribution of orbital period ratios of adjacent planets in extra-solar planetary systems discovered by the {\it Kepler} space telescope exhibits a peak near $\sim1.5$--$2$, a long tail of larger period ratios, and a steep drop-off in the number of systems with period ratios below $\sim1.5$. We find from this data that the dimensionless orbital separations have an approximately log-normal distribution. Using Hill's criterion for the dynamical stability of two planets, we find an upper bound on planet masses such that the most common planet mass does not exceed $10^{-3.2}m_*$, or about two-thirds Jupiter mass for solar mass stars. Assuming that the mass ratio and the dynamical separation (orbital spacings in units of mutual Hill radius) of adjacent planets are independent random variates, and adopting empirical distributions for these, we use Hill's criterion in a statistical way to estimate the planet mass distribution function from the observed distribution of orbital separations. We find that the planet mass function is peaked in logarithm of mass, with a peak value and standard deviation of $\log m/M_\oplus$ of $\sim(0.6-1.0)$ and $\sim(1.1-1.2)$, respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Debris Distribution in HD 95086 - A Young Analog of HR 8799

HD 95086 is a young early-type star that hosts (1) a 5 MJ planet at the projected distance of 56 AU revealed by direct imaging, and (2) a prominent debris disk. Here we report the detection of 69 um crystalline olivine feature from the disk using the Spitzer/MIPS-SED data covering 55-95 um. Due to the low resolution of MIPS-SED mode, this feature is not spectrally resolved, but is consistent with the emission from crystalline forsterite contributing 5% of the total dust mass. We also present detailed analysis of the disk SED and re-analysis of resolved images obtained by Herschel. Our results suggest that the debris structure around HD 95086 consists of a warm (175 K) belt, a cold (55 K) disk, and an extended disk halo (up to 800 AU), and is very similar to that of HR 8799. We compare the properties of the three debris components, and suggest that HD 95086 is a young analog of HR 8799. We further investigate and constrain single-planet, two-planet, three-planet and four-planet architectures that can account for the observed debris structure and are compatible with dynamical stability constraints. We find that equal-mass four-planet configurations of geometrically spaced orbits, with each planet of mass 5 MJ, could maintain the gap between the warm and cold debris belts, and also be just marginally stable for timescales comparable to the age of the system.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Planetary chaotic zone clearing: destinations and timescales

We investigate the orbital evolution of particles in a planet's chaotic zone to determine their final destinations and their timescales of clearing. There are four possible final states of chaotic particles: collision with the planet, collision with the star, escape, or bounded but non-collision orbits. In our investigations, within the framework of the planar circular restricted three body problem for planet-star mass ratio $μ$ in the range $10^{-9}$ to $10^{-1.5}$, we find no particles hitting the star. The relative frequencies of escape and collision with the planet are not scale-free, as they depend upon the size of the planet. For planet radius $R_p\ge0.001R_H$ where $R_H$ is the planet's Hill radius, we find that most chaotic zone particles collide with the planet for $μ\lesssim10^{-5}$; particle scattering to large distances is significant only for higher mass planets. For fixed ratio $R_p/R_H$, the particle clearing timescale, $T_{cl}$, has a broken power-law dependence on $μ$. A shallower power-law, $T_{cl}\sim μ^{-{1/3}}$, prevails at small $μ$ where particles are cleared primarily by collisions with the planet; a steeper power law, $T_{cl}\simμ^{-{3/2}}$, prevails at larger $μ$ where scattering dominates the particle loss. In the limit of vanishing planet radius, we find $T_{cl}\approx0.024μ^{-{3\over2}}$. The interior and exterior boundaries of the annular zone in which chaotic particles are cleared are increasingly asymmetric about the planet's orbit for larger planet masses; the inner boundary coincides well with the classical first order resonance overlap zone, $Δa_{cl,int}\simeq1.2μ^{0.28}a_p$; the outer boundary is better described by $Δa_{cl,ext}\simeq1.7μ^{0.31}a_p$, where $a_p$ is the planet-star separation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The inner solar system cratering record and the evolution of impactor populations

We review previously published and newly obtained crater size-frequency distributions in the inner solar system. These data indicate that the Moon and the terrestrial planets have been bombarded by two populations of objects. Population 1, dominating at early times, had nearly the same size distribution as the present-day asteroid belt, and produced the heavily cratered surfaces with a complex, multi-sloped crater size-frequency distribution. Population 2, dominating since about 3.8-3.7 Ga, has the same size distribution as near-Earth objects (NEOs), had a much lower impact flux, and produced a crater size distribution characterized by a differential -3 single-slope power law in the crater diameter range 0.02 km to 100 km. Taken together with the results from a large body of work on age-dating of lunar and meteorite samples and theoretical work in solar system dynamics, a plausible interpretation of these data is as follows. The NEO population is the source of Population 2 and it has been in near-steady state over the past ~3.7-3.8 gigayears; these objects are derived from the main asteroid belt by size-dependent non-gravitational effects that favor the ejection of smaller asteroids. However, Population 1 were main belt asteroids ejected from their source region in a size-independent manner, possibly by means of gravitational resonance sweeping during giant planet orbit migration; this caused the so-called Late Heavy Bombardment (LHB). The LHB began some time before ~3.9 Ga, peaked and declined rapidly over the next ~100 to 300 megayears, and possibly more slowly from about 3.8-3.7 Ga to ~2 Ga. A third crater population (Population S) consists of secondary impact craters that can dominate the cratering record at small diameters.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Do Centaurs preserve their source inclinations?

The Centaurs are a population of small, planet-crossing objects in the outer solar system. They are dynamically short-lived and represent the transition population between the Kuiper belt and the Jupiter family short-period comets. Dynamical models and observations of the physical properties of the Centaurs indicate that they may have multiple source populations in the trans-Neptunian region. It has been suggested that the inclination distribution of the Centaurs may be useful in distinguishing amongst these source regions. The Centaurs, however, undergo many close encounters with the giant planets during their orbital evolution; here we show that these encounters can substantially determine the inclination distribution of the Centaurs. Almost any plausible initial inclination distribution of a Kuiper belt source results in Centaurs having inclinations peaked near 10-20 degrees. Our studies also find that the Kuiper belt is an extremely unlikely source of the retrograde Centaur that has been observed.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

On the Non-uniform Distribution of the Angular Elements of Near-Earth Objects

We examine the angular distributions of near-Earth objects (NEOs) which are often regarded as uniform. The apparent distribution of the longitude of ascending node, $Ω$, is strongly affected by well-known seasonal effects in the discovery rate of NEOs. The deviation from the expected $π$-periodicity in the apparent distribution of $Ω$ indicates that its intrinsic distribution is slightly enhanced along a mean direction, $\barΩ=111^\circ$; approximately 53% of NEOs have $Ω$ values within $\pm90^\circ$ of $\barΩ$. We also find that each sub-group of NEOs (Amors, Apollos and Atens) has different observational selection effects which cause different non-uniformities in the apparent distributions of their arguments of perihelion $ω$, and longitudes of perihelion $\varpi$. For their intrinsic distributions, our analysis reveals that the Apollo asteroids have non-uniform $ω$ due to secular dynamics associated with inclination-eccentricity-$ω$ coupling, and the Amors' $\varpi$ distribution is peaked towards the secularly forced eccentricity vector. The Apollos' $ω$ distribution is axial, favoring values near 0 and $180^\circ$; the two quadrants centered at 0 and $180^\circ$ account for 55% of the Apollos' $ω$ values. The Amors' $\varpi$ distribution peaks near $\bar\varpi=4^\circ$; sixty-one percent of Amors have $\varpi$ within $\pm90^\circ$ of this peak. We show that these modest but statistically significant deviations from uniform random distributions of angular elements are owed to planetary perturbations, primarily Jupiter's. It is remarkable that this strongly chaotic population of minor planets reveals the presence of Jupiter in its angular distributions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Predictions for Shepherding Planets in Scattered Light Images of Debris Disks

Planets can affect debris disk structure by creating gaps, sharp edges, warps, and other potentially observable signatures. However, there is currently no simple way for observers to deduce a disk-shepherding planet's properties from the observed features of the disk. Here we present a single equation that relates a shepherding planet's maximum mass to the debris ring's observed width in scattered light, along with a procedure to estimate the planet's eccentricity and minimum semimajor axis. We accomplish this by performing dynamical N-body simulations of model systems containing a star, a single planet, and a disk of parent bodies and dust grains to determine the resulting debris disk properties over a wide range of input parameters. We find that the relationship between planet mass and debris disk width is linear, with increasing planet mass producing broader debris rings. We apply our methods to five imaged debris rings to constrain the putative planet masses and orbits in each system. Observers can use our empirically-derived equation as a guide for future direct imaging searches for planets in debris disk systems. In the fortuitous case of an imaged planet orbiting interior to an imaged disk, the planet's maximum mass can be estimated independent of atmospheric models.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Chaotic exchange of solid material between planetary systems: implications for lithopanspermia

We examine a low energy mechanism for the transfer of meteoroids between two planetary systems embedded in a star cluster using quasi-parabolic orbits of minimal energy. Using Monte Carlo simulations, we find that the exchange of meteoroids could have been significantly more efficient than previously estimated. Our study is relevant to astrobiology as it addresses whether life on Earth could have been transferred to other planetary systems in the solar system's birth cluster and whether life on Earth could have been transferred here from beyond the solar system. In the solar system, the timescale over which solid material was delivered to the region from where it could be transferred via this mechanism likely extended to several hundred million years (as indicated by the 3.8-4.0 Ga epoch of the Late Heavy Bombardment). This timescale could have overlapped with the lifetime of the Solar birth cluster (~100-500 Myr). Therefore, we conclude that lithopanspermia is an open possibility if life had an early start. Adopting parameters from the minimum mass solar nebula, considering a range of planetesimal size distributions derived from observations of asteroids and Kuiper Belt Objects and theoretical coagulation models, and taking into account Oort Cloud formation models, the expected number of bodies with mass > 10 kg that could have been transferred between the Sun and its nearest cluster neighbor could be of the order of 1E14-3E16, with transfer timescales of 10s Myr. We estimate that of the order of 3E8 x l(km) could potentially be life-bearing, where l(km) is the depth of the Earth crust in km that was ejected as the result of the early bombardment.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Planets Near Mean-Motion Resonances

The multiple-planet systems discovered by the Kepler mission exhibit the following feature: planet pairs near first-order mean-motion resonances prefer orbits just outside the nominal resonance, while avoiding those just inside the resonance. We explore an extremely simple dynamical model for planet formation, in which planets grow in mass at a prescribed rate without orbital migration or dissipation. We develop an analytic version of this model for two-planet systems in two limiting cases: the planet mass grows quickly or slowly relative to the characteristic resonant libration time. In both cases the distribution of systems in period ratio develops a characteristic asymmetric peak-trough structure around the resonance, qualitatively similar to that observed in the Kepler sample. We verify this result with numerical integrations of the restricted three-body problem. We show that for the 3:2 resonance, where the observed peak-trough structure is strongest, our simple model is consistent with the observations for a range of mean planet masses 20-100M_{\oplus}. This mass range is higher than expected, by at least a factor of three, from the few Kepler planets with measured masses, but part of this discrepancy could be due to oversimplifications in the dynamical model or uncertainties in the planetary mass-radius relation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The effect of orbital evolution on the Haumea (2003 EL61) collisional family

The Haumea family is currently the only identified collisional family in the Kuiper belt. We numerically simulate the long-term dynamical evolution of the family to estimate a lower limit of the family's age and to assess how the population of the family and its dynamical clustering are preserved over Gyr timescales. We find that the family is not younger than 100 Myr, and its age is at least 1 Gyr with 95% confidence. We find that for initial velocity dispersions of 50-400 m/s, approximately 20-45% of the family members are lost to close encounters with Neptune after 3.5 Gyr of orbital evolution. We apply these loss rates to two proposed models for the formation of the Haumea family, a graze-and-merge type collision between two similarly sized, differentiated KBOs or the collisional disruption of a satellite orbiting Haumea. For the graze-and-merge collision model, we calculate that >85% of the expected mass in surviving family members within 150 m/s of the collision has been identified, but that one to two times the mass of the known family members remains to be identified at larger velocities. For the satellite-break-up model, we estimate that the currently identified family members account for ~50% of the expected mass of the family. Taking observational incompleteness into account, the observed number of Haumea family members is consistent with either formation scenario at the 1 sigma level, however both models predict more objects at larger relative velocities (>150 m/s) than have been identified.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Inclination Mixing in the Classical Kuiper Belt

We investigate the long-term evolution of the inclinations of the known classical and resonant Kuiper belt objects (KBOs). This is partially motivated by the observed bimodal inclination distribution and by the putative physical differences between the low- and high-inclination populations. We find that some classical KBOs undergo large changes in inclination over gigayear timescales, which means that a current member of the low-inclination population may have been in the high-inclination population in the past, and vice versa. The dynamical mechanisms responsible for the time-variability of inclinations are predominantly distant encounters with Neptune and chaotic diffusion near the boundaries of mean motion resonances. We reassess the correlations between inclination and physical properties including inclination time-variability. We find that the size-inclination and color-inclination correlations are less statistically significant than previously reported (mostly due to the increased size of the data set since previous works with some contribution from inclination variability). The time-variability of inclinations does not change the previous finding that binary classical KBOs have lower inclinations than non-binary objects. Our study of resonant objects in the classical Kuiper belt region includes objects in the 3:2, 7:4, 2:1, and eight higher-order mean motion resonances. We find that these objects (some of which were previously classified as non-resonant) undergo larger changes in inclination compared to the non-resonant population, indicating that their current inclinations are not generally representative of their original inclinations. They are also less stable on gigayear timescales.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Long-Term Cycling of Kozai-Lidov Cycles: Extreme Eccentricities and Inclinations Excited by a Distant Eccentric Perturber

Kozai-Lidov oscillations of Jupiter-mass planets, excited by comparable planetary or brown dwarf mass perturbers were recently shown in numerical experiments to be slowly modulated and to exhibit striking features, including extremely high eccentricities and the generation of retrograde orbits with respect to the perturber. Here we solve this problem analytically for the case of a test particle orbiting a host star and perturbed by a distant companion whose orbit is eccentric and highly inclined. We give analytic expressions for the conditions that produce retrograde orbits and high eccentricities. This mechanism likely operates in various systems thought to involve Kozai-Lidov oscillations such as tight binaries, mergers of compact objects, irregular moons of planets and many others. In particular, it could be responsible for exciting eccentricities and inclinations of exo-planetary orbits and be important for understanding the spin-orbit (mis)alignment of hot Jupiters.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Secular resonance sweeping of the main asteroid belt during planet migration

We calculate the eccentricity excitation of asteroids produced by the sweeping $ν_6$ secular resonance during the epoch of planetesimal-driven giant planet migration in the early history of the solar system. We derive analytical expressions for the magnitude of the eccentricity change and its dependence on the sweep rate and on planetary parameters; the $ν_6$ sweeping leads to either an increase or a decrease of eccentricity depending on an asteroid's initial orbit. Based on the slowest rate of $ν_6$ sweeping that allows a remnant asteroid belt to survive, we derive a lower limit on Saturn's migration speed of $\sim0.15\AU\My^{-1}$ during the era that the $ν_6$ resonance swept through the inner asteroid belt (semimajor axis range 2.1--$2.8\AU$). This rate limit is for Saturn's current eccentricity, and scales with the square of Saturn's eccentricity; the limit on Saturn's migration rate could be lower if Saturn's eccentricity were lower during its migration. Applied to an ensemble of fictitious asteroids, our calculations show that a prior single-peaked distribution of asteroid eccentricities would be transformed into a double-peaked distribution due to the sweeping of the $ν_6$. Examination of the orbital data of main belt asteroids reveals that the proper eccentricities of the known bright ($H \leq10.8$) asteroids may be consistent with a double-peaked distribution. If so, our theoretical analysis then yields two possible solutions for the migration rate of Saturn and for the dynamical states of the pre-migration asteroid belt: a dynamically cold state (single-peaked eccentricity distribution with mean of $\sim0.05$) linked with Saturn's migration speed $\sim 4\AU\My^{-1}$, or a dynamically hot state (single-peaked eccentricity distribution with mean of $\sim0.3$) linked with Saturn's migration speed $\sim 0.8\AU\My^{-1}$.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Asymmetric impacts of near-Earth asteroids on the Moon

Recent lunar crater studies have revealed an asymmetric distribution of rayed craters on the lunar surface. The asymmetry is related to the synchronous rotation of the Moon: there is a higher density of rayed craters on the leading hemisphere compared with the trailing hemisphere. Rayed craters represent generally the youngest impacts. The purpose of this paper is to test the hypotheses that (i) the population of Near-Earth asteroids (NEAs) is the source of the impactors that have made the rayed craters, and (ii) that impacts by this projectile population account quantitatively for the observed asymmetry. We carried out numerical simulations of the orbital evolution of a large number of test particles representing NEAs in order to determine directly their impact flux on the Moon. The simulations were done in two stages. In the first stage we obtained encounter statistics of NEAs on the Earth's activity sphere. In the second stage we calculated the direct impact flux of the encountering particles on the surface of the Moon; the latter calculations were confined within the activity sphere of the Earth. A steady-state synthetic population of NEAs was generated from a debiased orbital distribution of the known NEAs. We find that the near-Earth asteroids do have an asymmetry in their impact flux on the Moon: apex-to-antapex ratio of 1.32 +/- 0.01. However, the observed rayed crater distribution's asymmetry is significantly more pronounced: apex-to-antapex ratio of 1.65 +/- 0.16. Our results suggest the existence of an undetected population of slower (low impact velocity) projectiles, such as a population of objects nearly coorbiting with Earth; more observational study of young lunar craters is needed to secure this conclusion.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Locating the planetesimals belts in the multiple-planet systems HD 128311, HD 202206, HD 82943 and HR 8799

In addition to the Sun, six other stars are known to harbor multiple planets and debris disks: HD 69830, HD 38529, HD 128311, HD 202206, HD 82943 and HR 8799. In this paper we set constraints on the location of the dust-producing planetesimals around the latter four systems. We use a radiative transfer model to analyze the spectral energy distributions of the dust disks (including two new Spitzer IRS spectra presented in this paper), and a dynamical model to assess the long-term stability of the planetesimals' orbits. As members of a small group of stars that show evidence of harboring a multiple planets and planetesimals, their study can help us learn about the diversity of planetary systems.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

A record of planet migration in the Main Asteroid Belt

The main asteroid belt lies between the orbits of Mars and Jupiter, but the region is not uniformly filled with asteroids. There are gaps, known as the Kirkwood gaps, in the asteroid distribution in distinct locations that are associated with orbital resonances with the giant planets; asteroids placed in these locations follow chaotic orbits and escape from the asteroid belt. Here we show that the observed distribution of main belt asteroids does not fill uniformly even those regions that are dynamically stable over the age of the solar system. We find a pattern of excess depletion of asteroids, particularly just outward of the Kirkwood Gaps associated with the 5:2, the 7:3, and the 2:1 jovian resonances. These features are not accounted for by planetary perturbations in the current structure of the solar system, but are consistent with dynamical ejection of asteroids by the sweeping of gravitational resonances during the migration of Jupiter and Saturn ~4 gigayears ago.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Chaotic Diffusion of Resonant Kuiper Belt Objects

We carried out extensive numerical orbit integrations to probe the long-term chaotic dynamics of the two strongest mean motion resonances of Neptune in the Kuiper belt, the 3:2 (Plutinos) and 2:1 (Twotinos). Our primary results include a computation of the relative volumes of phase space characterized by large- and small-resonance libration amplitudes, and maps of resonance stability measured by mean chaotic diffusion rate. We find that Neptune's 2:1 resonance has weaker overall long-term stability than the 3:2 -- only ~15% of Twotinos are projected to survive for 4 Gyr, compared to ~27% of Plutinos, based on an extrapolation from our 1-Gyr integrations. We find that Pluto has only a modest effect, causing a ~4% decrease in the Plutino population that survives to 4 Gyr. Given current observational estimates, and assuming an initial distribution of particles proportional to the local phase space volume in the resonance, we conclude that the primordial populations of Plutinos and Twotinos formerly made up more than half the population of the classical and resonant Kuiper Belt. We also conclude that Twotinos were originally nearly as numerous as Plutinos; this is consistent with predictions from early models of smooth giant planet migration and resonance sweeping of the Kuiper Belt, and provides a useful constraint for more detailed models.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Dynamical erosion of the asteroid belt and implications for large impacts in the inner solar system

The cumulative effects of weak resonant and secular perturbations by the major planets produce chaotic behavior of asteroids on long timescales. Dynamical chaos is the dominant loss mechanism for asteroids with diameters D > 10 km in the current asteroid belt. In a numerical analysis of the long term evolution of test particles in the main asteroid belt region, we find that the dynamical loss history of test particles from this region is well described with a logarithmic decay law. In our simulations the loss rate function that is established at t = 1 My persists with little deviation to at least t = 4 Gy. Our study indicates that the asteroid belt region has experienced a significant amount of depletion due to this dynamical erosion - having lost as much as ~50% of the large asteroids - since 1 My after the establishment of the current dynamical structure of the asteroid belt. Because the dynamical depletion of asteroids from the main belt is approximately logarithmic, an equal amount of depletion occurred in the time interval 10-200 My as in 0.2-4 Gy, roughly ~30% of the current number of large asteroids in the main belt over each interval. We find that asteroids escaping from the main belt due to dynamical chaos have an Earth impact probability of ~0.3%. Our model suggests that the rate of impacts from large asteroids has declined by a factor of 3 over the last 3 Gy, and that the present-day impact flux of D > 10 km objects on the terrestrial planets is roughly an order of magnitude less than estimates currently in use in crater chronologies and impact hazard risk assessments.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Two Dynamical Classes of Centaurs

The Centaurs are a transient population of small bodies in the outer solar system whose orbits are strongly chaotic. These objects typically suffer significant changes of orbital parameters on timescales of a few thousand years, and their orbital evolution exhibits two types of behaviors described qualitatively as random-walk and resonance-sticking. We have analyzed the chaotic behavior of the known Centaurs. Our analysis has revealed that the two types of chaotic evolution are quantitatively distinguishable: (1) the random walk-type behavior is well described by so-called generalized diffusion in which the rms deviation of the semimajor axis grows with time t as ~t^H, with Hurst exponent H in the range 0.22--0.95, however (2) orbital evolution dominated by intermittent resonance sticking, with sudden jumps from one mean motion resonance to another, has poorly defined H. We further find that these two types of behavior are correlated with Centaur dynamical lifetime: most Centaurs whose dynamical lifetime is less than ~22 Myr exhibit generalized diffusion, whereas most Centaurs of longer dynamical lifetimes exhibit intermittent resonance sticking. We also find that Centaurs in the diffusing class are likely to evolve into Jupiter-family comets during their dynamical lifetimes, while those in the resonance-hopping class do not.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Are Debris Disks and Massive Planets Correlated?

Using data from the Spitzer Space Telescope Legacy Science Program ``Formation and Evolution of Planetary Systems'' (FEPS), we have searched for debris disks around 9 FGK stars (2-10 Gyr), known from radial velocity (RV) studies to have one or more massive planets. Only one of the sources, HD 38529, has excess emission above the stellar photosphere; at 70 micron the signal-to-noise ratio in the excess is 4.7 while at wavelengths < 30 micron there is no evidence of excess. The remaining sources show no excesses at any Spitzer wavelengths. Applying survival tests to the FEPS sample and the results for the FGK survey published in Bryden et al. (2006), we do not find a significant correlation between the frequency and properties of debris disks and the presence of close-in planets. We discuss possible reasons for the lack of a correlation.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Assessing the massive young Sun hypothesis to solve the warm young Earth puzzle

A moderately massive early Sun has been proposed to resolve the so-called faint early Sun paradox. We calculate the time-evolution of the solar mass that would be required by this hypothesis, using a simple parametrized energy-balance model for Earth's climate. Our calculations show that the solar mass loss rate would need to have been 2-3 orders of magnitude higher than present for a time on the order of ~2 Gy. Such a mass loss history is significantly at variance (both in timescale and in the magnitude of the mass loss rates) with that inferred from astronomical observations of mass loss in younger solar analogues. While suggestive, the astronomical data cannot completely rule out the possibility that the Sun had the required mass loss history; therefore, we also examine the effects of the hypothetical historical solar mass loss on orbital dynamics in the solar system, with a view to identifying additional tests of the hypothesis. Planetary and satellite orbits provide a few tests, but these are weak or non-unique.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

Observational Limits on a Cold Distant Kuiper Belt

Almost all of the $>600$ known Kuiper belt objects (KBOs) have been discovered within 50 AU of the Sun. One possible explanation for the observed lack of KBOs beyond 50 AU is that the distant Kuiper belt is dynamically very cold, and thus thin enough on the sky to have slipped between previous deep survey fields. We have completed a survey designed to search for a dynamically cold distant Kuiper belt near the invariable plane of the Solar system. In $2.3 {\rm deg}^2$ we have discovered a total of 33 KBOs and 1 Centaur, but no objects in circular orbits beyond 50 AU. We find that we can exclude at 95% CL the existence of a distant disk inclined by $i\le1\arcdeg$ to the invariable plane and containing more than 1.2 times as many $D>185$ km KBOs between 50 and 60 AU as the observed inner Kuiper belt, if the distant disk is thinner than $σ=1\fdg75$.

0 citations0 reviewsNo code linkedOpen access