A unified theory of polynomial expansions and their applications involving Clebsch-Gordan type linearization relations and Neumann series

Motivated by their potential for applications in several diverse fields of physical, astrophysical, and engineering sciences, this paper aims at presenting a unified study of various classes of polynomial expansions and multiplication theorems associated with the general multivariable hypergeometric function (studied recently by A. W. Niukkanen and H. M. Srivastava), which provides an interesting and useful unifiation of numerous families of special functions in one and more variables, encoutered naturally (and rather frequently) in many physical, quantum chemical, and quantum mechanical situations. Several interesting applications of these general polynomial expansions are considered, not only in the derivations of various Clebsch-Gordan type linearization relations involving products of several Jacobi or Laguerre polynomials, but also to associated Neumann expansions in series of the Bessel functionsJ _v (z) andI _v (z) (and of their suitable products).     

Fourier techniques for an analysis of eclipsing binary light curves: I

The Fourier techniques developed so far for an analysis of eclipsing binary light curves have been re-discussed. The Fourier coefficients for the analysis have been derived in a simple form of series expansions, in terms of eclipse elements, valid for any type of eclipse (regardless of whetherr ₁r ₂).These coefficients may be utilized to solve the eclipse elements in terms of the observed characteristics of the light curves. A general relation between the observed quantitiesl and , and the eclipse elementsr _1,2,i andL ₁ has also been given in the form of series expansions which can be used for the synthesis of the light curves.     

Expansion theory for the elliptic motion of arbitrary eccentricity and semi-major axis

In this paper of the series, elliptic expansions in terms of the sectorial variables _j ⁽ⁱ⁾ introduced recently in Paper IV (Sharaf, 1982) to regularize highly oscillating perturbations force of some orbital systems will be established analytically and computationally for the fifth and sixth categories. For each of the elliptic expansions belonging to a category, literal analytical expressions for the coefficients of its trigonometric series representation are established. Moreover, some recurrence formulae satisfied by these coefficients are also established to facilitate their computations; numerical results are included to provide test examples for constructing computational algorithms.     

Expansion theory for the elliptic motion of arbitrary eccentricity and semi-major axis

In this paper of the series, elliptic expansions in terms of the sectorial variables _j ⁽ⁱ⁾ introduced by the author in Paper IV (Sharaf, 1982) to regularize the highly oscillating perturbation force of some orbital systems will be established analytically and computationally for the ninth, tenth, eleventh, and twelfth categories according to our adopted scheme of presentation drawn up in Paper V (Sharaf, 1983). For each of the elliptic expansions belonging to a category, literal analytical expressions for the coefficients of its trigonometric series representation are established. Moreover, some recurrence formulae satisfied by these coefficients are also established to facilitate their computation, and numerical results are included to provide test examples for constructing computational algorithms. Finally, the first collection of completed elliptic expansions in terms of _j ⁽ⁱ⁾ so explored will be given in Appendix A for the guidance of the reader.     

Expansion theory for elliptic motion of arbitrary eccentricity and semi-major axis

In this paper of the series, we arrive at the end of the second step of our regularization approach, and in which, elliptic expansions in terms of the sectorial variables _j ⁽ⁱ⁾ introduced by the author in Paper IV (Sharaf, 1982b) to regularize the highly oscillating perturbation force of some orbital systems will be established analytically and computationally for the thirteenth, fourteenth, fifteenth, and sixteenth categories according to our adopted scheme of presentation drawn up in Paper V (Sharaf, 1983). For each of the elliptic expansions belonging to a category, literal analytic expressions for the coefficients of its trigonometric series representation are established. Moreover, some recurrence formulae satisfied by these coefficients are also established to facilitate their computations, and numerical results are included to provide test examples for constructing computational algorithms. Finally, the second and the last collection of completed elliptic expansion will be given in Appendix B, such that, the materials of Appendix A of Paper VIII (Sharaf, 1985b) and those of Appendix B of the present paper provide the reader with the elliptic expansions in terms of _j ⁽ⁱ⁾ so explored for the second step of our regularization approach.     

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).     

Recurrence formulae for the Hansen's developments

In this paper we derive some recurrence formulae which can be used to calculate the Fourier expansions of the functions (r/a) ⁿ cosmv and (r/a) ⁿ sinmv in terms of the eccentric anomalyE or the mean anomalyM. We also establish a recurrence process for computing the series expansions for alln andm when the expansions of two basic series are known. These basic series were given in explicit form in the classical literature. The recurrence formulae are linear in the functions involved and thus make very simple the computation of the series.This work was supported by NASA contract No. NASr 54(06).—The paper was presented at the AIAA/AAS meeting, Princeton University, August 1969.     

Fourier techniques for an analysis of eclipsing binary light curves. VIa

The aim of the present paper is to deduce relations between the integral transformsA _2m, B_2m,andF _1,2 of the light curves of eclipsing binary systems. The integral transformsA _2m, B_2m,andF _1,2 have been related to one another by means of finite or rapidly converging infinite summations obtained by integrations of the series expansions of trigonometric functions.     

Analytical approach using KS elements to short-term orbit predictions including J 2

Analytical solutions using KS elements are derived. The perturbation considered is the Earth's zonal harmonic J ₂. The series expansions include terms of fourth power in the eccentricity. Only two of the nine KS element equations are integrated analytically due to the reasons of symmetry. The analytical solution is suitable for short-term orbit computations. Numerical studies show that reasonably good estimates of the orbital elements can be obtained in one step of 10 to 30 degrees of eccentric anomaly for near-Earth orbits of moderate eccentricity. For application purposes, the analytical solution can be effectively used for onboard computation in the navigation and guidance packages, where the modelling of J ₂ effect becomes necessary.     

Expansion theory for the elliptic motion of arbitrary eccentricity and semi-major axis

In this paper of the series, the expansions of the functionsH ₁,H ₂, andH ₃ will be established analytically and computationally form positive integer,q any real number and , are both positive <1. Full recursive computational algorithms with their numerical results will also be included.     

Analytical short-term orbit predictions withJ 3 andJ 4 in terms of KS elements

Analytical theory for short-term orbit motion of satellite orbits with Earth's zonal harmonicsJ ₃ andJ ₄ is developed in terms of KS elements. Due to symmetry in KS element equations, only two of the nine equations are integrated analytically. The series expansions include terms of third power in the eccentricity. Numerical studies with two test cases reveal that orbital elements obtained from the analytical expressions match quite well with numerically integrated values during a revolution. Typically for an orbit with perigee height, eccentricity and inclination of 421.9 km, 0.17524 and 30 degrees, respectively, maximum differences of 27 and 25 cm in semimajor axis computation are noted withJ ₃ andJ ₄ term during a revolution. For application purposes, the analytical solutions can be used for accurate onboard computation of state vector in navigation and guidance packages.     

Expansions of (r/a)m cos jv and (r/a)m sin jv to high eccentricities

Expansions of the functions (r/a)^cos jv and (r/a)m ^sin jv of the elliptic motion are extended to highly eccentric orbits, 0.6627 ... <e<1. The new expansions are developed in powers of (e–e*), wheree* is a fixed value of the eccentricity. The coefficients of these expansions are expressed in terms of the derivatives of Hansen's coefficients with respect to the eccentricity. The new expansions are convergent for values of the eccentricity such that |e–e*|<(e*), where the radius of convergence (e*) is the same of the extended solution of Kepler's equation. The new expansions are intrinsically related to Lagrange's series.     

Expansion theory for the elliptic motion of arbitrary eccentricity and semi-major axis

In this paper of the series, elliptic expansions in terms of the sectorial variables θ _j ⁽ⁱ⁾ introduced in Paper IV (Sharaf, 1982) to regularise highly oscillating perturbation force of some orbital systems will be explored for the first four categories. For each of the elliptic expansions belonging to a category, literal analytical expressions for the coefficients of its trigonometric series representation are established. Moreover, some recurrence formulae satisfied by these coefficients are also established to facilitate their computations, numerical results are included to provide test examples for constructing computational algorithms.     

Fourier analysis of the light curves of eclipsing variables,X

The aim of the present paper will be to develop methods for computation of the Fourier transforms of the light curves of eclipsing variables — due to any type of eclipses — as a function of a continuous frequency variablev. For light curves which are symmetrical with respect to the conjunctions (but only then) these transforms prove to be real functions ofv, and expressible as rapidly convergent expansions in terms of the momentsA _2m+1 of the light curves of odd orders. The transforms are found to be strongly peaked in the low-frequency domain (attaining a maximum forv=0), and become numerically insignificant forv>3. This is even more true of their power spectra.The odd momentsA _2m+1 — not encountered so far in our previous papers — are shown in Section 3 of the present communication to be expressible as infinite series in terms of the even momentsA _2m well known to us from Papers I–IV; and polynomial expressions are developed for approximating them to any desired degree of accuracy. The numerical efficiency of such expressions will be tested in Section 4, by application to a practical case, with satisfactory results.Lastly, in Section 5, an appeal to the Wiener-Khinchin theorem (relating the power spectra with autocorrelation function of the light curves) and Parseval's theorem on Fourier series will enable us to extend our previous methods for a specification of quadratic moments of the light curves in terms of the linear ones.     

Sur de nouvelles séries pour le problème de masses critiques de Routh dans le problème restreint plan des trois corps

The planar restricted 3-body problem, linearized in the neighborhood of Lagrangian equilibriaL ₄ andL ₅, has in general two distinct eigenvalues and their opposites. When they are pure imaginary and not multiples of each other, they generate two families of periodic solutions called long and short periodic families. This is essentially a consequence of the famous theorem of Liapunov (Siegel, 1956). We showed (Roels, 1971b) how to solve the problem when the eigenvalues are multiples of each other in building series with negative exponents instead of the integer expansions of Siegel (Roels and Lauterman, 1970). When the eigenvalues are equal, which is the case for the mass ratio of Routh, the problem was solved by Deprit and Henrard (1968) using formal series in ordinary unnormalized variables. That leads to very complicated series because of the use of variables that are not well adapted to the problem. The convergence of the series was proven by Meyer and Schmidt (1971). In this paper we solve the problem by using normalized variables. This brings us to build expansions with fractional exponents. So in summary, normalized variables generate integer series in the non-resonant cases, series with negative exponents in the case of resonancek≥3, and series with fractional exponents when the resonance is 1.     

The representation of the force functions of the Earth's and Moon's gravitation in ellipsoidal coordinates

This paper derives asymptotic expansions of ellipsoidal coordinates in Cartesian coordinates and an expansion in spherical harmonics of the dominant term for the solution of Laplace's equation corresponding to the gravitational force function for a two-dimensional finite body.On comparing the expansion of the dominant term derived here with known expansions of the force functions of the Earth's and Moon's gravitation the author obtains values for the semimajor axes and eccentricities of the singular ellipses of these bodies in terms of the second degree harmonic coefficientsc ₂₀ andc ₂₂.     

Astrophysical spectroscopy and neutron reactions: Integral transforms and Voigt functions

The Voigt functionsK(x, y) andL(x, y) which play an essential role in astrophysical spectroscopy and neutron physics are investigated and generalized from the viewpoint of integral operators. Unified representations and series expansions involving classical functions of mathematical physics and multivariable hypergeometric functions are established. From the delicate asymptotic analysis of Laplace and Hankel integral transforms we extract complete and rigorous asymptotic expansions of the generalized Voigt functions for large values of the variablesx andy which are of great value in the theory of spectral line profiles.     

Fourier analysis of the light curves of eclipsing variables

The aim of the present paper has been to construct analytic expressions for incomplete Fourier transforms underlying Kitamura's method for an analysis of the light curves of eclipsing binary systems (Kitamura, 1965).The expansions established for Kitamura's coefficientsc _n ands _n have been used as a basis for checking numerical accuracy of these coefficients tabulated by Kitamura in 1967; and the outcome bespeaks a high quality of his tables constructed by numerical quadratures.     

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).     

Expansions of elliptic motion based on elliptic function theory

New expansions of elliptic motion based on considering the eccentricitye as the modulusk of elliptic functions and introducing the new anomalyw (a sort of elliptic anomaly) defined byw=u/2K–/2,g=amu–/2 (g being the eccentric anomaly) are compared with the classic (e, M), (e, v) and (e, g) expansions in multiples of mean, true and eccentric anomalies, respectively. These (q,w) expansions turn out to be in general more compact than the classical ones. The coefficients of the (e,v) and (e,g) expansions are expressed as the hypergeometric series, which may be reduced to the hypergeometric polynomials. The coefficients of the (q,w) expansions may be presented in closed (rational function) form with respect toq, k, k=(1–k ²)^1/2,K andE, q being the Jacobi nome relatedk whileK andE are the complete elliptic integrals of the first and second kind respectively. Recurrence relations to compute these coefficients have been derived.on leave from Institute of Applied Astronomy, St.-Petersburg 197042, Russia