Generalized Hill's problem: Some cases of complete integrability

In this paper, we are investigating cases of integrability in the planar Hill's problem. The external potential U _extis supposed to be time independent in a given uniformly rotating frame. Cases of integrability of the relative motion of two interacting particles in the vicinity of an equilibrium solution of U _extare found. In all these cases, the form of the second integral is explicitly given, the first being the Jacobian one. Cases in which the interacting potential U between the two particles is of newtonian type are particularized.     

Stability of the classical type of relative equilibria of a rigid body in the J 2 problem

The motion of a point mass in the J ₂ problem is generalized to that of a rigid body in a J ₂ gravity field. The linear and nonlinear stability of the classical type of relative equilibria of the rigid body, which have been obtained in our previous paper, are studied in the framework of geometric mechanics with the second-order gravitational potential. Non-canonical Hamiltonian structure of the problem, i.e., Poisson tensor, Casimir functions and equations of motion, are obtained through a Poisson reduction process by means of the symmetry of the problem. The linear system matrix at the relative equilibria is given through the multiplication of the Poisson tensor and Hessian matrix of the variational Lagrangian. Based on the characteristic equation of the linear system matrix, the conditions of linear stability of the relative equilibria are obtained. The conditions of nonlinear stability of the relative equilibria are derived with the energy-Casimir method through the projected Hessian matrix of the variational Lagrangian. With the stability conditions obtained, both the linear and nonlinear stability of the relative equilibria are investigated in details in a wide range of the parameters of the gravity field and the rigid body. We find that both the zonal harmonic J ₂ and the characteristic dimension of the rigid body have significant effects on the linear and nonlinear stability. Similar to the classical attitude stability in a central gravity field, the linear stability region is also consisted of two regions that are analogues of the Lagrange region and the DeBra-Delp region respectively. The nonlinear stability region is the subset of the linear stability region in the first quadrant that is the analogue of the Lagrange region. Our results are very useful for the studies on the motion of natural satellites in our solar system.     

Analytical solutions for J 2-perturbed unbounded equatorial orbits

While solutions for bounded orbits about oblate spheroidal planets have been presented before, similar solutions for unbounded motion are scarce. This paper develops solutions for unbounded motion in the equatorial plane of an oblate spheroidal planet, while taking into account only the J ₂ harmonic in the gravitational potential. Two cases are distinguished: A pseudo-parabolic motion, obtained for zero total specific energy, and a pseudo-hyperbolic motion, characterized by positive total specific energy. The solutions to the equations of motion are expressed using elliptic integrals. The pseudo-parabolic motion unveils a new orbit, termed herein the fish orbit, which has not been observed thus far in the perturbed two-body problem. The pseudo-hyperbolic solutions show that significant differences exist between the Keplerian flyby and the flyby performed under the the J ₂ zonal harmonic. Numerical simulations are used to quantify these differences.     

The phase dependence of brightness and color of the lunar surface: a study based on integral photometric data

A procedure of an a posteriori correction of the available data on the integral photometry of the Moon is described. This procedure reduces the regular errors of the integral phase curves caused by variations of the libration parameters; the effect due to libration can reach 4%. A method allowing the integral measurements of the Moon to be compared correctly with the photometric measurements of the lunar areas or laboratory samples imitating the lunar soil has been developed. To approximate the phase curves of integral albedo in the phase-angle range from 6° to 120°, we proposed a simple empirical formula A _eq(α) = m _l e ^?ρα + m ₂ e ^?0.7α, where α is the phase angle, ρ is the factor of effective roughness, and m ₁ + m ₂ is the surface albedo at a zero phase angle. An empirical phase dependence of the slope of the lunar spectrum in the 360–1060 nm range has been obtained. The results may be used to test various theoretical models of the light scattering by the lunar surface and to calibrate the data of ground-based and space-borne spectrophotometric observations.     

Instability of the 2+2 body problem

The equations of motion of the 2+2 body problem (two interacting particles in the gravitational field of two much more massive primaries m₁ and m₂ in circular keplerian orbit) have an integral analogous to the Jacobi integral of the circular 2+1 body problem. We show here that with 2+2 bodies this integral does not give rise to Hill stability, i.e. to confinement for all time in a portion of the configuration space not allowing for some close approaches to occur. This is because all the level manifolds are connected and all exchanges of bodies between the regions surroundingm ₁,m ₂ and infinity do not contradict the conservation of the integral. However, it is worth stressing that some of these exchanges are physically meaningless, because they involve either unlimited extraction of potential energy from the binary formed by the small bodies (without taking into account their physical size) or significant mutual perturbations between the small masses without close approach, a process requiring, for the Sun-Jupiter-two asteroids system, timescales longer than the age of the Solar System.     

Motion in the generalized restricted three-body problem

This paper investigates the motion of an infinitesimal body in the generalized restricted three-body problem. It is generalized in the sense that both primaries are radiating, oblate bodies, together with the effect of gravitational potential from a belt. It derives equations of the motion, locates positions of the equilibrium points and examines their linear stability. It has been found that, in addition to the usual five equilibrium points, there appear two new collinear points L _n1, L _n2 due to the potential from the belt, and in the presence of all these perturbations, the equilibrium points L ₁, L ₃ come nearer to the primaries; while L ₂, L ₄, L ₅, L _n1 move towards the less massive primary and L _n2 moves away from it. The collinear equilibrium points remain unstable, while the triangular points are stable for 0<μ<μ _c and unstable for $\mu_{c} \le\mu\le\frac{1}{2}$ , where μ _c is the critical mass ratio influenced by the oblateness and radiation of the primaries and potential from the belt, all of which have destabilizing tendency. A practical application of this model could be the study of the motion of a dust particle near the oblate, radiating binary stars systems surrounded by a belt.     

Particle density in the plasmasphere during quiet periods

Density distribution of plasmaspheric particles in the equatorial plane is derived from a model of plasmaspheric streaming, which may produce S_q current system in the lower ionosphere, and from one integral of motion, which seems to be generally valid for steady-state magnetospheric convection. The results satisfy not only the observed features of S_q variation but also the observed pattern of the density distribution in the magnetospheric equatorial plane during quiet periods.     

Stablization and time transformation

We suggest a simple stabilization technique to reduce the along-track error in the numerical Integration of the Lagrange's equations of motion. We also investigate the equations of motion of the two-body problem after applying the Sundman transformation dt = r^αds, both with and without introducing the energy integral. In both cases, we show how the stability of the equations varies with α and in the case with the energy integral, we show that every solution is a quasi-periodic function of s with two frequencies.     

Schwarzschild-NUT solution in modified teleparallel gravity theory

Spherically symmetric tetrad field is applied to the field equation of modified teleparallel gravity theory. Some constraints are assumed on the resulting non-linear partial differential equations. Exact vacuum solution is derived with two constants of integration. The derived solution has a vanishing scalar torsion, i.e., T=T ^μ _νλ S _μ ^νλ =0 and its space-time is axially symmetric. It is shown that this solution can be rewritten as a product of three matrices, two local Lorentz transformations and one diagonal tetrad field. The local Lorentz transformations represent “so(3)” and a boost transformation and the diagonal tetrad reproduce the same metric field of the original tetrad.     

Asteroids,comets, meteoric matter: Proceedings of the IAU Colloquium No. 22, held at the University of Nice,France, April 4–6, 1972. Edited by C. Cristescu,W. J. Klepczynski,and B. Milet. Acad. Rep. Soc. Romania,Bucharest, 1974. 333 pp., 27 L

A model is proposed for single close encounters between two small masses, m₁and m₂, which orbit a much larger mass, M. The main new feature of the model is the assumption of conic motion of the center of mass of m₁and m₂ in the gravitational field of M. Comparisons of the model with the three-body equations of motion indicate that the model is a useful approximation for m₁, m₂ ? 10^?5M. The model is therefore applicable for encounters between bodies of the order of an earth mass or smaller in the presence of the sun. Comparisons are also made of outcomes obtained by the model with outcomes of numerical integration for a large variety of close encounters. The above comparisons reveal that for many purposes the model is an adequate approximation for those encounters with ? ≥ 4, where ? is the eccentricity of the hyperbolic orbit of m₁about m₂.     

Bianchi VI h with variable gravitational and cosmological constants

In order to study how the gravitational and the cosmological constants, G, Λ may vary, we consider two theoretical frameworks which are, a modification of the General Relativity and several scalar models (the standard, non-interacting and interacting models and their respective modifications to allow a G varying). We find exact self-similar solutions for the geometry Bianchi VI _h , (that is, the models: III, VI₀, and VI _h ,). Some physical and geometrical properties of the models are also discussed and we compare the obtained theoretical results with the current observational data. In the first of the theoretical models, we reach the conclusion that, from the structure of the field equations, the behaviour of Λ and G are related, but taking into account the observational data, we conclude that the Λ behaves as a positive decreasing time function while G is growing but in the long time regimen it tends to a constant value. In the scalar models, our solutions predict a “positive” dynamical cosmological constant in all the obtained solutions while the behaviour of G yields indeterminate, since its depends on a free parameter, G≈t ^2α , so it may be growing or decreasing as in the scalar-tensor theories.     

The Translational–Rotational Motion of an Earth Spheroid Satellite

In this paper we consider the translational–rotational motion of a spheroid satellite in the gravitational field, taking into account the asphericity of the earth. The harmonic coefficients of the earth’s gravitational field are taken up to J ₄. The equations of motion are obtained in terms of the canonical elements of Delaunay-Andoyer. A first order solution is obtained using the perturbing technique of Lie series.     

On the energy integral for first post-Newtonian approximation

The post-Newtonian approximation for general relativity is widely adopted by the geodesy and astronomy communities. It has been successfully exploited for the inclusion of relativistic effects in practically all geodetic applications and techniques such as satellite/lunar laser ranging and very long baseline interferometry. Presently, the levels of accuracy required in geodetic techniques require that reference frames, planetary and satellite orbits and signal propagation be treated within the post-Newtonian regime. For arbitrary scalar W and vector gravitational potentials \(W^j (j=1,2,3)\), we present a novel derivation of the energy associated with a test particle in the post-Newtonian regime. The integral so obtained appears not to have been given previously in the literature and is deduced through algebraic manipulation on seeking a Jacobi-like integral associated with the standard post-Newtonian equations of motion. The new integral is independently verified through a variational formulation using the post-Newtonian metric components and is subsequently verified by numerical integration of the post-Newtonian equations of motion.     

New exact solutions of Bianchi I,Bianchi III and Kantowski–Sachs spacetimes in scalar-coupled gravity theories via Noether gauge symmetries

The Noether symmetry approach is useful tool to restrict the arbitrariness in a gravity theory when the equations of motion are underdetermined due to the high number of functions to be determined in the ansatz. We consider two scalar-coupled theories of gravity, one motivated by induced gravity, the other more standard; in Bianchi I, Bianchi III and Kantowski–Sachs cosmological models. For these models, we present a full set of Noether gauge symmetries, which are more general than those obtained by the strict Noether symmetry approach in our recent work. Some exact solutions are derived using the first integrals corresponding to the obtained Noether gauge symmetries.     

Heat and mass transfer of an oscillatory flow with Hall current

In the first part of these notes new expressions—simpler than any previously obtained—are presented in integral form for the derivatives of the α _n ⁰ -functions (required for an interpretation of the observed light changes of eclipsing variables) with respect to the fractional radiir _{1, 2} and projected separation δ of their centres in terms of the modified Bessel functionsK _{0, 1} (x) of the second kind; and utilized for establishing new asymptotic formulae for the computation of ‘boundary integrals’ of the formJ _?1 ⁰ ,n(μ). In the second part of this paper, by a resort to bi-polar coordinates, we shall establish a new type of expansions for the α _n ⁰ -functions valid for any type of eclipses, and converging faster than the expansions of the cross-correlation integral of the form (1) for α _n ⁰ that have so far been established.     

Dust acoustic soliton and double layers with streaming dust and superthermal particles

Dust acoustic waves are investigated in plasma system containing dynamic and streaming dust, supertherrmal electrons and ions. Linear and nonlinear studies are carried out and elaborated with the help of parameters taken for Saturn’s F-ring. An energy integral equation is obtained by using the Sagdeev potential approach, and results are displayed by solving it analytically and numerically. The dependence of nonlinear structures on κ values, the ratio of electron to dust equilibrium densities μ _ed , Mach number M, and dust streaming speed v _d0 have been presented. The streaming speed appears as a destructive partner for the Mach number M in the pseudoenergy equation and hence plays a dominant modifying role in the formation of nonlinear structures. It plays a destructive role for some of the solitons and works as a source, for the emergence of new solitons (region). Formation of double layers are also investigated and shown that the amplitude, width and existence of double layers structures are predominantly affected by the presence of superthermal electrons, ions, and streaming dust beam.     

A Map for a Group of Resonant Cases in a quartic Galactic Hamiltonian

We present a map for the study of resonant motion in a potential made up of two harmonic oscillators with quartic perturbing terms. This potential can be considered to describe motion in the central parts of non-rotating elliptical galaxies. The map is based on the averaged Hamiltonian. Adding on a semi-empirical basis suitable terms in the unperturbed averaged Hamiltonian, corresponding to the 1:1 resonant case, we are able to construct a map describing motion in several resonant cases. The map is used in order to find thex − p _x Poincare phase plane for each resonance. Comparing the results of the map, with those obtained by numerical integration of the equation of motion, we observe, that the map describes satisfactorily the broad features of orbits in all studied cases for regular motion. There are cases where the map describes satisfactorily the properties of the chaotic orbits as well.     

On periodic orbits and resonance in extrasolar planetary systems

We study the evolution of an extrasolar planetary system with two planets, for planar motion, starting from an exact resonant periodic motion and increasing the deviation from the equilibrium solution. We keep the semimajor axes and the eccentricities of the two planets fixed and we change the initial conditions by rotating the orbit of the outer planet by Δω. In this way the resonance is preserved, but we deviate from the exact periodicity and there is a transition from order to chaos as the deviation increases. There are three different routes to chaos, as far as the evolution of (ω ₂ ? ω ₁) is concerned: (a) Libration → rotation → chaos, with intermittent transition from libration to rotation in between, (b) libration → chaos and (c) libration → intermittent interchange between libration and rotation → chaos. This indicates that resonant planetary systems where the angle (ω ₂ ? ω ₁) librates or rotates are not different, but are closely connected to the exact periodic motion.     

Exploring the nature of orbits in a galactic model with a massive nucleus

In the present article, we use an axially symmetric galactic gravitational model with a disk–halo and a spherical nucleus, in order to investigate the transition from regular to chaotic motion for stars moving in the meridian (r,z) plane. We study in detail the transition from regular to chaotic motion, in two different cases: the time independent model and the time evolving model. In both cases, we explored all the available range regarding the values of the main involved parameters of the dynamical system. In the time dependent model, we follow the evolution of orbits as the galaxy develops a dense and massive nucleus in its core, as mass is transported exponentially from the disk to the galactic center. We apply the classical method of the Poincaré (r,p_r) phase plane, in order to distinguish between ordered and chaotic motion. The Lyapunov Characteristic Exponent is used, to make an estimation of the degree of chaos in our galactic model and also to help us to study the time dependent model. In addition, we construct some numerical diagrams in which we present the correlations between the main parameters of our galactic model. Our numerical calculations indicate, that stars with values of angular momentum L_z less than or equal to a critical value L_zc, moving near to the galactic plane, are scattered to the halo upon encountering the nuclear region and subsequently display chaotic motion. A linear relationship exists between the critical value of the angular momentum L_zc and the mass of the nucleus M_n. Furthermore, the extent of the chaotic region increases as the value of the mass of the nucleus increases. Moreover, our simulations indicate that the degree of chaos increases linearly, as the mass of the nucleus increases. A comparison is made between the critical value L_zc and the circular angular momentum L_z0 at different distances from the galactic center. In the time dependent model, there are orbits that change their orbital character from regular to chaotic and vise versa and also orbits that maintain their character during the galactic evolution. These results strongly indicate that the ordered or chaotic nature of orbits, depends on the presence of massive objects in the galactic cores of the galaxies. Our results suggest, that for disk galaxies with massive and prominent nuclei, the low angular momentum stars in the associated central regions of the galaxy, must be in predominantly chaotic orbits. Some theoretical arguments to support the numerically derived outcomes are presented. Comparison with similar previous works is also made.