The total number of giant planets in debris discs with central clearings

Infrared spectra from the Spitzer Space Telescope ( SSC ) of many debris discs are well fit with a single blackbody temperature which suggest clearings within the disc. We assume that clearings are caused by orbital instability in multiple planet systems with similar configurations to our own. These planets remove dust-generating planetesimal belts as well as dust generated by the outer disc that is scattered or drifts into the clearing. From numerical integrations, we estimate a minimum planet spacing required for orbital instability (and so planetesimal and dust removal) as a function of system age and planet mass. We estimate that a 10⁸ yr old debris disc with a dust disc edge at a radius of 50 au hosted by an A star must contain approximately five Neptune mass planets between the clearing radius and the iceline in order to remove all primordial objects within it. We infer that known debris disc systems contain at least a fifth of a Jupiter mass in massive planets. The number of planets and spacing required is insensitive to the assumed planet mass. However, an order of magnitude higher total mass in planets could reside in these systems if the planets are more massive.     

干涉技术在天体测量中的应用

介绍了射电、光学和红外波段地面干涉的仪器以及近期的进展,如厘米波段的VLBA、EVLA、SKA等;毫米波段的ALMA、CARMA等;光学波段的Keck、VLTI/PRIMA,CHARA、MarkⅢ、NPOI、LBTI等.叙述了目前干涉技术在空间天体测量中应用的一些研究项目。如射电波段的VSOP、iARISE、RadioAstro;光学波段的HST/FGS、SIM PlanetQuest等.探讨了干涉技术在天体测量中的有关研究课题,并概述了其在我国的发展和应用情况.     

干支回推-天文年代学的一种新方法

“干支回推法”是针对中国历史时期所特有的干支材料而设计的新天文年代学方法．以一组假设的材料为例,系统地给出了该方法的完整求解过程．指出该方法的优势在于准确高效,能有效地处理某些年代学难题．此外,除材料的公历日期外还能将当时的历法细节一同迭代解出．     

On some applications of invariant manifolds

Taking transfer orbits of a collinear libration point probe, a lunar probe and an interplanetary probe as examples, some applications of stable and unstable invariant manifolds of the restricted three-body problem are discussed. Research shows that transfer energy is not necessarily conserved when invariant manifolds are used. For the cases in which the transfer energy is conserved, the cost is a much longer transfer time.     

卫星跟踪卫星模式中轨道参数需求分析

首次基于半解析法利用GRACE(Gravity Recovery and Climate Experiment)双星K波段星间速度误差、GPS接收机轨道误差和加速度计非保守力误差影响累计大地水准面精度的联合模型开展了卫星跟踪卫星模式中轨道参数的需求分析.建议我国将来首颗重力卫星的平均轨道高度设计为400 km和平均星间距离设计为220 km较优.此研究不仅为我国将来卫星重力测量计划中轨道参数的优化选取以及全球重力场精度的有效和快速估计提供了理论基础和计算保证,同时对将来国际GRACE Follow-On地球重力测量计划和GRAIL(Gravity Recovery and Interior Laboratory)月球重力探测计划的发展方向具有一定的指导意义.     

Non-linear fluid dynamics of eccentric discs

A new theory of eccentric accretion discs is presented. Starting from the basic fluid-dynamical equations in three dimensions, I derive the fundamental set of one-dimensional equations that describe how the mass, angular momentum and eccentricity vector of a thin disc evolve as a result of internal stresses and external forcing. The analysis is asymptotically exact in the limit of a thin disc, and allows for slowly varying eccentricities of arbitrary magnitude. The theory is worked out in detail for a Maxwellian viscoelastic model of the turbulent stress in an accretion disc. This generalizes the conventional alpha viscosity model to account for the non-zero relaxation time of the turbulence, and is physically motivated by a consideration of the nature of magnetohydrodynamic turbulence. It is confirmed that circular discs are typically viscously unstable to eccentric perturbations, as found by Lyubarskij, Postnov & Prokhorov, if the conventional alpha viscosity model is adopted. However, the instability can usually be suppressed by introducing a sufficient relaxation time and/or bulk viscosity. It is then shown that an initially uniformly eccentric disc does not retain its eccentricity as had been suggested by previous analyses. The evolutionary equations should be useful in many applications, including understanding the origin of planetary eccentricities and testing theories of quasi-periodic oscillations in X-ray binaries.     

An extension of Newton's apsidal precession theorem

Newton's apsidal precession theorem in Proposition 45 of Book I of the 'Principia' has great mathematical, physical, astronomical and historical interest. The lunar theory and the precession of the perihelion of the planet Mercury are but two examples of the applications of this theorem. We have examined the precession of orbits under varying force laws as measured by the apsidal angle θ( N , e ), where N is the index for the centripetal force law, for varying eccentricity e . The paper derives a general function for the apsidal angle, dependent only on e and N as the potential is spherically symmetric. Further, we explore approximate ways of the solution of this equation, in the neighbourhood of   N = 2  which happens to be the case of greatest historical interest. Exact solutions are derived where they are possible. The first derivatives  ∂θ/∂ N   and  ∂θ/∂ h   [where h ( N , e ) is the angular momentum] are analytically expressed in the neighbourhood of   N = 2  (case of the inverse square law). The value of  ∂θ/∂ N   is computed numerically as well for  1 ≤ N < 3  . The resulting integrals are interesting improper integrals with singularities at both limits. Some of the integrals, especially for   N = 2  , can be given in closed form in terms of generalized hypergeometric functions which are reducible in terms of algebraic and logarithmic functions. No evidence was found for isolated cases of zero precession as e was increased. The   N = 1  case of the logarithmic potential is also briefly discussed in view of its interest for the dynamics of eccentric orbits and its relevance to realistic galaxy models. The possibility of apsidal precession was also examined for a few cases of high-eccentricity asteroids and extrasolar planets. We find that these systems may provide interesting new laboratories for studies of gravity.     

Regular and chaotic motion in softened gravitational systems

The stability of the dynamical trajectories of softened spherical gravitational systems is examined, both in the case of the full N -body problem and that of trajectories moving in the gravitational field of non-interacting background particles. In the latter case, for   N 10 000  , some trajectories, even if unstable, had exceedingly long diffusion times, which correlated with the characteristic e-folding time-scale of the instability. For trajectories of   N ≈100 000  systems this time-scale could be arbitrarily large – and thus appear to correspond to regular orbits. For centrally concentrated systems, low angular momentum trajectories were found to be systematically more unstable. This phenomenon is analogous to the well-known case of trajectories in generic centrally concentrated non-spherical smooth systems, where eccentric trajectories are found to be chaotic. The exponentiation times also correlate with the conservation of the angular momenta along the trajectories. For N up to a few hundred, the instability time-scales of N -body systems and their variation with particle number are similar to those of the most chaotic trajectories in inhomogeneous non-interacting systems. For larger N (up to a few thousand) the values of the these time-scales were found to saturate, increasing significantly more slowly with N . We attribute this to collective effects in the fully self-gravitating problem, which are apparent in the time variations of the time-dependent Liapunov exponents. The results presented here go some way towards resolving the long-standing apparent paradoxes concerning the local instability of trajectories. This now appears to be a manifestation of mechanisms driving evolution in gravitational systems and their interactions – and may thus be a useful diagnostic of such processes.     

Chaotic zone boundary for low free eccentricity particles near an eccentric planet

We consider particles with low free or proper eccentricity that are orbiting near planets on eccentric orbits. Through collisionless particle integration, we numerically find the location of the boundary of the chaotic zone in the planet's corotation region. We find that the distance in semimajor axis between the planet and boundary depends on the planet mass to the 2/7 power and is independent of the planet eccentricity, at least for planet eccentricities below 0.3. Our integrations reveal a similarity between the dynamics of particles at zero eccentricity near a planet in a circular orbit and with zero free eccentricity particles near an eccentric planet. The 2/7th law has been previously explained by estimating the semimajor at which the first-order mean motion resonances are large enough to overlap. Orbital dynamics near an eccentric planet could differ due to first-order corotation resonances that have strength proportional to the planet's eccentricity. However, we find that the corotation resonance width at low free eccentricity is small; also the first-order resonance width at zero free eccentricity is the same as that for a zero-eccentricity particle near a planet in a circular orbit. This accounts for insensitivity of the chaotic zone width to planet eccentricity. Particles at zero free eccentricity near an eccentric planet have similar dynamics to those at zero eccentricity near a planet in a circular orbit.