Conversion from geocentric to geodetic coordinates

The geodetic latitude and the height of a satellite are obtained as power series in the ellipsoid's eccentricity, the terms being Fourier sums in the geocentric latitude with polynomials in 1/r as coefficients, with a view of determining the height to the fullest accuracy required by altimetry and geodesy from satellites.     

On the evaluation of an integral connected with the thermonuclear reaction rate in closed form

The standard expression of the reaction rate for low-energy, nonresonant nuclear reactions in nondegenerate plasma contains a parameter-dependent integral which in all previous calculations with physical or astrophysical background is considered as not capable of being evaluated in a closed form. So one usually resorts to approximation methods concerning large values of the parameter. At first we point out that CONSUL (1964) has given a series representation of the integral which was identified with a MEIJER 's G-function by MATHAI (1971). Next, in view of a physically more exact determination of the reaction rate formula, especially in connection with calculations concerning stellar energy generation, we consider a more general integral containing the mentioned one as special case and give an approximation-free representation by means of MEIJER 's G-function. The G-function so obtained may be conceived as complex-valued continuation of CONSUL 's series representation of a certain class of integrals contained in the considered one. From the series we extract a small parameter approximation of the special integral.     

The 27-Day Signal in Sunspot Number Series and the Solar Dynamo

We analyze the Wolf number daily series WN (1849 to present) as well as two other related series characterizing solar activity. Our analysis consists in computing the amplitude of a given Fourier component in a sliding time window and examining its long-term evolution. We start with the well-known 27.03- and 27.6-day periods and observe strong decadal variations of this amplitude as well as a sharp increase of the average value starting around 1905. We then consider a packet of 31 lines with periods from 25.743 to 28.453 days, which is shown to be a better representation of the synodic solar rotation. We first examine the temporal evolution of individual lines, then the energy of the packet. The energy of the packet increases sharply at the beginning of the 20th century, leading by more than two decades the well-known increase of the Wolf number. The nonaxisymmetry of sunspots increases before the total increase of activity and may be considered as a precursor. We discuss briefly and tentatively this observation in terms of solar dynamo theory.     

Note on the summation of Legendre series

A recursive procedure is established to evaluate series in themth derivatives of Legendre polynomials. It is applied to evaluate a gravitational potential, the components of its gradient and the elements of its Hessian.     

Photoelectric photometry of Z herculis

Two-colour light curves of Z Her were obtained observing the system from 1978 to 1981. The wave-like distortion at outside eclipses were derived. The amplitudes of the wave increase and its minimum migrates towards the decreasing orbital phase with a period of about 1.4 yr. An attempt was made to represent the wave-like distortions with a truncated Fourier series. The representation was found to be unsatisfactory since the distribution of the present observations is not regular but forms separate groups along the time axis. Therefore, further observations to cover the phases evenly in between the minima are needed.     

On the attitude dynamics of perturbed triaxial rigid bodies

Attitude dynamics of perturbed triaxial rigid bodies is a rather involved problem, due to the presence of elliptic functions even in the Euler equations for the free rotation of a triaxial rigid body. With the solution of the Euler–Poinsot problem, that will be taken as the unperturbed part, we expand the perturbation in Fourier series, which coefficients are rational functions of the Jacobian nome. These series converge very fast, and thus, with only few terms a good approximation is obtained. Once the expansion is performed, it is possible to apply to it a Lie-transformation. An application to a tri-axial rigid body moving in a Keplerian orbit is made.     

Compression of ephemerides

Approximations in the normL ₁ by Chebyshev polynomials are generated to represent astronomical ephemerides over large intervals of time.     

Predicted star counts and mean parallaxes down to the 23rd magnitude

Star counts and mean parallaxes as a function of B, V, R magnitudes down to 23 are presented. The data were computed by the use of two fundamental equations of stellar statistics. The assumed model considers the Galaxy as a symmetrical system with respect to its rotation axis and to its equatorial plane and as composed of the thin disk (main sequence and red giants), the thick disk and spheroid populations. Numbers of stars and mean parallaxes were derived in bins of galactic longitude and latitude of 30° and 10°, respectively. For the computation of the mean parallaxes depending on Galactic coordinates and magnitudes, series of products of Hermite and Legendre polynomials and of Fourier terms were used. The results of this paper may help in the planning of future survey missions and in the design of new telescopes. In addition, mean parallaxes can be used to derive corrections to absolute parallaxes and proper motions for any position in the sky.     

YORP torque as the function of shape harmonics

The second-order analytical approximation of the mean Yarkovsky–O'Keefe–Radzievskii–Paddack (YORP) torque components is given as an explicit function of the shape spherical harmonics coefficients for a sufficiently regular minor body. The results are based upon a new expression for the insolation function, significantly simpler than in previous works. Linearized plane-parallel model of the temperature distribution derived from the insolation function allows us to take into account a non-zero conductivity. Final expressions for the three average components of the YORP torque related with rotation period, obliquity and precession are given in a form of the Legendre series of the cosine of obliquity. The series have good numerical properties and can be easily truncated according to the degree of the Legendre polynomials or associated functions, with first two terms playing the principal role.     

Harmonic development of tide-generating potential of terrestrial planets

The aim of this study is to obtain high-accurate harmonic developments of the tide-generating potential (TGP) of Mercury, Venus and Mars. The planets’ TGP values have been first calculated on the base of DE/LE-406 numerical planetary/lunar ephemerides over a long period of time and then processed by a new spectral analysis method. According to this method the development is directly made to Poisson series where both amplitudes and arguments of the series’ terms are high-degree polynomials of time. A new harmonic development of Mars TGP is made over the time period 1900 AD–2100 AD and includes 767 second-order Poisson series’ terms of minimum amplitude equal to 10⁻⁷ m² s⁻². Analogous series composing both Mercury and Venus TGP harmonic models are built over the time period 1000 AD–3000 AD and include 1,061 and 693 terms, respectively. A modification of the standard HW95 format for representation of the terrestrial planets’ TGP is proposed. The number of terms in the planets’ TGP models transformed to the modified HW95 format is 650 for Mercury, 422 for Venus, and 480 for Mars. The quality of the new developments of the terrestrial planets’ TGP is better than that of the similar developments obtained earlier.

---

Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.      

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.     

An improved algorithm to compute circular functions of Poisson series

The efficiency in the computation of circular functions, such as cos(u) or sin(u), where u is a Poisson series, is important to derive accurate solutions of many problems of Celestial Mechanics, for instance, the orbital or rotational perturbed motion of natural or artificial bodies, since expansions in terms of Legendre functions and multiple Fourier series appear almost everywhere. Therefore, it is worth searching for alternative algorithms with lower computational cost. In this article, we propose a method based on the idea of elimination, which was originally applied to solve numerical problems, mainly in the case of matrix functions. Our comparisons with the traditional Taylor expansion prove that this new method can be more efficient when applied to compute the sine and cosine of a Poisson series.     

Updated Numerical Ephemerides of the Galilean Satellites of Jupiter

New versions of the ephemerides for the Galilean satellites of Jupiter (Io, Europa, Ganymede, and Callisto) constructed by numerically integrating the equations of motion of the satellites are presented. The satellite motionmodel takes into account the non-sphericity of Jupiter, the mutual perturbations of the satellites, and the perturbations from the Sun and major planets. The initial satellite motion parameters have been improved based on all the available series of ground-based optical observations spanning the interval 1891-2017, spacecraft observations, and radar observations. As a result, the coefficients of the expansion of the satellite coordinates and velocities in terms of Chebyshev polynomials in the interval 1891- 2025 have been obtained. The root-mean-square errors of the observations and the graphs of comparison of the constructed ephemerides both with the observations and with Lainey's numerical ephemerides are presented. The constructed ephemerides are publicly accessible.     

Electron temperature variations during solar-maximum conditions (200–3500 km)

The electron temperature variations are investigated above Arecibo, Jicamarca, Millstone Hill, St. Santin and a polar area—located at the meridian of Millstone Hill. The data analyzed represent quiet geomagnetic conditions (Kp ≤ 3) during a solar maximum (1967–1970). Between 200 and 600 km the electron temperature data stem from incoherent scatter measurements and above 600 km from the ISIS-1 observations. A simple analytical model which includes Fourier terms and cubic splines (for approximating the height dependence of the coefficients) describes the diurnal and seasonal pattern of the electron temperature in the altitude interval 200–3500 km. Three height regions are particularly striking, i.e. near 200 km where the diurnal variations show a sinusoidal pattern, the altitude interval up to approximately 1000 km which exhibits strong temperature gradients and a complex diurnal and seasonal structure, and the upper region beyond 1000 km which reflects again sinusoidal pattern but with a very pronounced latitudinal dependence.     

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

In this paper of the series, literal analytical expressions for the coefficients of the Fourier series representation ofG will be established for anyx _i; withn, N positive integers and |_i|<1 fori=1, 2, ... n. Moreover, the recurrence formulae satisfied by these coefficients will also be established. Illustrative analytical examples and a full recursive computational algorithm, with its numerical results, are included. The applications of the recurrence formulae are also illustrated by their stencils. As by-products of the analyses are two important periodic integrals developed analytically and computationally.     

Implications of Non-cylindrical Flux Ropes for Magnetic Cloud Reconstruction Techniques and the Interpretation of Double Flux Rope Events

Magnetic clouds (MCs) are a subset of interplanetary coronal mass ejections (ICMEs) which exhibit signatures consistent with a magnetic flux rope structure. Techniques for reconstructing flux rope orientation from single-point in situ observations typically assume the flux rope is locally cylindrical, e.g., minimum variance analysis (MVA) and force-free flux rope (FFFR) fitting. In this study, we outline a non-cylindrical magnetic flux rope model, in which the flux rope radius and axial curvature can both vary along the length of the axis. This model is not necessarily intended to represent the global structure of MCs, but it can be used to quantify the error in MC reconstruction resulting from the cylindrical approximation. When the local flux rope axis is approximately perpendicular to the heliocentric radial direction, which is also the effective spacecraft trajectory through a magnetic cloud, the error in using cylindrical reconstruction methods is relatively small (≈ 10^∘). However, as the local axis orientation becomes increasingly aligned with the radial direction, the spacecraft trajectory may pass close to the axis at two separate locations. This results in a magnetic field time series which deviates significantly from encounters with a force-free flux rope, and consequently the error in the axis orientation derived from cylindrical reconstructions can be as much as 90^∘. Such two-axis encounters can result in an apparent ‘double flux rope’ signature in the magnetic field time series, sometimes observed in spacecraft data. Analysing each axis encounter independently produces reasonably accurate axis orientations with MVA, but larger errors with FFFR fitting.     

Numerical Simulation of a Solar Active Region. I: Bastille Day Flare

We present three-dimensional unsteady modeling and numerical simulations of a coronal active region, carried out within the compressible single-fluid MHD approximation. We focus on AR 9077 on 14 July 2000, and the triggering of the X5.7 GOES X-ray class “Bastille Day” flare. We simulate only the lower corona, although we include a virtual photosphere and chromosphere below. The boundary conditions at the base of this layer are set using temperature maps from line intensities and line-of-sight magnetograms (SOHO/MDI). From the latter, we generate vector magnetograms using the force-free approximation; these vector magnetograms are then used to produce the boundary condition on the velocity field using a minimum energy principle (Longcope, Astrophys. J. 612, 1181, 2004). The reconnection process is modeled through a dynamical hyper-resistivity which is activated when the current exceeds a critical value (Klimas et al., J. Geophys. Res. 109, 2218, 2004). Comparing the time series of X-ray fluxes recorded by GOES with modeled time series of various mean physical variables such as current density, Poynting energy flux, or radiative loss inside the active region, we can demonstrate that the model properly captures the evolution of an active region over a day and, in particular, is able to explain the initiation of the flare at the observed time.     

Accurate analytical representation of Pluto modern ephemeris

An accurate development of the latest JPL’s numerical ephemeris of Pluto, DE421, to compact analytical series is done. Rectangular barycentric ICRF coordinates of Pluto from DE421 are approximated by compact Fourier series with a maximum error of 1.3 km over 1900–2050 (the entire time interval covered by the ephemeris). To calculate Pluto positions relative to the Sun, a development of rectangular heliocentric ICRF coordinates of the Solar System barycenter to Poisson series is additionally made. As a result, DE421 Pluto heliocentric positions by the new analytical series are represented to an accuracy of better than 5 km over 1900–2050.     

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

In this paper of the series, literal analytical expressions for the coefficients of the Fourier series representation ofF will be established for anyx _i ; withn, N positive integers 1 and | _i | fori=1, 2,...n. Moreover, the recurrence formulae satisfied by these coefficients will also be established. Illustrative analytical examples and a full recursive computational algorithm, with its numerical results, are included. The applications of the recurrence formulae are also illustrated by their stencils. As a by-product of the analyses is an integral which we may call a complete elliptic integral of thenth kind, in which the known complete elliptic integrals (1st, 2nd and 3rd kinds) are special cases of it.     

The Lambert <Emphasis Type="Italic">W</Emphasis> function and solutions to Kepler’s equation

This paper presents a method for the truncation of infinite Fourier–Bessel representations for functions requiring a solution to Kepler’s equation. Use is made of the Lambert W function to solve for the desired index that bounds the remainder terms of the series, within the prescribed tolerance. The enforcement of a maximum on the number of Bessel functions is also useful in truncating the Bessel functions themselves, resulting in an analytical representation of the solution to a desired tolerance, without the use of infinite series.