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.”     

A dynamical analysis of the 14 Herculis planetary system

Precision radial velocity measurements of the Sun-like dwarf 14 Herculis published by Naef et al., Butler et al. and Wittenmyer, Endl & Cochran reveal a Jovian planet in a 1760-d orbit and a trend indicating the second distant object. On the grounds of dynamical considerations, we test a hypothesis that the trend can be explained by the presence of an additional giant planet. We derive dynamical limits to the orbital parameters of the putative outer Jovian companion in an orbit within ∼13 au. In this case, the mutual interactions between the Jovian planets are important for the long-term stability of the system. The best self-consistent and stable Newtonian fit to an edge-on configuration of Jovian planets has the outer planet in 9-au orbit with a moderate eccentricity of ∼0.2 and confined to a zone spanned by the low-order mean motion resonances 5:1 and 6:1. This solution lies in a shallow minimum of (χ²_ν)^1/2 and persists over a wide range of the system inclination. Other stable configurations within 1σ confidence interval of the best fit are possible for the semimajor axis of the outer planet in the range of (6,13) au and the eccentricity in the range of (0, 0.3). The orbital inclination cannot yet be determined but when it decreases, both planetary masses approach ∼10 m _J and for i ∼ 30° the hierarchy of the masses is reversed.     

The potential for Earth-mass planet formation around brown dwarfs

Recent observations point to the presence of structured dust grains in the discs surrounding young brown dwarfs, thus implying that the first stages of planet formation take place also in the substellar regime. Here, we investigate the potential for planet formation around brown dwarfs and very low-mass stars according to the sequential core accretion model of planet formation. We find that, for a brown dwarf mass 0.05 M_⊙, our models predict a maximum planetary mass of  ∼5   M_⊕  , orbiting with semimajor axis ∼ 1 au. However, we note that the predictions for the mass–semimajor axis distribution are strongly dependent upon the models chosen for the disc surface density profiles and the assumed distribution of disc masses. In particular, if brown dwarf disc masses are of the order of a few Jupiter masses, Earth-mass planets might be relatively frequent, while if typical disc masses are only a fraction of Jupiter mass, we predict that planet formation would be extremely rare in the substellar regime. As the observational constraints on disc profiles, mass dependencies and their distributions are poor in the brown dwarf regime, we advise caution in validating theoretical models only on stars similar to the Sun and emphasize the need for observational data on planetary systems around a wide range of stellar masses. We also find that, unlike the situation around solar-like stars, Type II migration is totally absent from the planet formation process around brown dwarfs, suggesting that any future observations of planets around brown dwarfs would provide a direct measure of the role of other types of migration.     

New methods for large dynamic range problems in planetary formation

Modern N -body techniques for planetary dynamics are generally based on symplectic algorithms specially adapted to the Kepler problem. These methods have proven very useful in studying planet formation, but typically require the time-step for all objects to be set to a small fraction of the orbital period of the innermost body. This computational expense can be prohibitive for even moderate particle number for many physically interesting scenarios, such as recent models of the formation of hot exoplanets, in which the semimajor axis of possible progenitors can vary by orders of magnitude. We present new methods which retain most of the benefits of the standard symplectic integrators but allow for radial zones with distinct time-steps. These approaches should make simulations of planetary accretion with large dynamic range tractable. As proof-of-concept, we present preliminary science results from an implementation of the algorithm as applied to an oligarchic migration scenario for forming hot Neptunes.     

A planetary companion orbiting the intermediate-mass G Giant HD 173416

The discovery of a planetary companion to the intermediate-mass late-type giant star HD173416 from precise Doppler surveys of G and K giants at Xingiong station and Okayama Astrophysical Observatory (OAO) is presented in this letter. The planet has a minimum mass of 2.7 MJ an eccentricity of 0.21, a semimajor axis of 1.16 AU and an orbital period of 324 days.     

The effect on the Edgeworth–Kuiper Belt of a large distant tenth planet

We investigate the orbital evolution of both real and hypothetical Edgeworth–Kuiper Objects in order to determine whether any conclusions can be drawn regarding the existence, or otherwise, of the tenth planet postulated by Murray. We find no qualitative difference in the orbital evolution, and so conclude that the hypothetical planet has been placed on an orbit at such a large heliocentric distance that no evidence for the existence, or non-existence, can be found from a study of the known Edgeworth–Kuiper Objects.     

From the Oort cloud to observable short-period comets – I. The initial stage of cometary capture

We investigate the first stage of the dynamical evolution of Oort cloud comets entering the planetary region for the first time. To this purpose, we integrate numerically the motions of a large number of fictitious comets pertaining to two samples, both with perihelion distances up to 5.7 au and random inclinations; the first sample is composed of comets whose orbits have at least one node close to 5.2 au, while the second is not subject to this constraint. We examine the orbits when the comets come to aphelion after their first perihelion passage within the planetary region, and find that there is a clear statistical dependence of the energy perturbations on the Tisserand parameter. There appear to be two main processes, of comparable importance, governing the shortening of semimajor axes to values of less than 1000 au, i.e. planetary close encounters, especially with Jupiter, and indirect perturbations due to the shifting of the motion from barycentric to heliocentric and back; the former process mostly affects comets crossing the ecliptic at about 5.2 au, or on low-inclination orbits, while the latter mostly affects comets of small perihelion distance. This last result may help to understand the relative paucity of Halley-type comets with perihelion distances larger than about 1.5 au.     

Why are there so few hot Jupiters?

We use numerical simulations to model the migration of massive planets at small radii and compare the results with the known properties of 'hot Jupiters' (extrasolar planets with semimajor axes   a < 0.1  au). For planet masses   M _pl sin  i > 0.5 M _J  , the evidence for any 'pile-up' at small radii is weak (statistically insignificant), and although the mass function of hot Jupiters is deficient in high-mass planets as compared to a reference sample located further out, the small sample size precludes definitive conclusions. We suggest that these properties are consistent with disc migration followed by entry into a magnetospheric cavity close to the star. Entry into the cavity results in a slowing of migration, accompanied by a growth in orbital eccentricity. For planet masses in excess of 1 Jupiter mass we find eccentricity growth time-scales of a few ×10⁵ yr, suggesting that these planets may often be rapidly destroyed. Eccentricity growth appears to be faster for more massive planets which may explain changes in the planetary mass function at small radii and may also predict a pile-up of lower mass planets, the sample of which is still incomplete.     

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)     

Morphological and kinematic signatures of a binary central star in the planetary nebula Hu 2-1

We present H α , [N  ii ] and [O  iii ] ground-based and HST archive images, VLA–A 3.6-cm continuum and H92 α emission-line data and high-resolution long-slit [N  ii ] spectra of the planetary nebula Hu 2-1. A large number of structural components are identified in the nebula: an outer bipolar and an inner shell, two pairs of collimated bipolar structures at different directions, monopolar bow-shock-like structures, and an extended equatorial structure within a halo. The formation of Hu 2-1 appears to be dominated by anisotropic mass ejection during the late-AGB stage of the progenitor and by variable, 'precessing' collimated bipolar outflows during the protoplanetary nebula and/or early planetary nebula phases. Different observational results strongly support the existence of a binary central star in Hu 2-1, among them (1) the observed point-symmetry of the bipolar lobes and inner shell, and the departures from axial symmetry of the bipolar lobes, (2) the off-centre position of the central star, (3) the detection of mass ejection towards the equatorial plane, and (4) the presence of 'precessing' collimated outflows. In addition, (5) an analysis of the kinematics shows that the systemic velocity of the bipolar outflows does not coincide with the systemic velocity of the bipolar shell. We propose that this velocity difference is a direct evidence of orbital motion of the ejection source in a binary system. From a deduced orbital velocity of ∼10 km s⁻¹, a semimajor axis of ∼ 9–27 au and period of ∼ 25–80 yr are obtained, assuming a reasonable range of masses. These parameters are used to analyse the formation of Hu 2-1 within current scenarios of planetary nebulae with binary central stars.     

Detecting planets in planetary nebulae

We examine the possibility of detecting signatures of surviving Uranus/Neptune-like planets inside planetary nebulae. Planets that are not too close to the stars (orbital separation larger than ∼5 au) are likely to survive the entire evolution of the star. As the star turns into a planetary nebula, it has a fast wind and strong ionizing radiation. The interaction of the radiation and wind with a planet may lead to the formation of a compact condensation or tail inside the planetary nebula, which emits strongly in H α , but not in [O  iii ]. The position of the condensation (or tail) will change over a time-scale of ∼10 yr. Such condensations might be detected with currently existing telescopes.     

Motion of dust in mean motion resonances with planets

Effect of stellar electromagnetic radiation on the motion of spherical dust particle in mean motion orbital resonances with a planet is investigated. Planar circular restricted three-body problem with the Poynting–Robertson (P–R) effect yields monotonic secular evolution of eccentricity when the particle is trapped in the resonance. Planar elliptic restricted three-body problem with the P–R effect enables nonmonotonous secular evolution of eccentricity and the evolution of eccentricity is qualitatively consistent with the published results for the complicated case of interaction of electromagnetic radiation with nonspherical dust grain. Thus, it is sufficient to allow either nonzero eccentricity of the planet or nonsphericity of the grain and the orbital evolutions in the resonances are qualitatively equal for the two cases. This holds both for exterior and interior mean motion orbital resonances. Evolutions of argument of perihelion in the planar circular and elliptical restricted three-body problems are shown. Numerical integrations show that an analytic expression for the secular time derivative of the particle’s argument of perihelion does not exist, if only dependence on semimajor axis, eccentricity and argument of perihelion is admitted. Connection between the shift of perihelion and oscillations in secular eccentricity is presented for the planar elliptic restricted three-body problem with the P–R effect. Period of the oscillations corresponds to the period of one revolution of perihelion. Change of optical properties of the spherical grain with the heliocentric distance is also considered. The change of the optical properties: (i) does not have any significant influence on the secular evolution of eccentricity, (ii) causes that the shift of perihelion is mainly in the same direction/orientation as the particle motion around the Sun. The statements hold both for circular and noncircular planetary orbits.     

Kuiper带天体的轨道分布特性

1992年9月,夏威夷大学的D.Jewitt和加利福尼亚大学的J.Lun发现了海王星外绕太阳运行的第一个小天体1992QB1[1],开创了人类对于海王星外天体的实际观测的研究.近10年的接连不断发现,已经证实了海王星轨道外面存在着一个由大量的环绕太阳运动的小天体组成的环带[2].由于G.P.Kuiper曾在1951年的文章中提出过在冥王星的外边可能存在小天体的问题,因此人们一般把这个环带称为Kuiper带,你这些天体为“KuiperBelt　Objects”(KBOs),或从逻辑上称它们为“Trans-NeptunianObjects”(TNOs)[3]     

Dynamical behaviour of planetesimals temporarily captured by a planet from heliocentric orbits: basic formulation and the case of low random velocity

Planetesimals encountering with a planet cannot be captured permanently unless energy dissipation is taken into account, but some of them can be temporarily captured in the vicinity of the planet for an extended period of time. Such a process would be important for the origin and dynamical evolution of irregular satellites, short-period comets, and Kuiper-belt binaries. In this paper, we describe the basic formulation for the study of temporary capture of planetesimals from heliocentric orbits using three-body orbital integration, such as the definition of the duration and rate of temporary capture, and present results in the case of low random velocity of planetesimals. In the case of planetesimals initially on circular orbits, we find that planetesimals undergo a close encounter with the planet before they become temporarily captured. When planetesimals are scattered by the planet into the vicinity of one of periodic orbits around the planet, the duration of temporary capture tends to be extended. Typically, these capture orbits are in the retrograde direction around the planet. We evaluate the rate of temporary capture of planetesimals, and find that the ratio of this rate to their collision rate on to the planet increases with increasing semimajor axis of the planet. Similar results are obtained for planetesimals with non-zero but small random velocities, as long as Kepler shear dominates the relative velocity between the planet and planetesimals. For larger initial random velocities of planetesimals, temporary capture in both prograde and retrograde directions with much longer duration becomes possible.     

The Effect of Sweeping Secular Resonances on the Classical Kuiper Belt

Abstract The Kuiper Belt is a disk of small icy objects orbiting the Sun beyond Neptune. The region between 40-48AU in this disk is supposed to consist of dynamical “cold” objects on low-inclination orbits and is called the “Classical Kuiper Belt”. Recent observations reveal that there is a “hot” population with inclinations being as large as 30? residing in this region. Secular resonance sweeping, which took place in the late stage of formation of the planetary system when the residual nebula gas was dispersing, is a possible mechanism that can excite the orbits in this region. In this paper, we investigate in detail the excitation of orbital inclination by this mechanism. It is shown that the excitation depends sensitively on the angle δ between the midplane of the nebula gas and the invariable plane of the solar system. The excitation is very small when δ = 0?, but if the gas midplane coincides with the ecliptic, i.e. if δ ≈ 1.6?, then objects in the region of classical Kuiper belt can be excited to orbital inclinations as high as 30?, provided the nebula gas has the proper initial density and disperses at a proper rate. We also considered the orbital excitation by secular resonance sweeping with Jupiter on an inclined orbit and with migrating Jovian planets, and found the excitation is only slightly affected.     

Analysis of the Hipparcos Measurements of HD 10697: A Mass Determination of a Brown Dwarf Secondary

HD 10697 is a nearby main-sequence star around which a planet candidate has recently been discovered by means of radial velocity measurements (Vogt et al.). The stellar orbit has a period of about 3 yr, the secondary minimum mass is 6.35 Jupiter masses (MJ), and the minimum semimajor axis is 0.36 mas. Using the Hipparcos data of HD 10697 together with the spectroscopic elements of Vogt et al., we found a semimajor axis of 2.1+/-0.7 mas, implying a mass of 38+/-13 MJ for the unseen companion. We therefore suggest that the secondary of HD 10697 is probably a brown dwarf, orbiting around its parent star at a distance of 2 AU.     

Some properties of a two-body system under the influence of the Galactic tidal field

In this paper the effect of the Galactic tidal field on a Sun–comet pair will be considered when the comet is situated in the Oort cloud and planetary perturbations can be neglected. First, two averaged models were created, one of which can be solved analytically in terms of Jacobi elliptic functions. In the latter system, switching between libration and circulation of the argument of perihelion is prohibited. The non-averaged equations of motion are integrated numerically in order to determine the regions of the ( e ,  i ) phase space in which chaotic orbits can be found, and an effort is made to explain why the chaotic orbits manifest in these regions only. It is evident that for moderate values of semimajor axis, a ∼50 000 au , chaotic orbits are found for ( e <0.15 , 40°≤ i ≤140°) as determined by integrating the evolution of the comet over a period of 10⁴ orbits. These regions of chaos increase in size with increasing semimajor axis. The typical e-folding times for these orbits range from around 600 Myr to 1 Gyr and thus are of little practical interest, as the time-scales for chaos arising from passing stars are much shorter. As a result of Galactic rotation, the chaotic regions in ( e ,  i ) phase space are not symmetric for prograde and retrograde orbits.     

Core instability models of giant planet accretion – II. Forming planetary systems

We develop a simple model for computing planetary formation based on the core instability model for the gas accretion and the oligarchic growth regime for the accretion of the solid core. In this model several planets can form simultaneously in the disc, a fact that has important implications especially for the changes in the dynamic of the planetesimals and the growth of the cores since we consider the collision between them as a source of potential growth. The type I and type II migration of the embryos and the migration of the planetesimals due to the interaction with the disc of gas are also taken into account. With this model we consider different initial conditions to generate a variety of planetary systems and analyse them statistically. We explore the effects of using different type I migration rates on the final number of planets formed per planetary system such as on the distribution of masses and semimajor axis of extrasolar planets, where we also analyse the implications of considering different gas accretion rates. A particularly interesting result is the generation of a larger population of habitable planets when the gas accretion rate and type I migration are slower.     

A Bayesian Kepler periodogram detects a second planet in HD 208487

An automatic Bayesian Kepler periodogram has been developed for identifying and characterizing multiple planetary orbits in precision radial velocity data. The periodogram is powered by a parallel tempering Markov chain Monte Carlo (MCMC) algorithm which is capable of efficiently exploring a multiplanet model parameter space. The periodogram employs an alternative method for converting the time of an observation to true anomaly that enables it to handle much larger data sets without a significant increase in computation time. Improvements in the periodogram and further tests using data from HD 208487 have resulted in the detection of a second planet with a period of 909⁸²₋₉₂ d, an eccentricity of 0.37^0.26_−0.20, a semimajor axis of 1.87^0.13_−0.14 au and an M sin  i = 0.45^+0.11_−0.13 M _J. The revised parameters of the first planet are period = 129.8 ± 0.4 d, eccentricity = 0.20 ± 0.09, semimajor axis = 0.51 ± 0.02 au and M sin  i = 0.41 ± 0.05  M _J. Particular attention is paid to several methods for calculating the model marginal likelihood which is used to compare the probabilities of models with different numbers of planets.     

太阳多方模型对行星轨道要素变化的长期摄动影响

本文利用天体力学中的摄动理论和天体物理学中的气体星多方模型理论研究了太阳多方模型对行星轨道要素变化的长期摄动影响。文中给出了太阳日多方指数ｎ＝３的模型由于自转、扁度和内部密度分布等因素对行星轨道要素变化的长期摄动效应的理论结果。研究结果表明：行星轨道要素除长轴、偏心率和轨道倾角不受长期摄动外升交点经度、近日点经度以及平近点角均受长期摄动的影响。最后利用理论结果对行星轨道要素的长期投动效应做了数值计算，数值结果在表１中给出。