Tidal Evolution of Low-mass Kepler Planets

The planets with a radius < 4 R_⊕ observed by the Kepler mission exhibit a unique feature, and propose a challenge for current planetary formation models. The tidal effect between a planet and its host star plays an essential role in reconfiguring the final orbits of the short-period planets. In this work, based on various initial Rayleigh distributions of the orbital elements, the final semi-major axis distributions of the planets with a radius < 4 R_⊕ after suffering tidal evolutions are investigated. Our simulations have qualitatively revealed some statistical properties: the semi-major axis and its peak value all increase with the increase of the initial semi-major axis and eccentricity. For the case that the initial mean semi-major axis is less than 0.1 au and the mean eccentricity is larger than 0.25, the results of numerical simulation are approximately consistent with the observation. In addition, the effects of other parameters, such as the tidal dissipation coefficient, stellar mass and planetary mass, etc., on the final semi-major axis distribution after tidal evolution are all relatively small. Based on the simulation results, we have tried to find some clues for the formation mechanism of low-mass planets. We speculate that these low-mass planets probably form in the far place of protoplanetary disk with a moderate eccentricity via the type I migration, and it is also possible to form in situ.     

Slow modes in stellar systems with nearly harmonic potentials: I. Spoke approximation,radial orbit instability

Using a consistent perturbation theory for collisionless disk-like and spherical star clusters, we construct a theory of slow modes for systems having an extended central region with a nearly harmonic potential due to the presence of a fairly homogeneous (on the scales of the stellar system) heavy, dynamically passive halo. In such systems, the stellar orbits are slowly precessing, centrally symmetric ellipses (2: 1 orbits). Depending on the density distribution in the system and the degree of halo inhomogeneity, the orbit precession can be both prograde and retrograde, in contrast to systems with 1: 1 elliptical orbits where the precession is unequivocally retrograde. In the first paper, we show that in the case where at least some of the orbits have a prograde precession and the stellar distribution function is a decreasing function of angular momentum, an instability that turns into the well-known radial orbit instability in the limit of low angular momenta can develop in the system. We also explore the question of whether the so-called spoke approximation, a simplified version of the slow mode approximation, is applicable for investigating the instability of stellar systems with highly elongated orbits. Highly elongated orbits in clusters with nonsingular gravitational potentials are known to be also slowly precessing 2: 1 ellipses. This explains the attempts to use the spoke approximation in finding the spectrum of slow modes with frequencies of the order of the orbit precession rate. We show that, in contrast to the previously accepted view, the dependence of the precession rate on angular momentum can differ significantly from a linear one even in a narrow range of variation of the distribution function in angular momentum. Nevertheless, using a proper precession curve in the spoke approximation allows us to partially “rehabilitate” the spoke approach, i.e., to correctly determine the instability growth rate, at least in the principal (O(α _T^−1/2) order of the perturbation theory in dimensionless small parameter α _T, which characterizes the width of the distribution function in angular momentum near radial orbits.     

Tidal locking of habitable exoplanets

Potentially habitable planets can orbit close enough to their host star that the differential gravity across their diameters can produce an elongated shape. Frictional forces inside the planet prevent the bulges from aligning perfectly with the host star and result in torques that alter the planet’s rotational angular momentum. Eventually the tidal torques fix the rotation rate at a specific frequency, a process called tidal locking. Tidally locked planets on circular orbits will rotate synchronously, but those on eccentric orbits will either librate or rotate super-synchronously. Although these features of tidal theory are well known, a systematic survey of the rotational evolution of potentially habitable exoplanets using classic equilibrium tide theories has not been undertaken. I calculate how habitable planets evolve under two commonly used models and find, for example, that one model predicts that the Earth’s rotation rate would have synchronized after 4.5 Gyr if its initial rotation period was 3 days, it had no satellites, and it always maintained the modern Earth’s tidal properties. Lower mass stellar hosts will induce stronger tidal effects on potentially habitable planets, and tidal locking is possible for most planets in the habitable zones of GKM dwarf stars. For fast-rotating planets, both models predict eccentricity growth and that circularization can only occur once the rotational frequency is similar to the orbital frequency. The orbits of potentially habitable planets of very late M dwarfs ( Open image in new window ) are very likely to be circularized within 1 Gyr, and hence, those planets will be synchronous rotators. Proxima b is almost assuredly tidally locked, but its orbit may not have circularized yet, so the planet could be rotating super-synchronously today. The evolution of the isolated and potentially habitable Kepler planet candidates is computed and about half could be tidally locked. Finally, projected TESS planets are simulated over a wide range of assumptions, and the vast majority of potentially habitable cases are found to tidally lock within 1 Gyr. These results suggest that the process of tidal locking is a major factor in the evolution of most of the potentially habitable exoplanets to be discovered in the near future.     

Interaction of free-floating planets with a star–planet pair

The recent discovery of free-floating planets and their theoretical interpretation as celestial bodies, either condensed independently or ejected from parent stars in tight clusters, introduced an intriguing possibility. Namely, that some exoplanets are not condensed from the protoplanetary disk of their parent star. In this novel scenario a free-floating planet interacts with an already existing planetary system, created in a tight cluster, and is captured as a new planet. In the present work we study this interaction process by integrating trajectories of planet-sized bodies, which encounter a binary system consisting of a Jupiter-sized planet revolving around a Sun-like star. To simplify the problem we assume coplanar orbits for the bound and the free-floating planet and an initially parabolic orbit for the free-floating planet. By calculating the uncertainty exponent, a quantity that measures the dependence of the final state of the system on small changes of the initial conditions, we show that the interaction process is a fractal classical scattering. The uncertainty exponent is in the range (0.2–0.3) and is a decreasing function of time. In this way we see that the statistical approach we follow to tackle the problem is justified. The possible final outcomes of this interaction are only four, namely flyby, planet exchange, capture or disruption. We give the probability of each outcome as a function of the incoming planet’s mass. We find that the probability of exchange or capture (in prograde as well as retrograde orbits and for very long times) is non-negligible, a fact that might explain the possible future observations of planetary systems with orbits that are either retrograde (see e.g. Queloz et?al. Astron. Astrophys. 417, L1, 2010) or tight and highly eccentric.     

The tidal stripping of satellites

We present an improved analytic calculation for the tidal radius of satellites and test our results against N -body simulations.
The tidal radius in general depends upon four factors: the potential of the host galaxy, the potential of the satellite, the orbit of the satellite and the orbit of the star within the satellite . We demonstrate that this last point is critical and suggest using three tidal radii to cover the range of orbits of stars within the satellite. In this way we show explicitly that prograde star orbits will be more easily stripped than radial orbits; while radial orbits are more easily stripped than retrograde ones. This result has previously been established by several authors numerically, but can now be understood analytically. For point mass, power-law (which includes the isothermal sphere), and a restricted class of split power-law potentials our solution is fully analytic. For more general potentials, we provide an equation which may be rapidly solved numerically.
Over short times (≲1–2 Gyr ∼1 satellite orbit), we find excellent agreement between our analytic and numerical models. Over longer times, star orbits within the satellite are transformed by the tidal field of the host galaxy. In a Hubble time, this causes a convergence of the three limiting tidal radii towards the prograde stripping radius. Beyond the prograde stripping radius, the velocity dispersion will be tangentially anisotropic.     

Effect of the rotation and tidal dissipation history of stars on the evolution of close-in planets

Since 20 years, a large population of close-in planets orbiting various classes of low-mass stars (from M-type to A-type stars) has been discovered. In such systems, the dissipation of the kinetic energy of tidal flows in the host star may modify its rotational evolution and shape the orbital architecture of the surrounding planetary system. In this context, recent observational and theoretical works demonstrated that the amplitude of this dissipation can vary over several orders of magnitude as a function of stellar mass, age and rotation. In addition, stellar spin-up occurring during the Pre-Main-Sequence (PMS) phase because of the contraction of stars and their spin-down because of the torque applied by magnetized stellar winds strongly impact angular momentum exchanges within star–planet systems. Therefore, it is now necessary to take into account the structural and rotational evolution of stars when studying the orbital evolution of close-in planets. At the same time, the presence of planets may modify the rotational dynamics of the host stars and as a consequence their evolution, magnetic activity and mixing. In this work, we present the first study of the dynamics of close-in planets of various masses orbiting low-mass stars (from \(0.6~M_\odot \) to \(1.2~M_\odot \)) where we compute the simultaneous evolution of the star’s structure, rotation and tidal dissipation in its external convective envelope. We demonstrate that tidal friction due to the stellar dynamical tide, i.e. tidal inertial waves excited in the convection zone, can be larger by several orders of magnitude than the one of the equilibrium tide currently used in Celestial Mechanics, especially during the PMS phase. Moreover, because of this stronger tidal friction in the star, the orbital migration of the planet is now more pronounced and depends more on the stellar mass, rotation and age. This would very weakly affect the planets in the habitable zone because they are located at orbital distances such that stellar tide-induced migration happens on very long timescales. We also demonstrate that the rotational evolution of host stars is only weakly affected by the presence of planets except for massive companions.     

Tidal pattern instabilities on multi-moon planets

The equilibrium tide-generating forces in the lunar orbital plane of a planet of radius R are calculated for the case of N moons of mass Mi orbiting the planet at instantaneous polar coordinates (Di, αi). For the case of a single moon, there are only two high tides. For the case of two moons, it is found that there can exist a critical lunar orbital distance at which two high tides become unstable with respect to formation of three high tides. Bifurcation diagrams are presented which depict how the angular positions of the high and low tides on the planet vary with the lunar distances and lunar separation angle. Tidal stability diagrams, which illustrate the stability regions for various tidal patterns as a function of lunar distances and lunar separation angle, are presented for various values of D₂/D₁ and M₂/M₁. Generally speaking, the aforementioned tidal instability, and hence the propensity for formation of three high tides on a two-moon planet, exists over a significant range of lunar distances and separation angles provided that M₂/M₁∼³(D₂/D₁). For the case of N>2 moons, the tidal stability diagram becomes more complex, revealing a diversity of potential tidal patterns.     

Change in the orbital period of a binary system due to dynamical tides for main-sequence stars

We investigate the change in the orbital period of a binary system due to dynamical tides by taking into account the evolution of a main-sequence star. Three stars with masses of one, one and a half, and two solar masses are considered. A star of one solar mass at lifetimes t = 4.57 × 10⁹ yr closely corresponds to our Sun. We show that a planet of one Jupiter mass revolving around a star of one solar mass will fall onto the star in the main-sequence lifetime of the star due to dynamical tides if the initial orbital period of the planet is less than P _orb ≈ 2.8 days. Planets of one Jupiter mass with an orbital period P _orb ≈ 2 days or shorter will fall onto a star of one and a half and two solar masses in the mainsequence lifetime of the star.     

Terrestrial planet formation surrounding close binary stars

Most stars reside in binary/multiple star systems; however, previous models of planet formation have studied growth of bodies orbiting an isolated single star. Disk material has been observed around both components of some young close binary star systems. Additionally, it has been shown that if planets form at the right places within such disks, they can remain dynamically stable for very long times. Herein, we numerically simulate the late stages of terrestrial planet growth in circumbinary disks around ‘close’ binary star systems with stellar separations 0.05 AU?aB?0.4 AU and binary eccentricities 0?eB?0.8. In each simulation, the sum of the masses of the two stars is 1 M_⊙, and giant planets are included. The initial disk of planetary embryos is the same as that used for simulating the late stages of terrestrial planet formation within our Solar System by Chambers [Chambers, J.E., 2001. Icarus 152, 205-224], and around each individual component of the α Centauri AB binary star system by Quintana et al. [Quintana, E.V., Lissauer, J.J., Chambers, J.E., Duncan, M.J., 2002. Astrophys. J. 576, 982-996]. Multiple simulations are performed for each binary star system under study, and our results are statistically compared to a set of planet formation simulations in the Sun-Jupiter-Saturn system that begin with essentially the same initial disk of protoplanets. The planetary systems formed around binaries with apastron distances QB≡aB(1+eB)?0.2 AU are very similar to those around single stars, whereas those with larger maximum separations tend to be sparcer, with fewer planets, especially interior to 1 AU. We also provide formulae that can be used to scale results of planetary accretion simulations to various systems with different total stellar mass, disk sizes, and planetesimal masses and densities.     

Secular and tidal evolution of circumbinary systems

We investigate the secular dynamics of three-body circumbinary systems under the effect of tides. We use the octupolar non-restricted approximation for the orbital interactions, general relativity corrections, the quadrupolar approximation for the spins, and the viscous linear model for tides. We derive the averaged equations of motion in a simplified vectorial formalism, which is suitable to model the long-term evolution of a wide variety of circumbinary systems in very eccentric and inclined orbits. In particular, this vectorial approach can be used to derive constraints for tidal migration, capture in Cassini states, and stellar spin–orbit misalignment. We show that circumbinary planets with initial arbitrary orbital inclination can become coplanar through a secular resonance between the precession of the orbit and the precession of the spin of one of the stars. We also show that circumbinary systems for which the pericenter of the inner orbit is initially in libration present chaotic motion for the spins and for the eccentricity of the outer orbit. Because our model is valid for the non-restricted problem, it can also be applied to any three-body hierarchical system such as star–planet–satellite systems and triple stellar systems.     

Tides and angular momentum redistribution inside low-mass stars hosting planets: a first dynamical model

We introduce a general mathematical framework to model the internal transport of angular momentum in a star hosting a close-in planetary/stellar companion. By assuming that the tidal and rotational distortions are small and that the deposit/extraction of angular momentum induced by stellar winds and tidal torques are redistributed solely by an effective eddy-viscosity that depends on the radial coordinate, we can formulate the model in a completely analytic way. It allows us to compute simultaneously the evolution of the orbit of the companion and of the spin and the radial differential rotation of the star. An illustrative application to the case of an F-type main-sequence star hosting a hot Jupiter is presented. The general relevance of our model to test more sophisticated numerical dynamical models and to study the internal rotation profile of exoplanet hosts, submitted to the combined effects of tides and stellar winds, by means of asteroseismology are discussed.     

MHD simulation scenarios of the stellar wind interaction with Hot Jupiter magnetospheres

We present three-dimensional numerical simulations of the interaction between a Hot Jupiter and the stellar wind plasma of its host star in the framework of resistive magnetohydrodynamics (MHD). In a first step, we investigate the numerical realization of the plasma flow around the planet and the planetary magnetic field using a simplified model, before we simulate more realistic scenarios on the basis of the stellar wind model by Weber and Davis. A main goal is to understand the magnetic interaction between star and planet. In analogy to the well-known Jupiter Io scenario, we study the development of a magnetic field-aligned current system in different parameter regimes.     

Effect of radiation on the stability of a retrograde particle orbit in different stellar systems

We have numerically investigated the stability of retrograde orbits/trajectories around Jupiter and the smaller of the primaries in binary systems RW-Monocerotis (RW-Mon) and Krüger-60 in the presence of radiation. A trajectory is considered as stable if it remains around the smaller mass for at least few hundred binary periods. In case of circular binary orbit, we find that the third order resonance provides the basis for reduction of stability region of retrograde motion of particle in RW-Mon and Sun-Jupiter system both in the presence and absence of radiation. Considering finite ellipticity in Sun-Jupiter system we find that for distant retrograde orbits, radiation from the Sun increases the width of the stable region and covers a significant portion of the region obtained in the absence of solar radiation. Further, due to solar radiation pressure, the stable region in the neighborhood of Jupiter has been found to shift much below the characteristic asymptotic line for the periodic retrograde orbits. In case of Krüger-60 we observe the distant retrograde orbits around the smaller of the primaries get affected considerably with increase in radiation parameter β₁. Further the range of velocities for which stable motion may persist narrows down for distant retrograde orbits in this system.     

Dust on the Outskirts of the Jovian System

The outer region of the jovian system between ∼50 and 300 jovian radii from the planet is found to be the host of a previously unknown dust population. We used the data from the dust detector aboard the Galileo spacecraft collected from December 1995 to April 2001 during Galileo's numerous traverses of the outer jovian system. Analyzing the ion amplitudes, calibrated masses and speeds of grains, and impact directions, we found about 100 individual events fully compatible with impacts of grains moving around Jupiter in bound orbits. These grains have moderate eccentricities and a wide range of inclinations—from prograde to retrograde ones. The radial number density profile of the micrometer-sized dust is nearly flat between about 50 and 300 jovian radii. The absolute number density level (∼10 km⁻³ with a factor of 2 or 3 uncertainty) surpasses by an order of magnitude that of the interplanetary background. We identify the sources of the bound grains with outer irregular satellites of Jupiter. Six outer tiny moons are orbiting the planet in prograde and fourteen in retrograde orbits. These moons are subject to continuous bombardment by interplanetary micrometeoroids. Hypervelocity impacts create ejecta, nearly all of which get injected into circumjovian space. Our analytic and numerical study of the ejecta dynamics shows that micrometer-sized particles from both satellite families, although strongly perturbed by solar tidal gravity and radiation pressure, would stay in bound orbits for hundreds of thousands of years as do a fraction of smaller grains, several tenths of a micrometer in radius, ejected from the prograde moons. Different-sized ejecta remain confined to spheroidal clouds embracing the orbits of the parent moons, with appreciable asymmetries created by the radiation pressure and solar gravity perturbations. Spatial location of the impacts, mass distribution, speeds, orbital inclinations, and number density of dust derived from the data are all consistent with the dynamical model.     

Five new and three improved mutual orbits of transneptunian binaries

We present three improved and five new mutual orbits of transneptunian binary systems (58534) Logos-Zoe, (66652) Borasisi-Pabu, (88611) Teharonhiawako-Sawiskera, (123509) 2000 WK₁₈₃, (149780) Altjira, 2001 QY₂₉₇, 2003 QW₁₁₁, and 2003 QY₉₀ based on Hubble Space Telescope and Keck II laser guide star adaptive optics observations. Combining the five new orbit solutions with 17 previously known orbits yields a sample of 22 mutual orbits for which the period P, semimajor axis a, and eccentricity e have been determined. These orbits have mutual periods ranging from 5 to over 800 days, semimajor axes ranging from 1600 to 37,000 km, eccentricities ranging from 0 to 0.8, and system masses ranging from 2 × 10¹⁷ to 2 × 10²² kg. Based on the relative brightnesses of primaries and secondaries, most of these systems consist of near equal-sized pairs, although a few of the most massive systems are more lopsided. The observed distribution of orbital properties suggests that the most loosely-bound transneptunian binary systems are only found on dynamically cold heliocentric orbits. Of the 22 known binary mutual orbits, orientation ambiguities are now resolved for 9, of which 7 are prograde and 2 are retrograde, consistent with a random distribution of orbital orientations, but not with models predicting a strong preference for retrograde orbits. To the extent that other perturbations are not dominant, the binary systems undergo Kozai oscillations of their eccentricities and inclinations with periods of the order of tens of thousands to millions of years, some with strikingly high amplitudes.     

Injection of Oort Cloud comets: the fundamental role of stellar perturbations

We present Monte Carlo simulations of the dynamical evolution of the Oort cloud over the age of the Solar System, using an initial sample of one million test comets without any cloning. Our model includes perturbations due to the Galactic tide (radial and vertical) and passing stars. We present the first detailed analysis of the injection mechanism into observable orbits by comparing the complete model with separate models for tidal and stellar perturbations alone. We find that a fundamental role for injecting comets from the region outside the loss cone (perihelion distance q > 15 AU) into observable orbits (q < 5 AU) is played by stellar perturbations. These act in synergy with the tide such that the total injection rate is significantly larger than the sum of the two separate rates. This synergy is as important during comet showers as during quiescent periods and concerns comets with both small and large semi-major axes. We propose different dynamical mechanisms to explain the synergies in the inner and outer parts of the Oort Cloud. We find that the filling of the observable part of the loss cone under normal conditions in the present-day Solar System rises from <1% for a < 20 000 AU to about 100% for a ? 100 000 AU.     

An analytical method to compute comet cloud formation efficiency and its application

A quick analytical method is presented for calculating comet cloud formation efficiency in the case of a single planet or multiple-planet system for planets that are not too eccentric (e _p ≲ 0.3). A method to calculate the fraction of comets that stay under the control of each planet is also presented, as well as a way to determine the efficiency in different star cluster environments. The location of the planet(s) in mass-semi-major axis space to form a comet cloud is constrained based on the conditions developed by Tremaine (1993) together with estimates of the likelyhood of passing comets between planets; and, in the case of a single, eccentric planet, the additional constraint that it is, by itself, able to accelerate material to relative encounter velocity U ~ 0.4 within the age of the stellar system without sweeping up the majority of the material beforehand. For a single planet, it turns out the efficiency is mainly a function of planetary mass and semi-major axis of the planet and density of the stellar environment. The theory has been applied to some extrasolar systems and compared to numerical simulations for both these systems and the Solar System, as well as a diffusion scheme based on the energy kick distribution of Everhart (Astron J 73:1039–1052, 1968). The analytic results are in good agreement with the simulations.     

Orbit Orientation of the Planet in <Emphasis Type="Italic">τ</Emphasis> Boo (HD 120136)

We determined the rotation period and the parameters of the global magnetic field of τ Boo. This allowed us to estimate the inclination of the rotational axis of the star to the line of sight, as well as to obtain estimates of the inclination of the planet’s orbital plane to the stellar equator.     

Tidal evolution of hierarchical and inclined systems

We investigate the dynamical evolution of hierarchical three-body systems under the effect of tides, when the ratio of the orbital semi-major axes is small and the mutual inclination is relatively large (greater than 20°). Using the quadrupolar non-restricted approximation for the gravitational interactions and the viscous linear model for tides, we derive the averaged equations of motion in a vectorial formalism which is suitable to model the long-term evolution of a large variety of exoplanetary systems in very eccentric and inclined orbits. In particular, it can be used to derive constraints for stellar spin-orbit misalignment, capture in Cassini states, tidal-Kozai migration, or damping of the mutual inclination. Because our model is valid for the non-restricted problem, it can be used to study systems of identical mass or for the outer restricted problem, such as the evolution of a planet around a binary of stars. Here, we apply our model to various situations in the HD 11964, HD 80606, and HD 98800 systems.     

Model bond albedos of extrasolar giant planets

The atmospheres of extrasolar giant planets are modeled with various effective temperatures and gravities, with and without clouds. Bond albedos are computed by calculating the ratio of the flux reflected by a planet (integrated over wavelength) to the total stellar flux incident on the planet. This quantity is useful for estimating the effective temperature and evolution of a planet. We find it is sensitive to the stellar type of the primary. For a 5 M_Jup planet the Bond albedo varies from 0.4 to 0.3 to 0.6 as the primary star varies from A5V to G2V to M2V in spectral type. It is relatively insensitive to the effective temperature and gravity for cloud-free planets. Water clouds increase the reflectivity of the planet in the red, which increases the Bond albedo. The Bond albedo increases by an order of magnitude for a 13 M_Jup planet with an M2V primary when water clouds are present. Silicate clouds, on the other hand, can either increase or decrease the Bond albedo, depending on whether there are many small grains (the former) or few large grains (the latter).