首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 46 毫秒
1.
Most extrasolar planets discovered to date are more massive than Jupiter, in surprisingly small orbits (semimajor axes less than 3 AU). Many of these have significant orbital eccentricities. Such orbits may be the product of dynamical interactions in multiplanet systems. We examine outcomes of such evolution in systems of three Jupiter-mass planets around a solar-mass star by integration of their orbits in three dimensions. Such systems are unstable for a broad range of initial conditions, with mutual perturbations leading to crossing orbits and close encounters. The time scale for instability to develop depends on the initial orbital spacing; some configurations become chaotic after delays exceeding 108 y. The most common outcome of gravitational scattering by close encounters is hyperbolic ejection of one planet. Of the two survivors, one is moved closer to the star and the other is left in a distant orbit; for systems with equal-mass planets, there is no correlation between initial and final orbital positions. Both survivors may have significant eccentricities, and the mutual inclination of their orbits can be large. The inner survivor's semimajor axis is usually about half that of the innermost starting orbit. Gravitational scattering alone cannot produce the observed excess of “hot Jupiters” in close circular orbits. However, those scattered planets with large eccentricities and small periastron distances may become circularized if tidal dissipation is effective. Most stars with a massive planet in an eccentric orbit should have at least one additional planet of comparable mass in a more distant orbit.  相似文献   

2.
On the migration of a system of protoplanets   总被引：1，自引：0，他引：1  
The evolution of a system consisting of a protoplanetary disc with two embedded Jupiter-sized planets is studied numerically. The disc is assumed to be flat and non-self-gravitating; this is modelled by the planar (two-dimensional) Navier–Stokes equations. The mutual gravitational interaction of the planets and the star, and the gravitational torques of the disc acting on the planets and the central star are included. The planets have an initial mass of one Jupiter mass M Jup each, and the radial distances from the star are one and two semimajor axes of Jupiter, respectively.
During the evolution a joint wide annular gap is created by the planets. Both planets increase their mass owing to accretion of gas from the disc: after about 2500 orbital periods of the inner planet it has reached a mass of 2.3  M Jup, while the outer planet has reached a mass of 3.2  M Jup. The net gravitational torques exerted by the disc on the planets result in an inward migration of the outer planet on time-scales comparable to the viscous evolution time of the disc. The semimajor axis of the inner planet remains constant as there is very little gas left in its vicinity to induce any migration. When the distance of close approach eventually becomes smaller than the mutual Hill radius, the eccentricities increase strongly and the system may become unstable.
If disc depletion occurs rapidly enough before the planets come too close to each other, a stable system similar to our own Solar system may remain. Otherwise the orbits may become unstable and produce systems like υ And.  相似文献   

3.
赵佳  赵刚 《天文学进展》2012,30(1):48-63
自1995年第一颗类太阳恒星周围的系外行星发现以来,随着已发现的系外行星数目的增多,对系外行星性质的统计分析变得重要和有意义。截至2011年6月9日,共发现系外行星555颗。以这些系外行星的轨道参数为依据,对系外行星的性质进行统计分析,得出了一些有意义的结论。同时简要介绍现有的行星形成与演化模型并依据得出的行星统计性质对其进行检验,这对于系外行星的进一步探测具有一定的指导作用。  相似文献   

4.
Andrew W. Smith 《Icarus》2009,201(1):381-58
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.  相似文献   

5.
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−4 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 (∼10−4-10−2 times MMSN).  相似文献   

6.
In this article we investigate the masses, orbital periods, semimajor axis, eccentricities and radii of the existing exoplanets by comparing the first and second digit probabilities with Benford laws’s predictions. It is found that the masses, orbital periods and semimajor axis conform to Benfordós law quite well, but radii fail. It is also investigated the first digits occurrence corresponding to a given order of magnitude. We introduce a top function which can estimate all the probabilities of the first digit order-to-order.  相似文献   

7.
We investigate the dynamics of putative Earth-mass planets in the habitable zone (HZ) of the extrasolar planetary system OGLE-2006-BLG-109L, a close analogue of the Solar system. Our work is inspired by the work of Malhotra & Minton. Using the linear Laplace–Lagrange theory, they identified a strong secular resonance that may excite large eccentricity of orbits in the HZ. However, due to uncertain or unconstrained orbital parameters, the subsystem of Jupiters may be found in a dynamically active region of the phase space spanned by low-order mean-motion resonances. To generalize this secular model, we construct a semi-analytical averaging method in terms of the restricted problem. The secular orbits of large planets are approximated by numerically averaged osculating elements. They are used to calculate the mean orbits of terrestrial planets by means of a high-order analytic secular theory developed in our previous works. We found regions in the parameter space of the problem in which stable, quasi-circular orbits in the HZ are permitted. The excitation of eccentricity in the HZ strongly depends on the apsidal angle of jovian orbits. For some combinations of that angle, eccentricities and semimajor axes consistent with the observations, a terrestrial planet may survive in low eccentric orbits. We also study the effect of post-Newtonian gravity correction on the innermost secular resonance.  相似文献   

8.
We numerically investigate the stability of systems of 1 \({{\rm M}_{\oplus}}\) planets orbiting a solar-mass star. The systems studied have either 2 or 42 planets per occupied semimajor axis, for a total of 6, 10, 126, or 210 planets, and the planets were started on coplanar, circular orbits with the semimajor axes of the innermost planets at 1 AU. For systems with two planets per occupied orbit, the longitudinal initial locations of planets on a given orbit were separated by either 60° (Trojan planets) or 180°. With 42 planets per semimajor axis, initial longitudes were uniformly spaced. The ratio of the semimajor axes of consecutive coorbital groups in each system was approximately uniform. The instability time for a system was taken to be the first time at which the orbits of two planets with different initial orbital distances crossed. Simulations spanned virtual times of up to 1 × 108, 5 × 105, and 2 × 105 years for the 6- and 10-planet, 126-planet, and 210-planet systems, respectively. Our results show that, for a given class of system (e.g., five pairs of Trojan planets orbiting in the same direction), the relationship between orbit crossing times and planetary spacing is well fit by the functional form log(t c /t 0) = b β + c, where t c is the crossing time, t 0 = 1 year, β is the separation in initial orbital semimajor axis (in terms of the mutual Hill radii of the planets), and b and c are fitting constants. The same functional form was observed in the previous studies of single planets on nested orbits (Smith and Lissauer 2009). Pairs of Trojan planets are more stable than pairs initially separated by 180°. Systems with retrograde planets (i.e., some planets orbiting in the opposite sense from others) can be packed substantially more closely than can systems with all planets orbiting in the same sense. To have the same characteristic lifetime, systems with 2 or 42 planets per orbit typically need to have about 1.5 or 2 times the orbital separation as orbits occupied by single planets, respectively.  相似文献   

9.
We study the dynamics of planetary systems with two planets moving in the same plane, when frictional forces act on the two planets, in addition to the gravitational forces. The model of the general three-body problem is used. Different laws of friction are considered. The topology of the phase space is essential in understanding the evolution of the system. The topology is determined by the families of stable and unstable periodic orbits, both symmetric and non symmetric. It is along the stable families, or close to them, that the planets migrate when dissipative forces act. At the critical points where the stability along the family changes, there is a bifurcation of a new family of stable periodic orbits and the migration process changes route and follows the new stable family up to large eccentricities or to a chaotic region. We consider both resonant and non resonant planetary systems. The 2/1, 3/1 and 3/2 resonances are studied. The migration to larger or smaller eccentricities depends on the particular law of friction. Also, in some cases the semimajor axes increase and in other cases they are stabilized. For particular laws of friction and for special values of the parameters of the frictional forces, it is possible to have partially stationary solutions, where the eccentricities and the semimajor axes are fixed.  相似文献   

10.
In this work, we study the stability of hypothetical satellites of extrasolar planets. Through numerical simulations of the restricted elliptic three-body problem we found the borders of the stable regions around the secondary body. From the empirical results, we derived analytical expressions of the critical semimajor axis beyond which the satellites would not remain stable. The expressions are given as a function of the eccentricities of the planet, e P, and of the satellite, e sat. In the case of prograde satellites, the critical semimajor axis, in the units of Hill's radius, is given by a E≈ 0.4895   (1.0000 − 1.0305 e P− 0.2738 e sat). In the case of retrograde satellites, it is given by a E≈ 0.9309  (1.0000 − 1.0764 e P− 0.9812 e sat). We also computed the satellite stability region ( a E) for a set of extrasolar planets. The results indicate that extrasolar planets in the habitable zone could harbour the Earth-like satellites.  相似文献   

11.
Observations and results of orbit determination of the first known Mars Trojan asteroid (5261) Eureka are presented. We have numerically calculated the evolution of the orbital elements, and have analyzed the behavior of the motion during the next 2 Myr. Strong perturbations by planets other than Mars seem to stabilize the eccentricity of the asteroid by stirring the high order resonances present in the elliptic restricted problem. As a result, the orbit appears stable at least on megayear timescales. The difference of the mean longitudes of Mars and Eureka and the semimajor axis of the asteroid form a pair of variables that essentially behave in an adiabatic manner, while the evolution of the other orbital elements is largely determined by the perturbations due to other planets.  相似文献   

12.
We analyse the secular interactions of two coplanar planets which are not in mean motion resonances. The analysis is based on a high order (order 12) expansion of the perturbative potential in powers of the eccentricities. The model depends on only two parameters (the ratio of semi-major axis and the mass ratio of the planets) and can be reduced to a one degree of freedom system, allowing for an exhaustive parametric analysis. Following Pauwels [Pauwels T.: 1983, Celet. Mech. & Dyn. Astro. 30, 229–247] we map the phase space on a sphere, avoiding in this way the artificial singularities introduced by other mappings. We show that the 12 order expansion is able to describe correctly most of the exosolar planetary systems discovered so far, even if the eccentricities of these planets are considerably larger than the eccentricities of our own solar system. The expansion is even able to reproduce, at moderate eccentricities, the secular resonances discovered numerically by Michtchenko and Malhotra [Michtchenko, T. A. and Malhotra, R.: 2004, Icarus 168, 237–248] at moderate to large eccentricities. FNRS Research Fellow.  相似文献   

13.
Numerous studies in the past few years have analyzed possible effects of planetary migration on the small bodies of the Solar System (mainly asteroids and KBOs), with the double aim of explaining certain dynamical structures in these systems, as well as placing limits on the magnitude of the radial migration of the planets. Here we undertake a similar aim, only this time concentrating on the dynamical stability of planetary satellites in a migration scenario. However, different from previous works, the strongest perturbations on satellite systems are not due to the secular variation of the semimajor axes of the planets, but from the planetesimals themselves. These perturbations result from close approaches between the planetesimals and satellites.We present results of several numerical simulations of the dynamical evolution of real and fictitious satellite systems around the outer planets, under the effects of multiple passages of a population of planetesimals representing the large-body component of a residual rocky disk. Assuming that this component dominated the total mass of the disk, our results show that the present systems of satellites of Uranus and Neptune do not seem to be compatible with a planetary migration larger than even one quarter that suggested by previous studies, unless these bodies were originated during the late stage of evaporation of the planetesimal disk. For larger variations of the semimajor axes of the planets, most of the satellites would either be ejected from the system or suffer mutual collisions due to excitation in their eccentricities. For the systems of Jupiter and Saturn, these perturbations are not so severe, and even large migrations do not introduce large instabilities.Nevertheless, even a small number of 1000-km planetesimals in the region may introduce significant excitation in the eccentricities and inclinations of satellites. Adequate values of this component may help explain the present dynamical distribution of distant satellites, including the highly peculiar orbit of Nereid.  相似文献   

14.
The significant orbital eccentricities of most giant extrasolar planets may have their origin in the gravitational dynamics of initially unstable multiple planet systems. In this work, we explore the dynamics of two close planets on inclined orbits through both analytical techniques and extensive numerical scattering experiments. We derive a criterion for two equal mass planets on circular inclined orbits to achieve Hill stability, and conclude that significant radial migration and eccentricity pumping of both planets occurs predominantly by 2:1 and 5:3 mean motion resonant interactions. Using Laplace-Lagrange secular theory, we obtain analytical secular solutions for the orbital inclinations and longitudes of ascending nodes, and use those solutions to distinguish between the secular and resonant dynamics which arise in numerical simulations. We also illustrate how encounter maps, typically used to trace the motion of massless particles, may be modified to reproduce the gross instability seen by the numerical integrations. Such a correlation suggests promising future use of such maps to model the dynamics of more coplanar massive planet systems.  相似文献   

15.
We consider particles with low free or proper eccentricity that are orbiting near planets on eccentric orbits. Through collisionless particle integration, we numerically find the location of the boundary of the chaotic zone in the planet's corotation region. We find that the distance in semimajor axis between the planet and boundary depends on the planet mass to the 2/7 power and is independent of the planet eccentricity, at least for planet eccentricities below 0.3. Our integrations reveal a similarity between the dynamics of particles at zero eccentricity near a planet in a circular orbit and with zero free eccentricity particles near an eccentric planet. The 2/7th law has been previously explained by estimating the semimajor at which the first-order mean motion resonances are large enough to overlap. Orbital dynamics near an eccentric planet could differ due to first-order corotation resonances that have strength proportional to the planet's eccentricity. However, we find that the corotation resonance width at low free eccentricity is small; also the first-order resonance width at zero free eccentricity is the same as that for a zero-eccentricity particle near a planet in a circular orbit. This accounts for insensitivity of the chaotic zone width to planet eccentricity. Particles at zero free eccentricity near an eccentric planet have similar dynamics to those at zero eccentricity near a planet in a circular orbit.  相似文献   

16.
The Gliese 876 planetary system consists of two Jupiter-like planets having a nearly commensurate 2:1 orbital periods ratio. Because the semimajor axes of the planets are very small (of the order 0.1 au and 0.2 au, respectively), and the eccentricity of the inner companion is ≃0.3, the mutual perturbations are extremely large. However, many authors claim the long-term orbital stability of the system, at least over 500 Myr for initial conditions found by Rivera & Lissauer. Results of investigations of a migration of initially separated planets into the close 2:1 mean motion resonance lock from Lee & Peale also support the conclusion that the system should be stable for the lifetime of the parent star. Initial conditions of the system, found from non-linear N -body fits by Laughlin & Chambers and Rivera & Lissauer, to the radial velocity curve, formally allow for a variety of orbital configurations of the GJ 876 system, e.g. coplanar, with planetary inclinations in the range [≃30°, 90°], and with relative inclinations of orbital planes as high as 80°. Our work is devoted to the stability investigation of the systems originating from the fitted initial conditions. We study neighbourhoods of these initial states in the orbital parameter space. We found estimations of the 2:1 mean motion resonance width and dynamical limitations on the planetary masses. We also obtain a global representation of the domains of the orbital parameters space in which initial conditions leading to stable evolutions can be found. Our results can be useful in localization of the best, stable fits to the observational data. In our investigations we use the MEGNO technique (the Mean Exponential Growth factor of Nearby Orbits) invented by Cincotta & Simó. It allows us to distinguish efficiently and precisely between chaotic and regular behaviour of a planetary system.  相似文献   

17.
Rodney S Gomes 《Icarus》2003,161(2):404-418
I simulate the orbital evolution of the four major planets and a massive primordial planetesimal disk composed of 104 objects, which perturb the planets but not themselves. As Neptune migrates by energy and angular momentum exchange with the planetesimals, a large number of primordial Neptune-scattered objects are formed. These objects may experience secular, Kozai, and mean motion resonances that induce temporary decrease of their eccentricities. Because planets are migrating, some planetesimals can escape those resonances while in a low-eccentricity incursion, thus avoiding the return path to Neptune close encounter dynamics. In the end, this mechanism produces stable orbits with high inclination and moderate eccentricities. The population so formed together with the objects coming from the classical resonance sweeping process, originates a bimodal distribution for the Kuiper Belt orbits. The inclinations obtained by the simulations can attain values above 30° and their distribution resembles a debiased distribution for the high-inclination population coming from the real classical Kuiper Belt.  相似文献   

18.
Planetary, stellar and galactic physics often rely on the general restricted gravitational $N$ -body problem to model the motion of a small-mass object under the influence of much more massive objects. Here, I formulate the general restricted problem entirely and specifically in terms of the commonly used orbital elements of semimajor axis, eccentricity, inclination, longitude of ascending node, argument of pericentre, and true anomaly, without any assumptions about their magnitudes. I derive the equations of motion in the general, unaveraged case, as well as specific cases, with respect to both a bodycentric and barycentric origin. I then reduce the equations to three-body systems, and present compact singly- and doubly-averaged expressions which can be readily applied to systems of interest. This method recovers classic Lidov–Kozai and Laplace–Lagrange theory in the test particle limit to any order, but with fewer assumptions, and reveals a complete analytic solution for the averaged planetary pericentre precession in coplanar circular circumbinary systems to at least the first three nonzero orders in semimajor axis ratio. Finally, I show how the unaveraged equations may be used to express resonant angle evolution in an explicit manner that is not subject to expansions of eccentricity and inclination about small nor any other values.  相似文献   

19.
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 107 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 108 years, planets grow to Earth's size, but the final eccentricities are too high as in gas-free cases. If it is shorter than 106 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 106 and 108 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×106 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−4-10−3 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.  相似文献   

20.
We study orbits of planetary systems with two planets, for planar motion, at the 1/1 resonance. This means that the semimajor axes of the two planets are almost equal, but the eccentricities and the position of each planet on its orbit, at a certain epoch, take different values. We consider the general case of different planetary masses and, as a special case, we consider equal planetary masses. We start with the exact resonance, which we define as the 1/1 resonant periodic motion, in a rotating frame, and study the topology of the phase space and the long term evolution of the system in the vicinity of the exact resonance, by rotating the orbit of the outer planet, which implies that the resonance and the eccentricities are not affected, but the symmetry is destroyed. There exist, for each mass ratio of the planets, two families of symmetric periodic orbits, which differ in phase only. One is stable and the other is unstable. In the stable family the planetary orbits are in antialignment and in the unstable family the planetary orbits are in alignment. Along the stable resonant family there is a smooth transition from planetary orbits of the two planets, revolving around the Sun in eccentric orbits, to a close binary of the two planets, whose center of mass revolves around the Sun. Along the unstable family we start with a collinear Euler–Moulton central configuration solution and end to a planetary system where one planet has a circular orbit and the other a Keplerian rectilinear orbit, with unit eccentricity. It is conjectured that due to a migration process it could be possible to start with a 1/1 resonant periodic orbit of the planetary type and end up to a satellite-type orbit, or vice versa, moving along the stable family of periodic orbits.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号