Issue |
A&A
Volume 689, September 2024
|
|
---|---|---|
Article Number | A150 | |
Number of page(s) | 11 | |
Section | Astrophysical processes | |
DOI | https://doi.org/10.1051/0004-6361/202450548 | |
Published online | 10 September 2024 |
Dynamical friction in the quasi-linear formulation of modified Newtonian dynamics (QuMOND)
1
CNR-ISC, Via Madonna del Piano 17, 50022 Sesto Fiorentino, Italy
2
INAF-Osservatorio Astronomico di Arcetri, Largo Enrico Fermi 5, 50125 Firenze, Italy
3
Dipartimento di Fisica e Astronomia, Università di Firenze & INFN-Sezione di Firenze, Via Sansone 1, 50022 Sesto Fiorentino, Italy
4
Dipartimento di Fisica ”Giuseppe Occhialini”, Università di Milano Bicocca, Piazza della Scienza 3, 20126 Milano, Italy
e-mail: federico.re@unimib.it
5
INFN-Sezione di Milano, Via Celoria 15, 20133 Milano, Italy
6
Dipartimento di Scienze Fisiche, Informatiche e Matematiche, Università di Modena e Reggio Emilia, Via Campi 213/A, 41125 Modena, Italy
e-mail: caterina.chiari@unimore.it
7
CNR-NANO, Via Campi 213/A, 41125 Modena, Italy
Received:
29
April
2024
Accepted:
19
June
2024
Aims. We explore the dynamical friction on a test mass in gravitational systems in the quasi-linear formulation of modified Newtonian dynamics (QuMOND).
Methods. Exploiting the quasi-linearity of QuMOND, we derived a simple expression for the dynamical friction in akin to its Newtonian counterpart in the standard Chandrasekhar derivation. Moreover, adopting a mean field approach based on the Liouville equation, we were able to obtain a more rigorous (albeit in integral form) dynamical friction formula that can be evaluated numerically for a given choice of the QuMOND interpolation function.
Results. We find that our results are consistent with those of previous works. We observe that the dynamical friction is stronger in MOND with respect to a baryon-only Newtonian system with the same mass distribution. This amounts to a correction of the Coulomb logarithmic factor via additional terms that are proportional to the MOND radius of the system. Moreover, with the aid of simple numerical experiments, we confirm our theoretical predictions and those of previous works based on MOND.
Key words: gravitation / methods: analytical / stars: kinematics and dynamics / galaxies: kinematics and dynamics
© The Authors 2024
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.
1. Introduction
Modified Newtonian dynamics (hereafter MOND, Milgrom 1983) is an alternative theory of (classical) gravity that was originally introduced to solve the missing mass problem on astrophysical and cosmological scales without resorting to dark matter theory (hereafter DM). In a Lagrangian formulation (see Bekenstein & Milgrom 1984), MOND can be amended to a substitution of the Poisson equation relating the gravitational potential, Φ, and mass density, ρ, with the non-linear field equation:
where the acceleration, scaling as a0 ≈ 1.2 × 10−2ms−2, is a new universal constant; whereas the interpolating function, μ, (which is, in principle, unknown) is monotonic with the asymptotic behaviour, as follows:
A widely adopted form used for μ is:
According to the equations above, for ∇Φ ≫ a0, we can recover the Newtonian limit and Eq. (1) reduces to the Poisson equation. Vice versa, for ∇Φ ≪ a0 the system is in the so-called deep-MOND limit (hereafter, dMOND) and Eq. (1) becomes the p-Laplace equation:
For a given mass density, ρ, the right-hand-side of the Poisson equation and Eq. (1) are the same. We can therefore eliminate ρ, thereby obtaining the relation:
between the MOND and Newtonian force fields gM and gN, where S ≡ ∇ × h(ρ) is a density-dependent solenoidal field that zeros-out for systems in spherical, cylindrical, or planar symmetry, while it is typically non-zero for more general density profiles.
MOND has been rather successful in reproducing the kinematics of galaxies without the need for DM to be included (Bugg 2015; Zhu et al. 2023). In fact, it has proven to be a valid alternative to the ΛCDM paradigm where the latter is challenged, across a broad range of gravitational problems, such as the structural properties of dwarf galaxies (Milgrom 1994; Sánchez-Salcedo & Hernandez 2007; McGaugh & Milgrom 2013; Asencio et al. 2022); radial acceleration relation of galaxies (RAR, Lelli et al. 2017; Mistele et al. 2024b,a); profiles of low-surface-brightness galaxies (Sanders 2021); dynamics of wide binary stars (Hernandez et al. 2012, 2023; Chae 2024); evolution of globular clusters of ultra-diffuse galaxies (Bílek et al. 2021); and the orbital velocity of interacting galaxy pairs (Scarpa et al. 2022).
MOND has been extensively investigated numerically (Brada & Milgrom 1999; Knebe & Gibson 2004; Tiret & Combes 2007; Nipoti et al. 2007a,b,c; Sanders 2008; Wu et al. 2009; Malekjani et al. 2009; Roshan et al. 2021; Kroupa et al. 2022; Banik et al. 2022; Nagesh et al. 2023; Re & Di Cintio 2023; Sollima et al. 2023), for collisionless mean field processes; for example: violent relaxation and phase-mixing in monolithic collapse, galaxy merging, and the vertical dynamics of disk galaxies. Unfortunately, due to the non-linearity of the theory (cf. Eq. 1) and, hence, the absence of the superposition principle, there is much less information available on the collisional processes (see e.g. Bílek et al. 2021). For example, Ciotti & Binney (2004) used an approach based on the fluctuations of uniform field to estimate the MOND two-body relaxation (t2b, Chandrasekhar 1941a,b) and dynamical friction (tDF, Chandrasekhar 1943a,b,c) time scales in the dMOND limit. The outcome of this study (see also the numerical simulations in Nipoti et al. 2008) shows that a test mass, M, crossing a purely baryonic system would feel a stronger dynamical friction (hereafter, DF) force due to the encounters with the system’s stars in MOND than in the parent
Newtonian system without DM; however, it would undergo an only slightly more efficient DF than in the equivalent Newtonian system (ENS). The is comprised of the baryonic plus DM system constructed such that the potential is the same as the purely baryonic MOND model. The Newtonian and MOND DF time scales and
, are related as:
where ℛ is the ratio of the amount of DM to baryons in the ENS.
In MOND, the collisional direct N-body simulations are intrinsically unfeasible; thus, results such as those given by Eq. (6) cannot be exhaustively explored via numerical experiments. The reason for this is the absence of a MONDian expression for the force exerted by a point-like particle. To address this issue, Milgrom (1986) proposed the approximated expression for the force exchanged between two masses m1 and m2 placed at distance, r,
The latter is however valid only in the far field limit, thus making it unusable in a direct N-body code, as it lacks a regime bridging to the Newtonian 1/r2 limit.
More recently, Milgrom (2010) formulated a quasi-linear MOND theory (hereafter QuMOND), where the governing field equation is:
formally identical to Eq. (1) where now, ΦN is the Newtonian potential generated by the density ρ through the usual Poisson equation:
and the new interpolating function ν(y) is related to μ(x) via:
The asymptotic behaviours for ν in the Newtonian and in the dMOND limits are
We note that Eq. (8) can be rearranged as
that is, the QuMOND potential satisfies the Poisson equation for the auxiliary density . In practice, in QuMOND, we have to evaluate the Newtonian potential ΦN generated by ρ, evaluate its gradient and then via a non-linear algebraic step obtain the auxiliary MOND density
via Eq. (12) that acts as source for the potential ΦM.
At variance with the usual MOND formulation, in its quasi-linear formulation the interpolating function ν has a stronger effect on the form of the gravitational potential than its counterpart μ in Eq. (1). For the usual choice (3), we obtain:
We note that QuMOND contains a sort of superposition principle for the auxiliary densities . It is therefore tempting to try to derive an expression for the DF force in this context, making use of purely kinetic arguments as in the original derivation. In this paper, we first discuss a simple expression of the QuMOND DF summing up the contributions of the effective densities for point-particles of mass, m, interacting with a test mass, with some simple choices for ν. Moreover, we extend the mean field formulation of DF, originally derived in Newtonian gravity by Kandrup (1983) to the case of QuMOND and discuss the results of simple numerical experiments.
2. Dynamical friction: Classical approach
2.1. Newtonian case
Before tackling the QuMOND DF formula, we start by briefly recalling the ingredients of the derivation of its Newtonian expression (for a more detailed review, see e.g. Binney & Tremaine 2008 or Ciotti 2021). We consider a test mass, M, travelling at vM through an infinite background medium of particles of a mass, m, and (by now constant) number density, n, with velocity distribution function, f(v). In each encounter, M and a background particle exchange a force expressed as F ≈ GMm/b2, acting for an interval of time, Δt ≈ b/vM, where b is the impact parameter of the dynamical collision. The associated variation of momentum is Δp ≈ FΔt ≈ GMm/vMb. The rate of collisions of M is then
From this, the velocity diffusion coefficient along the trajectory of M can be written as:
In the equation above, the integral in the impact parameter b diverges for both b → 0 and b → ∞ (infrared and ultraviolet divergence, respectively). It is however possible to impose suitable cut-offs bmin and bmax. This factor, after integration takes the form of the Coulomb logarithm (in analogy with the same quantity in plasma physics) lnΛ, where Λ = bmax/bmin.
Let us now assume that f(v), defined such that
is the total number of particles, is position independent and isotropic for all velocities. For such a system, the DF force per unit mass acting on M is given as
where
is the so-called fractional velocity volume function (Ciotti 2021). Therefore, in the assumption of velocity isotropy, only particles moving slower than M contribute to its slowing down. Moreover for M ≫ m, in Eq. (17) the factor m(M + m)n becomes Mρ, where ρ = nm is the background mass density. In Newtonian systems with DM, two species of particles exert DF on the test mass M. Following Ciotti (2010), Eq. (17) can be rewritten for a model with a mass spectrum as
where is a velocity averaged Coulomb logarithm, ⟨m⟩ is the mean mass of the background particles, and
is the total fractional velocity volume function1 of the system. For the case of interest of a binary mass spectrum (i.e. stars and DM particles), we have
where ℳ = M/⟨m⟩ and
In the expressions above, x = nDM/n* and y = mDM/m* are the ratios of the number densities and masses of stars and DM particles, respectively and the ΨDM and Ψ* their fractional velocity volume functions. Stars and DM are not at energy equipartition, as the thermalization time scale would exceed the age of the universe of several order of magnitude due to the small mass of DM particles. Therefore, both species have the same velocity distribution dictated by their collective potential, hence ΨDM = Ψ* = ΨDM. As the DM is constituted by unknown elementary particles, with typical masses of the order of a few GeV, we have mDM ≪ m* and nDM ≫ n*. The average mass ⟨m⟩ is expressed as:
while
so that the total fractional velocity volume functions becomes:
For systems where the test mass M ≫ m* the latter further simplifies and and the DF formula can be simplified to:
In practice, the presence of a DM distribution has enhances the DF of a factor, ρTot/ρ*, with respect to the baryonic system with the same phase-space distribution. We note that Sánchez-Salcedo et al. (2006) also derived a simplified two-component DF formula for the cases of circular orbits in DM plus baryons King (1962) cores, concluding that DF is reduced by the presence of DM of a factor (1 + ℛ)−1 with respect to the case where the (same) mass is in one single component, following the light distribution.
2.2. QuMOND case
We first present a naive general expression of the MONDian DF, making use of the fact that in QuMOND, we have:
For each particle of the background system, we defined the so-called MOND radius as
With the same arguments of Sect. 2.1, we now find that the momentum variation of M is:
The form of the collision rate factor remains unchanged, as it is a geometric quantity. The diffusion coefficient Dv retains the same definition as in Eq. (15), but with a different integral over the impact parameter:
After substituting the expression in the DF formula, we obtain:
From Eq. (30) above, we notice immediately that the QuMOND DF is always larger than its Newtonian equivalent for the same baryonic mass distribution. This happens because ν2 > 1, thus returning a larger factor within the integral over the impact parameters. In addition, this is also in agreement with the fact that in MOND theories the far field behaviour of a given mass distribution falls off less rapidly than in Newtonian gravity, as implied by the logarithmic trend of the MOND potential.
We now evaluate Eq. (30) for a simple choice of ν. Knowing that ν has the asymptotic trends given in Eq. (11) we start with a simplified form:
which is essentially the sum of the Newtonian and the deep MOND régimes. In terms of μ, Eq. (31) corresponds to . For such a choice, the correction factor is
, and the integral on the impact parameters becomes:
where we assumed bmin ≪ rM in the last equality. We note that the maximum impact parameter bmax in the context of Newtonian gravity is usually assumed to be some (arbitrary) scale distance2 of the system at hand. In the original derivation of Chandrasekhar (1943a), where the system is infinite, the definition of bmax remains somewhat unclear. Some authors (see e.g. Van Albada & Szomoru 2020 and references therein) typically take bmax of the order of the average inter-particle distance. Here, we follow the same assumption, noting that for a star of a solar mass, the MOND radius rM is roughly 103 AU; this is much smaller than the typical scale radius of a galaxy, but comparable to the mean inter-particle distance among stars. This established, the QuMOND DF then becomes
In practice, at first order in QuMOND, each particle behaves as a point-like source exerting the usual Newtonian 1/r potential plus an infinitely extended ‘phantom DM’ halo whose contribution depends on the specific form of ν. Following this approach, the calculations as shown in Appendix A leads to the same expression of Eq. (33). We note that Eq. (33) differs from what we would obtain using the simple derivation of t2b in Sect. 2 of Ciotti & Binney 2004; namely, this derivation is proportional to the crossing time to evaluate the dynamical friction time scale with the Spitzer (1987) relation tDF = 2t2bm/(m + M). This is because their expression of t2b makes use of the force (7), implicitly assuming the dMOND regime.
Unlike the usual MOND interpolating function μ(x), the form of the QuMOND function ν(y) affects the behaviour of the MONDian force field. Adopting the commonly used Eq. (13) (see Kroupa et al. 2022), the DF expression becomes:
The explicit derivation is given in Appendix B.
Noteworthy, we observe in both cases that the MOND correction with respect to the Newtonian baryon only case is proportional to a function of the only bmax/rM
The same term always dominates the MOND correction of DF, regardless of the exact form of ν, while the smaller terms are affected by the particular choice. Recently, a number of authors (see e.g. Lelli et al. 2017) proposed an expression for ν supported by observational data proportional to the form of radial acceleration relation. However, its functional form is rather complex thus making the explicit evaluation of Eq. (30) particularly cumbersome.
As an example, in the top panel of Figure 1 we show for a test mass M = 103 M⊙ the ratio ηQuMOND/ηN of the QuMOND and Newtonian dynamical friction coefficient as a function of the stellar density, ρ*, for different values of the velocity dispersion σ (see e.g. Binney & Tremaine 2008) indicated in the figure with the colour map. In the corrective QuMOND term, of Eq. (34) and in the standard Coulomb logarithm, we assumed bmax = (ρ/M⊙)−1/3 and bmin = G(M + m)/σ2. It is evident that for an isolated baryon-only system (corresponding to a dark to luminous matter ratio, ℛ = 0) the DF in QuMOND is augmented (up to a factor ≈102) at low density, comparable to the typical intra-galactic density of 10−1 stars per pc3, with respect to the Newtonian case. Vice-versa (and as expected) in denser systems/regions the QuMOND correction is negligible, in particular for large values of the velocity dispersion, σ. Of course, a fair comparison between Newtonian and MONDian DF has to be made considering the presence of DM. In the bottom panel of Figure 1 we show ηQuMOND/ηN for the same values of ρ* and σ and ℛ = 10 (a typical DM to baryonic matter ratio). In this case, the Newtonian DF is given by Eq. (25). Remarkably, in this case the QuMOND DF is substantially lower than the Newtonian one including the dark mass contribution (substantially confirming the results of Sánchez-Salcedo et al. 2006 for a simplified circular orbit model and of Nipoti et al. 2008; Roshan et al. 2021 for the dynamics of bars in spiral galaxies), for almost all astrophysically relevant combinations of stellar density and velocity dispersion. In the low ρ* regime (say below ≈8 M⊙pc−3), the simplified QuMOND DF expression dominates over its Newtonian counterpart. This leads us to speculate that DM-dominated low ρ* systems, such as ultra-faint dwarf galaxies, could be used as a testing ground for MOND theories versus ΛCDM using the different behaviour of the DF (see e.g. Bílek et al. 2021; Oehm & Kroupa 2024). We stress the fact that in finite systems, both in Newtonian gravity and MOND, the orbital structure (i.e. the velocity anisotropy profile of the background particle distribution) has a relevant effect on the magnitude of DF. For example, a test mass, M, moving on a highly radial orbit in a radially anisotropic spherical system suffers less DF than in a systems with isotropic velocity distribution. In this case, in MOND, where only stars contribute to DF, the drag is further reduced when M reaches radii where no mass is present. A similar argument has been applied (among others) by Angus & Diaferio (2009) to explain in MOND the survival of the globular cluster systems of the Fornax dwarf spheroidal.
![]() |
Fig. 1. Ratio of the QuMOND and Newtonian DF coefficients as function of the stellar density, for ℛ = 0 (top) and 10 (bottom), for the case of a M = 103 M⊙ test mass. The dependence on the local velocity dispersion σ is colour coded. |
3. Mean field approach
A more rigorous evaluation of the DF expression can be carried out using the mean field formalism developed by Kandrup (1983) (see also Gilbert 1968 and Kandrup 1980) in the context of Newtonian gravity that can be applied to any kind of long-range force, not necessarily obeying the superposition principle.
We consider the usual system of N equal mass m (field) particles, with coordinates ri and momenta, pi = mvi, described by the (time-dependent) phase-space distribution function3, ℱ(r, p; t), a test particle M, with a coordinate, r0 = R, and momentum, p0 = P = MvM, which perturbs the initial phase-space distribution ℱ0. By virtue of the third law of dynamics, such perturbation corresponds to the force decelerating the test mass, M. We assume the field particles to be statistically uncorrelated in their initial state with Maxwellian velocity distribution. Moreover, we also take the limit of infinite, homogeneous distribution of field particles. In his original work, Kandrup (1983) initially assumed an external potential, Φ, confining the system (see Appendix C below). We note that MOND systems with self-consistent gravitational field, gin, embedded in an external gravitational field gext are prone to the so-called external field effect (hereafter EFE). The latter implies that for gin < a0 < gext the system is purely Newtonian, while for gin < gext < a0 the system behaves as a Newtonian model with rescaled gravitational constant, G′=Ga0/gext.
Under the assumptions given above, as usual, the average number density of particles, n = N/Vol, should be taken constant as N, Vol → ∞. The dynamical friction force can be therefore obtained simply as
.
In the same fashion of the special relativistic extensions of DF (see e.g. Syer 1994 and Chiari & Di Cintio 2023), we shift to the frame of reference of M such that for the new coordinates, , the linear evolution operator for ℱ takes the form
. See Eqs. (C.7) and (C.9) in the appendix below for its explicit expression.
The linear formulation of the evolution equation is allowed since MOND (and therefore QuMOND), although non-linear, is still a Lagrangian and local4 theory. At this stage, it is not necessary to consider explicitly the Bekenstein & Milgrom (1984) (or Milgrom 2010) Langrangian, since it is enough to consider the particle’s energies.
The significant difference with respect to the Newtonian case explored in Kandrup (1983) is that in MOND, we cannot substitute the potentials in a form:
due to the non-linearity of MOND implied by the absence of a Superposition principle. Such potentials however can be rather calculated in the QuMOND formalism (8).
A general expression for the DF force on the test particle can be written formally as
where Gℒ is the Greenian of the operator ℒ and dΣ = d3Nrd3Np is the differential element in phase space. The details of the derivation are discussed in Appendix C.
We now sketch the main simplifying hypotheses needed to perform practically the integrals in Eq. (37). Assuming that the system under consideration has a finite memory (as implied by its ergodicity), we can replace the time integration in τ on the finite interval [0; t] with an integration extended over the semi-infinite interval [0; ∞). Moreover, making use of the standard linear trajectory approximation (see e.g. Ter Haar 1977, see also Syer 1994), allows us to simplify Gℒ(τ → t)[Q]≅Q(t − τ). The latter is the most delicate step of the present calculation, since in this approximation, the effects of the different field particles decouple that is, in principle, not valid in non-linear theory.
Under the assumptions listed above, the DF formula becomes:
where , Ω is the phase space volume occupied by the system and only the field particles slower than VM induce the drag, so that we only for account their density, ρ|v < VM. At this stage, the problem is reduced to solving a three dimensional integral for a given choice of the QuMOND interpolation function, ν.
Eq. (38) in the Newtonian limit (i.e. ν → 1) can be solved by an integration, incrementally applying:
and
thus, we have:
Here, the term is the standard impact parameter (see again Binney & Tremaine 2008), thus, the integral
returns the usual Coulomb logarithm (see e.g. Kalnajs 1972), while the integral in d3v is rewritten as in the Newtonian case, therefore yielding Eq. (17).
The QuMOND case is, even for the simple form of ν given in Eq. (31), plagued by its intrinsic complexity. Analogously to the naive formula in Eq. (30) we notice the presence of two factors of ν. In this mean field formalism, the arguments for the two functions ν in (38) are still slightly different. From such difference , we recovers the logarithmic term in τ, as well as any additional QuMOND term. At variance with the Newtonian case, which requires some cut-off value only for the integral over dτ, in the QuMOND DF formula of Eq. (38), the integral over d3s is also diverging; therefore, it requires a cut-off at bmax. This yields terms of the form bmax/rM and
, which are vastly in agreement with the simplified approach described in Sect. 2.2.
In the deep MOND regime, where ν(y)→y−1/2, we can simplify Eq. (38) as:
We recall that Sánchez-Salcedo et al. (2006) (see also the numerical simulations of Bílek et al. 2021), used the circular orbit approximation and the estimates from Ciotti & Binney (2004) to formulate a local dMOND DF expression for a spherical system where the radial friction on M is expressed as:
In the expression above, is the typical Newtonian acceleration of the system with scale radius rs, and the circular velocity profile is given (see Milgrom 1983) as a function of the stellar radial mass profile M*(r), as follows:
Taking the same assumptions applied in Sánchez-Salcedo et al. 2006, we can simplify Eq. (43) as v and s are generally likely to be perpendicular, while v and are almost parallel. This allows us to then simplify the expression of the term
in scalar form and perform the integral in dτ before those in velocity, dv, and in space, ds, having imposed spherical symmetry. Standard (albeit tedious) algebra then yields
where the dMOND mean square velocity is given (Milgrom 1995) as
for a system with total baryonic mass M* and bmax (typically of the order of the local mean inter-particle distance) is fixed to a scale radius rs.
In Fig. 2, we show the DF acceleration in the limit of circular orbits on a mass of M = 2 × 105 M⊙, as function of the radial distance, r, as given in Eqs. (44) and (46), sinking in a cored spherical model (cfr. also Eq. 49 in the next section) with M* = 2 × 107 M⊙ and a core radius of rs = 0.3 kpc, corresponding to the model of a Fornax-like dwarf galaxy considered by Sánchez-Salcedo et al. 2006. With such a choice of parameters, the constant dMOND velocity dispersion given by Eq. (47) is ≈11km/s. We note that at large radii (e.g. above the core radius indicated in Figure by the vertical dashed line), the two expressions are somewhat comparable; however, in the inner part of the system, the DF is considerably lower in our approximated formulation. This is because for r → 0, the l.h.s. of Eq. (46) is mostly dominated by the contribution of Ψ(vc), that for the specific (but questionable) choice of a Maxwellian velocity distribution is expressed as:
![]() |
Fig. 2. Radial profiles of the DF force (per unit mass) in the circular orbit approximation as given by Eqs. (44,46); purple and cyan lines, respectively. |
then falls rapidly to zero, while in Eq. (44), it is partially compensated by the factor .
Interestingly, the QuMOND N-body simulations of Bílek et al. (2021) offer a hint that Eq. (44) overestimates the DF in the inner region of cored systems; whereas at large radii, the semianalytical estimate of Sánchez-Salcedo et al. (2006) is well matched by the simulations. The discrepancy is interpreted as an effect of the so-called dynamical buoyancy or core stalling (Banik & van den Bosch 2021) in a nearly harmonic core, preventing the further sinking of the test particle, M. We note that in MOND (and in particular in the dMOND limit assumed here), an almost constant density core does not exert a harmonic potential, thereby potentially undermining the core-stalling argument in MOND simulations. This leads to conjecture that the intrinsically less effective dMOND DF in the N-body simulation is indeed responsible for the discrepancy with the prediction of Eq. (44).
4. Numerical experiments
To clarify the different behaviour of DF in the QuMOND and Newtonian gravity model of a given stellar system, it is useful to perform a simple numerical experiment in a set-up where Newtonian gravity with DM the friction on the test object, M, would be negligible, while the gravitational field of the parent system would be in the MONDian regime, g ≪ a0, over a broad interval of radii.
We consider a M = 102 M⊙ star cluster moving through a dwarf galaxy with stellar mass of M* = 2 × 105 M⊙ and scale radius of rs = 0.8 kpc (parameters compatible to those of the Draco ultra-faint dwarf galaxy, Mashchenko et al. 2006), with a stellar distribution given by
corresponding to a Dehnen (1993)γ-model for γ = 0. The associated QuMOND gravitational potential given by Eq. (8) is (by construction) identical to the potential exerted by a Newtonian system with stellar density, as in Eq. (49), plus a DM halo with (cfr. Eq. 12). By contrast, the DF force experienced by M is different in the two paradigms: Newtonian and MOND, with Eq. (25) displaying the first case5 and Eq. (34) the second.
Using the numerical approach discussed in Appendix D, (see also Pasquato & Di Cintio (2020), Di Cintio et al. (2020, 2022)), we integrated different orbits under the effect of the same gravitational field and the two different DF expressions, under the assumption of local Maxwellian approximation (i.e. at all radii the velocity distribution is an isotropic Maxwell-Boltzmann with velocity dispersion approximated by ). In Figure 3, we show the evolution of the galactocentric distance, rg, of M over 12 Gyr for an orbit of initial ellipticity, e = 0.49. The green and blue curves mark the Newtonian and QuMOND simulations, respectively, while the black curve is the unperturbed orbit in the static potential of the model. We observe that while the Newtonian DF only minimally alters the orbit of M at around 9 Gyr, in QuMOND, the test particle suffers a rather strong orbital decay, with the pergalactic radius falling from ≈1.63 kpc down to ≈1.2 kpc at 12 Gyr. This is remarkably in agreement with the analytical estimates of Ciotti & Binney (2004), who predicted a MOND inspiral time of less than a Hubble time for an object that is considerably more massive than a star to the central regions of a Draco-like system. This time span is much shorter than its Newtonian analog and argued that this could be the reason behind the absence of star clusters in similar dwarf galaxies.
![]() |
Fig. 3. Evolution of the galactocentric radius, rg, for a star cluster orbiting through a dwarf galaxy in Newtonian gravity (green curve) and QuMOND (blue curve). The thin black line marks the case without dynamical friction. |
5. Conclusions
We investigated the dynamical friction in the Quasi linear formulation of MOND. Using a simple dimensional analysis, we find that the expression for the MONDian DF is augmented with respect to its Newtonian counterpart of a factor proportional to the MOND radius of the star, rM. This additional term becomes relevant when the maximum impact parameter (proportional to the average interparticle distance, Spitzer 1987) is bmax ≫ rM, as seen by Eqs. (33 and 34). In practice, in a dense star cluster or in galactic nuclei with a mean stellar density and where the typical inter-stellar distance is always smaller than rM, the DF is always purely Newtonian. Vice versa, in low-density systems, such as ultra-faint dwarf galaxies (Simon 2019) the bmax/rM correction is always of the order of the Coulomb logarithm, thus enhancing the (stellar) DF for a system (where, in principle, it should be negligible). In addition, we find that the explicit form of the QuMOND correction is dependent on the specific choice of the interpolating function ν. However, the strength of the enhancing term does not vary significantly for such different forms of ν. We also compared our results with a Newtonian plus DM expression obtained from the Ciotti (2010) formalism for DF in presence of a mass spectrum, confirming that in almost all astrophysically relevant regimes, the Newtonian DF (in presence of dark matter) dominates over its MOND counterpart.
Moreover, we have extended the mean field treatment of DF, pioneered by Kandrup (1983) in the Newtonian case to QuMOND. In this framework, alternatively to the fluctuation-based approach by Ciotti & Binney (2004), we recovered an integral expression that contains the usual dependence on ν. Unfortunately, the explicit evaluation of the DF force becomes rapidly cumbersome (at variance with the simpler Newtonian case) even in the dMOND limit; thus, we are forced to integrate Eq. (38) numerically (specific cases of this will be discussed in a future work). Working out in QuMOND the two-body relaxation time (and therefore the DF coefficient of Eq. 6 via the relation tDF = 2t2bm/M) with the fluctuation approach in Fourier space of Ciotti & Binney (2004) may, in principle, be possible; however, it would imply a double application of the classical Poisson equation and a single non-linear algebraic step is hindered by the implicit relation between the Fourier transform of the QuMOND potential and the associated auxiliary density . Using the mean field approach, we worked out an approximate formula, valid in the deep MOND regime for circular orbits analogous to the one presented in Sánchez-Salcedo et al. (2006). For the case of a globular cluster orbiting in a isolated Fornax-like dwarf galaxy, the two expressions are comparable, at least for radii larger than the core radius of the dwarf, below which our expression underestimates the DF. This result appears to be in agreement with the N-body simulation of Bílek et al. (2021).
Interestingly, simple numerical integration of a test particle in a QuMOND and its Newtonian equivalent system evidence that dynamical friction acts in a different way in the two paradigms, being considerably stronger in the first if the model is in a deep MOND regime (such as for DM-dominated dwarf galaxies). However, we stress the fact that in MOND, we cannot simply add a semi-analytic DF force to a mean field potential obtained solving numerically Eq. (1) for a given density (either imposed extrapolated from particles position) in the same fashion as in Alessandrini et al. (2014), Arca-Sedda & Capuzzo-Dolcetta (2014), Di Cintio et al. (2022). Therefore, our simple numerical estimates could potentially overestimate the MONDian DF. Using the quasi-linear formulation of MOND could, in principle, allow us to incorporate the contribution of DF (and possibly density fluctuations on a scale smaller than the particle-resolution) by adding it before evaluating to solve Eq. (8) (see Di Cintio 2023). In conclusion, it is also important to recall that as a consequence of the EFE, in systems with gint < gext < a0, the dynamics is essentially Newtonian with a rescaled G. In those cases, the DF force would be given by Eq. (17) augmented by the multiplicative factor (a0/gext)2. Therefore, a full mean field QuMOND treatment including the DF expression (discussed in Sect. 3) is only valid for an isolated system.
We note that Ciotti (2010) dubs with velocity volume function Ξ = 4πΨ.
Such a cut-off length in neutral plasma physics is unequivocally constrained by the Debye length, see Spitzer (1965).
MOND is local in the sense that the density at point ξ equals a differential operator acting on the potential in ξ. Vice versa, for the case of 1/rα forces with α ≠ 2, obeying the superposition principle this is not true (Stein 1970). In practice, for such general long-range interactions, density and potential are related only via integral relation over the whole domain occupied by the system (Di Cintio & Ciotti 2011; Di Cintio et al. 2013).
We recall that in a multimass system (e.g. stars and elementary particle sized DM), when the test mass, M, is much larger than the mean mass, ⟨m⟩, then the prefactor of Eq. (17) is dominated by M times the total mass density, as described in e.g. Ciotti (2021).
Acknowledgments
We express gratitude to Jan Pflamm-Altenburg, Pavel Kroupa and Luca Ciotti for the useful discussions at an early stage of this work. We also thank the anonymous Referee for his/her useful remarks that helped improving the presentation of our results.
References
- Alessandrini, E., Lanzoni, B., Miocchi, P., Ciotti, L., & Ferraro, F. R. 2014, ApJ, 795, 169 [NASA ADS] [CrossRef] [Google Scholar]
- Angus, G. W., & Diaferio, A. 2009, MNRAS, 396, 887 [CrossRef] [Google Scholar]
- Arca-Sedda, M., & Capuzzo-Dolcetta, R. 2014, ApJ, 785, 51 [Google Scholar]
- Asencio, E., Banik, I., Mieske, S., et al. 2022, MNRAS, 515, 2981 [CrossRef] [Google Scholar]
- Banik, U., & van den Bosch, F. C. 2021, MNRAS, 502, 1441 [NASA ADS] [CrossRef] [Google Scholar]
- Banik, I., Thies, I., Truelove, R., et al. 2022, MNRAS, 513, 129 [NASA ADS] [CrossRef] [Google Scholar]
- Bekenstein, J., & Milgrom, M. 1984, ApJ, 286, 7 [NASA ADS] [CrossRef] [Google Scholar]
- Bílek, M., Zhao, H., Famaey, B., et al. 2021, A&A, 653, A170 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Binney, J., & Tremaine, S. 2008, Galactic Dynamics: Second Edition (Princeton University Press) [Google Scholar]
- Brada, R., & Milgrom, M. 1999, ApJ, 519, 590 [NASA ADS] [CrossRef] [Google Scholar]
- Bugg, D. 2015, Canad. J. Phys., 93, 119 [NASA ADS] [Google Scholar]
- Chae, K.-H. 2024, ApJ, 960, 114 [NASA ADS] [CrossRef] [Google Scholar]
- Chandrasekhar, S. 1941a, ApJ, 93, 285 [NASA ADS] [CrossRef] [Google Scholar]
- Chandrasekhar, S. 1941b, ApJ, 93, 323 [NASA ADS] [CrossRef] [Google Scholar]
- Chandrasekhar, S. 1943a, ApJ, 97, 255 [Google Scholar]
- Chandrasekhar, S. 1943b, ApJ, 97, 263 [NASA ADS] [CrossRef] [Google Scholar]
- Chandrasekhar, S. 1943c, ApJ, 98, 54 [NASA ADS] [CrossRef] [Google Scholar]
- Chiari, C., & Di Cintio, P. 2023, A&A, 677, A140 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Ciotti, L. 2010, Am. Inst. Phys. Conf. Ser., 1242, 117 [NASA ADS] [Google Scholar]
- Ciotti, L. 2021, Introduction to Stellar Dynamics (Cambridge University Press) [CrossRef] [Google Scholar]
- Ciotti, L., & Binney, J. 2004, MNRAS, 351, 285 [NASA ADS] [CrossRef] [Google Scholar]
- Dehnen, W. 1993, MNRAS, 265, 250 [NASA ADS] [CrossRef] [Google Scholar]
- Di Cintio, P. 2023, https://doi.org/10.5281/zenodo.10397120 [Google Scholar]
- Di Cintio, P., & Casetti, L. 2022, in Multi-Scale (Time and Mass) Dynamics of Space Objects, eds. A. Celletti, C. Galeş, C. Beaugé, & A. Lemaître, 364, 152 [NASA ADS] [Google Scholar]
- Di Cintio, P., & Ciotti, L. 2011, Int. J. Bifur. Chaos, 21, 2279 [Google Scholar]
- Di Cintio, P., Ciotti, L., & Nipoti, C. 2013, MNRAS, 431, 3177 [NASA ADS] [CrossRef] [Google Scholar]
- Di Cintio, P., Ciotti, L., & Nipoti, C. 2020, IAU Symp., 351, 93 [NASA ADS] [Google Scholar]
- Gilbert, I. H. 1968, ApJ, 152, 1043 [NASA ADS] [CrossRef] [Google Scholar]
- Hernandez, X., Jiménez, M. A., & Allen, C. 2012, Eur. Phys. J. C, 72, 1884 [NASA ADS] [CrossRef] [Google Scholar]
- Hernandez, X., Verteletskyi, V., Nasser, L., & Aguayo-Ortiz, A. 2023, MNRAS, 528, 4720 [Google Scholar]
- Kalnajs, A. J. 1972, Astrophys. Space Sci. Libr., 31, 13 [NASA ADS] [CrossRef] [Google Scholar]
- Kandrup, H. E. 1980, Phys. Rep., 63, 1 [NASA ADS] [CrossRef] [Google Scholar]
- Kandrup, H. E. 1983, Ap&ss, 97, 435 [NASA ADS] [CrossRef] [Google Scholar]
- Kandrup, H. E., Pogorelov, I. V., & Sideris, I. V. 2000, MNRAS, 311, 719 [NASA ADS] [CrossRef] [Google Scholar]
- King, I. 1962, AJ, 67, 471 [Google Scholar]
- Knebe, A., & Gibson, B. K. 2004, MNRAS, 347, 1055 [NASA ADS] [CrossRef] [Google Scholar]
- Kroupa, P., Jerabkova, T., Thies, I., et al. 2022, MNRAS, 517, 3613 [CrossRef] [Google Scholar]
- Lelli, F., McGaugh, S. S., Schombert, J. M., & Pawlowski, M. S. 2017, ApJ, 836, 152 [Google Scholar]
- Malekjani, M., Rahvar, S., & Haghi, H. 2009, ApJ, 694, 1220 [CrossRef] [Google Scholar]
- Mannella, R. 2004, Phys. Rev. E, 69, 041107 [NASA ADS] [CrossRef] [Google Scholar]
- Mashchenko, S., Sills, A., & Couchman, H. M. 2006, ApJ, 640, 252 [NASA ADS] [CrossRef] [Google Scholar]
- McGaugh, S., & Milgrom, M. 2013, ApJ, 766, 22 [NASA ADS] [CrossRef] [Google Scholar]
- Milgrom, M. 1983, ApJ, 270, 365 [Google Scholar]
- Milgrom, M. 1986, ApJ, 302, 617 [NASA ADS] [CrossRef] [Google Scholar]
- Milgrom, M. 1994, ApJ, 429, 540 [NASA ADS] [CrossRef] [Google Scholar]
- Milgrom, M. 1995, ApJ, 455, 439 [NASA ADS] [CrossRef] [Google Scholar]
- Milgrom, M. 2010, MNRAS, 403, 886 [NASA ADS] [CrossRef] [Google Scholar]
- Mistele, T., McGaugh, S., Lelli, F., Schombert, J., & Li, P. 2024a, ApJ, 969, L3 [NASA ADS] [CrossRef] [Google Scholar]
- Mistele, T., McGaugh, S., Lelli, F., Schombert, J., & Li, P. 2024b, JCAP, 2024, 020 [Google Scholar]
- Nagesh, S. T., Kroupa, P., Banik, I., et al. 2023, MNRAS, 519, 5128 [NASA ADS] [CrossRef] [Google Scholar]
- Nipoti, C., Londrillo, P., & Ciotti, L. 2007a, ApJ, 660, 256 [NASA ADS] [CrossRef] [Google Scholar]
- Nipoti, C., Londrillo, P., & Ciotti, L. 2007b, MNRAS, 381, L104 [Google Scholar]
- Nipoti, C., Londrillo, P., Zhao, H., & Ciotti, L. 2007c, MNRAS, 379, 597 [NASA ADS] [CrossRef] [Google Scholar]
- Nipoti, C., Ciotti, L., Binney, J., & Londrillo, P. 2008, MNRAS, 386, 2194 [NASA ADS] [CrossRef] [Google Scholar]
- Oehm, W., & Kroupa, P. 2024, Universe, 10, 143 [NASA ADS] [CrossRef] [Google Scholar]
- Pasquato, M., & Di Cintio, P. 2020, A&A, 640, A79 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Re, F., & Di Cintio, P. 2023, A&A, 678, A110 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Roshan, M., Banik, I., Ghafourian, N., et al. 2021, MNRAS, 503, 2833 [NASA ADS] [CrossRef] [Google Scholar]
- Sánchez-Salcedo, F. J., Reyes-Iturbide, J., & Hernandez, X. 2006, MNRAS, 370, 1829 [CrossRef] [Google Scholar]
- Sánchez-Salcedo, F. J., & Hernandez, X. 2007, ApJ, 667, 878 [CrossRef] [Google Scholar]
- Sanders, R. H. 2008, MNRAS, 386, 1588 [NASA ADS] [CrossRef] [Google Scholar]
- Sanders, R. H. 2021, MNRAS, 507, 803 [NASA ADS] [CrossRef] [Google Scholar]
- Scarpa, R., Falomo, R., & Treves, A. 2022, MNRAS, 512, 544 [NASA ADS] [CrossRef] [Google Scholar]
- Simon, J. D. 2019, ARA&A, 57, 375 [NASA ADS] [CrossRef] [Google Scholar]
- Sollima, A., Nipoti, C., Calura, F., Pascale, R., & Baumgardt, H. 2023, MNRAS, 524, 5291 [Google Scholar]
- Spitzer, L. 1965, Physics of Fully Ionized Gases (Interscience Publishers) [Google Scholar]
- Spitzer, L. 1987, Dynamical Evolution of Globular Clusters (Princeton University Press) [Google Scholar]
- Stein, E. M. 1970, Singular Integrals and Differentiability Properties of Functions (PMS-30) (Princeton University Press) [Google Scholar]
- Syer, D. 1994, MNRAS, 270, 205 [NASA ADS] [CrossRef] [Google Scholar]
- Ter Haar, D. 1977, Lectures on Selected Topics in Statistical Mechanics [Google Scholar]
- Tiret, O., & Combes, F. 2007, A&A, 464, 517 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Van Albada, T. S., & Szomoru, A. 2020, Proc. Int. Astron. Union, 351, 532 [NASA ADS] [Google Scholar]
- Wu, X., Zhao, H., Wang, Y., Llinares, C., & Knebe, A. 2009, MNRAS, 396, 109 [NASA ADS] [CrossRef] [Google Scholar]
- Zhu, Y., Ma, H.-X., Dong, X.-B., et al. 2023, MNRAS, 519, 4479 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A: Effective dark matter halos
Combining Eqs. (9) and (8), we evaluate the equivalent matter distribution from Eq. (12) as
For a point particle, the density becomes ρ(r) = mδ3(r), so that gN(r) = − Gm/r3r, and
Substituting the simple interpolation (31), we then find
where the first term here is the real point particle, and the second one represents its phantom dark matter halo . Its cumulative mass function is
Following the approach of Chandrasekhar, the impulsive force generated by the particle and its effective DM halo is then
The associated momentum variation is and then the diffusion coefficient becomes
Keeping in mind that bm ≪ rM < bmax, the integral can be evaluated as . The QuMOND DF force is then
Appendix B: Explicit choice of the QuMOND interpolating function
If the usual choice (13) of the QuMOND interpolating function ν is used in Eq. (30), we have
Unfortunately, here we cannot take the usual approximation for small or large impact parameters b, since the integral runs from bmin < rM to bmax > rM. Vary large values of bmax (say of the order of the system size), return an important contribution already in the Newtonian case, in the form of large values of logΛ. This must also hold true a fortiori in MOND, as in the far-field the forces decay as 1/r rather than 1/r2. Therefore, we cannot expand ν to low-order terms in b/rM.
We must then solve the exact integral ∫ν(b)2db/b, where the square root term yields
where x = 2b/rM. Evaluating Eq. (B.2) for x = 2bmin/rM → 0 returns approximately
while for x = 2bmax/rM → ∞ becomes asymptotically
The integral (B.2) is
The integral over the impact parameters is then evaluated as
Appendix C: Further details on the mean field approach
For the system of N + 1 particles considered in Section 3 described by the distribution function ℱ(r, p; t) the average of a given phase-space observable Q, taken with respect to said distribution function is
where dΣ = d3Nrd3Np is the differential phase-space element.
The initial configuration ℱ0 is perturbed by the passage of the (N + 1)-th particle, (our usual test particle). In virtue of the Third Law of Dynamics, this corresponds to the force decelerating the test particle once the force due to the unperturbed configuration ℱ0 is subtracted:
The field particles are assumed statistically uncorrelated in their initial state . We also assume their velocity distribution to be a Maxwellian:
with temperature β and normalization 𝒩. The distribution of the spatial coordinates is chosen such that the field particles are limited within a certain volume; therefore, it can be also expressed as f1(ri) = e−βmiΦi(ri), with an efficient potential Φi.
For our purposes, we set f1 as the limit of an infinite, homogeneous distribution of field particles. Such a homogeneous distribution on a finite volume Vol ⊆ ℝ3 is given for , namely, when Φi has the profile of a rigid box. The average number density of particles can be then defined as n = N/|Vol|. In the limit of an infinite homogeneous system, this parameters should be taken constant as N, |Vol|→∞. The unperturbed force
vanishes, so that the dynamical friction
can be obtained simply as
.
The field particles evolves following the time-dependent Hamiltonian:
As discussed in Sect. (3), we cannot make the usual substitution of Wi = ∑j ≠ iWij(|ri − rj|) with Wij(r) = − Gmimj/r. This is because, in MOND, the superposition principle does not hold. The potentials should be rather expressed in using the QuMOND formalism (8), so that:
As a consequence, it is useless to define the forces between couples of particles F(j → i) = − ∂riWij(|ri − rj|) = −F(i → j), given that the total QuMOND force on the i−th particle is not simply . The force
can be calculated instead, only as
where Wi(r) is given by the non-linear QuMOND formula (C.5). The evolution equation is therefore expressed as
The sums here arise from those in the Hamiltonian (C.4). In practice, the meaning of Eq. (C.7) is just that the position of each field particle evolves according to the velocity at each instant, as well as the momentum follows the total force.
Following Kandrup (1983), it is useful to change the phase-space coordinates by taking the test particle’s reference:
The evolution operator for ℱ hence takes the form
where Ṗ acts as an apparent force on the i-th particle. Except for this, in this frame the forces do not change , since even in MOND they depend only on the respective particles positions. Moreover, it is useful to define the perturbation
of the distribution function so that
for any quantity Q. This will be used for evaluating Eq. (C.2), so that . The distribution function
evolves as
where again , and the term
vanishes in the infinite, homogeneous limit. In principle, keeping the contribution of
would allow us to investigate in a self-consistent fashion the MOND EFE, see Sect. 3. At this point, Eq. (C.11) can be rewritten in its simpler form:
Its solution can thus be formally expressed as
with the Greenian Gℒ of the operator ℒ. After all the substitution we have a formal, general expression for the dynamical friction on the test particle in the form
where becomes
in the infinite homogeneous system limit. Here we have exploited the fact that the full set of N + 1 particles is a closed, classic dynamical system, so that
.
To perform the integrals appearing in Eqs. (C.14), we now make some simplifying hypotheses. First of all, we perform the infinite homogeneous limit. Second, if the system has only a finite memory, we can replace the integration of τ on [0; t] with an integration of [0; ∞). Moreover, we make use of the standard linear trajectory approximation, which allows us to simplify Gℒ(τ → t)[Q]≅Q(t − τ). Under these approximations, the DF force becomes
This can be further simplified since ⟨FiFj⟩ = 0 for the off-diagonal combinations i ≠ j, since the field particles are by definition statistically uncorrelated. We are then left with
We note that the last step is justified only if the N field particles have all the same mass m, so that they are statistically equivalent in the initial phase-space distribution. Averaging ℱ0 for the i = 1 particle, practically means that only the f(r1, p1) distribution function is relevant, and hence
where is again the Maxwellian velocity distribution. The mean square speed is thus given by ⟨v2⟩ = 3/mβ that can be substituted to obtain
since and
, being
.
Appendix D: Numerical methods
The dynamics of a particle confined by a static potential Φ under the combined effect of force fluctuations and friction is given by the Langevin-type equations
In the most general case, the DF coefficient η depends explicitly on the position and velocity through the phase-space distribution of the embedding system (also sourcing the potential Φ) and δF is a fluctuating force per unit mass with distribution and amplitude connected to η via a fluctuation-dissipation relation. For the systems of interest (Newtonian or MOND) we solve Eq. (D.1) using the Mannella (2004) scheme
written here for simplicity for a 1D system and fixed time step Δt, where the coefficients c1, c2 and d1 are given by
In the equations above is sampled from a norm 1 Gaussian and ζ in the case of a delta correlated noise is fixed by the standard deviation of the distribution of F as
In practice, if such distribution is unknown, we would then, following Kandrup et al. (2000), assume ζ = ℰ where ℰ is the instantaneous relative (positive) energy per unit mass of the test particle along the orbit. In the simulations discussed in this paper the contribution of the fluctuating term is artificially set to 0 as it would be negligible for the type of gravitational systems under consideration. We note that if ν = ζ = 0 Eqs. (D.2) simply become the standard Verlet second order scheme in the drift-kick-drift form.
All Figures
![]() |
Fig. 1. Ratio of the QuMOND and Newtonian DF coefficients as function of the stellar density, for ℛ = 0 (top) and 10 (bottom), for the case of a M = 103 M⊙ test mass. The dependence on the local velocity dispersion σ is colour coded. |
In the text |
![]() |
Fig. 2. Radial profiles of the DF force (per unit mass) in the circular orbit approximation as given by Eqs. (44,46); purple and cyan lines, respectively. |
In the text |
![]() |
Fig. 3. Evolution of the galactocentric radius, rg, for a star cluster orbiting through a dwarf galaxy in Newtonian gravity (green curve) and QuMOND (blue curve). The thin black line marks the case without dynamical friction. |
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.