A grid search for families of periodic orbits in the restricted problem of three bodies

A method is described for the numerical determination of families of periodic orbits in the planar restricted problem of three bodies. The families are sought in their representation as curves in a two-dimensional space of parameters. A grid search is applied to the study of the evolution of satellite motion when the mass parameter is varied. Only that part of the space of parameters is investigated for which one of them, the relative energy constant, takes values larger than that corresponding to the inner Lagrangian pointL ₂. Critical values of the mass parameter are determined for which new families of simple or double periodic orbits appear inside the closed ovals of zero velocity.     

Third-order secular Solution of the variational equations of motion of a satellite in orbit around a non-spherical planet

We constructed an analytical theory of satellite motion up to the third order relative to the oblateness parameter of the Earth (J ₂). Equations of secular variations was developed for the first three orbital elements (a, e, i) of an artificial satellite. The secular variations are solved in a closed form.     

Nonuniqueness of families of periodic solutions in a four dimensional mapping

If we follow a family of periodic solutions along a closed path in a parameter space of two dimensions we may not return to the original solution when the parameters return to the original values. We study such nonuniqueness phenomena in simple and double period familis. Nonuniqueness appears if a closed path in the parameter space goes around a critical point. In some cases we find Riemann sheets in the same way as in multiply valued functions. In other cases the connections of various families change in a complicated way around the critical point. All these phenomena are explained analytically. At the critical point there is a collision of bifurcations. The changes of the connections of various families at such collisions of bifurcations are studied in some detail.     

Numerical continuation of families of frozen orbits in the zonal problem of artificial satellite theory

In the zonal problem of a satellite around the Earth, we continue numerically natural families of periodic orbits with the polar component of the angular momentum as the parameter. We found three families; two of them are made of orbits with linear stability while the third one is made of unstable orbits. Except in a neighborhood of the critical inclination, the stable periodic (or frozen) orbits have very small eccentricities even for large inclinations.     

Families of Asymmetric Periodic Orbits in the Restricted Three-body Problem

This paper studies the asymmetric solutions of the restricted planar problem of three bodies, two of which are finite, moving in circular orbits around their center of masses, while the third is infinitesimal. We explore, numerically, the families of asymmetric simple-periodic orbits which bifurcate from the basic families of symmetric periodic solutions f, g, h, i, l and m, as well as the asymmetric ones associated with the families c, a and b which emanate from the collinear equilibrium points L ₁, L ₂ and L ₃ correspondingly. The evolution of these asymmetric families covering the entire range of the mass parameter of the problem is presented. We found that some symmetric families have only one bifurcating asymmetric family, others have infinity number of asymmetric families associated with them and others have not branching asymmetric families at all, as the mass parameter varies. The network of the symmetric families and the branching asymmetric families from them when the primaries are equal, when the left primary body is three times bigger than the right one and for the Earth–Moon case, is presented. Minimum and maximum values of the mass parameter of the series of critical symmetric periodic orbits are given. In order to avoid the singularity due to binary collisions between the third body and one of the primaries, we regularize the equations of motion of the problem using the Levi-Civita transformations.     

Critical cases of 3-dimensional systems

We study the evolution of families of periodic orbits of simple 3-dimensional models representing the central parts of deformed galaxies. In some cases the evolution is non-unique, i.e. if we follow a closed path in the parameter space we do not return with the same periodic orbit. This happens when the path surrounds a critical point. We found that critical points are generated at particular collisions of bifurcations in limiting cases when the 3-D system is separated into a 2-D system and an independent oscillation along the third axis. The regions of stability and instability of some families of periodic orbits change in remarkable ways near the various collisions of bifurcations and around the critical points.     

Approximate General Solution Of The Two-Dimensional Duffing Dynamical System

This paper presents the approximate general solution of the triple well, double oscillator non-linear dynamical system. This system is non-integrable and the approximate general solution is calculated by application of the Last Geometric Theorem of Poincaré (Birkhoff, 1913, 1925). The original problem, known as the Duffing one, is a 1 degree of freedom system that, besides the conservative force component, includes dumping and external forcing terms (see details in the web site: http://www.uncwil.edu/people/hermanr/chaos/ted/chaos.html). The problem considered here is a 2 degree of freedom, autonomous and conservative one, without dumping, and of axisymmetric potential. The space of permissible motions is scanned for identification of all solutions re-entering after from one to nine oscillations and the precise families of periodic solutions are computed, including their stability parameter, covering all cases with periods T corresponding to 4osc/T. Seven sub-domains of the space of solutions were investigated in detail by zooming, an operation that proved the possibility to advance the accuracy of the approximate general solution to the level permitted by the integration routine. The approximation of the general solution, although impressive, provides clear evidence of the complexity of the problem and the need to proceed to larger period families. Nevertheless, it allows prediction of the areas where chaos and order regions in the Poincaré surfaces of section are to be expected. Examples of such surfaces of sections, as well as of types of closed solutions, are given. Two peculiar points of the space of solutions were identified as crossing, or source points from which infinite families of periodic solutions emanate. The morphology and stability of solutions of the problem are studied and discussed.     

Numerical investigation of the planar restricted three-body problem

In this article a method is described for the determination of families of periodic orbits, of the restricted problem of three bodies, as branchings of a given family of stable periodic orbits. Poincaré's method of successive crossings of a surface of section is applied for a value of the mass parameter corresponding to the Sun-Jupiter case of the restricted problem. New families are found, of the type of direct asteroids, having long periods and closing in space after many revolutions of the third body about the Sun. Their stability parameters are also given. The generating family, from which they branch, seems to have special significance for stability considerations.     

The manifolds of families of 3D periodic orbits associated to Sitnikov motions in the restricted three-body problem

This paper deals with the Sitnikov family of straight-line motions of the circular restricted three-body problem, viewed as generator of families of three-dimensional periodic orbits. We study the linear stability of the family, determine several new critical orbits at which families of three dimensional periodic orbits of the same or double period bifurcate and present an extensive numerical exploration of the bifurcating families. In the case of the same period bifurcations, 44 families are determined. All these families are computed for equal as well as for nearly equal primaries (μ = 0.5, μ = 0.4995). Some of the bifurcating families are determined for all values of the mass parameter μ for which they exist. Examples of families of three dimensional periodic orbits bifurcating from the Sitnikov family at double period bifurcations are also given. These are the only families of three-dimensional periodic orbits presented in the paper which do not terminate with coplanar orbits and some of them contain stable parts. By contrast, all families bifurcating at single-period bifurcations consist entirely of unstable orbits and terminate with coplanar orbits.     

Three-dimensional periodic oscillations generating from plane periodic ones around the collinear Lagrangian points

The three families of three-dimensional periodic oscillations which include the infinitesimal periodic oscillations about the Lagrangian equilibrium pointsL ₁,L ₂ andL ₃ are computed for the value =0.00095 (Sun-Jupiter case) of the mass parameter. From the first two vertically critical (|a _v |=1) members of the familiesa, b andc, six families of periodic orbits in three dimensions are found to bifurcate. These families are presented here together with their stability characteristics. The orbits of the nine families computed are of all types of symmetryA, B andC. Finally, examples of bifurcations between families of three-dimensional periodic solutions of different type of symmetry are given.     

Effects of solar activity on production rates of short‐lived cosmogenic radionuclides

The solar activity can be quantified by solar modulation parameter Φ that affects the heliospheric magnetic field. This activity influences the intensity of the galactic cosmic ray (GCR) particle flux within the solar system, and consequently, the differential primary particle spectra depend on the solar modulation parameter Φ (MeV). The modulation parameter Φ shows spatial and temporal variations (Leya and Masarik 2009). Some of the solar activity variations are cyclic and result in measurable effects as for example the 11‐year solar cycle. Variations in solar activity only induce small effects on the production of long‐lived cosmogenic radionuclides. This is due to the fact that activities measured in meteorites usually correspond to saturation values and represent long‐term average values. Long‐lived radionuclides often require millions of years of irradiation by GCR to reach saturation and therefore activity cycles average out. In contrast, one can expect strongly pronounced variations for saturation values caused by primary flux intensity variations, if short‐lived radionuclides with half‐lives ranging from days to a few years are investigated. Short‐lived cosmogenic nuclides were the subject of many experimental and theoretical investigations (e.g., Evans et al. 1982; Spergel et al. 1986; Neumann et al. 1997; Komura et al. 2002; Laubenstein et al. 2012). The aim of this work is to develop formulae for calculating production rates of radionuclides with short half‐life, taking into account temporal variations in the primary cosmic ray intensity. The developed formulae were applied to the Kosice and Chelyabinsk meteorites. The results for the Ko?ice meteorite were already published (Povinec et al. 2015). Here, we give a full explanation of underlying model.     

Periodic Orbits in the Restricted Three Body Problem with Radiation Pressure

This paper deals with the Restricted Three Body Problem (RTBP) in which we assume that the primaries are radiation sources and the influence of the radiation pressure on the gravitational forces is considered; in particular, we are interested in finding families of periodic orbits under theses forces. By means of some modifications to the method of numerical continuation of natural families of periodic orbits, we find several families of periodic orbits, both in two and three dimensions. As starters for our method we use some known periodic orbits in the classical RTBP. This revised version was published online in August 2006 with corrections to the Cover Date.     

Polarization of emerging resonance radiation: model profiles

The limiting polarization of a resonance line is examined for standard radiative transfer of polarized radiation in a semi-infinite scattering atmosphere with complete frequency redistribution. Two families of profiles of the line absorption coefficient, which are generalizations of Lorentz and Doppler profiles, are examined. It is shown that for both families this parameter approaches the Sobolev-Chandrasekhar limit when the fraction of absorption within the frequency interval (expressed in appropriate units) from −1 to 1 relative to the total absorption in the line approaches unity.     

Simple three-dimensional periodic orbits in a galactic-type potential

We study the families of simple periodic orbits in a three-dimensional system that represents the inner parts of a perturbed triaxial galaxy. The perturbations depend on two control parameters. We find the regions where each family is stable, simply unstable, doubly unstable, or complex unstable. the stable and simply unstable families produce other families by bifurcation. Several families reach a maximum (or minimum) perturbation and then are continued by other families. The bifurcations are direct or inverse. The transition from one type of bifurcation to the other is theoretically explained. Another important phenomenon is the splitting of one family into two, or the joining of two families into one. We do not have any complex instability in the limiting cases of two-dimensional motions (when one control parameter is zero).The two main families of periodic orbits are in most cases stable when the energy is smaller than the escape energy. Most high energy orbits are unstable. However, we found stable orbits even for energies about four times larger than the escape energy.     

Mid-infrared lightcurve of Vesta

We present a thermal mid-infrared lightcurve of Asteroid 4 Vesta and use this to infer variations in thermophysical properties over the surface. Vesta was observed over three nights during the May 2007 opposition with the Infrared Telescope Facility on Mauna Kea. Mid-infrared observations are compared to a model based on the Standard Thermal Model which is draped over a Vesta shape model derived from Hubble Space Telescope observations.A visible lightcurve with similar aspect was used to estimate the albedo as Vesta rotates. Shape and albedo can explain some of the features observed in the mid-infrared lightcurve. However, variations in the thermophysical properties, such as the “beaming parameter,” over Vesta’s surface are required to completely explain the observations.In order to match the mid-infrared magnitudes observed of Vesta, a beaming parameter of ∼0.862 is required which is higher than other Main Belt Asteroids such as Ceres and Pallas (0.756), indicating a smoother and/or rockier surface on Vesta. Variations in the beaming parameter with longitude are invoked to reproduce the observed thermal variations. Surface materials with relatively high beaming values, indicating a smoother and/or rockier surface, in the eastern hemisphere of Vesta coincide with locations where impact excavations may have produced surfaces that are younger and brighter relative to the western hemisphere.     

Evolution of families of double-and triple-periodic orbits

We study the evolution of the families of double-and triple-periodic orbits in a dynamical system that has closed zero velocity curves for arbitrarily large energies. We find three interesting features: (i) the characteristic x=x(h) of the family of double periodic orbits divides the (x,h)-plane into two unconnected parts; (i i) there is a sequence of sixteen closed characteristics, bifurcating from another one, each of them inside the previous one; (iii) inside the innermost characteristic of that sequence there is a sequence of eight pairs of close characteristics which are not connected with any of the previous characteristics.     

Restricted Three-Body Problem: An Approximation of its General Solution – Part One – the Manifold of Symmetric Periodic Solutions

This paper gives the results of a programme attempting to exploit ‘la seule bréche’ (Poincaré, 1892, p. 82) of non-integrable systems, namely to develop an approximate general solution for the three out of its four component-solutions of the planar restricted three-body problem. This is accomplished by computing a large number of families of ‘solutions précieuses’ (periodic solutions) covering densely the space of initial conditions of this problem. More specifically, we calculated numerically and only for μ = 0.4, all families of symmetric periodic solutions (1st component of the general solution) existing in the domain D:(x ₀ ∊ [−2,2],C ∊ [−2,5]) of the (x ₀, C) space and consisting of symmetric solutions re-entering after 1 up to 50 revolutions (see graph in Fig. 4). Then we tested the parts of the domain D that is void of such families and established that they belong to the category of escape motions (2nd component of the general solution). The approximation of the 3rd component (asymmetric solutions) we shall present in a future publication. The 4th component of the general solution of the problem, namely the one consisting of the bounded non-periodic solutions, is considered as approximated by those of the 1st or the 2nd component on account of the `Last Geometric Theorem of Poincaré' (Birkhoff, 1913). The results obtained provoked interest to repeat the same work inside the larger closed domain D:(x ₀ ∊ [−6,2], C ∊ [−5,5]) and the results are presented in Fig. 15. A test run of the programme developed led to reproduction of the results presented by Hénon (1965) with better accuracy and many additional families not included in the sited paper. Pointer directions construed from the main body of results led to the definition of useful concepts of the basic family of order n, n = 1, 2,… and the completeness criterion of the solution inside a compact sub-domain of the (x ₀, C) space. The same results inspired the ‘partition theorem’, which conjectures the possibility of partitioning an initial conditions domain D into a finite set of sub-domains D i that fulfill the completeness criterion and allow complete approximation of the general solution of this problem by computing a relatively small number of family curves. The numerical results of this project include a large number of families that were computed in detail covering their natural termination, the morphology, and stability of their member solutions. Zooming into sub-domains of D permitted clear presentation of the families of symmetric solutions contained in them. Such zooming was made for various values of the parameter N, which defines the re-entrance revolutions number, which was selected to be from 50 to 500. The areas generating escape solutions have being investigated. In Appendix A we present families of symmetric solutions terminating at asymptotic solutions, and in Appendix B the morphology of large period symmetric solutions though examples of orbits that re-enter after from 8 to 500 revolutions. The paper concludes that approximations of the general solution of the planar restricted problem is possible and presents such approximations, only for some sub-domains that fulfill the completeness criterion, on the basis of sufficiently large number of families.     

Bianchi type I two-fluid cosmological models with a variable G and Λ

This paper presents anisotropic, homogeneous two-fluid cosmological models in a Bianchi type I space–time with a variable gravitational constant G and cosmological constant Λ. In the two-fluid model, one fluid represents the matter content of the universe and another fluid is chosen to model the CMB radiation. We find a variety of solutions in which the cosmological parameter varies inversely with time t. We also discuss in detail the behavior of associated fluid parameters and kinematical parameters. This paper pictures cosmic history when the radiation and matter content of the universe are in an interactive phase. Here, Ω is closing to 1 throughout the cosmic evolution.      

The size distributions of asteroid families in the SDSS Moving Object Catalog 4

Asteroid families, traditionally defined as clusters of objects in orbital parameter space, often have distinctive optical colors. We show that the separation of family members from background interlopers can be improved with the aid of SDSS colors as a qualifier for family membership. Based on an ∼88,000 object subset of the Sloan Digital Sky Survey Moving Object Catalog 4 with available proper orbital elements, we define 37 statistically robust asteroid families with at least 100 members (12 families have over 1000 members) using a simple Gaussian distribution model in both orbital and color space. The interloper rejection rate based on colors is typically ∼10% for a given orbital family definition, with four families that can be reliably isolated only with the aid of colors. About 50% of all objects in this data set belong to families, and this fraction varies from about 35% for objects brighter than an H magnitude of 13 and rises to 60% for objects fainter than this. The fraction of C-type objects in families decreases with increasing H magnitude for H>13, while the fraction of S-type objects above this limit remains effectively constant. This suggests that S-type objects require a shorter timescale for equilibrating the background and family size distributions via collisional processing. The size distribution varies significantly among families, and is typically different from size distributions for background populations. The size distributions for 15 families display a well-defined change of slope and can be modeled as a “broken” double power-law. Such “broken” size distributions are twice as likely for S-type familes than for C-type families (73% vs. 36%), and are dominated by dynamically old families. The remaining families with size distributions that can be modeled as a single power law are dominated by young families (<1 Gyr). When size distribution requires a double power-law model, the two slopes are correlated and are steeper for S-type families. No such slope-color correlation is discernible for families whose size distribution follows a single power law. For several very populous families, we find that the size distribution varies with the distance from the core in orbital-color space, such that small objects are more prevalent in the family outskirts. This “size sorting” is consistent with predictions based on the Yarkovsky effect.     

RApid Temporal Survey (RATS) – I. Overview and first results

We present the aim and first results of the RApid Temporal Survey (RATS) made using the Wide Field Camera on the Isaac Newton Telescope. Our initial survey covers 3 square degrees, reaches a depth of   V ∼ 22.5  and is sensitive to variations on time-scales as short as 2 min: this is a new parameter space. Each field was observed for over 2 h in white light, with 12 fields being observed in total. Our initial analysis finds 45 targets which show significant variations. Around half of these systems show quasi-sinusoidal variations: we believe they are contact or short period binaries. We find four systems which show variations on a time-scale less than 1 h. The shortest period system has a period of 374 s. We find two systems which show a total eclipse. Further photometric and spectroscopic observations are required to fully identify the nature of these systems. We outline our future plans and objectives.