Numerical simulation of effects of the August 11, 1999 solar eclipse in the outer ionosphere

The results of a numerical simulation of such parameters of the topside ionosphere as concentration N _e and temperature T _e of electrons, and concentration n(H⁺) and fluxes along the magnetic field lines Φ(H⁺) of H⁺ ions at an altitude of ～2000 km for the conditions of the August 11, 1999 solar eclipse are presented. The calculations were performed using the Global Self-consistent Model of the Thermosphere, Ionosphere, and Protonosphere of the Earth (GSM TIP). It is shown that during the eclipse, in addition to a region of decreased values of T _e in the Northern Hemisphere and in the magnetically conjugate region in the Southern Hemisphere, regions of electron heating emerge in both hemispheres. Simultaneously, an extended region of decreased values of N _e comes into existence and moves behind the Moon’s shadow. Regions with decreased (down to ～30%) and enhanced (up to ～50%) concentrations of H⁺ ions are detected in the global distribution of these ions.     

A small parameter method in problems of maneuvering space vehicles with aerodynamic efficiency

The motion of a satellite with aerodynamic efficiency along a low near-circular orbit is considered in the paper. The controls of bank angle γ and lift coefficient C_y are used as control functions. The introduction of a small parameter ($\text{? = (}\text{?}{}_0\text{· S · g}{}_0\text{2G}\text{)}$) makes it possible to integrate an adjoint system of equations and to obtain an approximate solution to the complete problem in the class of piecewise-constant control functions. Maximum values for the coordinates of heading angle η and lateral derivation from the plane of a reference orbit ?, which are connected with orbit plane angle by the relation cos i = cos ? · cos h, are used as criteria of maneuvering capability for a satellite with aerodynamic efficiency. Optimal programs for bank angle and incidence variation are derived and the influence of lift-to-drag ratio on the vehicle maneuvering capabilities has been estimated.It is shown that the process of the optimal motion is a special kind of gravitational skipping similar to the Keplerian motion but with continuous descent.     

The optimization of the orbital Hohmann transfer

There are four bi-impulsive distinct configurations for the generalized Hohmann orbit transfer. In this case the terminal orbits as well as the transfer orbit are elliptic and coplanar. The elements of the initial orbit a₁, e₁ and the semi-major axis a₂ of the terminal orbit are uniquely given quantities. For optimization procedure, minimization is relevant to the independent parameter e_T, the eccentricity of the transfer orbit. We are capable of the assignment of minimum rocket fuel expenditure by using ordinary calculus condition of minimization for |ΔV_A|+|ΔV_B|=S.We exposed in detail the multi-steps of the optimization procedure. We constructed the variation table of S(e_T) which proved that S(e_T) is a decreasing function of e_T in the admissible interval [e_Tmin,e_Tmax]. Our analysis leads to the fact that e₂=1 for e_T=e_Tmax, i.e. the final orbit is a parabolic trajectory.     

An explanation of variable peculiarities in spectra of electrons and positrons of cosmic rays

It is shown that the excesses in the sum of fluxes (e ⁺ + e ^?) in cosmic rays in the energy range (10–1000) GeV and in the flux ratio e ⁺/(e ⁺ + e ^?) in the range > 10 GeV, observed both in recent and old experiments, can be explained by an accelerator of charged particles operating on the heliosphere periphery, in the region beyond the termination shock of the solar wind (～100 AU). Variations in the value and position of peculiarities in the spectra (e ⁺ + e ^?), as well as increasing ratio of fluxes e ⁺/(e ⁺ + e ^?), can be associated with variations of solar activity (and, as a consequence, of acceleration regimes) on different phases of the 11-years solar cycle. The results of numerical simulation of capture and acceleration of charged particles by packets of plasma waves in the heliospheric magnetic field are presented.     

The range of Alfvénic turbulence scales in the topside auroral ionosphere

The minimal scale of the Alfvénic turbulence transverse to the external magnetic field in the topside auroral ionosphere is investigated using electric field observations of the FAST spacecraft (the resolution 512 s^–1). The events in which the power law form of the electric fluctuation spectra with a 2.0–2.5 slope (typical of Alfvénic turbulence) remains unchanged down to acoustic gyroradius ρ_s or ion gyroradius ρ_i local values are illustrated for the first time. In this case, the character of spectrum variation does not change at the electron inertial length λ_e, which is much larger than ρ_s and ρ_i for FAST altitudes (apogee ~4000 km). We have tried to explain this experimental fact by consideration of the known scenarios of the appearance of a small transverse scale in an Alfvénic perturbation. It has been noted that the effects of front steepening in an inertial Alfvén wave with a finite amplitude, which propagates at an angle smaller than (m_e/m_i)^1/2 with respect to the transverse direction, can result in small transverse scales comparable with acoustic gyroradius appearing in a perturbation.     

Size quantization in high-temperature superconducting cuprates and a link to Einstein's diffusion law

Optimum doping of high-temperature superconductors (HTSC) defines a superconducting unit volume for each HTSC. For a single-mode HTSC, e.g., a cuprate with one CuO₂ plane, the volume is given by V_sc=cx², where c is the unit cell height and x the doping distance. The experimental resistivity at T_c is connected to the structure by ρ(exp)≈c×h/(2e²). Combining this result with the classical definition of resistivity leads to an equation similar to Einstein's diffusion law x²/(2τ)=h/(2M_eff)=D, where τ is the relaxation time, M_eff=2m_e and D the diffusion constant. It has also been shown that the mean free path d=x. The Einstein–Smoluchowski diffusion relation D=μk_BT_c provides a connection to T_c.     

Results of the early orbit manoeuvres of the ampte-uks

The Active Magnetospheric Particle Tracer Explorers (AMPTE) program consists of three satellites which were launched on 16th August 1984. The scientific aim of the mission is to inject lithium and barium tracer ions inside and outside the Earth's magnetosphere and to detect and monitor these ions as they diffuse through the inner magnetosphere. The first of these three satellites, the U.S. Charge Composition Explorer (CCE) was launched into an elliptical orbit of apogee 8 R_e. The other two satellites are the West German Ion Release Module (IRM) and the U.K. Subsatellite (UKS), both of which were launched on the same vehicle into a highly elliptical orbit of apogee 18 R_e. At discreet intervals during the mission the IRM will release ions into the solar wind, and the movement of these ions will be monitored by the UKS. Depending on the particular scientific requirement, the UKS has to be positioned accurately at a given distance behind the IRM. Initially the UKS has to be located 100 km behind the IRM, and held there for ～9 months. It will then be moved a distance of ～1 R_e behind the IRM. In order to manoeuvre the UKS around its orbit, a cold gas jet system is incorporated on the satellite, allowing impulses to be applied both along and perpendicular to the orbit velocity vector. The orbit control system also has to cater for relative orbit changes due to air drag at perigee, as the IRM and the UKS have different $\text{area}\text{mass}$ ratios. This paper presents an account of the orbit control system implemented on the UKS, together with the mathematical approach adopted, and results from manoeuvres made in the first weeks of the mission.     

Radial profile of pressure in a storm ring current as a function of D st

Using satellite data obtained near the equatorial plane during 12 magnetic storms with amplitudes from ?61 down to ?422 nT, the dependences of maximum in L-profile of pressure (L _m) of the ring current (RC) on the current value of D _st are constructed, and their analytical approximations are derived. It is established that function L _m(D _st ) is steeper on the phase of recovery than during the storm’s main phase. The form of the outer edge of experimental radial profiles of RC pressure is studied, and it is demonstrated to correspond to exponential growth of the total energy of RC particles on a given L shell with decreasing L. It is shown that during the storms’ main phase the ratio of plasma and magnetic field pressures at the RC maximum does not practically depend on the storm strength and L _m value. This fact reflects resistance of the Earth’s magnetic field to RC expansion, and testifies that during storms the possibilities of injection to small L are limited for RC particles. During the storms’ recovery phase this ratio quickly increases with increasing L _m, which reflects an increased fraction of plasma in the total pressure balance. It is demonstrated that function L _m(D _st ) is derived for the main phase of storms from the equations of drift motion of RC ions in electrical and magnetic fields, reflecting the dipole character of magnetic field and scale invariance of the pattern of particle convection near the RC maximum. For the recovery phase it is obtained from the Dessler-Parker-Sckopke relationship. The obtained regularities allow one to judge about the radial profile of RC pressure from ground-based magnetic measurements (data on the D _st variation).     

On Saturn’s rotation relative to a center of mass under the action of the gravitational moments of the Sun and Jupiter

Saturn’s rotation relative to a center of mass is considered within an elliptic restricted three-body problem. It is assumed that Saturn is a solid under the action of gravity of the Sun and Jupiter. The motions of Saturn and Jupiter are considered elliptic with small eccentricities e_S and e_J, respectively; the mean motion of Jupiter n_J is also small. We obtain the averaged Hamiltonian function for a small parameter of ε = n_J and integrals of evolution equations. The main effects of the influence of Jupiter on Saturn’s rotation are described: (α) the evolution of the constant parameters of regular precession for the angular momentum vector I₂; (β) the occurrence of new libration zones of oscillations I₂ near the plane of the celestial equator parallel to the plane of the Jupiter’s orbit; (γ) the occurrence of additional unstable equilibria of vector I₂ at the points of the north and south poles of the celestial sphere and, as a result, the existence of homoclinic trajectories; and (δ) the existence of periodic trajectories with arbitrarily large periods near the homoclinic trajectory. It is shown that the effects of (β), (γ), and (δ) are caused by the eccentricity e of the Jupiter’s orbit and are practically independent of Jupiter’s mass (within satellite approximation).     

A specific case of stability of cylindrical precession of a satellite

In a central Newtonian gravitational field, the motion of a dynamically symmetrical satellite along an elliptical orbit of arbitrary eccentricity is considered. The particular motion of the satellite is known when its axis of symmetry is perpendicular to the orbit plane, and the satellite rotates about this axis with a constant angular velocity (cylindrical precession). A nonlinear analysis of stability of this motion has been performed under the assumption that the geometry of the satellite mass corresponds to a thin plate. At small values of orbit eccentricity e the analysis is analytical, while numerical analysis is used for arbitrary values of e.     

Forecasting the velocity of quasi-stationary solar wind and the intensity of geomagnetic disturbances produced by it

A brief review is given of contemporary approaches to solving the problem of medium-term forecast of the velocity of quasi-stationary solar wind (SW) and of the intensity of geomagnetic disturbances caused by it. At the present time, two promising models of calculating the velocity of quasi-stationary SW at the Earth’s orbit are realized. One model is the semi-empirical model of Wang-Sheeley-Arge (WSA) which allows one to calculate the dependence V(t) of SW velocity at the Earth’s orbit using measured values of the photospheric magnetic field. This model is based on calculation of the local divergence f _S of magnetic field lines. The second model is semi-empirical model by Eselevich-Fainshtein-Rudenko (EFR). It is based on calculation in a potential approximation of the area of foot points on the solar surface of open magnetic tubes (sources of fast quasistationary SW). The new Bd-technology is used in these calculations, allowing one to calculate instantaneous distributions of the magnetic field above the entire visible surface of the Sun. Using predicted V(t) profiles, one can in EFR model calculate also the intensity of geomagnetic disturbances caused by quasi-stationary SW. This intensity is expressed through the K _p index. In this paper the EFR model is discussed in detail. Some examples of epignosis and real forecast of V(t) and K _p (t) are discussed. A comparison of the results of applying these two models for the SW velocity forecasting is presented.     

Luni-solar perturbations in the extended phase space representation of the Vinti problem

Perturbation theory is applied to the Vinti problem—motion about an oblate spheroid—to include the gravitational effects of the sun and moon. The problem is formulated using the extended phase space method which introduces a new independent variable similar to the true anomaly. The disturbing Hamiltonian H₁ for third bodies is of order J₂² (second order) and the final goal is a theory including second order short and long period terms and third order secular terms. The current paper however carries the development only to the second order in the secular terms and the first order in the periodic terms. Problems of including the higher orders are discussed. Therefore, in the development of H₁ all terms of order 10^?9 or larger are retained. The lunar emphemeris retains terms to e′² in the lunar eccentricity. The perturbation analysis is carried out by means of Lie series and is developed through the first order only which is consistent with the final accuracy desired. The generating function W₁ is obtained and separated into the long period, short period and secular terms. From W₁ the coordinates are defined from the Lie series by means of a transformation equation. These coordinates are non-singular for small eccentricity and inclination. Because of the complexity of the equations all algebraic computations were accomplished by means of a computerized Poisson series manipulator developed at the Naval Research Laboratory.     

Halo orbits in the vicinity of the L 2 libration point in the Sun-Earth system

Methods are proposed for constructing the orbits of spacecraft remaining for long periods of time in the vicinity of the L ₂ libration point in the Sun-Earth system (so-called halo orbits), and the trajectories of uncontrolled flights from low near-Earth orbits to halo orbits. Halo orbits and flight trajectories are constructed in two stages: A suitable solution to a circular restricted three-body problem is first constructed and then transformed into the solution for a restricted four-body problem in view of the real motions of the Sun, Earth, and Moon. For a halo orbit, its prototype in the first stage is a combination of a periodic Lyapunov solution in the vicinity of the L ₂ point and lying in the plane of large-body motion, with the solution for the linear second-order system describing small deviations of the spacecraft from this plane along the periodic solution. The desired orbit is found as the solution to the three-body problem best approximating the prototype in the mean square. The constructed orbit serves as a similar prototype in the second stage. In both stages, the approximating solution is constructed by continuation along a parameter that is the length of the approximation interval. Flight trajectories are constructed in a similar manner. The prototype orbit in the first stage is a combination of a solution lying in the plane of large-body motion and a solution for a linear second-order system describing small deviations of the spacecraft from this plane. The planar solution begins near the Earth and over time tends toward the Lyapunov solution existing in the vicinity of the L ₂ point. The initial conditions of both prototypes and the approximating solutions correspond to the spacecraft’s departure from a low near-Earth orbit at a given distance, perigee, and inclination.     

Results of measurements with the Planetary Fourier Spectrometer onboard Mars Express: Clouds and dust at the end of southern summer. A comparison with OMEGA images

We discuss the results of measurements made with the Planetary Fourier Spectrometer (PFS) onboard the Mars Express spacecraft. The data were obtained in the beginning of the mission and correspond to the end of summer in the southern hemisphere of Mars (L _s ～ 340°). Three orbits are considered, two of which passed through volcanoes Olympus and Ascraeus Mons (the height above the surface is about +20 km), while the third orbit intersects lowland Hellas (?7 km). The influence of the relief on the properties of the aerosol observed is demonstrated: clouds of water ice with a visual optical thickness of 0.1–0.5 were observed above volcanoes, while only dust was found during the observations (close in time) along the orbit passing through Hellas in low and middle latitudes. This dust is homogeneously mixed with gas and has a reduced optical thickness of 0.25±0.05 (at v = 1100 cm^?1). In addition to orographic clouds, ice clouds were observed in this season in the northern polar region. The clouds seen in the images obtained simultaneously by the mapping spectrometer OMEGA confirm the PFS results. Temperature inversion is discovered in the north polar hood below the level 1 mbar with a temperature maximum at about 0.6 mbar. This inversion is associated with descending movements in the Hadley cell.     

Collisions of stars by oscillating orbits of the second kind

Ordinary estimations of the number of star collisions in our galaxy—by simple kinematic considerations—lead to a very small number of such collisions: about one or even less every millions of years. However star collisions can occur through the following indirect way which has a much higher probability. (a) Binary stars are very common in our galaxy, about 30–50% of the stars. (b) If two binary stars meet a triple system can be formed by an ordinary exchange type motion. (c) A triple system is generally decomposed into the “inner orbit” (i.e. the relative orbit of the two nearest stars) and the “outer orbit” (i.e. the relative orbit of the third star with respect to the center of mass of the two nearest stars). The major axes of these two orbits have generally small perturbations and it is the same for the eccentricity of the outer orbit. On the contrary, if the relative inclination of the two orbits is large, the perturbations of the eccentricity of the inner orbit are important and can even in some cases lead to an eccentricity equal to one, that is to a collision of the two stars of the inner orbit.Such orbits can be called “oscillating orbits of the second kind”, indeed the first oscillating orbits—conceived by Khilmi and described for the first time in an example by Sitnikov—have unbounded mutual distances r_ij, but the system always come back to small sizes, it has an infinite number of very large expansions followed by strong contractions and, in the three-body case, an upper bound of lim inf (r_1.2 + r_1.3 + r_2.3) can be given in terms of the three masses and the integrals of motion. For the oscillating orbits of the second kind the mutual distances r_ij are bounded, but the velocities are unbounded (i.e. lim inf r_ij = 0 for at least one r_ij) and the system goes to a collision if the bodies have non-zero radius even small. The analytical study of the oscillating orbits of the second kind is a part of the general analytical study of the three-body problem, a part which must be valid for large eccentricities and large inclinations. The use of Delaunay's variables and of a Von Zeipel transformation lead to a first order integrable approximation, valid for any eccentricities and any inclinations, and giving the following results: (a) The oscillating orbits of the second kind occur when the angular momentum of the outer orbit has a modulus sufficiently close to the modulus of the total angular momentum of the three-body system. Hence these orbits occur for inclinations in the vicinity of 90°. (b) The oscillating orbits represent a set of positive measure of phase space and the first order study allows to give a rough estimation of the probability of collisions—even for stars of infinitely small radius. This probability, for given initial major axes and eccentricities and for isotropic arbitrary initial orientations, is generally of the order of m₃$\text{R}\text{M}$ (m₃ being the mass of the outer star, M the total mass and R the ratio of the period of the inner orbit to that of outer orbit).One question remains to be solved: how many collisions of stars are due to that phenomenon? That question is difficult because the probability of formation of a triple system by a random meeting of two binaries is very uneasy to estimate. However it seems that, compared to the usual evaluations based on pure kinematic considerations without gravitational effects, the number of collisions must be multiplied by a factor between one thousand and one million.     

Peculiarities of long-wave radio bursts from solar flares preceding strong geomagnetic storms

Radio bursts in the frequency range of 100–1500 kHz, recorded in 1997–2000 on the INTERBALL-1 satellite during the solar flares preceding the strong geomagnetic storms with D _st < ?100 nT, are analyzed in this paper. The observed long-wave III-type radio bursts of solar origin at frequencies of 1460 and 780 kHz were characterized by large values of the flux S _f = 10^?15 ?10^?17 W/m² Hz and duration longer than 10 min. The rapid frequency drift of a modulated radio burst continued up to a frequency of 250 kHz, which testified that the exciting agent (a beam of energetic electrons) propagated from the Sun to the Earth. All such flares were characterized by the appearance of halo coronal mass ejections, observed by the LASCO/SOHO, and by the presence of a southward B_z-component of the IMF, measured on the ACE and WIND spacecraft. In addition, shortly after radio bursts, the INTERBALL-1 satellite has recorded the fluxes of energetic electrons with E > 40 keV.     

Oxygen deficiency structure in iron-based high temperature superconductor GdFeAsO1−δ

Oxygen deficiency in the iron-based HTSC GdFeAsO_1?δ seems to create a parallelogram shaped Fe²⁺-ion/oxygen deficiency pattern in the Fe₂O₂ plane in c-direction. These two-dimensional nanostructures form superconducting current channels which are separated by h=0.828 nm. The doping distance in direction of the super-current shows a strong correlation to the transition temperature.     

Analytical solutions for the relative motion of spacecraft subject to Lorentz-force perturbations

A spacecraft capable of producing higher-than-natural electrostatic charges may achieve propellantless orbital maneuvering via the Lorentz-force interaction with a planetary magnetic field. Development of maneuver strategies for these propellantless vehicles is complicated by the fact that the perturbative Lorentz force acts along only a single line of action at any instant. Relative-motion dynamical models are developed that lead to approximate analytical solutions for the motion of charged spacecraft subject to the Lorentz force. These solutions indicate that the principal effects of the Lorentz force on a spacecraft in a circular orbit are to change the intrack position and to change the orbit plane. A rendezvous example is presented in which a spacecraft with a specific charge of ?3.81 $\times$ 10^?4 C/kg reaches a target vehicle initially 10 km away (on the same equatorial low-Earth orbit) in 1 day. Fly-around maneuvers may be achieved in low-Earth orbit with specific charges on the order of 0.001 C/kg.     

Low-energy Earth–Moon transfers involving manifolds through isomorphic mapping

Analysis and design of low-energy transfers to the Moon has been a subject of great interest for decades. Exterior and interior transfers, based on the transit through the regions where the collinear libration points are located, have been studied for a long time and some space missions have already taken advantage of the results of these studies. This paper is concerned with a geometrical approach for low-energy Earth-to-Moon mission analysis, based on isomorphic mapping. The isomorphic mapping of trajectories allows a visual, intuitive representation of periodic orbits and of the related invariant manifolds, which correspond to tubes that emanate from the curve associated with the periodic orbit. Two types of Earth-to-Moon missions are considered. The first mission is composed of the following arcs: (i) transfer trajectory from a circular low Earth orbit to the stable invariant manifold associated with the Lyapunov orbit at L₁ (corresponding to a specified energy level) and (ii) transfer trajectory along the unstable manifold associated with the Lyapunov orbit at L₁, with final injection in a periodic orbit around the Moon. The second mission is composed of the following arcs: (i) transfer trajectory from a circular low Earth orbit to the stable invariant manifold associated with the Lyapunov orbit at L₁ (corresponding to a specified energy level) and (ii) transfer trajectory along the unstable manifold associated with the Lyapunov orbit at L₁, with final injection in a capture (non-periodic) orbit around the Moon. In both cases three velocity impulses are needed to perform the transfer: the first at an unknown initial point along the low Earth orbit, the second at injection on the stable manifold, the third at injection in the final (periodic or capture) orbit. The final goal is in finding the optimization parameters, which are represented by the locations, directions, and magnitudes of the velocity impulses such that the overall delta-v of the transfer is minimized. This work proves how isomorphic mapping (in two distinct forms) can be profitably employed to optimize such transfers, by determining in a geometrical fashion the desired optimization parameters that minimize the delta-v budget required to perform the transfer.     

Energy dependence of the characteristic decay time of proton fluxes in solar cosmic ray events

Energetic solar proton events within the energy interval 1–48 MeV at the stage of their decay are considered over the period of 1974–2001. The dependence of the characteristic decay time on the proton energy in the assumed power-law representation τ(E) =E ^?n is analyzed for the events with an exponential decay form. The dependence of n on the heliolongitude of the flare (the particles source on the Sun) is studied.