Source author record

Yanqin Wu

Yanqin Wu 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

24works
2topics
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

24 published item(s)

preprint2021arXiv

Super-Earths and sub-Neptunes are Insensitive to Stellar Metallicity

Kepler planets (including super-Earths and sub-Neptunes, from 1 to 4 Earth radii) are likely formed before the gaseous proto-planetary disks have dissipated, as are the Jovian planets. If the metal content in these disks resembles that in the host stars, one might expect Kepler planets to occur more frequently, and to be more massive, around metal-rich stars. Contrary to these expectations, we find that the radii of Kepler planets (a proxy for mass) are independent of host metallicity. Previous claims that larger planets prefer more metal-rich stars can be adequately explained by the combined facts that more massive stars tend to host bigger planets, and that more massive stars are also more metal rich in the Kepler sample. We interpret this independence as that the mass of a Kepler planet is not determined by the availability of solids, but is instead regulated by an as yet unknown process. Moreover, we find that the occurrence rates of Kepler planets rise only weakly with stellar metallicity, a trend that is further flattened when the influence of close stellar binaries is accounted for. We explain this weak dependence, in contrast to the strong dependence exhibited by Jovian planets, using a phenomenological model, wherein the masses of proto-planetary disks have a much larger spread than the spread in stellar metallicity, and wherein the formation of Jovian planets requires disks that contain some $5$ times more solid than that needed to form Kepler planets. This model predicts that stars more metal-poor than half-solar should rarely host any Kepler planets.

preprint2020arXiv

The surprisingly low carbon mass in the debris disk around HD 32297

Gas has been detected in a number of debris disks. It is likely secondary, i.e. produced by colliding solids. Here, we report ALMA Band 8 observations of neutral carbon in the CO-rich debris disk around the 15--30 Myr old A-type star HD 32297. We find that C$^0$ is located in a ring at $\sim$110 au with a FWHM of $\sim$80 au, and has a mass of $(3.5\pm0.2)\times10^{-3}$ M$_\oplus$. Naively, such a surprisingly small mass can be accumulated from CO photo-dissociation in a time as short as $\sim$10$^4$ yr. We develop a simple model for gas production and destruction in this system, properly accounting for CO self-shielding and shielding by neutral carbon, and introducing a removal mechanism for carbon gas. We find that the most likely scenario to explain both C$^0$ and CO observations, is one where the carbon gas is rapidly removed on a timescale of order a thousand years and the system maintains a very high CO production rate of $\sim$15 M$_\oplus$ Myr$^{-1}$, much higher than the rate of dust grind-down. We propose a possible scenario to meet these peculiar conditions: the capture of carbon onto dust grains, followed by rapid CO re-formation and re-release. In steady state, CO would continuously be recycled, producing a CO-rich gas ring that shows no appreciable spreading over time. This picture might be extended to explain other gas-rich debris disks.

preprint2020arXiv

Tidal disruption versus planetesimal collisions as possible origins for the dispersing dust cloud around Fomalhaut

Recent analysis suggests that the faint optical point source observed around Fomalhaut from 2004-2014 (Fomalhaut b) is gradually fading and expanding, supporting the case that it may be a dispersing dust cloud resulting from the sudden disruption of a planetesimal. These types of disruptions may arise from catastrophic collisions of planetesimals, which are perturbed from their original orbits in the Fomalhaut dust ring by nearby giant planets. However, disruptions can also occur when the planetesimals pass within the tidal disruption field of the planet(s) that perturbed them in the first place, similar to the Shoemaker-Levy event observed in the Solar System. Given that a gravitationally focusing giant planet has a much larger interaction cross-section than a planetesimal, tidal disruption events can match or outnumber planetesimal collision events in realistic regions of parameter space. Intriguingly, the Fomalhaut dust cloud offers an opportunity to directly distinguish between these scenarios. A tidal disruption scenario leads to a very specific prediction of ephemerides for the planet causing the event. At a most probable mass of 66 Mearth, a semi-major axis of 117 AU, and a system age of 400-500 Myr, this planet would be readily detectable with the James Webb Space Telescope. The presence or absence of this planet at the specific, predicted position is therefore a distinctive indicator of whether the dispersing cloud originated from a collision of two planetesimals or from the disruption of a planetesimal in the tidal field of a giant planet.

preprint2016arXiv

Forming the Cold Classical Kuiper Belt in a light Disk

Large Kuiper Belt Objects are conventionally thought to have formed out of a massive planetesimal belt that is a few thousand times its current mass. Such a picture, however, is incompatible with multiple lines of evidence. Here, we present a new model for the conglomeration of Cold Classical Kuiper belt objects, out of a solid belt only a few times its current mass, or a few percent of the solid density in a Minimum Mass Solar Nebula. This is made possible by depositing most of the primordial mass in grains of size centimetre or smaller. These grains collide frequently and maintain a dynamically cold belt out of which large bodies grow efficiently: an order-unity fraction of the solid mass can be converted into large bodies, in contrast to the ~0.1% efficiency in conventional models. Such a light belt may represent the true outer edge of the Solar system, and it may have effectively halted the outward migration of Neptune. In addition to the high efficiency, our model can also produce a mass spectrum that peaks at an intermediate size, similar to the observed Cold Classicals, if one includes the effect of cratering collisions. In particular, the observed power-law break observed at ~30 km for Cold Classicals, one that has been interpreted as a result of collisional erosion, may be primordial in origin.

preprint2016arXiv

On the mass and origin of Chariklo's rings

Observations in 2013 and 2014 of the Centaur 10199 Chariklo and its ring system consistently indicated that the radial width of the inner, more massive ring varies with longitude. That strongly suggests that this ring has a finite eccentricity despite the fast differential precession that Chariklo's large quadrupole moment should induce. If the inferred apse alignment is maintained by the ring's self-gravity, as it is for the Uranian rings, we estimate a ring mass of a few times 10^16 g and a typical particle size of a few meters. These imply a short collisional spreading time of ~10^5 years, somewhat shorter than the typical Centaur dynamical lifetime of a few Myrs and much shorter than the age of the solar system. In light of this time constraint, we evaluate previously suggested ring formation pathways including collisional ejection and satellite disruption. We also investigate in detail a contrasting formation mechanism, the lofting of dust particles off Chariklo's surface into orbit via outflows of sublimating CO and/or N_2 triggered after Chariklo was scattered inward by giant planets. This latter scenario predicts that rings should be common among 100-km class Centaurs but rare among Kuiper belt objects and smaller Centaurs. It also predicts that Centaurs should show seasonal variations in cometary activity with activity maxima occurring shortly after equinox.

preprint2016arXiv

Warm Jupiters are less lonely than hot Jupiters: close neighbours

Exploiting the Kepler transit data, we uncover a dramatic distinction in the prevalence of sub-Jovian companions, between systems that contain hot Jupiters (periods inward of 10 days) and those that host warm Jupiters (periods between 10 and 200 days). Hot Jupiters, with the singular exception of WASP-47b, do not have any detectable inner or outer planetary companions (with periods inward of 50 days and sizes down to $2 R_{\rm Earth}$). Restricting ourselves to inner companions, our limits reach down to $1 R_{\rm Earth}$. In stark contrast, half of the warm Jupiters are closely flanked by small companions. Statistically, the companion fractions for hot and warm Jupiters are mutually exclusive, particularly in regard to inner companions. The high companion fraction of warm Jupiters also yields clues to their formation. The warm Jupiters that have close-by siblings should have low orbital eccentricities and low mutual inclinations. The orbital configurations of these systems are reminiscent of those of the low-mass, close-in planetary systems abundantly discovered by the Kepler mission. This, and other arguments, lead us to propose that these warm Jupiters are formed in-situ. There are indications that there may be a second population of warm Jupiters with different characteristics. In this picture, WASP-47b could be regarded as the extending tail of the in-situ warm Jupiters into the hot Jupiter region, and does not represent the generic formation route for hot Jupiters.

preprint2015arXiv

Atmospheres of low-mass planets: the "boil-off"

We show that, for a low-mass planet that orbits its host star within a few tenths of an AU (like the majority of the {\it Kepler} planets), the atmosphere it was able to accumulate while embedded in the proto-planetary disk may not survive unscathed after the disk disperses. This gas envelope, if more massive than a few percent of the core (with a mass below $10 M_\oplus$), has a cooling time that is much longer than the time-scale on which the planet exits the disk. As such, it could not have contracted significantly from its original size, of order the Bondi radius. So a newly exposed proto-planet would be losing mass via a Parker wind that is catalyzed by the stellar continuum radiation. This represents an intermediate stage of mass-loss, occurring soon after the disc has dispersed, but before the EUV/X-ray driven photoevaporation becomes relevant. The surface mass-loss induces a mass movement within the envelope that advects internal heat outward. As a result, the planet atmosphere rapidly cools down and contracts, until it has reached a radius of order $0.1$ Bondi radius, at which time the mass-loss effectively shuts down. Within a million years after the disk disperses, we find a planet that has only about ten percent of its original envelope, and a Kelvin-Helmholtz time that is much longer than its actual age. We suggest that this "boil-off" process may be partially responsible for the lack of planets above a radius of $2.5 R_\oplus$ in the {\it Kepler} data, provided planet formation results in initial envelope masses of tens of percent.

preprint2015arXiv

Conglomeration of kilometre-sized planetesimals

We study the efficiency of forming large bodies, starting from a sea of equal-sized planetesimals. This is likely one of the earlier steps of planet formation and relevant for the formation of the asteroid belt, the Kuiper belt and extra-solar debris disks. Here we consider the case that the seed planetesimals do not collide frequently enough for dynamical collisional to be important (the collisionless limit), using a newly constructed conglomeration code, and by carefully comparing numerical results with analytical scalings. In the absence of collisional cooling, as large bodies grow by accreting small bodies, the velocity dispersion of the small bodies ($u$) is increasingly excited. Growth passes from the well-known run-away stage (when $u$ is higher than the big bodies' hill velocity) to the newly discovered trans-hill stage (when $u$ and big bodies both grow, but $u$ remains at the big bodies' hill velocity). We find, concurring with the analytical understandings developed in Lithwick (2014), as well as previous numerical studies, that a size spectrum $dn/dR \propto R^{-4}$ results, and that the formation efficiency, defined as mass fraction in bodies much larger than the initial size, is $\sim {\mathrm{ a \,\, few}}\times R_\odot/a$, or $\sim 10^{-3}$ at the distance of the Kuiper belt. We argue that this extreme inefficiency invalidates the conventional conglomeration model for the formation of both our Kuiper belt and extra-solar debris disks. New theories, possibly involving direct gravitational collapse, or strong collisional cooling of small planetesimals, are required.

preprint2015arXiv

Know The Star, Know the Planet. IV. A Stellar Companion to the Host star of the Eccentric Exoplanet HD 8673b

HD 8673 hosts a massive exoplanet in a highly eccentric orbit (e=0.723). Based on two epochs of speckle interferometry a previous publication identified a candidate stellar companion. We observed HD 8673 multiple times with the 10 m Keck II telescope, the 5 m Hale telescope, the 3.63 m AEOS telescope and the 1.5m Palomar telescope in a variety of filters with the aim of confirming and characterizing the stellar companion. We did not detect the candidate companion, which we now conclude was a false detection, but we did detect a fainter companion. We collected astrometry and photometry of the companion on six epochs in a variety of filters. The measured differential photometry enabled us to determine that the companion is an early M dwarf with a mass estimate of 0.33-0.45 M?. The companion has a projected separation of 10 AU, which is one of the smallest projected separations of an exoplanet host binary system. Based on the limited astrometry collected, we are able to constrain the orbit of the stellar companion to a semi-major axis of 35{60 AU, an eccentricity ? 0.5 and an inclination of 75{85?. The stellar companion has likely strongly in uenced the orbit of the exoplanet and quite possibly explains its high eccentricity.

preprint2015arXiv

Spacing of Kepler Planets: Sculpting by Dynamical Instability

We study the orbital architecture of multi-planet systems detected by the Kepler transit mission using N-body simulations, focusing on the orbital spacing between adjacent planets in systems showing four or more transiting planets. We find that the observed spacings are tightly clustered around 12 mutual Hill radii, when transit geometry and sensitivity limits are accounted for. In comparison, dynamical integrations reveal that the minimum spacing required for systems of similar masses to survive dynamical instability for as long as a billion years is, ~10 if all orbits are circular and coplanar, and ~12 if planetary orbits have eccentricities ~0.02 (a value suggested by studies of planet transit-time-variations). This apparent coincidence, between the observed spacing and the theoretical stability threshold, leads us to propose that typical planetary systems were formed with even tighter spacing, but most, except for the widest ones, have undergone dynamical instability, and are pared down to a more anemic version of their former selves, with fewer planets and larger spacings. So while the high multiple systems (five or more transiting planets) are primordial systems that remain stable, the single or double planetary systems, abundantly discovered by the Kepler mission, may be the descendants of more closely packed high multiple systems. If this hypothesis is correct, we infer that the formation environment of Kepler systems should be different from that of the terrestrial planets.

preprint2015arXiv

The 3D Flow Field Around an Embedded Planet

3D modifications to the well-studied 2D flow topology around an embedded planet have the potential to resolve long-standing problems in planet formation theory. We present a detailed analysis of the 3D isothermal flow field around a 5 Earth-mass planet on a fixed circular orbit, simulated using our multi-GPU hydrodynamics code PEnGUIn. We find that, overall, the horseshoe region has a columnar structure extending vertically much beyond the Hill sphere of the planet. This columnar structure is only broken for some of the widest horseshoe streamlines, along which high altitude fluid descends rapidly into the planet's Bondi sphere, performs one horseshoe turn, and exits the Bondi sphere radially in the midplane. A portion of this flow exits the horseshoe region altogether, which we refer to as the "transient" horseshoe flow. The flow continues as it rolls up into a pair of up-down symmetric horizontal vortex lines shed into the wake of the planet. This flow, unique to 3D, affects both planet accretion and migration. It prevents the planet from sustaining a hydrostatic atmosphere due to its intrusion into the Bondi sphere, and leads to a significant corotation torque on the planet, unanticipated by 2D analysis. In the reported simulation, starting with a $Σ\sim r^{-3/2}$ radial surface density profile, this torque is positive and partially cancels with the negative differential Lindblad torque, resulting in a factor of 3 slower planet migration rate. Finally, we report 3D effects can be suppressed by a sufficiently large disk viscosity, leading to results similar to 2D.

preprint2014arXiv

A Statistical Reconstruction of the Planet Population Around Kepler Solar-Type Stars

Using the cumulative catalog of planets detected by the NASA Kepler mission, we reconstruct the intrinsic occurrence of Earth- to Neptune-size (1 - 4$R_{\oplus}$) planets and their distributions with radius and orbital period. We analyze 76,711 solar-type ($0.8<R_*/R_{\odot}<1.2 $) stars with 430 planets on 20-200~d orbits, excluding close-in planets that may have been affected by the proximity to the host star. Our analysis considers errors in planet radii and includes an "iterative simulation" technique that does not bin the data. We find a radius distribution that peaks at 2-2.8 Earth radii, with lower numbers of smaller and larger planets. These planets are uniformly distributed with logarithmic period, and the mean number of such planets per star is $0.46 \pm 0.03$. The occurrence is $\sim 0.66$ if planets interior to 20~d are included. We estimate the occurrence of Earth-size planets in the "habitable zone" (defined as 1-2$R_{\oplus}$, 0.99-1.7 AU for solar-like stars) as $6.4^{+3.4}_{-1.1} \%$. Our results largely agree with those of Petigura et al. (2013), although we find a higher occurrence of 2.8-4 Earth-radii planets. The reasons for this excess are the inclusion of errors in planet radius, updated Huber et al. (2014) stellar parameters, and also the exclusion of planets which may have been affected by proximity to the host star.

preprint2014arXiv

Frequency of Close Companions among Kepler Planets - a TTV study

A transiting planet exhibits sinusoidal transit-time-variations (TTVs) if perturbed by a companion near a mean-motion-resonance (MMR). We search for sinusoidal TTVs in more than 2600 Kepler candidates, using the publicly available Kepler light-curves (Q0-Q12). We find that the TTV fractions rise strikingly with the transit multiplicity. Systems where four or more planets transit enjoy roughly five times higher TTV fraction than those where a single planet transits, and about twice higher than those for doubles and triples. In contrast, models in which all transiting planets arise from similar dynamical configurations predict comparable TTV fractions among these different systems. One simple explanation for our results is that there are at least two different classes of Kepler systems, one closely packed and one more sparsely populated.

preprint2014arXiv

Stability of the Kepler-11 System and its Origin

A significant fraction of Kepler systems are closely-packed, largely coplanar and circular. We study the stability of a 6-planet system, Kepler-11, to gain insights on the dynamics and formation history of such systems. Using a technique called `frequency maps' as fast indicators for long-term stability, we explore the stability of Kepler-11 system by analyzing the neighbourhood space around its orbital parameters. Frequency maps provide a visual representation of chaos and stability, and their dependence on orbital parameters. We find that the current system is stable, but lies within a few percent of several dynamically dangerous 2-body mean-motion resonances. Planet eccentricities are restricted below a small value, $\sim 0.04$, for long-term stability, but planet masses can be more than twice their reported values (thus, allowing for the possibility of mass-loss by past photoevaporation). Based on our frequency maps, we speculate on the origin for instability in closely-packed systems. We then proceed to investigate how the system can have been assembled. The stability constraints on Kepler-11 (mainly, eccentricity constraints) suggest that if the system were assembled in-situ, a dissipation mechanism must have been at work to neutralize eccentricity excitation. On the other hand, if migration was responsible for assembling the planets, there has to be little differential migration among the planets, to avoid them either getting trapped into mean motion resonances, or crashing into each other.

preprint2013arXiv

Density and Eccentricity of Kepler Planets

We analyze the transit timing variations obtained by the Kepler mission for 22 sub-jovian planet pairs (17 published, 5 new) that lie close to mean motion resonances. We find that the TTV phases for most of these pairs lie close to zero, consistent with an eccentricity distribution that has a very low RMS value of e ~ 0.01; but about a quarter of the pairs possess much higher eccentricities, up to 0.1 - 0.4. For the low-eccentricity pairs, we are able to statistically remove the effect of eccentricity to obtain planet masses from TTV data. These masses, together with those measured by radial velocity, yield a best fit mass-radius relation M~3 M_E (R/R_E). This corresponds to a constant surface escape velocity of 20km/s. We separate the planets into two distinct groups, "mid-sized" (those greater than 3 R_E), and "compact" (those smaller). All mid-sized planets are found to be less dense than water and therefore contain extensive H/He envelopes, likely comparable in mass to that of their cores. We argue that these planets have been significantly sculpted by photoevaporation. Surprisingly, mid-sized planets are discovered exclusively around stars more massive than 0.8 M_sun. The compact planets, on the other hand, are often denser than water. Combining our density measurements with those from radial velocity studies, we find that hotter compact planets tend to be denser, with the hottest ones reaching rock density. Moreover, hotter planets tend to be smaller in sizes. These results can be explained if the compact planets are made of rocky cores overlaid with a small amount of hydrogen, < 1% in mass, with water contributing little to their masses or sizes. Photoevaporation has exposed bare rocky cores in cases of the hottest planets. Our conclusion that these planets are likely not water-worlds contrasts with some previous studies.

preprint2013arXiv

Kepler planets: a tale of evaporation

(Abridged) Inspired by the Kepler planet discoveries, we consider the thermal contraction of planets close to their parent star, under the influence of evaporation. The mass-loss rates are based on hydrodynamic models of evaporation that include both X-ray and EUV irradiation. We find that only low-mass planets with hydrogen envelopes are significantly affected by evaporation, with evaporation being able to remove massive hydrogen envelopes inward of 0.1 AU for Neptune-mass objects. We construct a theoretical population of planets with varying core masses, envelope masses, orbital separations, and stellar spectral types, and compare these against the sizes and densities measured for low-mass planets, both in the Kepler mission and from radial velocity surveys. This exercise leads us to conclude that evaporation is the driving force of evolution for close-in Kepler planets. In fact, some 50% of the Kepler planet candidates may have been significantly eroded. Evaporation explains two striking correlations observed in these objects: a lack of large radius/low density planets close to the stars, and a bimodal distribution in planet sizes with a deficit of planets around 2R_E. Planets that have experienced high X-ray exposures are generally smaller than this size, and those with lower X-ray exposures are typically larger. A bimodal distribution is naturally explained by the evaporation model, where, depending on their X-ray exposure, close-in planets can either hold on to hydrogen envelopes 1% in mass, or be stripped entirely. To quantitatively reproduce the observed features, we argue that not only do low-mass Kepler planets need to be made of rocky cores overlaid with hydrogen envelopes, but few of them should have initial masses above 20 M_E, and the majority of them should have core masses of a few Earth masses.

preprint2013arXiv

Secular chaos and its application to Mercury, hot Jupiters, and the organization of planetary systems

In the inner solar system, the planets' orbits evolve chaotically, driven primarily by secular chaos. Mercury has a particularly chaotic orbit, and is in danger of being lost within a few billion years. Just as secular chaos is reorganizing the solar system today, so it has likely helped organize it in the past. We suggest that extrasolar planetary systems are also organized to a large extent by secular chaos. A hot Jupiter could be the end state of a secularly chaotic planetary system reminiscent of the solar system. But in the case of the hot Jupiter, the innermost planet was Jupiter- (rather than Mercury-) sized, and its chaotic evolution was terminated when it was tidally captured by its star. In this contribution, we review our recent work elucidating the physics of secular chaos and applying it to Mercury and to hot Jupiters. We also present new results comparing the inclinations of hot Jupiters thus produced with observations.

preprint2012arXiv

Extracting Planet Mass and Eccentricity From TTV data

Most planet pairs in the Kepler data that have measured transit time variations (TTV) are near first-order mean-motion resonances. We derive analytical formulae for their TTV signals. We separate planet eccentricity into free and forced parts, where the forced part is purely due to the planets' proximity to resonance. This separation yields simple analytical formulae. The phase of the TTV depends sensitively on the presence of free eccentricity: if the free eccentricity vanishes, the TTV will be in phase with the longitude of conjunctions. This effect is easily detectable in current TTV data. The amplitude of the TTV depends on planet mass and free eccentricity, and it determines planet mass uniquely only when the free eccentricity is sufficiently small. We proceed to analyze the TTV signals of six short period Kepler pairs. We find that three of these pairs (Kepler-18,24,25) have TTV phase consistent with zero. The other three (Kepler-23,28,32) have small TTV phases, but ones that are distinctly non-zero. We deduce that the free eccentricities of the planets are small, < 0.01, but not always vanishing. Furthermore, as a consequence of this, we deduce that the true masses of the planets are fairly accurately determined by the TTV amplitudes, within a factor < 2. The smallness of the free eccentricities suggests that the planets have experienced substantial dissipation. This is consistent with the hypothesis that the observed pile-up of Kepler pairs near mean-motion resonances is caused by resonant repulsion. But the fact that some of the planets have non-vanishing free eccentricity suggests that after resonant repulsion occurred there was a subsequent phase in the planets' evolution when their eccentricities were modestly excited, perhaps by interplanetary interactions.

preprint2012arXiv

Ohmic Heating Suspends, not Reverses, the Cooling Contraction of Hot Jupiters

We study the radius evolution of close-in extra-solar jupiters under Ohmic heating, a mechanism that was recently proposed to explain the large observed sizes of many of these planets. Planets are born with high entropy and they subsequently cool and contract. We focus on two cases: first, that ohmic heating commences when the planet is hot (high entropy); and second, that it commences after the planet has cooled. In the former case, we use analytical scalings and numerical experiments to confirm that Ohmic heating is capable of suspending the cooling as long as a few percent of the stellar irradiation is converted into Ohmic heating, and the planet has a surface wind that extends to pressures of ~10 bar or deeper. For these parameters, the radii at which cooling is stalled are consistent with (or larger than) the observed radii of most planets. The only two exceptions are WASP-17b and HAT-P-32b. In contrast to the high entropy case, we show that Ohmic heating cannot significantly re-inflate planets after they have already cooled. This leads us to suggest that the diversity of radii observed in hot jupiters may be partially explained by the different epochs at which they are migrated to their current locations.

preprint2012arXiv

On the unusual gas composition in the Beta Pictoris debris disk

The metallic gas associated with the Beta Pic debris disk is not believed to be primordial, but arises from the destruction of dust grains. Recent observations have shown that carbon and oxygen in this gas are exceptionally overabundant compared to other elements, by some 400 times. We study the origin of this enrichment under two opposing hypothesis, preferential production, where the gas is produced with the observed unusual abundance (as may happen if gas is produced by photo-desorption from C/O-rich icy grains), and preferential depletion, where the gas evolves to the observed state from an original solar abundance (if outgassing occurs under high-speed collisions) under a number of dynamical processes. We include ... ... We find ... ...

preprint2012arXiv

Resonant Repulsion of Kepler Planet Pairs

Planetary systems discovered by the Kepler space telescope exhibit an intriguing feature. While the period ratios of adjacent low-mass planets appear largely random, there is a significant excess of pairs that lie just wide of resonances and a deficit on the near side. We demonstrate that this feature naturally arises when two near-resonant planets interact in the presence of weak dissipation that damps eccentricities. The two planets repel each other as orbital energy is lost to heat. This moves near-resonant pairs just beyond resonance, by a distance that reflects the integrated dissipation they experienced over their lifetimes. We find that the observed distances may be explained by tides if tidal dissipation is unexpectedly efficient (tidal quality factor ~10). Once the effect of resonant repulsion is accounted for, the initial orbits of these low mass planets show little preference for resonances. This could constrain their origin.

preprint2011arXiv

Planetesimals in Debris Disks of Sun-like Stars

Observations of dusty debris disks can be used to test theories of planetesimal coagulation. Planetesimals of sizes up to a couple thousand kms are embedded in these disks and their mutual collisions generate the small dust grains that are observed. The dust luminosities, when combined with information on the dust spatial extent and the system age, can be used to infer initial masses in the planetesimal belts. Carrying out such a procedure for a sample of debris disks around Sun-like stars, we reach the following two conclusions. First, if we assume that colliding planetesimals satisfy a primordial size spectrum of the form dn/ds ~ s^{-q}, observed disks strongly favor a value of q between 3.5 and 4, while both current theoretical expectations and statistics of Kuiper belt objects favor a somewhat larger value. Second, number densities of planetesimals are two to three orders of magnitude higher in detected disks than in the Kuiper belt, for comparably-sized objects. This is a surprise for the coagulation models. It would require a similar increase in the solid surface density of the primordial disk over that of the Minimum Mass Solar Nebula, which is unreasonable. Both of our conclusions are driven by the need to explain the presence of bright debris disks at a few Gyrs of age.

preprint2010arXiv

Secular Chaos and the Production of Hot Jupiters

In a planetary system with two or more well-spaced, eccentric, inclined planets, secular interactions may lead to chaos. The innermost planet may gradually become very eccentric and/or inclined, as a result of the secular degrees of freedom drifting towards equipartition of angular momentum deficit. Secular chaos is known to be responsible for the eventual destabilization of Mercury in our own Solar System. Here we focus on systems with three giant planets. We characterize the secular chaos and demonstrate the criterion for it to occur, but leave a detailed understanding of secular chaos to a companion paper (Lithwick & Wu, 2010). After an extended period of eccentricity diffusion, the inner planet's pericentre can approach the star to within a few stellar radii. Strong tidal interactions and ensuing tidal dissipation extracts orbital energy from the planet and pulls it inward, creating a hot Jupiter. In contrast to other proposed channels for the production of hot Jupiters, such a scenario (which we term "secular migration") explains a range of observations: the pile-up of hot Jupiters at 3-day orbital periods, the fact that hot Jupiters are in general less massive than other RV planets, that they may have misaligned inclinations with respect to stellar spin, and that they have few easily detectable companions (but may have giant companions in distant orbits). Secular migration can also explain close-in planets as low in mass as Neptune; and an aborted secular migration can explain the "warm Jupiters" at intermediate distances. In addition, the frequency of hot Jupiters formed via secular migration increases with stellar age. We further suggest that secular chaos may be responsible for the observed eccentricities of giant planets at larger distances, and that these planets could exhibit significant spin-orbit misalignment.

preprint2010arXiv

Theory of Secular Chaos and Mercury's Orbit

We study the chaotic orbital evolution of planetary systems, focusing on secular (i.e., orbit-averaged) interactions, because these often dominate on long timescales. We first focus on the evolution of a test particle that is forced by multiple massive planets. To linear order in eccentricity and inclination, its orbit precesses with constant frequencies. But nonlinearities can shift the frequencies into and out of secular resonance with the planets' eigenfrequencies, or with linear combinations of those frequencies. The overlap of these nonlinear secular resonances drive secular chaos in planetary systems. We quantify the resulting dynamics for the first time by calculating the locations and widths of nonlinear secular resonances. When results from both analytical calculations and numerical integrations are displayed together in a newly developed "map of the mean momenta" (MMM), the agreement is excellent. This map is particularly revealing for non-coplanar planetary systems and demonstrates graphically that chaos emerges from overlapping secular resonances. We then apply this newfound understanding to Mercury. Previous numerical simulations have established that Mercury's orbit is chaotic, and that Mercury might even collide with Venus or the Sun. We show that Mercury's chaos is primarily caused by the overlap between resonances that are combinations of four modes, the Jupiter-dominated eccentricity mode, the Venus-dominated inclination mode and Mercury's free eccentricity and inclination. Numerical integration of the Solar system confirms that a slew of these resonant angles alternately librate and circulate. We are able to calculate the threshold for Mercury to become chaotic: Jupiter and Venus must have eccentricity and inclination of a few percent. Mercury appears to be perched on the threshold for chaos.