Next Article in Journal
Origins and Natures of Inflation, Dark Matter and Dark Energy
Previous Article in Journal
On the Apparent Discretization of Spacetime and Its Connection with the Cosmological Constant
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Relevance of Dynamical Friction for the MW/LMC/SMC Triple System

1
scdsoft AG, Albert-Nestler-Str. 21, 76131 Karlsruhe, Germany
2
Helmholtz-Institut für Strahlen und Kernphysik (HISKP), Nussallee 14-16, 53115 Bonn, Germany
3
Astronomical Institute, Faculty of Mathematics and Physics, Charles University in Prague, V Holešovičkách 2, 18000 Prague, Czech Republic
*
Authors to whom correspondence should be addressed.
Universe 2024, 10(3), 143; https://doi.org/10.3390/universe10030143
Submission received: 4 December 2023 / Revised: 25 January 2024 / Accepted: 7 March 2024 / Published: 14 March 2024

Abstract

:
Simulations of structure formation in the standard cold dark matter cosmological model quantify the dark matter halos of galaxies. Taking into account dynamical friction between dark matter halos, we investigate the past orbital dynamical evolution of the Magellanic Clouds in the presence of the Galaxy. Our calculations are based on a three-body model of rigid Navarro–Frenk–White profiles for dark matter halos but were verified in a previous publication by comparison to high-resolution N-body simulations of live self-consistent systems. Under the requirement that the LMC and SMC had an encounter within 20 kpc between 1 and 4 Gyr ago in order to allow the development of the Magellanic Stream, using the latest astrometric data, the dynamical evolution of the MW/LMC/SMC system is calculated backwards in time. With the employment of the genetic algorithm and a Markov-Chain Monte-Carlo method, the present state of this system is unlikely, with a probability of < 10 9 ( 6 σ complement), because the solutions found do not fit into the error bars for the observed plane-of-sky velocity components of the Magellanic Clouds. This implies that orbital solutions that assume dark matter halos, according to cosmological structure formation theory, to exist around the Magellanic Clouds and the Milky Way are not possible with a confidence of more than 6 sigma.

1. Introduction

The standard ΛCDM cosmological model [1,2,3] requires about 81 per cent of the matter content of the Universe to comprise pressureless cold dark matter (CDM) particles, which are not described by the standard model of physics, and about 68 per cent of the total mass–energy density of the Universe to be composed of dark energy, represented by the cosmological constant Λ [4,5]. The related Λ WDM model is based on dark matter (DM) being made of warm dark matter particles that have a smaller mass than CDM particles but lead to very similar DM halos within which the observed galaxies reside (e.g., [6,7,8,9]). Detecting DM particles has become a major worldwide effort, which has so far yielded null results (e.g., [10,11,12,13,14,15] and references therein). Finding evidence for DM particles under the assumption that they have a finite interaction cross-section with standard-model particles may lead to null detection if the interaction cross-section is too small to be measurable. However, arguments based on the observed strong correlations between standard particles and dark matter suggest a significant cross-section [16]. This indicates an impasse in the experimental verification of the existence of DM particles.
An alternative and robust method for establishing the existence of DM particles is to develop techniques that rely solely on their gravitational interaction with ordinary matter. A suitable approach is given by Chandrasekhar dynamical friction (e.g., [5]): when a satellite galaxy with its own DM halo enters the DM halo of a host galaxy, its orbit decays as a result of dynamical friction, and the satellite merges with the host. The decay of the orbit does not depend on the mass of the DM particle but only on the mass density of DM particles, which is fixed by the cosmological parameters. This is the primary reason for the existence of the merger tree in standard cosmology (e.g., [17,18]) and for interacting galaxies to be thought of as merging galaxies. Effectively, a DM halo, being 10–20 times more extended and 50–100 times as massive as the observable part of a galaxy, works like a spider’s web. In contrast, if there were to be no DM halos around galaxies, then galaxy–galaxy encounters would be significantly less dissipative, galaxies would encounter each other multiple times, and mergers would be rare. In the eventuality that DM halos were not to exist, however, non-relativistic gravitational theory would need to be non-Newtonian, and given the correlations that galaxies are observed to obey, this theory would need to be effectively Milgromian [19,20,21,22]. Galaxy–galaxy encounters and mutual orbits in Milgromian dynamics and without DM halos have been studied (for example, the interacting Antennae galaxy pair [23], and the Milky Way–Andromeda binary explaining the mutually correlated planes of satellites around both hosts [24,25,26]). The law of gravitation fundamentally defines the formation and evolution of galaxies. For example, by mergers being much rarer in Milgromian gravitation (i.e., without DM), galaxies evolve largely in isolation [27,28,29], and the formation of cosmological structures proceeds differently than in the DM-based models [5,30,31]. It is thus of paramount importance to test for the existence of DM halos around galaxies.
The test based on Chandrasekhar dynamical friction for the existence of DM halos was introduced by Angus, Diaferio and Kroupa [32] by addressing the question of whether the present-day Galactocentric distances and motion vectors of observed satellite galaxies of the Milky Way (MW) conform to their putative infall many Gyr ago. The solutions for those satellite galaxies for which proper motions were available imply tension with the existence of DM halos since no infall solutions were found. In contrast, without the DM component, the satellite galaxies would be orbiting about the MW, having been most likely born as tidal dwarf galaxies during the MW–Andromeda encounter about 10 Gyr ago [25,26,33]. A further test for the presence of the DM component applying dynamical friction was achieved by Roshan et al. [34], who found the observed bars of disk galaxies to be too long and rotating too fast in comparison with the theoretically expected bars in the presence of DM halos that absorb the bar’s angular momentum. The reported discrepancy amounts to significantly higher than the 5-sigma threshold, such that the observations are incompatible with the existence of DM halos. Another independent application of the dynamical friction test is available on the basis of the observed distribution of matter within the M81 group of galaxies. The extended and connected tidal material implies multiple past close encounters of the group members. But the dynamical evolution of the M81 group of galaxies is difficult to understand theoretically if the galaxies are contained in the DM halos that are expected in standard cosmology ([35] and references therein). The problem is that, in the presence of DM halos, the group merges too rapidly to allow the tidal material to be dispersed as observed.
The interesting aspect of these results is that the three independent analyses of the orbital dynamics of the MW satellite galaxies, galactic bars and the M81 group agree in the conclusion that the data are difficult if not impossible to understand in the presence of DM halos. Because the implications of the non-existence of DM bears major implications for theoretical cosmology and galaxy formation and evolution, further tests are important to ascertain the above conclusions, or indeed to question them.
Given the high-accuracy and high-precision position and velocity data available today through the astrometric Gaia mission, the problem of how the existence of the massive and extended DM halos in which galaxies reside can be tested for is revisited here using the process of Chandrasekhar dynamical friction on the triple-galaxy system comprising the MW and the Large (LMC) and Small (SMC) Magellanic Clouds, including the Magellanic Stream. This system is comparable to the M81 system, as it too has a large gaseous structure stemming from the Magellanic Clouds, the Magellanic Stream, which constrains the past orbits of the components of the system. In the following, we assume that the standard cosmological model is correct and that the MW, the LMC and the SMC are each contained in DM halos that are consistent with those obtained from the theory of structure formation in the Λ CMD model of cosmology. That is, we associate the baryonic mass of each galaxy with a DM halo profile as predicted by the theory in order to test the theory.
Given that the LMC and SMC are about 50 and 60 kpc distant from the centre of the MW and about 20 kpc distant from each other and that the radii of the DM halos are such that all three galaxies are immersed in the dark matter halo of the MW and that the SMC is immersed in the dark matter halo of the LMC and vice versa, it is apparent that Chandrasekhar dynamical friction is likely to play a very significant role in establishing the orbits of the three galaxies relative to each other. Recent work [36,37,38,39,40,41] suggests the Magellanic Clouds to be on their first pericentre passage such that the bulk properties of their putative DM halos will not be significantly stripped. Detailed dynamical modelling of the Magellanic Clouds problem has, until now, not come to the conclusion that there is a problem.
For example, Besla et al. [37] performed simulations of the LMC and SMC in the first-infall scenario by assuming that dynamical friction from the MW DM halo can be neglected. This is not correct and significantly helps the LMC/SMC pair to exist longer before merging. Indeed, Lucchini et al. [42] studied the formation of the Magellanic Stream using hydrodynamical simulations in the first-infall scenario for the LMC and SMC. Their results, obtained in fully live DM halos of all components such that dynamical friction is correctly computed self-consistently, show that the LMC/SMC binary orbit shrinks significantly more rapidly (their Figure 1) than calculated by Besla et al. [37]. But Lucchini et al. adopted DM halo masses of the LMC ( 1.8 × 10 11 M ), of the SMC ( 1.9 × 10 10 M ) and of the MW ( 1.0 × 10 12 M ) that are significantly less massive than predicted by the Λ CDM model (see Table 1 below). Their calculations also do not consider whether the infalling LMC/SMC binary would have survived for sufficiently long before infall, because it is implausible for the present-day LMC/SMC binary to have formed during or shortly before infall, and whether the infall velocity vector is consistent with the Hubble flow 3.46 Gyr ago at the start of their simulation. Vasiliev [41] studied a scenario in which the LMC is on its second passage past the MW but ignored the SMC. The calculated orbits of the LMC decay due to dynamical friction consistently to the results presented here. Thus, this past work either neglects important contributions to dynamical friction or members in the problem or uses light-weight DM halos that are not consistent with the Λ CDM cosmological model. Essentially, past work has demonstrated that the MW/LMC/SMC triple system can broadly be understood in the context of orbital dynamics by effectively suppressing Chandrasekhar dynamical friction, but it does not demonstrate that the calculated orbits are consistent with the standard DM-based cosmological theory.
In order to assess whether the standard DM-based models can account for the existence of the observed MW/LMC/SMC triple system, the present work studies their orbital history by strictly constraining the DM halo of each of the three involved galaxies to be consistent with the Λ CDM-allowed DM halo masses. The dark matter component is thus assumed to be pressureless dust. The calculations furthermore impose the condition that the LMC and SMC had an encounter within 20 kpc of each other between 1 and 4 Gyr ago in order to allow the Magellanic Stream to form, and they take into account cosmological expansion to constrain the relative positions of the galaxies 5 and more Gyr ago. To check for the consistency of the solutions, two independent statistical methods were applied to search for all allowed initial conditions within the error bars for the transverse velocity components of the LMC and SMC. The two tests yielded indistinguishable results. The calculations based on rigid Navarro–Frenk–White (NFW) DM halo profiles are conservative because equivalent simulations with live DM halos lead to faster merging [35]. These calculations show there to be no orbital solutions and that the observed system comprising the MW/LMC/SMC plus Magellanic Stream cannot be understood in the presence of DM halos in the context of the Λ CDM model. The problem that gravitational theory leads to the observed system is left for a future contribution.

2. The Model

2.1. NFW Profiles

When testing a theory for consistency with data, it is of paramount importance to not mix the two: purely theoretically calculated properties need to be compared with empirical data that have not been modulated by the very theory to be tested, and vice versa. It is therefore not permissible to choose sub-massive DM halos that allow a solution in order to argue that Λ CDM or Λ WDM theory accounts for the MW/LMC/SMC system. Thus, here, we are adamant at requiring the use of the theoretically predicted NFW dark matter halo profiles [43]. These arise from self-consistent cosmological structure formation simulations such as those documented in [44] for the range of baryonic galaxy masses considered here. More generally, the NFW DM halo profile is a standardised outcome of structure formation simulations in Λ CDM theory (e.g., [45]), and cold or warm DM halos have been shown to have similar overall profiles and densities [6,7,8,9]. While the inner and outer DM halo profiles can differ from the NFW form, the latter well represents the distribution of theoretically predicted DM particles around the baryonic component of a galaxy, especially for galaxies that are on a first-infall orbit (e.g., [46,47,48,49]). In seeking orbital solutions, the stellar masses are allowed to vary within ± 30 per cent of their nominal values with the correspondingly different DM halo masses (Table 1 below). The range of theoretically predicted DM halo masses ( 41 to + 68 per cent for the MW, 42 to + 14 per cent for the LMC and 17 to + 16 per cent for the SMC) is thus accounted for.
The DM halo of either galaxy is treated as a rigid halo with a density profile according to [43] (NFW profile), truncated at the radius R 200 :
ρ ( r ) = ρ 0 r / R s 1 + r / R s 2 ,
where R s = R 200 / c , with R 200 denoting the radius yielding an average density of the halo of 200 times the cosmological critical density,
ρ c r i t = 3 H 2 8 π G ,
and the concentration parameter c [50],
log 10 c = 1.02 0.109 log 10 M vir 10 12 M .
The DM halo masses are derived from the stellar masses of the galaxies by means of Figure 7 of [44]. Note that using rigid DM halo profiles in such orbital computations is admissible because they have been verified to give conservative solutions (slower merging times) than self-consistent simulations with live DM halos [35] and, in particular, because the MW/LMC/SMC system has only a recent encounter history.

2.2. Dynamical Friction and Equations of Motion

Exploring the dynamics of bodies orbiting in the interior of DM halos implies that dynamical friction1 has to be taken into account in an appropriate manner [51]. Here, the formulation of Chandrasekhar’s dynamical friction as used in [35], their Equations (1) and (2), is applied.
Chandrasekhar’s formula gives a quick-to-compute estimate for orbital decay, which is needed for the sake of establishing statistical statements about merger rates between galaxies. High-resolution simulations of live self-consistent systems confirm our approach of employing this semi-analytical formula in our three-body calculations (see Section 7 as well as Figures 13 and 14 in [35]). Using computationally significantly more time-intensive live simulations, in fact, leads to more rapid orbital decay such that our here-employed semi-analytical approach is conservative by allowing a larger range of solutions than would be the case if live DM halos were to be used.
The equations of motion for the individual galaxies are as given in Appendix C of [35]2 and are here augmented with an additional term taking into account the Hubble flow as follows: Based on the assumption of a flat universe (curvature parameter k = 0 ), we extended the equations of motion by the cosmic acceleration term s ¨ s r i caused by the Hubble flow of the expanding universe. Here, s ( t ) is the scale factor of the Universe, and r i is the position of a galaxy in the centre-of-mass frame of the group. In Cartesian coordinates, the equations of motion are then given by ( k = 1 , 2 , 3 for the MW, LMC and SMC, respectively)
d 2 d t 2 ( r i ) k = m i · s ¨ s ( r i ) k + ( F i ) k ,
with
( F i ) k = j i ( r i ) k V i j + ( F i j D F ) k ( F j i D F ) k .
Here, F i is the total force acting on galaxy i, V i j is the potential energy between galaxies i and j, and F i j D F is the dynamical friction force acting on galaxy i caused by the overlap of the DM halos of galaxies i and j (according to actio est reactio, F j i D F is taken into account, too).
Assuming a flat dark energy-dominated cosmology ( Ω m , 0 + Ω Λ , 0 = 1 ), the second Friedmann equation becomes
s ¨ s = H 0 2 ( Ω Λ , 0 1 2 Ω m , 0 s 3 ) ,
where Ω m , 0 and Ω Λ , 0 are the matter and dark energy densities at the present time, respectively, scaled by ρ c r i t , and H 0 is the Hubble constant. Setting the conditions s ( 0 ) = 0 and s ˙ = H 0 at s ( t ) = 1 (present time) yields an analytical expression for the cosmic-scale factor:
s ( t ) = ( Ω m , 0 Ω Λ , 0 ) 1 3 sinh 2 3 ( 3 2 Ω Λ , 0 H 0 t ) ,
where t is the age of the Universe. For our calculations, we used the following values: Ω m , 0 = 0.315 , Ω Λ , 0 = 0.685 and H 0 = 67.3 km / s / Mpc . However, the cosmic acceleration term does not play a significant role for the time frame considered here. The order of magnitude contribution, compared to the forces between the DM halos in Equation (4), is in the range from 10 3 to 10 2 within the time range [ 5 Gyr, today]. All numerical computations were performed using SAP’s ABAP development workbench.

3. Observational Data

For the sake of easy accessibility, we list the observational constraints in Table 1, Table 2, Table 3, Table 4, Table 5, Table 6 and Table 7, and the references regarding the basic observational data in Table 2.
Table 1. Stellar masses of the galaxies (model (o)), varied by − 30 % (model (m)) and + 30 % (model (p)), and the derived DM halo masses according to Section 2.1.
Table 1. Stellar masses of the galaxies (model (o)), varied by − 30 % (model (m)) and + 30 % (model (p)), and the derived DM halo masses according to Section 2.1.
ObjectModelStellar MassDM Halo Mass
[ M ][ M ]
(o) 5 × 10 10 2.41 × 10 12
MW 30 % (m) 3.5 × 10 10 1.39 × 10 12
+ 30 % (p) 6.5 × 10 10 4.05 × 10 12
(o) 3.2 × 10 9 2.55 × 10 11
LMC 30 % (m) 2.24 × 10 9 1.47 × 10 11
+ 30 % (p) 4.16 × 10 9 2.90 × 10 11
(o) 5.3 × 10 8 1.07 × 10 11
SMC 30 % (m) 3.71 × 10 8 8.86 × 10 10
+ 30 % (p) 6.89 × 10 8 1.24 × 10 11
As mentioned in Section 2.1, we derived the DM halo masses of the MW, LMC and SMC according to [44] based on the stellar masses extracted from the references given in Table 2. Our approach, therefore, is in full correspondence with the predictions of ΛCDM regarding the DM halo masses. However, in order to overcome possible uncertainties in that regard and to achieve more independent results, we decided to apply our statistical evaluations in Section 4 individually to 27 mass combinations by varying the stellar mass of each galaxy by plus/minus 30 % , as displayed in Table 1. The 27 mass combinations are then denoted by the indicators “o” for the original masses, “m” for 30 % and “p” for + 30 % , for instance, “o-m-p” for the MW/LMC/SMC mass triple.
The relevant coordinates and velocities are presented in Table 3, Table 4 and Table 5, and the transformation to the Cartesian heliocentric equatorial system is shown in Table 6 and Table 7.
The Magellanic Stream consists of H I gas, which trails behind both the LMC and the SMC across a large fraction of the sky. It is thought to be the result of a combination of tidal forces and ram pressure stripping through the orbit of the LMC and SMC within the hot gaseous halo of the MW. Studies of the origin of the Magellanic Stream have shown that it was most likely created when the LMC and the SMC had a close encounter, during which some of the gas of the LMC and SMC became less bound to be subsequently removed from the pair through ram pressure stripping [39,40,62,63].

4. Statistical Methods

In order to cross-check the solutions and to enhance the confidence in them, we utilised two independent statistical methods, a Markov-Chain Monte-Carlo method (MCMC, see Section 4.1) and the genetic algorithm (GA, see Section 4.2). The reason why we do so is given in the introductory statement of Section 4.2.
This is the initial situation for each of the MCMC and GA methods: due to the error bars of the transverse velocity components in RA and DEC for the Magellanic Clouds (see Table 4), we have to consider four open parameters P i regarding the calculation of the three-body orbits of the dark matter halos, as displayed in Table 8.
This means that by employing the MCMC method and the GA method separately for each of the 27 mass combinations (see Section 3), we search for solutions of the three-galaxy orbits with best fits to the transverse velocity components of the Magellanic Clouds. It is important to note that our goal is to achieve results independent of the particular masses of the galaxies, not to find a best-fit mass combination.
Further, based on the radio-astronomical observational data concerning the Magellanic H I -Stream explained in Section 3, we specify the following broad condition, which we hereinafter refer to as the condition COND, regarding admissible past orbits of the Magellanic Clouds: the LMC and SMC encountered each other within the past time interval of [ 4 Gyr, 1 Gyr] at a pericentre distance of less than 20 kpc. Incorporating this condition, the algorithms searched for solutions by integrating Equation (4) backwards in time up to 5 Gyr.

4.1. Markov-Chain Monte Carlo (MCMC)

We followed a methodology proposed by [64] employing an affine-invariant ensemble sampler for the Markov-Chain Monte-Carlo method. A detailed guideline for implementation can be found in [65]. The basics of applying this formalism to our situation are outlined in Appendix D of [35].

4.1.1. Definition of the Posterior Probability Density

First of all, concerning the open parameters, we need to account for the error bars of the transverse velocity components v i of the Magellanic Clouds. This is ensured by an appropriate definition of the prior distribution p ( X ) , where X symbolises the parameter vector P 1 P 4 . With the values from Table 8, the four contributions to the prior distribution read ( i = 1 , , 4 )
p ( X ) i exp ( ( P i v i ) 2 2 σ i 2 ) .
Exploiting the minimal distance d 23 between the LMC and SMC within the time period [ 4 Gyr, 1 Gyr], the condition COND implies the likelihood function
P C ( X ) 1 , d 23 d p e r , exp d 23 d p e r 2 2 · d 0 2 , d 23 > d p e r ,
with d p e r = 20 kpc and d 0 = 5 kpc. The posterior probability density is then given by the product of Equations (8) and (9):
π ( X ) i = 1 4 p ( X ) i · P C ( X ) .
Here, the normalising constant for an overall probability of 1 is neglected because it is clear from the comparative search algorithm that the absolute values of π ( X ) are irrelevant.

4.1.2. The First Ensemble

Generating random integer numbers k { 0 , , 1000 } separately for each walker and for each open parameter of a given walker, the first ensembles are created according to ( i { 1 , , 4 } ):
P i = v i n σ i + k · 2 n σ i / 1000 ,
with varying n to account for extended error bars, as explained in Section 4.3.

4.2. Genetic Algorithm (GA)

The general aspects of the genetic algorithm are explained in detail in [66]. As pointed out there, the major advantage of this method is perceived to be its capability of avoiding local maxima in the process of searching for the global maximum of a given function (here, the fitness function). This feature makes it worthwhile to employ the GA method as a second independent statistical approach besides the MCMC method.
A precise description of the algorithm, especially instructions for its implementation, can also be found in [67]. To put it in a nutshell, the open parameters are related to genes, which are concatenated to genotypes. The set of a number of genotypes is considered to be a population, where the members pairs produce a follow-up generation with new features due to randomly performed cross-over and mutation. Winners per generation are determined by means of the fitness function, and by comparing the generations, the overall winner is found.
Comparing it to the more familiar MCMC method, we have the following parallels: walker <-> genotype, ensemble <-> population, set of ensembles <-> generations, posterior probability density <-> fitness function.

4.2.1. The GA Generations

Each open parameter from Table 8 is mapped to a 4-digit string (“gene”) a b c d i ( i { 1 , , 4 } ),
P i = ( v i n σ i ) + a b c d i · 2 n σ i / 10000 ,
again with varying n, like in Section 4.1.2, in order to account for extended error bars, as explained in Section 4.3. All genes together define the genotypes as 4 × 4-digit strings,
a b c d 1 a b c d 4 ,
which generate a population of N p o p genotypes. The first generation is established by randomly creating N p o p 4 × 4-digit strings.

4.2.2. The Definition of the Fitness Function

Our goal is to establish GA evaluations that are compatible with the evaluations by means of the Markov-Chain Monte-Carlo method. Therefore, we chose the definition of the fitness function to be identical to the definition of the posterior probability density (Equation (10)),
F ( X ) = i = 1 4 p ( X ) i · P C ( X ) ,
where the individual components are taken from Equations (8) and (9).

4.3. Results

Both statistical methods, MCMC and GA, deliver mutually consistent, and in fact indistinguishable, results.

4.3.1. Step I

As a first step, we tried to find solutions for which all four open parameters P i lie within the 1 σ error bars of the transverse velocity components of the Magellanic Clouds (see Table 8). Employing the methods MCMC and GA as search engines, using 1000 ensembles with 100 walkers (MCMC) and, accordingly, 1000 generations with a population of 100 genotypes (GA), we repeated the search 100 times for both statistical methods for each of the 27 mass combinations according to Table 1. This attempt to find a solution failed. That is, within the here-probed 1 σ uncertainty range of the transverse velocity components, neither algorithm found orbits that fulfil the condition COND. In other words, the MW/LMC/SMC plus Magellanic Stream system cannot exist in its observed configuration in the presence of dark matter halos.
To proceed further, we gradually extended the allowed intervals for the parameters P i to be v i ± n σ i with increasing integers n; see Equations (11) and (12). Neither the MCMC nor the GA statistical method found a solution for σ = 1 , 2 , 3 . For σ = 4 , both methods delivered solutions for a single mass combination only, namely, p-m-p. The details regarding all 27 mass combinations are given in Appendix A.

4.3.2. Step II

How to go about quantifying the first results obtained in Section 4.3.1? First of all, for a mass combination with the first solutions based on n σ error intervals, we chose the ( n + 1 ) σ intervals to be allowed for the open parameters. This is motivated by the thought that relaxing the overall n σ condition for all open parameters may result in solutions with individual smaller-than- 1 σ deviations of the parameters.
Furthermore, we constructed a “probability-sigma-grid” in the following sense: If a parameter P i lies within the error bar ( < 1 σ i ), then its probability weighting factor is 1, a conservative cautious setting. If a parameter P i lies outside the error bar, then its deviation from v i is weighted with the remainder of the probability function based on 0.1 σ intervals. For instance, if we have P i = v i + 3.74 σ i , then the corresponding probability weighting factor for this parameter is calculated according to the remainder of the probability function for 3.7 σ .
Based on the same search method as in Section 4.3.1, we searched for best-fit solutions by calculating an overall probability for the combination of the open parameters. The results are displayed in Table 9 and Table 10 for the methods MCMC and GA, respectively.
Regarding the MCMC method, it is interesting to note that, in some cases, individual parameters are not confined to the mentioned ( n + 1 ) σ intervals because the stretch moves can place walkers outside of the intervals specified by Equation (11), thus supporting the overall process of maximising the posterior probability density consisting of all open parameters. This is not possible for the GA method, as the genotypes are a priori confined to the preset intervals.
The result in Section 4.3.1 that the mass combination p-m-p provides the best-fit solution is confirmed by either statistical method, MCMC and GA. Furthermore, Figure 1 shows the trajectories, displayed as pairwise distances, of the individual best-fit solutions for the mass combinations o-o-o and p-m-p. Even for that detail, both methods deliver identical results.

4.3.3. Step III

Especially due to the MCMC method not being utilised as a pure search engine only, we established its validity by checking the results obtained in the previous Section 4.3.2 in terms of the autocorrelation function for the mass combinations o-o-o (original masses) and p-m-p (mass combination with the overall best-fit solutions according to Section 4.3.1 and Section 4.3.2) and creating sets of 1000 solutions in each case3. For the GA method, we undertook broader search runs, establishing 1000 solutions for the mentioned mass combinations.
The autocorrelation functions are presented in Figure 2 and Figure 3, demonstrating that convergence is achieved quickly for the MCMC method, and the improved probabilities are displayed in Table 11 and Table 12 for the MCMC and GA methods, respectively.

4.3.4. Interpretation

In Step I (Section 4.3.1), it was shown that no orbital solutions are possible if the observational quantities are allowed to span the ± 1 σ uncertainty range. In Step II (Section 4.3.2), the uncertainty range was therefore increased to ± ( n + 1 ) σ , with the result that values of n are needed to obtain viable orbital solutions that lie outside the 5 σ confidence range of the quantities such that the orbital solutions are likely with probabilities smaller than 10 9 (Table 9 and Table 10). Step III (Section 4.3.3) demonstrates that the MCMC and GA methods yield indistinguishable results, thus confirming the conclusion that orbital solutions do not exist for the observed MW/LMC/SMC plus Magellanic Stream system in the presence of the theoretically expected DM halos.
The dynamical behaviour of the MW/LMC/SMC system demonstrates the importance of dynamical friction in the context of interacting galaxies, as dynamical friction significantly influences the orbits. Figure 4 shows that the forces due to dynamical friction between the LMC and the SMC are comparable to the pure gravitational force between the overlapping DM halos.
As a thought experiment, we calculated the orbits back in time to 7 Gyr with and without dynamical friction for the overall best-fit solution derived for the mass combination p-m-p and for the best-fit solution for the original mass combination o-o-o. Dynamical friction can be turned off by omitting the terms F D F in Equation (5). This retains the gravitational pull of the DM halo but avoids the deceleration through Chandrasekhar dynamical friction. This Gedanken experiment is thus a rough approximation of the situation in Milgromian dynamics, according to which a galaxy generates a Milgromian gravitational potential that can be viewed as stemming from a Newtonian plus phantom DM halo that does not generate dynamical friction. With this approximation of MOND, the GA method as the search engine immediately delivers 17 mass combinations already within the 1 σ intervals (see also Appendix A). In other words, solutions matching the error intervals for the transverse velocities of the LMC and SMC are found readily without Chandrasekhar’s dynamical friction but with the potential generated by the DM halo. The result is displayed in Figure 5. With the absence of Chandrasekhar dynamical friction on DM halos, the Magellanic Clouds would have a long orbital lifetime as satellites of the MW, possibly being massive tidal dwarf galaxies formed during the MW–Andromeda encounter about 10 Gyr ago.
These results thus suggest that the MW/LMC/SMC plus Magellanic Stream system may have a straightforward orbital solution in Milgromian dynamics, and it is thus of much interest to simulate this system and its possible origin in this non-Newtonian framework.

5. Conclusions

The conclusions reached by this analysis are very robust, since both the MCMC and GA algorithms lead to indistinguishable results. Taking the observed configuration in the six-dimensional phase space of the MW/LMC/SMC plus Magellanic Stream system as the necessary boundary condition, it is impossible to find orbital solutions backwards in time that fulfil the very liberal condition COND (Section 4). The orbits accelerate too rapidly (backwards in time) such that the LMC and SMC could not have remained bound long enough to have the required encounter that is needed to have occurred to produce the Magellanic Stream. In other words, the system merges too rapidly forward in time to allow a close encounter as defined through COND and to still be visible today as two distinctly separated galaxies next to the MW. Solutions do not even appear if the DM halo masses of the three galaxies are allowed to be up to 42 per cent larger or smaller than those under the standard assumptions (see Table 1). The possibility that the LMC and SMC fell into the MW DM halo independently of each other but at a similar time in order to allow them to pair up to the observed binary is arbitrarily unlikely because the LMC and SMC would have to have had relative velocities to each other and to the MW that oppose the Hubble flow. Such an unlikely solution is, in any case, not possible because the condition COND requires the LMC and SMC to have had a close encounter at a time in the past such that the two would have merged today due to the mutual dynamical friction on their respective DM halos. In any case, neither the MCMC nor the GA method found such solutions. This work thus shows that the observed configuration of the MW/LMC/SMC plus Magellanic Stream system is not possible in the presence of the theoretically expected DM halos.
The results based on the Chandrasekhar dynamical friction test applied to the MW/LMC/SMC triple system arrived at here corroborate the previous evidence based on the same test but other systems, noted in the Introduction, that questions the existence of dark matter particles. The independently documented problems [5,68,69,70,71,72,73,74] of fitting the standard model of cosmology to the observed Universe on most probed scales are consistent with these results.
The thought experiment in Section 4.3.4, in which the potentials of the DM halos are kept but the Chandrasekhar dynamical friction term is set to zero, naturally leads to solutions. This experiment is an approximation to the situation in Milgromian dynamics and demonstrates that the origin and evolution of the MW/LMC/SMC plus Magellanic Stream system need to be studied in this non-Newtonian framework.

Author Contributions

Conceptualisation, W.O. and P.K.; methodology, W.O. and P.K.; software, W.O.; validation, W.O. and P.K.; formal analysis, W.O.; resources, W.O.; data curation, W.O.; writing—original draft preparation, W.O. and P.K.; writing—review and editing, W.O. and P.K.; visualisation, W.O. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

All data underlying this research are available in this article (see Section 3).

Acknowledgments

W.O. acknowledges the support of scdsoft AG in providing an SAP system environment for the numerical calculations. Without the support of scdsoft’s executives P. Pfeifer and U. Temmer, the innovative approach of programming the numerical tasks in SAP’s language ABAP would not have been possible.

Conflicts of Interest

Wolfgang Oehm was an employee of scdsoft AG. The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
ABAPSAP’s programming language
CDMCold dark matter
CONDCondition specified in Section 4
DECDeclination
DMDark matter
GAGenetic algorithm
ΛCDMDark-energy plus cold-dark-matter model of cosmology
ΛWDMDark-energy plus warm-dark-matter model of cosmology
LMCLarge Magellanic Cloud
MCMCMarkov-Chain Monte Carlo
MWThe Galaxy (Milky Way)
NEDNASA/IPAC extragalactic database
NFWNavarro, Frenk and White profile
RARight ascension
SMCSmall Magellanic Cloud

Appendix A. First Search Results

Confining the solutions to n σ intervals simultaneously for all transverse velocity components of the Magellanic Clouds, and based on 100 attempts for each n σ interval to find an appropriate solution for each of the 27 mass combinations, we obtained the following results for the MCMC method:
1 σ : none;
2 σ : none;
3 σ : none;
4 σ : p-m-p;
5 σ : o-m-p, m-m-p, p-m-o;
6 σ : o-m-o, o-m-m, o-p-m, m-m-o, m-m-m, p-m-m;
7 σ : o-o-o, o-o-m, o-o-p, o-p-o, o-p-p, m-o-o, m-o-m,
m-o-p, m-p-o, m-p-m, m-p-p, p-o-o;
8 σ : p-o-m, p-o-p, p-p-o, p-p-m, p-p-p.
For the GA method, we obtained the following:
1 σ : none;
2 σ : none;
3 σ : none;
4 σ : p-m-p;
5 σ : o-m-p, o-p-m, m-m-p, p-m-o;
6 σ : o-o-m, o-m-o, o-m-m, o-p-o, m-o-p, m-m-o,
m-m-m, p-m-m;
7 σ : o-o-o, o-o-p, o-p-p, m-o-o, m-o-m, m-p-o, m-p-m,
m-p-p, p-o-o, p-o-m, p-o-p, p-p-o, p-p-p;
8 σ : p-p-m.
However, when neglecting the effect of dynamical friction (by omitting the terms F D F in Equation (5)), the GA method already delivers 17 mass combinations with appropriate solutions within the 1 σ intervals:
o-m-p, m-o-o, m-o-m, m-o-p, m-m-o, m-m-m, m-m-p, m-p-o, m-p-m, m-p-p, p-o-o, p-o-m, p-o-p, p-m-p, p-p-o, p-p-m, p-p-p.

Notes

1
The term “dynamical dissipation” appears to be the more appropriate than “dynamical friction” because the physical process is not friction but the distortion of the orbits of the individual particles of the dark matter halos due to the long-distance gravitational forces. However, as the term “dynamical friction” is established in the community, we retain this term.
2
There is a printing mistake in the equation for Φ i ( s i ) , the third equation in Appendix C in [35]: the − sign in the first row is incorrect.
3
As explained in [35] the autocorrelation functions deliver an estimate for the convergence of the search algorithm for sufficiently large numbers of ensembles. Therefore, in order to establish full convergence for the approximate calculation of the autocorrelation functions, we created 50,000 ensembles for both mass combinations considered, o-o-o and p-m-p. The mentioned 1000 solutions fulfilling the conditions ( n + 1 ) σ intervals and COND were extracted from follow-up runs based on the set of walkers of ensemble 50,000 in either mass combination case.

References

  1. Efstathiou, G.; Sutherland, W.J.; Maddox, S.J. The cosmological constant and cold dark matter. Nature 1990, 348, 705–707. [Google Scholar] [CrossRef]
  2. Ostriker, J.P.; Steinhardt, P.J. The observational case for a low-density Universe with a non-zero cosmological constant. Nature 1995, 377, 600–602. [Google Scholar] [CrossRef]
  3. Peebles, P.J.; Ratra, B. The cosmological constant and dark energy. Rev. Mod. Phys. 2003, 75, 559–606. [Google Scholar] [CrossRef]
  4. Planck, Collaboration; Aghanim, N.; Akrami, Y.; Ashdown, M.; Aumont, J.; Baccigalupi, C.; Ballardini, M.; Banday, A.J.; Barreiro, R.B.; Bartolo, N.; et al. Planck 2018 results. VI. Cosmological parameters. arXiv 2018, arXiv:1807.06209. [Google Scholar]
  5. Kroupa, P.; Gjergo, E.; Asencio, E.; Haslbauer, M.; Pflamm-Altenburg, J.; Wittenburg, N.; Samaras, N.; Thies, I.; Oehm, W. The many tensions with dark-matter based models and implications on the nature of the Universe. arXiv 2023, arXiv:2309.11552. [Google Scholar] [CrossRef]
  6. Bullock, J.S.; Kravtsov, V.A.; Colín, P. Angular Momentum Profiles of Warm Dark Matter Halos. Astrophys. J. Lett. 2002, 564, L1–L4. [Google Scholar] [CrossRef]
  7. Knebe, A.; Arnold, B.; Power, C.; Gibson, B.K. The dynamics of subhaloes in warm dark matter models. Mon. Not. R. Astron. Soc. 2008, 386, 1029–1037. [Google Scholar] [CrossRef]
  8. Paduroiu, S. Warm Dark Matter in Simulations. Universe 2022, 8, 76. [Google Scholar] [CrossRef]
  9. Rose, J.C.; Torrey, P.; Villaescusa-Navarro, F.; Vogelsberger, M.; O’Neil, S.; Medvedev, M.V.; Low, R.; Adhikari, R.; Anglés-Alcázar, D. Inferring warm dark matter masses with deep learning. Mon. Not. R. Astron. Soc. 2024, 527, 739–755. [Google Scholar] [CrossRef]
  10. Baudis, L. Dark matter searches. Ann. Der Phys. 2016, 528, 74–83. [Google Scholar] [CrossRef]
  11. Zitzer, B.; VERITAS Collaboration. A Search for Dark Matter from Dwarf Galaxies using VERITAS. In Proceedings of the 38th International Conference on High Energy Physics (ICHEP2016), Chicago, IL, USA, 3–10 August 2016; p. 446. [Google Scholar]
  12. Hoof, S.; Geringer-Sameth, A.; Trotta, R. A global analysis of dark matter signals from 27 dwarf spheroidal galaxies using 11 years of Fermi-LAT observations. J. Cosmol. Astropart. Phys. 2020, 2020, 012. [Google Scholar] [CrossRef]
  13. Ferreira, E.G.M. Ultra-Light Dark Matter. arXiv 2020, arXiv:2005.03254v1. [Google Scholar] [CrossRef]
  14. Heisig, J.; Korsmeier, M.; Winkler, M.W. Dark matter or correlated errors? Systematics of the AMS-02 antiproton excess. arXiv 2020, arXiv:2005.04237v1. [Google Scholar] [CrossRef]
  15. Landim, R.G. Constraining dark photon dark matter with fast radio bursts. arXiv 2020, arXiv:2005.08621. [Google Scholar]
  16. Salucci, P. The distribution of dark matter in galaxies. Astron. Astrophys. Rev. 2019, 27, 2. [Google Scholar] [CrossRef]
  17. Bullock, J.S.; Johnston, K.V. Tracing Galaxy Formation with Stellar Halos. I. Methods. Astrophys. J. 2005, 635, 931–949. [Google Scholar] [CrossRef]
  18. Roper, W.J.; Thomas, P.A.; Srisawat, C. MEGA: Merger graphs of structure formation. Mon. Not. R. Astron. Soc. 2020, 494, 4509–4524. [Google Scholar] [CrossRef]
  19. Milgrom, M. A modification of the Newtonian dynamics as a possible alternative to the hidden mass hypothesis. Astrophys. J. 1983, 270, 365–370. [Google Scholar] [CrossRef]
  20. Bekenstein, J.; Milgrom, M. Does the missing mass problem signal the breakdown of Newtonian gravity? Astrophys. J. 1984, 286, 7–14. [Google Scholar] [CrossRef]
  21. Famaey, B.; McGaugh, S.S. Modified Newtonian Dynamics (MOND): Observational Phenomenology and Relativistic Extensions. Living Rev. Relativ. 2012, 15, 10. [Google Scholar] [CrossRef]
  22. Banik, I.; Zhao, H. From Galactic Bars to the Hubble Tension: Weighing Up the Astrophysical Evidence for Milgromian Gravity. Symmetry 2022, 14, 1331. [Google Scholar] [CrossRef]
  23. Renaud, F.; Famaey, B.; Kroupa, P. Star formation triggered by galaxy interactions in modified gravity. Mon. Not. R. Astron. Soc. 2016, 463, 3637–3652. [Google Scholar] [CrossRef]
  24. Zhao, H.; Famaey, B.; Lüghausen, F.; Kroupa, P. Local Group timing in Milgromian dynamics. A past Milky Way-Andromeda encounter at z > 0.8. Astron. Astrophys. 2013, 557, L3. [Google Scholar] [CrossRef]
  25. Bílek, M.; Thies, I.; Kroupa, P.; Famaey, B. MOND simulation suggests an origin for some peculiarities in the Local Group. Astron. Astrophys. 2018, 614, A59. [Google Scholar] [CrossRef]
  26. Banik, I.; Thies, I.; Truelove, R.; Candlish, G.; Famaey, B.; Pawlowski, M.S.; Ibata, R.; Kroupa, P. 3D hydrodynamic simulations for the formation of the Local Group satellite planes. Mon. Not. R. Astron. Soc. 2022, 513, 129–158. [Google Scholar] [CrossRef]
  27. Wittenburg, N.; Kroupa, P.; Famaey, B. The Formation of Exponential Disk Galaxies in MOND. Astrophys. J. 2020, 890, 173. [Google Scholar] [CrossRef]
  28. Eappen, R.; Kroupa, P.; Wittenburg, N.; Haslbauer, M.; Famaey, B. The formation of early-type galaxies through monolithic collapse of gas clouds in Milgromian gravity. Mon. Not. R. Astron. Soc. 2022, 516, 1081–1093. [Google Scholar] [CrossRef]
  29. Nagesh, S.T.; Kroupa, P.; Banik, I.; Famaey, B.; Ghafourian, N.; Roshan, M.; Thies, I.; Zhao, H.; Wittenburg, N. Simulations of star-forming main-sequence galaxies in Milgromian gravity. Mon. Not. R. Astron. Soc. 2023, 519, 5128–5148. [Google Scholar] [CrossRef]
  30. Katz, H.; McGaugh, S.; Teuben, P.; Angus, G.W. Galaxy Cluster Bulk Flows and Collision Velocities in QUMOND. Astrophys. J. 2013, 772, 10. [Google Scholar] [CrossRef]
  31. Wittenburg, N.; Kroupa, P.; Banik, I.; Candlish, G.; Samaras, N. Hydrodynamical structure formation in Milgromian cosmology. Mon. Not. R. Astron. Soc. 2023, 523, 453–473. [Google Scholar] [CrossRef]
  32. Angus, G.W.; Diaferio, A.; Kroupa, P. Using dwarf satellite proper motions to determine their origin. Mon. Not. R. Astron. Soc. 2011, 416, 1401–1409. [Google Scholar] [CrossRef]
  33. Bílek, M.; Thies, I.; Kroupa, P.; Famaey, B. Are Disks of Satellites Comprised of Tidal Dwarf Galaxies? Galaxies 2021, 9, 100. [Google Scholar] [CrossRef]
  34. Roshan, M.; Ghafourian, N.; Kashfi, T.; Banik, I.; Haslbauer, M.; Cuomo, V.; Famaey, B.; Kroupa, P. Fast galaxy bars continue to challenge standard cosmology. Mon. Not. R. Astron. Soc. 2021, 508, 926–939. [Google Scholar] [CrossRef]
  35. Oehm, W.; Thies, I.; Kroupa, P. Constraints on the dynamical evolution of the galaxy group M81. Mon. Not. R. Astron. Soc. 2017, 467, 273–289. [Google Scholar] [CrossRef]
  36. Besla, G.; Kallivayalil, N.; Hernquist, L.; Robertson, B.; Cox, T.J.; van der Marel, R.P.; Alcock, C. Are the Magellanic Clouds on Their First Passage about the Milky Way? Astrophys. J. 2007, 668, 949–967. [Google Scholar] [CrossRef]
  37. Besla, G.; Kallivayalil, N.; Hernquist, L.; van der Marel, R.P.; Cox, T.J.; Kereš, D. Simulations of the Magellanic Stream in a First Infall Scenario. Astrophys. J. Lett. 2010, 721, L97–L101. [Google Scholar] [CrossRef]
  38. D’Souza, R.; Bell, E.F. The infall of dwarf satellite galaxies are influenced by their host’s massive accretions. Mon. Not. R. Astron. Soc. 2021, 504, 5270–5286. [Google Scholar] [CrossRef]
  39. Hammer, F.; Yang, Y.B.; Flores, H.; Puech, M.; Fouquet, S. The Magellanic Stream System. I. Ram-Pressure Tails and the Relics of the Collision Between the Magellanic Clouds. Astrophys. J. 2015, 813, 110. [Google Scholar] [CrossRef]
  40. Wang, J.; Hammer, F.; Yang, Y.; Ripepi, V.; Cioni, M.R.L.; Puech, M.; Flores, H. Towards a complete understanding of the Magellanic Stream Formation. Mon. Not. R. Astron. Soc. 2019, 486, 5907–5916. [Google Scholar] [CrossRef]
  41. Vasiliev, E. Dear Magellanic Clouds, welcome back! arXiv 2023, arXiv:2306.04837. [Google Scholar] [CrossRef]
  42. Lucchini, S.; D’Onghia, E.; Fox, A.J. The Magellanic Stream at 20 kpc: A New Orbital History for the Magellanic Clouds. Astrophys. J. Lett. 2021, 921, L36. [Google Scholar] [CrossRef]
  43. Navarro, J.F.; Frenk, C.S.; White, S.D.M. The Structure of Cold Dark Matter Halos. Astrophys. J. 1996, 462, 563. [Google Scholar] [CrossRef]
  44. Behroozi, P.S.; Wechsler, R.H.; Conroy, C. The Average Star Formation Histories of Galaxies in Dark Matter Halos from z = 0–8. Astrophys. J. 2013, 770, 57. [Google Scholar] [CrossRef]
  45. Lukić, Z.; Reed, D.; Habib, S.; Heitmann, K. The Structure of Halos: Implications for Group and Cluster Cosmology. Astrophys. J. 2009, 692, 217–228. [Google Scholar] [CrossRef]
  46. Bullock, J.S.; Kolatt, T.S.; Sigad, Y.; Somerville, R.S.; Kravtsov, A.V.; Klypin, A.A.; Primack, J.R.; Dekel, A. Profiles of dark haloes: Evolution, scatter and environment. Mon. Not. R. Astron. Soc. 2001, 321, 559–575. [Google Scholar] [CrossRef]
  47. Dutton, A.A.; Macciò, A.V. Cold dark matter haloes in the Planck era: Evolution of structural parameters for Einasto and NFW profiles. Mon. Not. R. Astron. Soc. 2014, 441, 3359–3374. [Google Scholar] [CrossRef]
  48. Klypin, A.; Yepes, G.; Gottlöber, S.; Prada, F.; Heß, S. MultiDark simulations: The story of dark matter halo concentrations and density profiles. Mon. Not. R. Astron. Soc. 2016, 457, 4340–4359. [Google Scholar] [CrossRef]
  49. Heinze, F.M.; Despali, G.; Klessen, R.S. Not all subhaloes are created equal: Modelling the diversity of subhalo density profiles in TNG50. Mon. Not. R. Astron. Soc. 2024, 527, 11996–12015. [Google Scholar] [CrossRef]
  50. Macciò, A.V.; Dutton, A.A.; van den Bosch, F.C.; Moore, B.; Potter, D.; Stadel, J. Concentration, spin and shape of dark matter haloes: Scatter and the dependence on mass and environment. Mon. Not. R. Astron. Soc. 2007, 378, 55–71. [Google Scholar] [CrossRef]
  51. Chandrasekhar, S. Principles of Stellar Dynamics; The University of Chicago Press: Chicago, IL, USA, 1942. [Google Scholar]
  52. Pietrzyński, G.; Graczyk, D.; Gieren, W.; Thompson, I.B.; Pilecki, B.; Udalski, A.; Soszyński, I.; Kozłowski, S.; Konorski, P.; Suchomska, K.; et al. An eclipsing-binary distance to the Large Magellanic Cloud accurate to two per cent. Nature 2013, 495, 76–79. [Google Scholar] [CrossRef]
  53. Hilditch, R.W.; Howarth, I.D.; Harries, T.J. Forty eclipsing binaries in the Small Magellanic Cloud: Fundamental parameters and Cloud distance. Mon. Not. R. Astron. Soc. 2005, 357, 304–324. [Google Scholar] [CrossRef]
  54. McConnachie, A.W. The observed properties of dwarf galaxies in and around the local group. Astron. J. 2012, 144, 4. [Google Scholar] [CrossRef]
  55. Van der Marel, R.P.; Sahlmann, J. First gaia local group dynamics: Magellanic clouds proper motion and rotation. Astrophys. J. 2016, 832, L23. [Google Scholar] [CrossRef]
  56. Zivick, P.; Kallivayalil, N.; Linden, S.; Fritz, T.; Besla, G.; van der Marel, R.P.; Kozłowski, S.; Geha, M.C.; Kochanek, C.S.; Anderson, J.; et al. New Proper Motions of the Small Magellanic Cloud Using HST and Implications for Milky Way Mass. Proc. Int. Astron. Union 2017, 13, 394–395. [Google Scholar] [CrossRef]
  57. Binney, J.; Tremaine, S. Galactic Dynamics, 2nd ed.; Princeton University Press: Princeton, NJ, USA, 2008. [Google Scholar]
  58. Rubele, S.; Pastorelli, G.; Girardi, L.; Cioni, M.R.L.; Zaggia, S.; Marigo, P.; Bekki, K.; Bressan, A.; Clementini, G.; de Grijs, R.; et al. The VMC survey-XXXI: The spatially resolved star formation history of the main body of the Small Magellanic Cloud. Mon. Not. R. Astron. Soc. 2018, 478, 5017–5036. [Google Scholar] [CrossRef]
  59. Van der Marel, R.P.; Kallivayalil, N. Third-epoch Magellanic Cloud Proper Motions. II. The Large Magellanic Cloud Rotation Field in Three Dimensions. Astrophys. J. 2014, 781, 121. [Google Scholar] [CrossRef]
  60. McGaugh, S.S. A Precise Milky Way Rotation Curve Model for an Accurate Galactocentric Distance. Res. Notes AAS 2018, 2, 156. [Google Scholar] [CrossRef]
  61. Schönrich, R.; Binney, J.; Dehnen, W. Local kinematics and the local standard of rest. Mon. Not. R. Astron. Soc. 2010, 403, 1829–1833. [Google Scholar] [CrossRef]
  62. D’Onghia, E.; Fox, A.J. The Magellanic Stream: Circumnavigating the Galaxy. Annu. Rev. Astron. Astrophys. 2016, 54, 363–400. [Google Scholar] [CrossRef]
  63. Lucchini, S.; D’Onghia, E.; Fox, A.J.; Bustard, C.; Bland-Hawthorn, J.; Zweibel, E. The Magellanic Corona as the key to the formation of the Magellanic Stream. Nature 2020, 585, 203–206. [Google Scholar] [CrossRef]
  64. Goodman, J.; Weare, J. Ensemble samplers with affine invariance. Commun. Appl. Math. Comput. Sci. 2010, 5, 65–80. [Google Scholar] [CrossRef]
  65. Foreman-Mackey, D.; Hogg, D.W.; Lang, D.; Goodman, J. emcee: The MCMC Hammer. Publ. Astron. Soc. Pac. 2013, 125, 306. [Google Scholar] [CrossRef]
  66. Charbonneau, P. Genetic Algorithms in Astronomy and Astrophysics. Astrophys. J. Suppl. 1995, 101, 309. [Google Scholar] [CrossRef]
  67. Theis, C.; Kohle, S. Multi-method-modeling of interacting galaxies. I. A unique scenario for NGC 4449? Astron. Astrophys. 2001, 370, 365–383. [Google Scholar] [CrossRef]
  68. Di Valentino, E.; Melchiorri, A.; Silk, J. Cosmic Discordance: Planck and luminosity distance data exclude LCDM. arXiv 2020, arXiv:2003.04935. [Google Scholar]
  69. Di Valentino, E.; Melchiorri, A.; Silk, J. Investigating Cosmic Discordance. Astrophys. J. Lett. 2021, 908, L9. [Google Scholar] [CrossRef]
  70. Melia, F. A Candid Assessment of Standard Cosmology. Publ. Astron. Soc. Pac. 2022, 134, 121001. [Google Scholar] [CrossRef]
  71. Di Valentino, E.D. Challenges of the Standard Cosmological Model. Universe 2022, 8, 399. [Google Scholar] [CrossRef]
  72. Di Valentino, E. Cosmological tensions: Hints for a new concordance model? In Proceedings of the Sixteenth Marcel Grossmann Meeting, Online, 5–10 July 2021; On Recent Developments in Theoretical and Experimental General Relativity, Astrophysics, and Relativistic Field Theories. pp. 1770–1782. [Google Scholar] [CrossRef]
  73. Lopez-Corredoira, M. History and Problems of the Standard Model in Cosmology. arXiv 2023, arXiv:2307.10606. [Google Scholar] [CrossRef]
  74. Colgáin, E.Ó.; Pourojaghi, S.; Sheikh-Jabbari, M.M.; Sherwin, D. MCMC Marginalisation Bias and ΛCDM tensions. arXiv 2023, arXiv:2307.16349. [Google Scholar] [CrossRef]
Figure 1. Orbits, displayed as pairwise distances, calculated backwards in time up to 7 Gyr for the best-fit solutions for the mass combinations o-o-o and p-m-p with either statistical method, MCMC or GA. Left panel: Orbits obtained using the MCMC method. Right panel: Orbits obtained using the GA method.
Figure 1. Orbits, displayed as pairwise distances, calculated backwards in time up to 7 Gyr for the best-fit solutions for the mass combinations o-o-o and p-m-p with either statistical method, MCMC or GA. Left panel: Orbits obtained using the MCMC method. Right panel: Orbits obtained using the GA method.
Universe 10 00143 g001aUniverse 10 00143 g001b
Figure 2. The MCMC method: Autocorrelation functions related to the open parameters for the primordial mass combination o-o-o. The calculation is based on a set of 50,000 ensembles.
Figure 2. The MCMC method: Autocorrelation functions related to the open parameters for the primordial mass combination o-o-o. The calculation is based on a set of 50,000 ensembles.
Universe 10 00143 g002
Figure 3. The MCMC method: Same as for Figure 2 but for the overall best-fit mass combination p-m-p.
Figure 3. The MCMC method: Same as for Figure 2 but for the overall best-fit mass combination p-m-p.
Universe 10 00143 g003
Figure 4. Forces between the Magellanic Clouds based on today’s observational data with a distance of 22.5 kpc between the CMs of the dark matter halos and the present-day relative velocity of 89.8 pc/Myr ( 87.8 km/s), shown as the vertical dotted line.
Figure 4. Forces between the Magellanic Clouds based on today’s observational data with a distance of 22.5 kpc between the CMs of the dark matter halos and the present-day relative velocity of 89.8 pc/Myr ( 87.8 km/s), shown as the vertical dotted line.
Universe 10 00143 g004
Figure 5. Centre-of-mass distances, calculated backwards in time up to 7 Gyr for the best-fit solutions for the mass combinations o-o-o and p-m-p. Note that the best-fit solutions delivered by the MCMC and GA methods are indistinguishable (see Figure 1). Left panel: Full calculation including dynamical friction. Right panel: Based on identical initial conditions at present, orbits obtained by switching off dynamical friction.
Figure 5. Centre-of-mass distances, calculated backwards in time up to 7 Gyr for the best-fit solutions for the mass combinations o-o-o and p-m-p. Note that the best-fit solutions delivered by the MCMC and GA methods are indistinguishable (see Figure 1). Left panel: Full calculation including dynamical friction. Right panel: Based on identical initial conditions at present, orbits obtained by switching off dynamical friction.
Universe 10 00143 g005aUniverse 10 00143 g005b
Table 2. List of references for basic observational data.
Table 2. List of references for basic observational data.
RA and DEC for LMC, SMC:NASA Extragalactic Database
Distance LMC:[52]
Distance SMC:[53]
Radial velocities LMC, SMC:[54]
Transverse velocities LMC:[55]
Transverse velocities SMC:[56]
Stellar mass MW:[57]
Stellar mass LMC:[58]
Stellar mass SMC:[59]
Distance Galactic centre:[60]
Solar circular speed:[60]
Solar proper motion:[61]
Table 3. Observational data for LMC and SMC and in parts for the Galactic centre.
Table 3. Observational data for LMC and SMC and in parts for the Galactic centre.
ObjectRADECHeliocentricHeliocentric
(EquJ2000)(EquJ2000)DistanceRadial Velocity
LMC 80.894 69.756 49.97 kpc 262.2 km/s
SMC 13.187 72.829 60.6 kpc 145.6 km/s
MW 266.405 28.936 8.122 kpc
Table 4. Transverse velocity components for LMC and SMC.
Table 4. Transverse velocity components for LMC and SMC.
Object v RA v RA v DEC v DEC
[mas/yr][km/s][mas/yr][km/s]
LMC 1.872 ± 0.045 443.3 ± 10.7 0.224 ± 0.054 53.0 ± 12.8
SMC 0.820 ± 0.060 235.5 ± 17.2 1.230 ± 0.070 353.3 ± 20.1
Table 5. Circular velocity and proper motion of the Sun (Galactic coordinates).
Table 5. Circular velocity and proper motion of the Sun (Galactic coordinates).
v c v u v v v w
233.34 km/s 11.10 km/s 12.24 km/s 7.25 km/s
Table 6. Cartesian equatorial coordinates and heliocentric equatorial velocities for the Magellanic Clouds and the centre of the Galaxy.
Table 6. Cartesian equatorial coordinates and heliocentric equatorial velocities for the Magellanic Clouds and the centre of the Galaxy.
Objectxyz v x v y v z
[kpc][kpc][kpc][km/s][km/s][km/s]
LMC 2.736 17.073 46.883 415.6 208.4 227.8
SMC 17.419 4.081 57.899 340.5 162.1 243.4
MW 0.446 7.094 3.930 114.4 120.4 181.4
Table 7. Transformation of the error bars in RA and DEC for the transverse velocity components of the Magellanic Clouds to Cartesian equatorial coordinates.
Table 7. Transformation of the error bars in RA and DEC for the transverse velocity components of the Magellanic Clouds to Cartesian equatorial coordinates.
Velocity Component Δ v x Δ v y Δ v z
[km/s][km/s][km/s]
LMC RA ± 10.5 ± 1.69 ± 0
LMC DEC ± 1.90 ± 11.9 ± 4.43
SMC RA ± 3.93 ± 16.8 ± 0
SMC DEC ± 18.7 ± 4.39 ± 5.95
Table 8. Open parameters for the statistical methods and the corresponding 1 σ uncertainty values from the observational data.
Table 8. Open parameters for the statistical methods and the corresponding 1 σ uncertainty values from the observational data.
P 1 P 2 P 3 P 4
(LMC: v RA )(LMC: v DEC )(SMC: v RA )(SMC: v DEC )
v 1 v 2 v 3 v 4
443.3 km/s 53.0 km/s 235.5 km/s 353.3 km/s
σ 1 σ 2 σ 3 σ 4
10.7 km/s 12.8 km/s 17.2 km/s 20.1 km/s
Table 9. Results using the MCMC method: Probabilities regarding the evaluation of the error intervals of the plane-of-sky velocity components of the LMC and SMC for the individual best-fit solutions of each mass combination considered, as explained in Section 4.3.2.
Table 9. Results using the MCMC method: Probabilities regarding the evaluation of the error intervals of the plane-of-sky velocity components of the LMC and SMC for the individual best-fit solutions of each mass combination considered, as explained in Section 4.3.2.
CombinationDeviationDeviationBDeviationProbability
of Massesof P 1 of P 2 of P 3 of P 4
o-o-o> 6.8 σ > 3.9 σ < 1 σ > 3.6 σ 3.2 × 10 19
o-o-m> 6.3 σ > 3.1 σ < 1 σ > 1.7 σ 5.1 × 10 14
o-o-p> 6.8 σ < 1 σ < 1 σ > 6.8 σ 1.1 × 10 22
o-m-o> 5.0 σ < 1 σ < 1 σ > 5.7 σ 6.9 × 10 15
o-m-m> 3.4 σ > 1.0 σ < 1 σ > 6.8 σ 7.1 × 10 15
o-m-p> 4.3 σ > 1.2 σ > 1.1 σ > 5.9 σ 3.9 × 10 15
o-p-o> 6.5 σ > 2.4 σ > 1.4 σ > 4.3 σ 3.6 × 10 18
o-p-m> 5.9 σ > 2.3 σ < 1 σ > 2.1 σ 2.8 × 10 12
o-p-p> 6.7 σ > 1.8 σ > 4.3 σ > 6.9 σ 1.3 × 10 28
m-o-o> 5.7 σ > 1.2 σ > 3.2 σ > 6.9 σ 2.0 × 10 23
m-o-m> 6.8 σ > 1.0 σ < 1 σ > 6.9 σ 5.5 × 10 23
m-o-p> 6.2 σ > 2.4 σ < 1 σ > 6.9 σ 4.8 × 10 23
m-m-o> 5.1 σ < 1 σ < 1 σ > 6.0 σ 6.7 × 10 16
m-m-m> 5.1 σ > 1.0 σ < 1 σ > 6.3 σ 1.0 × 10 16
m-m-p> 5.4 σ < 1 σ > 1.0 σ > 5.8 σ 4.4 × 10 16
m-p-o> 6.9 σ > 1.5 σ > 2.8 σ > 6.9 σ 1.9 × 10 26
m-p-m> 6.8 σ < 1 σ > 4.8 σ > 6.8 σ 1.7 × 10 28
m-p-p> 6.9 σ > 2.8 σ > 1.1 σ > 6.8 σ 7.6 × 10 26
p-o-o> 6.2 σ < 1 σ < 1 σ > 7.9 σ 1.6 × 10 24
p-o-m> 6.5 σ > 1.9 σ > 2.1 σ > 7.9 σ 4.6 × 10 28
p-o-p> 5.9 σ > 1.9 σ < 1 σ > 7.7 σ 2.9 × 10 24
p-m-o> 4.9 σ > 1.4 σ < 1 σ > 5.8 σ 1.0 × 10 15
p-m-m> 5.4 σ < 1 σ < 1 σ > 5.9 σ 2.4 × 10 16
p-m-p> 4.7 σ < 1 σ > 2.0 σ > 3.5 σ 5.5 × 10 11
p-p-o> 7.7 σ > 1.2 σ < 1 σ > 7.9 σ 8.8 × 10 30
p-p-m> 7.3 σ < 1 σ > 1.4 σ > 8.4 σ 2.1 × 10 30
p-p-p> 6.9 σ > 2.0 σ > 1.0 σ > 7.8 σ 1.5 × 10 27
Table 10. Same as for Table 9, but for the GA method.
Table 10. Same as for Table 9, but for the GA method.
CombinationDeviationDeviationDeviationDeviationProbability
of Massesof P 1 of P 2 of P 3 of P 4
o-o-o> 6.8 σ > 3.0 σ > 1.1 σ > 4.0 σ 4.9 × 10 19
o-o-m> 6.4 σ < 1 σ < 1 σ > 3.1 σ 3.0 × 10 13
o-o-p> 6.5 σ > 1.2 σ < 1 σ > 6.9 σ 9.7 × 10 23
o-m-o> 4.5 σ > 1.6 σ < 1 σ > 5.9 σ 2.7 × 10 15
o-m-m> 4.7 σ > 1.0 σ < 1 σ > 6.3 σ 7.8 × 10 16
o-m-p> 4.4 σ < 1 σ > 1.9 σ > 5.9 σ 2.3 × 10 15
o-p-o> 6.4 σ > 3.5 σ < 1 σ > 3.5 σ 3.4 × 10 17
o-p-m> 5.9 σ > 3.0 σ < 1 σ > 1.9 σ 5.6 × 10 13
o-p-p> 5.9 σ < 1 σ > 6.0 σ > 6.9 σ 3.8 × 10 29
m-o-o> 6.3 σ > 1.5 σ > 1.4 σ > 6.9 σ 3.4 × 10 23
m-o-m> 6.9 σ > 1.5 σ > 1.6 σ > 6.8 σ 8.0 × 10 25
m-o-p> 6.8 σ > 1.0 σ < 1 σ > 6.7 σ 2.2 × 10 22
m-m-o> 4.6 σ > 1.7 σ > 1.0 σ > 6.1 σ 4.0 × 10 16
m-m-m> 5.2 σ < 1 σ > 2.1 σ > 5.8 σ 4.7 × 10 17
m-m-p> 5.3 σ > 1.0 σ > 1.4 σ > 5.6 σ 4.0 × 10 16
m-p-o> 6.9 σ > 1.4 σ > 3.5 σ > 6.8 σ 4.1 × 10 27
m-p-m> 6.9 σ > 1.3 σ > 4.0 σ > 6.9 σ 3.3 × 10 28
m-p-p> 6.8 σ < 1 σ > 3.0 σ > 6.9 σ 1.5 × 10 25
p-o-o> 6.9 σ > 1.9 σ > 6.0 σ > 6.9 σ 3.1 × 10 33
p-o-m> 7.3 σ < 1 σ > 1.1 σ > 7.8 σ 4.9 × 10 28
p-o-p> 6.3 σ > 1.0 σ > 4.9 σ > 6.9 σ 1.5 × 10 27
p-m-o> 4.9 σ < 1 σ > 1.2 σ > 5.8 σ 1.5 × 10 15
p-m-m> 3.4 σ > 2.0 σ < 1 σ > 6.9 σ 1.6 × 10 16
p-m-p> 5.5 σ > 1.0 σ > 1.2 σ > 2.7 σ 6.1 × 10 11
p-p-o> 7.5 σ < 1 σ > 2.9 σ > 7.9 σ 6.7 × 10 31
p-p-m> 7.8 σ > 4.1 σ > 4.6 σ > 7.7 σ 1.5 × 10 38
p-p-p> 7.4 σ > 1.9 σ < 1 σ > 7.5 σ 5.0 × 10 28
Table 11. Results from the MCMC method: Probabilities regarding the evaluation of the error intervals of the plane-of-sky velocity components of the LMC and SMC for the individual best-fit solutions of the mass combinations o-o-o and p-m-p based on sets of 1000 solutions, as explained in Section 4.3.3.
Table 11. Results from the MCMC method: Probabilities regarding the evaluation of the error intervals of the plane-of-sky velocity components of the LMC and SMC for the individual best-fit solutions of the mass combinations o-o-o and p-m-p based on sets of 1000 solutions, as explained in Section 4.3.3.
CombinationDeviationDeviationDeviationDeviationProbability
of Masses of P 1 of P 2 of P 3 of P 4
o-o-o > 6.8 σ > 3.9 σ < 1 σ > 3.6 σ 3.2 × 10 19
p-m-p > 5.3 σ > 2.8 σ > 1.0 σ > 1.0 σ 5.9 × 10 10
Table 12. Same as for Table 11, but for the GA method.
Table 12. Same as for Table 11, but for the GA method.
CombinationDeviationDeviationDeviationDeviationProbability
of Masses of P 1 of P 2 of P 3 of P 4
o-o-o > 6.9 σ > 3.2 σ < 1 σ > 3.6 σ 2.3 × 10 18
p-m-p > 5.5 σ > 2.2 σ < 1 σ > 1.6 σ 1.2 × 10 10
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Oehm, W.; Kroupa, P. The Relevance of Dynamical Friction for the MW/LMC/SMC Triple System. Universe 2024, 10, 143. https://doi.org/10.3390/universe10030143

AMA Style

Oehm W, Kroupa P. The Relevance of Dynamical Friction for the MW/LMC/SMC Triple System. Universe. 2024; 10(3):143. https://doi.org/10.3390/universe10030143

Chicago/Turabian Style

Oehm, Wolfgang, and Pavel Kroupa. 2024. "The Relevance of Dynamical Friction for the MW/LMC/SMC Triple System" Universe 10, no. 3: 143. https://doi.org/10.3390/universe10030143

APA Style

Oehm, W., & Kroupa, P. (2024). The Relevance of Dynamical Friction for the MW/LMC/SMC Triple System. Universe, 10(3), 143. https://doi.org/10.3390/universe10030143

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