Galaxies in the early Universe appear to have grown too big too fast assembling into massive monolithic objects more rapidly than anticipated in the hierarchical Lambda cold dark matter CDM structure formation paradigm. The available photometric data are consistent with there being a population of massive galaxies that form early z  10 and quench rapidly over a short 1 Gyr timescale consistent with the traditional picture for the evolution of giant elliptical galaxies. Similarly kinematic observations as a function of redshift show that massive spirals and their scaling relations were in place at early times. Explaining the early emergence of massive galaxies requires either an extremely efficient conversion of baryons into stars at z gt 10 or a more rapid assembly of baryons than anticipated in CDM. The latter possibility was explicitly predicted in advance by modified Newtonian dynamics MOND. We discuss some further predictions of MOND such as the early emergence of clusters of galaxies and early reionization.
['Galaxy evolution', 'Galaxy formation', 'Galaxy masses', 'High-redshift galaxies', 'Protogalaxies', 'High-redshift galaxy clusters', 'Cold dark matter', 'Modified Newtonian dynamics', '594', '595', '607', '734', '1298', '2007', '265', '1069', 'Astrophysics - Astrophysics of Galaxies', 'Astrophysics - Cosmology and Nongalactic Astrophysics']
Accelerated Structure Formation The Early Emergence of Massive Galaxies and Clusters of Galaxies
{'Accelerated Structure Formation: the Early Emergence of Massive Galaxies and Clusters of Galaxies': 'Stacy S. McGaugh, 1 James M. Schombert, 2 Federico Lelli, 3 and Jay Franck 4  \n1 Department of Astronomy, Case Western Reserve University, Cleveland, OH 44106, USA 2 Institute for Fundamental Science, University of Oregon, Eugene, OR 97403, USA 3 INAF - Arcetri Astrophysical Observatory, Largo Enrico Fermi 5, I-50125, Firenze, Italy 4 1108 Sherman St., Longmont, CO 80501, USA', 'ABSTRACT': 'Galaxies in the early universe appear to have grown too big too fast, assembling into massive, monolithic objects more rapidly than anticipated in the hierarchical ΛCDM structure formation paradigm. The available photometric data are consistent with there being a population of massive galaxies that form early ( z ≳ 10) and quench rapidly over a short ( ≲ 1 Gyr) timescale, consistent with the traditional picture for the evolution of giant elliptical galaxies. Similarly, kinematic observations as a function of redshift show that massive spirals and their scaling relations were in place at early times. Explaining the early emergence of massive galaxies requires either an extremely efficient conversion of baryons into stars at z > 10 or a more rapid assembly of baryons than anticipated in ΛCDM. The latter possibility was explicitly predicted in advance by MOND. We discuss some further predictions of MOND, such as the early emergence of clusters of galaxies and early reionization.  \nKeywords: Cold dark matter (265), Galaxy evolution (594), Galaxy formation (595), Galaxy masses (607), High-redshift galaxies (734), MOND (1069), Protogalaxies (1298)', '1. INTRODUCTION': 'The formation and evolution of galaxies has been a central concern of cosmology since Hubble (1929) demonstrated that spiral nebulae are external stellar systems of size comparable to the Milky Way. Ideas about the formation of galaxies have ranged from the monolithic collapse of giant gas clouds (Eggen et al. 1962) to assembly through the merger of numerous protogalactic fragments (Searle & Zinn 1978). With the launch of JWST, we now have the opportunity to directly observe the assembly of galaxies at early epochs, providing a direct test of these ideas.  \nIn Λ Cold Dark Matter (ΛCDM), galaxies form in dark matter halos that originate from primordial density fluctuations that start small and grow gradually (Schramm 1992; Peebles 1993). Massive halos, and the galaxies that they contain, assemble from the merger of smaller halos (Wechsler & Tinker 2018). This hierarchical formation of structure is often depicted as a merger tree (Somerville & Kolatt 1999). The objects that are giant galaxies today are the products of the assembly of many protogalactic fragments.  \nThere are two basic effects in play: (i) the assembly of mass and (ii) the emergence of an observable, luminous galaxy through the accretion of gas and its conversion  \ninto stars. The timeline of mass assembly (i) is well quantified by N-body simulations (De Lucia et al. 2006; Srisawat et al. 2013). The second step is highly uncertain, depending on many aspects of gas physics and star formation. However, the luminosity of an individual galaxy cannot outpace its assembly rate (Naab et al. 2009; van der Wel et al. 2009; Nipoti et al. 2009). Observations of the luminosities of galaxies thus test the predicted formation history.  \nWe describe galaxy formation models in section 2. Constraints on the evolution of high redshift galaxies predating the launch of JWST are discussed in section 3 and new insights from JWST data are examined in section 4. Complementary constraints from kinematic observations are discussed in section 5. Taken together, the data indicate that structure formed in the early universe at an accelerated pace relative to the predictions of ΛCDM(section 6). This result had been anticipated well in advance of the observations (Sanders 1998; McGaugh 2015) as discussed in section 7. Section 8 provides a succinct summary.  \nWe adopt a vanilla ΛCDM universe with Ω m = 0 . 3, Ω Λ = 0 . 7, and H 0 = 70 kms -1 Mpc -1 for cosmologydependent quantities.', '2. GALAXY FORMATION MODELS': 'It is important to have at least two distinct hypotheses to compare and contrast, so we consider both monolithic and hierarchical galaxy formation models. The passive evolution of a monolithic galaxy that forms early is motivated by traditional inferences about the evolutionary history of giant ellipticals (Thomas et al. 2005; Renzini 2006; Bregman et al. 2006; Schombert & Rakos 2009; Schombert 2016; Chiosi & Carraro 2002). In contrast, hierarchical galaxy formation is expected in ΛCDM (White & Frenk 1991), with testable predictions provided by both semi-analytic galaxy formation models (SAMs) and hydrodynamical simulations. Though something of a straw-man on its own, the monolithic case is useful as a proxy for the predictions (Sanders 1998) of MOND (Milgrom 1983).', '2.1. Hierarchical Λ CDM Models': "Hierarchical galaxy formation in ΛCDM is a combination of in situ star formation in the largest of many progenitors and ex situ growth via merging (Mo et al. 1998; Henriques et al. 2015; Rodriguez-Gomez et al. 2016; Angl'es-Alc'azar et al. 2017; Behroozi & Silk 2018). The largest galaxies are predicted to form latest, as they take the longest to assemble. That large galaxies are observed to contain the oldest stars (Bregman et al. 2006; Schombert 2016) may thus seem like a contradiction, but it may simply mean that many of the stars formed in protogalaxies at early times prior to them merging into the final giant galaxy (De Lucia et al. 2006; Cattaneo et al. 2008; van der Wel et al. 2009; Nipoti et al. 2009). We thus expect to see many small precursor galaxies at high redshift for every modern giant (Newman et al. 2012; Conselice et al. 2022).  \nThere are many papers in the literature discussing the evolution of galaxies, including theoretical works that attempt to predict how observable galaxies are associated with their parent dark matter halos. While important details vary from model to model (Knebe et al. 2015), the basic prediction of the hierarchical build up of mass is common to all as it is fundamental to the ΛCDM structure formation paradigm. This is well documented by many simulations, for example the IllustrisTNG suite of magnetohydrodynamical simulations (Springel et al. 2018; Pillepich et al. 2018; Marinacci et al. 2018; Naiman et al. 2018; Nelson et al. 2018).  \nFig. 1 illustrates hierarchical galaxy formation as realized in the high resolution TNG50 simulation (Pillepich et al. 2019; Nelson et al. 2019). The merger tree of one model galaxy (subhalo 6 of TNG50-1) is shown. Early, high redshift epochs are dominated by mergers, a feature that is common to all galaxies in the hierarchical  \nformation scenario. There is in general a great diversity of merger tree morphologies, with some continuing to experience significant mergers up to the present time. The illustrated example experiences its last merger and quenches relatively early ( z ∼ 1), but the build-up of stellar mass in its largest progenitor follows the typical evolutionary track fairly closely, being only a little ahead of the median (Rodriguez-Gomez et al. 2016) at z ≈ 1.  \nAfundamental aspect of hierarchical galaxy formation is that a massive galaxy at z = 0 is the sum of many parts. As we look to high redshift, we do not expect to see an early version of the modern galaxy, but rather its many precursor components. There is no single entity whose evolution we can trace. The closest thing to that is the largest progenitor, which takes time to assemble. For example, the typical L ∗ galaxy in the Illustris simulation (Rodriguez-Gomez et al. 2016) takes about half a Hubble time (until z ≈ 0 . 7) to assemble half of its stars (Fig.1). This is why the observation of massive ( M ∗ > 10 10 M ⊙ ) galaxies at z > 6 is surprising; such objects should be rare occurrences because there has not yet been sufficient time to assemble large objects from the many contributing protogalactic fragments (Fig. 1).  \nSimulated galaxy stellar mass functions are in tolerable agreement with the data for z ≤ 2 (Genel et al. 2014; Furlong et al. 2015). Above this redshift, they start to diverge from the data in the sense that there are more bright, high mass galaxies than expected (Franck & McGaugh 2017), with the discrepancy becoming clear at z > 4 (Steinhardt et al. 2016). From 6 < z < 10, there is a pronounced excess in the numbers of galaxies with M ∗ ≈ 10 10 M ⊙ over that predicted (McGaugh 2024). JWST observations extend this excess to even higher redshift (section 4). We will discuss possible paths to understand this in section 6 after exploring the evidence.", '2.2. Monolithic Models': "An important touchstone in galaxy evolution is the case of the passive evolution of a monolithic island universe (Eggen et al. 1962). The assumption implicit in this hypothesis is that practically all of the mass currently in a galaxy has always been part of it; it evolves as a closed box since a formation redshift z f . This provides a convenient picture, but is an unrealistic oversimplification. Galaxy-mass balls of gas do not magically appear in the early universe; mass must assemble from the initial condition of a nearly homogeneous early universe (Planck Collaboration et al. 2020). This picture nevertheless provides a useful null hypothesis, and a starting point for more realistic models that assemble mass rapidly if not instantaneously (Chiosi & Carraro 2002).  \nFigure 1. A merger tree for a model galaxy from the TNG50-1 simulation (Pillepich et al. 2019; Nelson et al. 2019, left panel) selected to have M ∗ ≈ 9 × 10 10 M ⊙ at z = 0; i.e., the stellar mass of a local L ∗ early type galaxy (Driver et al. 2022). Mass assembles hierarchically, starting from small halos at high redshift (bottom edge) with the largest progenitor traced along the left of edge of the merger tree. The size of the symbol is proportional to the halo mass and the color bar illustrates the specific star formation rate. The growth of stellar mass of the largest progenitor is shown in the right panel. This example (jagged line) is close to the median (dashed line) of comparable mass objects (Rodriguez-Gomez et al. 2016), and within the range of the scatter (the shaded band shows the 16th - 84th percentiles). A monolithic model that forms at z f = 10 and evolves with an exponentially declining star formation rate (eq. 1) with τ = 1 Gyr (purple line) is shown for comparison.  \n<!-- image -->  \nThe growth of the stellar mass of a monolithic galaxy is determined by its star formation history. A common prescription for the passive evolution of a predominantly old stellar population like that of a typical elliptical galaxy (Bregman et al. 2006; Rakos et al. 2008; Schombert & Rakos 2009; Schombert 2016) is an exponentially declining star formation history,  \nψ ( t ) = ψ 0 e -u . (1)  \nHere, ψ 0 is a star formation rate that sets the scale that leads to a final mass M f ∗ , and  \nu = t -t i τ (2)  \nwhere t i is the time after the Big Bang when star formation begins that we equate to the redshift of galaxy formation z f , and τ is the timescale over which star formation activity fades. A population can be said to be passively evolving if this timescale is much shorter than a Hubble time so that most of the stars form in the early universe and evolve passively thereafter. The stellar mass increases as  \nM ∗ ( t ) = M f ∗ (1 -e -u ) . (3)  \nWith this prescription for the build-up of stellar mass M ∗ ( t ) and the time-redshift relation of vanilla ΛCDM,  \nthe luminosity evolution L ( z ) can be calculated using the methods of stellar population synthesis (Renzini 2006).  \nAn exponential model with τ = 1 Gyr and a formation redshift z f = 10 ( t i = 464 Myr for vanilla ΛCDM) forms half of its stellar mass by z = 5 (Franck & McGaugh 2017). This provides a baseline against which to compare other models (Fig. 1). To relate light to mass, we adopt an empirical calibration of this model by equating the final mass M f ∗ to the characteristic stellar mass M ∗ ∗ = 9 × 10 10 M ⊙ of the Schechter fit to local early type galaxies (Driver et al. 2022) and by matching the corresponding luminosity L ∗ to the characteristic apparent magnitude m ∗ of cluster galaxies at z ≈ 1 (Mancone et al. 2010). This data-informed choice of the mass-tolight ratio is within the range expected for the stellar population models built by Franck & McGaugh (2017): evolutionary theory and data agree at the expected level.  \nTo improve on the exponential model without adding much complexity, we also consider a generalization that allows for a finite ramp-up of star formation before it quenches. Kelson et al. (2016) show that stochastic star formation leads to an average linear ramp up in the star formation rate ψ ∼ t as galaxies accrete gas. This will quickly build up stellar mass as M ∗ ∼ t 2 . This process is more monolithic than hierarchical, as it envisions in situ  \nFigure 2. The expansion (solid lines) and collapse (dotted lines) of spherical regions of various baryonic masses in MOND as a function of time (Sanders 1998). Globular clusters collapse first and form early. Large galaxies reach maximum expansion after ∼ 3 × 10 8 yr; rich clusters of galaxies do so after ∼ 3 × 10 9 yr. Laniakea (Tully et al. 2014) and Ho'oleilana (Tully et al. 2023) may be examples of vast structures near turnaround today. The redshift scale on the top axis assumes vanilla ΛCDM for reference, but the timeredshift relation may differ in MOND.  \n<!-- image -->  \nstar formation from gas accretion onto a single object. Indeed, star formation must quench rapidly in order not to overproduce stellar mass. Regardless of the precise mechanism by which quenching occurs (Peng et al. 2015; Kimmig et al. 2023), an obvious choice to model it is an exponential attenuation with a short e-folding time. Combining this with an initially linear ramp-up gives  \nψ ( u ) = ψ 0 ue -u (4)  \nThis is only a slight modification to the traditional exponential star formation history described above, providing it with a more realistic initial condition. In principle, we could consider separate timescales for the ramp up of star formation and for its quenching, and could also insert a time delay for the commencement of quenching. We eschew these details for now as unnecessary complications, as we seek only to quantify the approximate time scales relevant to explaining observations of massive galaxies at high redshift. We thus restrict ourselves to the two timescales t i and τ that are built into the definition of u (eq. 2).  \nIntegrating equation 4, the stellar mass grows as  \nM ∗ ( u ) = M f ∗ [1 -(1 + u ) e -u ] (5)  \nwhere M f ∗ is the final stellar mass. While no individual galaxy will have exactly this star formation history, eq.  \n5 provides a simple way to describe the formation and quenching timescales that characterize a population of galaxies as represented by the typical Schechter function L ∗ galaxy and its corresponding mass M ∗ ∗ .", '2.3. MOND Galaxy Formation Models': "MOND has been very successful as a theory of galaxy dynamics (Sanders & McGaugh 2002; Famaey & McGaugh 2012; Banik & Zhao 2022), but has no completely satisfactory cosmology (Famaey & McGaugh 2012; Wittenburg et al. 2023). Perhaps the most successful attempt to combine MOND and General Relativity to date is the Aether-Scalar-Tensor (AeST) theory of Skordis & Zglyph[suppress]lo'snik (2019). AeST has been shown to fit the power spectrum of both the cosmic microwave background and galaxies at low redshift (Skordis & Zglyph[suppress]lo'snik 2021), but many details of the theory remain to be explored. One possible shortcoming is an apparent tension between the parameters required to explain the data at small, intermediate, and large scales (Mistele et al. 2023).  \nWhile the deeper theory remains unknown, we know the universe is expanding, and that MOND often works to describe the dynamics of objects within it. In what follows, we consider what happens to a region within the expanding universe that is subject to the MOND force law (Sanders 1998, 2008).  \nMOND is inherently nonlinear. A growth factor of ∼ 10 5 from z = 1090 to z = 0 is achieved through nonlinear growth (Nusser 2002; Knebe & Gibson 2004; Llinares et al. 2008) rather than linear growth with cold dark matter. Calculating the nonlinear growth of structure in MOND is a nontrivial problem (Sanders 1998; Stachniewicz & Kutschera 2001; Sanders 2001; Nusser 2002; McGaugh 2004; Llinares et al. 2008; Feix 2016; Wittenburg et al. 2020) and considerable work remains to be done. Nevertheless, a common feature of these analyses is a period of rapid structure formation. Indeed, this seems unavoidable when considering a MONDian region within an expanding background.  \nSanders (1998) considered spherical regions in the MONDregime of low accelerations in an expanding Universe. After first showing that the usual early universe results (e.g., Big Bang Nucleosynthesis) are retained, Sanders (1998) showed that there is a characteristic length scale r c below which MOND dynamics should apply. The value of r c grows with time, with small regions entering the MOND regime first and larger ones later in a hierarchical sequence, albeit one greatly accelerated relative to ΛCDM. On scales larger than r c , the Universe remains homogeneous and isotropic, so the conventional Friedmann equation may continue to apply. On scales smaller than r c , the region becomes detached  \nfrom the Hubble expansion and is destined to recollapse under its own gravity. This process is inevitably inefficient: not all regions of the universe will be neat spheres that are cleanly in the MOND regime. Those which are must necessarily collapse rapidly, so there should be a portion of the galaxy population that forms early.  \nThe governing equation for an initially expanding spherical region in the MOND regime (Felten 1984; Sanders 1998) is  \n( ˙ r ) 2 = (˙ r i ) 2 -( 16 π 3 a 0 Gρ b r 3 0 ) 1 / 2 ln( r/r i ) , (6)  \nwhere r 0 is the comoving radius of the spherical region, r i and ˙ r i are the initial radius and expansion velocity thereof, ρ b is the baryonic mass density, and a 0 = 1 . 2 × 10 -10 ms -2 (Begeman et al. 1991; McGaugh et al. 2016). The introduction of the dimensional constant a 0 makes the problem scale-dependent (Felten 1984), so the initial conditions matter. Fortunately, r i and ˙ r i have an obvious interpretation: since this growth can only commence after radiation releases its grip on the baryons, the initial velocity is simply the cosmic expansion rate at that time while the initial radius specifies the mass of the object. The precise redshift when this occurs is sensitive to the cosmology (Sanders 1998), but happens early enough ( z ≳ 200) that the net result for galaxy mass objects is not particularly sensitive to the initial conditions: the start time is small compared to the subsequent evolution.  \nA region that evolves according to eq. 6 reaches a maximum radius and recollapses on a timescale comparable to that of the initial expansion (Fig. 2). In a pure MOND universe, such a region is destined to recollapse (Felten 1984) irrespective of its initial density (eq. 6). Small regions naturally collapse faster than large ones, so the process is inherently hierarchical, but the timescale is greatly accelerated relative to the linear case. The mass of the top-hat sets the timescale for decoupling from the Hubble flow (Sanders 1998) and recollapse (Stachniewicz & Kutschera 2001). Hydrodynamics is important at the scale of globular clusters, which briefly delays their collapse, but is less important at galaxy scales (Stachniewicz & Kutschera 2001).  \nGlobular cluster mass ( ∼ 10 5 M ⊙ ) objects collapse very quickly, in the first 100 Myr (Stachniewicz & Kutschera 2001). This is so fast that they will have ages that are practically indistinguishable from that of the universe itself (Ying et al. 2023). Massive ( ∼ 10 11 M ⊙ ) galaxies reach maximum expansion around 300 Myr (Sanders 1998), with collapse happening on a similar timescale (Stachniewicz & Kutschera 2001). Wittenburg et al. (2020) consider slightly lower mass objects, and  \nfound that the initial spheres collapse to form a thin, rotating disk after ∼ 500 Myr. This is the epoch of galaxy formation in MOND (see also Sanders 2001; McGaugh 2004; Sanders 2008; Llinares et al. 2008). Thus, even though structure formation in MOND remains intrinsically hierarchical, the timescales are so short that the monolithic models described in section 2.2 provide a reasonable first approximation.", '3. HIGH REDSHIFT GALAXIES BEFORE JWST': 'The discussion above outlines two distinct hypotheses: hierarchical and monolithic galaxy formation. Hierarchical galaxy formation is expected in ΛCDM while monolithic galaxy formation is a first approximation to the accelerated structure formation expected in MOND. We can test these hypotheses by observing the evolutionary development of galaxies over a wide range of redshifts. In this paper we focus on tracing the growth of galaxy stellar mass utilizing photometric observations over the range 0 < z < 15. We also consider the constraints on dynamical mass provided by kinematic observations over the past 11 Gyr. These complimentary lines of evidence provide a consistent picture in which a significant population of massive galaxies formed remarkably early.  \nAs we look to high redshift, the first objects we see are always the brightest beacons that have been lit at that time. We must therefore take care to consider how typical these objects are. To do so, we utilize Schechter (1976) function fits to quantify the characteristic stellar mass M ∗ ∗ of large number of galaxies at each redshift. Locally, giant ellipticals galaxies have M ∗ ∗ = 9 × 10 10 h -2 70 M ⊙ (Driver et al. 2022). We wish to know the evolution of the typical galaxy, M ∗ ∗ ( z ). Of course, we cannot see the evolution of a single galaxy over cosmic time, and must attempt to infer M ∗ ∗ ( z ) from snapshots at different redshifts. This is famously problematic, as we can never relate a particular galaxy to its progenitor at higher redshift (Bell et al. 2004). Nevertheless, we can compare the data to evolutionary tracks from models to see which might work and which do not.  \nThe first data from JWST has brought the formation epoch of galaxies and evolution of M ∗ ∗ ( z ) into sharp focus. While this is a story in progress, it is possible to place these data in the context of precursor work with deep fields observed by HST and Spitzer. Critically, we are not limited to photometric redshifts; there exist many spectroscopic redshifts for many sources (see, e.g., the compilations of Franck & McGaugh 2016a,b).', '3.1. Galaxies with Spectroscopic Redshifts': "Giant elliptical galaxies are routinely found in dense regions like clusters of galaxies (Dressler et al. 1997)  \nFigure 3. The redshift dependence of the Spitzer [4.5] apparent magnitude m ∗ of Schechter function fits to populations of galaxies in clusters and candidate protoclusters: each point represents all the galaxies in each cluster. Data from Mancone et al. (2010, black circles), Wylezalek et al. (2014, dark blue circles), and Franck & McGaugh (2017, blue circles) all have spectroscopic redshifts. The orange line is the prediction of the Munich galaxy formation model (Henriques et al. 2015) based on the Millennium simulations (Springel et al. 2005; Boylan-Kolchin et al. 2009). Open squares are mock observations of this model using the same algorithm that was applied to the data (Franck & McGaugh 2016a). The predicted characteristic magnitude is fainter than observed, diverging systematically for z > 2. The purple line is a model of a galaxy formed at z f = 10 with an exponential star formation history (eq. 1) with τ = 1 Gyr (Franck & McGaugh 2017) normalized at z ≈ 1. Galaxies like this apparently exist in the high redshift universe, before they were predicted to have assembled, and are common enough to dominate the Schechter fit for m ∗ .  \n<!-- image -->  \nthat typically have well-defined red sequences (Bell et al. 2004). This makes clusters and protoclusters convenient environments in which to find a sufficient number of galaxies at the same redshift to construct luminosity functions. This has been done to progressively higher redshifts by Mancone et al. (2010), Wylezalek et al. (2014), and Franck & McGaugh (2017). These studies are all informed by galaxies with spectroscopically observed redshifts, so there is no ambiguity about their cosmic distance as can happen with photometric redshifts. The structures identified as [proto]clusters are redshift spikes in the N ( z ) diagrams of redshift surveys (Franck & McGaugh 2016a,b), so they represent a population of galaxies at the same point in the history of the universe regardless of whether they are indeed a bound structure.  \nFig. 3 shows the dependence of the AB Spitzer [4.5] apparent magnitude m ∗ (Franck & McGaugh 2017) corresponding to the characteristic luminosity L ∗ obtained  \nfrom Schechter function fits to many dozens, and sometimes hundreds of galaxies in clusters and protocluster candidates. This is a characteristic quantity of the galaxy population, not just a few anecdotal examples. Franck & McGaugh (2017) found no significant difference between the characteristic luminosity of cluster luminosity functions and that of galaxies in the surrounding fields at z > 2, so there does not appear to be a strong environment bias at high redshift.  \nGalaxies become fainter with increasing redshift, as expected (Fig. 3). However, observed galaxies are brighter than anticipated by contemporaneous models, e.g., the Munich galaxy formation model (Henriques et al. 2015). The model behaves as expected: earlier galaxies are small protogalaxies, so their characteristic luminosity becomes progressively fainter with increasing redshift. This generic expectation of ΛCDM diverges progressively from the data at z > 2 (Fig. 3).  \nFigure 4. The characteristic stellar mass of the Schechter mass function M ∗ ∗ as a function of time with the corresponding redshift noted at top. The data from Fig. 3 (Mancone et al. 2010; Wylezalek et al. 2014; Franck & McGaugh 2017) are augmented with higher redshift data (Stefanon et al. 2021, light blue points). The purple line is the passively evolving monolithic model from Fig. 3 normalized to M ∗ ∗ = 9 × 10 10 M ⊙ for local elliptical galaxies (Driver et al. 2022, grey point). The dotted orange line shows the build-up of the most massive progenitor of a galaxy that reaches this mass by z = 0 in the Munich SAM (Henriques et al. 2015). This is indistinguishable from the result of Illustris (Rodriguez-Gomez et al. 2016, red dashed line). Stellar mass in the FIRE simulation (Angl'es-Alc'azar et al. 2017, green dot-dash line) grows faster, but only a little.  \n<!-- image -->  \nAs a check that the same quantity was being measured in both data and model, Franck (2018) made mock observations of lightcones from the Munich model (Henriques et al. 2015). The same algorithm was applied to the mock data that was used to identify protocluster candidates in the real data (Franck & McGaugh 2016a,b). The mock observations recover basically the same answer that is known directly from the model (squares in Fig. 3). If the real universe looked like the prediction of the Munich model, we could easily tell. While it is tempting to blame the details of star formation in this particular model, the primary problem is more fundamental. Model galaxies are faint because hierarchical assembly is incomplete at z > 3 (compare Figs. 1 and 3).  \nIn contrast, the data fall around the line representing a monolithic giant that formed at z f = 10 and followed an exponential star formation history (eq. 1). If massive galaxies form early and evolve passively, it would look like the characteristic magnitudes that are observed. In addition to capturing the general trend of the data at high redshift, the data are very well fit at low redshift ( z < 1 . 5). This represents passive stellar evolution over most of cosmic time ( ∼ 9 Gyr) after essentially all the stellar mass has been formed and early stochastic variations have had time to subside. These data look very much like the evolution of a massive monolith that was assembled into a single object already at high redshift, and not like ΛCDM models in which the largest progen-  \nr should have been much smaller and fainter at z > 2 (Figs. 1 and 3).  \nNo single model line will fit all the data in Fig. 3. The scatter at high redshift presumably reflects stochastic variations in star formation rates at early times (Kelson et al. 2016; Pallottini & Ferrara 2023) before the red sequence was established (Bell et al. 2004; Rakos et al. 2008; Schombert & Rakos 2009; Franck et al. 2015). Nevertheless, the evolutionary trend predicted by the monolithic model captures the essence of the data in a way that the nominal prediction of ΛCDM does not.", '3.2. Galaxies at z > 6 Before JWST': "The bulk of the data discussed above are for z < 4 with a couple of candidate protoclusters extending to z ≈ 6 (Franck & McGaugh 2016b). There have been many studies of galaxies to yet higher redshift that predate JWST (e.g., Finkelstein 2016; Grazian et al. 2015; Song et al. 2016; Stefanon et al. 2021; Santini et al. 2022; Weaver et al. 2023). These works provide M ∗ ∗ for field galaxies with stellar mass functions measured independently of any of the data described above, albeit at the cost of relying more, if not entirely, on photometric redshifts.  \nFig. 4 shows the characteristic stellar mass M ∗ ∗ as a function of redshift. The data from Fig. 3 are shown assuming an exponential star formation history to provide a mapping from magnitude to mass that preserves the distribution of the data. This compares well to the data at higher redshift for which the stellar mass  \nFig. 4 also shows the a priori predictions of several ΛCDM models. These include the Munich SAM (Henriques et al. 2015, as in Fig. 3) and the hydrodynamical simulations Illustris (Rodriguez-Gomez et al. 2016) and  \n<!-- image -->  \nFigure 5. The UV luminosity function (left) observed by Donnan et al. (2024, points) compared to that predicted for ΛCDM by (Yung et al. 2023, lines) as a function of redshift. Lines and points are color coded by redshift, with dark blue, light blue, green, orange, and red corresponding to z = 9, 10, 11, 12, and 14, respectively. There is a clear excess in the number density of galaxies that becomes more pronounced with redshift, ranging from a factor of ∼ 2 at z = 9 to an order of magnitude at z ≥ 11 (right).  \n<!-- image -->  \nestimates are entirely independent (Finkelstein 2016; Grazian et al. 2015; Song et al. 2016).  \nWe include in Fig. 4 the data of Stefanon et al. (2021) together with those from Fig. 3. This is the most conservative choice in the sense that the stellar mass estimates of Stefanon et al. (2021) are the lowest available at these redshifts. This happens in large part because Stefanon et al. (2021) make larger corrections for line emission from non-stellar sources. The data are all consistent; they simply attribute less of the observed luminosity to stars.  \nThe data are consistent with a population of massive galaxies that formed early and evolved passively. Despite its naive simplicity, the exponential star formation history (eq. 3 with τ = 1 Gyr) provides a remarkably reasonable depiction of the build up of the characteristic stellar mass seen in Fig. 4. This is highly non-trivial, as the mass build-up happens early while the luminosity evolution is most pronounced at late times (Fig. 3).  \nMore complex star formation histories are admissible. Indeed, stochastic star formation may well drive the scatter seen in the data around z ≈ 3 when the universe was only ∼ 2 Gyr old and stellar populations were necessarily still young. Regardless of the details of the early star formation history, it appears that there exists a population of massive galaxies that formed early and in which most of the stars were made long ago in a single object rather than the multiplicity of progenitors envisioned in hierarchical galaxy formation.  \nFIRE (Angl'es-Alc'azar et al. 2017). These all show basically the same thing. Galaxies are predicted to assemble gradually, with their most massive progenitor not reaching half the final stellar mass until half a Hubble time has passed ( z < 1). To give a specific example, the Munich SAM reaches half the final stellar mass at z = 0 . 68 when the universe is 7.3 Gyr old. The star formation prescription of FIRE makes more stars earlier, but it is only a small shift of the same basic result: mass assembles too slowly in ΛCDM models. Making star formation more efficient makes more stars earlier, but it does not assemble them into the massive individual galaxies that are observed.  \nThe shortfall of stellar mass in individual galaxies is especially severe at high redshift. At z ≈ 3, the largest progenitor of an L ∗ galaxy is only a tenth of its eventual z = 0 mass in FIRE. It is even less in Illustris, about 3%. At z ≈ 5, individual progenitor galaxies are not expected to have grown massive enough to even appear on Fig. 4, let alone to do so at z ≈ 10. Yet at these redshifts the data show that many galaxies with masses that already approach M ∗ ∗ ( z = 0) - enough galaxies to define L ∗ in a Schechter fit for every point in Fig. 4. This is the normal galaxy population, not just a few extreme individuals. Moreover, spectroscopic redshifts are required to be a part of the samples (Franck & McGaugh 2016a,b) that inform the Schechter function fits of Franck & McGaugh (2017), so there is no uncertainty due to photometric redshifts.", '4. HIGH REDSHIFT GALAXIES WITH JWST': 'JWST has made the observation of galaxies at z > 10 seem mundane, so it is worth recalling that this is a  \nrecent development. Franck & McGaugh (2017) built the model with z f = 10 shown in Figures 3 and 4 as an extreme upper limit in the context of the widespread presumption at the time that there was no possibility for massive galaxies to have formed that early. Consequently, early JWST results came as a surprise (e.g., Finkelstein et al. 2022a; Merlin et al. 2022; Melia 2023; Ferrara et al. 2023). However, they merely extend the trends already seen in earlier data, corroborating previous indications that galaxies grew too big too fast (Rocca-Volmerange et al. 2004; Steinhardt et al. 2016; Franck & McGaugh 2017; Merlin et al. 2019). The simple observation is that the high redshift universe contains a bounty of bright, morphologically mature galaxies (Ferreira et al. 2022, 2023) that are more luminous than had been anticipated by ΛCDM models (e.g., Yung et al. 2019a, 2022; Behroozi et al. 2020).', '4.1. The UV Luminosity Function': "There is a clear excess in the number density of M ∗ ≈ 10 10 M ⊙ galaxies over the predictions of contemporaneous ΛCDM models (Yung et al. 2019a,b) at z ≈ 8 (McGaugh 2024). This becomes more challenging to assess at z > 10 where much of the observed luminosity is in the ultraviolet where the short-lived nature of high mass stars makes it is difficult to assess the corresponding stellar mass. Bearing this caveat in mind, we can nevertheless compare (Fig. 5) observations (Donnan et al. 2024) with the updated predictions of Yung et al. (2023). There is a clear excess that is apparent at all luminosities. Similar results follow from comparison to other predictions (e.g., Fig. 1; Behroozi et al. 2020) and other analyses of the observations (Robertson et al. 2023). ΛCDM models did not anticipate the large number of relatively bright galaxies that are observed.  \nThe excess in the number density of UV-bright galaxies is not subtle, being an order of magnitude at z ≥ 11. This is true despite an upward adjustment of the density in the models by 0.3 to 0.4 dex from Yung et al. (2019a) to Yung et al. (2023). One can imagine a number of ways to further enhance the UV luminosity per unit mass (Finkelstein et al. 2023), but the salient observational fact is that the UV luminosity function barely evolves over the redshift range that the dark matter halo mass function is evolving rapidly. Consequently, any appeal to the efficiency of UV light production (or attenuation: Ferrara et al. 2023) must necessarily be fine-tuned to balance the barely-evolving UV luminosity function with the rapidly evolving dark matter halo mass function over a rather small window of cosmic time, there being only ∼ 10 8 years between z = 14 and 11.  \nFigure 6. Mass estimates for high redshift galaxies from JWST. Colored points based on photometric redshifts are from Adams et al. (2023a, dark blue triangles), Atek et al. (2023, green circles), Labb'e et al. (2023, open squares), Naidu et al. (2022, open star), Harikane et al. (2023, yellow diamonds), Casey et al. (2023, light blue left-pointing triangles), and Robertson et al. (2023, orange right-pointing triangles). Black points from Wang et al. (2023, squares), Carniani et al. (2024, triangles), Harikane et al. (2024, circles) and Castellano et al. (2024, star) have spectroscopic redshifts. The upper limit for the most massive galaxy in TNG100 (Springel et al. 2018) as assessed by Keller et al. (2023) is shown by the light blue line. This is consistent with the maximum stellar mass expected from the stellar mass-halo mass relation of Behroozi et al. (2020, solid blue line). These merge smoothly into the trend predicted by Yung et al. (2019b) for galaxies with a space density of 10 -5 dex -1 Mpc -3 (dashed blue line), though Yung et al. (2023) have revised this upwards by ∼ 0 . 4 dex (dotted blue line). This closely follows the most massive objects in TNG300 (Pillepich et al. 2018, red line). The light grey region represents the parameter space in which galaxies were not expected in ΛCDM. The dark grey area is excluded by the limit on the available baryon mass (Behroozi & Silk 2018; Boylan-Kolchin 2023).  \n<!-- image -->  \nIrrespective of the ultimate interpretation, Fig. 5 clearly illustrates the excess in the number of bright high redshift galaxies over that predicted a priori .", '4.2. The Mass and Evolution of Individual Galaxies': "Despite the clear excess in bright galaxies seen in Fig. 5, it is not yet possible to define a robust populationwide estimate of the corresponding stellar mass M ∗ ∗ as seen in Fig. 4. Nevertheless, there are many individual galaxies (Fig. 6) that have a reasonably persuasive claim to being high mass at high redshift (Adams et al. 2023a; Atek et al. 2023; Labb'e et al. 2023; Naidu et al. 2022; Harikane et al. 2023; Finkelstein et al. 2022b; Robertson et al. 2023). These are based on photometric red-  \nFigure 7. The stellar masses of quiescent galaxies from Nanayakkara et al. (2024). The inferred growth of stellar mass is shown for several cases, marking the time when half the stars were present (small green circles) to the quenching time (mid-size orange circles) to the epoch of observation (large red circles). Illustrative star formation histories following equation 5 are shown as dotted lines with parameters t i , τ in Gyr noted. We omit the remaining lines for clarity, as many cross. There is a wide distribution of formation times from very early ( t i = 0 . 2 Gyr) to relatively late ( > 1 Gyr), but all of the galaxies in this sample are inferred to build their stellar mass rapidly and quench early ( τ < 0 . 5 Gyr).  \n<!-- image -->  \ns, but similar examples with spectroscopic data are becoming available (Wang et al. 2023; Carniani et al. 2024; Harikane et al. 2024; Castellano et al. 2024; Price et al. 2024). The essential outstanding feature that they share is how bright they are: galaxies simply do not appear to be as faint as anticipated by ΛCDM models (Fig. 3). There are of course some interlopers (Arrabal Haro et al. 2023), but changes in the photometric redshift (as sometimes happens: Adams et al. 2023a) do not necessarily reduce the stellar mass, as a reduction in distance is accompanied by a shift of the flux to bands where the stellar mass-to-light ratio is larger than in the UV. The data discussed by Labb'e et al. (2023) are a case in point: all of the initial redshifts and masses were revised downwards after accounting for the on-sky calibration of JWST, but all of them remain problematic in terms of their mass for their redshift.  \nFigure 6 shows the stellar masses and redshifts of high redshift galaxies identified in JWST data. These are individual galaxies rather than the results of Schechter fits to many as in Fig. 4. Consequently, they are more plausibly subject to the concern of being rare instances of  \nextreme outliers. To assess the significance of these extreme values, the lines in Fig. 6 illustrate the maximum stellar mass found in various studies (Keller et al. 2023; Behroozi et al. 2020; Yung et al. 2019b, 2023).  \nTo make our own assessment, we queried available simulations to find the most massive model galaxy as a function of redshift. The most generous estimate emerges from the Illustris TNG300 (Pillepich et al. 2018) simulation (the red line Fig. 6), which closely tracks the estimate of Yung et al. (2023). These are the most massive objects in a large simulation volume, so observing them would require that JWST happened to observe regions of exceptional overdensity (McCaffrey et al. 2023; Kragh Jespersen et al. 2024). We therefore consider this line to be a conservative upper limit on stellar mass: objects near this line should be exceedingly rare, and nonexistent beyond it.  \nThere are a number of galaxies observed to exceed the limits illustrated in Fig. 6, including cases with spectroscopic redshifts (Carniani et al. 2024; Harikane et al. 2024). A few objects approach the limit where all the available baryons would need to already be formed into  \nstars (Boylan-Kolchin 2023). These galaxies are difficult to understand in ΛCDM, but they are consistent with a continuation in the trend already seen at intermediate redshift, so are less surprising from an empirical perspective.", '4.3. Quenched Galaxies': 'Another important observation is that of quenched galaxies at 3 < z < 4 (Schreiber et al. 2018; Merlin et al. 2019; Nanayakkara et al. 2024; Glazebrook et al. 2023). These galaxies have observed spectra that show the classic features of a stellar population aging after intense star formation at an earlier epoch. Not only do massive galaxies exist by z ≈ 4, but there are examples that have stellar populations that are old for the age of the universe at the redshift of observation. By modeling the observed spectra, it is possible to estimate the stellar mass at the time of quenching and roughly when when half the stellar mass was in place in addition to the mass at the observed redshift. This provides an approximate curve of growth for each galaxy, as illustrated in Fig. 7.  \nThe growth rates inferred by Nanayakkara et al. (2024) are consistent with the rapid rate of growth illustrated by the monolithic model (eq. 5). Individual galaxies vary in mass, but all are consistent with following a similar evolutionary trajectory. Further examples of such galaxies at yet higher redshift are discussed by Wang et al. (2024a). These galaxies grew too big too fast, well ahead of the expectation in ΛCDM models.  \nGalaxies that formed most of their stars early and then quenched are consistent with the traditional view of the evolution of monolithic early type galaxies. That the galaxies observed by Nanayakkara et al. (2024) have masses consistent with other high redshift galaxies and show the evolved spectra expected for stellar populations descended from earlier star formation suggests that the most obvious interpretation of the data is also the most likely: bright galaxies at high redshift are intrinsically luminous because they contain lots of stars. They appear to have formed as giant monoliths at early times.  \nThe depiction of evolutionary tracks in Fig. 7 implicitly assumes all the mass was assembled at an early time. It is also conceivable that the galaxies observed by Nanayakkara et al. (2024) at z ≈ 3 to 4 were not individual objects at the time of quenching or when half the stellar mass had formed. These events could instead have occurred in protogalactic fragments that subsequently merged to form the observed galaxies. Since the assembled stars are old for the epoch of observation, the assembly must occur as dry mergers devoid of star formation (Newman et al. 2012; Conselice et al. 2022).  \nFigure 8. The stellar masses and star formation rates of galaxies from Nanayakkara et al. (2024, red symbols). Downward pointing triangles are upper limits; some of these fall well below the edge of the plot so are illustrated as the line of points along the bottom. Also shown are objects selected from the TNG50 (Pillepich et al. 2019, solid squares) and TNG300 (Pillepich et al. 2018, open squares) simulations at z = 3 to cover the same range of stellar mass. Simulated objects with stellar masses comparable to real galaxies are mostly forming stars at a rapid pace. In the higher resolution TNG50, none have quenched as observed.  \n<!-- image -->  \nTo check what ΛCDM predicts, we have searched the TNG50 and TNG300 simulations (Nelson et al. 2018, 2019) for model galaxies at z = 3 with stellar masses in the same range as the data of Nanayakkara et al. (2024). Many examples of such objects exist in the simulations, but almost all are actively star forming (Fig. 8). Quenched galaxies are rare in the simulations at this redshift, and non-existent in the higher resolution TNG50, with all branches of the merger trees experiencing high specific star formation rates at z > 3 (Fig. 1). So while it is possible to find simulated objects of the observed stellar mass, their star formation histories are not a good match to those observed.', '5. GALAXY KINEMATICS': "The discussion heretofore has focussed on the photometric evidence. Kinematic observations provide an independent line of evidence that mature galaxies appeared early in the history of the Universe. Disk galaxies at intermediate redshift (1 < z < 3) are observed to have large rotation speeds (Neeleman et al. 2020), to be dynamically cold (Di Teodoro et al. 2016; Lelli et al. 2018, 2023; Rizzo et al. 2023), and to follow scaling relations like Tully-Fisher (Miller et al. 2012; Pelliccia et al. 2017). The Tully-Fisher relation persists up to at least z ≈ 2 . 5 when the universe was ∼ 2 . 5 Gyr  \nFig. 9 shows two scaling relations at both low and high redshift: the baryonic mass-circular speed relation (Tully & Fisher 1977; McGaugh et al. 2000) and the dark matter fraction-surface brightness relation (de Blok & McGaugh 1997; Starkman et al. 2018). Both relations are clearly present in the data of Nestor Shachar et al. (2023). There is little if any indication of evolution in either relation up to z ≈ 2 . 5. The good agreement between low and high redshift samples is remarkable given that we have made no attempt to reconcile the choice of circular velocity measure (Lelli et al. 2019) or the precise definition of baryonic mass, which is sensitive to the stellar population model and its evolution (Schombert et al. 2019). Any systematic differences between studies are apparently within the scatter induced by measurement uncertainties.  \n<!-- image -->  \nFigure 9. The baryonic Tully-Fisher (left) and dark matter fraction-surface brightness (right) relations. Local galaxy data (circles) are from Lelli et al. (2019) for the Tully-Fisher relation and for spiral disks (Hubble stage T < 9) from SPARC Lelli et al. (2016) for the dark matter fraction-surface brightness relation. Higher redshift data (squares) are from Nestor Shachar et al. (2023) in bins with equal numbers of galaxies color coded by redshift: 0 . 6 < z < 1 . 22 (blue), 1 . 22 < z < 2 . 14 (green), and 2 . 14 < z < 2 . 53 (red). Open squares with error bars illustrate the typical uncertainties. The relations known at low redshift also appear at higher redshift with no clear indication of evolution over a lookback time up to 11 Gyr.  \n<!-- image -->  \nold (Nestor Shachar et al. 2023). Individual galaxies with high circular speeds and relatively high rotation-tovelocity dispersion ratios are found up to z ≈ 5 (Rizzo et al. 2020, 2021; Lelli et al. 2021; Roman-Oliveira et al. 2023), barely one billion years after the Big Bang.  \nThe early appearance of massive, dynamically cold disks in the first few billion years after the Big Bang is contradictory to early ΛCDM predictions. For example, Mo et al. (1998) anticipated that 'present-day discs were assembled recently (at z ≤ 1).' Early disks are expected to be small and dynamically hot (Dekel & Burkert 2014; Zolotov et al. 2015; Krumholz et al. 2018; Pillepich et al. 2019). Kinematic scaling relations like Tully-Fisher are expected to emerge late and evolve significantly (e.g., Glowacki et al. 2021).  \nThe high rotation speeds observed in early disk galaxies are remarkable. These sometimes exceed 250 (Neeleman et al. 2020) or even 300 km s -1 (Nestor Shachar et al. 2023; Wang et al. 2024), comparable to the most massive local spirals (Noordermeer et al. 2007; Di Teodoro et al. 2021, 2023). Kinematics indicate large dynamical masses for these early galaxies. The problem is not limited to luminosity; the underlying dynamical mass is also larger than expected.  \nThe study of kinematics at still higher redshift is a nascent field, but there are already important individual cases. For example, the kinematics of ALESS 073.1 at z ≈ 5 indicate the presence of a massive stellar bulge as well as a rapidly rotating disk (Lelli et al. 2021). A similar case has been observed at z ≈ 6 (Tripodi et al. 2023). These kinematic observations indicate the presence of mature, massive disk galaxies well before they  \nwere expected to be in place (Pillepich et al. 2019; Wardlow 2021). Spiral galaxies are ubiquitous in JWST images up to z ∼ 6 (Ferreira et al. 2022, 2023; Kuhn et al. 2023).  \nThe dark matter fraction f DM = 1 -( V b /V c ) 2 , where V c is the observed circular velocity and V b is that due to the baryons at the same radius. Many high surface brightness galaxies are maximal in their inner regions, so V b → V c so f DM → 0 (Starkman et al. 2018). The precise value of f DM is very sensitive to the stellar population model, which determines the amplitude of V b . This is less critical for low surface brightness galaxies, which have long been known to be dark matter dominated (de Blok & McGaugh 1997). Fig. 9 shows that the relation between the dark matter fraction and surface brightness was already in place at intermediate redshift  \n(Nestor Shachar et al. 2023) and has not evolved much over cosmic time (Sharma et al. 2023).  \nKinematic observations to date show that dynamically cold, massive disks are already present in the universe at early times. These disk galaxies appear to follow the same kinematic scaling relations that are known locally. The presence of dynamically massive galaxies in settled kinematic scaling relations at early times does not follow naturally from the hierarchical galaxy formation picture.", '6. DISCUSSION: ΛCDM': "ΛCDMgalaxy formation models (e.g., Henriques et al. 2015; Rodriguez-Gomez et al. 2016; Angl'es-Alc'azar et al. 2017; Yung et al. 2019b; Behroozi et al. 2020) failed to anticipate the bright galaxies that are observed at both intermediate and high redshift. Forming enough stars is not the problem. The problem is assembling them into a single object. The presence of these quasimonolithic objects appears to violate the hierarchical assembly paradigm.", '6.1. Observed Numbers and the Mass Function': 'There are two basic issues, one theoretical and one observational. On the theory side, we must be careful about what the theory predicts, including the inherent uncertainties in doing so. The mass function of the dark matter halos is well predicted, but relating this to stellar mass requires a prescription for the star formation efficiency. Relating the mass of stars formed to the observed magnitudes requires a stellar population model. These are good, but never perfect. Stellar population models also play a role on the observational side in fitting the observed spectral energy distribution to obtain a simultaneous estimate of the stellar mass and redshift. This process can go amiss (e.g., Adams et al. 2023b), but enough spectroscopic observations now exist (Wang et al. 2023; Carniani et al. 2024; Harikane et al. 2024; Castellano et al. 2024; Price et al. 2024) that it seems unlikely that the whole problem can be one of misleading observations.  \nThe normalization of numerical simulations is not trivial. Yung et al. (2023) reassessed their predictions (Yung et al. 2019b) in the light of higher resolution simulations. This resulted in an upward shift in numbers at a given redshift by a factor of ∼ 2 (Fig. 6). This helps, but only a bit. There remain galaxies of larger mass at higher redshift than should be possible, including examples with spectroscopic redshifts (Wang et al. 2023; Carniani et al. 2024).  \nIn contrast, Katz et al. (2023) find a good match of the Sphinx simulation (Rosdahl et al. 2018, 2022) with the cumulative number counts at 9 < z < 13. In this case,  \nthe test may come at lower redshift to which it is not yet computationally feasibly to extend Sphinx. Fig. 34 of Katz et al. (2023) shows the surface density n ( z ) increasing rapidly with decreasing redshift, implying that this simulation may fit the high redshift data at the risk of overshooting the data at lower redshift (Merlin et al. 2019; Sabti et al. 2024). The evolutionary trajectories that are common in simulations do not have the right shape to explain the data (Fig. 10).', '6.2. Evolutionary Trajectories': "The expectation of hierarchical galaxy formation in ΛCDM is illustrated in Fig. 1. The stellar mass of a galaxy is a combination of in situ star formation in its largest progenitor and ex situ star formation in a multitude of protogalactic clumps that ultimately merge into it. This is a gradual process, with the median stellar mass reaching half its final value at z < 1 (De Lucia et al. 2006; De Lucia & Blaizot 2007; Rodriguez-Gomez et al. 2016; Angl'es-Alc'azar et al. 2017). The brightness of the typical galaxy is thus expected to diminish rapidly as we look to high redshift (Fig. 3), both because the largest progenitor is smaller and less mature, and because it should split into many precursor protogalaxies. Here the simple monolithic galaxy model provides a useful contrast in which all the baryons assemble promptly at an early time and the luminosity evolution is due entirely to in situ star formation.  \nFig. 10 combines the data from Figures 4 and 6 to show the stellar mass of galaxies observed over all of accessible cosmic time. These data are further augmented with recent spectroscopic observations by Xiao et al. (2023) and Carnall et al. (2024). The data begin to fill out the stellar mass-redshift diagram, with only the top right portion of Fig. 10 being empty. The remainder of the diagram is populated by galaxies of a wide range of masses and evolutionary states.  \nThe hierarchical ΛCDM model is represented in Fig. 10 by three evolutionary tracks from the Illustris simulation (Vogelsberger et al. 2014). These illustrate the median stellar mass growth of the largest progenitors that reach M f ∗ = 10 10 , 10 11 , and 10 12 M ⊙ at z = 0 (Rodriguez-Gomez et al. 2016). These trajectories parallel one another with the primary difference being that in normalization. These are typical objects; atypical objects also follow a similar trajectory with a higher or lower normalization (e.g., the maximum stellar mass line in Fig. 6). Other ΛCDM simulations (e.g., Angl'esAlc'azar et al. 2017; Katz et al. 2023) follow a similar trajectories with modest variations (Fig. 4). These variations stem from differences in the implementation of  \nFigure 10. The data from Figures 4 and 6 shown together using the same symbols. Additional JWST data with spectroscopic redshifts are shown from Xiao et al. (2023, green triangles) and Carnall et al. (2024). The data of Carnall et al. (2024) distinguish between star forming galaxies (small blue circles) and quiescent galaxies (red squares); the latter are in good agreement with the typical M ∗ ∗ determined from Schechter fits in clusters (large circles). The dashed red lines show the median growth predicted by the Illustris ΛCDM simulation (Rodriguez-Gomez et al. 2016) for model galaxies that reach final stellar masses of M ∗ = 10 10 , 10 11 , and 10 12 M ⊙ . The solid lines show monolithic models with M f ∗ = 9 × 10 10 M ⊙ and t i = τ = 0 . 3, 0.4, and 0.5 Gyr (equation 5) as might be appropriate for giant elliptical galaxies. The dotted line shows a model appropriate to a spiral galaxy with t i = 0 . 5 and τ = 13 . 5 Gyr.  \n<!-- image -->  \nbaryonic physics, not in the more fundamental assembly of mass.  \nMonolithic models that form early and reach a mass of M f ∗ = 9 × 10 10 M ⊙ at z = 0 are shown with three lines that follow eq. 5 with t i = τ = 0 . 3, 0.4, and 0.5 Gyr. These represent the traditional picture of a giant elliptical galaxy that is in place early, forms its stars in an initial burst, and evolves passively thereafter. A giant spiral model is also illustrated with t i = 0 . 5 and τ = 13 . 5 Gyr to show the effect of a more extended period of star formation.  \nThe evolutionary tracks of the monolithic models naturally explain the high redshift ( z > 6) region of the M ∗ -z plane. This region is not explored by the ΛCDM tracks, and is not accessible to typical ΛCDM models; we must appeal to rare outliers. The reason for this not merely a difference in star formation history; it is a consequence of the time required to hierarchically assemble mass (Fig. 1). Observed galaxies appear to have assembled more promptly than anticipated by ΛCDM.  \nA deviation from the predicted hierarchical assembly of mass is a considerably greater problem for ΛCDM than the details of star formation in the largest progenitor.  \nAnother problem is the shape of the simulated ΛCDM evolutionary tracks. The tracks are nearly linear in Fig. 10. This is fine for assembling a spiral galaxy with an extended star formation history, but it is the wrong shape to describe the upper envelope of the data. If a rare outlier is invoked to explain massive galaxies at high redshift, the growth tracks predict that the low redshift descendant of this object will be more massive than anything in the local universe.  \nIn contrast, the shape of monolithic model is a good match to the distribution of the data. It traces the envelope of the most massive galaxies at high redshift and arrives at a plausible mass at low redshift. The monolithic model is also consistent with quenched galaxies being common at intermediate redshifts. Objects that form quickly with short star forming timescale ( τ ≲ 1  \nGyr) form the bulk of their stars early (Fig. 4), so should appear quenched.  \nIt appears to us that the problem is the mass assembly history, not the star formation history. If correct, this implies that the hierarchical galaxy formation paradigm is broken; this would be a fundamental challenge to ΛCDM. However, we have considerable freedom to adjust the star formation history of protogalaxies at early times, so we consider this and other possibilities below.", '6.3. Possible Solutions': 'It is conceivable that we have cosmology wrong (Melia 2023; Li 2023). Short of that, there are a number of auxiliary hypotheses that we might invoke to attempt to save the phenomena. In a nutshell, we need to find a way to make the evolution of the brightest progenitor look like a monolithic model (Fig. 1).  \nThere is a strong selection effect to see the brightest objects at the limit of our observations, so perhaps the bright galaxies observed at high redshift are rare outliers (e.g., McCaffrey et al. 2023; Kragh Jespersen et al. 2024). This idea is challenged by the number of objects that exceed the extremal estimates for the single most massive object that can appear as a function of redshift (Fig. 6). It suffers even greater problems explaining the data at intermediate redshift ( z ≈ 3 - 4) where mature, massive ( M ∗ ≈ 10 11 M ⊙ ) galaxies exist (Schreiber et al. 2018; Merlin et al. 2019; Nanayakkara et al. 2024; Glazebrook et al. 2023). These objects were not anticipated to exist in ΛCDM, and cannot be attributed to rare outliers as they are common enough to define L ∗ in Schechter fits (Wylezalek et al. 2014; Steinhardt et al. 2016; Franck & McGaugh 2017). Rare objects may be part of the solution, but rarity by itself does not suffice.  \nAt the highest redshifts, much of the observed light is from the rest frame ultraviolet, which is subject to the considerable uncertainty associated with massive stars and their short lifetimes (Finkelstein et al. 2023). An obvious possibility is a top-heavy IMF (Harvey et al. 2024) that enhances the production of UV light per unit mass. This seems unlikely to maintain the nearly constant UV luminosity function that is observed as the underlying halo mass function evolves rapidly with redshift (section 4.1). This possibility is difficult to test, but it also does not sit well with the observation of quenched galaxies at intermediate redshift that appear to be the descendants of normal stellar populations.  \nSome of the bright sources could be AGN instead of galaxies. This possibility seems unlikely now that spectroscopic observations plainly reveal the spectra of normal stellar populations (Xiao et al. 2023; Price et al.  \n2024; Carnall et al. 2024). That excessively bright sources might be AGN was considered by Franck & McGaugh (2017) who found that the inference of numeroud bright galaxies persisted even if the brightest sources were ignored. Invoking AGN does not really help anyway, as it simply turns the problem of too many early stars into one of too many early supermassive black holes.  \nPerhaps the most obvious possibility is that star formation is more efficient at high redshift so that more of the available baryons form into stars (e.g., Xiao et al. 2023; Carnall et al. 2024; Wang et al. 2024b). This line of reasoning has a limit, as star formation cannot be more efficient than 100%: there comes a point when dark matter halos lack enough baryons to produce the observed stars (Boylan-Kolchin 2023). A few galaxies appear to challenge this disallowed region (Fig. 6), albeit not many. While these may prove illusory (Arrabal Haro et al. 2023), there do seem to be a number of high redshift galaxies for which the stellar fraction approaches the limit of the available baryons ( f ∗ → 1: Xiao et al. 2023; Carnall et al. 2024). This is apparently necessary to maintain consistency with ΛCDM (Carnall et al. 2024; Wang et al. 2024b).  \nOne does not simply turn all the available baryons into stars. Local star formation is nowhere near that efficient (Leroy et al. 2008), and we find it difficult to imagine that it can be. The net result must be to form enough stars early on to mimic a monolithic model as seen in Fig. 10. This is a big ask, as the star formation must not only be highly efficient, but then it must rapidly quench (Kimmig et al. 2023). Such models must also respect the limit on the baryon content of galaxies, which is typically low ( f ∗ < 0 . 5) at low redshift (McGaugh et al. 2010).  \nSome sort of super-efficient star formation (SESF) appears necessary at early times to maintain consistency with ΛCDM. This is a mode of star formation that is utterly unfamiliar in the local universe. We briefly speculate on the conditions in the early universe that may lead to SESF. First, it seems likely that SESF, if it happens, contributes to the formation of giant elliptical galaxies. These may be associated with somewhat rare peaks in the initial power spectrum that preferentially reside in dense environments. Since structure correlates with structure, we imagine that these protogalaxies have a large amount of substructure. In effect, the bottom tier of the merger tree in Fig. 1 starts in a dense configuration. We then imagine that there is some threshold for SESF that is crossed once a sufficient density of substructure is obtained, and the evolution proceeds rapidly  \nin a manner reminiscent of violent relaxation (LyndenBell 1967).  \nThis crude outline is not guaranteed to happen, much less to drive star formation so efficient that f ∗ → 1. To achieve something like SESF, the mode of star formation needs to change dramatically once an arbitrary threshold distinguishing proto-ellipticals from proto-spirals is crossed. For this to work, it is necessary for feedback to be suppressed so that star formation can completely consume the available baryons. This seems outlandish, but is perhaps the least outlandish of the available options.  \nIt is necessary to invoke auxiliary hypotheses like SESF to avoid the conclusion that observations of the high redshift universe are genuinely problematic for ΛCDM (Haslbauer et al. 2022b). Indeed, the usual linear growth rate cannot reconcile the new JWST results with previous results from HST (Sabti et al. 2024) without invoking such extreme hypotheses. Perhaps the data themselves indicate nonlinearity.', '7. DISCUSSION: MOND': "That structure could and should form at an accelerated pace was anticipated well in advance by Sanders (1998), Stachniewicz & Kutschera (2001), McGaugh (2004, 2018), Sanders (2008), Llinares et al. (2008), and others - see McGaugh (2015) and references therein. The new physics driving the prediction of early structure formation is MOND (Milgrom 1983). MOND has a lengthy track record of predictive success (Milgrom 2014), many aspects of which are not satisfactorily explained by dark matter (McGaugh 2020). The early formation of massive galaxies is another predictive success.  \nSanders (1998) was the first to explicitly predict that 'Objects of galaxy mass are the first virialized objects to form (by z = 10).' Contrast this with the contemporaneous ΛCDM statement by Mo et al. (1998): 'presentday discs were assembled recently (at z ≤ 1).' One of these a priori predictions is consistent with the data.  \nThe work of Sanders (2008) anticipates the success of the monolithic model seen above. The observation of bright, massive galaxies at high redshift is precisely what is expected in MOND (Fig. 2). The assembly of mass is greatly accelerated by the nonlinearity of MOND (Nusser 2002); there is no need to invoke SESF or other unlikely effects.  \nTo a good approximation, galaxies in MOND evolve as monolithic island universes after a chaotic early period of hierarchical assembly. This hierarchical assembly is greatly accelerated by the nonlinearity of MOND relative to the linear case of ΛCDM. In situ growth dom-  \nnates and galaxies follow an evolutionary path dictated by the usual astrophysics of gas accretion, cooling, star formation, and feedback (Combes 2014). Mergers certainly continue to happen at late times, but over longer timescales and with less frequency than in ΛCDM (Tiret & Combes 2008).  \nObservations at z ≳ 15 may reveal the rapid early assembly of the first giant galaxies. The transition from hierarchical to monolithic behavior should be sudden, and likely lies in the range 15 < z < 20. However, the precise timing is highly dependent on the underlying cosmology, which need not follow ΛCDM exactly at these redshifts where there are few empirical constraints on the expansion history of the universe and the corresponding time-redshift relation (McGaugh 2018). This depends on the underlying theory, which is may be a generalization of General Relativity along the lines discussed by Skordis & Zglyph[suppress]lo'snik (2019). Exactly what this theory is remains a profound question.", '7.1. Clusters': "A further prediction of Sanders (1998) is that 'larger structure develops rapidly.' For example, massive clusters of galaxies should form early in MOND. A region destined to become a cluster of galaxies will reach maximum expansion after 2 - 3 Gyr (Fig. 2), so clusters should emerge as recognizable objects fairly early in the development of the universe.  \nThe predicted emergence of clusters in MOND is certainly earlier than anticipated in ΛCDM. Kravtsov & Borgani (2012) show that a ∼ 10 15 M ⊙ cluster is barely getting started at z = 3 (see their Fig. 6). Note that in MOND there is no cold dark matter, so the equivalently massive clusters is of order 10 14 M ⊙ . Sanders (1998) predicts 'that by z = 3 not only do massive galaxies exist but they are also significantly clustered (the density of the 10 14 M ⊙ region would be enhanced by a factor of 6.5 over the mean at this redshift).'  \nFranck & McGaugh (2016a,b) identify dozens of protocluster candidates at redshifts 2 < z < 6 . 6 (Fig. 11). Of these, the sixteen most reliable candidates have N ≥ 10 spectroscopically confirmed members with overdensities ranging from 5 to 20 with a median δ = 9 . 5. Similar structures have recently been identified by Laporte et al. (2022), who identify a protocluster candidate at z = 7 . 66, (Morishita et al. 2023) who find another protocluster at z = 7 . 88, and Shah et al. (2024), who identify six massive protoclusters around z ≈ 3. This is consistent with the findings of Franck & McGaugh (2016b) and the predictions of MOND. Clusters should not appear this early in ΛCDM (Mortonson et al. 2011; Kravtsov & Borgani 2012).  \nFigure 11. Measured velocity dispersions of protocluster candidates (Franck & McGaugh 2016a,b) as a function of redshift. Point size grows with the assessed probability that the identified overdensities correspond to a real structure: all objects are shown as small points, candidates with P > 50% are shown as light blue mid-size points, and the large dark blue points meet this criterion and additionally have at least ten spectroscopically confirmed members. The MOND mass for an equilibrium system in the low acceleration regime is noted at right; these are comparable to cluster masses at low redshift.  \n<!-- image -->  \nA number of the best candidate clusters of Franck & McGaugh (2016a,b) are at z ∼ 3. The median velocity dispersion of candidate clusters is ∼ 600 kms -1 (Fig. 11). This is about twice that of equivalent systems found in lookback cones in ΛCDM simulations (Franck 2018). That is, protocluster candidates identified as spikes in N ( z ) have larger velocity dispersions in the data than in simulations at the same redshift. This is another indication that the discrepancy is one in mass, not just luminosity. These systems should not yet be bound in ΛCDM, but it is possible that they have already formed in MOND (Fig. 2).  \nThe median observed velocity dispersion of intermediate redshift cluster candidates corresponds to a mass of ∼ 10 14 M ⊙ for systems in dynamical equilibrium in the deep MOND regime of low acceleration ( ≪ a 0 ). These numbers are very much in line with the long-standing prediction of Sanders (1998), and are similar to the masses inferred for clusters in MOND at low redshift (McGaugh 2015; Li et al. 2023; Tian et al. 2024). Galaxy clusters are known to display a residual mass discrepancy in MOND of about a factor of two (Sanders 2003, 2007; Angus et al. 2008), so there is an apparent need for an extra mass component, possibly in the form of undetected baryons (Milgrom 2008; Kelleher & Lelli 2024). We are thus in the curious situation that the masses of galaxy clusters are problematic for MOND, but their formation time is potentially problematic for ΛCDM, as are other properties like cluster collision speeds (Angus &McGaugh 2008; Katz et al. 2013) and the largest mass objects (Asencio et al. 2021).", '7.2. Larger Structures': "The morphology of cosmic structure is similar in ΛCDM and MOND (Llinares et al. 2008; McGaugh 2015), but the accelerated structure formation of MOND pushes all benchmarks for structure formation to earlier times in cosmic history. This provides a natural explanation for the strong clustering of high-redshift quasars (e.g., Shen et al. 2007; Clowes et al. 2013; Timlin et al. 2018; Eilers et al. 2024; Pizzati et al. 2024) and other very large features (Horv'ath et al. 2014, 2015; Bal'azs et al. 2015; Lopez et al. 2022; Colg'ain et al. 2024). These features are larger than expected in ΛCDM, but quite natural in MOND.  \nAnother indication of early structure is that clusters already appear themselves to be clustered at high redshift (Franck & McGaugh 2016b). This is visible in Fig. 11 as the multiplicity of objects at similar redshift, e.g., at z ≈ 5 . 7 and z ≈ 6 . 6. In addition to explaining unexpectedly overdense regions, MOND is also effective at making large, empty voids (Llinares et al. 2008; McGaugh 2015), a persistent puzzle in ΛCDM (Nusser et al. 2005; Peebles & Nusser 2010).  \nSanders (1998) notes that the largest structures nearing turnaround today would be superclusters (Sankhyayan et al. 2023). These are not expected to be simple objects, as they are inevitably far from the simple spherical approximation of eq. 6. Indeed, it was recognized early (Felten 1984) that MOND would induce anisotropy on unexpectedly large scales. These considerations anticipate the size and complex kinematics of objects like Laniakea (Tully et al. 2014) and Ho'oleilana (Tully et al. 2023).  \nAt late times, it becomes difficult for the entire universe to remain isotropic (Felten 1984) in MOND. This motivates searches for anisotropy in the expansion rate (Colin et al. 2019; Rahman et al. 2022; Kumar Aluri et al. 2023; Bengaly et al. 2024), and would go some way to explain the tension between the dipole anisotropy of number counts of distant sources and the cosmic microwave background (Secrest et al. 2021, 2022; Dom'enech et al. 2022).", '7.3. Early Reionization': "McGaugh (2004) predicted that 'the most obvious signature of MOND-induced structure formation is an early onset of reionization.' This is a natural consequence of early structure formation. The reionization of the universe is achieved by the ultraviolet radiation of conventional sources like Population II stars. It will be patchy and require an extended time to complete, starting around z ≈ 17 (McGaugh 2004, 2018). The uncertainty in the precise redshift of onset is large thanks to the inverse relation between time and redshift. The onset of patchy reionization earlier than z ≳ 12 favors MOND over ΛCDM.  \nJWST observations are beginning to show signs of the predicted early reionization. Tang et al. (2024) find indications of patchy reionization, with some lines of sight being transparent to Ly α photons to surprisingly high redshift ( z ∼ 8). These are the bubbles of early reionization, as expected in MOND.  \nThere is also an apparent crisis in the budget of UV photons at high redshift (Mu˜noz et al. 2024). JWST observations are in tension with the optical depth due to Thomson scattering estimated in fits to Planck data ( τ = 0 . 058 ± 0 . 006 Tristram et al. 2024). This quantity is covariant with other parameters in such fits; a higher optical depth τ ≈ 0 . 17 provides a good fit to the Planck data at ℓ < 600 in the absence of cold dark matter (McGaugh 2004). The high density of UV photons observed by JWST is as expected in MOND. A further consequence of early reionization is an enhanced ISW effect (McGaugh 2004), for which there is some evidence (Granett et al. 2008; Nadathur et al. 2012; Kov'acs et al. 2020).  \nThe early onset of structure formation has further consequences for the high redshift universe. Predictions for 21 cm absorption at cosmic dawn and during the dark ages are discussed by McGaugh (2018). The depth of the absorption signal can be deeper than in ΛCDM at both epochs ( z ≈ 17 and z ≈ 100). The redshift-dependent power spectrum contains further clues. Due to the nonlinear growth of structure, one expects little power in fluctuations entering the dark ages ( z ∼ 150) but more  \nby their end ( z ∼ 50) in MOND than in ΛCDM. These would be clear signs of nonlinear structure formation and a departure from the ΛCDM cosmology.", '8. SUMMARY': 'We have examined evidence concerning the evolution of galaxies across a large range of redshifts for which data are available. There appears to be a population of bright galaxies that formed early and grew rapidly. There is copious photometric evidence that indicates the existence of this population, and important corroborative evidence from kinematics.', '8.1. Photometric Evidence': 'The galaxy population that grew too big too fast has photometric properties that suggest it is  \n- · luminous, with apparent magnitudes considerably brighter than anticipated by contemporaneous ΛCDM models (Fig. 3);\n- · massive, with typical examples approaching the mass of a local L ∗ galaxy already by z ≈ 3 when the universe was only ∼ 2 Gyr old (Fig. 4);\n- · old, with stellar populations consistent with forming at high redshift z f ≳ 10 (Figs. 10 and 7) and quenching early (Fig. 8);\n- · common in clusters at z ≲ 3 and already clustered in protoclusters at z ≲ 6 (Fig. 11); and\n- · consistent with the population of bright galaxies observed by JWST at z > 10 (Fig. 10) where bright galaxies are more common than anticipated (Figs. 5 and 6).  \nThese observed properties do not sit well with hierarchical ΛCDM models, which predict that local massive galaxies were divided into many progenitor protogalaxies at the observed redshifts (Fig. 1). A more natural interpretation is that galaxies are bright at high redshift because they had already grown large. This provides an evolutionary trajectory that maps nicely between observations at low, intermediate, and high redshift. Note that it is not necessary for all galaxies to form early and follow such a trajectory, just enough of them to define a bright L ∗ in early ( z ≈ 3) clusters (Franck & McGaugh 2017). This population was not anticipated by ΛCDM models and is not easily reconciled with them: these galaxies grew too big too fast.', '8.2. Kinematic Evidence': 'Photometric observations are complemented by kinematic observations that trace the mass, not just the  \nlight. JWST observations show that morphologically mature spiral galaxies are common at early times (Ferreira et al. 2022, 2023). Kinematic observations to date (section 5) show that rotationally supported galaxies  \n- · formed early, by z ≳ 6 (Smit et al. 2018; Rowland et al. 2024; Xu et al. 2024) when z ≤ 1 had been the nominal expectation (Mo et al. 1998);\n- · rotate fast, with circular speeds in excess of 250 kms -1 (Rizzo et al. 2021; Lelli et al. 2021), comparable to massive local spirals;\n- · obey the Baryonic Tully-Fisher relation (Lelli et al. 2018, 2019; Nestor Shachar et al. 2023); and\n- · obey the dark matter fraction-surface brightness relation (Fig. 9).  \nThese observations show that at least some spiral galaxies formed early and became massive rapidly. Note that kinematic observations imply a large dynamical mass: it is not just a matter of stars producing more light per unit dark matter halo mass.  \nHigh redshift galaxies appear remarkably mature. Kinematic scaling relations were established early and appear to have evolved little over most of cosmic time (the past ∼ 11 Gyr back to z ≈ 2 . 5). This does not sit comfortably with the gradual assembly predicted by hierarchical galaxy formation (Fig. 1).', '8.3. Structure Formation in MOND': 'The early formation of massive galaxies was explicitly predicted a quarter of a century ago by Sanders (1998). The new physics driving this prediction of accelerated structure formation is MOND, a theory that has had many other predictive successes (e.g., Sanders & Verheijen 1998; McGaugh & de Blok 1998; McGaugh 2011; McGaugh & Milgrom 2013a,b; Milgrom 2015; Sanders 2019a; McGaugh 2020; Mistele et al. 2024a,b). The nonlinearity of MOND causes growth to occur at a much higher rate early than expected with linear growth in ΛCDM (Sanders & McGaugh 2002).  \nMOND makes a number of long-standing predictions about early structure formation:  \n- · early reionization at z ≳ 12 (McGaugh 2004);\n- · massive galaxies at z ≳ 10 (Sanders 1998, 2008);\n- · early emergence of the cosmic web (by z ≈ 5: Llinares et al. 2008; McGaugh 2015);\n- · rich clusters of galaxies form by z ≈ 2 (Sanders 1998);  \n- · an enhanced ISW effect (McGaugh 2004);\n- · large voids swept clear by low redshift (McGaugh 2015);\n- · the largest scales forming at the present time (Sanders 1998) were anticipated to be comparable to Laniakea (Tully et al. 2014);\n- · the universe may depart from the isotropic ideal of the cosmological principle at late times (Felten 1984; Sanders 1998).  \nA number of puzzling observations in cosmology were anticipated by MOND, including the early formation of massive galaxies. The predictive power of MOND is not limited to the dynamics of individual galaxies.  \nDespite the predictive successes of MOND, we do not yet know how to construct a cosmology based on it. In contrast, ΛCDM provides a good fit to a wide range of cosmological observables, but does not provide a satisfactory explanation of the many phenomena that were predicted by MOND (Famaey & McGaugh 2012), nor is it clear that it can do so (Kroupa 2015; Sanders 2019b; McGaugh 2020; Merritt 2020, 2021; Roshan et al. 2021; Haslbauer et al. 2022a,b; Kroupa et al. 2024; Oehm & Kroupa 2024). We find ourselves caught between two very different theories that seem irreconcilable despite applying to closely related yet incommensurate lines of evidence (McGaugh 2015). The simple force law hypothesized by MOND has made enough successful a priori predictions that it cannot be an accident: it must be telling us something. What that is remains as mysterious as the composition of dark matter.  \nThis research has made use of data from a wide variety of sources, themselves enabled by great observatories like the Hubble Space Telescope, Spitzer, and JWST. We gratefully acknowledge the community and society that make these investigations possible. We thank the referee for a detailed reading and many useful suggestions. We also acknowledge conversations with many others, in particular Sara Tosi, Francis Duey, Konstantin Haubner, Tiffany Visgaitis, Tobias Mistele, Pengfei Li, and Marcel Pawlowski.', 'REFERENCES': "Wittenburg, N., Kroupa, P., Banik, I., Candlish, G., & Samaras, N. 2023, MNRAS, 523, 453,  \n- doi: 10.1093/mnras/stad1371  \nWittenburg, N., Kroupa, P., & Famaey, B. 2020, ApJ, 890, 173, doi: 10.3847/1538-4357/ab6d73  \n- Wylezalek, D., Vernet, J., De Breuck, C., et al. 2014, ApJ, 786, 17, doi: 10.1088/0004-637X/786/1/17  \nXiao, M., Oesch, P., Elbaz, D., et al. 2023, arXiv e-prints, arXiv:2309.02492, doi: 10.48550/arXiv.2309.02492  \n- Xu, Y., Ouchi, M., Yajima, H., et al. 2024, arXiv e-prints, arXiv:2404.16963, doi: 10.48550/arXiv.2404.16963  \n- Ying, J. M., Chaboyer, B., Boudreaux, E. M., et al. 2023, AJ, 166, 18, doi: 10.3847/1538-3881/acd9b1\n- Yung, L. Y. A., Somerville, R. S., Finkelstein, S. L., Popping, G., & Dav'e, R. 2019a, MNRAS, 483, 2983, doi: 10.1093/mnras/sty3241\n- Yung, L. Y. A., Somerville, R. S., Finkelstein, S. L., Wilkins, S. M., & Gardner, J. P. 2023, MNRAS, doi: 10.1093/mnras/stad3484\n- Yung, L. Y. A., Somerville, R. S., Popping, G., et al. 2019b, MNRAS, 490, 2855, doi: 10.1093/mnras/stz2755\n- Yung, L. Y. A., Somerville, R. S., Ferguson, H. C., et al. 2022, MNRAS, 515, 5416, doi: 10.1093/mnras/stac2139\n- Zolotov, A., Dekel, A., Mandelker, N., et al. 2015, MNRAS, 450, 2327, doi: 10.1093/mnras/stv740"}