Chaotic gravitational zones around a regularly shaped complex rotating body

In preparation for the Rosetta mission, the location and widths of gravitational resonances surrounding a regularly shaped and possibly complex rotating body are mapped following the second fundamental model of resonance. It is found that for uniaxial rotation of the central body, the surrounding resonances are widest for prograde orbits. If the figure axis is tilted with respect to the spin axis of the central body, an additional number of wide resonances appear with a preference for prograde and inclined orbits, and the occurrence of initial conditions which lie in the globally connected chaotic web is significantly increased. For larger rotational excitations, it is seen how these new additional resonances overlap internally at low eccentricity for very large semi-major axes. However, with exceptions for some excited short-axis rotational modes of the central body, it is argued that most resonances vanish for retrograde orbits lying in the plane normal to the body spin, and that resonant or non-resonant stability therefore can be expected for a wide range of mean orbit eccentricities.     

Long term evolution and chaotic diffusion of the insolation quantities of Mars

As the obliquity of Mars is strongly chaotic, it is not possible to give a solution for its evolution over more than a few million years. Using the most recent data for the rotational state of Mars, and a new numerical integration of the Solar System, we provide here a precise solution for the evolution of Mars' spin over 10 to 20 Myr. Over 250 Myr, we present a statistical study of its possible evolution, when considering the uncertainties in the present rotational state. Over much longer time span, reaching 5 Gyr, chaotic diffusion prevails, and we have performed an extensive statistical analysis of the orbital and rotational evolution of Mars, relying on Laskar's secular solution of the Solar System, based on more than 600 orbital and 200,000 obliquity solutions over 5 Gyr. The density functions of the eccentricity and obliquity are specified with simple analytical formulas. We found an averaged eccentricity of Mars over 5 Gyr of 0.0690 with standard deviation 0.0299, while the averaged value of the obliquity is 37.62° with a standard deviation of 13.82°, and a maximal value of 82.035°. We find that the probability for Mars' obliquity to have reached more than 60° in the past 1 Gyr is 63.0%, and 89.3% in 3 Gyr. Over 4 Gyr, the position of Mars' axis is given by a uniform distribution on a spherical cap limited by the obliquity 58.62°, with the addition of a random noise allowing a slow diffusion beyond this limit. We can also define a standard model of Mars' insolation parameters over 4 Gyr with the most probable values 0.068 for the eccentricity and 41.80° for the obliquity.     

Orbit and spin evolution of synchronous binary stars on the main sequence(a theoretical improvement to the analytical method)

This paper provides a method to study the solution of equations for synchronous binary stars with large eccentricity on the main sequence.The theoretical results show that the evolution of the eccentricity is linear with time or follows an exponential form,and the semi-major axis and spin vary with time in an exponential form that are different from the results given in a previous paper.The improved method is applicable in both cases of large eccentricity and small eccentricity.In addition,the number of terms in the expansion of a series with small eccentricity is very long due to the series converging slowly.The advantage of this method is that it is applicable to cases with large eccentricity due to the series converging quickly.This paper chooses the synchronous binary star V1143 Cyg that is on the main sequence and has a large eccentricity(e = 0.54) as an example calculation and gives the numerical results.Lastly,the evolutionary tendency including the evolution of orbit and spin,the time for the speed up of spin,the circularization time,the orbital collapse time and the life time are given in the discussion and conclusion.The results shown in this paper are an improvement on those from the previous paper.     

Lunar shape does not record a past eccentric orbit

The Moon has long been known to have an overall shape not consistent with expected past tidal forces. It has recently been suggested (Garrick-Bethell, I., Wisdom, J., Zuber, M.T. [2006]. Science 313, 652-655) that the present lunar moments of inertia indicate a past high-eccentricity orbit and, possibly, a past non-synchronous spin-orbit resonance. Here I show that the match between the lunar shape and the proposed orbital and spin states is much less conclusive than initially proposed. Garrick-Bethell et al. (Garrick-Bethell, I., Wisdom, J., Zuber, M.T. [2006]. Science 313, 652-655) spin and shape evolution scenarios also completely ignore the physics of the capture into such resonances, which require prior permanent deformation, as well as tidal despinning to the relevant resonance. If the early lunar orbit was eccentric, the Moon would have been rotating at an equilibrium non-synchronous rate determined by it eccentricity. This equilibrium supersynchronous rotation would be much too fast to allow a synchronous spin-orbit lock at e = 0.49, while the capture into the 3:2 resonance is possible only for a very constrained lunar eccentricity history and assuming some early permanent lunar tri-axiality. Here I show that large impacts in the early history of the Moon would have frequently disrupted this putative resonant rotation, making the rotation and eccentricity solutions of Garrick-Bethell et al. (Garrick-Bethell, I., Wisdom, J., Zuber, M.T. [2006]. Science 313, 652-655) unstable. I conclude that the present lunar shape cannot be used to support the hypothesis of an early eccentric lunar orbit.     

Dynamics and gravitational waveforms in spinning compact eccentric binaries

The effects of eccentricity on the Hamiltonian dynamics of post-Newtonian spinning compact binaries and gravitational radiation from eccentric orbits are discussed. The simulation results of scans for chaos show that the eccentricity has a great effect on the dynamics without considering dissipation due to gravitational radiation. Moreover, both the dynamics behavior and the orbital eccentricity jointly modulate the gravitational waveforms, and the spin–spin coupling effect play an important role during the gravitational radiation of inspiral and coalescence. Moreover, the imprint of characteristic of the original system can be deduced from the time-domain and frequency-domain waveforms.     

Orbital evolution of small binary asteroids

About 15% of both near-Earth and main-belt asteroids with diameters below 10 km are now known to be binary. These small asteroid binaries are relatively uniform and typically contain a fast-spinning, flattened primary and a synchronously rotating, elongated secondary that is 20-40% as large (in diameter) as the primary. The principal formation mechanism for these binaries is now thought to be YORP (Yarkovsky-O’Keefe-Radzievskii-Paddack) effect induced spin-up of the primary followed by mass loss and accretion of the secondary from the released material. It has previously been suggested (?uk, M. [2007]. Astrophys. J. 659, L57-L60) that the present population of small binary asteroids is in a steady state between production through YORP and destruction through binary YORP (BYORP), which should increase or decrease secondary’s orbit, depending on the satellite’s shape. However, BYORP-driven evolution has not been directly modeled until now. Here we construct a simple numerical model of the binary’s orbital as well the secondary’s rotational dynamics which includes BYORP and selected terms representing main solar perturbations. We find that many secondaries should be vulnerable to chaotic rotation even for relatively low-eccentricity mutual orbits. We also find that the precession of the mutual orbit for typical small binary asteroids might be dominated by the perturbations from the prolate and librating secondary, rather than the oblate primary. When we evolve the mutual orbit by BYORP we find that the indirect effects on the binary’s eccentricity (through the coupling between the orbit and the secondary’s spin) dominate over direct ones caused by the BYORP acceleration. In particular, outward evolution causes eccentricity to increase and eventually triggers chaotic rotation of the secondary. We conclude that the most likely outcome will be reestablishing of the synchronous lock with a “flipped” secondary which would then evolve back in. For inward evolution we find an initial decrease of eccentricity and secondary’s librations, to be followed by later increase. We think that it is likely that various forms of dissipation we did not model may damp the secondary’s librations close to the primary, allowing for further inward evolution and a possible merger. We conclude that a merger or a tidal disruption of the secondary are the most likely outcomes of the BYORP evolution. Dissociation into heliocentric pairs by BYORP alone should be very difficult, and satellite loss might be restricted to the minority of systems containing more than one satellite at the time.     

On the rotation of co-orbital bodies in eccentric orbits

We investigate the resonant rotation of co-orbital bodies in eccentric and planar orbits. We develop a simple analytical model to study the impact of the eccentricity and orbital perturbations on the spin dynamics. This model is relevant in the entire domain of horseshoe and tadpole orbit, for moderate eccentricities. We show that there are three different families of spin–orbit resonances, one depending on the eccentricity, one depending on the orbital libration frequency, and another depending on the pericenter’s dynamics. We can estimate the width and the location of the different resonant islands in the phase space, predicting which are the more likely to capture the spin of the rotating body. In some regions of the phase space the resonant islands may overlap, giving rise to chaotic rotation.     

Secular resonance,solar spin down,and the orbit of Mercury

A mechanism capable of accounting for the large mean eccentricity (0.175) and inclination (7°.2) of Mercury is discussed. Provided the gravitational field of the rapidly rotating primordial Sun had a sufficiently large second degree harmonic (i.e., J₂ ? order 10^?3), subsequent solar spin down would drive the orbit of Mercury through two secular resonances with Venus, one involving the precession of the line of apsides, the other one involving the regression of the nodal line. Resonance passage generates contributions to the eccentricity and inclination that are proportional to the square root of the characteristic solar spin down time. We find that an initial solar rotation l period of $\text{P ? 5}\text{1}\text{2}\text{hr}$ guarantees passage through resonance and that a spin down time of $\text{τ = Ω|dΩ/dt|}{}^{\mathrm{?1}}$ of order 10⁶ years could have produced the observed eccentricity and inclination. Such a primordial rotation rate is comparable to the measured rotations of very young stars and the spin down time appears consistent with the time scale derived for magnetic braking of the Sun's rotation by an intense solar wind during a T-Tauri stage of solar evolution.     

Behavior of the Rayleigh-Taylor mode in a dusty plasma with rotational and shear flows

The stability of a dusty plasma with sheared rotational flows is investigated. Using the fluid model together with the Bayly nonmodal approach, the inhomogeneous partial differential equations governing short-wavelength perturbations at the center of a rotational flow field or vortex structure are obtained. The effects of flow eccentricity, strength of the flow shear, as well as concentration of dust grains on the stability of the perturbations are investigated numerically. It is found that flow shear can cause secondary Rayleigh-Taylor instability of a rotational flow.     

Rotational Properties of Cometary Nuclei

We review several techniques used to retrieve rotational parameters from observations. The spin period of a dozen of comets retrieved with these techniques are summarized. We describe how the spin period of comet Hale-Bopp (C/1995 O1) has been calculated with a high accuracy (11.30–11.34 h). Although several authors converged to a spin axis orientation at (α,δ) = (275 ± 15°, -55 ± 5°), detailed studies indicate that the dust jets morphology in 1996–1997 may be incompatible with this orientation. Comet 19P/Borrelly has been recently observed by the Deep Space 1 spacecraft. At the same time, its spin axis orientation and period have been determined by several authors to be respectively (α,δ) = (225 ± 15°, -10 ± 10°)and 26h. These two comets are likely to be in (or close to) a principal axis spin state. We discuss new modeling of the spin state of comet 46P/Wirtanen, the target of the Rosetta mission. The model involves a three-dimensional shape and thermal model, from which the torque of the non gravitational force is calculated at each time step. The moments of inertia are computed for each irregular shape. The results from numerical integrations show that this comet can remain in a principal axis spin state during more than 10 orbits if the spin period does not get above～6 h. If the spin period increases, its nucleus gets rapidly into excited spin states. It shows that even small and very active short-period comets are not necessarily in non principal axis spin states. In the last section, the consequences of recent observations and modeling of the rotational parameters of comet nuclei are discussed, and unsolved problems are presented.     

Supernova kicks and misaligned Be star binaries

Be stars are rapidly spinning B stars surrounded by an outflowing disc of gas in Keplerian rotation. Be star/X-ray binary systems contain a Be star and a neutron star. They are found to have non-zero eccentricities and there is evidence that some systems have a misalignment between the spin axis of the star and the spin axis of the binary orbit. The eccentricities in these systems are caused by a kick to the neutron star during the supernova that formed it. Such kicks would also give rise to misalignments. In this paper, we investigate the extent to which the same kick distribution can give rise to both the observed eccentricity distribution and the observed misalignments. We find that a Maxwellian distribution of velocity kicks with a low velocity dispersion,  σ_k≈ 15 km s⁻¹  , is consistent with the observed eccentricity distribution but is hard to reconcile with the observed misalignments, typically   i ≥ 25°  . Alternatively, a higher velocity kick distribution,  σ_k= 265 km s⁻¹  , is consistent with the observed misalignments but not with the observed eccentricities, unless post-supernova circularization of the binary orbits has taken place. We discuss briefly how this might be achieved.     

The equilibrium configuration of rapid rotating compact stars and some gravitational effects

In this paper, the equilibrium configurations of rapid rotating compact stars and some gravitational effects are studied within the general relativity by use of the Harrison-Wheeler equation of state and by the self-consistent field method. Numerical calculations show that the equilibrium configuration of a rotating star is a spheroid. For large spin velocities, say, ω > 3.0 × 10² sec⁻¹ the eccentricity and mass increase very rapidly as the angular velocity increases, for the critical angular velocity of the rotating star, the eccentricity is about 0.7, the increase in mass is about 10–35%. The difference of the gravitational redshifts at the surface of the star caused by rotation, and the difference of the light bending when the beam moves in the direction of rotation or in the opposite direction are obvious.     

Non-gravitational force modeling of Comet 81P/Wild 2: II. Rotational evolution

In this paper, we have studied both the dynamical and the rotational evolution of an 81P/Wild 2-like comet under the effects of the outgassing-induced force and torque. The main aim is to study if it is possible to reproduce the non-gravitational orbital changes observed in this comet, and to establish the likely evolution of both orbital and rotational parameters. To perform this study, a simple thermophysical model has been used to estimate the torque acting on the nucleus. Once the torque is calculated, Euler equations are solved numerically considering a nucleus mass directly estimated from the changes in the orbital elements (as determined from astrometry). According to these simulations, when the water production rate and changes in orbital parameters for 1997, as well as observational rotational parameters for 2004 are imposed as constraints, the change in the orbital period of 81P/Wild 2, , will decrease so that to , which is similar to the actual tendency observed from 1988 up to 1997. This nearly constant decreasing can be explained as due to a slight drift of the spin axis orientation towards larger ecliptic longitudes. After studying the possible spin axis orientations proposed for 1997, simulations suggest that the spin obliquity and argument (I,Φ)=(56°,167°) is the most likely. As for rotational evolution, changes per orbit smaller than 10% of the actual spin velocity are probable, while the most likely value corresponds to a change between 2 and 7% of the spin velocity. Equally, net changes in the spin axis orientation of 4°-8° per orbit are highly expected.     

Invariant rotational curves in Sitnikov's Problem

The Sitnikov's Problem is a Restricted Three-Body Problem of Celestial Mechanics depending on a parameter, the eccentricity,e. The Hamiltonian,H(z, v, t, e), does not depend ont ife=0 and we have an integrable system; ife is small the KAM Theory proves the existence of invariant rotational curves, IRC. For larger eccentricities, we show that there exist two complementary sequences of intervals of values ofe that accumulate to the maximum admissible value of the eccentricity, 1, and such that, for one of the sequences IRC around a fixed point persist. Moreover, they shrink to the planez=0 ase tends to 1.     

PSR J1753−2240: a mildly recycled pulsar in an eccentric binary system

We report the discovery of PSR J1753−2240 in the Parkes Multibeam Pulsar Survey data base. This 95-ms pulsar is in an eccentric binary system with a 13.6-d orbital period. Period derivative measurements imply a characteristic age in excess of 1 Gyr, suggesting that the pulsar has undergone an episode of accretion-induced spin-up. The eccentricity and spin period are indicative of the companion being a second neutron star, so that the system is similar to that of PSR J1811−1736, although other companion types cannot be ruled out at this time. The companion mass is constrained by geometry to lie above 0.48 solar masses, although long-term timing observations will give additional constraints. If the companion is a white dwarf or a main-sequence star, optical observations may yield a direct detection of the companion. If the system is indeed one of the few known double neutron star systems, it would lie significantly far from the recently proposed spin-period/eccentricity relationship.     

Extension of Cassini's Laws

We investigate the Cassini's laws which describe the rotational motion in a 1:1 spin-orbit resonance. When this rotational motion follows the conventional Cassini's laws, the figure axis coincides with the angular momentum axis. In this case we underline the differences between the rotational Hamiltonian for a 'slow rotating' body like the Moon and for a 'fast rotating' body like Phobos. Then, we study a more realistic rotational Hamiltonian where the angle J between the figure axis and the angular momentum axis could be different from zero. This Hamiltonian has not been studied before. We have found a new particular solution for this Hamiltonian which could be seen as an extension of the Cassini's laws. In this new solution the angle J is constant, which is not zero, and the precession of the angular momentum plane is equal to the mean motion of the argument of pericenter of the rotating body. This type of rotational motion is only possible when the orbital eccentricity of the rotating body is not zero. This new law enables describing in particular, the Moon mean rotational motion for which the mean value of the angle J is found to be equal to 103.9±0.7 s of arc.     

Rotational dynamics of subsolar sublimating triaxial comets

A model with subsolar water sublimation on a triaxial, ellipsoidal comet nucleus is presented for the calculation of reactive torques. The resulting differential equations describing the comet's rotation are then Hamiltonian, and gravity-gradients are trivial to include. While effects derived from a weak perturbing function are neither able to change the rotational excitation nor the spin magnitude of the nucleus, it is shown how the spin orientation of comets can change significantly over an orbital run. However, of the four comets studied, 1P, 19P, 46P and the Rosetta target 67P, 19P and 46P were the only objects clearly exhibiting this feature, thereby confirming a technique used to derive the consequences of a more elaborate model of sublimation induced torques. In particular, the rotational parameters of 67P were seen to be very stable, indicating that a highly kinematical model of its rotation for the mapping of the comet's gravitational field during the Rosetta mission can be used. The model's hierarchy with 1P/Halley as the object with highest excitation probability, is consistent with observations.     

Tidal locking of habitable exoplanets

Potentially habitable planets can orbit close enough to their host star that the differential gravity across their diameters can produce an elongated shape. Frictional forces inside the planet prevent the bulges from aligning perfectly with the host star and result in torques that alter the planet’s rotational angular momentum. Eventually the tidal torques fix the rotation rate at a specific frequency, a process called tidal locking. Tidally locked planets on circular orbits will rotate synchronously, but those on eccentric orbits will either librate or rotate super-synchronously. Although these features of tidal theory are well known, a systematic survey of the rotational evolution of potentially habitable exoplanets using classic equilibrium tide theories has not been undertaken. I calculate how habitable planets evolve under two commonly used models and find, for example, that one model predicts that the Earth’s rotation rate would have synchronized after 4.5 Gyr if its initial rotation period was 3 days, it had no satellites, and it always maintained the modern Earth’s tidal properties. Lower mass stellar hosts will induce stronger tidal effects on potentially habitable planets, and tidal locking is possible for most planets in the habitable zones of GKM dwarf stars. For fast-rotating planets, both models predict eccentricity growth and that circularization can only occur once the rotational frequency is similar to the orbital frequency. The orbits of potentially habitable planets of very late M dwarfs ( Open image in new window ) are very likely to be circularized within 1 Gyr, and hence, those planets will be synchronous rotators. Proxima b is almost assuredly tidally locked, but its orbit may not have circularized yet, so the planet could be rotating super-synchronously today. The evolution of the isolated and potentially habitable Kepler planet candidates is computed and about half could be tidally locked. Finally, projected TESS planets are simulated over a wide range of assumptions, and the vast majority of potentially habitable cases are found to tidally lock within 1 Gyr. These results suggest that the process of tidal locking is a major factor in the evolution of most of the potentially habitable exoplanets to be discovered in the near future.     

Periodic perturbations in close binary systems due to rotational distortion

This paper deals with the perturbations which rotation can produce to the orbital elements of a close binary system. The rectangular components,R, S andW of the disturbing accelerations due to rotation have been substituted to the Gauss form of Lagrange's planetary equations to yield the first order approximation. The results obtained are exact for any value of orbital eccentricity between the values 0<e<1 and for arbitrary inclinations of the rotational axes to the orbital plane.First and second order approximations are given for the special case when equators are coplanar to the orbit.     

Spin Axis Direction of Comet 19p/Borrelly Based on Observations from 2000 and 2001

We calculate the direction of the rotational angular momentum vector,M, of comet 19P/Borrelly based on rotational lightcurve data from 2000, groundbased imaging of the coma during the Deep Space 1encounter, and the basic near-nucleus coma morphology as revealed by the Deep Space 1 spacecraft. For the most likely direction, we derivea family of solutions (with center at RA = 221°, Dec = -7°) if the direction of M is towards the sunward hemisphere during the Deep Space 1 encounter, whereas if the rotation is of opposite sense, the diametrically opposite family of solutions (with center at RA = 41°, Dec = 7°) would result. We argue that the coma morphology in September 2001 is consistent with the nucleus being a principal axis rotator or one observationally indistinguishable from it. Therefore, for all practical purposes, the direction of the rotational angular momentum vector coincides with the spin axis. We also discuss why the determination of the spin axis direction based on observations from the last apparition is in disagreement with the current result.