Simulations of stellar disks in galaxies with non-Newtonian interactions

We describe gravitationalN-body simulations to investigate whether various non-Newtonian interactions between the stars of a system could explain the flat rotational curves which are characteristic of actual isolated spiral galaxies. It is shown that replacing the standard Newtonian interaction by the models of Sanders (1984), Kuhn and Kruglyak (1987) and Milgrom (1983), no massive halo (or dark matter) is required to produce the flat rotational curves of the systems under consideration. All models also generate the exponential surface mass density distribution which is in agreement with that observed in disk-shaped galaxies. In relation to the spiral structure of galaxies, we present the evidence that the non-Newtonian interactions can reproduce the multiple armed patterns in stellar disks without dark matter.     

Treatment of close approaches in the numerical integration of the gravitational problem ofN bodies

One of the main difficulties encountered in the numerical integration of the gravitationaln-body problem is associated with close approaches. The singularities of the differential equations of motion result in losses of accuracy and in considerable increase in computer time when any of the distances between the participating bodies decreases below a certain value. This value is larger than the distance when tidal effects become important, consequently,numerical problems are encounteredbefore the physical picture is changed. Elimination of these singularities by transformations is known as the process of regularization. This paper discusses such transformations and describes in considerable detail the numerical approaches to more accurate and faster integration. The basic ideas of smoothing and regularization are explained and applications are given.     

Earth–Mars transfers with ballistic escape and low-thrust capture

In this paper novel Earth–Mars transfers are presented. These transfers exploit the natural dynamics of n-body models as well as the high specific impulse typical of low-thrust systems. The Moon-perturbed version of the Sun–Earth problem is introduced to design ballistic escape orbits performing lunar gravity assists. The ballistic capture is designed in the Sun–Mars system where special attainable sets are defined and used to handle the low-thrust control. The complete trajectory is optimized in the full n-body problem which takes into account planets’ orbital inclinations and eccentricities. Accurate, efficient solutions with reasonable flight times are presented and compared with known results.     

Coupled-oscillators

Modal coupling oscillation models for the stellar radial pulsation and coupled-oscillators are reviewed. Coupled-oscillators with the second-order and third-order terms seemed to behave non-systematically. Using the equation by Schwarzschild and Savedoff (1949) with the dissipation term of van del Pol's type which is third-order, we demonstrate the effect of each term. The effects can be understood by the terms of the nonlinear dynamics, which is recently developing, that is. phase-locking, quasi-periodicity, period doubling, and chaos. As the problem of stellar pulsation, especially of double-mode cepheids on the period-ratio, we examine the dependence on the stellar structure from which the coupling constants in the second-order terms are derived. Eigen functions for adiabatic pulsations had been used for the calculation of the constants. It is noted that only two set of the constants are available, that is, for the polytrope model withn = 3 and a cepheid model without convection. Some examples of nonlinear dynamical effects will be shown.It is shown that if the constants were suitable values, the period-ratio of double-mode cepheids is probably realized. The possibility is briefly suggested.     

Wintner's collinear and flat solutions are nowhere dense

It is shown that in then-body problem with generalized attraction law, the sets of initial conditions which lead to Wintner's collinear, respectively flat, motion are nowhere dense relative to the set of initial conditions that define solutions inR ³.     

Galaxies in the connection machine

Many problems in galactic dynamics require computer speeds orders of magnitude larger than what can be achieved on current single-processor machines. In the near future such speeds are likely to become available through computer architectures based on large-scale, fine-grained parallelism. An example of a highly parallel computer is the Connection Machine, with up to 65,636 processors. We have benchmarked gravitationalN-body algorithms on the Connection Machine, and compared those with similar benchmarks which we have obtained on more traditional vector supercomputers. Our conclusions are: (1) The direct summation algorithm, with of orderN ² interactions forN particles, can be made to run with high efficiency on either type of computer. As a result, the Connection Machine clearly wins in speed over all supercomuters tested, with the exception of an 8-processor ETA, which shows a comparable performance. (2) A more efficient tree algorithm reduces the growth of the number of interactions fromN ² toN logN. However, the greater complexity of this algorithm causes a considerable degradation of efficiency, by a factor which is larger on the Connection Machine than on vector supercomputers. As a result, out tree code runs at comparable speeds on both types of machines, with the notable exception of the 8-processor ETA, which has an extrapolated speed for running our tree code which is higher than any of the other machines we have tested.     

Conjugate Points in the Gravitational n-Body Problem

In an effort to understand the nature of almost periodic orbits in the n-body problem (for all time t) we look first to the more basic question of the oscillatory nature of solutions of this problem (on a half-line, usually taken as R ⁺). Intimately related to this is the notion of a conjugate point(due to A. Wintner) of a solution. Specifically, by rewriting the mass unrestricted general problem of n-bodies in a symmetric form we prove that in the gravitational Newtonian n-body problem with collisionless motions there exists arbitrarily large conjugate points in the case of arbitrary (positive) masses whenever the cube of the reciprocal of at least one of the mutual distances is not integrable at infinity. The implication of this result is that there are possibly many Wintner oscillatorysolutions in these cases (some of which may or may not be almost periodic). As a consequence, we obtain sufficient conditions for all continuable solutions (to infinity) to be either unbounded or to allow for near misses (at infinity). The results also apply to potentials other than Newtonian ones. Our techniques are drawn from results in systems oscillation theory and are applicable to more general situations. Dedicated to the memory of Robert M. (Bob) Kauffman, formerly Professor of the University of Alabama in Birmingham     

Dynamics of very rich open clusters

The oldest open clusters in our Galaxy set the lower limit to the age of the Galactic Disk (9–10 Gyr). Although they appear to be very rich now, it is clear that their primordial populations were much larger. Often considered as transitional objects, these populous open clusters show structural differences with respect to globular clusters so their dynamics and characteristic evolutionary time scales can also be different. On the other hand, their large membership lead to different dynamical evolution as compared with average open clusters. In this paper, the differential features of the evolution of rich open clusters are studied using N-body simulations, including several of the largest (10⁴ stars) published direct collisional N-body calculations so far, which were performed on a CRAY YMP. The disruption rate of rich open clusters is analysed in detail and the effect of the initial spatial distribution of the stars in the cluster on its dynamics is studied. The results show that cluster life-time depends on this initial distribution, decreasing when it is more concentrated. The effect of stellar evolution on the dynamical evolution of rich clusters is an important subject that also has been considered here. We demonstrate that the cluster's life-expectancy against evaporation increases because of mass loss by evolving high-mass stars. This revised version was published online in July 2006 with corrections to the Cover Date.     

Direct N-body modelling of stellar populations: blue stragglers in M67

We present a state-of-the-art N -body code which includes a detailed treatment of stellar and binary evolution as well as the cluster dynamics. This code is ideal for investigating all aspects relating to the evolution of star clusters and their stellar populations. It is applicable to open and globular clusters of any age. We use the N -body code to model the blue straggler population of the old open cluster M67. Preliminary calculations with our binary population synthesis code show that binary evolution alone cannot explain the observed numbers or properties of the blue stragglers. On the other hand, our N -body model of M67 generates the required number of blue stragglers and provides formation paths for all the various types found in M67. This demonstrates the effectiveness of the cluster environment in modifying the nature of the stars it contains, and highlights the importance of combining dynamics with stellar evolution. We also perform a series of N =10 000 simulations in order to quantify the rate of escape of stars from a cluster subject to the Galactic tidal field.     

On the number of central configurations in the N-body problem

Central configurations are critical points of the potential function of the n-body problem restricted to the topological sphere where the moment of inertia is equal to constant. For a given set of positive masses m ₁,..., m _n we denote by N(m ₁, ..., m _n, k) the number of central configurations' of the n-body problem in ^k modulus dilatations and rotations. If m _{n 1},..., m _n, k) is finite, then we give a bound of N(m ₁,..., m _n, k) which only depends of n and k.     

How the Method of Minimization of Action Avoids Singularities

The method of minimization of action is a powerful technique of proving the existence of particular and interesting solutions of the n-body problem, but it suffers from the possible interference of singularities. The minimization of action is an optimization and, after a short presentation of a few optimization theories, our analysis of interference of singularities will show that:(A) An n-body solution minimizing the action between given boundary conditions has no discontinuity: all n-bodies have a continuous and bounded motion and thus all eventual singularities are collisions;(B) A beautiful extension of Lambert's theorem shows that, for these minimizing solutions, no double collision can occur at an intermediate time;(C) The proof can be extended to triple and to multiple collisions. Thus, the method of minimization of action leads to pure n-body motions without singularity at any intermediate time, even if one or several collisions are imposed at initial and/or final times.This method is suitable for non-infinitesimal masses only. Fortunately, a similar method, with the same general property with respect to the singularities, can be extended to n-body problems including infinitesimal masses.     

Improbability of collinear solutions in the n-body problem with generalized attraction law

It is shown that in the n-body problem with generalized attraction law (inverse (α + l )-power of the distance, α > 0) the set of initial conditions which lead to collinear motion is of Lebesgue measure zero and nowhere dense with respect to the set of initial conditions that define solutions in R³ or R².     

Structural stability of homothetic solutions of the collinearn-body problem

We investigate specific homothetic solutions of then-body problem which both begin and end in a simultaneous collision of all of the particles. Under a suitable change of variables, these solutions become heteroclinic orbits, i.e., they lie in the intersection of the stable and unstable manifolds of distinct equilibrium points. Our main result is that these manifolds intersect transversely along these orbits. This proves that the homothetic solutions are structurally stable.Partially supported by NSF Grant MCS 77-00430.     

Periodic solutions of circular-elliptic type in the planarN-body problem

Letn2 mass points with arbitrary masses move circularly on a rotating straight-line central-configuration; i.e. on a particular solution of relative equilibrium of then-body problem. Replacing one of the mass points by a close pair of mass points (with mass conservation) we show that the resultingN-body problem (N=n+1) has solutions, which are periodic in a rotating coordinate system and describe precessing nearlyelliptic motion of the binary and nearlycircular collinear motion of its center of mass and the other bodies; assuming that also the mass ratio of the binary is small.     

On the variational approach to the periodic n-body problem

This expository paper gathers some of the results obtained by the author in recent works in collaboration with Davide Ferrario and Vivina Barutello, focusing on the periodic n-body problem from the perspective of the calculus of variations and minimax theory. These researches were aimed at developing a systematic variational approach to the equivariant periodic n-body problem in the two and three-dimensional space. The purpose of this paper is to expose the main problems and achievements of this approach. The material here was exposed in the talk that given at the Meeting CELMEC IV promoted by SIMCA (Società italiana di Meccanica Celeste).     

Photoionization of cool MHD disk winds

Recent ultraviolet observations point out that there is hot, dense plasma associated with the optical jet in some T Tauri stars. In this contribution, cool MHD disk wind physics is reviewed by means of a self-similar analytical model to analyze whether hot (Te ? 80,000 K) and dense (n_e ? 10⁹ cm^-3) plasma can be produced in disk winds. It is shown that these high densities can only be achieved at the base of the wind where the stellar X-rays radiation field is strong. The propagation of the X-rays radiation through the disk wind is analyzed: a cocoon of photoionized gas is generated around the star. However, it is difficult to foresee how temperatures as high as ～ 5 × 10⁴ can be reached unless a significant fraction of the X-rays radiation is produced by magnetic reconnection at the boundary between the stellar magnetosphere and the accretion disk.     

New stacked central configurations for the planar 5-body problem

A stacked central configuration in the n-body problem is one that has a proper subset of the n-bodies forming a central configuration. In this paper we study the case where three bodies with masses m ₁, m ₂, m ₃ (bodies 1, 2, 3) form an equilateral central configuration, and the other two with masses m ₄, m ₅ are symmetric with respect to the mediatrix of the segment joining 1 and 2, and they are above the triangle generated by {1, 2, 3}. We show the existence and non-existence of this kind of stacked central configurations for the planar 5-body problem.     

Invariant sets and polhodes in the rigid body problem

In this work we will describe the sets in the rigid body phase space where the energy and angular momentum are constant, and it will turn out that in nontrivial cases they will simply take the form of cartesian products of the polhodes byS ¹. These sets are important for the global study of said geodesic mechanical system for being invariant under Euler's equations (energy and momentum are constant along their solutions).To motivate from something more familiar in celestial mechanics, we will begin to relate the problem to Smale's study of the planarn-body problem (Smale, 1970) and Easton's study of the planar 3-body problem (Easton, 1971), exemplifying in particular with the central force problem.In the last Sections 4 and 5, we extent our methods to give results for generalized solids on Lie groups, mentioning the further extensions to transitive mechanical systems.Proceedings of the Sixth Conference on Mathematical Methods in Celestial Mechanics held at Oberwolfach (West Germany) from 14 to 19 August, 1978.This work was partially supported by the Consejo Nacional de Ciencia y Tecnología (México) under grant PNCB-049.     

On the instability of weakly radially anisotropic star clusters

We discuss contradictions existing in the literature in the problem on the stability of collisionless spherical stellar systems, which are the simplest anisotropic generalization of the well-known polytropic models. On the one hand, calculations of the growth rates within the framework of a linear stability theory and N-body simulations suggest that these systems should become stable when the parameter s characterizing the degree of anisotropy of the stellar velocity distribution becomes lower than some critical value s _crit > 0. On the other hand, according to Palmer and Papaloizou, the growth rate should be nonzero up to the isotropic limit s = 0. Using our method of determining the eigenmodes of stellar systems, we show that even though the mode growth rates in weakly radially anisotropic systems of this type are nonzero, they are exponentially small, i.e., decrease as γ ∝ exp(−a/s) when s → 0. For slightly radially anisotropic systems with a finite lifetime, this actually implies stability.     

Partial Reduction in the N-Body Planetary Problem using the Angular Momentum Integral

We present a new set of variables for the reduction of the planetary n-body problem, associated to the angular momentum integral, which can be of any use for perturbation theory. The construction of these variables is performed in two steps. A first reduction, called partial is based only on the fixed direction of the angular momentum. The reduction can then be completed using the norm of the angular momentum. In fact, the partial reduction presents many advantages. In particular, we keep some symmetries in the equations of motion (d'Alembert relations). Moreover, in the reduced secular system, we can construct a Birkhoff normal form at any order. Finally, the topology of this problem remains the same as for the non-reduced system, contrarily to Jacobi's reduction where a singularity is present for zero inclinations. For three bodies, these reductions can be done in a very simple way in Poincaré's rectangular variables. In the general n-body case, the reduction can be performed up to a fixed degree in eccentricities and inclinations, using computer algebra expansions. As an example, we provide the truncated expressions for the change of variable in the 4-body case, obtained using the computer algebra system TRIP.