A boundary equation for uniformly rotating polytropes

A linear approximating equation exists for the boundary of a uniformly rotating polytrope. The equation in η=(ξ₁−ξ)/ξ₁ permits rapid calculation of the polytrope radius for any latitude, and is accurate for angular velocities of rotation nearly to critical rotation. Data in this paper apply to a polytrope indexn=3.0.     

Rotational distortion of polytropes by analytic approximation

This paper modifies the first-order perturbation theory of Chandrasekhar, for rotational distortion of polytropes. Comparison with numerical integrations by other authors demonstrates that the present analytic theory is as accurate as other published first-order theories. The present theory is in a form permitting rapid calculation of boundary shapes as a function of the rotation parameter,v, and the polytrope index,n. Results are presented for a critically rotating polytrope, for the casen=3.     

A second-order perturbation theory for differentially rotating gaseous polytropes

The aim of the present investigation has been to specify the structure of a differentially rotating gaseous polytrope, by extending Chandrasekhar's method to include second-order terms in the perturbation parameter. The formal results obtained permit the determination of the polytropic structure for all types of differential rotation of cylindrical symmetry. The formalism developed works efficiently in the limiting case of uniform rotation and the results obtained are compared with those of previous investigations of uniformly rotating polytropes.     

Possible solution to the problem of the extreme mass ratio W UMa-type binaries

When the total angular momentum of a binary system is at a critical (minimum) value, a tidal instability occurs (Darwin's instability), eventually forcing the stars to merge into a single, rapidly rotating object. The instability sets in at some critical separation which in the case of contact binaries corresponds to a minimum mass ratio depending on dimensionless gyration radius k ₁. If one considers   n = 3  polytrope (fully radiative primary with  Γ₁= 4/3  ),   k ²₁= 0.075  and   q _min≈ 0.085–0.095  . There appears to be, however, some W UMa-type binaries with q values very close, if not below these theoretical limits, implying that primary in these systems is probably more centrally condensed. We try to solve the discrepancy between theory and observations by considering rotating polytropes. We show by deriving and solving a modified Lane–Emden equation for   n = 3  polytrope that including the effects of rotation does increase the central concentration and could reduce   q _min  to as low as 0.070–0.074, more consistent with the observed population.     

The general solution of the HENON–HEILES problem

The general solution of the Henon–Heiles system is approximated inside a domain of the (x, C) of initial conditions (C is the energy constant). The method applied is that described by Poincaré as ‘the only “crack” permitting penetration into the non-integrable problems’ and involves calculation of a dense set of families of periodic solutions that covers the solution space of the problem. In the case of the Henon–Heiles potential we calculated the families of periodic solutions that re-enter after 1–108 oscillations. The density of the set of such families is defined by a pre-assigned parameter ε (Poincaré parameter), which ascertains that at least one periodic solution is computed and available within a distance ε from any point of the domain (x, C) for which the approximate general solution computed. The approximate general solution presented here corresponds to ε = 0.07. The same solution is further improved by “zooming” into four square sub-domain of (x, C), i.e. by computing sufficient number of families that reduce the density parameter to ε = 0.003. Further zooming to reduce the density parameter, say to ε = 10⁻⁶, or even smaller, although easily performable in both areas occupied by stable as well as unstable solutions, was found unnecessary. The stability of all members of each and all families computed was calculated and presented in this paper for both the large solution domain and for the sub-domains. The correspondence between areas of the approximate general solution occupied by stable periodic solutions and Poincaré sections with well-aligned section points and also correspondence between areas occupied by unstable solutions and Poincaré sections with randomly scattered section points is shown by calculating such sections. All calculations were performed using the Runge-Kutta (R-K) 8th order direct integration method and the large output received, consisting of many thousands of families is saved as “Atlas of the General Solution of the Henon–Heiles Problem,” including their stability and is available at request. It is concluded that approximation of the general solution of this system is straightforward and that the chaotic character of its Poincaré sections imposes no limitations or difficulties.     

Numerical exploration of the photogravitational restricted five-body problem

We study numerically the restricted five-body problem when some or all the primary bodies are sources of radiation. The allowed regions of motion as determined by the zero-velocity surface and corresponding equipotential curves, as well as the positions of the equilibrium points are given. We found that the number of the collinear equilibrium points of the problem depends on the mass parameter β and the radiation factors q _i , i=0,…,3. The stability of the equilibrium points are also studied. Critical masses associated with the number of the equilibrium points and their stability are given. The network of the families of simple symmetric periodic orbits, vertical critical periodic solutions and the corresponding bifurcation three-dimensional families when the mass parameter β and the radiation factors q _i vary are illustrated. Series, with respect to the mass (and to the radiation) parameter, of critical periodic orbits are calculated.     

New parametric classes of exact solutions in general relativity and polytropic nature fluid ball with constant sound speed

In this paper we have presented a method of obtaining varieties of new parametric classes of spherically symmetric analytic solutions of the general relativistic field equations in canonical coordinates. A number of previously known classes of solutions has been rediscovered which describe perfect fluid balls with infinite central pressure and infinite central density though their ratio is positively finite and less then one. From the solutions of one of the class we have constructed a causal model of polytrope with constant sound speed Corresponding to the polytrope model we have maximized the Neutron star mass 3.26 M_⊙ with the linear dimensions 32.27 kms with surface red shift 0.7355 and for other class we have constructed a causal model in which outmarch of pressure and density is monotonically decreasing and pressure–density ratio is positive and less than 1 throughout with in the ball. Corresponding to this model we have maximized the Neutron star mass 3.09 M_⊙ with the linear dimensions 30.55 kms with surface red shifts 0.5811.     

Rotational dynamics of celestial bodies

In this paper the ‘class of near homoaxial rotations’ is defined, being suitable for treatment of problems of nonuniform rotation of stars. This class is represented by a proper form of the so-called ‘velocity tensor’, the latter describing efficiently the motion of a deformable finite material continuum in the common frame. The ‘class of particular near homoaxial rotations’ is then defined, characterized by simple transformation equations of the velocity tensor in two noninertial frames; namely, in a ‘frame rotating uniformly’ relative to the common frame, and in a ‘frame rotating nonuniformly’ relative to it. A sufficient condition is also derived so that a near homoaxial rotation be reducible to a particular one. ‘Preferred frames’ are then defined in the sense that they preserve a near homoaxial rotation in its class when referring thismotion; that is, such frames keep invariant the intertial class of the motion. Finally, a method is proposed for constructing a nonuniformly rotating preferred frame, to which a near homoaxial rotation is referred simply as ‘radial distortion’.     

Effect of perturbations on the non linear stability of triangular points in the restricted three-body problem with variable mass

The effect of small perturbations ε and ε ^′ in the Coriolis and the centrifugal forces, respectively on the nonlinear stability of the triangular points in the restricted three-body problem with variable mass has been studied. It is found that, in the nonlinear sense, the triangular points are stable for all mass ratios in the range of linear stability except for three mass ratios, which depend upon ε, ε ^′ and β, the constant due to the variation in mass governed by Jeans’ law.     

Structure-distortion of polytropic stars by differential rotation: Third-order perturbation theory

In order to specify the structure of a rapidly and differentially rotating gaseous polytrope, we extend Chandrasekhar's perturbation theory to include third-order terms in the perturbation parameter. In the present paper, the theory developed is required for a subsequent numerical treatment of the structure-determination.This research was supported by the Research Development Project of the University of Patras, Greece.     

Evolution of astrophysical jets generated by a vortex mechanism

The nonlinear dynamics of a rotating jet is examined following its ejection from a compact gravitating object by a vortex mechanism. A scenario is described in which a dense stream expands and is subsequently transformed into a nonstationary vortex consisting of a cylindrical core and a “sheath.” At this stage of development, a converging radial flow of matter in the differentially rotating nonuniform sheath collimates the jet and speeds up the rotation of the core, as well as the flow of matter along the jet, in accordance with a power law or “explosive” instability, until the velocity discontinuity at the surface of the core approaches the sound speed. Flows of this type have low energy dissipation and can serve as unique channels for the acceleration and collimation of jet eruptions from young stars, active galactic nuclei, and quasars. Translated from Astrofizika, Vol. 52, No. 1, pp. 135–145 (February 2009).     

Minimal Energy and Gravitation

Quantum theory in Robertson – Walker spacetime suggests the existence of a minimal energy ε of the order of 10⁻⁴⁵ erg. Reasonable forms for ε give the expansion factor R=R(t)(t= the cosmic time) with no need of gravitational field equations.Einstein's theory should be modified in gravitational fields of strength less than ε c/ħ ∼ 10⁻⁸ cm/s² where c is the speed of light and ħ is the reduced Planck constant. The cosmological term λ is expected to decrease as the universe expands.In the Appendix, ε is derived from a big bang – big crunch Newtonian cosmology.     

On a uniformly rotating magnetic polytrope

The equilibrium general magnetic field inside a magneto-rotating star, assumed to be a polytrope, has been determined more accurately, for large general magnetic field. Furthermore the effect of such field on the structure and oscillations of a slowly rotating polytrope has been studied forn=1,0, 1.5, 2.0, and 3.0.     

Bifurcation points in the planar problem of three bodies

An essential parameter in the planar problem of three bodies is the product of the square of the angular momentum and of the total energy (c ² H). The role of this parameter, which may be called abifurcation parameter, in establishing regions of possible motions has been shown by Marchal and Saari (1975) and Zare (1976a). There exist critical values of this parameter below which exchange between bodies cannot occur. These critical values may be calledbifurcation points.This paper gives an analytical criterion to obtain these bifurcation points for any given masses of the participating bodies.     

Turbulent compressible convection with rotation–penetration below a convection zone

We examine the behaviour of penetrative turbulent compressible convection under the influence of rotation by means of three dimensional numerical simulations. We estimate the extent of penetration below a stellar-type rotating convection zone in an f-plane configuration. Several models have been computed with a stable-unstable-stable configuration by varying the rotation rate (Ω), the inclination of the rotation vector and the stability of the lower stable layer. The spatial and temporal average of kinetic energy flux (F_k) is computed for several turnover times after the fluid has thermally relaxed and is used to estimate the amount of penetration below the convectively unstable layer. Our numerical experiments show that with the increase in rotational velocity, the downward penetration decreases. A similar behaviour is observed when the stability of the lower stable layer is increased in a rotating configuration. Furthermore, the relative stability parameter S shows an S ^−1/4 dependence on the penetration distance implying the existence of a thermal adjustment region in the lower stable layer rather than a nearly adiabatic penetration region.     

Rosette Central Configurations, Degenerate Central Configurations and Bifurcations

In this paper we find a class of new degenerate central configurations and bifurcations in the Newtonian n-body problem. In particular we analyze the Rosette central configurations, namely a coplanar configuration where n particles of mass m₁ lie at the vertices of a regular n-gon, n particles of mass m₂ lie at the vertices of another n-gon concentric with the first, but rotated of an angle π /n, and an additional particle of mass m₀ lies at the center of mass of the system. This system admits two mass parameters μ = m₀/m₁ and ε = m₂/m₁. We show that, as μ varies, if n > 3, there is a degenerate central configuration and a bifurcation for every ε > 0, while if n = 3 there is a bifurcation only for some values of ε.     

On the nature of compact galactic nuclei

Spherically symmetric stellar clusters (compact galactic nuclei and globular clusters), far advanced toward the state of complete statistical equilibrium in the course of evolution, are investigated. The equation of state of such systems (a polytrope with an index k = 0.5) is derived and their main characteristics are calculated. It is shown that compact galactic nuclei must consist mainly of rapidly rotating neutron stars and white dwarfs. It is demonstrated that pulsars may be created by the evaporation of neutron stars from the nucleus of our Galaxy. The number of such pulsars is ~3.10⁶. Translated from Astrofizika, Vol. 41, No. 1, pp. 41–50, January-March, 1998.     

Fundamental energy and closed space

The existence of the universal quantization law E=n ε E =any energy; n = an integer, ε = the fundamental energy ∼ħ c/R with ħ = the reduced Planck constant, c = the speed of light, R = the curvature radius of the closed cosmological space) is advocated and discussed. A possible connection between ε and the mass of elementary particles is pointed out.     

Hamiltonian Stability of Spin–Orbit Resonances in Celestial Mechanics

The behaviour of ‘resonances’ in the spin-orbit coupling in celestial mechanics is investigated in a conservative setting. We consider a Hamiltonian nearly-integrable model describing an approximation of the spin-orbit interaction. The continuous system is reduced to a mapping by integrating the equations of motion through a symplectic algorithm. We study numerically the stability of periodic orbits associated to the above mapping by looking at the eigenvalues of the matrix of the linearized map over the full cycle of the periodic orbit. In particular, the value of the trace of the matrix is related to the stability character of the periodic orbit. We denote by ε_* (p/q) the value of the perturbing parameter at which a given elliptic periodic orbit with frequency p/q becomes unstable. A plot of the critical function ε_* (p/q) versus the frequency at different orbital eccentricities shows significant peaks at the synchronous resonance (for low eccentricities) and at the synchronous and 3:2 resonances (at higher eccentricities) in good agreement with astronomical observations. This revised version was published online in July 2006 with corrections to the Cover Date.