Empirical evidence for tidal evolution in transiting planetary systems

Most transiting planets orbit very close to their parent star, causing strong tidal forces between the two bodies. Tidal interaction can modify the dynamics of the system through orbital alignment, circularization, synchronization and orbital decay by exchange of angular moment. Evidence for tidal circularization in close-in giant planet is well known. Here, we review the evidence for excess rotation of the parent stars due to the pull of tidal forces towards spin-orbit synchronization. We find suggestive empirical evidence for such a process in the present sample of transiting planetary systems. The corresponding angular momentum exchange would imply that some planets have spiralled towards their star by substantial amounts since the dissipation of the protoplanetary disc. We suggest that this could quantitatively account for the observed mass–period relation of close-in gas giants. We discuss how this scenario can be further tested and point out some consequences for theoretical studies of tidal interactions and for the detection and confirmation of transiting planets from radial velocity and photometric surveys.     

Instability of the 2+2 body problem

The equations of motion of the 2+2 body problem (two interacting particles in the gravitational field of two much more massive primaries m₁ and m₂ in circular keplerian orbit) have an integral analogous to the Jacobi integral of the circular 2+1 body problem. We show here that with 2+2 bodies this integral does not give rise to Hill stability, i.e. to confinement for all time in a portion of the configuration space not allowing for some close approaches to occur. This is because all the level manifolds are connected and all exchanges of bodies between the regions surroundingm ₁,m ₂ and infinity do not contradict the conservation of the integral. However, it is worth stressing that some of these exchanges are physically meaningless, because they involve either unlimited extraction of potential energy from the binary formed by the small bodies (without taking into account their physical size) or significant mutual perturbations between the small masses without close approach, a process requiring, for the Sun-Jupiter-two asteroids system, timescales longer than the age of the Solar System.     

Tidal torque constants as a critical test for an evolution of close eclipsing binaries

Individual tidal torque λ ₂,E ₂ and apsidal-motion k ₂ constants were calculated for 112 close eclipsing binaries (CEBs) with Detached components belonging to the Main Sequence (DMS-type) from the catalogue by Svechnikov and Perevozkina (Catalogue of orbital elements, masses and luminosities of variable stars of DMS-type and some results of its statistical treatment, Ural State University Press, Yekaterinburg, pp. 1–5, 1999) and for 95 detached binaries taken from the catalogue by Torres et al. (Astron. Astrophys. Rev. 18:67, 2010) on the base of theoretical evolutionary stellar models including tidal torque constants by Claret (Astron. Astrophys. 424:919, 2004). A method of the inversion of model track grid into isochrones was formulated as a complex interpolation procedure for DMS-binaries data. Sets of isochrones were computed in k ₂–M, k ₂–R, λ ₂–M, λ ₂–R, E ₂–M, and E ₂–R planes. Calculated tidal torque constants allow to test stellar structure theory by comparing observed and estimated values of apsidal motion period and analyzing the correlation between timescales of synchronization, circularization, magnetic braking, as well as nuclear burning of DMS-components.     

Gravitational radiation from approaching neutron stars and the rotational explosion mechanism for collapsing supernovae

We consider the evolution of a neutron star binary system under the effect of two factors: gravitational radiation and mass transfer between the components. Gravitational radiation is specified under the justified assumption of a circular orbit and point masses and in the approximation of a weak gravitational field at nonrelativistic velocities of the binary components. During the first evolutionary phase determined only by gravitational radiation, the neutron stars approach each other according to a simple analytical solution. The second evolutionary phase begins at the time of Roche-lobe filling by the low-mass component, when the second factor, mass transfer as a result of mass loss by the latter, also begins to affect the evolution. Under the simplest assumptions of conservative mass transfer and exact equality between the Roche-lobe radius and the radius of the low-mass neutron star, it is still possible to extend the analytical solution of the problem of evolution to its second phase. We present this complete solution at both phases and, in particular, give theoretical light curves of gravitational radiation that depend only on two dimensionless parameters (m _t and δ ₀). Based on the solution found, we analyze the theoretical gravitational signals from SN 1987A; this analysis includes the hypothesis about the rotational explosion mechanism for collapsing supernovae.     

Numerical integration of dynamical systems with Lie series

The integration of the equations of motion in gravitational dynamical systems—either in our Solar System or for extra-solar planetary systems—being non integrable in the global case, is usually performed by means of numerical integration. Among the different numerical techniques available for solving ordinary differential equations, the numerical integration using Lie series has shown some advantages. In its original form (Hanslmeier and Dvorak, Astron Astrophys 132, 203 1984), it was limited to the N-body problem where only gravitational interactions are taken into account. We present in this paper a generalisation of the method by deriving an expression of the Lie terms when other major forces are considered. As a matter of fact, previous studies have been done but only for objects moving under gravitational attraction. If other perturbations are added, the Lie integrator has to be re-built. In the present work we consider two cases involving position and position-velocity dependent perturbations: relativistic acceleration in the framework of General Relativity and a simplified force for the Yarkovsky effect. A general iteration procedure is applied to derive the Lie series to any order and precision. We then give an application to the integration of the equation of motions for typical Near-Earth objects and planet Mercury.     

Asteroids,comets, meteoric matter: Proceedings of the IAU Colloquium No. 22, held at the University of Nice,France, April 4–6, 1972. Edited by C. Cristescu,W. J. Klepczynski,and B. Milet. Acad. Rep. Soc. Romania,Bucharest, 1974. 333 pp., 27 L

A model is proposed for single close encounters between two small masses, m₁and m₂, which orbit a much larger mass, M. The main new feature of the model is the assumption of conic motion of the center of mass of m₁and m₂ in the gravitational field of M. Comparisons of the model with the three-body equations of motion indicate that the model is a useful approximation for m₁, m₂ ? 10^?5M. The model is therefore applicable for encounters between bodies of the order of an earth mass or smaller in the presence of the sun. Comparisons are also made of outcomes obtained by the model with outcomes of numerical integration for a large variety of close encounters. The above comparisons reveal that for many purposes the model is an adequate approximation for those encounters with ? ≥ 4, where ? is the eccentricity of the hyperbolic orbit of m₁about m₂.     

Equilibrium velocities of a planetesimal population

The random velocities v(m) of planetesimal populations specified by maximum and minimum masses and a number density n(m) ∝ m^?q are calculated interatively based on two different physical models involving ratios of rates: (1) excitation of kinetic energy by gravitational perturbation and elastic collision equal to damping of kinetic energy by inelastic collisions; and (2) excitation of kinetic energy a ratio b (～3 usually) to doubling of mass. Model (2) follows the theory of Safronov closely. Both physical models are developed approximately, using averaged factors for collision dissipation and velocity ratios, and then more precisely, allowing for reference orbit differences and for plausible variations in collisional energy dissipation with impact velocity and planetesimal mass ratios. The approximate model (2) agrees reasonably well with the results of Safronov. Both precise models are applied to populations approximating those generated by the calculations of Greenberg, Cox, and Wetherill, producing qualitatively similar velocities v(m). These results encourage analytic models of planetesimal population growth, incrementing masses and calculating velocity distributions in alternate steps. The principal improvement needed in the models is more realistic collision energy partitioning.     

Bodily tides near spin–orbit resonances

Spin–orbit coupling can be described in two approaches. The first method, known as the “MacDonald torque”, is often combined with a convenient assumption that the quality factor Q is frequency-independent. This makes the method inconsistent, because derivation of the expression for the MacDonald torque tacitly fixes the rheology of the mantle by making Q scale as the inverse tidal frequency. Spin–orbit coupling can be treated also in an approach called “the Darwin torque”. While this theory is general enough to accommodate an arbitrary frequency-dependence of Q, this advantage has not yet been fully exploited in the literature, where Q is often assumed constant or is set to scale as inverse tidal frequency, the latter assertion making the Darwin torque equivalent to a corrected version of the MacDonald torque. However neither a constant nor an inverse-frequency Q reflect the properties of realistic mantles and crusts, because the actual frequency-dependence is more complex. Hence it is necessary to enrich the theory of spin–orbit interaction with the right frequency-dependence. We accomplish this programme for the Darwin-torque-based model near resonances. We derive the frequency-dependence of the tidal torque from the first principles of solid-state mechanics, i.e., from the expression for the mantle’s compliance in the time domain. We also explain that the tidal torque includes not only the customary, secular part, but also an oscillating part. We demonstrate that the lmpq term of the Darwin–Kaula expansion for the tidal torque smoothly passes zero, when the secondary traverses the lmpq resonance (e.g., the principal tidal torque smoothly goes through nil as the secondary crosses the synchronous orbit). Thus, we prepare a foundation for modeling entrapment of a despinning primary into a resonance with its secondary. The roles of the primary and secondary may be played, e.g., by Mercury and the Sun, correspondingly, or by an icy moon and a Jovian planet. We also offer a possible explanation for the “improper” frequency-dependence of the tidal dissipation rate in the Moon, discovered by LLR.     

Families of periodic orbits in the restricted four-body problem

In this paper, families of simple symmetric and non-symmetric periodic orbits in the restricted four-body problem are presented. Three bodies of masses m ₁, m ₂ and m ₃ (primaries) lie always at the apices of an equilateral triangle, while each moves in circle about the center of mass of the system fixed at the origin of the coordinate system. A massless fourth body is moving under the Newtonian gravitational attraction of the primaries. The fourth body does not affect the motion of the three bodies. We investigate the evolution of these families and we study their linear stability in three cases, i.e. when the three primary bodies are equal, when two primaries are equal and finally when we have three unequal masses. Series, with respect to the mass m ₃, of critical periodic orbits as well as horizontal and vertical-critical periodic orbits of each family and in any case of the mass parameters are also calculated.     

Circularization of Binary Orbits

In a previous paper Mayer and Hanna (1991) discussed the process of binary orbit circularization for detached eclipsing binaries. An agreement of observational data with Tassoul's circularization time formula for binaries with masses greater than 2.5 M was found. In the present paper we also examine orbital circularization in evolved spectroscopic binaries.     

Equations of motion of the restricted problem of (2 + 2) bodies when primaries are magnetic dipoles and minor bodies are taken as electric dipoles

In this problem of the restricted (2 + 2) bodies we have considered two magnetic dipoles of masses M ₁ and M ₂(M ₁ > M ₂) moving in circular Keplarian orbit about their centre of mass. Two minor bodies of masses m ₁, m ₂(m _j< M ₂) are taken as electric dipoles in the field of rotating magnetic dipoles. These minor bodies interact with each other but do not perturb the primaries.We have found equations of motions which differ from that of Goudas and Petsagouraki's (1985).     

The restricted planar isosceles three-body problem with non-negative energy

We consider a restricted three-body problem consisting of two positive equal masses m ₁ = m ₂ moving, under the mutual gravitational attraction, in a collision orbit and a third infinitesimal mass m ₃ moving in the plane P perpendicular to the line joining m ₁ and m ₂. The plane P is assumed to pass through the center of mass of m ₁ and m ₂. Since the motion of m ₁ and m ₂ is not affected by m ₃, from the symmetry of the configuration it is clear that m ₃ remains in the plane P and the three masses are at the vertices of an isosceles triangle for all time. The restricted planar isosceles three-body problem describes the motion of m ₃ when its angular momentum is different from zero and the motion of m ₁ and m ₂ is not periodic. Our main result is the characterization of the global flow of this problem.     

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.     

The Gyldén‐type problem revisited: More refined analytical solutions

We resume and consistently extend our previous researches concerning the Gyldén‐type problem (a two‐body problem with time‐dependent equivalent gravitational parameter). To approach most of the concrete astronomical situations to be modelled in this way, we consider a periodic small perturbation. For the nonresonant case, we present a second‐order analytical solution. For the resonant case, we adopt the most realistic astronomical situation: only one dominant term of the Hamiltonian. In this case we point out a fundamental model of resonance, common to every resonant situation, and, moreover, identical to the first fundamental model of resonance . Considering the simplest model of periodic change of the equivalent gravitational parameter, we .nd that all possible resonances are con.ned to the first fundamental model. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical explorations of the rectilinear problem of three bodies

We present some results of a numerical exploration of the rectilinear problem of three bodies, with the two outer masses equal. The equations of motion are first given in relative coordinates and in regularized variables, removing both binary collision singularities in a single coordinate transformation. Among our most important results are seven periodic solutions and three symmetric triple collision solutions. Two of these periodic solutions have been continued into families, the outer massm ₃ being the family parameter. One of these families exists for all masses while the second family is a branch of the first at a second-kind critical orbit. This last family ends in a triple collision orbit.Proceedings of the Sixth Conference on Mathematical Methods in Celestial Mechanics held at Oberwolfach (West Germany) from 14 to 19 August, 1978.     

Dynamical History of the Earth–Moon System

Differential equations describing the tidal evolution of the earth's rotation and of the lunar orbital motion are presented in a simple close form. The equations differ in form for orbits fixed to the terrestrial equator and for orbits with the nodes precessing along the ecliptic due to solar perturbations. Analytical considerations show that if the contemporary lunar orbit were equatorial the evolution would develop from an unstable geosynchronous orbit of the period about 4.42 h (in the past) to a stable geosynchronous orbit of the period about 44.8 days (in the future). It is also demonstrated that at the contemporary epoch the orbital plane of the fictitious equatorial moon would be unstable in the Liapunov's sense, being asymptotically stable at early stages of the evolution. Evolution of the currently near-ecliptical lunar orbit and of the terrestrial rotation is traced backward in time by numerical integration of the evolutional equations. It is confirmed that about 1.8 billion years ago a critical phase of the evolution took place when the equatorial inclination of the moon reached small values and the moon was in a near vicinity of the earth. Before the critical epoch t _cr two types of the evolution are possible, which at present cannot be unambiguously distinguished with the help of the purely dynamical considerations. In the scenario that seems to be the most realistic from the physical point of view, the evolution also has started from a geosynchronous equatorial lunar orbit of the period 4.19 h. At t < t _cr the lunar orbit has been fixed to the precessing terrestrial equator by strong perturbations from the earth's flattening and by tidal effects; at the critical epoch the solar perturbations begin to dominate and transfer the moon to its contemporary near-ecliptical orbit which evolves now to the stable geosynchronous state. Probably this scenario is in favour of the Darwin's hypothesis about originating the moon by its separation from the earth. Too much short time scale of the evolution in this model might be enlarged if the dissipative Q factor had somewhat larger values in the past than in the present epoch. Values of the length of day and the length of month, estimated from paleontological data, are confronted with the results of the developed model.     

The geometry of the roche coordinates and zero-velocity curves in the photogravitational three-body problem

The aim of the paper is to study the geometry of the Roche curvilinear coordinates (, , ) in the photogravitational circular restricted three-body problem, with varying radiation pressure, and special attention is given to the geometry of zero-velocity curves specified by the coordinate. The radiation pressure exerted by the primary bodies on the infinitesimal third body is considered the same (q ₁ =q ₂), and the primaries are taken to have equal masses (m ₁ =m ₂). The full range of values of the common radiation factor is explored, from the valueq ₁ =q ₂ = 1 (the gravitational three-body problem) down toq ₁ =q ₂ 0. It is found that radiation has a strong influence on the geometry of the Roche coordinates and the zero-velocity curves.     

Quasinormal modes of a black hole with quintessence-like matter and a deficit solid angle: scalar and gravitational perturbations

In the previous paper (Li et al. in Phys. Lett. B 666:125–130, 2008), we show the solutions of Einstein equations with static spherically-symmetric quintessence-like matter surrounding a global monopole. Furthermore, this monopole become a black hole with quintessence-like matter and a deficit solid angle when it is swallowed by an ordinary black hole. We study its quasinormal modes by WKB method in this paper. The numerical results show that both the real part of the quasinormal frequencies and the imaginary part decrease as the state parameter w, for scalar and gravitational perturbations. And we also show variations of quasinormal frequencies of scalar and gravitational fields via different ε (deficit solid angel parameter) and different ρ ₀ (density of static spherically-symmetric quintessence-like matter at r=1), respectively.     

Closed form Vaidya-Tikekar type charged fluid spheres with pressure

Recently, Bijalwan (Astrophys. Space Sci. doi:, 2011) discussed all important solutions of charged fluid spheres with pressure and Gupta et al. (Astrophys. Space Sci. doi:, 2010) found first closed form solutions of charged Vaidya-Tikekar (V-T) type super-dense star. We extend here the approach evolved by Bijalwan (Astrophys. Space Sci. doi:, 2011) to find all possible closed form solutions of V-T type super-dense stars. The existing solutions of Vaidya-Tikekar type charged fluid spheres considering particular form of electric field intensity are being used to model massive stars. Infact at present maximum masses of the star models are found to be 8.223931M _Θ and 8.460857M _Θ subject to ultra-relativistic and non-relativistic conditions respectively. But these stars with such are large masses are not well behaved due to decreasing velocity of sound in the interior of star. We present new results concerning the existence of static, electrically charged perfect fluid spheres that have a regular interior. It is observed that electric intensity used in this article can be used to model superdense stars with ultrahigh surface density of the order 2×10¹⁴ gm/cm³ which may have maximum mass 7.26368240M _Θ for ultra-relativistic condition and velocity of sound found to be decreasing towards pressure free interface. We solve the Einstein-Maxwell equations considering a general barotropic equation of state with pressure. For brevity we don’t present a detailed analysis of the derived solutions in this paper.     

Families of periodic orbits in the planar three-body problem

A periodic orbit of the restricted circular three-body problem, selected arbitrarily, is used to generate a family of periodic motions in the general three-body problem in a rotating frame of reference, by varying the massm ₃ of the third body. This family is continued numerically up to a maximum value of the mass of the originally small body, which corresponds to a mass ratiom ₁:m ₂:m ₃?5:5:3. From that point on the family continues for decreasing massesm ₃ until this mass becomes again equal to zero. It turns out that this final orbit of the family is a periodic orbit of the elliptic restricted three body problem. These results indicate clearly that families of periodic motions of the three-body problem exist for fixed values of the three masses, since this continuation can be applied to all members of a family of periodic orbits of the restricted three-body problem. It is also indicated that the periodic orbits of the circular restricted problem can be linked with the periodic orbits of the elliptic three-body problem through periodic orbits of the general three-body problem.