Arxiv:1802.01493V7 [Astro-Ph.GA] 10 Feb 2021 Ceeainbtti Oe Eoe Rbeai Tselrsae.I Scales

Home , Galactic disc

arXiv:1802.01493v7 [astro-ph.GA] 10 Feb 2021 ceeainbtti oe eoe rbeai tselrsae.I scales. stellar acceleration at MOND problematic the becomes require model this but acceleration ysatn ihacnomlyivratgaiyato,poue an scale gravit produces phenomenological the some modifies action, below [14] gravity is Milgrom modifyin invariant by acceleration of (MOND) conformally Dynamics example Newtonian a An dy Modified an with quantities. of [10]-[11] large lack starting experiments at the by particles and EDGES theorem” dark the galaxy for co from dwarf search it ”Dual temperature the hydrogen successful of [7 w low failing If formation up the galaxy DM. in come [8], Catastrophe” A without to Momentum attempt ”Angular requires, the general including discrepancy in mass theories (BEC) gravity where Condensate a Modified Bose-Einstein San DM speeds and a baryons. [4] rotational proposes Milgrom the [6] [3], between Khoury sca [2]-Lelli front, connection observed McGaugh non-Newtonian see small a galaxies, the for rich However, gas particles. of (BTFR) (DM) matter dark cold niil asi eurdo nices fteMN acceleration MOND the galax of increase of an region or X-ray required is temperature mass high invisible of treatment acceleration po1Gc[3.Gvntemn ucse fGnrlRlvt,teq the Relavity, General of successes many the Given [33]. Gpc 1 upto fe eobnto,Ddlo 2]Segl[3.I ag cl stru scale smaller large significantly In a [23]. cosm with Spergel MOND non-bary in [22] strong overdensity Dodelson a o the recombination, requires T use after further of the spectrum gravity. acceleration limits power not however of matter peculiar worth epoch The speed a pre-recombination is acceptable the of It in an stands, that acceleration provide it [16],[20]. to to As continues MOND shown origin of is relativistic scales. [21] a galactic Skordis for by at ma search of MOND the order of as 7 relevance over emer galaxies the should spherical strengthen MOND for that observations a view MOND FJR the recent the reinforces Comparing [18] Milgrom[17]. observations (FJR), curve relation Jackson Faber the 2] MN 3] mretGaiy[1 n O 3] tlredistances larger large a At accommodate include can [32]. These which MOG acceleration [28]. and like A1689 MOND [31] cluster a Gravity have in Emergent function [30], r DM GMOND central the the [29], match in with to profile formation used temperature cluster are the for galaxy accounting Modelling in especially temperature. observed the omnrsos otennNwoinrttoa uv faga a of curve rotational non-Newtonian the to response common A aibeMdfidNwoinMcaisI:Broi ul Fi Tully Baryonic II: Mechanics Newtonian Modified Variable hc pr rmteNwoinaclrto,gvsrs oa to acceleration rise MOND gives the acceleration, which Newtonian the in from apart which ASnmes ?? numbers: PACS TRprit u otecnevto faglrmmnu an momentum value angular l acceleration of scale MOND conservation few phenomenological the a the to to to shrinks due R shell realisti persists Fisher mass Within BTFR Tully outer Baryonic the dominat As the growth. acceleration acceleration. obeys overdensity MOND-like velocity the the rotational where of shell value end a at an in growing acceleration, to this lead in cloud can protogalactic size Way Milky eetyw n igemti ouinfrapitms res mass point a for solution single-metric a find we Recently eateto lcrcladEetoi niern,Univ Engineering, Electronic and Electrical of Department a 0 ovr nsrnt codn oteevrnetweems discr mass where environment the to according strength in vary to a 0 a erdc ag cl tutr iia oteCMrsls In results. CDM the to similar structure scale large a reproduce can a Dtd eray1,2021) 11, February (Dated: 0 INTRODUCTION srpae ytecsooia ceeain esuya study We acceleration. cosmological the by replaced is a .C Wong C. C. 0 erdcsteBF,seas aayMGuh1][6,and Famaey-McGaugh[15]-[16], also see BTFR, the reproduces , a 0 ncptnilt apnbroi cutcoscillations acoustic baryonic dampen to potential onic tlt time. late at a ntdsi assae uaoe.a.[9,further [19], al. et. Durazo scale, mass in gnitudes dtebro as ntemdfidD theory DM modified the On mass. baryon the nd 0 go,Nse 2] te oie rvt theories gravity Modified Other [27]. Nusser egion, tr omto iuain usr[4 hw that shows [24] Nusser simulation, formation cture hc h rwho nua momentum angular of growth the which 2] nEMOND, In [28]. hs hr Mcnitrc oesrnl with strongly more interact can DM where phase umr,aval ODlk hoysesto seems theory MOND-like viable a summary, n eto swehrsc ceeaincnb found. be can acceleration such whether is uestion lto BF)wt mle MOND smaller a with (BTFR) elation lse 2][6,i sfudta ihrmore either that found is it [25]-[26], cluster y tro miia ayncTlyFse Relation Tully-Fisher Baryonic empirical of tter nt itne,terttoa velocity rotational the distances, ength efo n ibemdfidgaiyter.The theory. gravity modified viable any from ge ODt xli h ihaosi ek [16]. peaks acoustic high the explain to MOND f a ayi htis”etna”dnmc requires dynamics ”Newtonian” its that is laxy diinlMN-ieacceleration MOND-like additional n rvt sCnomlGaiy[2-1] which [12]-[13], Gravity Conformal is gravity g problems modelling CDM some mitigate uld 0 ODuigsmlto ean problematic remains simulation using MOND otipratyteicesnl difficult increasingly the importantly most d soe etnadteotrmass outer the and Newton over es elc farltvsi hoylmt h use the limits theory relativistic a of lack he t oiatnnNwoinacceleration, non-Newtonian dominant a ith yafco f4wl encsayt produce to necessary be will 4 of factor a by dn na xadn nvre[1], universe expanding an in iding ,Cpaa wr aaisaon Centaurus around galaxies dwarf Coplanar ], ODrmisapeoeooia theory phenomenological a remains MOND esit,teoedniystops overdensity the redshifts, c diinltnbe”osat acceleration. ”constant” tunable additional aia rcini aayfrain[] the [9], formation galaxy in friction namical riyo ogKn.H.K. Kong. Hong of ersity lgclstain.Tehg gravitational high The situations. ological MN,aGnrlsdvrin fMOND of versions Generalised a EMOND, dD rdcin ihglci rotational galactic with predictions DM nd h ODaclrto grows acceleration MOND the d n htacaso ee hoyproposed theory TeVes of class a that ing toa ceeainwe h Newtonian the when acceleration ational es[] otne omtvt h hunt the motivate to continues [5], ders ohEODadEegn Gravity Emergent and EMOND both , a 0 nrae tlreNewtonian large at increases hrRelation sher pnisappear epancies MOND a 2

Emerging studies of Early type galaxy (ETG) rotational curves open up a new perspective on gravitational acceleration of large structures during cosmic evolution. Recently, Genzel et.al. [34] and Lang [35] show that early type galaxies have Newtonian rotational curves outside its scale length r0, which in fact fall oﬀ faster than the (Newtonian) exponential disk from 2r 5r , where 5r is the galactic disk size. By modelling DM density proﬁle 0 ∼ 0 0 based on the NSW model, the authors of [35] obtain the following equation for the rotational speed vrot, speed due to disk mass vdisk and speed due to dark matter vDM

2 2 2 2 r vrot(r)= vdisk(r)+ vDM (r) 2σ ( ), (1) − r0

2 2 where the (1-dimensional) σ is the velocity dispersion usually taken to be vrot = 3σ and the last term in the Eq. (1) represents pressure support inside the disk. Using simulations, the authors in [35] reproduce the observations by using maximal baryonic density (no DM allowed) inside the disk and including pressure support term. This work suggests that ” at z 2 there is no need for dark matter inside the galactic disk upto 5r0 ”. Work by Simons et.al. [36] demonstrates that∼ from z 2 to z 0.2, galactic rotational speed starts growing at a rate correlating with galactic mass growth and a rapid∼ decrease∼ of velocity dispersion σ. This is independently confirmed by Kassin et. al. [37]. Pure disk galaxies [38] observations also show that from z 1, a central exponential disk profile is maintained without major mergers. Their central surface brightness remain ∼unchanged over 1 >z> 0.4, but over 0.4 >z> 0.02 the stellar disk mass increases upto 40% while the central brightness and scale length remain unchanged. This indicates that after z < 0.4 a significant amount of gas from large radii has moved onto the central disk. However, this gas cannot get through into the stellar disk and is kept in the region r0 2r0. Continuous cold gas accretion is known to be a crucial ingredient in the evolution of spiral galaxies, see for∼ example [38]-[40]. The gas accretion time is almost an order of magnitude larger than the free fall time under a Newtonian acceleration, it is generally believed that this gas comes from cold streaming of the cosmic filaments. However, if the gas accretion comes from the same protogalactic cloud, one would expect that the dynamics of the accretion gas would be very different to that of Newtonian collapse.

The collapse of HI gas from large radii onto a stellar disk is also observed in the DiskMass Survey [41]-[43]. At early time the mass profile of HI gas is low inside 5r0, but increases to a peak outside 5r0 and then becomes flattened. The rotational curve also flattens outside 5r0 where the HI mass profile is flat. At late times, the HI mass profile is seen to have moved inside the 5r0 distances and peaks strongly around 2r0.

Even more specifically, [44] finds that for ETG at small radii outside r0, the amount of DM required is small ( 10%). However, for large radii outside 5r0, the rotational curves follows a BTFR with very tight scatter [45].∼ This BTFR matchs the late time BTFR of spiral galaxies at 60 70% of maximal mass to light ratio. − In MOND term, this is equivalent to requiring a smaller MOND accleration a0 at higher reshift and at larger distances. The above observations show that for Early type spiral galaxies, a nearly Newtonian central stellar disk develops first while its outer MOND-like rotating gas at radii larger than 5r0 moves very slowly towards the stellar disk.

In studies involving galaxies, the background metric is generally assumed to be the Einstein-Straus (Swiss-cheese) model [46] or the McVittie metric [47], in which a central mass region is described by the Schwarzschild metric, while outside of a junction the relevant metric is that of the Friedmann Robertson Walker (FRW) metric. Neither of these solutions is without criticism. It is known that a Einstein-Straus (Swiss cheese) solution does not exist in cylindrical coordinates, and the McVittie metric has a pressure singularity at small distances. A review of the difficulties in both approaches can be found in [48]. In strong lensing applications for galaxies, the debates between these two metrics remain ongoing [49]-[53]. However, these solutions are not necessarily the only choices. In a previous work [1], we find a single ”inhomogeneous” metric from solving Einstein equations for a point mass residing in an expanding FRW background, in which a MOND-like potential (acceleration) at slow speeds arises. This potential has no tunable parameter and depends strictly on the central baryonic density and the cosmological background density. In [1] we find that combining Newtonian and non-Newtonian acceleration (which we name ”VMOND” acceleration) amounts to the MOND scheme that interpolates between Newtonian and non-Newtonian acceleration. In our case, the MOND acceleration a0 is replaced by the cosmological acceleration. The VMOND acceleration can play a similar role of DM in density perturbation growth and lift the higher CMB acoustic peaks. In [54], Turnaround radius bounds from the Einstein-Straus and McVittie metric appear to be violated by very large structures, while we find that [55] the metric with VMOND potential provides a larger turnaround radius bound that can accommodate these structures. We note that the late time VMOND acceleration is negligible outside scale length distance of Milky way like galaxy, but is comparable to MOND only at a distance scale much larger distances. This is in fact desirable since the early time 3 galaxy rotational curves outside the scale length require a pure Newtonian potential as discussed above. Therefore our question is whether the late time BTFR has its origin from the dynamics of a protogalaxy in a VMOND potential.

In [19], the velocity dispersions-mass relations show that Ellpitical galaxies typically has a central Newtonian central region characterised by a radius scale RM above which the velocity dispersions become flat and follow a non-Newtonian FJR. Galaxy formation modelling starting with a spherical overdense cloud at high redshifts in a MOND acceleration produces analytically a FJR, Milgrom[17]. In simulation, MOND produces a characteristic radius for a central Newtonian region and above which the velocity dispersions-mass relation is that of FJR, Sanders[56]. Before the advent of DM, one of the most favoured models for galaxy formation in Newtonian potential is the Mestel theory for an initial spherical uniform density (or constant mass) cloud with constant angular velocity to collapse into a disk. Typically, there is a central core which collapses earlier to form an expontial disk with a constant angular speed while a relatively small amount of late arrival gas, which comes from a constant mass cloud, collapses to a disk with constant rotational speed. The model is found to work well by Crampin and Hoyle [57] and Freeman [58], but to be valid upto multiple scale length distances an unusually high mass to light ratio (M/L) is required. Without DM, the observed extended flat rotational curves of spiral galaxies require the evolution of an outer cloud in a non-Newtonian aceleration. We envisage that for a large protogalactic cloud, there is a central region which obeys a nearly Newtonian gravity and collapse following the Mes- tel theory, while the outer cloud can follow a different and much slower path of collapse under a VMOND acceleration.

In section 2, we recall the main features of the new metric and calculate the orbit of a single particle in the VMOND potential. In section 3, we study in VMOND potential a model protogalaxy that breaks away from its cosmological background at a large enough angular momentum. In section 4, we examine whether an overdensity having a power-law dependence on distances can lead to a central disk and the late-time state of the outer non-Newtonian mass shells. The last section is a discussion and conclusion.

2.1: THE MODEL

In a previous work, we ﬁnd a Tolman-Lemaˆitre metric

∂β 2 ds2 = e2β(̺,τ)d̺2 e2β(̺,τ)dΩ2 + c2dτ 2 (2) − ∂̺ − that describes the motion of test particles in the vicinity of a central point-mass M placed in an expanding ﬂat cosmological background, with proper time τ, comoving distance ̺ and dΩ2 = dθ2 + sin2θdϕ2. The physical radius is given by

τ 2/3 3 r = eβ(̺,τ) = a(τ) C ̺ C a(y)−3/2dy ; C = √2GM, (3) 1 − 1 1 2 Zτ0 a(τ) is the scale factor and for a(τ) = 1 one recovers the Lemaˆitre form of the Schwarzschild metric. Here r a(τ)̺2/3 for large ̺. ∼

In curvature coordinates this metric takes a familiar form dr2 ds2 = Zc2dt2 r2dφ2, (4) − Z − where

2GM H(z)r 2 Z =1 , (5) − rc2 − c2 r where G, M, c, r are the Newton’s constant, the central point mass, the speed of light and the radial distance respectively, H(z) is the Hubble constant at red-shift z. The slow speed energy equation is given by

1 h2 E = r˙2 + U H(z), r . (6) 2 2r2 − 4

1 2GM 2 GM 1 U H(z), r = H(z)r = H(z)√2GMr + H(z)2r2, (7) 2 r − r − 2 r where E is the total energy of the particle, h is the angular momentum per unit mass. The radial acceleration has the form h2 GM GM a¨ r¨ = 3 2 H r. (8) r − r − r 2r − a We call the sum of the Newtonian and the non-Newtonian gravitational acceleration as the ”VMOND” acceleration. These equations of motion can be directly obtained from the Tolman-Lemaˆitre metric [1], the metric in curvature coordinates provides a convenient perspective from the Schwarzschild metric view point.

We can read oﬀ some notable features speciﬁc to this metric.

In a particle free fall, at speciﬁc angular momentum h = 0 and energy E = 0, the energy equation reduces to the simple relation

2GM r˙ = Hr . (9) − r r which depicts that the particle will follow the Hubble expansion at large distances but will follow a Newtonian path at small distances. For the Schwarzschild de-Sitter metric, the above relation is not automatic but requires a swiss cheese construction. At the distance from centre

2GM 1/3 r = (10) H2 the metric becomes Minkowskian. This is also the radius wherer ˙ = 0 from Eq. (9) where the metric is isotropic and homogeneous. In the Lemaître-Tolman metric in [1] with time τ, and r(τ) as defined in [1], at rc, we have dH/dτ =0 where the Junction conditions [59] dr(τ) = 0; P =0 (11) dτ is satisfied, where P is the pressure which is zero at (Minkowskian) empty space.

The Hubble parameter at redshift z is given by 8πG 8πG 8πG H2(z)= ρ = ρ + ρ + ρ = H2 Ω (1 + z)4 +Ω (1 + z)3 +Ω , H2 = ρ , (12) 3 H 3 r m Λ 0 r m Λ 0 3 c where for i = r, m, Λ, ρi are cosmological background densities of radiation, pressureless matter which we assume is essentially baryonic matter and dark energy respectively. H0 and ρc are the Hubble parameter and the critical density of the present epcoh. For ρi,0 denotes the density of i compoent at the present epoch, we have the density parameters Ω = ρi,0 . Also the acceleration equation i ρc a¨ 4πG 8πG = 2ρ + ρ + ρ . (13) a − 3 r m 3 Λ M 4π 3 We can rewrite Eq.(8) in term of local ”central” baryonic matter density ρb(r)= V (r) , V (r)= 3 r

h2 4πGr a¨ r¨ = ρ (r)+ ρ (r)√ρ + r. (14) r3 − b b H 3 a p As we are interested in a protogalaxy in the redshift range 15

ρb δ(r)ρm = ρb(r), or δ(r)= . (15) ρm 5

As the uniform radiation energy density ρ and temperature T take the form ρ T 4 and the adiabaticity of the γ r ∝ perturbations on superhorizon scale δ( nr ) = 0, n ,n which are photon and baryon number respectively, which leads nb γ b to δρr δT 4 δρm 4 δγ = =4 = = δ. (16) ρr T 3 ρm 3

Observationally, temperature ﬂuctuations from recombination is related to δb by δT 1 = δ (17) T 3 This δ is taken as the initial value of baryon overdensity δint in its subsequent evolution after recombination.

From Baker [60] Eq.(11). the mass density in an enclosed spherical comoving volume radius ̺ is given by ∂ Λ ∂β e3β β˙2 =8πT 4 e3β. (18) ∂̺ − 3 4 ∂̺ 4 where T4 = Gρ(̺) is the time-time component of the momentum tensor, ρ(̺) is the mass density in comoving radius ̺. Calculate in terms of coordinate r using Eq.(3) gives local mass density, see [1]

1/2 ρ = 1 + ∆(δ) ρm, ∆(δ) = (δ + δ ). (19) To obtain the baryonic overdensity growth equation in comoving coordinates, instead of usual total local mass density

ρ = (1+ δ)ρm, (20) in this metric, the local total mass density is now given by Eq.(19). The overdensity growth is given by the familiar equation with δ replaced by ∆(δ), a˙ 1 ∆(¨ δ)+2 ∆(˙ δ)= 2 Φ=4πG∆(δ)ρ . (21) a a2 ▽̺ m For δ 1, Eq.(21) has a growth solution ∆ a(t) so that δ a(t)2. This result can also be seen in a top hat ≪ ∝ ∝ 1/2 2 spherical collapse model where the dominant acceleration is from the term δ ρm [1]. This δ a growth rate can account for the eﬀects attributed to DM to improve δ a(t) growth. We note that in a MOND∝ acceleration approach, the overdensity δ also grows as a2, Nusser [24].∝ This same growth rate is also used in elliptical galaxy formation simulation under MOND from an initial overdensity, Sanders [56]. It is important to stress that in this inhomogeneous metric, the matter perturbation in principle is allowed to grow beyond a linear regime. In the case of a protogalactic cloud experiencing no external angular momentum disturbances, δ is expected to grow near unity as a2 and after which it will continue to grow as a(t).

When perturbed by suﬃciently strong angular momentum, the overdensity δ can stop growing. After that point, the central mass M(r) is ﬁxed and the central mass density ρb(r) = M(r)/V (r) is no longer linked with the cosmological background but increases as r decreases and is only used to contrast with the cosmological density ρH in which the dark energy component will become important at low redshifts.

For a central point-mass epoch, from Eq.(8) the Non-Newtonian radial acceleration is given as GM GM H(z)2r GM r¨V M = H(z) = 2 = 2 a0(z, r) (22) − r 2r −r r r 2 −r r p where a0(z, r) is the cosmological acceleration at distance r from the mass centre. This is reminiscient of the −10 2 MOND accelerationr ¨ = √gN a0 where gN is the Newtonian acceleration and a0 = 1.2 10 m/s is the ﬁxed phenomenological MOND− acceleration. ×

Similar to MOND, if a0(z, r) is larger than the Newtonian acceleration, we will enter into a non-Newtonian −1 −1 acceleration. For late time galaxies, a quick check using a high value of H0 = 75 kms (Mpc) and a large distance outside the stellar disk of the Milky Way at r = 50kpc, the cosmological acceleration is a (H , r)=4.6 10−15m/s2. 0 0 × This a0(H0, r) is too weak to match the MOND acceleration. However, as the above ETG discussion show, the non-Newtonian rotating gas appears at an epoch with larger H(z) and at a distance that is larger than r > 50kpc, where a0(z, r) is much larger than a0(H0, r). To see whether it is plausible to reproduce the phenomenological MOND acceleration, we consider a protogalactic cloud. 6

2.2: The particle orbit in the VMOND potential outside the overdensity growth region

Subsequent to a protogalaxy break-away, we will need to deal with orbits of a point mass embedded in an expanding background. Before we proceed, it is thus useful to study of the orbit of a slow speed particle of unit mass around a fixed central mass M with specific angular momentum h under the influence of a central potential with polar coordinate (r, ϕ) (e.g. see [61]]) that satisfies r r−1dr ϕ = (23) 0 2 2 r 2 (E V )r 1 Z h − − q where E is the total energy. Eq.(23) becomes r r−1dr r r−1dr ϕ = = . (24) r0 2 2 2 2 r0 2 (GM/r + E H(z)√2GMr + H(z) r /2)r 1 2 GM ρH (z) 2 Z h Z 2 ( (1 ) + E) 1 − − h r − ρb(r) − q r q where ρH = ρm + ρΛ is the background cosmic density ranging from the matter dominant to the dark enery dominant epoch and ρb(r) is the baryonic density at r given by M ρ (r)= . (25) b 4π 3 3 r We stress that we are no longer considering a protogalaxy growth, but only a point particle in a VMOND potential starting from matter dominant epoch to the present time.

We set

H2r3 2 ρ (z) 2 α = 1 = 1 H , 0 α 1. (26) − 2GM − ρ (r) ≤ ≤ r s b 2 h El 2 We follow the usual parametrisations by taking a length scale l with 2GM = l, E0 = GM and eccentricity ǫ = 1+4E0.

Rearranging Eq. (24) we obtain r 2lr−2dr r 2lr−2dr ϕ = = (27) 2 2 r0 r0 Z 1 2 1 2l Z 1 1 2l 2l 2α ǫ 1+ 2 (α 1) 2 α ǫ 1 2 2 + ǫ − − ǫ − r − ǫ − ǫ r ǫr ǫ s s Setting 2l ǫx = (28) − r we obtain dx ϕ = (29) Z 1 2α 2 1 2 x x − ǫ − ǫ − s Although α(r) depends on r and z, its variation is based on a cosmological time scale. For an orbital path in which α(r) changes very little, we can approximate (29) by taking α = α0 = α(rkep) (constant) in the integral, and the linear term in x is ﬁnite as long as ǫ> 0. We obtain

α2 1 1/2 2l 1+ 0 − ǫ sinϕ α , (30) ǫ2 ≈ 0 − r That is, in Keplerian form 2l′ r = (31) 1 ǫ′sinϕ − 7 where

l ǫ α2 1 1/2 l′ = , ǫ′ = 1+ 0 − . (32) α α ǫ2 0 0 For a particle with positive energy E > 0 at high redshift and distance r where α is large, the particle orbit can start from an elliptical orbit and moves into a hyperpolic orbit.

In the dark energy (cosmological constant) dominant epoch, H(z) = H0 is nearly constant (α 1). We consider speciﬁcally the Earth’s orbit around the Sun with radius at 1AU =1.496 1011m, we have ∼ × −26 −3 −26 ρH ρc 1.1 10 kgm 1.1 10 −23 = = × 4π = × −4 7.71 10 , (33) ρ (r) ρ (r) M⊙/ (AU)3 1.42 10 ∼ × b b 3 ×

ρ H =8.78 10−12, α 1 1.75 10−11. (34) ρ (r) × ≈ − × r b VMOND correction to the closest approach is an O(10−11) eﬀect, while the Viking observation uncertainties for Earth and Mars orbits are 100m and 150m respectively, which are at O(10−9) [62]-[63].

3.1: INITIAL CONDITIONS OF A LARGE PROTOGALACTIC CLOUD

After recombination, at T = 3000K, under Newtonian gravity (with dark matter) the Jean Mass MJ estimate is 6 5 6 MJ < 10 M⊙ so that high density objects like globular clusters with MJ = 10 10 M⊙ in principle can form at z 100. In a Newtonian potential, it is generally assumed that when the baryonic∼ density perturbation reaches O(1) the∼ protogalaxy will start to break away from the expanding universe. It is expected that massive galaxies begin to form at redshﬁt z 15 to 50 Barkana-Loeb [64]. The earliest observed galaxy is GN-z11 at z 11.09 [65] with mass 9 ∼ ∼ 10 M⊙ and the detection of redshifted 21cm absorption at z 17 provide important signposts on the plausible ∼ 9 11 ≈ break-away epoch for galaxies with mass 10 M⊙ 10 M⊙. ∼ 9 11 To simulate the formation of galaxy with mass 10 M⊙ 10 M⊙, the authors in [66] use an initial uniform −23 −3 ∼ cloud density ρint = 1.11 10 kgm in a MOND background potential. A Milky Way formation simulation × −23 3 under a pure Newtonian potential is done using initial cloud density at ρint =7 10 kg/m [67]. We want to pick as guard-rail a Milky Way-like protogalaxy with its observed masses M(r), assuming× that the overall strucure is not altered by a past encounter with Andromeda in the local group. Within the Milky-Way scale length r0 = 8.5kpc, 10 the observed stellar mass is found to be M(r0)∗ = 6.07 10 M⊙ [68]-[69]. Assuming pure Newtonian potential 2 × 10 and using the relation h = 1.1GMr0, the observed rotational speeds at r0 gives M(r0)=8.4 10 M⊙. In [70], outside the Milky Way stellar disk, an exponential disk of HI gas is observed up to 60kpc and [71]× uses in simulation 10 10 10 a dark gas ring of 2.4 10 M⊙ (baryonic gas 1.2 10 M⊙) and thus total baryonic mass is 9.5 10 M⊙. We ∼ × × ∼ × therefore take the Milky Way baryonic mass estimate Mdistance at to be

10 10 M =8.4 10 M⊙,M = M =9.5 10 M⊙. r0 × 50kpc 50 ×

The protogalaxy mass is then chosen to be that of M50kpc. Our primary interest is the BTFR in the region of 2r 6r , where the present galactic baryonic mass density ρ at r and r = 50kpc respectively are given by 0 ∼ 0 0 50 ρ =2.42 10−21kg/m3, ρ =1.34 10−23kg/m3. (35) r0 × 50 × GM Qualitatively speaking, we note that the VMOND scheme involves comparing the Newtonian acceleration r2 with 1 2 the cosmological acceleration 2 H(z) r at distance r, which is equivalent to comparing a central baryonic denisty ρb(r) at distance r with an uniform cosmological density ρH (z). Given ρ50, if the initial protogalactic cloud radius is twice the size of its late time radius, at that point the protogalactic cloud density is O(10−24)kg/m3 which would be surpassed by the early time cosmological density say at z 20. Above this redshift, we will expect the non-Newtonian acceleration dominates in the gravitational acceleration.∼ 8

We use the simple analytic treatment of gas cloud at constant uniform density which is represented by spherical shells moving under the eﬀect of cosmic expansion, Gunn and Gott[76], White[77] -Nusser[79]. The mass shells move independently without friction between neighbouring shells. One point on the mass shell can be represented by a point particle at a distance ri from the centre in a VMOND potential.

We start by considering the largest spherical shell, from Eq.(6)-Eq.(7)

1 1 2 E = r˙2 √δ 1 H2(z)r2. (36) 2 − 2 − where M(r) is the enclosed mass inside a sphere at radius r. For δ = 0, E = 0 and a shell at radius r will simply follow the Hubble ﬂow. As δ grows, E is not necessarily constant. For example, in a Newtonian potential (keeping only the δ term), E grows as GM(r)/r. However, Eq.(36) admits a E = 0 solution, so thatr ˙ 0 as δ 1. Since the startingr> ˙ 0 is outward moving,− typically whenr ˙ reaches 0 and reverses direction, the overdensity→ breaks→ away from its cosmological background. We recall that with angular momentum, the radial acceleration equation Eq.(8) in terms of densities in matter dominant epoch is

h2 4πGρ (z)r r¨ = 1+ δ + √δ m . (37) r3 − 3 Keeping h2 = 0 for spherical overdensity, we note that δ represents the Newtonian contribution and √δ the non New- tonian contribution to the acceleration, so that when δ √δ the contribution from non Newtonian term dominates over the Newtonian and if δ √δ the contribution we≪ are back to Newtonian dominance. Therefore, in a growing overdensity, the size of density≫ perturbation δ determines whether or not the acceleration is Newtonian dominant.

3.2: Eﬀect of a large Angular momentum

During the overdensity growth, the overdense region can receive an angular momentum due to tidal fields from its neighbours. A common assumption is that angular momentum behaves as a small perturbation and does not affect the overdensity growth (in this case h2 = 0 is kept in Eq.(21). In the small δ regime, the angular momentum growth based on Tidal fields of neighbours is modelled by Peebles [72] in which the time dependent angular momentum per unit mass is given by the second order perturbation which grows as

h(t) t5/3. (38) ∝ so that the angular speed grows as

h4 a(t), (39) r4 ∝ while White [73] obtains a ﬁrst order perturbation from a ﬂatten sphere which grows as

h(t) t. (40) ∝ that the angular speed grows as

h4 1 . (41) r4 ∝ a(t) Casuso and Burkert in [74] -[75] propose that an expanding void can also provide a source of the protogalaxy’s angular momentum. In this scenario, the angular momentum transfer from the void causes the gravitational collapse and the break-away of the overdensity. In this study we consider the possibility that the growing angular momentum can catch up to the gravitational acceleration, and thus plays a critical role in ending the overdensity growth.

We consider mass shells of spherical overdense cloud which intersects the mid-plane perpendicular to the rotational axis. These mass rings are applicable to an overdensity in either spherical shape or as a thin disk. A typical mass ring has the equation of motion Eq.(37) with h > 0. At some δ, the protogalaxy picks up an angular speed h2r−4 which for second order perturbation starts to grow either as a or for a larger ﬁrst order perturbation falls 9

−1 −3 as a . As the background density ρm a continues to fall so that the gravitational acceleration due to central density perturbation in Eq. (37) grows as∝ 4πG 4πG δ + δ1/2 ρ = δ a−1 + δ1/2a−2 ρ , (42) 3 m int int 3 int where δint and ρint are initial overdensity and background density at recombination respectively. Initially, at small a(t) the central gravitational acceleration dominates over the angular speed. As a(t) grows, the large gravitational acceleration (a−2) term drops faster than the angular speed (either a or a−1). While the angular speed from the Peeble model (h2r−4 a) can catch up to the gravitational acceleration, for the White model (h2r−4 a−1), a ∝ ∝ signiﬁcantly large initial angular speed comparing to δint is required for the angular speed to catch up to the central gravitational acceleration at large a. In this work we assume this is the case and at some distance ri for the largest shell and zi, one will reach the point h2 4πGρ (z ) ϕ˙ 2 = i = δ + √δ m i , (43) r4 3 i here δ then takes on a ”ﬁxed” value.

For a spherical cloud, if angular velocity is uniform, the angular velocity is the same for mass shells at diﬀer- ent radius r from the rotating axis. Eq.(43) is then true for every horizontal slice of the sphere in cylindrical coordinates. From the overdensity evolution equation Eq.(21), the overdensity mass in the spherical mass shell also 4π 3 freezes at M(ri)= 3 δρm(zi)ri .

3.3: The turnaround point of the largest mass shell

After the overdensity stops growing, the energy equation for the mass shell changes to

1 h2 1 2GM 2 GM(r) ρ (z) 2 r˙2 + = H(z)r + E = α + E, α = 1 H , (44) 2 2r2 2 r − r − ρ (r) r s b where E and M(r) are now fixed by δ and H(zi), also ρb(r) no longer linked to the background matter density. The turnaround is the point ρH (z)= ρb(r) (α = 0) which represents a point where there is no gravitational potential. This is the point the VMOND metric becomes Minkowskian. Initiallyr ˙ > 0 remains outgoing from the mass centre and − 1 α = (1 δ 2 )2 is large. As ρ (z) decreases with z and α also decreases. From Eq.(44), bothr ˙ and h2/r2 decrease. − H There exists a solution with a maximum radius rta where h2 α =0, r˙ta =0, 2 = E > 0. (45) rta This turnaround point is also regarded as the break-away point of the overdensity from the Hubble flow. As redshift reduces further, α increases again, the mass shell starts to move towards the mass centre so that r < rta. However, we stress thatr ˙ is small near the turnaround point, the change in α depends on the reduction of cosmological background density ρH (z), which is very slow, the mass shell can take a long (cosmic) time to come back to the mass centre. If the system eventually approaches viralisation, the mass shells will cross one another and mix to reach an equilibrium state over a long relaxion time. In this work, we assume that the large non-Newtonian mass shells remain in collisionless collapse and has not come to significant mixing.

3.4: Early time Baryonic Tully Fisher Relation

Next we focus on the outer-most mass shell. From the conservation of angular momentum, the angular momentum 2 h (ri) at the point of overdensity freeze remains unchanged. For simplicity, we consider that this mass shell returns to radius ri after reaching the turnaround point, where M(ri) is the overdensity mass and Eq.(43) remains true. For the mass shell intersecting the mid-plane, Eq.(43) in the δ √δ non-Newtonian region leads to the rotational speed ≪ V (ri) h4 4π V (r )4 = i = ρ (z )r (1 + √δ)2 GM(r )= a (z , r )GM(r ). (46) i r4 3 m i i i 0 i i i i 10

Here the ρm(zi) is set at the growth-freeze redshift value zi due to angular momentum conservation and the Newtonian acceleration part is represented by the factor √δ. Eq.(46) already shows that the mass shell rotational speed follows the BTFR, with an eﬀective MOND acceleration given by 4π a (z , r )= ρ (z )r (1 + √δ)2. (47) 0 i i 3 m i i We note that in Eq.(46), a small variation of δ does not change the overall power law behaviour between V and M. It is also possble to write Eq.(46) in the form

1 GM(r ) V (r )2 = 1+ i . (48) i √ r δ i For large δ, a variation of δ maintains the Newtonian nature of the acceleration with some varying invisble matter eﬀect. However, for small δ, the δ−1/2 term dominates the bracket and disrupts the Newtonian power-law behaviour, i.e. with a strongly varying dark matter mass proﬁle. Thus variations of a small δ is better described by the power law in Eq.(46). As δ varies from a small value to a large value, we expect to see V x M with n varies from x 4 to −1 ∝ ∼ x 2. As δρm a , small δ is associated with a larger mass. For larger δ one could obtain an exponent n which is ∼ ∝ 9 slightly smaller than 4. In Lelli [80], for smaller mass (< 10 M⊙) galaxies, the observed BTFR exponent is tending towards x< 4.

Consider the largest mass shell at ri which is far away from r50 but eventually shrinks to galactic disk scale r50 at 50kpc by a contracting factor n such that

ri = nr50 (49) for n> 1. From the conservation of the angular momentum of the mass shell, the rotational speed will increase as its radius shrinks. After mass shell radius contracts to 50kpc, from Eq.(46) the rotational speed continues to obey the BTFR but with an acceleration larger than a0(zi, ri) having a higher rotational speed

4 4 4 4 V M V (r50)= n V (ri)= n a0(zi, ri)GM(ri)= a0 (r50)GM(ri), (50) so that 4πG aV M (r )= a (z , r )n4; aV M (r )=(1+ √δ)2 ρ (z ) 50kpc n5. (51) 0 50 0 i i 0 50 3 m i × × where aV M (r ) should be comparable with the phenomenological MOND acceleration a =1.2 0.3 10−10ms−2. 0 50 0 ± ×

We note that after the mass shell radius ri shrinks to the size r50, we have 4π M = ρ r3 , ρ = δn3ρ (z ). (52) 50kpc 3 50 50 50 m i

ρ50 for our system is given in Eq.(35), which is an observable density. Using the present day critical density ρc = 1.11 10−26kg/m3, the background matter density at redshift z , where Ω =0.05 is given by × i m 3 ρm(zi)=(0.05)(1 + zi) ρc. (53)

This is the same cosmological assumption used in Sanders [56] to simulate elliptical galaxy formation in MOND acceleration, where the overdensity is assumed to stop growing when the ”non-Newtonian acceleration” is twice the size of the cosmological acceleration.

−4 ∆T −5 At combination with z = 1080, we consider a believable δ = δint 1.65 10 ( T 5.5 10 ). Since δ 2 ∼ × ∼ × grows as a(t) , we can calculate the overdensity δ at redshift zi by the scale factor growth to obtain

1081 2 13.88 δ = δ , √δ = (54) int 1+ z 1+ z i i V M so that we can ﬁx n and obtain a0 (r50) from Eq.(51). 11

(We should point out that the Hubble tension and the KBC void question are under investigation and could incur upto 10% uncertainty in H0 value but it should not alter the qualitative picture of this analysis.)

We focus on an observationally viable (from 21 cm radiation) redshift range 15

1 1 c/H a = cH = H2r γ; γ = 0 =2.67 104, (55) 0 6 0 2 0 50 3 50kpc × × where the factor γ speciﬁes the MOND observed value. Eq.(51) takes the form

1 aV M (r )= H2r γ(z ), γ(z )=(1+ √δ)2Ω (1 + z )3n5 = 1.45 102 (1 + z )4/3 (1 + √δ)2. (56) 0 50 2 0 50 i i m i × × i −24 3 If δ freezes at redshift zi = 17, we have √δ = 0.771 (δ = 0.594), n = 1.90 and ρm(zi)=3.23 10 kg/m so that γ(z )=2.14 104 × i × V M a0 (r50)=0.803 a0. (57)

−23 3 For δ freezes at higher redshift zi = 38, √δ = 0.355 (δ = 0.126), n = 1.47 and ρm(zi)=3.29 10 kg/m , γ(z )=3.51 104 so that × i × V M a0 (r50)=1.31 a0. (58)

At zi = 23, we have

V M a0 (r50)= a0. (59)

The above demonstrates that given a realistic δint value at recombination, within a large range of redshifts (17

4.1 Central disk

Similar to the nearly spherical galaxies, spiral galaxies also has a Newtonian core which appears as an exponential disk. We recall that the gravitational collapse of an uniform density, uniform angular speed spherical cloud in a Newtonian potential can be analytically described by Mestel [81] and Freeman [58] to produce an approximately exponential central disk. This is also supported by simulation [67]. However, the VMOND acceleration analysis above based on an uniform density overdense sphere does not automatically produce a central Newtonian dominant region where overdensity reaches O(1). From Eq.(44) that δ 1 leads to a large ρb(r) after turnaround which ensures that α > 0 in subsequent redshﬁts and the total energy can≥ take the value E < 0 which guarantees a Newtonian gravitational collapse. Although our primary aim is to understand BTFR at large radii, it is interesting to see how a central Newtonian region may arise.

A plausible scenario is suggested by Nusser in [27] that the size of overdensity has a power law dependence on r.

r S δ(r)= δ i (60) r where ri is the maximum size of the cloud, and the exponent S > 0 is a constant. Here δ(r) plays the role of an eﬀective δ at r. The acceleration equation Eq.(37) becomes

h2 δ δ1/2 4πG r¨ = + ρ r. (61) r3 − (r/r )S (r/r )S/2 3 m i i 12

S The factor δ(r)= δ(ri/rM ) > 1 leads to a central Newtonian dominant region inside a radius rM region, where

r δ1/S r . (62) M ≤ i

For ﬁxed S > 0 and ri, we see that RM starts to increase as δ increases from δint. This could lead to Newtonian gravitational collapse in the small r region at very high redshifts. We consider the scenario that when outer overdensity stops growing δ < 1, the overall overdensity also stop growing. The Newtonian region inside RM has (δ(r) > 1) can however remain at its slowr ˙ region while the angular speed is transported from the outer region. The mass density of the protogalactic cloud is no longer uniform. Insider rM the mass density is given by

r S δ(r)ρ (z )= δρ i (63) m i m r and the central mass upto shell at r is

r r 4πδρ M(r)=4π r2ρ(r)dr =4πρ δ(r)r2dr = m rS r(3−S). (64) m (3 S) i Z0 Z0 − We can estimate S using the observation that a Newtonian central mass M(r ) is 90% of the total protogalaxy M ∼ mass at ri, so that

3−S M(ri) 10 ri − 3−S = = = δ S (65) M(r ) 9 r M M

For z = 17.8, δ =0.5 we obtain S 2.6 , the central Newtonian sphere radius is given by rM =0.768 ri = 74.88kpc. The outer non-Newtonian gas shell≈ is about 23% of the total spherical cloud radius.

The overdensity’s power-law behavour in r−S is one assumption that can lead to a central Newtonian region. Other decreasing functions such as an exponential function

δ(r)= δAe−r/rB , (66) where A and rB are constants, could lead to a large portion of the central Newtonian sphere having nearly uniform density, which can make the Mestel model calculation more meaningful. We see that for an overdensity δ(r) with a reasonable decreasing function in r, one could obtain at high reshifts a central region already turnarounded and follows the Newtonian dynamics, while the outer region follows a non-Newtonian path heading towards a slow turnaround. The proposal of a power-law decreasing function was based on the observed galactic matter density profiles. The primordial black-holes or cosmological defects [82]-[83] proposed to account for the Supermassive Black Holes (SMBHs) at high redshifts (z 6) could also affect the density profile of the overdensity. ∼ In the case that the bulk of mass distribution is close enough to uniform density central Newtonian cloud, the Mestel [81] analysis could be an useful approximate solution. The free fall time tff from a z = 17 central break- away cloud with uniform density ρ 3.23 10−24kgm−3 is t 1.1Gyr, where the end of collapse redshift is z 7.7. ≥ × ff ∼ ∼ From the discussion of Early Type Spiral Galaxy observations in the introduction, nearly Newtonian central disks are seen to be accompanied by non-Newtonian rotating gas at large radii. A recent observation is reported, [84] that, a dense central protogalaxy at z 2.279 is clearly seen with its external cloud spaning upto a distance of 130kpc. ∼

4.2 Mass Shells inside ri

Outside rM , the non-Newtonian mass shells at rj inside ri, which will reach inside 50kpc, the angular speed from the Eq.(43) is

h2 4πGρ (z ) ϕ˙ 2 = j = δ + √δ m i (67) r4 3 j 13 whereϕ ˙ is the angular speed of the mass shell at ri. We assume that the mass shells at r < ri has the same angular speed and behave similarly as the largest mass shell at ri.

As a mass shell in the outer region will spend a number of redshifts to reach its turnaround and slowly rec- ollapses, it take a signiﬁcantly longer time than Newtonian free fall time to move towards the centre. This could account for the continuous infalling hydrogen gas to the Milky way over 8 Gyrs. Before the mass rings on the mid-plane starts to return from turnaround radius, the spherical shells will have enough time to collapse in the z direction onto a thin disk along the mid-plane where a ﬂat disk potential is more appropriate.

2 2 For constantϕ ˙ , h(rj ) =ϕr ˙ j is ﬁxed at the break-away. Mestel [81] suggests that if an uniform density sphere with uniform angular velocity collapses onto a thin spheroid without turbulent mixing, the distribution of angular momentum remains the same as that of the sphere. This idea is supported by ﬁndings in a group of disk galaxies in Crampin and Hoyle (1964). Following this idea, the observed h(rj ) on the disk remains that of the spherical cloud at rj .

For a ﬂat disk, ΦV is the VMOND gravitational potential, the acceleration of a ring is given by r¨ 1 ∂Φ (r) =ϕ ˙ 2 + V . (68) r r ∂r z=0

Using the thin disk potential in Mestel [81], at some r = rj , we have

h2 r¨j j 2 √ 4πG 2 Gρ(rj ) 4πG 8πG = 4 π Gαρ(rj ) = (δ + δ) ρm(zi) π Gρ(rj ) H(z)π ρm(z)+ ρΛ (69) rj rj − 3 − − r 2 − 3 3

2 3M where π Gρ(ri) is the disk Newtonian gravitational acceleration at rj and ρ(rj ) = δρm = 3 . The last part of 4πrj Eq.(69) comes from writing out the acceleration terms of the VMOND acceleration.

As the central Newtonian region collapses early, the outer mass shell will remain unmixed until equilibrium. From Eq.(69), we see that for Newtonian gravity (α = 1) [81], the initial angular acceleration of the sphere at rj is less than the disk gravitational acceleration, the acceleration r¨j < 0 is directed towards the mass centre. As r shrinks, rj 2 −4 −3 ϕ˙ r grows faster than ρ(rj ) r and the ring shell will reach equilibrium at about 2/5rj. The dynamics is qualitatively∝ similar here except for∝ the additional accelerations due to VMOND and the background accelerations.

After the turnaround, the eﬀective acceleration of a simple mass shell remains attractive toward the mass centre until it reaches equilibrium where acceleration goes to zero or when there is signiﬁcant mixing of shells. As these outer rings of disk progress towards an equilibrium over time, we consider the rotational velocity at r = 35kpc at late time, which is given by 4πG V 4(r )= ρ (z )(r )(1 + √δ)2n5GM = aV M (r )GM . (70) 35 3 m i 35 35 0 35 35

The mass lying between Mr0 M50 is about 10% of the overall mass and M35 0.98M50. In selecting galactic mass M in a Tully-Fisher plot, M =∼ M M is used in practical approximation. ∼ 35 ∼ 50

Consider an overdensity freezed at δ = 0.381 (zi = 21.5) where ri = 93kpc (n = 1.86), the outer shell will have travelled a distance of 43kpc to reach 50kpc. We assume a nearly uniform shell radial speed for different shells on the disk. The inner shells at rj = 40kpc would collapse onto the stellar disk and attains a dynamical process of inflow-outflow equilibrium [39], [84]-[88]. The disk shell at r35 comes from a shell at rj = 78kpc so that the V M contraction factor is nj = 78/35 and a0 at r35 = 35kpc on the disk is given by

n 5 r aV M (r )= aV M (r ) j 35 =1.20 10−10ms−2. (71) 0 35 0 50 n r × 50

This is also within the error bar of the phenomenological a0 in the literature. In a large central mass spiral galaxy such as our Milky Way, the Tully Fisher relation plot V 4 M is obtained where rotational speed V is taken as the V M ∝ average speed over the ﬂat speed range [80], so that a0 (r35) will be closer to the average observed value chosen than V M a0 (r50). 14

DISCUSSION AND CONCLUSION

With no modification to General Relativity, a metric (inhomogeneous cosmological solution) is found to possess a slow speed VMOND acceleration which includes both the Newtonian acceleration and a MOND-like acceleration where the phenomenologically a0 is replaced by the cosmological acceleration. At late time and at solar system scale, the VMOND acceleration is effectively Newtonian and satisfies a very stringent observation limit. In a previous work we find that in matter-radiation equality epoch, the non-Newtonian gravitational acceleration at high redshifts provides an external driving potential which can lead to higher odd peaks and could also dampen the baryonic acoustic oscillations after recombination. In the matter dominant epoch, the dominant non-Newtonian acceleration drives the overdensity to grow as δ a2 which provides the fast growth of overdensity as required without DM. The Early Type Galaxy (z 2) observations∝ that rotational curve is Newtonian at small scale length radii but obeys BTFR at large radii indicate∼ that it is the rotating gas at large radii which remains under the influence of the non-Newtonian acceleration over cosmic time.

Galactic overdensity growing under Newtonian gravity (with or without DM) is expected to reach values O(1) where the angular speed remains subdominant. In this study of a Milky Way scale protogalaxy, we follow a stronger assumption in the literature that over time an uniform angular acceleration catches up to the overdensity’s gravitational acceleration leading to the end of overdensity growth where overdensity value freezes. The angular speed is thus directly connected to the overdensity freeze value. At high redshifts, the overdensity freezes at a value much smaller than unity, where the MOND-like acceleration dominates. For the largest mass shell, BTFR holds automatically with a VMOND acceleration a0(zi, ri) deﬁned by the cosmological acceleration at overdensity freeze. As this mass shell’s radius contracts to galactic disk size, the VMOND acceleration a0(zi, r) grows towards the phenomenological MOND value a0 due to the conservation of angular momentum. From a realistic initial δ value at recombination and an observationally acceptable range of redshifts, one can predict the a0 value for the largest mass shell of an overdensity at late time. The predicted a0 matches the phenomenological MOND acceleration value a0 within its observable bounds. This is the main result of this work that the phenomenological MOND acceleration a0 scale has a natural and cosmological origin.

Observations requires that a central region of an overdense cloud collapses under Newtonian acceleration, which can be achieved if the overdensity has a power law (or simply a decreasing function behaviour) dependence on its radius. Mass shells in the outer layer having similar angular speed as the largest mass shell follow the same dynamics. These spherical shells collapse in the z direction without turbulent mixing in a thin disk potential. We follow the Mestel assumption that a speciﬁc angular momentum of a mass shell is invariant under contraction, the outer mass shells remain at a non-equilbrium state while moving slowly towards an equilibrium. We ﬁnd that after the outer disk contracts to within a distance of 50kpc, the VMOND acceleration upto 50kpc remains within the range of the phenomenological value of a0.

In Nusser [24], the overdensity under MOND can also grow as a(t)2 whch can lead to a Large Scale Struc- tures similar to a CDM universe, except that the MOND acceleration needs to be reduced by a ”signicant” factor (up to a factor of 10). In the VMOND scenario, the overdensity only grows as a(t)2 in the high redshift regime (z 12), which can lead to early structures appearance, but the growth would slow down to a(t) for z < 12. It is interesting≫ to see whether VMOND’s weaker growth rate is sufficient to mitigate the problem encountered by MOND. In galaxy formation, we notice that strong non-Newtonian gravitational acceleration only appears at high redshifts and at large distances. As redshift drops, the gravitational acceleration around a central mass eventually becomes Newtonian dominant, it is the angular momentum attained at high redshifts that remain non-Newtonian. For galaxies, instead of being an equilibrium system of gases where the force law is modified in the weak Newtonian MOND regime, in VMOND the non-Newtonian outer shells are relics of dynamics frozen from the high redshift epoch. For Elliptical galaxies, random angular momentum attained in non-Newtonian region by large mass shells could show up in the Faber-Jackson relations. Satellite galaxy that is formed from a small central mass at high redshifts will retain its high velocity dispersions. However, the Tidal Dwarf Galaxies which are formed by debris at late time naturally should have only Newtonian acceleration [89]. The high speed collision of 2 subclusters in El Gordo is accounted for by MOND [90] and could also be accounted for by VMOND if the collision starts at high redshifts. As the mass density around a central mass also includes an accompanied mass density which is cosmic time-varying, this could also lead to both strong and weak lensing effects. 15

References

[1] C. C. Wong, ”Variable Modified Newtonian Mechanics I: Single metric universe”, arXiv: 1601.00376. [2] S. S. McGaugh, APJ, 609, 652 (2004); ApJ. 632, 859 (2005); AJ, 143, 40 (2012). [3] F. Lelli, S. S. McGaugh, J. M. Schombert, M. S. Pawlowski, ”One Law to Rule Them All: The Radial Acceleration Relation of Galaxies” M. S, 2017, ApJ, 836, 152. [4] M Milgrom, R.H. Sanders,”MOND rotation curves of very low mass spiral galaxies”, arXiv:0611494, AJ Lett. 658, L17 (2007). [5] R. H. Sanders, ”The prediction of rotation curves in gas-dominated dwarf galaxies with modified dynamics”, arXiv:1811.0526. [6] Lasha Berezhiani, J. Khoury, ” Theory of Dark Matter Superfluidity”, Phys. Rev. D 92, 103510 (2015); arXiv:1507.01019. [7] Van den Bosch, Burket and Swaters, MBRAS, 326, 1205-1215 (2001). [8] O. Müller, et. al. ”A whirling plane of satellite galaxies around Centaurus A challenges cold dark matter cosmology”, arXiv:1802.00081; Müller et al., Science 359, 534, 2018. [9] P. Kroupa, ” Galaxies as simple dynamical systems: observational data disfavor dark matter and stochastic star formation”, Canadian Journal of Physics, 2015, 93(2): 169-202. [10] J. D. Bowman, A. E. Rogers, R. A. Monslave, T. J. Mozden, & N. Mahesh, Nature, 555, 67 (2018);arXiv: 1810.05912. [11] S. McGaugh, ”Strong Hydrogen Absorption at Cosmic Dawn: the Signature of a Baryonic Universe”, arXiv:1803.02365; To appear in RNAAS. [12] P. D. Mannheim, ”Alternatives to Dark Matter and Dark Energy”, arXiv:astro-ph/0505266v2, Prog. Part. Nuvl. Phys. 56 (2006) 340-445. [13] P.D Mannheim, D. Kazanas, AJ, 324,635 (1989). [14] M. Milgrom, APJ, vol. 270, 371-383 (1983). [15] B. Famaey, S. McGaugh, ”Challenges for Lambda-CDM and MOND”, arXiv: 1301.0623, Proceedings of the Meeting of the International Association for Relativistic Dynamics, IARD 2012, Florence. [16] S. McGaugh, ”A tale of two paradigms: the mutual incommensurability of ΛCDM and MOND. arXiv:1404.7525, Canadian Journal of Physics 93,250 (2015). [17] M. Milgrom, ApJ, 287, 571 (1984). [18] S. McGaugh, ”Predictions and Outcomes for the Dynamics of Rotating Galaxies”, arXiv:2004.14402. [19] R.Durazo R, et. al. ”A Universal Velocity Dispersion Profile for Pressure supported systems:Evidence for MONDian Gravity across Seven Orders of Magnitude in Mass”, APJ, 837: 179 (2017), arXiv:1506.04099. [20] B. Famaey, S. McGaugh, Living Reviews in Relativity, 15, 2012, 10; arXiv:1112.3960. [21] C. Skordis, T. Zlosnik, ” A general class of gravitational theories as alternatives to dark matter where the speed of gravity always equals the speed of light” , arXiv:1905.09465; ”A new relativistic theory for Modified Newtonian Dynamics”, arXiv:2007.00082. [22] , S. Dodelson, ”The Real Problem with MOND”, arXiv:1112.1320. [23] K. Pardo, D. N. Spergel, ”What is the price of abandoning dark matter? Cosmological constraints on alternative gravity theories”, arXiv:2007.00555, accepted for publication in PRL. [24] A. Nusser, ”Modified Newtonian dynamics of large-scale structure”, arXiv:0109016, MNRAS, 331,909, (2002). [25] A. Aguirre, J. Schaye, E. Quataert, ApJ, 561, 550 (2001). [26] R. H. Sanders, MNRAS, 342, 901, 2003. [27] A. Nusser, E. Pointecouteau, ”Modeling the formation of galaxy clusters in MOND”, arXiv:0505095, MNRAS. 366, 969-976 (2006). [28] A. Hodson, H. Zhao, ”Galaxy Cluster A1689 in Modifed MOND, MOG and Emergent Gravity”, arXiv:1703.10219. [29] H. Zhao, B. Famaey, Phys. Rev. D, 86, 067301 (2012). [30] J. Khoury, Phys. Rev. D, 91, 024022 (2015). [31] E. P. Verlinde, ”Emergent Gravity and the Dark Universe”, SciPost Phys. 2, 016 (2017), arXiv:1611.02269. [32] J. W. Moffat, ”Acceleration in Modified Gravity (MOG) and the Mass-Discrepancy Baryonic Relation”, arXiv:1610.06909v2. ′ [33] M. Haslbauer, I. Banik, P. Kroupa, ”The KBC void and Hubble tension contradict CDM on a Gpc scale a Milgromian dynamics as a possible solution”, arXiv:2009.11292v1, Published in the Monthly Notices of the Royal Astronomical Society. [34] R. Genzel et al., Nature 543, 397-401 (2017). [35] P. Lang et al., arXiv: 1703.05491, submitted to Astrophysical Journal. [36] Simons et. al., ”An Epoch of Disk Assembly” arXiv:1705.03474 ”z ∼ 2. [37] Kassin et. al., ”The Epoch of disk settling z ∼ 1 to now” arXiv:1207.7072, submitted to APJ. [38] S. Sachdeva, K. Saha, ”Survival of the pure disk galaxies over the last 8 billion years, arXiv:1602.08942v1. [39] G. Pezzulli, F. Fraternali, ”Accretion, radial flow and abundance gradients in spiral galaxies”, arXiv:1510.04289. [40] M. Mollá et. al.,”The role of gas infall in the evolution of disc galaxies”, MNRAS,1-12, 2016. [41] R. Swaters et. al., ”The link between the Baryonic Mass Distrtibution and the Rotation Curve Shape”, arXiv: 1207.2729. 16

accepted for publication in MNRAS. [42] T. Martinsson et. al., ” The DiskMass Survey, X. Radio synthesis imaging of spiral galaxies” arXiv:1510.7666, A & A, 585, A99(2016). [43] G. Angus et.al., ” Mass Models of disk galaxies from the DiskMass Survey in MOND”, arXiv:1505.05522, accepted for publication in MNRAS. [44] M. Cappellari, Ann. Rev. Astron. Astrophys. 54:597-665, 2016; arXiv:1602.04267v1. [45] M. Den Heijer et. al., A&A, 581 (2015) A98; arXiv:1509.05236. [46] A. Einstein, E. G. Straus, ”The Influence of the Expansion of Space on the Gravitation Fields Surrounding the Invidual Stars.”, Reviews of Modern Physics, 17(2-3), 120124. (1945); A. Einstein, E. G. Straus, ”Corrections and Additional Remarks to our Paper: The Influence of the Expansion of Space on the Gravitation Fields Surrounding the Invidual Stars.”, Reviews of Modern Physics, 18(1), 148149 (1946). [47] G. C. McVittie, ”The Mass-particle in an Expanding Universe.”, Monthly Notices of the Royal Astronomical Society, 93(5), 325339 (1933). [48] A. Oeftiger, ”On the Effect of Global Cosmological Expansion on Local Dynamics”, http://www.oeftiger.net/global-cosmic-expansion.pdf. [49] W. Rindler, M. Ishak, ” The Contribution of the Cosmological Constant to the Relativistic Bending of Light Revisited”, arXiv: 0709.2948, Phys.Rev.D76:043006,2007. [50] D. Lebedev, K Lake, ”On the influence of the cosmological constant on trajectories of light and associated measurements in Schwarzschild de Sitter space”, arXiv:1308.4931. [51] V. Faraoni, M. Lapierre-Leonard,”Beyond lensing by the cosmological constant” arXiv:1608.03164. [52] Oliver Plattella,” On the effect of the cosmological expansion on the gravitational lensing by a point mass”, Universe, 2016, 2(4), 25; arXiv:1609.00270. [53] Ali and Bhattancharya, ” Light bending, static dark energy and related uniqueness of Schwarzschild-de Sitter spacetime”, arXiv:1705.08602v2, Phys. Rev. D 97 (2018) 024029. [54] D. Tanoglidis, V. Pavlidou, T. N. Tomaras, ”Testing Lambda CDM cosmology at turnaround: where to look for violations of the bound?”, arXiv:1412.6671, JCAP, 12, 060 (2015). [55] C. C. Wong, ”Variable modified Newtonian Mechanics III: Turnaround radius”, arXiv:1910.10477 . [56] R. H. Sanders, ”Forming galaxies with MOND”, arXiv:arXiv:0712.2576, accepted for publication in MNRAS. [57] D. Crampin, F. Hoyle, APJ, 74, 186, 1964. [58] K. Freeman, ”On the disks of spiral and SO galaxy” APJ, vol 160, June 1970. [59] C. Leibovitz’s, Phys. rev. D Vol. 4, No. 10 1969 (2949-2955), Appendix A. [60] G. A, Baker Jr., ”Effects on the structure of the universe of an accelerating expansion” ”General Relativity and Gravitation” June 2002,volume 34, issue 6, pp 767-791. [61] J. Daboul, M. Nieto, ”Exact, E=0, solutions for general power-law potentials. I Classical orbits”. arXiv:hep-th/9408057, Phys. Rev. E52, 4430 (1995). [62] R. D. Reasenberg et.al. Astrophys. J. 234, L219 (1979). [63] J. D. Anderson, J. K. Campbell, R. F. Jurgens, E. L. Lau, X X Newhall, M. A. Slade III, and E. M. Standish, Jr., in: Proceedings of the Sixth Marcel Grossmann Meeting on General Relativity, Part A, ed. H. Sato and T. Nakamura, (World Scientific, Singapore, 1992), p. 353. [64] R. Barkana, A. Loeb, ” In the Beginning: The First Sources of Light and the Reionization of the Universe”, Phys. Rept. 349, 125-238, 2001; arXiv: 0010468. [65] P. A. Oesch, G. Brammer, P. van Dokkum; et al., ”A Remarkably Luminous Galaxy at z=11.1 Measured with Hubble Space Telescope Grism Spectroscopy”. APJ. 819 (2) 129 (2016); arXiv:1603.00461. [66] N. Wittenburg, P. Kroupa, B. Famaey,”The formation of exponential disk galaxies in MOND”, arXiv:2002.10941, accepted for publication in ApJ 18.01.2020. [67] J. I. Davies, ”A Heavy Baryonic Galactic Disc”, arXiv:1204.4649. [68] T.Licquia, J. Newman, APJ 806, 96 2(015); arXiv:1407.1078v3. [69] T. Licquia, J. Newman, ApJ 831, 71 (2016); arXiv:1607.05281. [70] P. M. W. Kalberla, L. Dedes, ”Global properties of the HI distribution in the outer Milky Way”, arXiv:0804.4831. [71] P. M. W Kalberla et. al. ” Dark matter in the Milky Way, II. the HI gas distribution as a tracer of the gravitational potential ”, A&A 2006-6326 (2018); arXiv 0704.3925. [72] P. J. E. Peebles, APJ. 155, 393 (1969). [73] S. D. M. White, ”Angular Momentum Growth in Protogalaxies”, APJ. 296:38-41 (1984). [74] E. Casuso, J. E. Beckman, ”On the Origin of the Angular momentum of Galaxies: Cosmological tidal torques supplemented by the Coriolis force”, MNRAS, 1-9, 2015. [75] A. Burkert, E D’Onghia, ”Galaxy formation and the cosmological angular momentum problem”, arXiv:0409540. [76] Gunn and Gott. AJ,176;1-19,1972 August 15. [77] S. White, D. Zaritsky, AJ, 394:1-6, (1992). [78] A. Fillmore, P. Goldreich, AJ, 281:1-8, (1984). [79] A. Nusser, MNRAS, 325, 1397-1401(2001). [80] F. Lelli, et. al, ”The baryonic Tully-Fisher relation for different velocity definitions and implications for galaxy angular momentum” MNRAS Vol. 484, 32673278 (2019); arXiv:1901.05966. [81] L. Mestel,”On the Galactic law of rotation”, MNRAS, 126, 553, 1963. [82] S. F. Bramberger et.al.”Cosmic String Loops as the Seeds of Super-Massive Black Holes”, arXiv:1503.02317, 17

JCAP06(2015)007. [83] R. Brandenberger, H. Jiao, ”Cosmic Textures and Global Monopoles as seeds for Super-Massive Black Holes”, arXiv:1908.04585. [84] D. C. Martin et. al., ”A giant protogalactic disk linked to the cosmic web”, Nature, 524, 192-195 (2015). [85] E. Pitts, R. Taylor 1989, MNRAS, 240, 373. [86] E. Pitt, R.Taylor 1996, MNRAS, 280, 1101. [87] Mayor and Vigroux (1981). [88] Lacey and Fall (1985). [89] F. Lelli et. al., ”Tidal Dwarf Galaxies” Disc Formation at z = 0, Galaxies 2015, 2, 1-x; arXiv:1511.00689. [90] E. Asencio, I. Banik, P. Kroupa, ”A massive blow for CDM the high redshift, mass, and collision velocity of the interacting galaxy cluster El Gordo contradicts concordance cosmology”, arXiv:2012.03950, MNRAS 500, 5249.