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

Yoram Lithwick

Yoram Lithwick 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-ph.SRphysics.flu-dynastro-ph.COastro-phastro-ph.GAastro-ph.HEnlin.CDphysics.ao-ph

Catalog footprint

What is connected

30works
9topics
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

30 published item(s)

preprint2021arXiv

The Criterion for Chaos in Three-Planet Systems

We establish the criterion for chaos in three-planet systems, for systems similar to those discovered by the Kepler spacecraft. Our main results are as follows: (i) The simplest criterion, which is based on overlapping mean motion resonances MMR's), only agrees with numerical simulations at a very crude level. (ii) Much greater accuracy is attained by considering neighboring MMR's that do not overlap. We work out the width of the chaotic zones around each of the neighbors, and also provide simple approximate expressions for the widths. (iii) Even greater accuracy is provided by the overlap of three-body resonances (3BR's), which accounts for fine-grained structure seen in maps from N-body simulations, and also predicts the Lyapunov times. Previous studies conflict on whether overlap of MMR's or of 3BR's drive interplanetary chaos. We show that both do, and in fact they are merely different ways of looking at the same effect. (iv) We compare both criteria with high-resolution maps of chaos from N-body simulations, and show that they agree at a high level of detail.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Convection with Misaligned Gravity and Rotation: Simulations and Rotating Mixing Length Theory

We present numerical simulations, using two complementary setups, of rotating Boussinesq thermal convection in a three-dimensional Cartesian geometry with misaligned gravity and rotation vectors. This model represents a small region at a non-polar latitude in the convection zone of a star or planet. We investigate the effects of rotation on the bulk properties of convection at different latitudes, focusing on determining the relation between the heat flux and temperature gradient. We show that our results may be interpreted using rotating mixing length theory (RMLT). The simplest version of RMLT (due to Stevenson) considers the single mode that transports the most heat. This works reasonably well in explaining our results, but there is a systematic departure from these predictions (up to approximately $30\%$ in the temperature gradient) at mid-latitudes. We develop a more detailed treatment of RMLT that includes the transport afforded by multiple modes, and we show that this accounts for most of the systematic differences. We also show that convectively-generated zonal flows and meridional circulations are produced in our simulations, and that their properties depend strongly on the dimensions of the box. These flows also affect the heat transport, contributing to departures from RMLT at some latitudes. However, we find the theoretical predictions of the multi-mode theory for the mid-layer temperature gradient, the root-mean-square (RMS) vertical velocity, the RMS temperature fluctuation, and the spatial spectrum of the heat transport at different latitudes, are all in reasonably good agreement with our numerical results when zonal flows are small.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Inner Boundary Condition in Quasi-Lagrangian Simulations of Accretion Disks

In simulations of viscously evolving accretion disks, the inner boundary condition is particularly important. If treated incorrectly, it induces incorrect behavior very quickly, because the viscous time is shortest near the inner boundary. Recent work has determined the correct inner boundary in Eulerian simulations. But in quasi-Lagrangian simulations (e.g., SPH, moving mesh, and mesh-less), where the inner boundary is modeled by removing mass within a finite zone, the inner density profile typically becomes anomalously depleted. Here we show how the boundary condition should be applied in such codes, via a simple modification of the usual approach: when one removes mass, one must speed up the remaining material so that the disk's angular momentum is unchanged. We show with both 1D and 2D moving-mesh (AREPO) simulations that this scheme works as desired in viscously evolving disks. It produces no spurious density depletions and is independent of the mass removal rate, provided that the disk is adequately resolved and that the mass removal rate is not so extreme as to trigger instabilities. This "torque-free" mass removal technique permits the use of quasi-Lagrangian codes to simulate viscously evolving disks, while including a variety of additional effects. As an example, we apply our scheme to a 2D simulation of an accretion disk perturbed by a very massive planet, in which the disk is evolved to viscous steady state.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Long-Lived Eccentric Modes in Circumbinary Disks

Hydrodynamical simulations show that circumbinary disks become eccentric, even when the binary is circular. Here we demonstrate that, in steady state, the disk's eccentricity behaves as a long-lived free mode trapped by turning points that naturally arise from a continuously truncated density profile. Consequently, both the disk's precession rate and eccentricity profile may be calculated via the simple linear theory for perturbed pressure-supported disks. By formulating and solving the linear theory we find that (i) surprisingly, the precession rate is roughly determined by the binary's quadrupole, even when the quadrupole is very weak relative to pressure; (ii) the eccentricity profile is largest near the inner edge of the disk, and falls exponentially outwards; and (iii) the results from linear theory indeed agree with what is found in simulations. Understanding the development of eccentric modes in circumbinary disks is a crucial first step for understanding the long term (secular) exchange of eccentricity, angular momentum and mass between the binary and the gas. Potential applications include the search for a characteristic kinematic signature in disks around candidate binaries and precession-induced modulation of accretion over long timescales.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Pileups and Migration Rates for Planets in Low Mass Disks

We investigate how planets interact with viscous accretion disks, in the limit that the disk is sufficiently low mass that the planet migrates more slowly than the disk material. In that case, the disk's surface density profile is determined by the disk being in viscous steady state (VSS), while overflowing the planet's orbit. We compute the VSS profiles with 2D hydrodynamical simulations, and show that disk material piles up behind the planet, with the planet effectively acting as a leaky dam. Previous 2D hydrodynamical simulations missed the pileup effect because of incorrect boundary conditions, while previous 1D models greatly overpredicted the pileup due to the neglect of non-local deposition. Our simulations quantify the magnitude of the pileup for a variety of planet masses and disk viscosities. We also calculate theoretically the magnitude of the pileup for moderately deep gaps, showing good agreement with simulations. For very deep gaps, current theory is inadequate, and we show why and what must be understood better. The pileup is important for two reasons. First, it is observable in directly imaged protoplanetary disks, and hence can be used to diagnose the mass of a planet that causes it or the viscosity within the disk. And second, it determines the planet's migration rate. Our simulations determine a new Type-II migration rate (valid for low mass disks), and show how it connects continuously with the well-verified Type-I rate.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Dynamical constraints on the origin of hot and warm Jupiters with close friends

Gas giants orbiting their host star within the ice line are thought to have migrated to their current locations from farther out. Here we consider the origin and dynamical evolution of observed Jupiters, focusing on hot and warm Jupiters with outer friends. We show that the majority of the observed Jupiter pairs (twenty out of twenty-four) will be dynamically unstable if the inner planet was placed at >~1AU distance from the stellar host. This finding is at odds with formation theories that invoke the migration of such planets from semi-major axes >~1AU due to secular dynamical processes (e.g., secular chaos, Lidov-Kozai oscillations) coupled with tidal dissipation. In fact, the results of N-body integrations show that the evolution of dynamically unstable systems does not lead to tidal migration but rather to planet ejections and collisions with the host star. This and other arguments lead us to suggest that most of the observed planets with a companion could not have been transported from further out through secular migration processes. More generally, by using a combination of numerical and analytic techniques we show that the high-e Lidov-Kozai migration scenario can only account for less than 10% of all gas giants observed between 0.1-1 AU. Simulations of multi-planet systems support this result. Our study indicates that rather than starting on highly eccentric orbits with orbital periods above one year, these "warm" Jupiters are more likely to have reached the region where they are observed today without having experienced significant tidal dissipation.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Secular dynamics of multiplanet systems: implications for the formation of hot and warm Jupiters via high-eccentricity migration

Hot Jupiters (HJs) are Jupiter-like planets that reside very closely to their host star, within $\sim 0.1\,\mathrm{AU}$. Their formation is not well understood. It is generally believed that they cannot have formed in situ, implying that some form of migration must have occurred after their initial formation. We study the production of HJs through secular evolution in multiplanet systems with three to five planets. In this variant of high-$e$ migration, the eccentricity of the orbit of the innermost planet is excited on secular time-scales, triggering orbital migration due to tidal dissipation. We use a secular dynamics code and carry out a population synthesis study. We find that HJs are only produced if the viscous time-scale is short ($\approx 0.014$ yr). In contrast, in up to $\approx 0.3$ of systems, the innermost planet is tidally disrupted. The orbital period distribution is peaked around 5 d, consistent with observations. The median HJ mass is $1\,M_\mathrm{J}$ with a maximum of $\approx 2 \, M_\mathrm{J}$, similar to observed HJs. Approximately 0.1 of the HJs have retrograde orbits with respect to the stellar spin. We do not find any warm Jupiters in our simulations, i.e. planets with semimajor axes between 0.1 and 1 AU.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Numerical and Analytical Modelling of Transit Time Variations

We develop and apply methods to extract planet masses and eccentricities from observed transit time variations (TTVs). First, we derive simple analytic expressions for the TTV that include the effects of both first- and second-order resonances. Second, we use N-body Markov chain Monte Carlo (MCMC) simulations, as well as the analytic formulae, to measure the masses and eccentricities of ten planets discovered by Kepler that have not previously been analyzed. Most of the ten planets have low densities. Using the analytic expressions to partially circumvent degeneracies, we measure small eccentricities of a few percent or less.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Densities and Eccentricities of 139 Kepler Planets from Transit Time Variations

We extract densities and eccentricities of 139 sub-Jovian planets by analyzing transit time variations (TTVs) obtained by the Kepler mission through Quarter 12. We partially circumvent the degeneracies that plague TTV inversion with the help of an analytical formula for the TTV. From the observed TTV phases, we find that most of these planets have eccentricities of order a few percent. More precisely, the r.m.s. eccentricity is 0.018^{+0.005}_{-0.004}, and planets smaller than 2.5R_\earth are around twice as eccentric as those bigger than 2.5R_\earth. We also find a best-fit density-radius relationship ρ~3 g/cm^3 \times (R/3R_\earth)^{-2.3} for the 56 planets that likely have small eccentricity and hence small statistical correction to their masses. Many planets larger than 2.5R_\earth are less dense than water, implying that their radii are largely set by a massive hydrogen atmosphere.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Theory and simulations of rotating convection

We study thermal convection in a rotating fluid in order to better understand the properties of convection zones in rotating stars and planets. We first derive mixing-length theory for rapidly-rotating convection, arriving at the results of Stevenson (1979) via simple physical arguments. The theory predicts the properties of convection as a function of the imposed heat flux and rotation rate, independent of microscopic diffusivities. In particular, it predicts the mean temperature gradient; the rms velocity and temperature fluctuations; and the size of the eddies that dominate heat transport. We test all of these predictions with high resolution three-dimensional hydrodynamical simulations of Boussinesq convection in a Cartesian box. The results agree remarkably well with the theory across more than two orders of magnitude in rotation rate. For example, the temperature gradient is predicted to scale as the rotation rate to the 4/5th power at fixed flux, and the simulations yield $0.75\pm 0.06$. We conclude that the mixing length theory is a solid foundation for understanding the properties of convection zones in rotating stars and planets.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

After Runaway: The Trans-Hill Stage of Planetesimal Growth

When planetesimals begin to grow by coagulation, they enter an epoch of runaway, during which the biggest bodies grow faster than all the others. The questions of how runaway ends and what comes next have not been answered satisfactorily. Here we show that runaway is followed by a `trans-hill stage' that commences once the bodies become trans-hill, i.e. once the Hill velocity of the bodies that dominate viscous stirring matches the random speed of the small bodies they accrete. Subsequently, the small bodies' speed grows in lockstep with the big bodies' sizes, such that the bodies remain in the trans-hill state. Trans-hill growth is crucial for determining the efficiency of growing big bodies, as well as their growth timescale and size spectrum. We work out the properties of trans-hill growth analytically and confirm these numerically. Trans-hill growth has two sub-stages. In the earlier one, collisional cooling is irrelevant, in which case the efficiency of forming big bodies remains very low (0.1% in the Kuiper belt) and their mass spectrum is flat. This explains results from previous coagulation simulations for both the Kuiper belt and the asteroid belt. The second sub-stage commences when small bodies begin to collide with one another. Collisional cooling then controls the evolution, in which case the efficiency of forming big bodies rises, and their size spectrum becomes more top-heavy. Trans-hill growth terminates in one of two ways, depending on parameters. First, mutual accretion of big bodies can become significant and conglomeration proceeds until half the total mass is converted into big bodies. This mode of growth may explain the size distributions of minor bodies in the Solar System, and is explored in forthcoming work. Second, if big bodies become separated by their Hill radius, oligarchy commences. This mode likely precedes the formation of fully-fledged planets.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Nonlinear evolution of the elliptical instability in the presence of weak magnetic fields

We investigate whether the elliptical instability is important for tidal dissipation in gaseous planets and stars. In a companion paper, we found that the conventional elliptical instability results in insufficient dissipation because it produces long-lived vortices that then quench further instability. Here, we study whether the addition of a magnetic field prevents those vortices from forming, and hence leads to enhanced dissipation. We present results from magnetohydrodynamical simulations that evolve the elliptical instability in a local patch of a rotating planet or star, in the presence of a weak magnetic field. We find that magnetic fields do indeed prevent vortices from forming, and hence greatly enhance the steady state dissipation rate. In addition, the resulting turbulence acts as a small-scale dynamo, amplifying the initially weak field. The inferred tidal dissipation is potentially important at short orbital periods. For example, it can circularise hot Jupiters with orbital periods shorter than 2.5 days, and synchronise their spins with their orbits out to 6 days. However, it appears unable to account for the hot Jupiters that appear to have been circularised out to six to ten day orbital periods. It also cannot account for the inferred circularisation of many close binary stars.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Nonlinear evolution of the tidal elliptical instability in gaseous planets and stars

Tidally distorted rotating stars and gaseous planets are subject to a well-known linear fluid instability -- the elliptical instability. It has been proposed that this instability might drive enough energy dissipation to solve the long-standing problem of the origin of tidal dissipation in stars and planets. But the nonlinear outcome of the elliptical instability has yet to be investigated in the parameter regime of interest, and the resulting turbulent energy dissipation has not yet been quantified. We do so by performing three dimensional hydrodynamical simulations of a small patch of a tidally deformed fluid planet or star subject to the elliptical instability. We show that when the tidal deformation is weak, the nonlinear outcome of the instability leads to the formation of long-lived columnar vortices aligned with the axis of rotation. These vortices shut off the elliptical instability, and the net result is insufficient energy dissipation to account for tidal dissipation. However, further work is required to account for effects neglected here, including magnetic fields, turbulent convection, and realistic boundary conditions.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Secular Dynamics in Hierarchical Three-Body Systems

The secular approximation for the evolution of hierarchical triple configurations has proven to be very useful in many astrophysical contexts, from planetary to triple-star systems. In this approximation the orbits may change shape and orientation, on time scales longer than the orbital time scales, but the semi major axes are constant. For example, for highly inclined triple systems, the Kozai-Lidov mechanism can produce large-amplitude oscillations of the eccentricities and inclinations. Here we revisit the secular dynamics of hierarchical triple systems. We derive the secular evolution equations to octupole order in Hamiltonian perturbation theory. Our derivation corrects an error in some previous treatments of the problem that implicitly assumed a conservation of the z-component of the angular momentum of the inner orbit (i.e., parallel to the total angular momentum of the system). Already to quadrupole order, our results show new behaviors including the possibility for a system to oscillate from prograde to retrograde orbits. At the octupole order, for an eccentric outer orbit, the inner orbit can reach extremely high eccentricities and undergo chaotic flips in its orientation. We discuss applications to a variety of astrophysical systems, from stellar triples to merging compact binaries and planetary systems. Our results agree with those of previous studies done to quadrupole order only in the limit in which one of the inner two bodies is a massless test particle and the outer orbit is circular;our results agree with previous studies at octupole order for the eccentricity evolution, but not for the inclination evolution.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Adding particle collisions to the formation of asteroids and Kuiper belt objects via streaming instabilities

Modelling the formation of super-km-sized planetesimals by gravitational collapse of regions overdense in small particles requires numerical algorithms capable of handling simultaneously hydrodynamics, particle dynamics and particle collisions. While the initial phases of radial contraction are dictated by drag forces and gravity, particle collisions become gradually more significant as filaments contract beyond Roche density. Here we present a new numerical algorithm for treating momentum and energy exchange in collisions between numerical superparticles representing a high number of physical particles. We adopt a Monte Carlo approach where superparticle pairs in a grid cell collide statistically on the physical collision time-scale. Collisions occur by enlarging particles until they touch and solving for the collision outcome, accounting for energy dissipation in inelastic collisions. We demonstrate that superparticle collisions can be consistently implemented at a modest computational cost. In protoplanetary disc turbulence driven by the streaming instability, we argue that the relative Keplerian shear velocity should be subtracted during the collision calculation. If it is not subtracted, density inhomogeneities are too rapidly diffused away, as bloated particles exaggerate collision speeds. Local particle densities reach several thousand times the mid-plane gas density. We find efficient formation of gravitationally bound clumps, with a range of masses corresponding to contracted radii from 100 to 400 km when applied to the asteroid belt and 150 to 730 km when applied to the Kuiper belt, extrapolated using a constant self-gravity parameter. The smaller planetesimals are not observed at low resolution, but the masses of the largest planetesimals are relatively independent of resolution and treatment of collisions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Hot Jupiters from Secular Planet--Planet Interactions

About 25 per cent of `hot Jupiters' (extrasolar Jovian-mass planets with close-in orbits) are actually orbiting counter to the spin direction of the star. Perturbations from a distant binary star companion can produce high inclinations, but cannot explain orbits that are retrograde with respect to the total angular momentum of the system. Such orbits in a stellar context can be produced through secular (that is, long term) perturbations in hierarchical triple-star systems. Here we report a similar analysis of planetary bodies, including both octupole-order effects and tidal friction, and find that we can produce hot Jupiters in orbits that are retrograde with respect to the total angular momentum. With distant stellar mass perturbers, such an outcome is not possible. With planetary perturbers, the inner orbit's angular momentum component parallel to the total angular momentum need not be constant. In fact, as we show here, it can even change sign, leading to a retrograde orbit. A brief excursion to very high eccentricity during the chaotic evolution of the inner orbit allows planet-star tidal interactions to rapidly circularize that orbit, decoupling the planets and forming a retrograde hot Jupiter.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Self-Similar Solutions of Triaxial Dark Matter Halos

We investigate the collapse and internal structure of dark matter halos. We consider halo formation from initially scale-free perturbations, for which gravitational collapse is self-similar. Fillmore and Goldreich (1984) and Bertschinger (1985) solved the one dimensional (i.e. spherically symmetric) case. We generalize their results by formulating the three dimensional self-similar equations. We solve the equations numerically and analyze the similarity solutions in detail, focusing on the internal density profiles of the collapsed halos. By decomposing the total density into subprofiles of particles that collapse coevally, we identify two effects as the main determinants of the internal density structure of halos: adiabatic contraction and the shape of a subprofile shortly after collapse; the latter largely reflects the triaxiality of the subprofile. We develop a simple model that describes the results of our 3D simulations. In a companion paper, we apply this model to more realistic cosmological fluctuations, and thereby explain the origin of the nearly universal (NFW-like) density profiles found in N-body simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Eccentric Kozai Mechanism for a Test Particle

We study the dynamical evolution of a test particle that orbits a star in the presence of an exterior massive planet, considering octupole-order secular interactions. In the standard Kozai mechanism (SKM), the planet's orbit is circular, and so the particle conserves vertical angular momentum. As a result, the particle's orbit oscillates periodically, exchanging eccentricity for inclination. However, when the planet's orbit is eccentric, the particle's vertical angular momentum varies and its Kozai oscillations are modulated on longer timescales---we call this the eccentric Kozai mechanism (EKM). The EKM can lead to behavior that is dramatically different from the SKM. In particular, the particle's orbit can flip from prograde to retrograde and back again, and it can reach arbitrarily high eccentricities given enough time. We map out the conditions under which this dramatic behavior (flipping and extreme eccentricities) occurs, and show that when the planet's eccentricity is sufficiently high, it occurs quite generically. For example, when the planet's eccentricity exceeds a few percent of the ratio of semi-major axes (outer to inner), around half of randomly oriented test particle orbits will flip and reach extreme eccentricities. The SKM has often been invoked for bringing pairs of astronomical bodies (star-star, planet-star, compact-object pairs) close together. Including the effect of the EKM will enhance the rate at which such matchmaking occurs.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Origin of Dark Matter Halo Profiles

A longstanding puzzle of fundamental importance in modern cosmology has been the origin of the nearly universal density profiles of dark matter halos found in N-body simulations -- the so-called NFW profile. We show how this behavior may be understood, simply, by applying adiabatic contraction to peaks of Gaussian random fields. We argue that dynamical friction acts to reduce enormously the effect of random scatter in the properties of initial peaks, providing a key simplification. We compare our model predictions with results of the ultra-high resolution Via Lactea-II N-body simulation, and find superb agreement. We show how our model may be used to predict the distribution of halo properties like concentration. Our results suggest that many of the basic properties of halo structure may be understood using extremely simple physics.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

The Intrinsic Size of Sagittarius A* from 0.35 cm to 6 cm

We present new high-resolution observations of Sagittarius A* at wavelengths of 17.4 to 23.8 cm with the Very Large Array in A configuration with the Pie Town Very Long Baseline Array antenna. We use the measured sizes to calibrate the interstellar scattering law and find that the major axis size of the scattering law is smaller by ~6% than previous estimates. Using the new scattering law, we are able to determine the intrinsic size of Sgr A* at wavelengths from 0.35 cm to 6 cm using existing results from the VLBA. The new law increases the intrinsic size at 0.7 cm by ~20% and <5% at 0.35 cm. The intrinsic size is 13^{+7}_{-3} Schwarzschild radii at 0.35 cm and is proportional to lambda^gamma, where gamma is in the range 1.3 to 1.7.

0 citations0 reviewsNo code linkedOpen access