Eccentricities of Planets in Binary Systems

The most puzzling property of the extrasolar planets discovered by recent radial velocity surveys is their high orbital eccentricities, which are very difficult to explain within our current theoretical paradigm for planet formation. Current data reveal that at least 25% of these planets, including some with particularly high eccentricities, are orbiting a component of a binary star system. The presence of a distant companion can cause significant secular perturbations in the orbit of a planet. At high relative inclinations, large-amplitude, periodic eccentricity perturbations can occur. These are known as “Kozai cycles” and their amplitude is purely dependent on the relative orbital inclination. Assuming that every planet host star also has a (possibly unseen, e.g., substellar) distant companion, with reasonable distributions of orbital parameters and masses, we determine the resulting eccentricity distribution of planets and compare it to observations? We find that perturbations from a binary companion always appear to produce an excess of planets with both very high (?0.6) and very low (e ? 0.1) eccentricities. The paucity of near-circular orbits in the observed sample implies that at least one additional mechanism must be increasing eccentricities. On the other hand, the overproduction of very high eccentricities observed in our models could be combined with plausible circularization mechanisms (e.g., friction from residual gas) to create more planets with intermediate eccentricities (e? 0.1–0.6).     

A keplerian method to estimate perturbations in the restricted three-body problem

We develop a new and fast method to estimate perturbations by a planet on cometary orbits. This method allows us to identify accurately the cases of large perturbations in a set of fictitious orbits. Hence, it can be used in constructing perturbation samples for Monte Carlo simulations in order to maximize the amount of information. Furthermore, the estimated perturbations are found to yield a good approximation to the real perturbation sample. This is shown by a comparison of the perturbations obtained by the new estimator with the results of numerical integration of regularized equations of motion for the same orbits in the same dynamical model: the three-dimensional elliptic restricted three-body problem (Sun-Jupiter-comet).     

Planet formation around M-dwarfs: the moving snow line and super-Earths

Planets result from a series of processes within a circumstellar disk. Evidence comes from the near planar orbits in the Solar System and other planetary systems, observations of newly formed disks around young stars, and debris disks around main-sequence stars. As planet-hunting techniques improve, we approach the ability to detect systems like the Solar System, and place ourselves in context with planetary systems in general. Along the way, new classes of planets with unexpected characteristics are discovered. One of the most recent classes contains super Earth-mass planets orbiting a few AU from low-mass stars. In this contribution, we outline a semi-analytic model for planet formation during the pre-main sequence contraction phase of a low-mass star. As the star contracts, the “snow line”, which separates regions of rocky planet formation from regions of icy planet formation, moves inward. This process enables rapid formation of icy protoplanets that collide and merge into super-Earths before the star reaches the main sequence. The masses and orbits of these super-Earths are consistent with super-Earths detected in recent microlensing experiments.     

Close encounters in young stellar clusters: implications for planetary systems in the solar neighbourhood

The stars that populate the solar neighbourhood were formed in stellar clusters. Through N -body simulations of these clusters, we measure the rate of close encounters between stars. By monitoring the interaction histories of each star, we investigate the singleton fraction in the solar neighbourhood. A singleton is a star which formed as a single star, has never experienced any close encounters with other stars or binaries, or undergone an exchange encounter with a binary. We find that, of the stars which formed as single stars, a significant fraction is not singletons once the clusters have dispersed. If some of these stars had planetary systems, with properties similar to those of the Solar System, the planets' orbits may have been perturbed by the effects of close encounters with other stars or the effects of a companion star within a binary. Such perturbations can lead to strong planet–planet interactions which eject several planets, leaving the remaining planets on eccentric orbits. Some of the single stars exchange into binaries. Most of these binaries are broken up via subsequent interactions within the cluster, but some remain intact beyond the lifetime of the cluster. The properties of these binaries are similar to those of the observed binary systems containing extrasolar planets. Thus, dynamical processes in young stellar clusters will alter significantly any population of Solar System-like planetary systems. In addition, beginning with a population of planetary systems exactly resembling the Solar System around single stars, dynamical encounters in young stellar clusters may produce at least some of the extrasolar planetary systems observed in the solar neighbourhood.     

Simple three-dimensional periodic orbits in a galactic-type potential

We study the families of simple periodic orbits in a three-dimensional system that represents the inner parts of a perturbed triaxial galaxy. The perturbations depend on two control parameters. We find the regions where each family is stable, simply unstable, doubly unstable, or complex unstable. the stable and simply unstable families produce other families by bifurcation. Several families reach a maximum (or minimum) perturbation and then are continued by other families. The bifurcations are direct or inverse. The transition from one type of bifurcation to the other is theoretically explained. Another important phenomenon is the splitting of one family into two, or the joining of two families into one. We do not have any complex instability in the limiting cases of two-dimensional motions (when one control parameter is zero).The two main families of periodic orbits are in most cases stable when the energy is smaller than the escape energy. Most high energy orbits are unstable. However, we found stable orbits even for energies about four times larger than the escape energy.     

Dynamical relaxation and massive extrasolar planets

Following the suggestion of Black that some massive extrasolar planets may be associated with the tail of the distribution of stellar companions, we investigate a scenario in which 5 N 100 planetary mass objects are assumed to form rapidly through a fragmentation process occuring in a disc or protostellar envelope on a scale of 100 au. These are assumed to have formed rapidly enough through gravitational instability or fragmentation that their orbits can undergo dynamical relaxation on a time-scale of ∼100 orbits.
Under a wide range of initial conditions and assumptions, the relaxation process ends with either (i) one potential 'hot Jupiter' plus up to two 'external' companions, i.e. planets orbiting near the outer edge of the initial distribution; (ii) one or two 'external' planets or even none at all; (iii) one planet on an orbit with a semi-major axis of 10 to 100 times smaller than the outer boundary radius of the inital distribution together with an 'external' companion. Most of the other objects are ejected and could contribute to a population of free-floating planets. Apart from the potential 'hot Jupiters', all the bound objects are on orbits with high eccentricity, and also with a range of inclination with respect to the stellar equatorial plane. We found that, apart from the close orbiters, the probability of ending up with a planet orbiting at a given distance from the central star increases with the distance. This is because of the tendency of the relaxation process to lead to collisions with the central star. The scenario we envision here does not impose any upper limit on the mass of the planets. We discuss the application of these results to some of the more massive extrasolar planets.     

On the migration of a system of protoplanets

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.     

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.     

Tidal tails of star clusters

Based on recent findings of a formation mechanism of substructure in tidal tails by Küpper et al., we investigate a more comprehensive set of N -body models of star clusters on orbits about a Milky Way like potential. We find that the predicted epicyclic overdensities arise in any tidal tail no matter which orbit the cluster follows as long as the cluster lives long enough for the overdensities to build up.
The distance of the overdensities along the tidal tail from the cluster centre depends for circular orbits only on the mass of the cluster and the strength of the tidal field, and therefore decreases monotonically with time, while for eccentric orbits the orbital motion influences the distance, causing a periodic compression and stretching of the tails and making the distance oscillate with time. We provide an approximation for estimating the distance of the overdensities in this case.
We describe an additional type of overdensity which arises in extended tidal tails of clusters on eccentric orbits, when the acceleration of the tidal field on the stellar stream is no longer homogeneous. Moreover, we conclude that a pericentre passage or a disc shock is not the direct origin of an overdensity within a tidal tail. Escape due to such tidal perturbations does not take place immediately after the perturbation but is rather delayed and spread over the orbit of the cluster. All observable overdensities are therefore of the mentioned two types. In particular, we note that substructured tidal tails do not imply the existence of dark matter substructures in the haloes of galaxies.     

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.     

Eccentric Extrasolar Planets: The Jumping Jupiter Model

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 10⁸ 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.     

Stability of inclined orbits of terrestrial planets in habitable zones

In long-term stability studies of terrestrial planets moving in the habitable zone (HZ) of a sun-like star, we distinguish four different configurations: (i) planets moving in binary star systems, (ii) the inner type (where the gas giant moves outside the HZ), (iii) the outer type (where the gas giant is closer to the star, than the HZ) and (iv) the Trojan type (where the gas giant moves in the HZ). Since earlier calculations indicated, that the stability of the motion in the HZ also depends on the inclination of the terrestrial planet orbits, we present a detailed numerical investigation to show correlations between the eccentricity, the mass and the distance of the giant planet for various inclinations of the terrestrial planets. The orbital stability of the HZ was examined for all four configurations stated above. While we could find hardly any stable orbits for the first three types for inclinations higher than 40^°, the Trojan planets can be stable up to an inclination of 60^°. Additionally, we could also find some stabilizing effects of the inclination for the first three types. As dynamical model we used the elliptic restricted three-body problem, which consists of two massive and one mass-less body. This allows an application to all detected and future extrasolar single planet systems.     

Analytical direct planetary perturbations on the orbital motion of satellites

An analytical expansion of the disturbing function arising from direct planetary perturbations on the motion of satellites is derived. As a Fourier series, it allows the investigation of the secular effects of these direct perturbations, as well as of every argument present in the perturbation. In particular, we construct an analytical model describing the evection resonance between the longitude of pericenter of the satellite orbit and the longitude of a planet, and study briefly its dynamic. The expansion developed in this paper is valid in the case of planar and circular planetary orbits, but not limited in eccentricity or inclination of the satellite orbit.     

Planetary microlensing signals from the orbital motion of the source star around the common barycentre

With several detections, the technique of gravitational microlensing has proven useful for studying planets that orbit stars at Galactic distances, and it can even be applied to detect planets in neighbouring galaxies. So far, planet detections by microlensing have been considered to result from a change in the bending of light and the resulting magnification caused by a planet around the foreground lens star. However, in complete analogy to the annual parallax effect caused by the revolution of the Earth around the Sun, the motion of the source star around the common barycentre with an orbiting planet can also lead to observable deviations in microlensing light curves that can provide evidence for the unseen companion. We discuss this effect in some detail and study the prospects of microlensing observations for revealing planets through this alternative detection channel. Given that small distances between lens and source star are favoured, and that the effect becomes nearly independent of the source distance, planets would remain detectable even if their host star is located outside the Milky Way with a sufficiently good photometry (exceeding present-day technology) being possible. From synthetic light curves arising from a Monte Carlo simulation, we find that the chances for such detections are not overwhelming and appear practically limited to the most massive planets (at least with current observational set-ups), but they are large enough for leaving the possibility that one or the other signal has already been observed. However, it may remain undetermined whether the planet actually orbits the source star or rather the lens star, which leaves us with an ambiguity not only with respect to its location, but also to its properties.     

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)     

The photometric method of detecting other planetary systems

The photometric method detects planets orbiting other stars by searching for the reduction in the light flux or the change in the color of the stellar flux that occurs when a planet transits a star. A transit by Jupiter or Saturn would reduce the stellar flux by approximately 1% while a transit by Uranus or Neptune would reduce the stellar flux by 0.1%. A highly characteristic color change with an amplitude approximately 0.1 of that for the flux reduction would also accompany the transit and could be used to verify that the source of the flux reduction was a planetary transit rather than some other phenomenon. Although the precision required to detect major planets is already available with state-of-the-art photometers, the detection of terrestrial-sized planets would require a precision substantially greater than the state-of-the-art and a spaceborne platform to avoid the effects of variations in sky transparency and scintillation. Because the probability is so small of observing a planetary transit during a single observation of a randomly chosen star, the search program must be designed to continuously monitor hundreds or thousands of stars. The most promising approach is to search for large planets with a photometric system that has a single-measurement precision of 0.1%. If it is assumed that large planets will have long-period orbits, and that each star has an average of one large planet, then approximately 10⁴ stars must be monitored continuously. To monitor such a large groups of stars simultaneously while maintaining the required photometric precision, a detector array coupled by a fiber-optic bundle to the focal plane of a moderate aperture (≈ 1 m), wide field of view (≈50°) telescope is required. Based on the stated assumptions, a detection rate of one planet per year of observation appears possible.     

The adapted solar wind system as cause for a momentum transfer to the sun and its consequences for the orbital motions of keplerian objects

In the following paper we argue that each wind-driving star in relative motion with respect to the ambient interstellar medium experiences a force exerted on its central wind-generating body. The exact magnitude of this force depends on the actual geometry of the counterflow configuration of stellar and interstellar winds for a particular kinematic situation which is especially sensitive to whether the interstellar flow is subsonic or supersonic. It will, however, be demonstrated here that this force is of an accelerating nature, i.e., it operates like a rocket-motor, as long as the peculiar motion of the wind-driving star with respect to the ambient interstellar medium remains subsonic.Here we use a specific analytical model to describe theoretically the specific counterflow configuration for the case of the solar system in a subsonic peculiar motion with respect to the local interstellar medium assuming irrotational and incompressible flows. We can work out a quantitative number for the accelerating force governing the Sun's motion at present. The net reaction force exerted on the solar body is then mediated by the asymmetric boundary conditions to which the distant solar wind field has to adapt.Next we study the indirect action of such a force on orbiting Keplerian objects like planets, planetesimals and comets. Since this force only influences the central solar body, but not the planets themselves, the problem is different from the treatment of a constant perturbation force perturbing the Keplerian orbits. We present a perturbation analysis treating the action of a corresponding position-dependent perturbation force resulting in secular changes of the orbital elements of Keplerian objects. It is found that changes are accumulating more rapidly in time the closer to the sun the orbiting bodies are. Main axis and perihelion distances are systematically increasing. Especially pronounced are changes in the perihelion position angle of the objects. For solar wind mass losses larger than the Sun's present value by a factor of 1000 (T-Tauri phase of the Sun,) the migration periods calculated for the planet Mercury are of the same order of magnitude as that for corresponding general relativistic migration.     

Debris disc stirring by secular perturbations from giant planets

Detectable debris discs are thought to require dynamical excitation ('stirring'), so that planetesimal collisions release large quantities of dust. We investigate the effects of the secular perturbations of a planet, which may lie at a significant distance from the planetesimal disc, to see if these perturbations can stir the disc, and if so over what time-scale. The secular perturbations cause orbits at different semimajor axes to precess at different rates, and after some time   t _cross  initially non-intersecting orbits begin to cross. We show that   t _cross∝ a ^9/2_disc/( m _pl e _pl a ³_pl)  , where   m _pl, e _pl  and   a _pl  are the mass, eccentricity and semimajor axis of the planet, and   a _disc  is the semimajor axis of the disc. This time-scale can be faster than that for the growth of planetesimals to Pluto's size within the outer disc. We also calculate the magnitude of the relative velocities induced among planetesimals and infer that a planet's perturbations can typically cause destructive collisions out to 100 s of au. Recently formed planets can thus have a significant impact on planet formation in the outer disc which may be curtailed by the formation of giant planets much closer to the star. The presence of an observed debris disc does not require the presence of Pluto-sized objects within it, since it can also have been stirred by a planet not in the disc. For the star ε Eridani, we find that the known radial velocity planet can excite the planetesimal belt at 60 au sufficiently to cause destructive collisions of bodies up to 100 km in size, on a time-scale of 40 Myr.     

Limits on the Orbits of Possible Eccentric and Inclined Moons of Extrasolar Planets Orbiting Single Stars

Limits are placed on the range of orbits and masses of possible moons orbiting extrasolar planets which orbit single central stars. The Roche limiting radius determines how close the moon can approach the planet before tidal disruption occurs; while the Hill stability of the star–planet–moon system determines stable orbits of the moon around the planet. Here the full three-body Hill stability is derived for a system with the binary composed of the planet and moon moving on an inclined, elliptical orbit relative the central star. The approximation derived here in Eq. (17) assumes the binary mass is very small compared with the mass of the star and has not previously been applied to this problem and gives the criterion against disruption and component exchange in a closed form. This criterion was applied to transiting extrasolar planetary systems discovered since the last estimation of the critical separations (Donnison in Mon Not R Astron Soc 406:1918, 2010a) for a variety of planet/moon ratios including binary planets, with the moon moving on a circular orbit. The effects of eccentricity and inclination of the binary on the stability of the orbit of a moon is discussed and applied to the transiting extrasolar planets, assuming the same planet/moon ratios but with the moon moving with a variety of eccentricities and inclinations. For the non-zero values of the eccentricity of the moon, the critical separation distance decreased as the eccentricity increased in value. Similarly the critical separation decreased as the inclination increased. In both cases the changes though very small were significant.     

Stability of S-type Orbits in Binaries

This numerical investigation is concerned with the stability of planets moving around one component of a double star system. Since the discovery of four extra solar planets moving in such orbits, there is a growing interest of stability studies thereto. We determined the stable regions in the elliptic restricted three body problem, for the whole range of mass-ratios from 0.1 to 0.9, by means of the Fast Lyapunov Indicators. The computations have been carried out for eccentricities of the binary and of the planet in the range 0–0.5. Therefore, we present for the first time the variation of the stable regions when the initial eccentricity of the planet is increased. We have found a correlation between the reduction of the stable zones if the eccentricity of the planet or of the binary is increased — of course the latter one is more effective.This revised version was published online in October 2005 with corrections to the Cover Date.