Kinematic and Thermodynamic Properties of the Galaxy


We present the first detailed physically-based thermodynamic model of the Galaxy as obtained from high-resolution kinematical data from the inner stellar halo. We interpret the observed distribution of Hyper Velocity Stars (HVS) as a physical manifestation of the Maxwell-Boltzmann (M-B) probability distribution expected for a fully virialized galactic system existing in quasi-equilibrium with its surroundings. The conventional view is that the HVS sample originates from chance gravitational encounters that have attained enough speed to escape the Galaxy. We counter that the HVS population is created by thermodynamic effects and link the observed mid-disk velocity peak with the M-B most probable velocity VPeak ≈ Vmp=432 km/s. Most of the HVS population originates at the Galactic virialization radius of 23 kpc with the current observed sample reproducing the M-B distribution with fidelity.


Hyper Velocity Stars (HVS) have become a valuable tool to constrain the total dynamic mass of the Milky Way. We investigate a thermodynamic origin for the HVS population and compare it against popular ejection (slingshot) mechanisms. We perform a combined analysis (HVS sample and stellar halo kinematics) and comfortably fit a Maxwell-Boltzmann distribution to these data sets, a result consistent with virialized Galaxy models.  In this paper, we advance a physically-based approach conveniently termed the Rotation Curve-Spin Parameter (RC-SP) model to distinguish it from ɅCDM.

The Virial Galaxy

The RC-SP solution employs a baryonic-based interpretation based on the spin parameter equation, and includes angular momentum and total energy associated with the Galactic “state” (Peebles 1971) (La Fortune 2016). In addition to obtaining galactic dynamics from extended rotation velocity profiles, this paper strengthens the RC-SP model via precision measures of inner stellar halo component dispersion velocities and HVS “escape velocity” analysis. The RC-SP approach is based on two classical equations, the Virial Theorem and Newton’s second law for circular motion. The Virial Theorem is expressed below and includes the constraint which limits the theorem to isolated, self-gravitating systems in “equilibrium” and is equally applicable to dark matter halos or baryonic disks. Only two parameters are required to determine the global “system” properties of galaxies, RVirial and VEsc (where RVirial ≈ RDisk):

Although Newtonian dynamics ensures VCirc = VEsc/√2 at RVirial, it cannot provide a defined VCirc profile as a function of radius. ɅCDM theory removes this difficulty as dark matter halos conveniently have constant VCirc (flat virial halo rotation) permitting use of Newton’s second law for circular motion:

This substitution effectively decouples galactic VEsc (and by association, VCirc) from any significant baryonic influence. Rather than constraining VCirc to a theoretical value throughout RVirial, we treat VCirc as a direct observable, now possible with the availability of accurate component velocities measured between 6 and 30 kpc, a region spanning both the inner and outer Galactic disks (King III 2015). In following sections we explore the implications of non-flat circular velocity against the latest, most sophisticated ɅCDM simulations. We leverage the unbound Hyper Velocity Star (HVS) population to estimate Galactic VEsc and provide a physical (classical) origin and explanation for King’s recently discovered kinematic feature.

Galactic Rotation – RC-SP versus ɅCDM Models

In this section, we examine a well-cited rotation curve from Bhattacharjee augmented with data from Bajkova and Bobylev. This composite rotation curve is reproduced in Figure 2. Included are three ɅCDM models, labeled 1, 2, and 3 from the original figure (Bajkova 2016). The RC-SP Galactic rotation curve fit is shown by the black dash. This curve fit is based on observed velocities within the disk (R≤40 kpc) and the Keplerian decline beyond. Into the original figure, we have inserted a model of King’s velocity peak (gray dash) where observations are entirely missing. As shown below, all three ɅCDM rotation curves smoothly span this range, perhaps unaware of the recent discovery of this kinematic feature.

Figure 1: Galactic rotation curve from Bhattacharjee, Bobylev and Bajkova. Panel (a) illustrates the Obs/RC-SP integrated rotation profile (black dash). Prominent kinematic features include a spike in rotation to escape velocity (gray dashed) and a Keplerian decline beyond the Galactic disk. ɅCDM models are identified 1, 2 and 3 (see Bajkova ref. for details). Panel (b) is a semi-log version of panel (a). This view provides the specific R-V combination (open circle) corresponding to the RC-SP Galactic dynamic mass. This mass is obtained from the fit (Keplerian decline) beyond 40 kpc. The average rotation velocity (blue dash) and the Sun’s velocity and position (red dot) is provided for reference. Image source – Fig.1 (Bajkova 2016)

In the above figure, all three ɅCDM models conflict with the Galaxy data especially at outer radii where it is evident that a conventional Keplerian decline provides a better fit than the flat rotation prediction (Sofue 2015) (Huang 2016). The gray dashed curve in the region of Bajkova’s missing data is based on a recently measured kinematical feature in the inner stellar halo. We challenge the notion this feature is a perturbative/ transient halo artifact, contending it is long-lived and thermodynamic in origin.

Observed Kinematics of the Galactic Stellar Halo

The kinematics of the stellar halo serves as a sensitive probe of Galactic dynamics especially within ɅCDM cosmology where the stellar and dark matter halos share the same space. We examine in detail the properties of the stellar halo recently obtained from King’s high precision component velocity dispersion survey spanning 6 to 30 kpc.

Our focus is a previously identified kinematic feature termed the “tangential dip.” New observations have resulted in this dip becoming a significant trough, creating more tension between galactic kinematics and ɅCDM and MOND model expectations. The main take-away is that this this stellar feature should not be discounted or ignored within any truly accurate model of the Galaxy.

Figure 2 shows King’s component velocity dispersion results (R, φ, and θ). We have annotated King’s data with isotropic dispersion (black dash) corresponding to a central baryonic mass MBar = 0.085x1012Mʘ and average Galactic rotation velocity (blue dash). The lower right-hand panel provides the anisotropy coefficient β, demonstrating the extremity of this “dip” that translates to high velocity circular orbits.

Figure 2: a., b., and c. panels provide observed velocity dispersions for radial, polar and azimuthal components. Included with the original data is the MW average disk rotation velocity VCirc = 230 km/s (blue dash) and the isotropic baryonic ‘point mass’ curve (black dash). In the lower right panel, Extreme depression in anisotropy coefficient, β = [1-(σ θ 2 φ 2)/2σ R 2] denotes a spike in highly tangential orbits near mid-disk. To date, this kinematic substructure has not been incorporated into ɅCDM models of the Galaxy. Image source - Fig. 10 (King III 2015)

The dynamic mass distribution within the Galaxy is roughly traced as the net positive difference in velocity between the baryonic isotropic curve (black dash) and particular dispersion components. We find very little velocity support from the radial, with azimuthal and vertical components being dominant. This particular kinematic substructure cannot be reproduced within the context of the theoretical properties of dark matter halos.

In the next section, we next construct a dynamical model that relies on this complex but subtle kinematic substructure. We combine King’s results with those obtained from Hyper Velocity Star (HVS) surveys to construct a physically consistent (kinematical, dynamical, and thermodynamical) model of the Galaxy. Note we emphasis a single RC-SP Galaxy model based on observation rather than simulation ad hoc “best fits.”

Hyper Velocity Star Orbital Parameters

In this section, we advance a thermodynamic origin for the observed population distribution of Hypervelocity Stars (HVS). This unbound stellar population of stars is receiving attention as a method to quantify the dynamic (or dark matter halo) mass of the Galaxy. Currently, the HVS population is thought to acquire extreme velocities through intense gravitational ejection mechanisms deep inside the Galactic core (Tauris 2015) (Fragione 2016a) (Rossi 2016) (Fragione 2016b). Dynamic masses obtained by ejection mechanisms rely on “chance” encounters (constrained to the Galactic center) and complex three-body gravitational interactions to obtain a HVS model population. Figure 3 shows the HVS sample space against four dark matter halo mass models (see inset). In this figure, Fragione regarded stars VObs >275 km/s as “unbound” and a HVS candidate. To this original figure, we include King’s summed component velocities in quadrature (blue dash) with its measured peak velocity of 432 km/s at 23 kpc. Just beyond the peak, we find Galactic velocities plunge into a conventional Keplerian decline beyond the baryonic disk equivalent to MDyn = 0.5 x1012Mʘ and not the excessively high dark matter halo masses depicted.

Figure 3: Hypervelocity star sample population R-V space for the MW compared against RC-SP radial escape velocity profile (bold blue dotted) and  four dark matter halo models from Fragione (various black - see inset for details). King III quadrature summed velocity dispersion components (blue dashed spike).As the MW is a virialized system, we assign the MW’s 432 km/s peak as the Maxwell-Boltzmann “most probable” velocity. RC-SP parameterized average disk and Keplerian decline (both black dashed) as labeled. Image source – Fig 2 (Fragione 2016a)

The above figure highlights an issue which has been plaguing ɅCDM cosmology since inception, the tremendous insensitivity between halo properties and observation. Due to this lack of connection and high uncertainty between the dark matter and baryonic constituents, a particular halo model with a halo mass between 1.2 – 1.7 x1012Mʘ could only be “favored” over the others. In effect, halo mass uncertainty is equivalent to the RC-SP dynamic mass of the entire Galaxy. From Figure 3, we certainly observe a link between VPeak at 23 kpc as the virial radius of origin for the HVS population. As such, the HVS sample should be distributed based on thermodynamic considerations. The next step is to assign a physical mechanism responsible for this particular profile for the HVS sample population – the Maxwell-Boltzmann probability distribution.

A Thermodynamic Solution to Explain HVS Sample/Population Statistics

In this section we focus on the Maxwell-Boltzmann form and define peak velocity equivalent to the most probable velocity Vmp = 432 km/s within the distribution with (Wu 2014):

Figure 5 below compares Rossi’s expectation for the HVS velocity distributions (red dash and solid black) based on the gravitational ejection model. As it appears, the dearth of data beyond the peak indicates a narrow distribution at high velocities directly attributable to the very deep gravitational well of the dark matter halo. The M-B distribution (blue dash) shows a significant high velocity tail should be present.

Figure 4: HVS probability density functions per black hole ejection schemes (solid black and red dash). Observed “unbound” late B-star (~3.5 Mʘ) HVS distribution (gray histograms) from Brown 2014. The alternative RC-SP/Maxwell-Boltzmann probability distribution (blue dash) is based on most probable velocity Vmp=432 km/s (black dash) with arbitrary vertical scale. Image source - Fig. 3 (Rossi 2016)

As shown above, Rossi contends that the linear decline in HVS distribution in the low velocity tail is expected. We find the M-B distribution (blue dash) in this region is actually more linear than either ejection model, but no true discrimination between models is possible < 432 kms-1.

At the high velocity end of the distribution, Rossi contends the steep decline in the model is real with HVS becoming increasing rare at higher velocities. We contend it is the high velocity tail of the HVS distribution that distinguishes the M-B solution over chance gravitational encounters. In Figure 6 below, we expand the HVS sample distribution beyond ≈400 km/s by including lower stellar mass G/K dwarfs (Tauris 2015). In the figure below, the absence of HVS < 350 km/s is due to arbitrary truncation of the data.

Figure 5: Velocity distribution of HVS late B-type (blue histograms) and G/K-dwarfs (green histograms) with respect to the Galactic rest frame. The combined data traces the M-B probability distribution (blue dash) for a most probable velocity, Vmp= 432 km/s. The M-B vertical scale is matched to the data for comparative purposes. The M-B average velocity is VAvg = 487 km/s and root-mean-squared Vrms=529 km/s. This latter value is near the oft cited global Galactic escape velocity of 550 km/s. Image source - Fig. 1 (Tauris 2015)

We find HVS stellar mass tracks well with the M-B solution, with lighter G/K dwarfs exhibit greater net velocity than heavier late B-stars, accurately tracing the overall M-B distribution and the magnitude of escape velocities in relation to King’s velocity dispersion results. Of course, the stellar universal Initial Mass Function (IMF) needs to be considered as it directly influences the mass of the star that could become a HVS candidate, i.e., the IMF exhibits a peak in stellar mass between 0.2Mʘ to 4Mʘ (Baldry 2003) (Offner 2014). We would expect the M-B distribution would become fully “occupied” but this (to date) is not the case. Either the missing data is due to severe under-sampling or a more subtle effect not yet fully understood.

This thermodynamic interpretation is physically consistent with a virialized, quasi- equilibrium system in highly ordered motion (Struck 2016). Struck interpreted this ordered motion as “free energy” that will be thermalized in the future. We contend ordered galactic motion is “potential energy” contributing to total dynamic mass today.


Recent data suggests the Galaxy is an open thermodynamic “system” in quasi-equilibrium with its external “surroundings.” Under this model, we employ observed kinematics of the inner stellar halo and the latest sampling of Hyper Velocity Stars to indicate a thermodynamic origin without imposing any deviation from classical mechanics. Thanks to the Winnower for open access publishing, the research community and private communication that make this information available to the larger audience. This paper is dedicated to my dad.

Appendix – King’s Stellar Inner Halo Velocity Dispersion


Bajkova, A.T., Bobylev, V.V. "Rotation Curve and Mass Distribution in the Galaxy from the Velocities of Objects at Distances up to 200 kpc." arXiv, Jul 27, 2016:


Fragione, G., Capuzzo-Dolcetta, R., Kroupa, P. "Hypervelocity stars from young stellar clusters in the Galactic Centre." arXiv, Sep 17, 2016b:

Fragione, G., Loeb, A. "Constraining Milky Way mass with Hypervelocity Stars." arXiv, Aug 4, 2016a:

Huang, Y., Liu, X.-W., Yuan, M.-S. et al. "The Milky Way’s rotation curve out to 100 kpc and its constraint on." arXiv, Apr 5, 2016:

King III, C., Brown W.R., Geller, M.J., Kenyon, S.J. "STELLAR VELOCITY DISPERION AND ANISOTROPY OF THE MILKY WAY INNER HALO." The Astrophysical Journal, 813:89, Nov 10, 2015:

La Fortune, J.M. "The Baryonic Tully-Fisher Relationship: Empirical Evidence for a Baryonic Solution to the Galactic “Missing Mass” Problem." The Winnower, Jul 15, 2016:

Offner, S.S.R, Clark, P.C., Hennebelle, P., et al. "The Origin and Universality of the Stellar Initial Mass Function." arXiv, May 19, 2014:

Peebles, P.J.E. "Rotation of Galaxies and the Gravitational Instability Picture." A&A, 11, 377-386, Dec 1970 rec'd 23, 1971:

Rossi, E.M., Marchetti, T., Cacciato, M., Kuiack, K., Sari, R. "Joint constraints on the Galactic dark matter halo and Galactic Centre from hypervelocity stars." arXiv, Aug 5, 2016:

Sofue, Y. "Dark Halos of M31 and the Milky Way." arXiv, Apr 12, 2015:

Struck, C., Elmegren, B.G. "Near-exponential surface densities as hydrostatic, nonequilibrium profiles in galaxy discs." arXiv, Sep 28, 2016:

Tauris, T.M. "Maximum speed of hypervelocity stars ejected from binaries." MNRAS 448, 1 L6-L10, Mar 21, 2015:

Wu, Y-C., Y. "An Exploration of the Limits of the Maxwell-Boltzmann Distribution." arXiv, Oct 25, 2014:


This article and its reviews are distributed under the terms of the Creative Commons Attribution 4.0 International License, which permits unrestricted use, distribution, and redistribution in any medium, provided that the original author and source are credited.