General relativistic dynamics applied to the rotation curves of galaxies

We extend our general relativistic analysis of galactic rotation curves (see Cooperstock and Tieu in Int. J. Mod. Phys. A 22:2293, 2007) with galaxies NGC 2841, NGC 2903 and NGC 5033. As before, we employ the solution of the Einstein field equations of general relativity with an expansion in Bessel functions. As in our earlier studies, the fits to the data are found to be very precise and the calculated baryonic masses are lower than those based upon Newtonian gravity. Also as in our previous studies, the galactic radii at which the optical luminosities terminate are seen to correlate with densities near 10^−21.75 kg m⁻³. This concordance lends further support to the correctness of the procedure as well as providing a potentially valuable piece of information in the understanding of galactic evolution.     

Equilibrium structure of stars obeying a generalized differential rotation law

In the present paper we have considered the problem of determining the equilibrium structure of differentially rotating stars in which the angular velocity of rotation varies both along the axis of rotation and in directions perpendicular to it. For this purpose, a generalized law of differential rotation of the type ² =b ₀+b ₁ s ²+b ₂ s ⁴+b ₃ z ²+b ₄ z ⁴+b ₅ z ² s ² (here is a nondimensional measure of the angular velocity of a fluid element distants from the axis of rotation andz from the plane through the centre of the star perpendicular to the axis of rotation, andb's are suitably chosen parameters) has been used. Whereas Kippenhahn and Thomas averaging approach has been used to incorporate the rotational effects in the stellar structure equations, Kopal's results on Roche equipotentials have been used to obtain the explicit form of the stellar structure equations, which incorporate the rotational effects up to second order of smallness in the distortion parameters. The method has been used to compute the equilibrium structure of certain differentially rotating polytropes. Certain differentially rotating polytropes. Certain differentially rotating models of the Sun have also been computed by using this approach.     

The dynamics of S0 galaxies and their Tully–Fisher relation

This paper investigates the detailed dynamical properties of a relatively homogeneous sample of disc-dominated S0 galaxies, with a view to understanding their formation, evolution and structure. By using high signal-to-noise ratio long-slit spectra of edge-on systems, we have been able to reconstruct the complete line-of-sight velocity distributions of stars along the major axes of the galaxies. From these data, we have derived both model distribution functions (the phase density of their stars) and the approximate form of their gravitational potentials.
The derived distribution functions are all consistent with these galaxies being simple disc systems, with no evidence for a complex formation history. Essentially no correlation is found between the characteristic mass scalelengths and the photometric scalelengths in these galaxies, suggesting that they are dark-matter dominated even in their inner parts. Similarly, no correlation is found between the mass scalelengths and asymptotic rotation speed, implying a wide range of dark matter halo properties.
By comparing their asymptotic rotation speeds with their absolute magnitudes, we find that these S0 galaxies are systematically offset from the Tully–Fisher relation for later-type galaxies. The offset in luminosity is what one would expect if star formation had been suddenly switched off a few Gyr ago, consistent with a simple picture in which these S0s were created from ordinary later-type spirals which were stripped of their star-forming interstellar medium when they encountered a dense cluster environment.     

Modeling the Newtonian dynamics for rotation curve analysis of thin-disk galaxies

We present an efficient,robust computational method for modeling the Newtonian dynamics for rotation curve analysis of thin-disk galaxies.With appropriate mathematical treatments,the apparent numerical difficulties associated with singularities in computing elliptic integrals are completely removed.Using a boundary element discretization procedure,the governing equations are transformed into a linear algebra matrix equation that can be solved by straightforward Gauss elimination in one step without further ...     

The theory of canonical perturbations applied to attitude dynamics and to the Earth rotation. Osculating and nonosculating Andoyer variables

In the method of variation of parameters we express the Cartesian coordinates or the Euler angles as functions of the time and six constants. If, under disturbance, we endow the “constants” with time dependence, the perturbed orbital or angular velocity will consist of a partial time derivative and a convective term that includes time derivatives of the “constants”. The Lagrange constraint, often imposed for convenience, nullifies the convective term and thereby guarantees that the functional dependence of the velocity on the time and “constants” stays unaltered under disturbance. “Constants” satisfying this constraint are called osculating elements. Otherwise, they are simply termed orbital or rotational elements. When the equations for the elements are required to be canonical, it is normally the Delaunay variables that are chosen to be the orbital elements, and it is the Andoyer variables that are typically chosen to play the role of rotational elements. (Since some of the Andoyer elements are time-dependent even in the unperturbed setting, the role of “constants” is actually played by their initial values.) The Delaunay and Andoyer sets of variables share a subtle peculiarity: under certain circumstances the standard equations render the elements nonosculating. In the theory of orbits, the planetary equations yield nonosculating elements when perturbations depend on velocities. To keep the elements osculating, the equations must be amended with extra terms that are not parts of the disturbing function [Efroimsky, M., Goldreich, P.: J. Math. Phys. 44, 5958–5977 (2003); Astron. Astrophys. 415, 1187–1199 (2004); Efroimsky, M.: Celest. Mech. Dyn. Astron. 91, 75–108 (2005); Ann. New York Acad. Sci. 1065, 346–374 (2006)]. It complicates both the Lagrange- and Delaunay-type planetary equations and makes the Delaunay equations noncanonical. In attitude dynamics, whenever a perturbation depends upon the angular velocity (like a switch to a noninertial frame), a mere amendment of the Hamiltonian makes the equations yield nonosculating Andoyer elements. To make them osculating, extra terms should be added to the equations (but then the equations will no longer be canonical). Calculations in nonosculating variables are mathematically valid, but their physical interpretation is not easy. Nonosculating orbital elements parameterise instantaneous conics not tangent to the orbit. (A nonosculating i may differ much from the real inclination of the orbit, given by the osculating i.) Nonosculating Andoyer elements correctly describe perturbed attitude, but their interconnection with the angular velocity is a nontrivial issue. The Kinoshita–Souchay theory tacitly employs nonosculating Andoyer elements. For this reason, even though the elements are introduced in a precessing frame, they nevertheless return the inertial velocity, not the velocity relative to the precessing frame. To amend the Kinoshita–Souchay theory, we derive the precessing-frame-related directional angles of the angular velocity relative to the precessing frame. The loss of osculation should not necessarily be considered a flaw of the Kinoshita–Souchay theory, because in some situations it is the inertial, not the relative, angular velocity that is measurable [Schreiber, K. U. et al.: J. Geophys. Res. 109, B06405 (2004); Petrov, L.: Astron. Astrophys. 467, 359–369 (2007)]. Under these circumstances, the Kinoshita–Souchay formulae for the angular velocity should be employed (as long as they are rightly identified as the formulae for the inertial angular velocity).     

A new scalar tensor theory for gravity and the flat rotation curves of spiral galaxies

In this paper we reproduce flat rotation curves of spiral galaxies by discussing a scalar tensor theory of gravity which includes the Higgs field as scalar field. The galaxy density distribution is assumed to be homogeneous or polytropic and in the galaxy core there exists a large concentrated mass.     

Cosmological constant and the flat rotation curves of galaxies

It is shown that Einstein's cosmological constant can explain of the flat rotation curves of spiral galaxies.     

A new approach to a relativistic model of galaxies

The observation that stars in galaxies have constant period for a large region from the centre is used to place a constraint on the form of the line element in that region.     

Magnetic fields and the dynamics of spiral galaxies

We investigate the dynamics of magnetic fields in spiral galaxies by performing 3D magnetohydrodynamics simulations of galactic discs subject to a spiral potential using cold gas, warm gas and a two-phase mixture of both. Recent hydrodynamic simulations have demonstrated the formation of interarm spurs as well as spiral arm molecular clouds, provided the interstellar medium model includes a cold H  i phase. We find that the main effect of adding a magnetic field to these calculations is to inhibit the formation of structure in the disc. However, provided a cold phase is included, spurs and spiral arm clumps are still present if β≳ 0.1 in the cold gas. A caveat to the two-phase calculations though is that by assuming a uniform initial distribution, β≳ 10 in the warm gas, emphasizing that models with more consistent initial conditions and thermodynamics are required. Our simulations with only warm gas do not show such structure, irrespective of the magnetic field strength.
Furthermore, we find that the introduction of a cold H  i phase naturally produces the observed degree of disorder in the magnetic field, which is again absent from simulations using only warm gas. Whilst the global magnetic field follows the large-scale gas flow, the magnetic field also contains a substantial random component that is produced by the velocity dispersion induced in the cold gas during the passage through a spiral shock. Without any cold gas, the magnetic field in the warm phase remains relatively well ordered apart from becoming compressed in the spiral shocks. Our results provide a natural explanation for the observed high proportions of disordered magnetic field in spiral galaxies and we thus predict that the relative strengths of the random and ordered components of the magnetic field observed in spiral galaxies will depend on the dynamics of spiral shocks.     

Flat rotation curves of spiral galaxies and the dark matter particles

It is shown that a gas of massive bosonic particles ( eV), e.g. Higgs particles, surrounding the disk galaxies is able to generate flat rotation curves without internal difficulties with physical principles.     

Dark matter particles and the flat rotation curves of spiral galaxies

It is shown that a gas of massive bosonic particles (m 60 eV), e.g. Higgs particles, surrounding disk galaxies is able to generate the flat rotation curves of spiral galaxies without internal difficulties with physical principles. We have analyzed 36 galaxies and find good agreement with the empirical data.     

cD galaxies of apparent supergiant sizes due to the curvature of space

By combining two two-dimensional subspaces, closed into themselves due to curvature, it is possible to create a model of three-dimensional space of the same properties. If the Universe is a space of this type, its effect is that of a monstrous lens. Close objects are observed to diminish according to the normal law of perspective; however, the remote galaxies are seen to be very highly magnified.The apparent angular size² of a galaxy is more than the size¹ in flat space according to relation:² =¹ cosec , where is the angular distance from the observer to the galaxy. The diameter² d of a galaxy in curved space must be in the same relation to a diameter¹ d with no curvature of space:² d=¹ d cosec . The apparent angular size² and diameter² d are distorted shapes in consequence of an optical illusion caused by the spatial curvature.It is necessary to distribute the multitude of galaxies into two parts in accordance with their location on the close or reverse hemihypersphere of the Universe. The minimum of apparent angular size² of a galaxy of diameter¹ d is at the equatorial zone.The most likely candidates for location in the reverse hemi-hypersphere are cD's of apparent supergiant sizes due, probably, to the curvature of space. The existence of supergiant sizes of galaxies is the second indirect proof, besides superluminal velocities, that the Universe is closed into itself through curvature. The third indirect evidence, i.e., inductive confirmation of the same fact, is the superposition of galaxies which need not inevitably be a new alternative to the present theories of collisions, cannibalism, merger, etc.The fourth indirect proof of the positive curvature of the Universe is the occurrence of background radiation, because that must vanish in hyperbolic space irrespective of its origin. The gravitational lens effect acquires another theoretical form, as usual, in the case of remote galaxies, because it is impossible to distinguish between gravitator and lensing image.     

Coronal expansion and the transfer of solar angular momentum to the interplanetary space

We discuss the question of loss of angular momentum through coronal expansion. From a large volume of data on Type-1 cometary tails we have confirmed the presence of a tangential component in the coronal expansion, which has not only a stochastic component but also a constant component of 9.8 km/s. Through coronal expansion the Sun has lost 80% of its angular momentum since it evolved on to the main sequence and the angular velocity of the Sun is decreasing exponentially. This result should have a large effect on the dynamical evolution of the Sun.     

Model rotation curves of disk galaxies oriented at an arbitrary angle to the plane of the sky

We consider the effects of projection, absorption, and velocity dispersion on the shape of the rotation curve for a galaxy as a function of its disk inclination to the plane of the sky. We conclude that for galaxies with a fairly massive compact bulge these effects lead to a marked discrepancy between the rotation curve and the circular velocity curve even if their disks are viewed far from edge-on, especially for the rotation curves constructed from observations of the stellar component.     

Application of a newN-body code to the encounters of galaxies and dissipationless formation of elliptical galaxies

A new code has been used to simulate the formation of elliptical galaxies via a dissipationless collapse by a numerical experiment. Preliminary results of numerical simulations are presented for initial 2T/|W| conditions ranging from 10^–2 to 10^–5. The importance of inhomogeneity in the initial density distribution and of merging is outlined.Paper presented at the 11th European Regional Astronomical Meetings of the IAU on New Windows to the Universe, held 3–8 July, 1989, Tenerife, Canary Islands, Spain.     

Using space manifold dynamics to deploy a small satellite constellation around the Moon

The possibility of communicating with the far side of the Moon is essential for keeping a continuous radio link with lunar orbiting spacecraft, as well as with manned or unmanned surface facilities in locations characterized by poor coverage from Earth. If the exploration and the exploitation of the Moon must be sustainable in the medium/long term, we need to develop the capability to realize and service such facilities at an affordable cost. Minimizing the spacecraft mass and the number of launches is a driving parameter to this end. The aim of this study is to show how Space Manifold Dynamics can be profitably applied in order to launch three small spacecraft onboard the same launch vehicle and send them to different orbits around the Moon with no significant difference in the Delta-V budgets. Internal manifold transfers are considered to minimize also the transfer time. The approach used is the following: we used the linearized solution of the equations of motion in the Circular Restricted Three Body Problem to determine a first–guess state vector associated with the Weak Stability Boundary regions, either around the collinear Lagrangian point L1 or around the Moon. The resulting vector is then used as initial state in a numerical backward-integration sequence that outputs a trajectory on a manifold. The dynamical model used in the numerical integration is four-body and non-circular, i.e. the perturbations of the Sun and the lunar orbital eccentricity are accounted for. The trajectory found in this way is used as the principal segment of the lunar transfer. After separation, with minor maneuvers each satellite is injected into different orbits that lead to ballistic capture around the Moon. Finally, one or more circularization maneuvers are needed in order to achieve the final circular orbits. The whole mission profile, from launch to insertion into the final lunar orbits, is modeled numerically with the commercial software STK.     

Solution to the rotation of the elastic earth by method of rigid dynamics

Hamiltonian mechanics is applied to the problem of the rotation of the elastic Earth. We first show the process for the formulation of the Hamiltonian for rotation of a deformable body and the derivation of the equations of motion from it. Then, based on a simple model of deformation, the solution is given for the period of Euler motion, UT1 and the nutation of the elastic Earth. In particular it is shown that the elasticity of the Earth acts on the nutation so as to decrease the Oppolzer terms of the nutation of the rigid Earth by about 30 per cent. The solution is in good agreement with results which have been obtained by other, different approaches.