Families of orbits in planar anisotropic fields

The aim of the planar inverse problem of dynamics is: given a monoparametric family of curves f(x, y) = c, find the potential V (x, y) under whose action a material point of unit mass can describe the curves of the family. In this study we look for V in the class of the anisotropic potentials V(x, y) = v(a²x² + y²), (a=constant). These potentials have been used lately in the search of connections between classical, quantum, and relativistic mechanics. We establish a general condition which must be satisfied by all the families produced by an anisotropic potential. We treat special cases regarding the families (e. g. families traced isoenergetically) and we present certain pertinent examples of compatible pairs of families of curves and anisotropic potentials. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Programmed motion in the presence of homogeneity

In the framework of the inverse problem of dynamics, we face the following question with reference to the motion of one material point: Given a region T_orb of the xy plane, described by the inequality g (x, y) ≤ c₀, are there potentials V = V (x, y) which can produce monoparametric families of orbits f (x, y) = c (also to be found) lying exclusively in the region T_orb? As the relevant PDEs are nonlinear, an answer to this question (generally affirmative, but not with assurance) can be given by the procedure of the determination of certain constants specifying the pertinent functions. In this paper we ease the mathematics involved by making certain simplifying assumptions referring to the homogeneity of both the function g (x, y) (describing the boundary of T_orb) and of the slope function γ(x, y) = f_y/f_x (representing the required family f (x, y) = c). We develop the method to treat the so formulated problem and we show that, even under these restrictive assumptions, an affirmative answer is guaranteed provided that two algebraic equations have in common at least one solution (© 2009 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

On the adelphic potentials compatible with a set of planar orbits

Given a planar potentialB=B(x, y), compatible with a monoparametric family of planar orbitsf(x, y)=c, we face the problem of producing potentialsA=A(x, y), adelphic toB(x, y), i.e. nontrivial potentials which have in common withB(x, y) the given set of orbits. We establish a linear, second order partial differential equation for a functionP(x, y) and we prove that, to any definite positive solution of this equation, there corresponds a potentialA(x, y) adelphic toB(x, y).     

Special Families of Orbits in the Direct Problem of Dynamics

The direct problem of dynamics in two dimensions is modeled by a nonlinear second-order partial differential equation, which is therefore difficult to be solved. The task may be made easier by adding some constraints on the unknown function = f _y /f _x , where f(x, y) = c is the monoparametric family of orbits traced in the xy Cartesian plane by a material point of unit mass, under the action of a given potential V(x, y). If the function is supposed to verify a linear first-order partial differential equation, for potentials V satisfying a differential condition, can be found as a common solution of certain polynomial equations.The various situations which can appear are discussed and are then illustrated by some examples, for which the energy on the members of the family, as well as the region where the motion takes place, are determined. One example is dedicated to a Hénon—Heiles type potential, while another one gives rise to families of isothermal curves (a special case of orthogonal families). The connection between the inverse/direct problem of dynamics and the possibility of detecting integrability of a given potential is briefly discussed.This revised version was published online in October 2005 with corrections to the Cover Date.     

Motion of rigid bodies in a set of redundant variables

In this article we study the conditions for obtaining canonical transformationsy=f(x) of the phase space, wherey(y ₁,y ₂,...,y _2n ) andx(x ₁,x ₂,...,x _2m ) in such a way that the number of variables is increased. In particular, this study is applied to the rotational motion in functions of the Eulerian parameters (q ₀,q ₁,q ₂,q ₃) and their conjugate momenta (Q ₀,Q ₁,Q ₂,Q ₃) or in functions of complex variables (z ₁,z ₂,z ₃,z ₄) and their conjugate momenta (Z ₁,Z ₂,Z ₃,Z ₄) defined by means of the previous variables. Finally, our article include some properties on the rotational motion of a rigid body moving about a fixed point.     

Families of planar orbits generated by homogeneous potentials

The second order partial differential equation which relates the potentialV(x,y) to a family of planar orbitsf(x,y)=c generated by this potential is applied for the case of homogeneousV(x,y) of any degreem. It is shown that, if the functionf(x,y) is also homogeneous, there exists, for eachm, a monoparametric set of homogeneous potentials which are the solutions of an ordinary, linear differential equation of the second order. Iff(x,y) is not homogeneous, in general, there is not a homogeneous potential which can create the given family; only if =f _y /f _x satisfies two conditions, a homogeneous potential does exist and can be determined uniquely, apart from a multiplicative constant. Examples are offered for all cases.     

Compatible potentials and orbits in synodic systems

For a given family of orbits f(x,y) = c ^* which can be traced by a material point of unit in an inertial frame it is known that all potentials V(x,y) giving rise to this family satisfy a homogeneous, linear in V(x,y), second order partial differential equation (Bozis,1984). The present paper offers an analogous equation in a synodic system Oxy, rotating with angular velocity . The new equation, which relates the synodic potential function (x,y), = –V(x, y) + % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaSqaaSqaai% aaigdaaeaacaaIYaaaaaaa!3780!\[\tfrac{1}{2}\]²(x ² + y ²) to the given family f(x,y) = c ^*, is again of the second order in (x,y) but nonlinear.As an application, some simple compatible pairs of functions (x,y) and f(x, y) are found, for appropriate values of , by adequately determining coefficients both in and f.     

Family boundary curves for autonomous dynamical systems

The notion of the family boundary curves (FBC), introduced recently for two-dimensional conservative systems, is extended to account for, generally, nonconservative autonomous systems of two degrees of freedom. Formulae are found for the force componentsX (x, y),Y (x, y) which produce a preassigned family of orbitsf(x, y)=c lying inside a preassigned, open or closed, regionB(x, y)0 of the xy plane.     

Stability of motion in the Sitnikov 3-body problem

We study the stability of motion in the 3-body Sitnikov problem, with the two equal mass primaries (m ₁ = m ₂ = 0.5) rotating in the x, y plane and vary the mass of the third particle, 0 ≤ m ₃ < 10⁻³, placed initially on the z-axis. We begin by finding for the restricted problem (with m ₃ = 0) an apparently infinite sequence of stability intervals on the z-axis, whose width grows and tends to a fixed non-zero value, as we move away from z = 0. We then estimate the extent of “islands” of bounded motion in x, y, z space about these intervals and show that it also increases as |z| grows. Turning to the so-called extended Sitnikov problem, where the third particle moves only along the z-axis, we find that, as m ₃ increases, the domain of allowed motion grows significantly and chaotic regions in phase space appear through a series of saddle-node bifurcations. Finally, we concentrate on the general 3-body problem and demonstrate that, for very small masses, m ₃ ≈ 10⁻⁶, the “islands” of bounded motion about the z-axis stability intervals are larger than the ones for m ₃ = 0. Furthermore, as m ₃ increases, it is the regions of bounded motion closest to z = 0 that disappear first, while the ones further away “disperse” at larger m ₃ values, thus providing further evidence of an increasing stability of the motion away from the plane of the two primaries, as observed in the m ₃ = 0 case.     

Motion on two-dimensional manifolds in rotating coordinates

The two degree-of-freedom system in rotating coordinates: \.u – 2nv = V _x, \.v + 2nu = V _y, \.x = u, \.y = v and its Jacobi integral define a time-dependent velocity field on a differentiable, two-dimensional manifold of integral curves. Explicit time dependence is determined by the dynamical system, coordinate frame, and initial conditions. In the autonomous cases, orbits are level curves of an autonomous function satisfying a second-order, quasi-linear, partial differential equation of parabolic type. Important aspects of the theory are illustrated for the two-body problem in rotating coordinates.     

Study of faint galaxies in the field of GRB021004

We present an analysis of BV R _c I _c observations of the field sized around 4′ × 4′ centered at the host galaxy of the gamma-ray burst GRB021004 with the 6-m BTA telescope of the Special Astrophysical Observatory of the Russian Academy of Sciences. We measured the magnitudes and constructed the color diagrams for 311 galaxies detected in the field (S/N>3). The differential and integral counts of galaxies up to the limit, corresponding to 28.5 (B), 28.0 (V), 27.0 (R _c ), 26.5 (I _c ) were computed. We compiled the galaxy catalog, consisting of 183 objects, for which the photometric redshifts up to the limiting magnitudes 26.0 (B), 25.5 (V), 25.0 (R _c ), 24.5 (I _c ) were determined using the HyperZ code. We then examined the radial distribution of galaxies based on the z estimates. We have built the curves expected in the case of a uniform distribution of galaxies in space, and obtained the estimates for the size and contrast of the possible super-large-scale structures, which are accessible with the observations of this type.     

The Parker Problem and the Theory of Coronal Heating

To illustrate his theory of coronal heating, Parker initially considers the problem of disturbing a homogeneous vertical magnetic field that is line-tied across two infinite horizontal surfaces. It is argued that, in the absence of resistive effects, any perturbed equilibrium must be independent of z. As a result random footpoint perturbations give rise to magnetic singularities, which generate strong Ohmic heating in the case of resistive plasmas. More recently these ideas have been formalized in terms of a magneto-static theorem but no formal proof has been provided. In this paper we investigate the Parker hypothesis by formulating the problem in terms of the fluid displacement. We find that, contrary to Parker's assertion, well-defined solutions for arbitrary compressibility can be constructed which possess non-trivial z-dependence. In particular, an analytic treatment shows that small-amplitude Fourier disturbances violate the symmetry ∂_z = 0 for both compact and non-compact regions of the (x, y) plane. Magnetic relaxation experiments at various levels of gas pressure confirm the existence and stability of the Fourier mode solutions. More general footpoint displacements that include appreciable shear and twist are also shown to relax to well-defined non-singular equilibria. The implications for Parker's theory of coronal heating are discussed.     

The magnetic structures of electron phase-space holes formed in the electron two-stream instability

It is well known that the parallel cuts of the parallel and perpendicular electric field in electron phase-space holes (electron holes) have bipolar and unipolar structures, respectively. Recently, electron holes in the Earth’s plasma sheet have been observed by THEMIS satellites to have detectable fluctuating magnetic field with regular structures. Du et al. (2011) investigated the evolution of a one-dimensional (1D) electron hole with two-dimensional (2D) electromagnetic particle-in-cell (PIC) simulations in weakly magnetized plasma (Ω _e <ω _pe , where Ω _e and ω _pe are the electron gyrofrequency and electron plasma frequency, respectively), which initially exists in the simulation domain. The electron hole is unstable to the transverse instability and broken into several 2D electron holes. They successfully explained the observations by THEMIS satellites based on the generated magnetic structures associated with these 2D electron holes. In this paper, 2D electromagnetic particle-in-cell (PIC) simulations are performed in the x–y plane to investigate the nonlinear evolution of the electron two-stream instability in weakly magnetized plasma, where the background magnetic field (B₀ = B₀[(e)\vec] _x)(\mathbf{B}_{0} =B_{0}\vec{\mathbf{e}} _{x}) is along the x direction. Several 2D electron holes are formed during the nonlinear evolution, where the parallel cuts of E _x and E _y have bipolar and unipolar structures, respectively. Consistent with the results of Du et al. (2011), we found that the current along the z direction is generated by the electric field drift motion of the trapped electrons in the electron holes due to the existence of E _y , which produces the fluctuating magnetic field δB _x and δB _y in the electron holes. The parallel cuts of δB _x and δB _y in the electron holes have unipolar and bipolar structures, respectively.     

Intergalactic medium: X-ray background and compton parameter

A model of intergalactic medium heated by QSOs and cooled by the expansion of the universe and Compton cooling is studied in the framework of a Friedmann-Robertson-Walker universe. Cosmological evolution functions of the comoving density of QSO's as well as the case of no evolution are considered. The theoretical X-ray background spectrum (through thermal bremsstrahlung) and Comptony parameter are calculated including relativistic corrections in the electron-electron, electron-proton and electron-photon interactions. The observed X-ray background and the upper limit of the Compton parametery cobe given by the COBE satellite are used to adjust, for each value of reheating redshiftsz _c ranging from 0.1 to 5.0, the present values of the temperatureT ₀ and densityn ₀ of the intergalactic gas. Forz _c > 0.25, when the theoretical X-ray spectrum fits the observed one, the adjusted values ofT ₀ andn ₀ imply iny >y cobe. On the other hand, whenT ₀ andn ₀ are consistent withy cobe, the calculated X-ray spectrum is lower than the observed one. Unless 100% of the observed X-ray background is due to discrete sources and if the intergalactic medium contributes more than 2.5% to such background we come to the interesting result that the medium must have been heated atz _c < 1. In this case we shall have to explain the high energy rates necessary to heat the intergalactic medium. Forz _c 0.25, it is possible to find values ofT ₀ andn ₀ such that both the calculated X-ray background and the y parameter simultaneously reproduce the corresponding observed values. However, in this case, unless it could be shown to be otherwise by future observations or theoretical studies, it seems that the model of hot intergalactic medium is not plausible because of the high energies required to heat the intergalactic gas.     

“Stickiness” in mappings and dynamical systems

We present results of a study of the so-called “stickiness” regions where orbits in mappings and dynamical systems stay for very long times near an island and then escape to the surrounding chaotic region. First we investigated the standard map in the form x_i+1 = x_i+y_i+1 and y_i+1 = y_i+K/2π · sin(2πx_i) with a stochasticity parameter K = 5, where only two islands of regular motion survive. We checked now many consecutive points—for special initial conditions of the mapping—stay within a certain region around the island. For an orbit on an invariant curve all the points remain forever inside this region, but outside the “last invariant curve” this number changes significantly even for very small changes in the initial conditions. In our study we found out that there exist two regions of “sticky” orbits around the invariant curves: A small region I confined by Cantori with small holes and an extended region II is outside these cantori which has an interesting fractal character. Investigating also the Sitnikov-Problem where two equally massive primary bodies move on elliptical Keplerian orbits, and a third massless body oscillates through the barycentre of the two primaries perpendicularly to the plane of the primaries—a similar behaviour of the stickiness region was found. Although no clearly defined border between the two stickiness regions was found in the latter problem the fractal character of the outer region was confirmed.     

Volume and mass distribution in selected asteroid families

The main focus of this paper is calculation of the diameters of asteroids belonging to five families (Vesta, Eos, Eunomia, Koronis, and Themis). To do that, we used the HCM algorithm applied for a data set containing 292,003 numbered asteroids, and a numerical procedure for choosing the crucial parameter of the HCM, called “the cutting velocity” v_cut. It was established with a precision as high as 1 m s^?1. Thereafter, we used the WISE (Wide‐field Infrared Survey Explorer) catalog to set a range of albedo for the largest members of each family considered. The albedo data were supported by the data concerning color classification (SDSS MOC4). The asteroids with albedo out of this range were classified as interlopers and were therefore disqualified as family members. Sizes were calculated for the asteroids with albedo within the acceptable range. For the other asteroids (those chosen by means of the HCM, but with albedo not listed in the WISE), the value of albedo of the largest member of the family was adopted. Results are given in a set of figures showing the families on the planes (a, e), (a, i), (e, i). Diameters and volumes of the asteroids that are the individual members of a family were calculated on the basis of their known or assumed albedo and on their absolute magnitude. Volumes of the parent bodies of the families were found on the basis of the cumulative volume distribution of these families. We also studied the secular resonances of the family members. We have shown that the locations of members of the considered asteroid families are related to the lines of secular resonances z₁, z₂, and z₃ with Saturn.     

The geometry of the Roche coordinates

The exact geometry of the Roche curvilinear coordinates (, , ) in which corresponds to the zero-velocity surfaces is investigated numerically in the plane, as well as in the spatial, case for various values of the mass-ratio between the two point-masses (m ₁,m ₂) constituting a binary system.The geometry of zero-velocity surfaces specified by -values at the Lagrangian points are first discussed by taking their intersections with various planes parallel to thexy-, xz- andyz-planes. The intersection of the zero-velocity surface specified by the -value at the Lagrangian equilateral-triangle pointsL _4,5 with the planex=1/2 discloses two invariable curves passing through the pointsL _4,5 and situated symmetrically with respect to thexy-plane whose form is independent of the mass-ratio.The geometry of the remaining two coordinates (, ) orthogonal to the zero-velocity surfaces is investigated in thexy- andxz-planes from extensive numerical integrations of differential equations generated from the orthogonality relations among the coordinates. The curves (x, y)=constant in thexy-plane are found to be separated into three families by definite envelopes acting as boundaries whose forms depend upon the mass-ratio only: the inner -constant curves associated with the masspointm ₁, the inner -constant curves associated with the mass-pointm ₂ and the outer -constant curves. All the -constant curves in thexy-plane coalesce at either of the Lagrangian equilateraltriangle pointsL _4,5, except for a limiting case coincident with thex-axis. The curves (x, z)=constant in thexz-plane are also separated by definite envelopes depending upon the mass-ratio into different families: the inner -constant curves associated with the mass-pointm ₁, the inner -constant curves associated with the mass-pointm ₂ and the outer -constant curves on both sides out of the envelopes. For larger values ofz, the curves =constant tend asymptotically to the line perpendicular to thex-axis and passing through the centre of mass of the system, except for a limiting case coincident with thex-axis. The geometrical aspects of the envelopes for the curves (x, y)=constant in thexy-plane and the curves (x, z)=constant in thexz-plane are also discussed independently.In the three-dimensional space, the Roche coordinates can be conveniently defined in such a way as to correspond to the polar coordinates in the immediate neighbourhood of the origin, and to the cylindrical coordinates at great distances. From numerical integrations of simultaneous differential equations generating spatial curves orthogonal to the zero-velocity surfaces, the surfaces (x, y, z)=constant and the surfaces (x, y, z)=constant are constructed as groups of such spatial curves with common values of some parameters specifying the respective surfaces.On leave of absence from the University of Tokyo as an Honorary Fellow of the Victoria University of Manchester.     

Filling factors and scale heights of the diffuse ionized gas in the Milky Way

The combination of dispersion measures of pulsars, distances from the model of Cordes & Lazio (2002) and emission measures from the WHAM survey enabled a statistical study of electron densities and filling factors of the diffuse ionized gas (DIG) in the Milky Way. The emission measures were corrected for absorption and contributions from beyond the pulsar distance. For a sample of 157 pulsars at |b | > 5. and 60° < ℓ < 360°, located in mainly interarm regions within about 3 kpc from the Sun, we find that: (1) The average volume filling factor along the line of sight and the mean density in ionized clouds are inversely correlated: ( ) = (0.0184 ± 0.0011) ^{–1.07 ± 0.03} for the ranges 0.03 < < 2 cm^–3 and 0.8 > > 0.01. This relationship is very tight. The inverse correlation of and causes the well‐known constancy of the average electron density along the line of sight. As (z ) increases with distance from the Galactic plane |z |, the average size of the ionized clouds increases with |z |. (2) For |z| < 0.9 kpc the local density in clouds n _c(z ) and local filling factor f (z ) are inversely correlated because the local electron density n _e(z ) = f (z )n _c(z ) is constant. We suggest that f (z ) reaches a maximum value of >0.3 near |z | = 0.9 kpc, whereas n _c(z ) continues to decrease to higher |z |, thus causing the observed flattening in the distribution of dispersion measures perpendicular to the Galactic plane above this height. (3) For |z | < 0.9 kpc the local distributions n _c(z ), f (z ) and (z ) have the same scale height which is in the range 250 < h ≲ 500 pc. (4) The average degree of ionization of the warm atomic gas (z ) increases towards higher |z | similarly to (z ). Towards |z | = 1 kpc, (z ) = 0.24 ± 0.05 and (z ) = 0.24 ± 0.02. Near |z | = 1 kpc most of the warm, atomic hydrogen is ionized. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Distribution of gaps in emission line redshifts of QSOs

A gap in the distribution of a parameter is simply the absence of the parameter for the values corresponding to the gap. The gap in the emission line redshift (z) of QSOs thus represents absence of QSOs with emission line redshift values corresponding to the gap region. Gaps in emission line redshifts of QSOs have been analysed statistically with updated data consisting of 1549 values. The study indicates: (i) There is a critical redshiftz _c =2.4, which separates two distinct phases in the creation of QSOs. Forz>z _c , the creation appears to have been a slow process. Atz?z _c there was a triggering action which produced a burst of QSOs simultaneously. Forz _c , the rate of production of QSOs have been fast. (ii) The distribution of gaps atz _c ; appear to be consequence of periodicities, provided the periodicities involved are perfect and the redshift values are accurate. (iii) The distribution of gaps atz>z _c are not random, but follow a definite trend.