On theJ-integrals associated with close binary systems: New recursion formulae

The integralsJ _, ^m were introduced by Kopal for the numerical evaluation of the light changes exhibited by eclipsing binaries when both the tidal and rotational distrotions are taken into account.This paper is a sequel to a previous one to appear in this journal and aims at ascertaining some recursion formulae for these integrals to alleviate the computational complexity of the problem.Using a relationship existing between theJ-integrals and the Appell hypergeometric series of the first kind, we have been able to obtain recursion formulae affecting all three parametersm , of these integrals. The present stage of development has also allowed for a complete enumeration of all independent recursion formulae applicable to the case of partial eclipses.Various recursion formulae, given here for the first time, generalize previous results by Kopal which were valid form=0 ory=0.     

Evaluation of the light changes for distorted eclipsing binaries

This paper discusses a method for improving on the numerical evaluation of the light changes exhibited by a distorted eclipsing binary system.In the theory formulated by Kopal (1959), certain boundary integrals due to the distortion of both components have been calculated in terms of the Appell hypergeometric series of the first kind. The values of the four parameters appearing in these series differ according as to whether one is dealing with a partial or an annular eclipse.To accelerate the numerical evaluation of the light changes one should avoid recomputing such infinite series for contiguous values of the parameters. This can be achieved by making use of certain recursion formulae which hold for the foregoing series.We have provided here a procedure that yields forty-eight recursion formulae for the Appell hypergeometric series and have specifically calculated four new independent recursion formulae relevant to the astrophysical problem.     

Fourier analysis of the light curves of eclipsing variables,XV

A new general expression for the theoretical momentsA _2m of the light curves of eclipsing systems has been presented in the form of infinite series expansion. In this expansion, the terms have been given as the product of two different polynomials which satisfy certain three-term recursion formulae, and the coefficients diminish rapidly with increasing number of terms. Thus, the numerical values of the theoretical momentsA _2m can be generated recursively up to four significant figures for any given set of eclipse elements. This can be utilized to solve the eclipse elements in two ways: (i) with an indirect method (for the procedures see Paper XIV, Kopal and Demircan, 1978), (ii) with a direct method as minimization to the observational momentsA _2m (area fitting). The procedures given in Paper XIV for obtaining the elements of any eclipsing system consisting of spherical stars have been automated by making use of the new expression for the momentsA _2m of the light curves. The theoretical functionsf ₀,f ₂,f ₄,f ₆,g ₂ andg ₄ which are the functions ofa andc ₀, have been used to solve the eclipse elements from the observed photometric data. The closed-form expressions for the functionsf ₂,f ₄ andf ₆ have also been derived (Section 3) in terms of Kopal'sI-integrals.The automated methods for obtaining the eclipse elements from one minimum alone have been tested on the light curves of YZ (21) Cassiopeiae under the spherical model assumptions. The results of these applications will be given in Section 5 which follows a brief introduction to the procedure we followed.     

Fourier analysis of the light curves of eclipsing variables,XXII

The theoretical values of the momentsA _2m for any type of eclipses, expressed in terms of the elementsL ₁,a andc ₀, have been derived in the simple forms of rapidly convergent expansions to the series of Chebyshev polynomials, Jacobi polynomials and KopalJ-integrals (Kopal, 1977c) and hold good for any real (not necessarily integral) value ofm0.The aim of the present paper has been to establish explicit expressions for the Jacobian and its fast enough computation in the light changes of close eclipsing systems, arising from the partial derivative of different pairs ofg-functions (Kopal and Demircan, 1978, Paper XIV) with respect toa andc ₀ ² , for any type of eclipses (be these occultations or transit, partial, total or annular) and for any arbitrary degreel of the adopted law of limb-darkening. The functional behaviour of this Jacobian would determine the reasonable light curve in connection with geometrical determinacy of the parametersa andc ₀. In the expansion of Jacobian, the terms consist of two polynomials which satisfy certain three-term recursion relations having the eclipse parametersa andc ₀, as their arguments.Closed form expressions forf-functions, as well as of the Jacobian (e.g.,m=1, 2, 3), obtaining in the case of total eclipses, are given for a comparative discussion with the theoretical values of Jacobian derived from partial derivative of different pairs ofg-functions.The numerical magnitude of Jacobian would determine the best combination of the momentsA _2m in the different pairs ofg-functions and definite results would follow in the subsequent paper of this series (Edalati, 1978c, Paper XXIV).     

Fourier analysis of the light curves of eclipsing variables. XIII

New expressions for the fractional loss of light _l ⁰ have been derived in the simple forms of rapidly converging expansions to the series of Chebyshev polynomials, Jacobi polynomials, and Kopal'sJ-integrals. In these expansions, which are a supplement to those given by Kopal (1977b), variablesk andh occur in different products that simplify the numerical computation. The treatment follows the new definition of _l ⁰ which has been recently developed by Kopal (1977a).     

Determination of the elements of eclipsing variables,RW Gem and AY Cam,by Fourier analysis of their light changes

The aim of the present paper has been to present an analysis of the light curve of two eclipsing systems RW Gem and AY Cam by Fourier analysis of the light changes in the frequency domain which was developed by Kopal (1975a, b, c, d, e; 1976).In Section 1, the subject is introduced in a general way, with the intention of laying the foundation of the light curve analysis. Section 2 contains the evaluation of the empirical values of the theoretical momentA _2m is demonstrated, with the equation of the condition given. Then the equations forA _2m in terms of the elements of the total and the annular eclipses, including partial and annular phase of transit eclipse, follow.The analysis of the light curves of the two eclipsing binaries (RW Gem and AY Cam), the results and the discussion of our solution, are outlined in Section 3.     

Partial derivatives of eclipse functions αoc and αtr

Expressions are given for partial derivatives of eclipse functions with respect to geometrical depth,p, and the ratio of radii,k. The derivatives are evaluated for critical combinations ofp andk at which indeterminacies occur and the resulting expressions are listed. All expressions are given in a form suitable for numerical evaluation. Notation employed is that of Merrill.     

Coronal emission line polarization

A discussion of a program for the computation of coronal emission line polarization is presented. The starting point is a general formulation of the scattering function for magnetic dipole transitions between any two total angular momentum levels, J J, J ± 1. Illustration of the behavior of the scattering function for different transitions is given. The integration of the scattering function over the solar disk and along the line of sight accounting for arbitrary distribution of magnetic fields as well as an inhomogeneous temperature and density structure of the corona is considered next.Sample results are presented for the numerical computation of the angle of maximum polarization and the degree of maximum polarization to be expected from idealized magnetic field configurations such as radial and dipole. A computation is included for a realistic field configuration predicted to exist at the time of the 1966 eclipse. The magnetic field input to the scattering calculation is based upon the potential field extension of photospheric magnetic fields. It is the purpose of the sample calculations to demonstrate how the measurement of emission polarization measurements can be interpreted in terms of the direction of coronal magnetic fields. Factors which lend ambiguity to such interpreations are clearly illustrated from the examples. These include the Hanle-effect depolarization and the depolarization at the Van Vleck angle.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

A study of secondary cosmic ray flux variation during the annular eclipse of 15 January 2010 at Rameswaram, India

The gamma ray flux was measured during the annular solar eclipse on January 15, 2010 at Rameswaram, India using NaI (Tl) scintillator detectors and the variation in charged particles and gamma rays was monitored using Geiger Mu¨ller (GM) counters. The gamma ray flux variation was studied in energy range of 0.1–3.5 MeV. Weather parameters (temperature and humidity) were continuously recorded throughout the duration of the eclipse at the site and correlation between gamma ray flux and weather parameters was examined. Analysis of the secondary cosmic ray flux using Fast Fourier Transform (FFT) was carried out to study the impact of the eclipse on the flux modulation. An overall decrease in flux was observed by both GM counters and scintillator detectors. A relative enhancement observed for short time during the eclipse which could be associated with the presence of counter electrojet observed at Rameswaram. This is suggestive of an increase in secondary cosmic ray flux at the geomagnetic equator during every counter electrojet due to decrease in geomagnetic rigidity.     

Long-term orbit computations with KS uniformly regular canonical elements with oblateness

This paper concerns with the study of KS uniformly regular canonical elements with Earth's oblateness. These elements, ten in number, are all constant in the unperturbed motion and even in the perturbed motion, the substitution is straightforward and elementary due to the transformation laws being explicit and closed expression. By utilizing the recursion formulas of Legendre's polynomials, we are able to include any number of Earth's zonal harmonics J _n in the package and also economize the computations. A fixed step-size fourth-order Runge-Kutta-Gill method is employed for numerical integration of the canonical equations.Utilizing 5 test cases covering a large range of semimajor axis and eccentricity, we have carried out computations to study the effects of Earth's zonal harmonics (up to J ₃₆) and integration step-size variation. Bilinear relations and energy equation are used for checking the accuracies of numerical integration. From the application point of view, the package is utilized to study the behaviour of 900 km height near-circular sun-synchronous satellite orbit over a longer duration of 220 days time (nearly 3078 revolutions) and the necessity of including more number of Earth's zonal harmonic terms is noticed. The package is also used to study the effect of higher zonal harmonics on three 900 km height near-circular orbits with inclinations of 60, 63.2, and 65 degrees, by including Earth's zonal harmonics up to J ₂₄. The mean eccentricity (e _m) is found to have long-periods of 459.6, 6925.1 and 1077.6 days, respectively. Sharp changes in the variation of _m near the minima to e_m are noticed. The values of _m are found to be very near to +-90 degrees at the extrema of e_m. The same orbit is employed to study the effect of variation of inclination from 0 to 180 degrees on long-period (T) of eccentricity with J ₂ to J ₂₄ terms. T is found to increase rapidly as we proceed towards the critical inclinations.     

Motion of artificial satellites in the set of Eulerian redundant parameters,II

In this paper the full recurrent power series solution is established for J ₂-gravity perturbed motion in terms of the Eulerian redundant parameters. Applications of the method for the problem of the final state prediction are illustrated by numerical examples of some typical ballistic missiles, a final state of very high accuracy is obtained for each case study.     

Prediction of satellite motions in the Earth's gravitational field with axial symmetry by the KS regularized theory

In this paper, economical and stable recurrence formulae for the Earth's zonal potential and its gradient for the KS regularized theory will be established for any numberN of the zonal harmonic coefficient. A general recursive computational algorithm based on these formulae is also established for the initial value problem of the KS theory for the prediction of artificial satellites in the Earth's gravitational field with axial symmetry. Applications of the algorithm for the problem of the final state prediction are illustrated by numerical examples of three test orbits each for two geopotential models corresponding toN=2 andN=36. A final state of any desired accuracy is obtained for each case study, a result which shows the flexibility of the algorithm.     

Variation of Spectral Radiance of the Zenith Sky During the Annular Eclipse on May 21, 2012

The absolute brightness of the zenith sky was measured using a simple calibrated spectrometer during the annular solar eclipse event on May 21, 2012 in Fujioka City, Japan (36.2924°N, 139.0823°E). The sensitivity of the spectrometer was calibrated as a function of wavelength between 400 and 700 nm with an integral sphere. The brightness of the sky decreased to 6 % of its usual condition at the maximum magnitude of the annular eclipse of 0.95 for all wavelengths. The curve describing the variation of sky brightness accords well with the total luminosity of the solar disk estimated by a simple model that accounts for the limb darkening effect. This study provides zenith sky radiance as a function of wavelength and solar elevation angle, which is useful for the investigation of new optical instruments for atmospheric studies.     

Prediction of satellites in Earth's gravitational field with axial symmetry using Burdet's regularized theory

In this paper, economical and stable recurrence formulae for the Earth's zonal potential and its gradient for Burdet's regularized theory will be established for any number N of the zonal harmonic coefficients. A general recursive computational algorithm based on these formulae is also established for the initial value problem of Burdet oscillator for the prediction of artificial satellites in the Earth's gravitational field with axial symmetry. Applications of the algorithm for the problem of the final state prediction are illustrated by numerical examples of three test orbits each for two geopotential models corresponding to N = 2 and N = 36.A final state of any desired accuracy is obtained for each case study, a result which shows the flexibility of the algorithm.Dept. of Astronomy, KAU     

Using eclipse observations to test scintillation models

Near second and third contact during a solar eclipse the spatial spectrum of the solar illumination changes as the relative power at high spatial frequencies increases strongly. Since groundlevel atmospheric scintillation depends on a weighted integral of the image power spectrum, we can expect to see a measureable time dependence to solar scintillation during an eclipse. This effect was observed during an annular solar eclipse and quantitatively compared with a scintillation model.Operated by the Association of Universities for Research in Astronomy, Inc. (AURA) under cooperative agreement with the National Science Foundation.     

Frequency domain light curve analysis of partially eclipsing binary systems: With an application to β Persei (Algol)

Kopal's new iterative method for analysing partially eclipsing binary light curves in the frequency domain has been put into a form suitable for applications, and explicit forms for basic expressions developed. To facilitate the computations in practice, the required newS _m -functions, which are expressed in terms ofJ _, ⁰ -integrals, have been extensively tabulated by making use of hypergeometric series. The automated method has been tested successfully on the light curves of Persei (Algol). Short information on the system and the revised sets of elements including a new determination of the limb-darkening coefficient in the ultraviolet have been presented.     

Derivatives of theH-function

The first two derivatives and the mixed derivative of theH-function with respect to the angular variable and the albedo of single scattering are studied in the case of isotropic scattering using the kernel approximation method. This allows to obtain simple formulae the accuracy of which is estimated to be very good in rather low orders of approximation. Some samples of numerical results are given with eight decimal accuracy.     

Analytical solution to the perturbedJ 2 motion of Artificial Satellite in terms of Euler parameters

In this paper an analytical solution for the differential equations governing the motion of Artificial Satellite under the oblateness of the earth will be developed. We start with the differential equations in terms of Euler parameters. To compact algebra, we introduce the Cayley-Klein type complex variables. Comparison between numerical and analytical final states withJ ₂ will be given for a test case.     

Frozen orbits for satellites close to an Earth-like planet

We say that a planet is Earth-like if the coefficient of the second order zonal harmonic dominates all other coefficients in the gravity field. This paper concerns the zonal problem for satellites around an Earth-like planet, all other perturbations excluded. The potential contains all zonal coefficientsJ ₂ throughJ ₉. The model problem is averaged over the mean anomaly by a Lie transformation to the second order; we produce the resulting Hamiltonian as a Fourier series in the argument of perigee whose coefficients are algebraic functions of the eccentricity — not truncated power series. We then proceed to a global exploration of the equilibria in the averaged problem. These singularities which aerospace engineers know by the name of frozen orbits are located by solving the equilibria equations in two ways, (1) analytically in the neighborhood of either the zero eccentricity or the critical inclination, and (2) numerically by a Newton-Raphson iteration applied to an approximate position read from the color map of the phase flow. The analytical solutions we supply in full to assist space engineers in designing survey missions. We pay special attention to the manner in which additional zonal coefficients affect the evolution of bifurcations we had traced earlier in the main problem (J ₂ only). In particular, we examine the manner in which the odd zonalJ ₃ breaks the discrete symmetry inherent to the even zonal problem. In the even case, we find that Vinti's problem (J ₄+J ₂ ² =0) presents a degeneracy in the form of non-isolated equilibria; we surmise that the degeneracy is a reflection of the fact that Vinti's problem is separable. By numerical continuation we have discovered three families of frozen orbits in the full zonal problem under consideration; (1) a family of stable equilibria starting from the equatorial plane and tending to the critical inclination; (2) an unstable family arising from the bifurcation at the critical inclination; (3) a stable family also arising from that bifurcation and terminating with a polar orbit. Except in the neighborhood of the critical inclination, orbits in the stable families have very small eccentricities, and are thus well suited for survey missions.     

Complete Zonal Problem of the Artificial Satellite: Generic Compact Analytic First Order in Closed Form

This paper is a contribution to the Theory of the Artificial Satellite, within the frame of the Lie Transform as canonical perturbation technique (elimination of the short period terms). We consider the perturbation by any zonal harmonic J _n (n ≥ 2) of the primary on the satellite, what we call here the complete zonal problem of the artificial satellite. This is quite useful for primaries with symmetry of revolution. We give an analytical formula to compute directly the first order averaged Hamiltonian. The computation is carried out in closed form for all terms, avoiding therefore tedious expansions in the eccentricity or in any anomaly; this feature makes the averaging process, not only valid for all kind of elliptic trajectories but at the same time it yields the averaged Hamiltonian in a very short and compact way. The formula allows us to now skip the averaging process, which means an asymptotic gain of a factor 3n/2 regarding the computational cost of the n ^th zonal. Our analytical formulae have been widely checked, by comparison on one hand with published works (Brouwer, 1959) (which contained results for particular zonal harmonics, let’s say typically from J ₂ to J ₈), and on the other hand with the results of 3 symbolic manipulation software, among which the MM (standing for ‘Moon’s series Manipulator’), which has already been used and described in (De Saedeleer B., 2004). Additionally, the first order generator associated with this transformation is given into the same closed form, and has also been validated.