Dynamics of rotation of super-Earths

We numerically investigate the dynamics of rotation of several close-in terrestrial exoplanet candidates. In our model, the rotation of the planet is disturbed by the torque of the central star due to the asymmetric equilibrium figure of the planet. We model the shape of the planet by a Jeans spheroid. We use surfaces of section and spectral analysis to explore numerically the rotation phase space of the systems adopting different sets of parameters and initial conditions close to the main spin–orbit resonant states. One of the parameters, the orbital eccentricity, is critically discussed here within the domain of validity of orbital circularization timescales given by tidal models. We show that, depending on some parameters of the system like the radius and mass of the planet, eccentricity etc., the rotation can be strongly perturbed and a chaotic layer around the synchronous state may occupy a significant region of the phase space. 55 Cnc e is an example.     

Refinement of the phase of a spherical planet located at a small distance from the Sun

Formulas for refining the phase of a spherical planet located a small distance from the Sun are derived. Finite heliocentric distance of the planet results in the formation on its visible disk of the geometric terminator, which is not coincident with the orthographic terminator. The visible disk is assumed to be observed from Earth in orthographic projection. We suggest introducing linear and surface phases for the geometric terminator in accordance with two existing definitions of the phase of a planet. Linear and surface phases of a planet are shown to be given by different sets of formulas. An example of the computation of the phase of Mercury is given.     

Capture of grains into resonances through Poynting-Robertson drag

We review here some relevant problems connected to the evolution of circumstellar dust grains, subjected to Poynting-Robertson (PR) drag, and perturbed by first-order resonances with a planet on a circular orbit. We show that only outer mean motion resonances are able to counteract the damping effect of PR drag. However, the high orbital eccentricities reached by the particle lead to orbit crossings with the planet. This is a serious difficulty for a permanent trapping to be achieved. In any case, we show that the time spent in the resonance is long enough for statistical effects (accumulation at the resonant radius) to be significant. We underline some difficulties associated with this problem, namely, the non-adiabaticity of motion in the resonance phase space and the existence of close encounters with the planet at high eccentricities.     

Investigation of some of the principal geometric effects on planetary polarization

Three major geometric factors which are likely to influence theoretical interpretation of planetary polarization measurements, viz., observer—planet distance, horizontal inhomogeneity of planetary disk, and deviation from a spherical body, are investigated.The distance effect is examined for regional as well as global polarizations. For convenience of analysis, the expressions for zenith and azimuth angles of incident and emergent light appropriate for a snap-shot observation are derived as explicit functions of distance between observer and planet. Sample computations for Venus indicate that regional polarization near the planetary limb is significantly affected by the observer's distance. This effect should be particularly noticeable when an observation is made at a phase angle around which the single scattering polarization of atmospheric scattering agents exhibits a steep variation. The global polarization at large phase angles (measured at disk-center) is gradually moved toward smaller phase angles, as the observer approaches the planet. Any narrow polarization features such as rainbow and glory at small phase angles are heavily smoothed out.The effects of horizontal inhomogeneity are investigated with a planetary disk having highly polarizing regions at high latitudes. Comparison of theoretical global polarization computed for such a disk with the Pioneer Venus OCPP measurements shows a possible change in cloud-haze stratification approximately at 50° latitude, consistent with other imaging observations. An approximate analytical representation of residual polarization at zero phase angle is then derived to compare to the numerical results for Venus. An attempt is also made to explain the relatively large magnitude of residual polarization observed on Jupiter.Finally, to study the effects of nonsphericity of planetary body, the global polarizations are computed for a spheroidal planet. The global polarization tends to increase as the planet's oblateness increases. However, for Jupiter and Saturn, such effect may be of secondary importance.     

Allowance for the phase of Jupiter when determining the contact times in some phenomena of its Galilean satellites

We suggest a new method for predicting the phenomena observed in Jovian system of Galilean satellites that takes into account the planet’s phase effect. The method allows one to determine the geocentric times of the contacts of the satellite and its shadow with the illuminated part of the planet’s visible disk that occur near its inferior geocentric and inferior heliocentric conjunctions, respectively. The calculation is performed in the orthographic approximation for the geometric center of the satellite and its shadow by taking into account the curvature of the satellite’s orbit and the visible flattening of Jovian disk. The correction for the phase to the satellite’s contact time is determined from the phase shift of the center of the planet’s disk.     

The impact of comet SL9 with Jupiter: IR observations at TIRGO

Infrared images of Jupiter have been obtained on 5 nights before, during and shortly after the period of the impacts of the fragments of comet Shoemaker-Levy 9 (1993e) with the giant planet. Long lived bright spots produced by the impacts have been measured and analyzed. By measuring the intensity variation of the spots as a function of Jupiter rotation we show that these spots are likely constituted by large and thin clouds of dust located above the methane layer. The IR relative albedos has been also measured for some of these spots.On leave from Center for Astrophysics - Cambridge (USA)     

The origin of dark inclusions in Allende: New evidence from lithium isotopes

Abstract— Aqueous and thermal processing of primordial materials occurred prior to and during planet formation in the early solar system. A record of how solid materials were altered at this time is present in the carbonaceous chondrites, which are naturally delivered fragments of primitive asteroids. It has been proposed that some materials, such as the clasts termed “dark inclusions” found in type III chondrites, suggest a sequence of aqueous and thermal events. Lithium isotopes (⁶Li and ⁷Li) can reveal the role of liquid water in dark inclusion history. During aqueous alteration, ⁷Li passes preferentially into solution leaving ⁶Li behind in the solid phase and, consequently, any relatively extended periods of interaction with ⁷Li‐rich fluids would have left the dark inclusions enriched in the heavier isotope when compared to the meteorite as a whole. Our analyses of lithium isotopes in Allende and its dark inclusions reveal marked isotopic homogeneity and no evidence of greater levels of aqueous alteration in dark inclusion history.     

Habitable zones of host stars during the post-MS phase

A star will become brighter and brighter with stellar evolution, and the distance of its habitable zone will become larger and larger. Some planets outside the habitable zone of a host star during the main sequence phase may enter the habitable zone of the host star during other evolutionary phases. A terrestrial planet within the habitable zone of its host star is generally thought to be suitable for the existence of life. Furthermore, a rocky moon around a giant planet may be also suitable for life to survive, provided that the planet–moon system is within the habitable zone of its host star. Using Eggleton’s code and the boundary flux of the habitable zone, we calculate the habitable zone of our Solar system after the main sequence phase. It is found that Mars’ orbit and Jupiter’s orbit will enter the habitable zone of the Solar system during the subgiant branch phase and the red giant branch phase, respectively. And the orbit of Saturn will enter the habitable zone of Solar during the He-burning phase for about 137 million years. Life is unlikely at any time on Saturn, as it is a giant gaseous planet. However, Titan, the rocky moon of Saturn, may be suitable for biological evolution and become another Earth during that time. For low-mass stars, there are similar habitable zones during the He-burning phase as our Solar, because there are similar core masses and luminosities for these stars during that phase.     

Asteroidal jet streams

Detailed studies of the asteroidal belt are of importance for clarifying whether the asteroids are fragments of an exploded planet or represent an intermediate state in the accretion of planets.A study of the Hirayama family Flora shows that it contains three groups of bodies travelling in almost identical orbits, thus constituting three jet streams.It is shown that the formation of jet streams is difficult to reconcile with the exploded planet view.In order to decide whether the formation of jet streams is a general phenomenon in the asteroidal belt, it is necessary to determine the orbits of a large number of small asteroids.     

Deciphering the origin of the regular satellites of gaseous giants - Iapetus: The Rosetta ice-moon

Ever since their discovery the regular satellites of Jupiter and Saturn have held out the promise of providing an independent set of observations with which to test theories of planet formation. Yet elucidating their origins has proven elusive. Here we show that Iapetus can serve to discriminate between satellite formation models. Its accretion history can be understood in terms of a two-component gaseous subnebula, with a relatively dense inner region, and an extended tail out to the location of the irregular satellites, as in the SEMM model of Mosqueira and Estrada (2003a,b) (Mosqueira, I., Estrada, P.R. [2003a]. Icarus 163, 198-231; Mosqueira, I., Estrada, P.R. [2003b]. Icarus 163, 232-255). Following giant planet formation, planetesimals in the feeding zone of Jupiter and Saturn become dynamically excited, and undergo a collisional cascade. Ablation and capture of planetesimal fragments crossing the gaseous circumplanetary disks delivers enough collisional rubble to account for the mass budgets of the regular satellites of Jupiter and Saturn. This process can result in rock/ice fractionation as long as the make up of the population of disk crossers is non-homogeneous, thus offering a natural explanation for the marked compositional differences between outer solar nebula objects and those that accreted in the subnebulae of the giant planets. For a given size, icy objects are easier to capture and to ablate, likely resulting in an overall enrichment of ice in the subnebula. Furthermore, capture and ablation of rocky fragments become inefficient far from the planet for two reasons: the gas surface density of the subnebula is taken to drop outside the centrifugal radius, and the velocity of interlopers decreases with distance from the planet. Thus, rocky objects crossing the outer disks of Jupiter and Saturn never reach a temperature high enough to ablate either due to melting or vaporization, and capture is also greatly diminished there. In contrast, icy objects crossing the outer disks of each planet ablate due to the melting and vaporization of water-ice. Consequently, our model leads to an enhancement of the ice content of Iapetus, and to a lesser degree those of Titan, Callisto and Ganymede, and accounts for the (non-stochastic) compositions of these large, low-porosity outer regular satellites of Jupiter and Saturn. For this to work, the primordial population of planetesimals in the Jupiter-Saturn region must be partially differentiated, so that the ensuing collisional cascade produces an icy population of ?1 m size fragments to be ablated during subnebula crossing. We argue this is likely because the first generation of solar nebula ∼10 km planetesimals in the Jupiter-Saturn region incorporated significant quantities of ²⁶Al. This is the first study successfully to provide a direct connection between nebula planetesimals and subnebulae mixtures with quantifiable and observable consequences for the bulk properties of the regular satellites of Jupiter and Saturn, and the only explanation presently available for Iapetus’ low density and ice-rich composition.     

Planetesimal growth in turbulent discs before the onset of gravitational instability

It is difficult to imagine a planet formation model that does not at some stage include a gravitationally unstable disc. Initially unstable gas–dust discs may form planets directly, but the high surface density required has motivated the alternative that gravitational instability occurs in a dust sub-layer only after grains have grown large enough by electrostatic sticking. Although such growth up to the instability stage is efficient for laminar discs, previous research concluded that realistic disc turbulence catastrophically increases the settling time, thereby requiring additional processes to facilitate planet formation on the needed time scales. We develop a different model for the influence of turbulence on the collisional velocity of grains and on the scale height of the dust layer and find that the earlier conclusions must be revisited. The model produces a disc-radius dependent time scale to reach a gravitationally unstable phase of planet formation. For a range of dust sticking and disc parameters, we find that for viscosity parameters α ⩽ 10⁻³, this time scale is short enough over a significant range in radii R that turbulence does not catastrophically slow the early phases of planet formation, even in the absence of agglomeration enhancement agents like vortices.     

Calculations of the effects of angular momentum on the early evolution of Jupiter

The planet Jupiter is assumed to have formed as a subcondensation in the solar nebula. The initial phase of its evolution is one of hydrostatic contraction with radiative energy transport. Calculations of evolutionary sequences through this phase are presented, including the effects of angular momentum. The calculations are carried out in two space dimensions under the assumptions of axial symmetry, constancy of angular velocity on cylindrical surfaces about the rotation axis, a pressure-density relation given by the polytrope of index 3, conservation of angular momentum, and a homogeneous composition. The results show that under certain physically reasonable initial distributions of density and angular momentum the formation of a central planet and a rotating circumplanetary envelope is possible, while under assumptions a point of instability is reached that probably results in the breakup of the condensation by fission into two or more parts. The models are discussed with reference to the present angular momenta of Jupiter and its regular satellites.     

On the co-orbital motion in the planar restricted three-body problem: the quasi-satellite motion revisited

In the framework of the planar and circular restricted three-body problem, we consider an asteroid that orbits the Sun in quasi-satellite motion with a planet. A quasi-satellite trajectory is a heliocentric orbit in co-orbital resonance with the planet, characterized by a nonzero eccentricity and a resonant angle that librates around zero. Likewise, in the rotating frame with the planet, it describes the same trajectory as the one of a retrograde satellite even though the planet acts as a perturbator. In the last few years, the discoveries of asteroids in this type of motion made the term “quasi-satellite” more and more present in the literature. However, some authors rather use the term “retrograde satellite” when referring to this kind of motion in the studies of the restricted problem in the rotating frame. In this paper, we intend to clarify the terminology to use, in order to bridge the gap between the perturbative co-orbital point of view and the more general approach in the rotating frame. Through a numerical exploration of the co-orbital phase space, we describe the quasi-satellite domain and highlight that it is not reachable by low eccentricities by averaging process. We will show that the quasi-satellite domain is effectively included in the domain of the retrograde satellites and neatly defined in terms of frequencies. Eventually, we highlight a remarkable high eccentric quasi-satellite orbit corresponding to a frozen ellipse in the heliocentric frame. We extend this result to the eccentric case (planet on an eccentric motion) and show that two families of frozen ellipses originate from this remarkable orbit.     

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.     

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)     

Theoretical interpretation of photometric properties of the Martian surface and atmosphere

Earth-based UBV photometry, high-quality photographs from the Lowell Observatory collection, and Mariner 9 data have been combined with a new radiative transfer theory to derive physical parameters for the Martian surface and atmosphere, both before and during the 1971 dust storm. We find that the dust particles of the storm had a single-scattering albedo of 0.84 ± 0.02 and an asymmetry factor of 0.35 ± 0.10 in green (V) light. The geometric albedo of Mars was 0.15 and the phase integral 1.83, which yield 0.27 for the Bond albedo. The mean optical thickness of the “clear” atmosphere averaged over the whole planet was 0.15 ± 0.05 and was not detectably dependent on wavelength. Geometric albedos for the surface are 0.25 (light areas) and 0.17 (dark areas) in V, 0.095 in B (both areas), and 0.060 in U (both areas). The soil particles are moderately backward scattering with an asymmetry factor of ?0.20, indicating them to be rather opaque. The mean surface roughness, on a scale larger than that of individual dust particles and therefore large compared with the wavelength, is 0.57. This represents the depth/radius ratio of an average hole and it is only one-half as large as values typical for the Moon and asteroids.     

Dynamical relaxation and the orbits of low-mass extrasolar planets

We consider the evolution of a system containing a population of massive planets formed rapidly through a fragmentation process occurring on a scale on the order of 100 au and a lower mass planet that assembles in a disc on a much longer time-scale. During the formation phase, the inner planet is kept on a circular orbit owing to tidal interaction with the disc, while the outer planets undergo dynamical relaxation. Interaction with the massive planets left in the system after the inner planet forms may increase the eccentricity of the inner orbit to high values, producing systems similar to those observed.     

Consequences of tidal interaction between disks and orbiting protoplanets for the evolution of multi-planet systems with architecture resembling that of Kepler 444

We study orbital evolution of multi-planet systems with masses in the terrestrial planet regime induced through tidal interaction with a protoplanetary disk assuming that this is the dominant mechanism for producing orbital migration and circularization. We develop a simple analytic model for a system that maintains consecutive pairs in resonance while undergoing orbital circularization and migration. This model enables migration times for each planet to be estimated once planet masses, circularization times and the migration time for the innermost planet are specified. We applied it to a system with the current architecture of Kepler 444 adopting a simple protoplanetary disk model and planet masses that yield migration times inversely proportional to the planet mass, as expected if they result from torques due to tidal interaction with the protoplanetary disk. Furthermore the evolution time for the system as a whole is comparable to current protoplanetary disk lifetimes. In addition we have performed a number of numerical simulations with input data obtained from this model. These indicate that although the analytic model is inexact, relatively small corrections to the estimated migration rates yield systems for which period ratios vary by a minimal extent. Because of relatively large deviations from exact resonance in the observed system of up to 2 %, the migration times obtained in this way indicate only weak convergent migration such that a system for which the planets did not interact would contract by only \({\sim }1\,\%\) although undergoing significant inward migration as a whole. We have also performed additional simulations to investigate conditions under which the system could undergo significant convergent migration before reaching its final state. These indicate that migration times have to be significantly shorter and resonances between planet pairs significantly closer during such an evolutionary phase. Relative migration rates would then have to decrease allowing period ratios to increase to become more distant from resonances as the system approached its final state in the inner regions of the protoplanetary disk.     

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.     

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

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