Round-off error in long-term orbital integrations using multistep methods

Techniques for reducing roundoff error are compared by testing them on high-order Störmer and summetric multistep methods. The best technique for most applications is to write the equation in summed, function-evaluation form and to store the coefficients as rational numbers. A larger error reduction can be achieved by writing the equation in backward-difference form and performing some of the additions in extended precision, but this entails a larger cpu cost.     

Numerical methods of star formation history measurement and applications to seven dwarf spheroidals

A comprehensive study of the measurement of star formation histories from colour–magnitude diagrams (CMDs) is presented, with an emphasis on a variety of subtle issues involved in the generation of model CMDs and maximum likelihood solution. Among these are the need for a complete sampling of the synthetic CMD, the use of proper statistics for dealing with Poisson-distributed data (and a demonstration of why χ ² must not be used), measuring full uncertainties in all reported parameters, quantifying the goodness-of-fit, and questions of binning the CMD and incorporating outside information. Several example star formation history measurements are given. Two examples involve synthetic data, in which the input and recovered parameters can be compared to locate possible flaws in the methodology (none were apparent) and measure the accuracy with which ages, metallicities and star formation rates can be recovered. Solutions of the histories of seven Galactic dwarf spheroidal companions (Carina, Draco, Leo I, Leo II, Sagittarius, Sculptor and Ursa Minor) illustrate the ability to measure star formation histories given a variety of conditions – numbers of stars, complexity of star formation history and amount of foreground contamination. Significant measurements of ancient >8 Gyr star formation are made in all seven galaxies. Sculptor, Draco and Ursa Minor appear entirely ancient, while the other systems show varying amounts of younger stars.     

New mass formula for spherically symmetric stellar configurations

A new formula is derived for the mass of spherically symmetric stellar configurations. An expression will be given where the mass square is connected to an integral over the pressure of gravitating matter. This formula turns out to hold for Newtonian gravity, for Einstein's GRT, for projective and bimetric scalar-tensor and further theories of gravitation.     

Motion in the field of a rotationally symmetric potential: Exact use of an approximate equation for the derivative of the field of directions along the normal to a trajectory

The previously derived equation (Agekyan 1974) for the derivative ?f/?n of the field of directions along the normal to a trajectory is approximate, because differentiating along the normal takes the point out of the orbit and changes the third integral of motion. However, on the envelope of the trajectory, i.e., on the contour of an orbit or a fold, ?f/?n undergoes a discontinuity of the second kind. Many authors have used this property to find points of the contours of orbits and folds. Although the integrable equation is approximate, the envelope points are determined accurately.     

Noise reduction methods in the analysis of near infrared lunar occultation light curves for high angular resolution measurements

A lunar occultation （LO） technique in the near-infrared （NIR） provides angular resolution down to milliarcseconds for an occulted source, even with ground- based 1 m class telescopes. LO observations are limited to brighter objects because they require a high signal-to-noise ratio （S/N ~40） for proper extraction of angular diameter values. Hence, methods to improve the S/N ratio by reducing noise using Fourier and wavelet transforms have been explored in this study. A sample of 54 NIR LO light curves observed with the IR camera at Mt Abu Observatory has been used. It is seen that both Fourier and wavelet methods have shown an improvement in S/N compared to the original data. However, the application of wavelet transforms causes a slight smoothing of the fringes and results in a higher value for angular diameter. Fourier transforms which reduce discrete noise frequencies do not distort the fringe. The Fourier transform method seems to be effective in improving the S/N, as well as improving the model fit, particularly in the fainter regime of our sample. These methods also provide a better model fit for brighter sources in some cases, though there may not be a significant imorovement in S/N.     

天体轨道长期数值积分的误差估计方法

数值积分方法是进行天体力学研究的重要工具, 尤其对于行星历表的研究工作而言. 由于在使用数值方法计算天体轨道时, 最终误差通常是难以预知的, 所以在面对精度要求较高或者积分时间较长的工作时具体积分方案的设计---尤其是当使用定步长方法时的步长选择---需要十分谨慎, 因为这将意味着是否能在时间成本可以被接受的范围内使解的精度达到要求. 因此, 在使用数值方法解决实际问题时如何快速寻找效率与精度之间的最佳平衡点是每一个数值积分方法的设计者与使用者都会面临的难题. 为解决这一问题, 在定步长条件下对数值积分方法的舍入误差概率分布函数以及截断误差积累量对步长的依赖关系和随时间的增长关系进行了深入研究. 基于所得结论, 提出了一种仅需较少的数值实验资料即可对选择任意时间步长积分至任意积分时刻时的舍入误差概率分布函数与截断误差积累量进行准确估计的方法, 并使用Adams-Cowell方法对该误差估计方法在圆周期轨道条件下进行了验证. 该误差估计方法在未来有望用于不同数值算法的性能对比研究, 同时也可以对数值积分方法求解实际轨道问题时的决策工作带来重要帮助.     

Degeneracy of semi-analytic solution for axially symmetric lens equations

We derive a simple semi-analytical approximation for lens equations with an arbitrary radially symmetric mass density ρ( r ), when   r /ξ₀≪ 1  and ξ₀ is the scalelength of the density profile. At the strong lensing regime, which is mostly constrained by the inner part of the mass density profile, we assume ρ∝ r ^α.
A dark matter (DM) haloes (GNFW model) are parametrized through a shape parameter α, a concentration parameter c ₁ and the total mass M . We apply our semi-analytical model to show how the solutions of the axially symmetric lens equations are degenerated in respect to the parameters α and c ₁.
In the case of an asymmetric dual image lens system, similar effective degeneracy is produced when the geometry of the lens is relaxed. Because it is impossible to determine the exact location of the source image, a family of solutions is acquired when the mass of the lens object and location of the observed images are fixed.
Our results indicate that the amount of degeneration is only weakly affected by the asymmetry in the lensing geometry set-up, e.g. the observational effective degeneracy is very close to the true physical degeneracy of the Einstein ring solutions. Basically with high-enough values for the concentration parameter, the degeneracy spawns the whole range for the shape parameter  α=[−2.0, −1.0]  .     

Numerical simulations of the Lyman α forest – a comparison of gadget-2 and enzo

We compare simulations of the Lyman α forest performed with two different hydrodynamical codes, gadget-2 and enzo . A comparison of the dark matter power spectrum for simulations run with identical initial conditions show differences of 1–3 per cent at the scales relevant for quantitative studies of the Lyman α forest. This allows a meaningful comparison of the effect of the different implementations of the hydrodynamic part of the two codes. Using the same cooling and heating algorithm in both codes, the differences in the temperature and the density probability distribution function are of the order of 10 per cent. The differences are comparable to the effects of box size and resolution on these statistics. When self-converged results for each code are taken into account, the differences in the flux power spectrum – the statistics most widely used for estimating the matter power spectrum and cosmological parameters from Lyman α forest data – are about 5 per cent. This is again comparable to the effects of box size and resolution. Numerical uncertainties due to a particular implementation of solving the hydrodynamic or gravitational equations appear therefore to contribute only moderately to the error budget in estimates of the flux power spectrum from numerical simulations. We further find that the differences in the flux power spectrum for enzo simulations run with and without adaptive mesh refinement are also of the order of 5 per cent or smaller. The latter require 10 times less CPU time making the CPU time requirement similar to that of a version of gadget-2 that is optimized for Lyman α forest simulations.     

Continuation of normal doubly symmetric orbits in conservative reversible systems

In this paper we introduce the concept of a quasi-submersive mapping between two finite-dimensional spaces, we obtain the main properties of such mappings, and we introduce “normality conditions” under which a particular class of so-called “constrained mappings” are quasi-submersive at their zeros. Our main application is concerned with the continuation properties of normal doubly symmetric orbits in time-reversible systems with one or more first integrals. As examples we study the continuation of the figure-eight and the supereight choreographies in the N-body problem.     

Numerical integration of the variational equations of satellite orbits

A recurrent power series (RPS) method is constructed for the numerical integration of the equations of motion together with the variational equations of N point masses orbiting around an oblate spheroid. By the term “variational equations” we mean the equations of the partial derivatives of the bodies’ position and velocity components with respect to the initial conditions, the relative masses and the spheroid's oblateness coefficients J₂ and J₄. The construction of recursive relations for the partial derivatives involved in the variational equations is based on partial differentiation of the corresponding recursive relations for the integration of the equations of motion. Since the number of the auxiliary variables needed for this complex system becomes tremendously large when N>1, special care must be taken during computer implementation, so as to minimize the amount of computer memory needed as well as the cost in CPU time. The RPS method constructed in this way is tested for N=1,…,4 using real initial conditions of the Saturnian satellite system. For various sets of satellites, we monitor the behaviour of all the corresponding partial derivatives. The results show a prominent difference in the behaviour of the partial derivatives between resonant and non-resonant orbital systems.     

Multiple methods for estimating the bispectrum of the cosmic microwave background with application to the MAXIMA data

We describe different methods for estimating the bispectrum of cosmic microwave background data. In particular, we construct a minimum-variance estimator for the flat-sky limit and compare results with previously studied frequentist methods. Application to the MAXIMA data set shows consistency with primordial Gaussianity. Weak quadratic non-Gaussianity is characterized by a tunable parameter   f _NL  , corresponding to non-Gaussianity at a level of  ∼10⁻⁵ f _NL  (the ratio of non-Gaussian to Gaussian terms), and we find limits of   f _NL= 1500 ± 950  for the minimum-variance estimator and   f _NL= 2700 ± 1650  for the usual frequentist estimator. These are the tightest limits on primordial non-Gaussianity, which include the full effects of the radiation transfer function.     

绝热非径向太阳振动的数值计算方法

本文采用改进后的平行Shooting方法，即连续正交归一化方法，通过增加系统的秩和把微分方程从线性变成非线性的方法，使得病态方程变成非病态方程，克服了平行Shooting方法（不连续正交归一化方法）的主要困难：一是选择径向网格点数目的问题，如果被检验振动模的径向节点数大或小，所需要的网格点数则多或少；二是在远离阈值的时候，尚未有合适的方法提高计算精度的问题．数值计算出了标准太阳模型的绝热非径向太阳振动的理论本征频率谱．从数值计算的结果发现改进后的平行Shooting方法比线性化的Henyey方法和平行Shooting方法收敛性更好且精度更高．其误差可以达到1％的精度．