Precision of the ephemerides of outer planetary satellites

Ephemerides of planetary satellites are needed to address many problems. These ephemerides are used for subsequent observations. A comparison of the available ephemerides with new observations allows the accuracy of the former to be assessed. However, the precision of the ephemerides must be known a priori when solving the tasks. In this paper we formulate and solve the problem of estimating the precision of the ephemerides of outer planetary satellites derived from observations when applied up to the future moments.The methods of assessing the precision of ephemerides involve producing a set of samples of the same ephemeris inferred from observations with different samples of Monte Carlo generated random errors (RO) superimposed onto it. The statistical parameters of simulated observational errors are based on the results of the reduction of real satellite observations. We compute the deviations of the samples of the ephemeris from the standard ephemeris inferred from real observations and adopt the root-mean-square deviation of the apparent coordinates as the precision of the ephemeris. We also use alternative methods: one based on the matrix of covariances of parameter errors (RP), and another one based on bootstrap samples of observations (BS).We use three methods (RO, RP, and BS) to estimate the precision of the ephemerides of all the 107 outer planetary satellites over the 2010-2020 time interval. The precision of the ephemerides of different satellites varies from 0.05 to 4.0 arcsec. For a number of satellites new observations are of vital importance for maintaining the precision of the ephemerides at a level that would allow identification of satellites during the reduction of observations. For some satellites the precision of their ephemerides is of the order of the sizes of their orbits and such satellites can be considered to have been lost. We show that the method of bootstrap samples (BS) can give doubtful results in the cases where there are few observations, which covered a time interval that is shorter than the orbital period of the satellite.Our results suggest obtaining more precise ephemeris making new observations at the times of maximum estimated errors of the ephemeris.All the inferred estimates of the precision of ephemerides are available from the MULTI-SAT ephemeris server: www.imcce.fr/sat (IMCCE), www.sai.msu.ru/neb/nss/index.htm (SAI).     

A theory of integration for the extended Delaunay method

In this paper a method for the integration of the equations of the extended Delaunay method is proposed. It is based on the equations of the characteristic curves associated with the partial differential equation of Delaunay-Poincaré. The use of the method of characteristics changes the partial differential equation for higher order approximations into a system of ordinary differential equations. The independent variable of the equations of the characteristics is used instead of the angular variables of the Jacobian methods and the averaging principle of Hori is applied to solve the equations for higher orders. It is well known that Jacobian methods applied to resonant problems generally lead to the singularity of Poincaré. In the ideal resonance problem, this singularity appears when higher order approximations of the librational motion are considered. The singularity of Poincaré is non-essential and is caused by the choice of the critical arguments as integration variables. The use of the independent variable of the equation of the characteristics in the place of the critical angles eliminates the singularity of Poincaré.     

The solution of Jacobi's virial equation for celestial bodies

A relationship between the potential energy and the moment of the inertia for celestial bodies is heuristically discovered. This relationship consists in the constancy of the product of formfactors for the potential energy and the moment of the inertia. The product is independent of the body mass and its radial mass distribution.We find the exact solution of Jacobi's virial equation for a gravitating spherical body based on the relationship obtained. This solution represents the unharmonic radial oscillations of the body. The solution is valid for a wide class of celestial bodies including variable stars and relativistic objects for which a relativistic analog of Jacob's equation is derived.The period of the radial oscillations of the planets is estimated with the help of the solution found. We note the coincidence of the experimental data and our theoretical calculations for the Sun.We show the important role of the Coulomb forces in the formation of the planets. It is demonstrated that the Coulomb forces result in the relation between the planet masses and their average molecular weight.     

The 1997 Campaign of Observations of Mutual Occultations and Eclipses in the System of Jupiter's Galilean Satellites: Final Results

This paper presents the final results of the campaign of photometric observations of the Galilean satellites of Jupiter at the time of their mutual occultations and eclipses at the epoch of 1997 at observatories of Kazakhstan, Russia, and Ukraine. These results essentially contribute to the worldwide database of this type. They will help to refine the models of motion of the Galilean satellites of Jupiter. A comparison is made of the obtained data with the results of the worldwide campaign of observations of mutual occultations and eclipses in 1997. From the results of these observations, the differences of planetocentric coordinates of two satellites are calculated for a sequence of instants of time. Processing is based on the theory and the model of mutual occultations and eclipses of natural planetary satellites, which were developed by one of the authors in previous papers (Emel'yanov, 1999, 2000). In this paper, some new elements of the technique were used. We estimated the accuracy of observations. The paper contains information on the observing conditions and describes briefly the instruments employed.     

Two-integral distribution functions for axisymmetric systems

Some formulae are presented for finding two-integral distribution functions (DFs) which depend only on the two classical integrals of the energy and the magnitude of the angular momentum with respect to the axis of symmetry for stellar systems with known axisymmetric densities. They come from a combination of the ideas of Eddington and Fricke and they are also an extension of those shown by Jiang and Ossipkov for finding anisotropic DFs for spherical galaxies. The density of the system is required to be expressed as a sum of products of functions of the potential and of the radial coordinate. The solution corresponding to this type of density is in turn a sum of products of functions of the energy and of the magnitude of the angular momentum about the axis of symmetry. The product of the density and its radial velocity dispersion can be also expressed as a sum of products of functions of the potential and of the radial coordinate. It can be further known that the density multiplied by its rotational velocity dispersion is equal to a sum of products of functions of the potential and of the radial coordinate minus the product of the density and the square of its mean rotational velocity. These formulae can be applied to the Binney and the Lynden-Bell models. An infinity of the odd DFs for the Binney model can be also found under the assumption of the laws of the rotational velocity.     

Solar oscillations and the equation of state

Summary Accurate measurements of observed frequencies of solar oscillations are providing a wealth of data on the properties of the solar interior. The frequencies depend on solar structure, and on the properties of the plasma in the Sun. Here we consider in particular the dependence on the thermodynamic state. From an analysis of the equations of stellar structure, and the relevant aspects of the properties of the oscillations, we argue that in the convection zone one can isolate information about the equation of state which is relatively unaffected by other uncertainties in the physics of the solar interior. We review the different treatments that have been used to describe the thermodynamics of stellar plasmas. Through application of several of these to the computation of models of the solar envelope we demonstrate that the sensitivity of the observed frequencies is in fact sufficient to distinguish even quite subtle features of the physics of solar matter. This opens up the possibility of using the Sun as a laboratory for statistical mechanics, under conditions that are out of reach in a terrestrial laboratory.     

Dynamic Structure of the GLONASS and GPS Orbital Space: Problem of Disposal of Retired Objects

The paper presents the results of a study of the dynamic structure of the orbital space of the navigation systems GLONASS and GPS. It is shown that the dynamic structure of the GLONASS region is determined by the action of one stable Lidov–Kozai secular resonance. The motion of almost all the retired objects of the GLONASS system is stable throughout the 100-year study period. In the GPS region, there is an orbital resonance and a large number of secular resonances. Their combined influence leads to a rapid increase in the eccentricity of the orbits of the retired objects of the system. Features of the dynamic structure of the orbital space are used to find the graveyard (parking) orbits of the retired objects of navigation systems.     

The accuracy of observations of the Galilean satellites of Jupiter taken at the Abastumani Astrophysical Observatory of the Academy of Sciences of Georgia

The sets of photographic observations of the Galilean satellites of Jupiter taken at the Abastumani Astrophysical Observatory of the Academy of Sciences of Georgia are analyzed here. Positional observations of the system of Jupiter were made in the period from 1985 to 1994 with the use of the double Zeiss astrograph in order to determine the exact coordinates of Jupiter and its satellites. The accurate positions of the satellites and Jupiter itself, as well as their stellar (equatorial) coordinates relative to the stars of the currently available catalogs and the relative ??satellite ?? satellite?? coordinates were obtained from the observations. From the comparison of the observation results with the modern theories of motion of satellites, the accuracy in determining the positions of the satellites and Jupiter was analyzed. The results of observations are presented in the Pulkovo database of observations of Solar System bodies that is accessible to users at http://www.puldb.ru.     

On the formation of pulsar radiation diagrmas

A model of pulsars is discussed in which formation of a polar diagram of the radiation is influenced by the motion of the source around a neutron star with a velocity close to that of light. For a power-law frequency-spectrum of the radiation and isotropy of the diagram in a system of coordinates rotating with the source, the width of the observed pulse is shown to be independent of frequency.The proposed explanation of the second period characteristics of type CP 1919 pulsars is based on the effect of relativistic motion of the radiation source. The positions are established (relative to the axis of rotation of the star) of the local sources of radiation in the optical and in the radio ranges for the pulsar NP 0532. It is shown how the polarization characteristics of the optical radiation of this pulsar may be connected with the effects of relativistic orbiting of the source of radiation about the star.     

Closed solution of the Hill differential equation for short arcs and a local mass anomaly in the central body

The Hill differential equation describes the relative motion of a satellite w.r.t. a circular reference orbit. The deviations in the orbit are caused by a residual acceleration, which is small compared to the effect of the central gravitational field. In this paper, the acceleration of a local mass anomaly in the central body is considered, which rotates w.r.t. the inertial frame with a constant angular velocity. The mass anomaly is modeled by a superposition of radial base functions. The potential and the gradient of each base function are represented in the orbit by the Keplerian elements, the rotation rate of the central body and the parameters of the base function, i.e. the position of its center, the shape and a scaling factor. The inhomogeneous solution of the Hill differential equation for short arcs is found by means of the Laplace transform. A few lower orders of the solution require an additional Laplace transform, to consider the so-called resonance cases. The final deviations are described in a closed and differentiable formula of the Keplerian elements and the parameters of the base function.     

Simple Model of a Stochastically Excited Solar Dynamo

The aim of this paper is to investigate the dynamical nature of the complexity observed in the time evolution of the sunspot number. We report a detailed analysis of the sunspot number time series, and use the daily records to build the phase space of the underlying dynamical system. The observed features of the phase space prompted us to describe the global behavior of the solar cycle in terms of a noise-driven relaxation oscillator. We find the equations whose solutions best fit the observed series, which adequately describe the shape of the peaks and the oscillations of the system. The system of equations obtained from this fitting procedure is shown to be equivalent to a truncation of the dynamo equations. A linear transformation maps the phase space of these equations into the phase space reconstructed from the observations. The irregularities of the solar cycle were modeled through the introduction of a stochastic parameter in the equations to simulate the randomness arising in the process of eruption of magnetic flow to the solar surface. The mean values and deviations obtained for the periods, rise times and peak values, are in good agreement with the values obtained from the sunspot time series.     

Gasdynamical Treatment of Mass Transfer in Symbiotic Stars

The results of two-dimensional, non-adiabatic gas dynamical simulations of gas flows in symbiotic stars are presented. It is shown that for a binary system with components that do not fill their Roche lobes, the structure of the gaseous stream is determined not only by the flow from the vicinity of the inner Lagrangian point, but also by the flow of matter caused by the orbital motion of the accretor through the gas of the stellar wind.The calculated gaseous flow structure consists of a set of shocks and tangential discontinuities. It is found that for all considered cases two bow shocks exist. One of the bow shocks is located in front of the accretor in the path of its orbital motion, and the second one is between the components.It is found that matter is accreted in a spiral fashion, forming an accretion disk. In steady-state simulations, a structure with shocks (both trailing and leading) in the disk was observed. Gas dynamical perturbations of the disk's outer edge by the accreting gas inflow lead to the formation of two intensive trailing shocks propagating inward. Gas dynamical perturbations of the inner part of the disk result in the occurrence of leading shocks. One possible mechanism for the generation of spiral shocks is dissipative instability. The existence of spiral shocks results in a redistribution of the angular momentum of the disk material and an increase of the accretion rate.     

Reflectance model for densely packed media: Estimates of the surface properties of the high-albedo satellites of Saturn

Interpretation of photometric and polarimetric observations of atmosphereless celestial bodies faces the problems connected with both the insufficient accuracy and level of details in groundbased observations and the current state of the theory of the multiple scattering of light. In application to sparse media, where the electromagnetic waves, propagating between the scatterers, can be considered as spherical (the socalled far-field approximation), this theory is rather well developed for both the diffuse and coherent components of the scattered radiation. In this paper, we show that this theory can be also successfully applied to the measurements of polarization of light scattered by densely packed, though nonabsorbing or weakly absorbing, media. For this purpose, we calculated the models for a semi-infinite layer of the medium composed of randomly oriented clusters of spherical particles and compared them with the data of laboratory and astronomical measurements. The potential of the present approach is illustrated by an example of the interpretation of the polarization measurements of the ice satellites of Saturn—Rhea and Enceladus—which allowed some properties of the surface of these celestial bodies to be estimated. In particular, the ratio of the surface area that makes no contribution to the negative polarization of light reflected at small phase angles to the area producing the negative polarization branch was found. Under the assumption of the same albedo of these areas, this ratio turned out to be 3.31–3.66 and 1.7–3.8 for Rhea and Enceladus, respectively. For Enceladus, it is difficult to obtain a sufficiently narrow range of the estimated parameters, since the number of measurement points in the phase dependence of polarization of this satellite is small. For the surface of Rhea, the estimated packing density of particles, participating in the opposition effects, is approximately 15%, while their smallest size is of the order of the wavelength of visible light.     

Precise Determination of the Orientation of the Solar Image

Accurate heliographic coordinates of objects on the Sun have to be known in several fields of solar physics. One of the factors that affect the accuracy of the measurements of the heliographic coordinates is the accuracy of the orientation of a solar image. In this paper the well-known drift method for determining the orientation of the solar image is applied to data taken with a solar telescope equipped with a CCD camera. The factors that influence the accuracy of the method are systematically discussed, and the necessary corrections are determined. These factors are as follows: the trajectory of the center of the solar disk on the CCD with the telescope drive turned off, the astronomical refraction, the change of the declination of the Sun, and the optical distortion of the telescope. The method can be used on any solar telescope that is equipped with a CCD camera and is capable of taking solar full-disk images. As an example to illustrate the method and its application, the orientation of solar images taken with the Gyula heliograph is determined. As a byproduct, a new method to determine the optical distortion of a solar telescope is proposed.     

Photometric relief of the previously uninvestigated surface of mercury

The results of investigations of the surface relief of Mercury with a classical photometric method are reported. A subject of the photometric method is the intensity of light reflected by the surface of the planet. The main data for the photometric study are the high-resolution images of Mercury received from the MESSENGER space station during its first flyby over Mercury. The images of the surface of Mercury were downloaded from the NASA web-site (http://messenger.jhuapl.edu) and converted to a digital form for photometric measurements. The reflectance characteristics of the surface were calculated according to the model of a three-dimensional scattering phase function (Shevchenko, 1979; 2004a; 2006). From the photometric processing of the space-borne images, the reflectance of four morphologic types of the surface structure of Mercury was determined. With the Hapke model of the bidirectional reflectance (Hapke, 2001), the structural inhomogeneity of the regions of the morphologic forms with centimeter-sized rough-ness was estimated.     

Modeling of H II region radiation surrounding the starburst knot taking into account the evolution of structures formed by the superwind

The method for the multicomponent photoionization modeling of the H II region radiation surrounding the starburst knot is presented. The internal structure of the H II region has been determined using the evolutionary model of the superwind bubble from the starburst center. Models of Chevalier and Clegg (1985) and Weaver et al. (1977) have been used to determine the radial distribution of the gas density, the velocity of gas layers, and the temperature in the region of the superwind free expansion and in the cavity, respectively. The chemical content of the internal components of the bubble has been set by the results of the modeling of the evolutionary population synthesis. External components of our models describe a high-density, thin layer of gas formed by the shock wave of stellar superwind from the surrounding gas and a typical H II region, respectively. Input model parameters have been taken from the precalculated evolutionary starburst models based on three types of evolutionary tracks. Evolutionary grids of multicomponent low-metallicity models are calculated. A comparative analysis of the results of their calculation with the observational data has been carried out.     

On the storage of high-energy protons in the solar corona: The cyclotron instability

We consider the problem of long-time storage of high-energy protons, accelerated in the process of a flare, in coronal magnetic traps. From the viewpoint of the storage, one of the most important plasma instabilities is the kinetic cyclotron instability of the Alfvén waves. We carry out a detailed theoretical analysis of the instability for typical conditions of the solar corona. It is the refraction of the Alfvén waves in combination with a drastic decrease of the instability growth rate with an increase of the angle between the directions of the wave vector and the stationary magnetic field that leads to the possibility of the long-term storage of the flare protons. Sufficient conditions of the storage are determined.     

A mechanism for the growth phase of magnetospheric substorms

A mechanism is presented whereby the rate of energy dissipation in the magnetosphere is controlled by the particle density in the plasma sheet in the near geomagnetic tail. The mechanism is based on a model in which the plasma sheet is sustained by injection of solar-wind particles into the dayside magnetosphere. The efficiency of the injection is controlled by solarwind parameters, in particular, the north-south component of the interplanetary magnetic field; the maximum injection rate occurs when the interplanetary field is northward. During geomagnetically quiet times, this source balances the loss of particles from the edges of the tail current sheet. If the dayside source rate is reduced (e.g. by a southward-turning interplanetary magnetic field), then the plasma sheet is depleted and the rate of magnetic merging is enhanced in the earthward portion of the tail current sheet. This period of steadily-enhanced merging is associated with the growth phase, i.e. the period of enhanced magnetospheric convection for about one hour preceding the breakup of a polar magnetic or auroral substorm. The breakup can be understood as the result of the collapse of a portion of the tail current sheet following the local depletion of the plasma sheet.     

The struve dynasty in the history of astronomy in Ukraine

The impact of the Struve astronomical dynasty on the development of astronomy in Ukraine in the 19th–20th centuries is studied. First of all, the role of F.G.W. Struve and O.V. Struve in the formation of astronomical research programs at the observatories at the Kharkiv, Kyiv, Odesa, and Mykolayiv, in equipping the observatories with instruments, in practical training of astronomers as well as in the organization of astronomy-geodetic expeditions (19th century). Particular attention is paid to the activity of L.O. Struve as a director of the Astronomical observatory of Kharkiv University and his works conducted together with G.A. Shajn and B.P. Gerasimovich (20th century) as well as to the impact of his scientific and public activity, including one he made as a President of IAC, on the development of astronomy in the Soviet Union and Ukraine. A range of important documents from the archives of Russian Academy of Sciences, Institute of astronomy and State Archive of Ukraine are cited. The article was translated by the authors.     

On the Formation of Planetesimals: Radial Contraction of the Dust Layer Interacting with the Protoplanetary Disk Gas

Radial contraction of the dust layer in the midplane of a gas–dust protoplanetary disk that consists of large dust aggregates is modeled. Sizes of aggregates vary from centimeters to meters assuming the monodispersion of the layer. The highly nonlinear continuity equation for the solid phase of the dust layer is solved numerically. The purpose of the study is to identify the conditions under which the solid matter is accumulated in the layer, which contributes to the formation of planetesimals as a result of gravitational instability of the dust phase of the layer. We consider the collective interaction of the layer with the surrounding gas of the protoplanetary disk: shear stresses act on the gas in the dust layer that has a higher orbital velocity than the gas outside the layer, this leads to a loss of angular momentum and a radial drift of the layer. The stress magnitude is determined by the turbulent viscosity, which is represented as the sum of the α-viscosity associated with global turbulence in the disk and the viscosity associated with turbulence that is localized in a thin equatorial region comprising the dust layer and is caused by the Kelvin–Helmholtz instability. The evaporation of water ice and the continuity of the mass flux of the nonvolatile component on the ice line is also taken into account. It is shown that the accumulation of solid matter on either side of the ice line and in other regions of the disk is determined primarily by the ratio of the radii of dust aggregates on either side of the ice line. If after the ice evaporation the sizes (or density) of dust aggregates decrease by an order of magnitude or more, the density of the solid phase of the layer’s matter in the annular zone adjacent to the ice line from the inside increases sharply. If, however, the sizes of the aggregates on the inner side of the ice line are only a few times smaller than behind the ice line, then in the same zone there is a deficit of mass at the place of the modern asteroid belt. We have obtained constraints on the parameters at which the layer compaction is possible: the global turbulence viscosity parameter (α < 10^?5), the initial radial distribution of the surface density of the dust layer, and the distribution of the gas surface density in the disk. Restrictions on the surface density depend on the size of dust aggregates. It is shown that the timescale of radial contraction of a dust layer consisting of meter-sized bodies is two orders of magnitude and that of decimeter ones, an order of magnitude greater than the timescale of the radial drift of individual particles if there is no dust layer.