Close encounters of a rotating star with planets in parabolic orbits of varying inclination and the formation of hot Jupiters

In this paper we extend the theory of close encounters of a giant planet on a parabolic orbit with a central star developed in our previous work (Ivanov and Papaloizou in MNRAS 347:437, 2004; MNRAS 376:682, 2007) to include the effects of tides induced on the central star. Stellar rotation and orbits with arbitrary inclination to the stellar rotation axis are considered. We obtain results both from an analytic treatment that incorporates first order corrections to normal mode frequencies arising from stellar rotation and numerical treatments that are in satisfactory agreement over the parameter space of interest. These results are applied to the initial phase of the tidal circularisation problem. We find that both tides induced in the star and planet can lead to a significant decrease of the orbital semi-major axis for orbits having periastron distances smaller than 5?C6 stellar radii with tides in the star being much stronger for retrograde orbits compared to prograde orbits. Assuming that combined action of dynamic and quasi-static tides could lead to the total circularisation of orbits this corresponds to observed periods up to 4?C5 days. We use the simple Skumanich law to characterise the rotational history of the star supposing that the star has its rotational period equal to one month at the age of 5 Gyr. The strength of tidal interactions is characterised by circularisation time scale, t _ev , which is defined as a typical time scale of evolution of the planet??s semi-major axis due to tides. This is considered as a function of orbital period P _obs , which the planet obtains after the process of tidal circularisation has been completed. We find that the ratio of the initial circularisation time scales corresponding to prograde and retrograde orbits, respectively, is of order 1.5?C2 for a planet of one Jupiter mass having P _obs ~ 4 days. The ratio grows with the mass of the planet, being of order five for a five Jupiter mass planet with the same P _orb . Note, however, this result might change for more realistic stellar rotation histories. Thus, the effect of stellar rotation may provide a bias in the formation of planetary systems having planets on close orbits around their host stars, as a consequence of planet?Cplanet scattering, which favours systems with retrograde orbits. The results reported in the paper may also be applied to the problem of tidal capture of stars in young stellar clusters.     

Obliquity evolution of extrasolar terrestrial planets

We have investigated the obliquity evolution of terrestrial planets in habitable zones (at ∼1 AU) in extrasolar planetary systems, due to tidal interactions with their satellite and host star with wide varieties of satellite-to-planet mass ratio (m/Mp) and initial obliquity (γ₀), through numerical calculations and analytical arguments. The obliquity, the angle between planetary spin axis and its orbit normal, of a terrestrial planet is one of the key factors in determining the planetary surface environments. A recent scenario of terrestrial planet accretion implies that giant impacts of Mars-sized or larger bodies determine the planetary spin and form satellites. Since the giant impacts would be isotropic, tilted spins (sinγ₀∼1) are more likely to be produced than straight ones (sinγ₀∼0). The ratio m/Mp is dependent on the impact parameters and impactors' mass. However, most of previous studies on tidal evolution of the planet-satellite systems have focused on a particular case of the Earth-Moon systems in which m/Mp?0.0125 and γ₀∼10° or the two-body planar problem in which γ₀=0° and stellar torque is neglected. We numerically integrated the evolution of planetary spin and a satellite orbit with various m/Mp (from 0.0025 to 0.05) and γ₀ (from 0° to 180°), taking into account the stellar torques and precessional motions of the spin and the orbit. We start with the spin axis that almost coincides with the satellite orbit normal, assuming that the spin and the satellite are formed by one dominant impact. With initially straight spins, the evolution is similar to that of the Earth-Moon system. The satellite monotonically recedes from the planet until synchronous state between the spin period and the satellite orbital period is realized. The obliquity gradually increases initially but it starts decreasing down to zero as approaching the synchronous state. However, we have found that the evolution with initially tiled spins is completely different. The satellite's orbit migrates outward with almost constant obliquity until the orbit reaches the critical radius ∼10-20 planetary radii, but then the migration is reversed to inward one. At the reversal, the obliquity starts oscillation with large amplitude. The oscillation gradually ceases and the obliquity is reduced to ∼0° during the inward migration. The satellite eventually falls onto the planetary surface or it is captured at the synchronous state at several planetary radii. We found that the character change of precession about total angular momentum vector into that about the planetary orbit normal is responsible for the oscillation with large amplitude and the reversal of migration. With the results of numerical integration and analytical arguments, we divided the m/Mp-γ₀ space into the regions of the qualitatively different evolution. The peculiar tidal evolution with initially tiled spins give deep insights into dynamics of extrasolar planet-satellite systems and discussions of surface environments of the planets.     

The Effect of Tidal Interaction with a Gas Disk on Formation of Terrestrial Planets

We have performed N-body simulation on final accretion stage of terrestrial planets, including the effect of damping of eccentricity and inclination caused by tidal interaction with a remnant gas disk. As a result of runway and oligarchic accretion, about 20 Mars-sized protoplanets would be formed in nearly circular orbits with orbital separation of several to ten Hill radius. The orbits of the protoplanets would be eventually destabilized by long-term mutual gravity and/or secular resonance of giant gaseous planets. The protoplanets would coalesce with each other to form terrestrial planets through the orbital crossing. Previous N-body simulations, however, showed that the final eccentricities of planets are around 0.1, which are about 10 times higher than the present eccentricities of Earth and Venus. The obtained high eccentricities are the remnant of orbital crossing. We included the effect of eccentricity damping caused by gravitational interaction with disk gas as a drag force (“gravitational drag”) and carried out N-body simulation of accretion of protoplanets. We start with 15 protoplanets with 0.2M⊕ and integrate the orbits for 10⁷ years, which is consistent with the observationally inferred disk lifetime (in some runs, we start with 30 protoplanets with 0.1M⊕). In most runs, the damping time scale, which is equivalent to the strength of the drag force, is kept constant throughout each run in order to clarify the effects of the damping. We found that the planets' final mass, spatial distribution, and eccentricities depend on the damping time scale. If the damping time scale for a 0.2M⊕ mass planet at 1 AU is longer than 10⁸ years, planets grow to Earth's size, but the final eccentricities are too high as in gas-free cases. If it is shorter than 10⁶ years, the eccentricities of the protoplanets cannot be pumped up, resulting in not enough orbital crossing to make Earth-sized planets. Small planets with low eccentricities are formed with small orbital separation. On the other hand, if it is between 10⁶ and 10⁸ years, which may correspond to a mostly depleted disk (0.01-0.1% of surface density of the minimum mass model), some protoplanets can grow to about the size of Earth and Venus, and the eccentricities of such surviving planets can be diminished within the disk lifetime. Furthermore, in innermost and outermost regions in the same system, we often find planets with smaller size and larger eccentricities too, which could be analogous to Mars and Mercury. This is partly because the gravitational drag is less effective for smaller mass planets, and partly due to the “edge effect,” which means the innermost and outermost planets tend to remain without collision. We also carried out several runs with time-dependent drag force according to depletion of a gas disk. In these runs, we used exponential decay model with e-folding time of 3×10⁶ years. The orbits of protoplanets are stablized by the eccentricity damping in the early time. When disk surface density decays to ?1% of the minimum mass disk model, the damping force is no longer strong enough to inhibit the increase of the eccentricity by distant perturbations among protoplanets so that the orbital crossing starts. In this disk decay model, a gas disk with 10⁻⁴-10⁻³ times the minimum mass model still remains after the orbital crossing and accretional events, which is enough to damp the eccentricities of the Earth-sized planets to the order of 0.01. Using these results, we discuss a possible scenario for the last stage of terrestrial planet formation.     

Orbital stability of systems of closely-spaced planets

An investigation of the stability of systems of 1 M_⊕ (Earth-mass) bodies orbiting a Sun-like star has been conducted for virtual times reaching 10 billion years. For the majority of the tests, a symplectic integrator with a fixed timestep of between 1 and 10 days was employed; however, smaller timesteps and a Bulirsch-Stoer integrator were also selectively utilized to increase confidence in the results. In most cases, the planets were started on initially coplanar, circular orbits, and the longitudinal initial positions of neighboring planets were widely separated. The ratio of the semimajor axes of consecutive planets in each system was approximately uniform (so the spacing between consecutive planets increased slowly in terms of distance from the star). The stability time for a system was taken to be the time at which the orbits of two or more planets crossed. Our results show that, for a given class of system (e.g., three 1 M_⊕ planets), orbit crossing times vary with planetary spacing approximately as a power law over a wide range of separation in semimajor axis. Chaos tests indicate that deviations from this power law persist for changed initial longitudes and also for small but non-trivial changes in orbital spacing. We find that the stability time increases more rapidly at large initial orbital separations than the power-law dependence predicted from moderate initial orbital separations. Systems of five planets are less stable than systems of three planets for a specified semimajor axis spacing. Furthermore, systems of less massive planets can be packed more closely, being about as stable as 1 M_⊕ planets when the radial separation between planets is scaled using the mutual Hill radius. Finally, systems with retrograde planets can be packed substantially more closely than prograde systems with equal numbers of planets.     

The extreme physical properties of the CoRoT-7b super-Earth

The search for rocky exoplanets plays an important role in our quest for extra-terrestrial life. Here, we discuss the extreme physical properties possible for the first characterised rocky super-Earth, CoRoT-7b (R_pl = 1.58 ± 0.10 R_Earth, M_pl = 6.9 ± 1.2 M_Earth). It is extremely close to its star (a = 0.0171 AU = 4.48 R_st), with its spin and orbital rotation likely synchronised. The comparison of its location in the (M_pl, R_pl) plane with the predictions of planetary models for different compositions points to an Earth-like composition, even if the error bars of the measured quantities and the partial degeneracy of the models prevent a definitive conclusion. The proximity to its star provides an additional constraint on the model. It implies a high extreme-UV flux and particle wind, and the corresponding efficient erosion of the planetary atmosphere especially for volatile species including water. Consequently, we make the working hypothesis that the planet is rocky with no volatiles in its atmosphere, and derive the physical properties that result. As a consequence, the atmosphere is made of rocky vapours with a very low pressure (P ? 1.5 Pa), no cloud can be sustained, and no thermalisation of the planet is expected. The dayside is very hot (2474 ± 71 K at the sub-stellar point) while the nightside is very cold (50-75 K). The sub-stellar point is as hot as the tungsten filament of an incandescent bulb, resulting in the melting and distillation of silicate rocks and the formation of a lava ocean. These possible features of CoRoT-7b could be common to many small and hot planets, including the recently discovered Kepler-10b. They define a new class of objects that we propose to name “Lava-ocean planets”.     

Mass-radius curve for extrasolar Earth-like planets and ocean planets

By comparison with the Earth-like planets and the large icy satellites of the Solar System, one can model the internal structure of extrasolar planets. The input parameters are the composition of the star (Fe/Si and Mg/Si), the Mg content of the mantle (Mg# = Mg/[Mg + Fe]), the amount of H₂O and the total mass of the planet. Equation of State (EoS) of the different materials that are likely to be present within such planets have been obtained thanks to recent progress in high-pressure experiments. They are used to compute the planetary radius as a function of the total mass. Based on accretion models and data on planetary differentiation, the internal structure is likely to consist of an iron-rich core, a silicate mantle and an outer silicate crust resulting from magma formation in the mantle. The amount of H₂O and the surface temperature control the possibility for these planets to harbor an ocean. In preparation to the interpretation of the forthcoming data from the CNES led CoRoT (Convection Rotation and Transit) mission and from ground-based observations, this paper investigates the relationship between radius and mass. If H₂O is not an important component (less than 0.1%) of the total mass of the planet, then a relation (R/REarth)=ab(M/MEarth) is calculated with (a,b)=(1,0.306) and (a,b)=(1,0.274) for ¹⁰⁻²MEarth<M<MEarth and MEarth<M<10MEarth, respectively. Calculations for a planet that contains 50% H₂O suggest that the radius would be more than 25% larger than that based on the Earth-like model, with (a,b)=(1.258,0.302) for ¹⁰⁻²MEarth<M<MEarth and (a,b)=(1.262,0.275) for MEarth<M<10MEarth, respectively. For a surface temperature of 300 K, the thickness of the ocean varies from 150 to 50 km for planets 1 to 10 times the Earth's mass, respectively. Application of this algorithm to bodies of the Solar System provides not only a good fit to most terrestrial planets and large icy satellites, but also insights for discussing future observations of exoplanets.     

Stability of a planet in the HD 41004 binary system

The Hill stability criterion is applied to analyse the stability of a planet in the binary star system of HD 41004 AB, with the primary and secondary separated by 22 AU, and masses of 0.7 M_⊙ and 0.4 M_⊙, respectively. The primary hosts one planet in an S‐type orbit, and the secondary hosts a brown dwarf (18.64 M_J) on a relatively close orbit, 0.0177 AU, thereby forming another binary pair within this binary system. This star‐brown dwarf pair (HD 41004 B+Bb) is considered a single body during our numerical calculations, while the dynamics of the planet around the primary, HD 41004 Ab, is studied in different phase‐spaces. HD 41004 Ab is a 2.6 M_J planet orbiting at the distance of 1.7 AU with orbital eccentricity 0.39. For the purpose of this study, the system is reduced to a three‐body problem and is solved numerically as the elliptic restricted three‐body problem (ERTBP). The Hill stability function is used as a chaos indicator to configure and analyse the orbital stability of the planet, HD 41004 Ab. The indicator has been effective in measuring the planet's orbital perturbation due to the secondary star during its periastron passage. The calculated Hill stability time series of the planet for the coplanar case shows the stable and quasi‐periodic orbits for at least ten million years. For the reduced ERTBP the stability of the system is also studied for different values of planet's orbital inclination with the binary plane. Also, by recording the planet's ejection time from the system or collision time with a star during the integration period, stability of the system is analysed in a bigger phase‐space of the planet's orbital inclination, ≤ 90°, and its semimajor axis, 1.65–1.75 AU. Based on our analysis it is found that the system can maintain a stable configuration for the planet's orbital inclination as high as 65° relative to the binary plane. The results from the Hill stability criterion and the planet's dynamical lifetime map are found to be consistent with each other. (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Contribution of tidal dissipation to lunar thermal history

The possible contributions of tidal heating to lunar thermal history are investigated. Analytic determinations of tidal dissipation in a homogeneous, incompressible Moon and in a two-layer Moon with a soft core and rigid mantle are given as a function of position in the Moon and as a function of Earth-Moon separation. The most recent information on the historical values of the lunar obliquity is employed, and we present results for the constant values of orbital eccentricity of e = 0.0 and e = 0.055. For a simplified orbital evolution and a dissipation factor Q = 100, the total increase in the mean lunar temperature for the homogeneous case does not exceed several tens of degrees. For the two-layer models the local dissipation may be enhanced over that of the homogeneous Moon by a factor of 5 for a core radius of 0.5 lunar radii and by a factor of 100 for a core radius of 0.95 lunar radii. The corresponding factors for the total dissipation are 3 and 15 for the two values of core radii, respectively. We conclude that tidal contributions to lunar thermal history are probably not important. But under special circumstances the enhanced dissipation in a two-layer Moon could have led to a spectacular thermal event.     

How common are Earth-Moon planetary systems?

The Earth’s comparatively massive moon, formed via a giant impact on the proto-Earth, has played an important role in the development of life on our planet, both in the history and strength of the ocean tides and in stabilizing the chaotic spin of our planet. Here we show that massive moons orbiting terrestrial planets are not rare. A large set of simulations by Morishima et al. (Morishima, R., Stadel, J., Moore, B. [2010]. Icarus. 207, 517-535), where Earth-like planets in the habitable zone form, provides the raw simulation data for our study. We use limits on the collision parameters that may guarantee the formation of a circumplanetary disk after a protoplanet collision that could form a satellite and study the collision history and the long term evolution of the satellites qualitatively. In addition, we estimate and quantify the uncertainties in each step of our study. We find that giant impacts with the required energy and orbital parameters for producing a binary planetary system do occur with more than 1 in 12 terrestrial planets hosting a massive moon, with a low-end estimate of 1 in 45 and a high-end estimate of 1 in 4.     

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.     

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.     

Accretion of terrestrial planets from oligarchs in a turbulent disk

We have investigated the final accretion stage of terrestrial planets from Mars-mass protoplanets that formed through oligarchic growth in a disk comparable to the minimum mass solar nebula (MMSN), through N-body simulation including random torques exerted by disk turbulence due to Magneto-Rotational Instability. For the torques, we used the semi-analytical formula developed by Laughlin et al. [Laughlin, G., Steinacker, A., Adams, F.C., 2004. Astrophys. J. 608, 489-496]. The damping of orbital eccentricities (in all runs) and type-I migration (in some runs) due to the tidal interactions with disk gas is also included. Without any effect of disk gas, Earth-mass planets are formed in terrestrial planet regions in a disk comparable to MMSN but with too large orbital eccentricities to be consistent with the present eccentricities of Earth and Venus in our Solar System. With the eccentricity damping caused by the tidal interaction with a remnant gas disk, Earth-mass planets with eccentricities consistent with those of Earth and Venus are formed in a limited range of disk gas surface density (∼10⁻⁴ times MMSN). However, in this case, on average, too many (?6) planets remain in terrestrial planet regions, because the damping leads to isolation between the planets. We have carried out a series of N-body simulations including the random torques with different disk surface density and strength of turbulence. We found that the orbital eccentricities pumped up by the turbulent torques and associated random walks in semimajor axes tend to delay isolation of planets, resulting in more coagulation of planets. The eccentricities are still damped after planets become isolated. As a result, the number of final planets decreases with increase in strength of the turbulence, while Earth-mass planets with small eccentricities are still formed. In the case of relatively strong turbulence, the number of final planets are 4-5 at 0.5-2 AU, which is more consistent with Solar System, for relatively wide range of disk gas surface density (∼¹⁰⁻⁴-¹⁰⁻² times MMSN).     

There are ten,not eight planets

In 1946, E. Sevin postulated the global vibrations of the Sun with a period P ₀ = 1/9 day and a “wavelength” L ₀ = c × P ₀ = 19.24 AU and predicted the tenth planet at a mean distance of 4.0 × L ₀ ≈ 77.0 AU from the Sun (c is the speed of light). The global vibrations of the Sun, precisely with the period of 1/9 day, were actually detected in 1974. Recently, the largest Kuiper Bell object 2003 UB₃₁₃, or Eris, with an orbital semimajor axis ≈ 3.5 × L ₀ ≈ 67.5 AU was discovered. We adduce arguments for the status of Eris as our tenth planet: (i) the object is larger and farther from the Sun than Pluto and (ii) the semimajor axis of Eris agrees well with the sequence of planetary distances that follows from the resonance spectrum of the Solar system dimensions (with the scale L ₀ and for all 11 orbits, including those of Pluto, Eris, and the asteroid belt). We point to a mistake of the Prague (2006) IAU Assembly, which excluded Pluto from the family of planets by introducing a new, highly controversial class of objects—“dwarf planets.”     

Could we identify hot ocean-planets with CoRoT, Kepler and Doppler velocimetry?

Planets less massive than about 10 MEarth are expected to have no massive H-He atmosphere and a cometary composition (∼50% rocks, 50% water, by mass) provided they formed beyond the snowline of protoplanetary disks. Due to inward migration, such planets could be found at any distance between their formation site and the star. If migration stops within the habitable zone, this may produce a new kind of planets, called ocean-planets. Ocean-planets typically consist in a silicate core, surrounded by a thick ice mantle, itself covered by a 100 km-deep ocean. The possible existence of ocean-planets raises important astrobiological questions: Can life originate on such body, in the absence of continent and ocean-silicate interfaces? What would be the nature of the atmosphere and the geochemical cycles? In this work, we address the fate of hot ocean-planets produced when migration ends at a closer distance. In this case the liquid/gas interface can disappear, and the hot H₂O envelope is made of a supercritical fluid. Although we do not expect these bodies to harbor life, their detection and identification as water-rich planets would give us insight as to the abundance of hot and, by extrapolation, cool ocean-planets. The water reservoir of these planets seems to be weakly affected by gravitational escape, provided that they are located beyond some minimum distance, e.g. 0.04 AU for a 5-Earth-mass planet around a Sun-like star. The swelling of their water atmospheres by the high stellar flux is expected not to significantly increase the planets' radii. We have studied the possibility of detecting and characterizing these hot ocean-planets by measuring their mean densities using transit missions in space—CoRoT (CNES) and Kepler (NASA)—in combination with Doppler velocimetry from the ground—HARPS (ESO) and possible future instruments. We have determined the domain in the [stellar magnitude, orbital distance] plane where discrimination between ocean-planets and rocky planets is possible with these instruments. The brightest stars of the mission target lists and the planets closest to their stars are the most favorable cases. Full advantage of high precision photometry by CoRoT, and particularly Kepler, can be obtained only if a new generation of Doppler instruments is built.     

Spin states and climates of eccentric exoplanets

The known extrasolar planets exhibit a wide range of orbital eccentricities e. This has a profound influence on their rotations and climates. Because of tides in their interiors, mostly solid exoplanets are expected eventually to despin to a state of spin-orbit resonance, where the orbital period is some integer or half-integer times the rotation period. The most important of these resonances is the synchronous state, where the planet's spin period exactly equals its orbital period (like Earth's Moon, and indeed most of the regular satellites in the Solar System). Such planets seem doomed to roast on one side and freeze on the other. However, synchronous planets rock back and forth by an angle of ∼2Arcsine with respect to the sub-stellar point. For e=0.055 (as for the Moon), this optical libration amounts to only ∼6°; but for a synchronous planet with e=0.50, for example, it would rise to ∼59°. This greatly expands the temperate “twilight zone” near the terminator and considerably improves the planet's prospects for habitability. For e?0.72389, the optical libration exceeds 90°; for such planets, the sector of permanent night vanishes, while the sunniest region splits in two. Furthermore, the synchronous state is not the only possible spin resonance. For example, Mercury (with e≈0.206) has an orbital period exactly 1.5 times its rotation period. A terrestrial exoplanet with e=0.40, say, is liable to have an orbital period of 2.0, 2.5, or 3.0 times its spin period. The corresponding insolation patterns are generally complicated, and all different from the synchronous state. Yet these non-synchronous resonances also protect certain longitudes from the worst extremes of temperature and solar radiation, and improve the planet's habitability, compared to non-resonant rotation. These results also have implications for the direct detectability of extrasolar planets, and the interpretation of their thermal emissions.     

Exoplanet's Figure and Its Interior

Along with the development of the observing technology, the observation and study on the exoplanets’ oblateness and apsidal precession have achieved significant progress. The oblateness of an exoplanet is determined by its interior density profile and rotation period. Between its Love number k₂ and core size exists obviously a negative correlation. So oblateness and k₂ can well constrain its interior structure. Starting from the Lane-Emden equation, the planet models based on different polytropic indices are built. Then the flattening factors are obtained by solving the Wavre's integro-differential equation. The result shows that the smaller the polytropic index, the faster the rotation, and the larger the oblateness. We have selected 469 exoplanets, which have simultaneously the observed or estimated values of radius, mass, and orbit period from the NASA (National Aeronautics and Space Administration) Exoplanet Archive, and calculated their flattening factors under the two assumptions: tidal locking and fixed rotation period of 10.55 hours. The result shows that the flattening factors are too small to be detected under the tidal locking assumption, and that 28% of exoplanets have the flattening factors larger than 0.1 under the fixed rotation period of 10.55 hours. The Love numbers under the different polytropic models are solved by the Zharkov's approach, and the relation between k₂ and core size is discussed.     

Orbital resonances in the inner neptunian system: II. Resonant history of Proteus, Larissa, Galatea, and Despina

We investigate the orbital history of the small neptunian satellites discovered by Voyager 2. Over the age of the Solar System, tidal forces have caused the satellites to migrate radially, bringing them through mean-motion resonances with one another. In this paper, we extend our study of the largest satellites Proteus and Larissa [Zhang, K., Hamilton, D.P., 2007. Icarus 188, 386-399] by adding in mid-sized Galatea and Despina. We test the hypothesis that these moons all formed with zero inclinations, and that orbital resonances excited their tilts during tidal migration. We find that the current orbital inclinations of Proteus, Galatea, and Despina are consistent with resonant excitation if they have a common density . Larissa's inclination, however, is too large to have been caused by resonant kicks between these four satellites; we suggest that a prior resonant capture event involving either Naiad or Thalassa is responsible. Our solution requires at least three past resonances with Proteus, which helps constrain the tidal migration timescale and thus Neptune's tidal quality factor: 9000<QN<36,000. We also improve our determination of Qs for Proteus and Larissa, finding 36<QP<700 and 18<QL<200. Finally, we derive a more general resonant capture condition, and work out a resonant overlap criterion relevant to satellite orbital evolution around an oblate primary.     

基于无量纲的轨道参数开展卫星性质统计研究

目前已发现了285颗围绕太阳系八大行星公转的卫星, 它们的轨道和物理性质呈现了丰富多样性. 目前为止, 几乎所有的卫星研究工作都基于单个卫星系统或者卫星群, 似乎缺少统一的研究. 提出了一个新的与行星性质无关、只与恒星半径有关的轨道参数n, 定义为以太阳半径为单位的轨道半长轴的自然对数. 不同行星的卫星的n值都存在双极分布, 绝大部分卫星在$n\gtrsim2$区间, 其次在$n\lesssim-1$区间, 位于中间区域的行星则很少. 从卫星物理参数和轨道参数与n的关系中发现, 分属六大行星的卫星有明显的共同特征. 首先, 轨道偏心率和轨道倾角偏大的卫星的n值都在3.5左右, 它们都是巨行星的不规则卫星. 其次, n值在-1和1之间的卫星绝大部分体积大、质量大、反照率高、自转速度慢. 从文献中找到11颗系外卫星候选体, 获得了它们轨道n值和卫星质量, 发现后者也是在-1< n< 1区间最大,其他区间偏小.这些统一的 规律暗示,太阳系内不同行星的卫星形成机制以及太阳系外卫星的形成机制可能一样或类似.