Patterns of fracture and tidal stresses on Europa

the hypothesis that lineaments on Europa are fractures produced by tidal distortion and planetary volume change is examined by comparing the orientations of dark bands, triple bands, and cuspate ridges to fracture patterns predicted for tidal distortion due to orbital recession and orbital eccentricity. If short, reticulate dark band nnear the anti-Jove point are tension cracks which formed in response to tidal distortion, they could only have been produced by orbital eccentricity. Long, arcuate dark band and triple bands peripheral to the anti-Jove point orientations which suggest that they are strike-slip faults which formed in response to orbital recession. If cuspate ridges are compressional features, their orientations and distribution suggest that they formed in response to combined orbital recession and a decrease in planetary volume. Stresses due to orbital eccentricit could have produced tension cracks near the anti-Jove point only if tensile failure occurred either prior to the accumulation of orbital recession stresses or after they had relaxed. Surface fracturing, if a consequence of tidal deformation, places important constraints on the orbital evolution of Europa.     

Periodic perturbations in close binary systems due to tidal distortion (non-lagging tides)

In a previous paper (Zafiropoulos and Kopal, 1982; hereafter referred to as Paper I) we have studies the effects of rotational distortion on the orbital elements. The aim of the present paper is to investigate the secular and periodic perturbations of the orbital elements due to tidal distortion. For tidal distortion when tides do not lag, the Gaussian form of Lagrange's planetary equations has been employed to yield the first- and second-order approximations. The results obtained include the effects produced by the second, third and fourth harmonic distortions. The first order approximation for non-lagging tides has been expressed by means of Hansen coefficients.     

The relativistic factor in the orbital dynamics of point masses

There is a growing population of relativistically relevant minor bodies in the Solar System and a growing population of massive extrasolar planets with orbits very close to the central star where relativistic effects should have some signature. Our purpose is to review how general relativity affects the orbital dynamics of the planetary systems and to define a suitable relativistic correction for Solar System orbital studies when only point masses are considered. Using relativistic formulae for the N body problem suited for a planetary system given in the literature we present a series of numerical orbital integrations designed to test the relevance of the effects due to the general theory of relativity in the case of our Solar System. Comparison between different algorithms for accounting for the relativistic corrections are performed. Relativistic effects generated by the Sun or by the central star are the most relevant ones and produce evident modifications in the secular dynamics of the inner Solar System. The Kozai mechanism, for example, is modified due to the relativistic effects on the argument of the perihelion. Relativistic effects generated by planets instead are of very low relevance but detectable in numerical simulations.     

天体运动的轨道偏心率研究概述---从恒星系统到行星系统

偏心率是描述天体运动轨道的重要参数之一, 能够为揭示天体的动力学演化提供重要线索, 进而帮助理解天体形成与演化的过程及背后的物理机制. 随着天文观测技术的不断发展, 人们对于天体运动轨道的研究已经走出太阳系, 包含的系统也从大质量端的恒星系统延伸到了低质量端的行星系统. 聚焦天体轨道偏心率研究, 回顾了目前在恒星系统(包括主序恒星、褐矮星以及致密星)和行星系统(包括太阳系外巨行星以及``超级地球''、``亚海王星''等小质量系外行星)方面取得的进展, 总结了不同尺度结构下偏心率研究的一些共同之处和待解决的问题. 并结合当下和未来的相关天文观测设备和项目, 对未来天体轨道偏心率方面的研究工作进行了展望.     

Giant planet migration through the action of disk torques and planet-planet scattering

This paper presents a parametric study of giant planet migration through the combined action of disk torques and planet-planet scattering. The torques exerted on planets during Type II migration in circumstellar disks readily decrease the semi-major axes a, whereas scattering between planets increases the orbital eccentricities ?. This paper presents a parametric exploration of the possible parameter space for this migration scenario using two (initial) planetary mass distributions and a range of values for the time scale of eccentricity damping (due to the disk). For each class of systems, many realizations of the simulations are performed in order to determine the distributions of the resulting orbital elements of the surviving planets; this paper presents the results of ∼8500 numerical experiments. Our goal is to study the physics of this particular migration mechanism and to test it against observations of extrasolar planets. The action of disk torques and planet-planet scattering results in a distribution of final orbital elements that fills the a-? plane, in rough agreement with the orbital elements of observed extrasolar planets. In addition to specifying the orbital elements, we characterize this migration mechanism by finding the percentages of ejected and accreted planets, the number of collisions, the dependence of outcomes on planetary masses, the time spent in 2:1 and 3:1 resonances, and the effects of the planetary IMF. We also determine the distribution of inclination angles of surviving planets and the distribution of ejection speeds for exiled planets.     

Radial velocity searches for other planetary systems:Current status and future prospects

Measurement of variations in the radial velocities of stars due to the reflex orbital motion of the star around the planetary-system barycenter constitutes a powerful method of searching for substellar or planetary mass companions. After several years of patient data acquisition, radial-velocity searches for planetary systems around other stars are now beginning to bear fruit. In late 1995 and early 1996, three candidate systems were announced with Jovian-mass planets around solar-type stars. The current paradigm for low-mass star formation suggests that planetary systems should be able to form in the circumstellar disks surrounding young stellar objects. These newly discovered systems, and other discoveries which will soon follow them, will test critically our understanding of the processes of star- and planet-formation. We review the techniques used in these radial-velocity searches and their results to date. We then discuss planned improvements in the surveys, and the prospects for the next 20 years.     

On the stability of circumbinary planetary systems

The dynamics of circumbinary planetary systems (the systems in which the planets orbit a central binary) with a small binary mass ratio discovered to date is considered. The domains of chaotic motion have been revealed in the “pericentric distance–eccentricity” plane of initial conditions for the planetary orbits through numerical experiments. Based on an analytical criterion for the chaoticity of planetary orbits in binary star systems, we have constructed theoretical curves that describe the global boundary of the chaotic zone around the central binary for each of the systems. In addition, based on Mardling’s theory describing the separate resonance “teeth” (corresponding to integer resonances between the orbital periods of a planet and the binary), we have constructed the local boundaries of chaos. Both theoretical models are shown to describe adequately the boundaries of chaos on the numerically constructed stability diagrams, suggesting that these theories are efficient in providing analytical criteria for the chaoticity of planetary orbits.     

Possible patterns in the distribution of planetary formation regions

Eris, an object larger than Pluto, is known to reside in the transneptunian region further away than Pluto. One can wonder whether its semimajor orbital axis fits in a generalized Titius–Bode law, in the same way as Pluto does. We performed a new least-squares fit to a generalized Titius–Bode law including Eris and found that not only does Eris fit in the trend, but also the correlation coefficient improves. In addition, there is a remarkable symmetry of the location of the planetary formation regions with respect to Jupiter when the natural logarithm of the heliocentric distance is used as the metric. The issue of whether the observed patterns have some physical meaning or are due to mere chance is addressed using a Monte Carlo approach identical to that of Lynch. Although the probability of chance occurrence is highly dependent on the way in which the random configurations of synthetic planetary systems are selected, we find that in all reasonable scenarios of random planetary systems the probability of chance occurrence of the observed patterns is small (below 1 per cent in most cases). If the trend were used as a prediction tool, one might expect another planet or dwarf planet or a swarm of bodies with semimajor orbital axis of 120 ± 20 au. Simple calculations show that the protoplanetary nebula most likely had enough mass to allow the accretion of at least a dwarf planet at that distance. We also found that if the surface density of the nebula decayed with heliocentric distance ( r ) as a power of −2, the regular spacing in ln  r in the Solar system could be a natural consequence of the existence of a threshold mass for planetary formation.     

Orbital evolution of a planet with tidal dissipation in a restricted three-body system

The angle between planetary spin and the normal direction of an orbital plane is supposed to reveal a range of information about the associated planetary formation and evolution. Since the orbit's eccentricity and inclination oscillate periodically in a hierarchical triple body and tidal friction makes the spin parallel to the normal orientation of the orbital plane with a short timescale in an isolated binary system, we focus on the comprehensive effect of third body perturbation and tidal mechanism on the angle. Firstly, we extend the Hut tidal model(1981) to the general spatial case, adopting the equilibrium tide and weak friction hypothesis with constant delay time, which is suitable for arbitrary eccentricity and any angle ? between the planetary spin and normal orientation of the orbital plane. Furthermore, under the constraint of angular momentum conservation, the equations of orbital and ratational motion are given. Secondly, considering the coupled effects of tidal dissipation and third body perturbation, and adopting the quadrupole approximation as the third body perturbation effect, a comprehensive model is established by this work. Finally, we find that the ultimate evolution depends on the timescales of the third body and tidal friction. When the timescale of the third body is much shorter than that of tidal friction, the angle ? will oscillate for a long time,even over the whole evolution; when the timescale of the third body is observably larger than that of the tidal friction, the system may enter stable states, with the angle ? decaying to zero ultimately, and some cases may have a stable inclination beyond the critical value of Lidov-Kozai resonance. In addition, these dynamical evolutions depend on the initial values of the orbital elements and may aid in understanding the characteristics of the orbits of exoplanets.     

Numerical Instruments for the Analysis of Secular Dynamics of Exoplanetary Systems

A review is given of modern numerical methods for the analysis of resonant and chaotic dynamics: calculation of the Lyapunov characteristic exponents, the MEGNO method, and the maximum eccentricity method. These methods are used to construct stability diagrams for the planetary systems γ Cep, HD 196885, and HD 41004. The diagrams are analyzed to determine the most probable values taken by the orbital parameters of the exoplanets and obtain estimates for the Lyapunov time of their orbital dynamics. The stability diagrams constructed using the different methods are compared to analyze their effectiveness in the study of secular dynamics of exoplanetary systems.     

Metallicity of Sun-like G-stars that have Exoplanets

By considering the physical and orbital characteristics of G type stars and their exoplanets, we examine the association between stellar mass and its metallicity that follows a power law. Similar relationship is also obtained in case of single and multiplanetary stellar systems suggesting that, \(\hbox {Sun}^{\prime }\)s present mass is about 1% higher than the estimated value for its metallicity. Further, for all the stellar systems with exoplanets, association between the planetary mass and the stellar metallicity is investigated, that suggests planetary mass is independent of stellar metallicity. Interestingly, in case of multiplanetary systems, planetary mass is linearly dependent on the stellar absolute metallicity, that suggests, metal rich stars produce massive (\(\ge \)1 Jupiter mass) planets compared to metal poor stars. This study also suggests that there is a solar system planetary missing mass of \({\sim }\)0.8 Jupiter mass. It is argued that probably 80% of missing mass is accreted onto the Sun and about 20% of missing mass might have been blown off to the outer solar system (beyond the present Kuiper belt) during early history of solar system formation. We find that, in case of single planetary systems, planetary mass is independent of stellar metallicity with an implication of their non-origin in the host star’s protoplanetary disk and probably are captured from the space. Final investigation of dependency of the orbital distances of planets on the host stars metallicity reveals that inward migration of planets is dominant in case of single planetary systems supporting the result that most of the planets in single planetary systems are captured from the space.     

Global dynamics of the Gliese 876 planetary system

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

Determination of the properties of Mercury's magnetic field by the MESSENGER mission

The MESSENGER mission to Mercury, to be launched in 2004, will provide an opportunity to characterize Mercury's internal magnetic field during an orbital phase lasting one Earth year. To test the ability to determine the planetary dipole and higher-order moments from measurements by the spacecraft's fluxgate magnetometer, we simulate the observations along the spacecraft trajectory and recover the internal field characteristics from the simulated observations. The magnetic field inside Mercury's magnetosphere is assumed to consist of an intrinsic multipole component and an external contribution due to magnetospheric current systems described by a modified Tsyganenko 96 model. Under the axis-centered-dipole approximation without correction for the external field the moment strength is overestimated by ∼4% for a simulated dipole moment of , and the error depends strongly on the magnitude of the simulated moment, rising as the moment decreases. Correcting for the external field contributions can reduce the error in the dipole term to a lower limit of ∼1-2% without a solar wind monitor. Dipole and quadrupole terms, although highly correlated, are then distinguishable at the level equivalent to an error in the position of an offset dipole of a few tens of kilometers. Knowledge of the external magnetic field is therefore the primary limiting factor in extracting reliable knowledge of the structure of Mercury's magnetic field from the MESSENGER observations.     

Study of dynamical processes at the final stages of planetary system formation: Resonance motion of giant planets

According to current observational data, planets of many exoplanetary systems have resonant motion. The formation of resonance configurations is studied within a unified model of planetary migration. Planets in the observed systems 24 Sex, HD 37124, HD 73526, HD 82943, HD 128311, HD 160691, Kepler 9, NN Ser, which are moving in the 2: 1 resonance, could have been captured into this resonance due to both the Type I and II migration with a wide range of parameters. The migration conditions are defined for the formation of HD 45364 and HD 200964 that are in the 3: 2 and 4: 3 first-order resonances, correspondingly. The results obtained for HD 200964 show that planets can be captured in the first-order resonances, when the outer-to-inner orbital period ratios for the planets are less than 3: 2, only if Type I migration rates are large, and the mass of at least one planet is substantially less than the modern masses of the observed giant planets. The formation of the HD 102272, HD 108874, HD 181433 and HD 202206 systems with planets in high-order resonances is considered. The capture into these resonances can be realized with very slow Type II migration. Possible bounds for migration parameters are considered. In particular, it has been found that the capture of HD 108874 into the 4: 1 resonance is possible only if the angle between the plane of planetary orbits and the plane of sky is appreciably less than 90°, i.e., the planetary masses are a few times larger than the minimum values. The capture of HD 202206 into the 5: 1 resonance is possible at low migration rates; however, another mechanism is required to explain the high observed eccentricity of the inner planet (for example, strong gravitational interaction between the planets). Resonant configurations can be disrupted due to the interaction between planets and remaining fragments of the planetesimal disk as, for example, may occur in the three-planet system 47 UMa. The specific orbital features observed for this system are explained.     

Computer simulations of planetary accretion dynamics: Sensitivity to initial conditions

We have tested the implications and limitations of Program ACRETE, a scheme based merely on Newtonian physics and accretion with unit sticking efficiency, devised by Dole in 1970 to simulate the origin of the planets. The dependence of the results on a variety of radial and vertical density distribution laws, on the ratio of gas to dust in the solar nebula, on the total nebula mass, and on the orbital eccentricity, ?, of the accreting grains are explored. Only for a small subset of conceivable cases are planetary systems closely like our own generated. Many models have tendencies toward one of two preferred configurations: multiple-star systems, or planetary systems in which Jovian planets either have substantially smaller masses than in our system or are absent altogether. But for a wide range of cases recognizable planetary systems are generated, ranging from multiple-star systems with accompanying planets, to systems with Jovian planets at several hundred astronomical units, to single stars surrounded only by asteroids. Many systems exhibit planets like Pluto and objects of asteroidal mass, in addition to usual terrestrial and Jovian planets. No terrestrial planets were generated more massive than five Earth masses. The number of planets per system is for most cases of order 10, and, roughly, inversely proportional to ?. All systems generated obey a relation of the Titius-Bode variety for relative planetary spacing. The case with which planetary systems are generated using such elementary and incomplete physical assumptions supports the idea of abundant and morphologically diverse planetary systems throughout the Galaxy.     

Effect of Orbital Distance on the Atmospheric Escape of Exoplanets

Atmospheric escape is an important sector in the evolution of planetary atmosphere, and its energy is mainly originated from the radiation of the host star at the high energy band. The radiation flux drops dramatically with the increase of orbital distance, there is a large difference of planetary atmospheric escape in different orbits, so it is necessary to study the impact of orbital distance on the atmospheric escape of an exoplanet. We consider the radiation transfer and the photochemical reactions of multiple kinds of particles to study the variation of planetary atmospheric escape with the orbital distance by using a 1-D hydrodynamic model. Due to the large differences of the spectra of host stars in different evolution stages, the Astrophysical Plasma Emission Code (APEC) in the X-Ray Spectral Fitting Package (XSPEC) is used to obtain the spectra of solar-type stars with different ages as the input spectra of the model. The results indicate that the escape rates of the exoplanets in different orbits are different significantly, and the escape mechanism is converted from the drastic hydrodynamic escape into the moderate Jeans escape as the orbital distance increases, the smaller the planetary gravitational potential, the younger the star-planet system, the larger the distance of this conversion. The correlation between the escape rate and the radiation flux decreases for the short-period exoplanets in a younger star-planet system. It is shown that the classical energy-limited escape theory is not suitable for this kind of exoplanets. These results have enriched the studies on the atmospheric escape of exoplanets, especially, extended the studies on the escape mechanism and energy conversion under different orbital distances and stellar ages.     

Fresip: A mission to determine the character and frequency of extra-solar planets around solar-like stars

FRESIP (FRequency of Earth-Sized Inner Planets) is a mission designed to detect and characterize Earth-sizes planets around solar-like stars. The sizes of the planets are determined from the decrease in light from a star that occurs during planetary transits, while the orbital period is determined from the repeatability of the transits. Measurements of these parameters can be compared to theories that predict the spacing of planets, their distribution of size with orbital distance, and the variation of these quantities with stellar type and multiplicity. Because thousands of stars must be continually monitored to detect the transits, much information on the stars can be obtained on their rotation rates and activity cycles. Observations of p-mode oscillations also provide information on their age and composition. These goals are accomplished by continuously and simultaneously monitoring 500 solar-like stars for evidence of brightness changes caused by Earth-sized or larger planetary transits. To obtain the high precision needed to find planets as small as the Earth and Venus around solar-like stars, a wide field of view Schmidt telescope with an array of CCD detectors at its focal plane must be located outside of the Earth's at mosphere. SMM (Solar Maximum Mission) observations of the low-level variability of the Sun (1:100,000) on the time scales of a transit (4 to 16 hours), and our laboratory measurements of the photometric precision of charge-coupled devices (1:100,000) show that the detection of planets as small as the Earth is practical. The probability for detecting transits is quite favorable for planets in inner orbits. If other planetary systems are similar to our own, then approximately 1% of those systems will show transits resulting in the discovery of 50 planetary systems in or near the habitable zone of solar-like stars.     

Canonical planetary perturbation equations for velocity-dependent forces,and the Lense-Thirring precession

A form of planetary perturbation theory based on canonical equations of motion, rather than on the use of osculating orbital elements, is developed and applied to several problems of interest. It is proved that, with appropriately selected initial conditions on the orbital elements, the two forms of perturbation theory give rise to identical predictions for the observable coordinates and velocities, while the orbital elements themselves may be strikingly different. Differences between the canonical form of perturbation theory and the classical Lagrange planetary perturbation equations are discussed. The canonical form of perturbation theory in some cases has advantages when the perturbing forces are velocity-dependent, but the two forms of perturbation theory are equivalent if the perturbing forces are dependent only on position and not on velocity. The canonical form of the planetary perturbation equations are derived and applied to the Lense Thirring precession of a test body in a Keplerian orbit around a rotating mass source.     

Analytical method for the effects of the asteroid belt on planetary orbits

Analytic expressions are derived for the perturbation of planetary orbits due to a thick constant density asteroid belt. The derivations include extensions and adaptations of Plakhov's analytic expressions for the perturbations in five of the orbital elements for closed orbits around Saturn's rings. The equations of Plakhov are modified to include the effect of ring thickness and additional equations are derived for the perturbations in the sixth orbital element, the mean anomaly. The gravitational potential and orbital perturbations are derived for the asteroid belt with and without thickness, and for a hoop approximation to the belt. The procedures are also applicable to Saturn's rings and the newly discovered rings of Uranus.The effects of the asteroid belt thickness on the gravitational potential coefficients and the orbital motions are demonstrated. Comparisons between the Mars orbital perturbations obtained using the analytic expressions and those obtained using numerical integration are discussed. The effects of the asteroid belt on the Earth based ranging to Mars are also demonstrated.