Next Article in Journal
Mixtures of Semi-Parametric Generalised Linear Models
Previous Article in Journal
Propagation Prediction of Body Waves in Fluid-Saturated Soils with Flow-Independent Viscosity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Simplification of Galactic Dynamic Equations

School of Mathematical Science, Fudan University, Shanghai 200433, China
Symmetry 2022, 14(2), 407; https://doi.org/10.3390/sym14020407
Submission received: 10 January 2022 / Revised: 5 February 2022 / Accepted: 16 February 2022 / Published: 18 February 2022
(This article belongs to the Section Physics)

Abstract

:
Galactic dynamics is the foundation for simulating galactic structure and for solving other problems. However, the traditional dynamic equations include some unreasonable assumptions and are therefore scientifically invalid. In this paper, by introducing the following three working assumptions, we established the galactic dynamics of high precision and convenient formalism. 1. In the research of large-scale structure, the retarded potential of the gravitational field should be taken into account, and the weak field and low velocity approximation of Einstein’s field equation should be adopted. 2. The stars in a fully developed galaxy should be zero-pressure and inviscid fluid, and the equation of motion is different from that of ordinary continuum mechanics. Stars move along geodesics. 3. The structure of the galaxy is only related to the total mass density distribution. The equation of state of dark halo is different from that of ordinary luminous interstellar matter, so their trajectories are also very different. In a galaxy, the dark halo and the ordinary matter are automatically separated. The total mass density distribution can be presupposed according to the observation data, and then it can be determined by comparing the solution of the equations with the observed data. These assumptions and treatments are supported by theory and observation. The variables of the equations of simplified galactic dynamics are separated from each other, and the equations are well-posed and can be solved according to a definite procedure. The solution explains the Tully–Fisher relation. Therefore, this simplified dynamic equation system provides a more reasonable and practical framework for the further study of galactic structure, and can solve many practical problems. In addition, it is closely related to the study of dark matter halo in galaxy.

1. Introduction

Most fully developed galaxies have vivid spiral structures that have a profound effect on the generation and evolution of stars. In the past 70 years, a great deal of research has been done to reveal the essence of spiral structures. Although great progress has been made in understanding the mechanism of spiral arm formation, there are still many mysteries to be revealed. We do not even know whether the spiral is a long-standing phenomenon, or whether it is a transient and multiple regeneration phenomenon in the evolution of galaxies [1,2,3]. The prevailing view now is that the galactic spiral structure should be a quasi-stationary density wave [4]. According to the concept of density wave proposed by B. Lindblad, C. C. Lin and F. H. Shu analyzed and solved the quasi-stable solution of planetary fluid dynamics in gravitational field. However, there are other ideas that the spiral structure is a short-lived but recurring spiral pattern formed by the vibration amplification in the rotating galaxy disk [5], or caused by external forces [6]. In the literature [7,8], by comparing the distant galaxies observed by the Hubble Space Telescope with those observed on the ground, the results suggest that the spiral structure in the range of 0 < z < 1 in the universe may be a long-standing quasi-stationary pattern.
In density wave theory, Lin and Shu assume that the spiral structure is quasi-stationary in the form of density waves on the disk of a galaxy. The key points of this assumption are as follows: 1. The spiral galaxy has an invisible background gravitational field, which has a regular spiral structure, and the visible optical spiral arm is only the external manifestation of this spiral force field. 2. The galactic spiral pattern is made up of flowing materials, which is not the fixed same materials on the spiral arm, but changes over time. 3. The entire spiral pattern is quasi-steady and rotates at constant angular velocity Ω like a rigid body. The above assumption is called the quasi-stable spiral structure assumption.
A spiral galaxy consists of at least three parts: disk, bulge and halo. Because of the large dispersion velocity of the two parts of the bulge and halo, it is not easy to produce density disturbance, so it is usually treated as rigid body, and its function is to form a part of the fundamental gravitational field. Only the galactic disk is part of the density wave. In terms of composition, most of the matter in the galactic disk is stars, but there is also a small amount of gas. Because of the great distance between stars, the average ‘collision interval’ is several orders of magnitude larger than the age of the galaxy, so the collision can be omitted.
If the spiral arm of a spiral galaxy is made of fixed matter, and the matter near the center of the galaxy rotates faster than the matter at the edge of the galaxy, then the spiral arm becomes tighter and tighter. After only a few cycles, the arms become entangled and indistinguishable. Density wave theory has succeeded in avoiding this contradiction, with the arm matter not made up of fixed planets but of regions with higher star densities, similar to traffic jams on highways. For dynamic reasons, luminous matter such as a star moves slowly within the spiral arm, lingering for more time and gathering together, but the stars outside the spiral arm move fast and have lower densities.
Despite the great success of density wave theory, there are still some defects and problems. 1. nonlinear effects; 2. analysis of the open arm and barred spiral system; 3. detailed study of the near-central region; 4. dynamical classification of galaxies and so on. When using density wave theory to make galaxy arm pattern, for some cases, the rotary arm pattern is in good agreement with the actual one. However, some other galaxies are not the case, and the barred spiral structure is not clearly explained. The better coincident part is mainly on the outer side of the spiral galaxy, but the results are often far from the actual observations when approaching the center.
The orbit of a star in a non-central potential field is usually a non-closed and complicated spatial curve, which is closely related to the initial velocity and its direction [1]. Therefore, it is more reasonable to take the star system in a galaxy as fluid rather than mass points. Although the dynamics of both models are essentially the same Newtonian mechanics, the initial and boundary conditions are set differently. In the case of fluid, the streamlines are consistent fields, which are a natural result of stars generating from nebula.
W. Dehnen and J. Binney fitted the mass distribution within the galaxy with the observed data (such as the rotation curve and the Oort constant [9]), but found that the mass distribution of the fitting was uncertain. In [10,11], the authors numerically simulated the two-dimensional stellar fluid dynamics embedded in a planar galactic disk in dark matter halo and solved the Boltzmann equation to the second-order moment [1]. The references [12] propose a kinetic approach to explain the formation of spirals and rings in galaxy clusters based on spiral arms, rings, and pseudo-rings driven near the unstable equilibrium point of a given rotating barred potential, and so on, associated with the invariant manifold composed of periodic orbitals around the equilibrium point. By adjusting the kinetic parameters of the main galaxy, the spiral and ringed structures are obtained. The dark matter is beyond the standard model, how to disclose its nature is one of the most challenging problem in astrophysics.
In this paper, we propose a more reasonable and practical galactic dynamics to study the structure and properties of galaxies, which is closely linked to the study of dark matter halo. This framework has the following advantages: 1. The galactic dynamic equations are highly accurate and convenient to solve for the solution, but the previously used equations is not suitable due to unreasonable assumptions. 2. We derived the structure equations for a stable non-warped galaxy. 3. The galactic dynamics is the foundation to analyze and explain other properties of galaxies, on which other problems can be easily solved. We derived and confirmed the experimental law—the Tully–Fisher relation according to the model. 4. By calculating the high precision total density distribution ρ , we can study the dynamic properties of dark matter halo such as the state functions. In simplified dynamical equations, the physical variables are separated from each other. The equations can be solved by fixed procedure, and give a clearer explanation for the spiral and barred structure of galaxies. Thus, this simplified system of equations provides a more reasonable and practical framework for further study of the galactic structure and properties of the dark matter halo.
The dark matter halo provides a dynamical background for the generation and development of stars and planets in a galaxy, and the relationship is similar to that of the amniotic fluid and the fetus. Therefore, the mechanical properties of dark matter are very different from those of ordinary luminous matter, and their moving trajectories are also very different. Technical schemes for detecting dark matter particles by means of detecting ordinary particles are not feasible, because dark matter particles do not interact with ordinary particles at all and any other interaction will result in instable structure of a galaxy. Dark matter detection can only be performed by dynamic methods.

2. Simplification of Galactic Dynamic Equation

In the current theoretical analysis, the commonly used galactic dynamics is similar to the following one [13]
2 Φ = κ ρ ,
( t + V · ) V + 2 Ω × V = Φ 1 ρ P ,
t P + · ( P V ) = P ( γ 1 ) · V Λ ,
t ρ + · ( ρ V ) = 0 ,
in which κ = 4 π G , Ω is the angular speed of the galactic disk, γ adiabatic index, Λ the cooling term reflecting the gravitational disturbance. Some more complicated models also introduce viscous term.
Next we take the dynamic Equations (1)–(4) as example to analyze several problems in current galactic dynamics. 1. The Einstein’s field equation of gravity is essentially a wave equation. In [14], the measurement discloses that gravitational waves travel at c to a precision of 10 15 . The diameters of galaxies are over 100,000 light-years, so the effect of retarded potential of gravity should be taken into account in solving the dynamic equations in such cases. The stars near the center of a galaxy turn around the center several times before the change of the corresponding gravitational field reaches the edge of the galaxy. Whereas the Newton’s gravity is an action at a distance and ignores this time delay. So there is a hidden danger to use (1) to calculate the gravitational field in such an occasion, and the relativistic effect must be considered in the large-scale structure.
The weak field and low velocity approximation of the Einstein’s field equation is as follows (see Appendix A)
α α Φ = κ ρ ,
in which ρ is the total equivalent mass density. The solution of (5) is retarded potential
Φ ( t , x ) = κ 4 π c 2 r c t ρ ( c t r , y ) r d 3 y , r = | x y | .
In contrast with the Newtonian gravitational potential of (1), the retarded potential has finite propagating speed c for perturbation. Not only in the cases of high-speed we need relativity, but also in the cases of large distance or large structure we need relativity. The galactic dynamics must be discussed by general relativity, because the characteristic time of gravity propagating from the center to the edge is much larger than the moving period of the stars near the center, and this time difference has strong effect on the structure of the galaxy.
2. In the early stages of galaxy formation, luminous matter consists of hot plasma. The evolution process is complex at this time, and the dynamics should indeed take into account the influence of pressure, temperature, viscosity, and even magneto fluid effect. However as galaxies mature, luminous matter forms stars, and the motion of the stellar system becomes more and more orderly, which is similar to what happened in the solar system. The stars and ordinary luminous matter in the galaxy are mainly moving along geodesics, and these materials form zero-pressure and non-viscous fluid, so the energy-momentum tensor of the fluid should be as follows
T s μ ν = ρ s U μ U ν ,
where ρ s is the comoving mass density of stars and ordinary luminous particles, U μ is the 4-vectors speed distribution of ordinary matter. The ordinary matter satisfies the energy conservation law T s ; ν μ ν = 0 independent of dark halo. In the form of continuity equations and equations of motion, we obtain the dynamic equations of the stellar system as
U μ μ ρ s + ρ s U ; μ μ = 0 ,
U ν U ; ν μ = 0 .
For convenience, we take c = 1 as unit of speed. The Equations (2) and (3) derived from continuum theory and thermodynamics are not consistent with the actual conditions as mentioned above. The variables in such equations are strongly coupled, which brings great inconvenience and complexity to the analysis and solution of the equations.
3. The equation of state of the dark halo is different from that of ordinary luminous interstellar matter. For dark halo with self-action nonlinear potential, its energy-momentum tensor has the following form [15,16]
T μ ν = ( ρ tot + P ) U μ U ν + ( W P ) g μ ν ,
where W = W ( ρ tot ) corresponds to the nonlinear potential of the dark particles, which has the effect of negative pressure and causes the moving trajectories of the dark halo to deviate from geodesics. For dark spinors, we have W ρ P , so it also acts like dark energy. According to the energy-momentum conservation law T ; ν μ ν = 0 , we obtained the continuity equation U μ T ; ν μ ν = 0 and the equation of motion of the dark halo
( ρ tot + P ) U ν U ; ν μ = ( g μ ν U μ U ν ) ν ( P W ) ,
which is quite different from the geodesic Equation (8). Thus the dark halo is automatically separated from ordinary matter during galaxy formation. So near the solar system, there are few traces of the dark matter that dominates galaxies. Dark halo consists of diffuse gases, and the equation of motion (10) is a complete 1 + 3 dimensional nonlinear field equation. To analyze or solve such equation system is extremely difficult. However, the equation system (6) can be regarded as a planar dynamics in the domain that deviates slightly from the center of the galaxy, and that greatly simplifies the analysis and solving process for the equations. The galactic dynamics currently used confuses the two different kinds of concepts in (6) and (9). For example, in Equations (1)–(4), the mass density ρ in each equation should be different concept. Therefore, the generality and effectiveness of the traditional galactic dynamic equations are quite weak, and it is difficult to obtain the right results in accordance with the actual situation.
Taking into account all above factors, we obtain the galactic dynamic equations in Minkowski space-time under weak field and low speed approximations of relative error ∼ 10 6 (see Appendix A)
α α Φ = κ ρ ,
d d t V ( t + V · ) V = Φ ,
d d t ρ s ( t + V · ) ρ s = ρ s · V ,
where V ( t , x ) is the velocity field of ordinary matter, and the total equivalent mass density is given by
ρ = ˙ ρ tot 2 W + 3 P .
In the above galactic dynamic equations, the variables ( ρ , V , ρ s ) are relatively separated from each other, and the equations for Φ and ρ s are linear and easy to be solved. These features are significant advantages over traditional dynamical Equations (1)–(4). If the total mass density ρ is known, the gravitational potential of Φ ( t , x ) can be integrated, and the velocity field of V or stellar orbits can be obtained from (12)
d 2 d t 2 X ( t ) = Φ ( t , X ) .
If the gravitational potential Φ is independent of time, the mechanical energy of each star is conserved. Integrating the above equation, we obtained
1 2 V 2 + Φ ( X ) = E ( X 0 ) .
Thus, the inappropriate approximation in traditional galactic dynamics brings unnecessary complexity of the equations.
In theory, the complete dynamic equation should include the continuity equation and the dynamic equation of the dark halo. However these equations rely on the equations of state for dark matter and dark energy, which is now an unsolved puzzle. The following analysis shows, the lack of this knowledge can be compensated by observational data. For example, by the rotational velocity curves of the stars, we can then conversely derive some important information about the total density distribution and spiral structure of the dark halo. In some sense, we study the inverse problem of the literature [9,12].

3. Structural Equations of Non-Warped Stationary Galaxies

In the context of general relativity, the structure of galaxies depends on the distribution and properties of dark halo, which is unknown in the present situation [15,16,17]. Because the gravity in the galactic disk is very weak, except for the region near the center, the single gravitational potential Φ can describe the spiral structure precisely enough. We express equations in spherical coordinate system ( t , r , θ , ϕ ) . In this case, we have g k l = d i a g ( 1 , r 2 , r 2 sin 2 θ ) and the non-zero Christoffel symbols
Γ θ θ r = r , Γ ϕ ϕ r = r sin 2 θ , Γ r θ θ = Γ r ϕ ϕ = 1 r , Γ ϕ ϕ θ = sin θ cos θ , Γ θ ϕ ϕ = cot θ .
By 3-d tensor calculus, (11) and (12) becomes
c 2 t 2 Φ g k l ( k l Γ k l m m ) Φ = κ ρ ,
t V m + V k ( k V m + Γ k l m V l ) = g m l l Φ .
Substituting (17) into the above equations, we obtained the dynamic system in spherical coordinate system as
c 2 t 2 Φ ( r 2 + 2 r r + 1 r 2 L ^ 2 ) Φ + κ ρ = 0 ,
D t V r r V θ 2 r sin 2 θ V ϕ 2 + r Φ = 0 ,
D t V θ + 2 r V r V θ sin θ cos θ V ϕ 2 + θ Φ r 2 = 0 ,
D t V ϕ + 2 r V r V ϕ + 2 cot θ V θ V ϕ + ϕ Φ ( r sin θ ) 2 = 0 ,
where the derivative operators are defined as
L ^ 2 = θ 2 + cot θ θ + 1 sin 2 θ φ 2 , D t = t + V r r + V θ θ + V ϕ ϕ ,
and the velocity of the stars is defined in the form of a three-dimensional contra-variant vector
V r d d t r ( t ) , V θ d d t θ ( t ) , V ϕ d d t ϕ ( t ) .
Equations (21)–(23) is the Newton’s second law for stellar motion.
For the stars moving in the galactic disk, we have
θ = π 2 , V θ 0 .
In this case (22) holds automatically, so we obtained the two-dimensional equation of motion and continuity equation for stellar fluid in non-warped disk as [3,10,15]
( t + V r r + V ϕ ϕ ) V r r V ϕ 2 + r Φ = 0 ,
( t + V r r + V ϕ ϕ ) V ϕ + 2 r V r V ϕ + 1 r 2 ϕ Φ = 0 ,
( t + V r r + V ϕ ϕ ) Σ + ( r V r + ϕ V ϕ + 1 r V r ) Σ = 0 ,
where Σ 0 is the mass surface density of stars in the disk.
For quasi-stationary spiral structures, the galactic structure moves as a whole at constant angular velocity Ω around the z-axis. Therefore, under the coordinate transformation φ = ϕ Ω t , the solution of the structural equation in the new coordinate system ( t , r , θ , φ ) is static, that is, the solution is independent of time t. So we obtained the following structural equations for non-warped galaxies
( r 2 + 2 r r + 1 r 2 L ^ 2 Ω 2 c 2 φ 2 ) Φ = κ ρ ,
( V r r + V φ φ ) V r r ( V φ + Ω ) 2 + r Φ = 0 ,
( V r r + V φ φ ) V φ + 2 r V r ( V φ + Ω ) + 1 r 2 φ Φ = 0 ,
( V r r + V φ φ ) Σ + ( r V r + φ V φ + 1 r V r ) Σ = 0 .
The terms related with Ω in the above formulas have obvious physical significance, r ( V φ + Ω ) 2 is centrifugal force, and 2 r V r Ω Coriolis force. Because of the background gravitational force, Equations (29)–(31) is represented in spherical coordinate system, but (32) is described in polar coordinate system.

4. Second-Order Approximation of Structural Equation

The Equations (29)–(32) can be solved according to the following procedure: At first, a reasonable total mass density ρ is assumed, which can be expressed as series expansion of spherical harmonic function Y l m ( θ , φ ) . The distribution of the dark halo can be determined inversely by comparing the calculated solutions with the observed data. Second, for given ρ , (29) is a linear equation with respect to Φ , which can be easily solved with natural boundary conditions. Third, by (30) and (31), the velocity function can be represented by the trigonometric series of φ . Due to symmetry, the trigonometric series has only the terms with an integer multiple of a fixed base frequency, usually only even terms. The fourth, Equation (32) for the surface mass density is a first-order linear equation, which is also easy to solve. In this way, the stellar distribution and movement are perfectly determined for non-warped stationary galaxies.
In what follows, we demonstrate the solving process of using the triangular series expansion according to this procedure. Consider the following second-order approximation of mass and gravitational potential distributions
ρ = ρ 0 + ( ρ 1 + ρ 2 cos 2 φ + ρ 3 sin 2 φ ) sin 2 θ ,
Φ = Φ 0 + ( Φ 1 + Φ 2 cos 2 φ + Φ 3 sin 2 φ ) sin 2 θ ,
in which all ( ρ n , Φ n ) are functions of r and satisfy ρ 0 . The relative error of (33) and (34) is about 3 2 s < 12 % , where s 1 reflects the degree of mass distribution similar to (33). The corresponding second order velocity of the stars in the disk is given by
V r = v 1 cos 2 φ + v 2 sin 2 φ , V θ = 0 ,
V φ = ω 0 + ω 1 cos 2 φ + ω 2 sin 2 φ ,
in which all ( v n , ω n ) are functions of r.
Substituting (33) and (34) into (29), we obtained the equation of Φ n as
Φ 0 + 2 r Φ 0 + 4 r 2 Φ 1 = κ ρ 0 ,
Φ 1 + 2 r Φ 1 6 r 2 Φ 1 = κ ρ 1 ,
4 Ω 2 c 2 Φ k + Φ k + 2 r Φ k 6 r 2 Φ k = κ ρ k ,
where k = 2 , 3 . The solution is given by
Φ 1 = κ 5 [ 1 r 3 0 r ρ 1 ( χ ) χ 4 d χ + r 2 r ρ 1 ( χ ) χ d χ ] ,
Φ 0 = 1 r 0 r ( χ 2 κ ρ 0 4 Φ 1 ) d χ r 1 χ ( χ 2 κ ρ 0 4 Φ 1 ) d χ
and
Φ k = κ c 5 32 Ω 5 r 2 [ A ( r ) 0 r B ( χ ) ρ k ( χ ) d χ + B ( r ) r A ( χ ) ρ k ( χ ) d χ ] ,
in which
A ( r ) = 6 Ω c cos ( 2 Ω r c ) + 4 Ω 2 r 2 3 c 2 r c 2 sin ( 2 Ω r c ) , B ( r ) = 4 Ω 2 r 2 3 c 2 r c 2 cos ( 2 Ω r c ) 6 Ω c sin ( 2 Ω r c ) .
Substituting (34)–(36) into (30) and (31), and restricting θ = 1 2 π , we obtained
0 = d d r [ Φ 0 + Φ 1 + 1 4 ( v 1 2 + v 2 2 ) ] 1 2 R [ 2 ( ω 0 + Ω ) 2 + ω 1 2 + ω 2 2 ] + ω 1 v 2 ω 2 v 1 + [ 2 ω 0 v 2 2 r ( ω 0 + Ω ) ω 1 + Φ 2 ] cos 2 φ [ 2 ω 0 v 1 + 2 r ( ω 0 + Ω ) ω 2 Φ 3 ] sin 2 φ +
as well as
0 = 1 2 [ ( r ω 1 + 2 ω 1 ) v 1 + ( r ω 2 + 2 ω 2 ) v 2 ] + [ ( r ω 0 + 2 ω 0 + 2 Ω ) v 1 + 2 r ω 0 ω 2 + 2 r Φ 3 ] cos 2 φ + [ ( r ω 0 + 2 ω 0 + 2 Ω ) v 2 2 r ω 0 ω 1 2 r Φ 2 ] sin 2 φ .
Let the coefficients of the terms ( sin 2 φ , cos 2 φ ) in (43) and (44) be 0, we obtained
v 1 = K [ ω 0 r 2 Φ 3 + 2 r ( ω 0 + Ω ) Φ 3 ] ,
v 2 = K [ ω 0 r 2 Φ 2 + 2 r ( ω 0 + Ω ) Φ 2 ] ,
ω 1 = K [ ( 1 2 r 2 ω 0 + r ( ω 0 + Ω ) ) Φ 2 + 2 ω 0 Φ 2 ] ,
ω 2 = K [ ( 1 2 r 2 ω 0 + r ( ω 0 + Ω ) ) Φ 3 + 2 ω 0 Φ 3 ] ,
K { r 2 [ r ( ω 0 + Ω ) ω 0 + 2 Ω ( 2 ω 0 + Ω ) ] } 1 .
Substituting the above results into (44), that is, into the equation
( r r ω 1 + 2 ω 1 ) v 1 + ( r r ω 2 + 2 ω 2 ) v 2 = 0 ,
we obtained a first-order linear ordinary differential equation for ω 0
ω 0 = 1 r ( 3 ω 0 + Ω ) + ω 0 r 2 · Φ 3 Φ 2 Φ 2 Φ 3 Φ 3 Φ 2 Φ 2 Φ 3 + ( ω 0 + Ω ) Φ 3 Φ 2 Φ 2 Φ 3 Φ 3 Φ 2 Φ 2 Φ 3 ,
and another solution with no physical meaning ω 0 = 2 r ( ω 0 + Ω ) . It is important to note that, while Φ 2 0 or Φ 3 0 , ω 0 cannot be determined by above Equation (51). In this case (50) holds automatically. Substituting (45)–(49) into (43), i.e., substituting them into
d d r [ Φ 0 + Φ 1 + 1 4 ( v 1 2 + v 2 2 ) ] = ω 1 v 2 + ω 2 v 1 + 1 2 r [ 2 ( ω 0 + Ω ) 2 + ω 1 2 + ω 2 2 ] ,
we obtained a constraint between the mass density distribution and Ω , so Ω is determined by the total mass density. From (45)–(51) we learn that, the stellar speed V is completely determined by ( ρ 2 , ρ 3 , Ω ) . After the velocity distribution is determined, the surface mass density Σ ( r , φ ) of the stars can be solved from (32). At this point, the second-order approximation of galactic dynamics is completely solved. The above process only involves solving ordinary differential equations.

5. Simple Solutions to Stationary Galaxy Structure

The calculations in the preceding section show that the structure of a stationary galaxy is determined by ( ρ 2 , ρ 3 , Ω ) . If the total mass density distribution is known, it is easy to determine the structure of the galaxy according to the above procedure. However, it is difficult to directly determine the density distribution. The data observed are mainly the moving speeds of the stars. High-precision observation show that the stellar speeds in galaxies in a larger range is approximately equal. To illustrate the above solving procedure in detail and obtained some enlightening conclusions, we concretely solve two of the simplest examples. As a working hypothesis, we assume [18,19],
ω 0 = v r Ω , ( r [ R 0 , R 1 ] ) ,
where v is a constant velocity with a typical value of | v | = 200 ∼400 km/s, [ R 0 , R 1 ] is the effective domain of assumption (53). Usually we have R 0 = 100 ∼500 pc and R 1 = 10 ∼60 kpc, which are about the radius of the bulge and the radius of the visible range, respectively. In the bulge region r < R 0 , the mass distribution can be approximated by a rigid sphere. No loss of generality, we can set ( ω 0 > 0 , v > 0 ) , otherwise we can make a transformation ϕ ϕ to obtained this condition. So we always have v > r Ω in the effective domain. As shown below, using the star’s velocity curve v = v ( r ) to calculate the total mass density is a good method to study the properties of dark halo [15,20,21,22,23].

5.1. Barred Spiral Galaxy

We consider the simple case ρ 2 0 , ρ 3 = 0 . In this case we have P 3 = 0 , and P 2 ( r ) can be calculated by (42). By (45)–(49), we obtained
v 1 = 0 , v 2 = r ( v r Ω ) P 2 + 2 v P 2 2 ( v r Ω ) 2 v 2 ,
ω 2 = 0 , ω 1 = r v P 2 + 4 ( v r Ω ) P 2 2 r [ 2 ( v r Ω ) 2 v 2 ] .
The effective region of the above formulas is r < ( 2 2 ) v 2 Ω . Substituting (53)–(55) into (32), we obtained the surface density of stars as
Σ = [ Σ 0 + C 2 r v 2 cos ( 2 φ ) ] 1 ,
in which
Σ 0 = r v 2 e 2 ω 1 v 2 d r [ C 1 + 2 C 2 v r Ω r v 2 e r 2 ω 1 ( χ ) v 2 ( χ ) d χ d r ] ,
( C 1 , C 2 ) are integration constants.
The above results show the structural characteristics of a barred spiral galaxy. For given total mass density ρ and boundary values, all the above formulas can be calculated concretely, but the analytical expression is a little complicated. If a restriction v 2 = 0 is added, that is V r 0 , then the stars move around circles, the results are relatively simple. In this case, by (54), we have
Φ 2 = q r 2 ( v r Ω ) 2 , Φ 3 = 0 ,
in which q 0 is an integral constant. Substituting (57) into (45)–(49) we obtained
V r = 0 , V φ = 1 r 3 ( v r Ω ) ( r 2 + q cos 2 φ ) .
Again by (32), we obtained the surface mass density of stars
Σ = ϱ ( r ) r 2 + q cos 2 φ ,
in which ϱ ( r ) is determined by boundary condition. Substituting the above results into (52), we obtained a restriction between Ω and potential
Φ 0 + Φ 1 = v 2 ln r r 0 q 2 24 r 4 ( 3 v 2 8 v r Ω + 6 r 2 Ω 2 ) ,
in which r 0 > 0 . In the case of q = 0 , the solution (57)–(60) is exact, which correspond to elliptical galaxy. However, as long as the conditions deviate a little, the elliptical galaxy will evolve into a spiral one. The spiral structure should be an inevitable result for the evolution of ordinary galaxies.
As an application of (60), we derive the Tully–Fisher relation L v 4 [24]. In astrophysics, we have the following empirical laws for a galaxies
L M R 2 ,
in which L is the luminosity, M the total mass and R the visible radius. For an elliptical galaxy, in (60) we have
q = Ω = 0 , | Φ 1 | | Φ 0 | .
Then, we obtained
4 π G ρ = Φ + 2 r Φ = ˙ v 2 r 2 .
So we have relation
M = 4 π 0 R ρ r 2 d r = ˙ v 2 G R .
Taking v as known parameter and solving (61) and (64), we obtain relations
R v 2 , L M v 4 .
In general case, if the perturbative terms Ω , V r and L ^ 2 in (29) and (30) can be omitted and v is almost a constant, then we have equations
r Φ = ˙ v 2 r , 4 π G ρ = r 2 Φ + 2 r r Φ = ˙ v 2 r 2 .
Equation (66) implies that the Tully–Fisher relation also approximately holds for all stable galaxies.

5.2. Spiral Galaxy

Now, we calculate the following simple case of Ω = 0 ,
ω 0 = v r , Φ 2 = P 1 cos ( ξ r ) , Φ 3 = P 1 sin ( ξ r ) ,
where ξ > 0 is a constant, P ( r ) is a function of r. Substituting (67) into (51), we obtained a linear differential eguation for P ( r )
P 4 r P + ( ξ 2 + 8 r 2 ) P = 0 .
The solution is given by
P = r 5 [ C 1 J α ( ξ r ) + C 2 J α ( ξ r ) ] , α = 1 2 7 ,
in which ( C 1 , C 2 ) are constants, ( J α , J α ) is Bessel function of imaginary parameter
J ν ( x ) = x 2 ν k = 0 1 k ! 1 Γ ( ν + k + 1 ) x 2 2 k , ( ν = ± α ) .
Substituting the above relations into (45)–(49), we have
Φ = Φ 0 + [ Φ 1 + P 1 cos γ ] sin 2 θ , V r = 1 v P 2 [ ξ r P cos γ ( r P 2 P ) sin γ ] , r V φ = 1 2 v P 2 [ ( r P 4 P ) cos γ + ξ r P sin γ ] + v , Σ = C 3 P 2 ξ r 5 ( r P 2 P ) cos γ + ξ r P sin γ + 2 r 2 v 2 P 2 r 3 d r ,
in which γ = ξ r 2 φ . All of these variables are characterized by the structure of a spiral galaxy, which are closer to the Archimedes spiral rather than the logarithmic one. Since the pitch angle of galaxy is approximately 10 40 [25], we have estimation ξ = ( 2 11 ) R 1 1 .

6. Discussion and Conclusions

To explain the galaxy structure and flat velocity curves, [26] declared that it is unnecessary to introduce ‘hidden mass hypothesis’, and the following modified Newtonian dynamics was proposed,
M g μ ( a / a 0 ) a = F .
In order to obtain the Tully–Fisher relation, in [27] the centripetal acceleration a = V 2 / r and the following gravity were used,
g N = M G r 2 .
Since the visible mass distribution in galaxies is disc-like, if the dark matter is absent, then the main part of the Newtonian gravitational potential of a galaxy should be similar to that of a disc mass distribution [28]
Φ d = 2 G M R 1 n = 0 C 1 / 2 n r R 1 2 n P 2 n ( cos θ ) + 2 G M r R 1 2 | cos θ | , r R 1 ,
Φ d = 2 G M R 1 n = 0 C 1 / 2 n + 1 R 1 r 2 n + 1 P 2 n ( cos θ ) , r > R 1 ,
where C 1 / 2 n is defined by 1 + x = C 1 / 2 n x n . In the disc ( r R 1 , θ = 1 2 π ) , we have cos θ = 0 and
Φ d = 2 G M R 1 [ 1 1 4 r R 1 2 + ] , g N = G M r R 1 3 + .
How can we obtain (72)? Equation (72) is actually a proof of the existence of dark matter, because it is just the Newtonian gravity of the spherically symmetric distribution of mass. The reasonable explanation for the Tully–Fisher relation L v 4 should be (63)–(66).
In addition, Equation (71) is also not a good model in physics, because it obviously violates some fundamental principles of physics such as covariance and the Poisson equation or wave equation. For problems as complex as galactic dynamics, a systematic research based on the first principles is needed. Ref. [26] only used 5 simple algebraic formulas without any logical relation between them: partial assumptions, partial fitting to the observed data and most of the content is stating facts, opinions and conjectures. For complex problems, this style of research is difficult to effectively solve the problem, and is easily disoriented. One without sharp logic can hardly check the validity of the results or find out where a bug exists. From the review paper [29], we can find a number of studies misled by the model, much of time and resource was wasted by the professional barriers. The simple treatment is feasible for organizing experimental data and laws, but the galactic dynamics and structure involve large scale space-time geometry and partial differential equations of gravity and matter distribution, so it is impossible to obtain profound results by simple formula and conjecture. In this case, open mind and multidisciplinary collaboration is important, because some problems that are difficult for the experts in the subject but simple for others.
For the mature galaxies, in the paper three working hypotheses are introduced to simplify the dynamic equations: 1. To research the Large-scale structure of galaxy, we should consider the retarded potential of gravity, which takes longer time to propagate to border of galaxy than the revolution period of the stars near the center. In this case, to use Newtonian gravity directly is unreasonable and it should be replaced by wave equation, i.e., the weak field and low velocity approximation of Einstein’s field Equation (5).
2. The stellar system of mature galaxies should be zero-pressure and viscosity–Free fluid having different equations of motion from ordinary continuum mechanics. The pressure of such system has no physical meaning and cannot be defined, and the introduction of pressure will only bring errors and troubles. The star moves along geodesic and the dynamic equation of stellar system should be (7) and (8).
3. The equation of state of dark halo is different from that of ordinary luminous interstellar matter, so their trajectories are very different, and the dark halo is automatically separated from ordinary matter in galaxies. The structure of galaxies is only determined by the total mass density, which can be assumed beforehand based on observations. Then, the total mass density can be determined by comparing the solution of structural equations with the observed data.
By these processes, we find that the variables in the dynamical equations of galaxies are separated from each other. The equations are well-posed and can be solved according to a fixed procedure. Traditional equations include unreasonable assumptions and are too complicated to be studied in depth. In fact, the basic objects of Nature are designed elaborately, and the corresponding equations are simple and symmetrical. If our research is trapped in complex computation and incomprehensible dilemma, this situation reflects that the model used has gone wrong with big probability.
The dynamic equation system derived above provides a more reasonable and practical framework for simulating the structure of galaxies and is closely linked to the study of dark matter halo. The physical significance of the above framework is as follows: 1. A system of highly accurate and convenient galactic dynamic Equations (11)–(13) is established, but the previously used equations such as (1)–(4) is not suitable for galactic structure. 2. For a mature galaxy, its structure is stable, and (29)–(32) gives the structure equations for a stable non-warped galaxy. 3. Equations (11)–(13) is the foundation to analyze and explain other properties of galaxies, on which other problems can be easily solved, such as the Tully–Fisher relation derived by (66). It is clear that the discreteness of the Tully–Fisher relation between different types of galaxies can also be derived by calculating the errors of (66) for various galaxies. 4. By calculating the fine structure of the total density distribution ρ in (11), we can study the dynamic properties of dark matter halo, such as the state function W ( ρ ) and so on. In conclusion, galactic dynamics is the theoretical foundation for the study of galactic structure and properties, and a prerequisite for solving other problems. The previous dynamic equations included unreasonable assumptions, while the framework in the paper provides a reasonable dynamic model of higher accuracy and more convenience.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

I am grateful to an anonymous reviewer; this paper has been modified according to their suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Linearization of Einstein Field Equation

In this section, by weak-field and low-speed approximation of Einstein field equation, we will establish the dynamic equations describing the galactic system. Some fundamental contents can be found in [30], but here we make a more systematic and detailed derivation for galactic dynamics.
In the context of general relativity, the whole system of the dynamic equations for galactic evolution should be Einstein field equation
G μ ν R μ ν 1 2 g μ ν R = κ ¯ T μ ν , ( κ ¯ 8 π G c 4 ) ,
combined with the law of energy-momentum conservation and the equation of state of the gravitating source. In most cases, the classical approximation for the total energy-momentum tensor takes the following form [16]
T μ ν = ( ρ + P ) U μ U ν + ( W P ) g μ ν ,
where W corresponds to the nonlinear potentials, which acts as negative pressure and leads to the deviation from geodesic. According to the law of energy-momentum conservation or the Bianchi identity T ; ν μ ν = 0 , we can derive the continuity equation U μ T ; ν μ ν = 0 and the equation of motion for the source as follows
U μ μ ( ρ + W ) = ( ρ + P ) U ; μ μ ,
( ρ + P ) U ν U ; ν μ = ( g μ ν U μ U ν ) ν ( P W ) .
In the case of W ρ P , such as for the nonlinear spinor, we find that the stream lines are quite different from the geodesics U ν U ; ν μ = 0 . Therefore, a fully relativistic simulation of the galactic evolution should include such terms. However, in the following non-relativistic approximation, the effects of ( P , W ) are merged into an effective mass-energy density ρ , which becomes much simpler.
For convenience, we take c = 1 as the unit of velocity. Noting the facts that collisions between stars rarely occur, the trajectories of the ordinary matter, such as atoms, are almost geodesics, so for stars, the following zero-pressure and inviscid energy-momentum tensor holds
T s μ ν = ρ s U μ U ν ,
in which ρ s is the comoving mass density of stars, and U μ is the 4-vector speed of the stellar flow. For energy-momentum tensor of a compound system, we have the following useful theorem [31], which means the energy-momentum of any independent subsystem is conserved, respectively.
Theorem A1.
Assume matter consists of two subsystems I and II, namely L m = L I ( ϕ ) + L I I ( ψ ) , then we have
T μ ν = T I μ ν + T I I μ ν .
If the subsystems I and II have not interaction with each other, namely,
δ δ ψ L I ( ϕ ) = δ δ ϕ L I I ( ψ ) = 0 ,
then the two subsystems have independent energy-momentum conservation laws, respectively,
T I ; ν μ ν = 0 , T I I ; ν μ ν = 0 .
According to Theorem A1, ordinary matter satisfies the law of energy-momentum conservation independent of that of the dark halo, so we have T s ; ν μ ν = 0 . Expressing it in the form of equations of continuity and motion, we obtain the dynamic equations for the stars
U μ μ ρ s + ρ s U ; μ μ = 0 , U ν U ; ν μ = 0 .
The total energy-momentum tensor of the galaxy is still given by (A2), and satisfies the dynamic Equations (A3) and (A4). Using (A1) and (A2), we obtain
R = κ ¯ ( ρ + 4 W 3 P ) ,
where R = g μ ν R μ ν is the scalar curvature. Substituting (A9) into (A1), we obtain
R μ ν = κ ¯ ( ρ + P ) U μ U ν + 1 2 κ ¯ ( ρ + 2 W P ) g μ ν ,
where U μ is the average 4-vector speed of all gravitating source.
In order to make a weak-field approximation, we choose the harmonic coordinate system, which leads to the usual Cartesian coordinate system when linearizing of metric. Then we have the de Donder coordinate condition
Γ μ g α β Γ α β μ = 1 g ν ( g g μ ν ) = 0 ,
where g = | det ( g ) | . Denote the Minkowski metric by η μ ν = η μ ν = d i a g ( 1 , 1 , 1 , 1 ) . For the weak-field approximation, we have the linearization for the metric
g μ ν η μ ν + h μ ν , g μ ν = ˙ η μ ν h μ ν , h μ ν = η μ α η ν β h α β , h = h μ μ = η μ ν h μ ν , g = ˙ 1 + h , g = ˙ 1 + 1 2 h .
In the above approximation we have a precision O ( | h μ ν | 2 ) Φ 2 + V 2 / c 2 10 6 , which is sufficient for calculations of galactic dynamics. Next, we directly use = to replace = ˙ for convenience. By straightforward calculation, we obtain the linearization for other parameters
Γ α β μ = 1 2 η μ ν ( α h ν β + β h α ν ν h α β ) , Γ μ = ν ( h μ ν 1 2 η μ ν h ) , R μ ν = 1 2 α α h μ ν 1 2 ( η μ α ν Γ α + η ν α μ Γ α ) , R μ ν = 1 2 α α h μ ν 1 2 ( η μ α α Γ ν + η ν α α Γ μ ) , R = 1 2 α α h α Γ α .
In which α α = t 2 2 is the d’Alembert operator.
In the harmonic coordinate system, we have
Γ μ = ν ( h μ ν 1 2 η μ ν h ) = 0 ,
R μ ν = 1 2 α α h μ ν , R μ ν = 1 2 α α h μ ν ,
R = 1 2 α α h , G μ ν = 1 2 α α ( h μ ν 1 2 η μ ν h ) .
From (A11) and (A13), we find that if Γ μ = 0 , t Γ μ = 0 at any given time t = t 0 , it will always hold due to the Bianchi identity G ; ν μ ν = 0 .
To compare with electromagnetism and to understand the physical meaning of the parameters, denote
Φ = 1 2 h t t = 1 2 h t t , A = ( h t x , h t y , h t z ) = ( h t x , h t y , h t z ) ,
H = ( h a b ) = ( h a b ) , ( { a , b } { 1 , 2 , 3 } ) , B = × A .
In the International System of Units, we have the order of magnitude of the metric components
c 2 | h a b | c | A k | | Φ | 1 , ( a b ) ,
which means | h a b | | A k | | Φ | 1 if we take c = 1 as the unit.
For the present purpose, we define the stellar speed V by
V 1 U 0 ( U 1 , U 2 , U 3 ) ,
which is approximately equivalent to the usual definition. For galaxies, we have the following order of magnitude
| V | 300 km / s = 10 3 c , A κ ¯ V , h a b κ ¯ | V | 2 , ( a b ) ,
in which the coefficient κ ¯ is also a number of small value. Then, according to
1 = g μ ν U μ U ν = ( 1 + 2 Φ 2 A · V + g a b V a V b ) 1 2 U 0 ,
by omitting the O ( V 2 ) terms, the low-speed assumption gives
U 0 = 1 Φ + A · V .
Substituting (A17) and (A18) into (A8) and omitting the higher-order terms, we obtain the continuity equation and motion equation for the stars
( t + V · ) ρ s = ρ s [ · V + ( t Φ + · A ) ] ,
( t + V · ) V = Φ + ( t A + V t Φ ) + V × B + V · t H .
In (A19), we used the de Donder condition Γ 0 = 0 in the form
1 2 t ( h x x + h y y + h z z ) = ( t Φ + · A ) .
The equation of motion (A20) is similar to electrodynamics. From it we learn that, Φ gives the Newtonian gravitational potential, and A leads to the gravimagnetic field B .
By (A9) and (A13), we have
α α h = 2 κ ¯ ( ρ + 4 W 3 P ) .
By (A10), (A12) and (A22), we obtain the dynamic equations for h μ ν
α α h μ ν = 2 κ ¯ ( ρ + P ) U μ U ν + κ ¯ ( ρ + 2 W P ) η μ ν ,
α α χ μ ν = 2 κ ¯ [ ( ρ + P ) U μ U ν + ( W P ) η μ ν ] ,
where χ μ ν = h μ ν 1 2 η μ ν h . If the average speed of the dark halo is also small, omitting O ( U 2 ) from (A23), we obtain h x x = h y y = h z z 2 Ψ , h a b = 0 , ( a b ) and
α α Φ = 1 2 α α h 00 = 4 π G ρ ,
α α Ψ = 1 2 α α h k k = 4 π G ρ ˜ ,
where ρ and ρ ˜ are the effective mass densities with little difference. Their zeroth-order approximation gives
ρ = ρ [ 2 ( U 0 ) 2 1 ] 2 W + P [ 2 ( U 0 ) 2 + 1 ]
= ˙ ρ 2 W + 3 P .
ρ ˜ = ˙ ρ + 2 W P .
In the following discussion, only the zeroth-order approximation of (A19), (A20) and (A25) are involved. The precision of these equations is sufficient for calculations of galactic dynamics in the effective domain r [ R 0 , R 1 ] .
In order to obtain the the Hamiltonian formalism of Dirac equation in curved soace-time, we need the natural coordinate system(NCS) [32]
d s 2 = g t t d t 2 g ¯ k l d x k d x l , d τ = g t t d t , d V = g ¯ d 3 x .
In which d s is the 4-d length of line element, d τ is the Newton’s absolute cosmic time element and d V is the absolute volume element of space at time t. The NCS generally exists and the global simultaneity is unique. Only in NCS can we clearly establish the Hamiltonian formalism and calculate the Noether charges.
The NCS is a wonderful coordinate system different from the Gaussian coordinate system that is valid only in the neighborhood of the initial Cauchy hypersurface. The isodistant translating hypersurface will deform soon, so that the metric d s 2 = d t 2 + g k l d x k d x l becomes invalid. The NCS is also different from the Einstein’s lift moving along a geodesic, namely the co-moving coordinate system, as this requires the lift to be an infinitesimal volume in curved space-time. While NCS holds unconditionally and globally, and its time is objective cosmic time. g t t represents gravity and cannot be merged into the time coordinate t.
However, the NCS is not suitable for the linearization of the Einstein’s field equations, because for a rotating galaxy the coordinate transformation from harmonic coordinate system to NCS is similar to
r = r , θ = θ , φ = ω ( r , θ ) t + φ ,
and such transformation with rotation cannot be linearized.

References

  1. Binney, B.; Tremaine, S. Galactic Dynamics, 2nd ed.; Princeton University Press: Princeton, NJ, USA, 2008. [Google Scholar]
  2. Orlova, N.; Korchagin, V.; Kikuchi, N.; Miyama, S.M.; Moiseev, A. Application of the global modal approach to the spiral galaxies. Astrophys. Space Sci. 2003, 284, 739. [Google Scholar] [CrossRef]
  3. Meidt, S.E.; Rand, R.J.; Merrifield, M.R.; Debattista, V.P.; Shen, J. Tests of the Radial Tremaine-Weinberg Method. Astrophys. J. 2008, 676, 899. [Google Scholar] [CrossRef]
  4. Lin, C.C.; Shu, F.H. On the Spiral Structure of Disk Galaxies. Astrophys. J. 1964, 140, 646. [Google Scholar] [CrossRef]
  5. Goldreich, P.; Lynden-Bell, D., II. Spiral Arms as Sheared Gravitational Instabilities. Mon. Not. R. Astron. Soc. 1965, 130, 125. [Google Scholar] [CrossRef]
  6. Julian, W.H.; Toomre, A. Non-Axisymmetric Responses of Differentially Rotating Disks of Stars. Astrophys. J. 1966, 146, 810. [Google Scholar] [CrossRef]
  7. Somerville, R.S.; Barden, M.; Rix, H.W.; Bell, E.F.; Beckwith, S.V.; Borch, A.; Wolf, C. An Explanation for the Observed Weak Size Evolution of Disk Galaxies. Astrophys. J. 2008, 672, 776. [Google Scholar] [CrossRef] [Green Version]
  8. Kendall, S.; Kennicutt, R.C.; Clarke, C.; Thornley, M.D. Tracing Spiral Density Waves in M81. Mon. Not. R. Astron. Soc. 2008, 387, 1007. [Google Scholar] [CrossRef]
  9. Dehnen, W.; Binney, J. Mass Models of the Milky Way. Mon. Not. R. Astron. Soc. 1998, 294, 429. [Google Scholar] [CrossRef] [Green Version]
  10. Vorobyov, E.I.; Theis, C. Boltzmann Moment Equation Approach for the Numerical Study of Anisotropic Stellar Discs. Mon. Not. R. Astron. Soc. 2006, 373, 197. [Google Scholar] [CrossRef]
  11. Vorobyov, E.I.; Theis, C. Shape and Orientation of Stellar velocity Ellipsoids in Spiral Galaxies. Mon. Not. R. Astron. Soc. 2008, 383, 817. [Google Scholar] [CrossRef] [Green Version]
  12. Romero-Gómez, M.; Athanassoula, E.M.; Masdemont, J.J.; García-Gómez, C. The Formation of Spiral Arms and Rings in Barred Galaxies. Astron. Astrophys. 2007, 472, 63. [Google Scholar] [CrossRef] [Green Version]
  13. Riols, A.; Latter, H.N. Spiral Density Waves and Vertical Circulation in Protoplanetary Discs. Mon. Not. R. Astron. Soc. 2018, 476, 5115. [Google Scholar] [CrossRef] [Green Version]
  14. Virgo & LIGO Collaborations. GW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral. Phys. Rev. Lett. 2017, 119, 161101. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Gu, Y.Q. Stationary Spiral Structure and Collective Motion of the Stars in a Spiral Galaxy. In Old Problems and New Horizons in World Physics; Christianto, V., Smarandache, F., Eds.; Nova Science Publishers: New York, NY, USA, 2019; p. 323. [Google Scholar]
  16. Gu, Y.Q. Clifford Algebra and Unified Field Theory; LAP LAMBERT Academic Publishing: Chisinau, Moldova, 2020. [Google Scholar]
  17. Mori, M.; Umemura, M. The Evolution of Galaxies from Primeval Irregulars to Present-Day Ellipticals. Nature 2006, 440, 644. [Google Scholar] [CrossRef] [PubMed]
  18. Sofue, Y.; Tutui, Y.; Honma, M.; Tomita, A.; Takeda, Y. Central Rotation Curves of Spiral Galaxies. Astrophys. J. 2009, 523, 136. [Google Scholar] [CrossRef] [Green Version]
  19. Bownstein, J.R.; Moffat, J.W. Galaxy Rotation Curves without Nonbaryonic Dark Matter. Astrophys. J. 2006, 636, 721. [Google Scholar] [CrossRef] [Green Version]
  20. Salucci, P.; Martins, C.F. The Mass Distribution in Spiral Galaxies. EAS Publ. 2009, 36, 133. [Google Scholar] [CrossRef] [Green Version]
  21. Persic, M.; Salucci, P.; Stel, F. Rotation Curves of 967 Spiral Galaxies: Implications for Dark Matter. AIP Conf. Proc. 1995, 336, 137. [Google Scholar]
  22. Persic, M.; Salucci, P.; Stel, F. The Universal Rotation Curve of Spiral Galaxies—I. The Dark Matter Connection. Mon. Not. R. Astron. Soc. 1996, 281, 27. [Google Scholar] [CrossRef] [Green Version]
  23. Widrow, L.M.; Pym, B.; Dubinski, J. Dynamical Blueprints for Galaxies. Astrophys. J. 2008, 679, 1239. [Google Scholar] [CrossRef] [Green Version]
  24. Tully, R.B.; Fisher, J.R. A new method of determining distance to galaxies. Astron. Astrophys. 1977, 54, 661–673. [Google Scholar]
  25. Yu, S.Y.; Ho, L.C. On the Connection between Spiral Arm Pitch Angle and Galaxy Properties. Astrophys. J. 2019, 871, 194. [Google Scholar] [CrossRef] [Green Version]
  26. Milgrom, M. A Modification of the Newtonian Dynamics as a Possible Alternative to the Hidden mass Hypothesis. Astrophys. J. 1983, 270, 365. [Google Scholar] [CrossRef]
  27. Milgrom, M. A Modification of the Newtonian Dynamics—Implications for Galaxies. Astrophys. J. 1983, 270, 371. [Google Scholar] [CrossRef]
  28. Croxton, C.A. Introductory Eigenphysics; John Wiley & Sons Ltd.: New York, NY, USA, 1974. [Google Scholar]
  29. Banik, I.; Zhao, H.S. From galactic bars to the Hubble tension—Weighing up the astrophysical evidence for Milgromian gravity. arXiv 2021, arXiv:2110.06936. [Google Scholar]
  30. Weinberg, S. Gravitation and Cosmology; (Ch5.4, Ch7, Ch.9.8, Ch.10.1); Wiley: New York, NY, USA, 1972. [Google Scholar]
  31. Gu, Y.Q. Dynamical Reason for a Cyclic Universe. Symmetry 2021, 13, 2272. [Google Scholar] [CrossRef]
  32. Gu, Y.Q. Natural Coordinate System in Curved Space-time. J. Geom. Symmetry Phys. 2018, 47, 51–62. [Google Scholar] [CrossRef] [Green Version]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gu, Y.-Q. Simplification of Galactic Dynamic Equations. Symmetry 2022, 14, 407. https://doi.org/10.3390/sym14020407

AMA Style

Gu Y-Q. Simplification of Galactic Dynamic Equations. Symmetry. 2022; 14(2):407. https://doi.org/10.3390/sym14020407

Chicago/Turabian Style

Gu, Ying-Qiu. 2022. "Simplification of Galactic Dynamic Equations" Symmetry 14, no. 2: 407. https://doi.org/10.3390/sym14020407

APA Style

Gu, Y. -Q. (2022). Simplification of Galactic Dynamic Equations. Symmetry, 14(2), 407. https://doi.org/10.3390/sym14020407

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop