Aims. We extend the analysis of a physical model within the standard cosmology that robustly predicts a high starformation efficiency SFE in massive galaxies at cosmic dawn due to feedbackfree starbursts FFBs. This model implies an excess of bright galaxies at z  10 compared to the standard models based on the low SFE at later epochs an excess that is indicated by JWST observations. Methods. Here we provide observable predictions of galaxy properties based on the analytic FFB scenario. These can be compared with simulations and JWST observations. We use the model to approximate the SFE as a function of redshift and mass assuming a maximum SFE of SUBmaxSUB  0.2  1 in the FFB regime. From this we derive the evolution of the galaxy mass and luminosity functions as well as the cosmological evolution of stellar and starformation densities. We then predict the starformation history SFH galaxy sizes outflows gas fractions metallicities and dust attenuation all as functions of mass and redshift in the FFB regime. Results. The major distinguishing feature of the model is the occurrence of FFBs above a mass threshold that declines with redshift. The luminosities and star formation rates in bright galaxies are predicted to be in excess of extrapolations of standard empirical models and standard cosmological simulations an excess that grows from z  9 to higher redshifts. The FFB phase of 100 Myr is predicted to show a characteristic SFH that fluctuates on a timescale of 10 Myr. The stellar systems are compact RSUBeSUB  0.3 kpc at z  10 and declining with z. The galactic gas consists of a steady wind driven by supernovae from earlier generations with high outflow velocities FWHM  1400  6700 km sSUP1SUP low gas fractions lt 0.1 low metallicities 0.1 ZSUBSUB and low dust attenuation ASUBUVSUB  0.5 at z  10 and declining with z. We make tentative comparisons with current JWST observations for initial insights anticipating more complete and reliable datasets for detailed quantitative comparisons in the future. The FFB predictions are also offered in digital form.
['galaxies: evolution', 'galaxies: formation', 'galaxies: halos', 'galaxies: high-redshift', 'galaxies: ISM', 'galaxies: starburst', 'Astrophysics - Astrophysics of Galaxies', 'Astrophysics - Cosmology and Nongalactic Astrophysics']
Feedbackfree starbursts at cosmic dawn Observable predictions for JWST
{'Observable predictions for JWST': 'Zhaozhou Li 1, ⋆ , Avishai Dekel 1 , 2, ⋆⋆ , Kartick C. Sarkar 1 , 3 , 4 , Han Aung 1 , Mauro Giavalisco 5 , Nir Mandelker 1 , and Sandro Tacchella 6 , 7  \n- 1 Center for Astrophysics and Planetary Science, Racah Institute of Physics, The Hebrew University, Jerusalem, 91904, Israel\n- 2 Santa Cruz Institute for Particle Physics, University of California, Santa Cruz, CA 95064, USA\n- 3 School of Physics and Astronomy, Tel Aviv University, Tel Aviv, 6997801, Israel\n- 4 Department of Space, Planetary & Astronomical Sciences and Engineering, Indian Institute of Technology Kanpur, 208016, India\n- 5 University of Massachusetts Amherst, Amherst, MA 01003-9305, USA\n- 6 Kavli Institute for Cosmology, University of Cambridge, Cambridge, CB3 0HA, UK\n- 7 Cavendish Laboratory, University of Cambridge, Cambridge, CB3 0HE, UK  \nJune 4, 2024', 'ABSTRACT': 'Aims. Weextend the analysis of a physical model within the standard cosmology that robustly predicts a high star-formation e ffi ciency (SFE) in massive galaxies at cosmic dawn due to feedback-free starbursts (FFBs). This model implies an excess of bright galaxies at z ≳ 10 compared to the standard models based on the low SFE at later epochs, an excess that is indicated by JWST observations.  \nMethods. Here we provide observable predictions of galaxy properties based on the analytic FFB scenario. These can be compared with simulations and JWST observations. We use the model to approximate the SFE as a function of redshift and mass, assuming a maximum SFE of ϵ max = 0 . 2 -1 in the FFB regime. From this, we derive the evolution of the galaxy mass and luminosity functions as well as the cosmological evolution of stellar and star-formation densities. We then predict the star-formation history (SFH), galaxy sizes, outflows, gas fractions, metallicities, and dust attenuation, all as functions of mass and redshift in the FFB regime.  \nResults. The major distinguishing feature of the model is the occurrence of FFBs above a mass threshold that declines with redshift. The luminosities and star formation rates in bright galaxies are predicted to be in excess of extrapolations of standard empirical models and standard cosmological simulations, an excess that grows from z ∼ 9 to higher redshifts. The FFB phase of ∼ 100 Myr is predicted to show a characteristic SFH that fluctuates on a timescale of ∼ 10 Myr. The stellar systems are compact ( R e ∼ 0 . 3 kpc at z ∼ 10 and declining with z ). The galactic gas consists of a steady wind driven by supernovae from earlier generations, with high outflow velocities (FWHM ∼ 1400 -6700 km s -1 ), low gas fractions ( < 0 . 1), low metallicities ( ≲ 0 . 1 Z ⊙ ), and low dust attenuation ( A UV ∼ 0 . 5 at z ∼ 10 and declining with z ). We make tentative comparisons with current JWST observations for initial insights, anticipating more complete and reliable datasets for detailed quantitative comparisons in the future. The FFB predictions are also o ff ered in digital form.  \nKey words. galaxies: evolution - galaxies: formation - galaxies: halos - galaxies: high-redshift - galaxies: ISM - galaxies: starburst', '1. Introduction': 'It has been argued from first principles that within the standard cosmological paradigm, the star-formation e ffi ciency (SFE) in massive galaxies at cosmic dawn is expected to be exceptionally high due to feedback-free starbursts (FFBs; Dekel et al. 2023, hereafter D23). SFE at cosmic dawn is predicted to be significantly higher than the low SFE that is valid at later epochs and lower masses, where star formation is assumed to be suppressed by stellar and supernova feedback and / or AGN feedback. The SFE is defined here as the integrated fraction of the accreted gas onto a dark-matter halo that is converted to stars:  \nϵ = M s f b M v , (1)  \nwhere M s is the galaxy stellar mass, M v is the halo virial mass, and f b = Ω b / Ω m = 0 . 16 is the cosmic baryon fraction. An alternative, instantaneous SFE can be defined by relating the star  \nformation rate (SFR) to the gas-accretion rate, ϵ \' = SFR / f b ˙ M v. While at lower redshifts, typically ϵ ∼ 0 . 1 or lower (Behroozi et al. 2019, hereafter B19; see also e.g., Behroozi et al. 2013; Rodríguez-Puebla et al. 2017; Moster et al. 2018; Tacchella et al. 2013, 2018), it is predicted to be of order unity for the most massive galaxies at z ∼ 10 and above.  \nThe simple key idea behind the FFB scenario is that in starforming clouds, there is a natural window that lasts about two million years between a starburst and the onset of e ff ective stellar wind and supernova feedback. Thus, if the gas density in a star-forming cloud is above a threshold of ∼ 3 × 10 3 cm -3 , the free-fall time is ∼ 1 Myr or shorter, allowing the star formation to proceed to near completion in a manner that is free of wind and supernova feedback. D23 showed that several other necessary conditions for FFB are likely to be fulfilled under the specified characteristic high-density conditions. For example, if these starbursts occur in ∼ 10 6 -7 M ⊙ Jeans or Toomre clouds of sizes larger than ∼ 10 pc, the surface density is su ffi ciently high for the ∼ 1 Myr period to also be free of radiative feedback such as radiative pressure or photoionization (Menon et al. 2023; Grudi\'c  \net al. 2021). By an interesting coincidence, at the FFB densities, the cooling time to star-forming temperatures of ∼ 10K is shorter than the free-fall time, as long as the metallicity is not negligible, allowing the short starbursts required for an e ff ective time gap between the burst and the onset of feedback. A unique feature of FFB is that, at metallicities significantly lower than solar, the cooling time becomes shorter than the free-fall time only after the cloud has contracted to a density comparable to the FFB density, delaying star formation until the onset of FFB.  \nFurthermore, it has been shown that the clouds undergoing FFB are shielded against SN-driven winds and UV radiation from post-FFB clusters of earlier generations, provided that the clouds are more massive than ∼ 10 4 M ⊙ . For a galaxy with a stellar mass of ∼ 10 10 M ⊙ in a halo of ∼ 10 11 M ⊙ at z ∼ 10, the FFB is predicted to occur in about 10 4 globular-cluster-like clusters populating a compact, subkiloparsec-sized galaxy, divided into about ten generations during the global halo inflow time of ∼ 80 Myr. The FFB scenario has been analyzed in two extreme configurations: clouds in a spherical collection of shells, or a disk, associated with low or moderate angular momentum in the cold streams that feed the galaxy, respectively, with qualitatively similar results. While D23 suspected that a low metallicity in the star-forming clouds may help to diminish the potentially disruptive e ff ects of stellar winds during the first free-fall times, we argue below that a subsolar metallicity is actually not a necessary condition for FFB. Indeed, while the supernovae in postFFB clusters generate metals and dust, the newly forming FFB clouds, which are shielded against the supernova-driven winds, are not significantly contaminated by metals or dust.  \nD23 evaluated the regime in redshift and halo mass where FFBs with a high SFE are expected. These authors crudely predicted FFBs above the threshold line  \nM v , ff b 10 10 . 8 M ⊙ ∼ 1 + z 10 ! -6 . 2 . (2)  \nThis should be modulated with quenching above a threshold halo mass of ∼ 10 12 M ⊙ or more because of suppression of the cold gas supply through the halos and / or AGN feedback (Dekel & Birnboim 2006; Dekel et al. 2009, 2019). The quenching mass, combined with eq. (2), only permits FFBs above a redshift of z ∼ 6, and then only below the quenching halo mass.  \nThe purpose of this paper is to present predictions of the FFB scenario based on the analytic toy model of FFBs for various observable galaxy properties. These predictions are to be compared with current and especially new JWST observations, as well as future simulations. We compare the FFB predictions with a highz extrapolation of the \'standard\' empirical model, UniverseMachine (UM, B19), as well as with predictions from a \'standard" semi-analytic model (SAM, Yung et al. 2019a,b, 2024a) that was calibrated to match z ∼ 0 observations, which resemble current cosmological simulations where the conditions for FFB are not resolved by the subgrid models. We add tentative, incomplete JWST observation data in certain figures to provide an initial qualitative impression of the compatibility between the FFB scenario and observations. In order to help researchers perform comparisons to other models, simulations, and observations, we provide digital tables and codes for reproducing the predictions online. 1  \nThe paper is organized as follows. In §2 we describe how we represent the FFB scenario with a predicted SFE as a function  \nof redshift and halo mass. In §3 we present the predicted mass function and luminosity function at di ff erent redshifts, and the cosmological evolution of stellar density and SFR density. In §4 we address the predicted bursty star-formation history (SFH) and how it can be constrained by observations. In §5 we introduce a steady wind model for the gas profile in an FFB galaxy, which is then used in the following three sections. In §6 we describe the predicted radii of galaxies at cosmic dawn. In §7 we estimate the gas fraction and metallicity in FFB galaxies. In §8 we evaluate the expected dust attenuation in such galaxies. In §9 we provide a discussion of certain issues. In §10 we summarize our results and outline our conclusions.  \nWe adopt a flat Λ CDM cosmology hereafter with parameter values Ω m = 0 . 3, f b ≡ Ω b / Ω m = 0 . 16, H 0 = 70 km s -1 Mpc -1 , σ 8 = 0 . 82, and n s = 0 . 95, which is close to the Planck cosmology (e.g., Planck Collaboration et al. 2020).', '2. Star-formation efficiency': 'As a basis for deriving observable predictions, we start with formulating the SFE as a function of halo mass and redshift according to the FFB analytic model.', '2.1. Star formation rate': "The distinction between the FFB regime and the nonFFB regime is based on the threshold line M v , ff b( z ) given in eq. (2). Well above this threshold, we assume that all the star-forming galaxies are FFB galaxies, while well below the threshold we assume that they are all nonFFB galaxies, with a smooth transition across the vicinity of the threshold line. We thus take the fraction of FFB galaxies among the star-forming galaxies to be  \nf ff b( M v , z ) = S log[ M v / M v , ff b( z )] ∆ log M ! , (3)  \nwhere the smooth transition is arbitrarily approximated by a sigmoid function, S ( x ) = (1 + e -x ) -1 , varying smoothly from zero to unity, and ∆ log M is assumed quite arbitrarily to be 0.15 dex (with no qualitative e ff ects on the results).  \nThe mean SFR in an FFB galaxy is taken to be  \n⟨ SFR ff b ⟩ = ϵ max f b ˙ M v , (4)  \nwhere ϵ max is the maximum average SFE valid in the FFB regime, and ˙ M v is the halo growth rate, which is approximately  \n˙ M v( M v , z ) / M v ≃ s M β v , 12 (1 + z ) 2 . 5 , (5)  \nwith s = 0 . 03Gyr -1 , β = 0 . 14 and M v , 12 ≡ M v / 10 12 M ⊙ (Neistein &Dekel 2008; Dekel et al. 2013).  \nWe consider the global average maximum SFE ϵ max as a free parameter ranging from 0 . 2 to unity. As demonstrated below, predictions with ϵ max ∼ 0 . 2 appear to fit most of the current observational estimates, while ϵ max = 1 represents the theoretical upper limit for reference. FFB galaxies are predicted below to have a periodic bursty SFH (§4), where the SFE fluctuates about the mean value ϵ max, with peaks of SFE ∼ 1 followed by periods of very low SFE. A value of ϵ below unity can be interpreted as the associated duty cycle, given by the ratio of the characteristic duration of the peak ( ∼ 2 Myr) to the typical duration of each cycle ( ∼ 10 Myr).  \n1  \nFig. 1. Star formation rate and SFE adopted for the FFB model as a function of halo mass and redshift ( M v , z ). Left: SFR. Middle: Instantaneous SFE ϵ ' = SFR / ( f b ˙ M v). Right: Integrated SFE ϵ = M s / ( f b M v). The SFR and the instantaneous SFE are interpolated across the FFB threshold of eq. (2) (solid black curve), from the low SFE based on the standard UM model (B19) well below the threshold, to ϵ max = 1 well above the threshold. A smaller ϵ max gives similar results but with correspondingly lower values for FFB galaxies. The colored area is limited to halos more abundant than 6 σ peaks, with the halo masses of 2 σ and 4 σ peaks shown for reference (dashed lines). The integrated and instantaneous SFEs are extremely similar.  \n<!-- image -->  \n<!-- image -->  \nFig. 2. Star formation e ffi ciency adopted for the FFB model, as in Fig. 1, but via curves of SFE versus M v for a given z (left) and SFE versus z for a given M v (right). The SFE is made to coincide with the standard UM model (B19) (dashed lines) at low masses and redshifts and to approach ϵ max above the FFB thresholds of mass and redshift. The shaded areas refer to a maximum SFE that ranges from ϵ max = 0 . 2 (bottom) to 1 (top). The FFB threshold is at roughly 0 . 5 ϵ max. The SFE does not reach ϵ max for M v ≥ 10 12 . 5 M ⊙ where SFR quenching is assumed.  \n<!-- image -->  \nWe estimate below in §4 that the variation in halo accretion rate and the bursty SFH induce a scatter of σ ∼ 0 . 3 dex in the SFR. Assuming SFR ff b to follow a lognormal distribution, the median log SFR is then log ⟨ SFR ff b ⟩ -0 . 5 σ 2 ln 10.  \nFor nonFFB galaxies, we adopt the empirical model U ni -verse M achine (UM, B19). The UM has been calibrated with pre-JWST observations, including stellar mass functions and quenched fractions (at z < 4), specific and cosmic SFRs ( z < 10), UV luminosity functions (4 < z < 10), and clustering ( z ∼ 0), where each observable is reproduced over the available redshift range. The SFR of star-forming and quenched galaxies, SFRum , sf and SFRum , q, are modeled separately in UM as functions of halo mass and redshift, along with a quiescent fraction f q( M v , z ). As a crude approximation, we extrapolate the UM fitting formulae from Appendix H of B19 to nonFFB galaxies in the relevant mass and redshift ranges at cosmic dawn.  \nAt a given M v and z , the probability distribution of the SFR in star-forming galaxies is assumed to be  \np (SFR | M v , z ) = f ff b p ff b(SFR) + (1 -f ff b) p um , sf (SFR) , (6)  \nwhere the quantities on the right-hand side are all at a given ( M v , z ). As in B19, the SFR of di ff erent types, SFR ff b and SFRum , sf , are assumed to follow separately lognormal distributions, p ff b and p um , sf . The mean SFR for a given halo mass is then ⟨ SFR | M v , z ⟩ = R SFR p (SFR | M v , z ) dSFR.  \nWhen referring to the whole galaxy population, the expression in eq. (6) should be multiplied by (1 -f q), and an additional term, f q p um , q(SFR), should refer to the quenched galaxies. The latter is assumed to dominate above a certain threshold quenching mass, where the SFR is suppressed by a hot CGM and / or the onset of e ff ective AGN feedback (e.g., Dekel et al. 2019). While this quenching mass threshold is very uncertain in B19 due to the lack of constraints at z > 4, we find that the FFB predictions at z > 6 are not sensitive to the quenching details, as most  \nhighz galaxies are well below the quenching threshold and thus expected to be star-forming.", '2.2. Stellar mass': 'In the original UM model, the stellar mass was computed selfconsistently by integrating the in situ SFR and the mass accreted through mergers. However, this is only workable with provided halo merger trees. For simplicity, instead of resorting to detailed merger trees, we use here the median halo growth history. We compute the median stellar mass by starting from the stellar mass of the UM model and adding the time-integrated di ff erence between the mean SFR of the UM + FFB model (SFRum +ff b, §2.1) and the mean SFR of the original UM model (SFRum), namely, 2  \nM s( M v , z ) = M s , um + Z t ( z ) 0 ⟨ SFRum +ff b -SFRum ⟩ d t , (7)  \nwhere the quantities on the right-hand side are at given ( M v , z ) and where M s , um( M v , z ) is taken from the median of the stellarto-halo mass relation in the original UM (B19). 3 We note that the integration on the right-hand side is performed for the typical halo assembly histories M \' v ( z \' | M v , z ). Integrating eq. (5), approximating 1 + z ≃ ( t / t 1) -2 / 3 at z ≫ 1 with t 1 = (2 / 3) Ω -1 / 2 m H -1 0 ≃ 17 . 3 Gyr (Dekel et al. 2013), the mass-growth history of a halo of M v at z is  \nM \' v , 12 ( z \' | M v , z ) ≃ " M -β v , 12 + 3 2 s t 1 β ( z \' -z ) # -1 /β . (8)  \nWe assume p (log M s | M v) to follow a normal distribution. While B19 do not provide a fitting formula for the scatter, we find that it can be approximated by  \nσ s = 0 . 1 + 0 . 3 S (13 . 9 -0 . 3 z ) -log M v 0 . 2 ! , (9)  \nvarying from 0 . 4dex to 0.1dex from the low mass end to the high mass end. 4  \nWe show the SFR and the instantaneous and cumulative SFE (SFR / f b ˙ M v and M s / f b M v) as functions of halo mass and redshift in Fig. 1, and in an alternative way in Fig. 2, respectively. The SFE is made to match the fiducial UM model (B19) in low mass and low redshift halos and it approaches ϵ max in massive halos and high redshifts above the FFB threshold.  \nThe FFB model predicts a strong redshift dependence in the SFE and the stellar-to-halo mass ratio following the evolving FFB threshold halo mass (eq. 2). A cautionary note that should be mentioned is that certain pre-JWST data (Stefanon et al. 2021) indicate only a weaker evolution of the SFE with redshift at z < 10.', '3.1. Stellar mass function': "The stellar mass function is given by  \nΦ (log M s) = Z p (log M s | M v) n ( M v) d M v , (10)  \nwhere the conditional probability distribution in the first term is assumed to be a lognormal distribution with the median given by eq. (7) and the standard deviation by eq. (9), and n ( M v) is the halo mass function in the Λ CDM cosmology, adopted from Watson et al. (2013, with implementation by Diemer 2018) based on simulations that are reliable at z ∼ 10 and beyond (unlike most other similar work). The cumulative number density of galaxies with stellar mass above a threshold value is then N ( > M s) = R ∞ log M s Φ (log M ' s ) d log M ' s .  \nFig. 3 shows the predicted stellar mass function at z = 6-9, 12 . 5, and 16. As an alternative way to present the galaxy abundance, Fig. 4 addresses the stellar-mass distribution via the cumulative number density of galaxies with stellar mass above a threshold value, M s > 10 9 M ⊙ or 10 10 M ⊙ , as a function of redshift. 5 The FFB predictions cover the range of ϵ max = 0 . 2 -1 (areas between the orange and blue lines). They are compared to the standard UM model (B19) and the standard SAM (Yung et al. 2024a). The number density of 10 10 M ⊙ galaxies at z ≥ 8 predicted by the FFB model with ϵ max = 0 . 2 is 1dex higher than the standard models, and even more so at higher redshifts or when assuming a larger ϵ max. Also shown are the first observational results from JWST, tentatively indicating general consistency with the FFB model for massive galaxies at high redshifts for ϵ max ≳ 0 . 2.  \nThe uncertainty in the estimates of stellar masses makes the comparison using the SMF as in Fig. 3 particularly tentative, as indicated by the big scatter among the di ff erent datasets, especially at z = 9. This is because the present stellar libraries, photoionization templates, and the wavelength range covered by the photometry are likely inadequate. For example, the inclusion of the mid-infra bands (MIRI) can make a significant di ff erence in the stellar mass inferred from the spectral energy distributions (SED) (Papovich et al. 2023; Wang et al. 2024). We note that using a flexible bursty SFH in the SED fitting tends to yield more massive galaxies (Harvey et al. 2024), which makes the data agree better with the FFB model with ϵ max = 0 . 2. Despite the uncertainties concerning the SFHs of high-z galaxies (Tacchella et al. 2022b), certain studies suggested that allowing flexible nonparametric SFH may reproduce the true masses more reliably by accounting for outshining (Giménez-Arteaga et al. 2024; Narayanan et al. 2024). Theoretically, the FFB model indeed expects the highz massive galaxies to have bursty SHFs (§4).", '3.2. Luminosity function': 'The UV luminosity function (UVLF) is derived from the halo mass function as  \nΦ ( M UV) = Z p ( M UV | M v) n ( M v) d M v , (11)  \nFig. 3. Stellar mass function. Shown are the predictions of the FFB model at z = 6-9, 12 . 5, and 16. The blue and orange lines refer to ϵ max = 1 and 0.2, respectively. Shown for comparison are the extrapolated standard UM model (B19, solid gray line) and the results of a standard semianalytic model (Yung et al. 2024a, solid dashed line). The FFB model predicts a greater abundance of massive galaxies, particularly evident from M s ≳ 10 9 M ⊙ at z ∼ 9. The excess for lower masses also becomes larger at higher z . The open symbols show pre-JWST observational constraints (Kikuchihara et al. 2020; Stefanon et al. 2021) and the filled symbols refer to tentative JWST constraints (Navarro-Carrera et al. 2024; Wang et al. 2024; Harvey et al. 2024; Weibel et al. 2024). The FFB predictions with ϵ max ∼ 0 . 2 seem to be consistent with the observed mass function at the massive end, though the evidence for deviations of the data from the standard UM model is only marginal for this statistic, especially at z ∼ 9. This is compared to the stronger evidence provided below by other observables. Indeed, we should note that SMF measurements su ff er from large uncertainties in both observation and modeling. For example, assuming a bursty SFH (as expected by the FFB scenario, §4) rather than a smooth SFH during SED fitting results in an increased abundance of massive galaxies (gray squares; Harvey et al. 2024).  \n<!-- image -->  \nFig. 4. Distribution of stellar masses. Shown are the FFB predictions for the number density of galaxies with stellar masses above 10 9 M ⊙ (left panel) or 10 10 M ⊙ (right panel). The blue and orange lines refer to ϵ max = 1 and 0.2, respectively. The standard UM model (B19) is shown for comparison (gray line). The FFB predictions are above the UM model for z ≥ 6, and more so at higher redshifts. Shown for tentative comparison are the observational measurements compiled by Chworowsky et al. (2023), including JWST data (Labbé et al. 2023; Chworowsky et al. 2023) as well as earlier HST data (Song et al. 2016; Stefanon et al. 2021; Weaver et al. 2023). The JWST measurements are systematically higher than the previous measurements and they are above the standard UM model. The observations seem to favor the FFB prediction with ϵ max ≳ 0 . 2.  \n<!-- image -->  \nwhere M UV is the absolute UV magnitude (in the AB magnitude system). Here p ( M UV | M v), the probability distribution of the UV magnitude at a given halo mass and redshift, is assumed to be a Gaussian distribution. For the median, we adopt a simplified relation between the UV magnitude and the stellar mass obtained  \nfrom stellar population synthesis in Yung et al. (2024a, fig. 7), 6  \nmedian { M UV | M v } = -2 . 3 log( M s / 10 9 M ⊙ ) -20 . 5 , (12)  \n2  \nFig. 5. Luminosity function. Shown are the FFB predictions for the rest-frame UVLFs of galaxies at z = 9-13 and 16. The colored bands represent dust-corrected FFB predictions for ϵ max = 0 . 2 (orange) or 1 (blue). The dust correction assumes either the shell or the disk version of the FFB model at the bottom or top bounds of the shaded area, respectively. The UVLFs uncorrected for dust are shown by dashed lines. Shown for comparison is the fiducial UM model (B19) (gray), corrected for dust (solid) or uncorrected (dashed). Also shown is the fiducial SAM simulation without dust attenuation (Yung et al. 2024a, dotted gray). The FFB predictions are well above the standard models, e.g., at z ∼ 9 for M UV ≥-19 and at z ∼ 12 for M UV ≥-18. The open gray data points are pre-JWST constraints (McLeod et al. 2016; Oesch et al. 2018; Morishita et al. 2018; Stefanon et al. 2019; Bowler et al. 2020; Bouwens et al. 2021; Harikane et al. 2022; Finkelstein et al. 2022a; only for z ≲ 10). The colored open symbols are measurements based on galaxies photometrically selected with JWST (Adams et al. 2024; Bouwens et al. 2023a,b; Castellano et al. 2022; Donnan et al. 2023; Finkelstein et al. 2022b; Harikane et al. 2023a; McLeod et al. 2024; Morishita & Stiavelli 2023; Naidu et al. 2022; Pérez-González et al. 2023; Yan et al. 2023; partly complied by Shen et al. 2023). In particular, bright-filled symbols are based on the JWST NGDEEP Epoch 1 (Leung et al. 2023, red), the complete CEERS sample (Finkelstein et al. 2023a, blue), and the JWST spectroscopically confirmed galaxies (Harikane et al. 2024, black). At z ∼ 11, the FFB predictions seem to match the tentative bright-end JWST observations with ϵ max ∼ 0 . 2 -0 . 5. We note that the photometric constraints at z ∼ 16 are highly uncertain and possibly subject to redshift interlopers (Harikane et al. 2024).  \n<!-- image -->  \n]  \nFig. 6. Evolution of the UV luminosity function at M UV = -17 (left) and -20 (right). Similar to Fig. 5, the colored bands represent dust-corrected FFB predictions for ϵ max = 0 . 2 (orange) or 1 (blue) and the gray line shows the fiducial UM model (B19). The UVLFs uncorrected for dust are shown by dashed lines. The FFB predictions at M UV = -20 are higher than the UM model starting at z ∼ 7 and more so at higher redshifts. The open gray data points in the right panel are JWST measurements (Adams et al. 2024; Finkelstein et al. 2023a; McLeod et al. 2024). The tentative comparison indicates a fit with ϵ max ∼ 0 . 2 at z = 8 -12 and possibly higher e ffi ciencies at higher redshifts.  \n<!-- image -->  \nwhere M s is taken from the median stellar-to-halo mass relation in eq. (7), taking into account the enhanced luminosity in the FFB regime. The variance of { M UV | M v } is assumed to be σ 2 M UV = 0 . 3 2 + (2 . 3 σ s) 2 , combining an intrinsic scatter of 0.3 mag in { M UV | M s } (cf. Yung et al. 2024a, fig. 7) and an uncertainty propagated from the scatter in the stellar-to-halo mass relation  \n(eq. 9), 2 . 3 σ s, where the factor 2.3 comes from eq. (12). We note that the high-luminosity end of the UVLF depends on the assumed variance (Shen et al. 2023).  \nThe UV magnitudes can be corrected for dust attenuation A UV( M v , z ), by replacing median { M UV } with median { M UV } + A UV in eq. (11). The attenuation of FFB and nonFFB galaxies are  \n3  \nFig. 7. Cumulative surface number density of galaxies. Shown is the surface number density per unit angular area at redshifts greater than z , for magnitude limits m F277W < 28 . 5 , 29 . 5, and 31.5 in the three panels, respectively. The FFB predictions are for ϵ max = 1 . 0 (blue) and 0.2 (orange), with (solid) or without (dashed) dust correction. They are well above the standard UM predictions (gray). JWST observations are shown in the left and middle panels (gray shaded area), referring to the 68% confidence interval for z > 8 . 5 galaxies from CEERS (Finkelstein et al. 2023a) and NGDEEP (Leung et al. 2023) respectively. The horizontal dotted line indicates the surface density corresponding to one galaxy per JWST NIRCam field of view (9 . 7arcmin 2 ). The tentative comparison indicates a fit with ϵ max = 0 . 2 -1. Future JWST observations (e.g., DESTINY proposed for Cycle 3) can possibly reach down to m F277W ∼ 31 . 5.  \n<!-- image -->  \nFig. 8. Cosmological evolution of comoving densities. Left: Stellar mass. Middle: UVluminosity. Right: SFR. The densities are computed based on galaxies brighter than M UV = -17 . 5 or -17 as indicated in the labels. Shown are the FFB predictions for the bounding values of ϵ max = 0 . 2 and 1 (orange and blue respectively). Shown in comparison is the standard UM model (B19, gray). The shaded bands represent rough uncertainty estimates that correspond to varying the magnitude limits by 0.5 mag up or down. The FFB predictions are well above the UM model, starting with a small excess at z ∼ 8 and growing to an order of magnitude at z ∼ 12. Shown for tentative comparison are observational data (symbols). The FFB model seems to match the tentative data for UV and SFR density with ϵ max ∼ 0 . 2. The references include measurements of stellar density (Madau & Dickinson 2014; Papovich et al. 2023; Santini et al. 2023), UV density, and SFR density estimated from UV (McLeod et al. 2016; Oesch et al. 2018; Adams et al. 2023; Donnan et al. 2023; Harikane et al. 2023a; Pérez-González et al. 2023; McLeod et al. 2024; Finkelstein et al. 2023a; see also e.g., Coe et al. 2013; Ellis et al. 2013; Bouwens et al. 2020, 2023b,a; Robertson et al. 2023a).  \n<!-- image -->  \nincorporated separately, using attenuation models from §8 and B19 (eq 23), respectively.  \nFig. 5 shows the FFB predicted UV luminosity function at z = 9-13 and 16, for ϵ max = 0 . 2 and 1. Shown are both the uncorrected LF and the LF corrected for dust attenuation based on §8. The LF from the FFB model is compared to the LF based on the standard UM model (B19), with and without their correction for dust. It is also compared to the LF in the standard SAM (Yung et al. 2024a) with no correction for dust. The FFB LF at the bright end is significantly higher than the standard predictions, more so at higher redshifts. For example, at z ∼ 12 the excess is higher by more than an order of magnitude already at M UV ∼-20. Current observational results from pre-JWST and from JWST are displayed on top of the theory predictions. While the scatter in the observational estimates is still large, the FFB model seems to provide a qualitative fit to the data, with an SFE that increases with redshift, e.g., indicating ϵ max = 0 . 2 -0 . 5 at z ∼ 11.  \nTo highlight the excess of the luminosity function at high redshifts, Fig. 6 shows the redshift evolution of the UVLF at a given absolute magnitude, M UV = -17 or -20. For example, at  \nM UV ∼-20, the FFB LF is higher than the UM model by more than an order of magnitude already at z ∼ 10, and more so at higher redshifts. Current observational results from JWST are displayed on top. The fit of the FFB model to the data indicates ϵ max ∼ 0 . 2 at z = 8 -12 and possibly higher e ffi ciencies at higher redshifts.  \nFig. 7 shows the FFB-predicted surface number density of galaxies above a given redshift for di ff erent magnitude limits. In order to compute this quantity, we converted the absolute magnitude M UV to the JWST apparent band m F277W using the luminosity distance and the K-correction. A linear relation is assumed between the rest-frame M UV and the observed-frame IR m F277W, whose slope and intercept are computed by fitting to a sample of galaxies in a semi-analytic model simulation (Yung et al. 2019a,b), in which spectral energy distributions (SEDs) were constructed based on the star-formation histories. A tentative comparison to observational results from the complete JWST CEERS sample (Finkelstein et al. 2023a) indicates consistency with the FFB model for ϵ max = 0 . 2 -1. This is also consistent with the peak SFE estimated from the UVLF (e.g., Sun & Furlanetto 2016; Sipple & Lidz 2024).  \n3.3. Evolution of the cosmic stellar mass and SFR densities  \nGiven the stellar mass function from eq. (10) and the luminosity function from eq. (11), the cosmic density of stellar mass and UV luminosity are derived as functions of redshift,  \nρ s( z ) = Z ⟨ M UV | M s ⟩ < -17 . 5 M s Φ (log M s | z ) d log M s , (13)  \nρ UV( z ) = Z M UV < -17 . 5 L UV Φ ( M UV | z ) d M UV , (14)  \nwhere L UV is the luminosity corresponding to M UV. Similarly, the cosmic density of SFR can be estimated as  \nρ SFR( z ) = Z ⟨ M UV | M v ⟩ < -17 ⟨ SFR | M v , z ⟩ n ( M v | z ) d M v , (15)  \nwhere ⟨ SFR | M v , z ⟩ is the mean SFR for a given halo mass and n ( M v | z ) is the halo mass function (Watson et al. 2013). For comparison with observation, the above integrations are limited to galaxies brighter than M UV = -17 . 5 or -17. In observation, ρ UV is commonly used as a proxy for ρ SFR, through ρ SFR ≃ κ UV ρ UV with a constant factor κ UV ∼ 1 . 15 × 10 -28 M ⊙ yr -1 erg -1 s Hz that depends on the assumed initial mass function (IMF) and SFH (Madau & Dickinson 2014).  \nFig. 8 shows the evolution with redshift of the average densities of stellar mass, UV luminosity, and SFR. The FFB model predictions are shown for the bounding values of ϵ max = 0 . 2 and 1. They are compared to the predictions of the standard UM model (B19), showing a significant over-density for z ≥ 8, which is growing with redshift. Current observational results from JWST are displayed on top of the predictions, indicating general consistency with the FFB model, tentatively indicating ϵ max ∼ 0 . 2. We note that the comparison between di ff erent densities and corresponding observations may each lead to a somewhat di ff erent result. This is partly because the di ff erent densities as derived from observations may involve di ff erent assumptions concerning quantities such as the IMF and the SFH, as well as the metallicity and dust obscuration, all a ff ecting the conversion between SFR and UV emissivity.', '4. Star formation history': 'As discussed in §7 of D23, the main phase of star formation in an FFB galaxy is expected to last for a timescale comparable to the halo virial time t vir ∼ 80 Myr (1 + z ) -3 / 2 10 , which is at z ∼ 10 in the same ball park as the time for accretion of most of the baryons (eqs. 31 and 32 there). This long phase is divided to N gen ∼ 10 consecutive generations, each of a period t gen ∼ t vir / N gen. Each generation starts with a peak episode that consists of FFB starbursts in many clusters and consumes gas fast. This peak episode lasts for a duration ∆ t of a few cloud free-fall times ( ∼ a few Myr) and is expected to be followed by a low-SFR period until the newly accumulated gas is su ffi cient for a new episode of fragmentation, bringing the onset of the next generation. The peak SFR is expected to be a factor t gen / ∆ t higher than the average SFR. Based on the accretion rate in eq. 31 of D23, the latter is the SFE ( ϵ ) times the average baryonic accretion rate,  \n⟨ SFR ⟩ ≃ 65 M ⊙ yr -1 ϵ M 1 . 14 v , 10 . 8 (1 + z ) 5 / 2 10 . (16)  \nThe number of generations for the two extreme toy models of shells and disks that were considered in D23 are  \nN gen ≃        7 . 1 M 1 / 3 v , 10 . 8 (1 + z )10 n -1 / 2 3 . 5 (shell) 10 . 6 λ -5 / 6 . 025 M -0 . 05 v , 10 . 8 (1 + z ) -1 / 3 10 (disc) . (17)  \nFor disks this is based on the generation mass given in eq. 60 of D23, assuming Q = 0 . 67, and the total mass f b M v. For shells the generation mass is given in eq. 46 of D23, with R str from eq. 35, assuming T = 10 4 K and stream spin λ s = 0 . 025, and with R v from eq. 23 there. At the threshold for FFB, eq. (2), with n = 10 3 . 5 cm -3 , this becomes  \nN gen ≃ ( 7 . 1 (1 + z ) -1 . 07 10 (shell) 10 . 6 λ -5 / 6 . 025 (disc) . (18)  \nThen, the period of each generation is  \nt gen ≃ t vir N gen ≃ ( 11 . 4 Myr (1 + z ) -0 . 43 10 (shell) 7 . 7 Myr (1 + z ) -1 . 5 10 λ 5 / 6 . 025 (disc) . (19)  \nAt z ∼ 10, the generation period is thus expected to be on the order of 10 Myr, so we consider periods in the range t gen ≃ 5 -20 Myr.  \nThe peak duration reflects the spread in starburst time among the FFB clusters of a given generation. A lower limit would thus be the FFB time ∼ 1 Myr. An upper estimate could be the cloudcrushing time of a typical FFB cloud under the cumulative supernova wind from the clusters of earlier generations, which based on eq. 17 of D23 is ∼ 5 Myr for a 10 6 M ⊙ cluster. We therefore consider peak durations in the range ∆ t ≃ 2 -5 Myr.  \nThe SFH, with the two free parameters t gen and ∆ t , can be constrained by the distribution of the ratio of two di ff erent SFR indicators that are sensitive to di ff erent timescales, based on a sample of FFB galaxies. Here we use the SFR based on common SFR indicators such as H α (or other Balmer lines) and FUV (1350-1750Å), which are mostly sensitive to stars with ages of a few Myr and tens of Myr, respectively.  \nThe upper panels of Fig. 9 show the toy model SFH for different choices of t gen and ∆ t and the associated H α and FUV signals as a function of time. We note that the SFH as indicated by any Balmer line, such as H α , H β , etc., is the same at first order (Tacchella et al. 2022a). The responses of the SFR using Balmer emission lines and FUV emission follow the approach in Caplar & Tacchella (2019). We predict the luminosity evolution for a Simple Stellar Population (SSP) using the Flexible Stellar Population Synthesis code (FSPS; Conroy et al. 2009; Conroy & Gunn 2010). 7 We adopt the MILES stellar library and the MIST isochrones, with a metallicity of 0 . 1 Z ⊙ for the stars and for the gas phase (without a qualitative e ff ect on the results).  \nAs can be seen in the upper panels of Fig. 9, the Balmer peak height is suppressed and its width is extended compared to the intrinsic SFR peak, and these two e ff ects are stronger for the FUV peak. The extension is because the photons are emitted from stars that live longer than the burst width. The duration for the Balmer lines is shorter because they require H-ionizing photons ( λ < 912Å) which are produced by massive short-lived stars, while the FUV photons are also produced by less massive, older stars. The bottom panels of Fig. 9 then show the corresponding distributions of the ratio log SFRBalmer / SFRUV, assuming that the galaxies are observed in a random phase of their periodic SFH. The distribution is narrow at a high ratio above unity when the duty cycle ∆ t / t gen is high. In all cases there is a peak at high SFRBalmer / SFRUV, corresponding to the early times within a generation when H α is higher than FUV. If ∆ t / t gen is su ffi ciently low, and especially when t gen is high, at late times within a generation H α drops to low values while FUV is still  \n1  \nFig. 9. Star formation history in an FFB galaxy at z ∼ 10. The intrinsic SFH consists of several generations, separated by t gen ∼ 10 Myr. Each generation consists of a peak of nearly simultaneous FFB starbursts that last for ∆ t = 2 -5 Myr, followed by a quiet period of gas accumulation until the onset of the following generation of bursts. Shown are three cases where ∆ t = 2 , 5 and t gen = 10 , 20 Myr. Top: Simplified representation of the periodic bursty SFH (gray histogram) and as deduced from two di ff erent SFR indicators, using a Balmer line (H α , H β , etc.; blue) or FUV ( ∼ 1500Å; orange), which probe the SFR at di ff erent time-scales. Bottom: Distribution of the ratio between the SFRs as indicated by Balmer lines and by FUV for galaxies with SFH in the corresponding top panels. Quoted are the standard deviation and the skewness of the distribution, which can be used to constrain the characteristic ∆ t and t gen by observations. For a high duty cycle, ∆ t / t gen, the distribution is dominated by a single skewed peak at a ratio larger than unity, while a long quiet period t gen -∆ t broadens the distribution toward ratios smaller than unity. This characteristic bursty SFH may bias the observed luminosities to even higher values at the bright end.  \n<!-- image -->  \nsignificant, leading to a second peak of the distribution at low SFRBalmer / SFRUV values. As illustrated by the numbers in the lower panels of Fig. 9, one can use two parameters to quantify the shape of the distribution, e.g., the overall standard deviation ( σ ) and the skewness ( γ 1). For example, while σ = 0 . 14 and γ 1 = -0 . 56 for the case ( ∆ t , t gen) = (5 , 10) Myr, they are σ = 0 . 44 and γ 1 = 0 . 09 for (2 , 20) Myr.  \nThe SFH illustrated here is obviously oversimplified, and is meant for qualitative purposes. For example, the stochasticity in gas accretion can introduce an additional scatter of ∼ 0 . 3 dex (Dekel et al. 2013), which may dominate the variations in the SFH. More realistic predictions for the SFH variability will be obtained using cosmological simulations, semi-analytical or hydrodynamical. These will allow predictions for the temporal power spectrum of the SFH (e.g., Caplar & Tacchella 2019; Iyer et al. 2020; Tacchella et al. 2020) and the corresponding observable characteristics.  \nWe recall that the bursty SFH, which is naturally predicted for the FFB phase, is likely to contribute to the enhanced luminosity function at the bright end, where the mass function is dropping steeply, causing low-mass galaxies to preferably scatter to high luminosities (Shen et al. 2023; Sun et al. 2023a,b). This will add to the main e ff ect of enhanced SFE by FFB.', '5. Steady wind and gas profile': 'An important element of an FFB galaxy is the outflowing wind driven by the cumulative supernova ejecta from the earlier generations of star clusters within the galaxy while the outer shell is forming new FFB star clusters. In D23 we crudely modeled this wind based on the mass in each generation of FFB, and used it to evaluate the condition for shielding of the new FFB clouds against crushing by this wind and for estimating the shell radius in balance between the forces exerted by the outflowing wind and the inflowing streams. Here, we model the wind somewhat di ff erently as a steady wind that depends on the average SFR and the associated SFE and not explicitly on the generation mass, thus simplifying the expressions and reducing the dependence on free parameters. This leads to a prediction of the expected observed line width as a function of the SFE, and the steady-state gas density profile as functions of the SFE, SFR, and galaxy radius. We use this model in the following sections to re-evaluate the galaxy radius in the shell version of the FFB model, and to estimate the observable gas fraction, metallicity, and dust attenuation, all as functions of mass and redshift.  \nIn the following analysis, we assume a constant SFE with the mean value as determined by the halo mass and redshift. As discussed in §4, FFB galaxies are expected to have bursty star formation histories with fluctuating SFE, which could influence the results presented below such as the estimates of the galaxy size and dust attenuation. For the purpose of the current qualitative estimates, we tentatively ignore these e ff ects here and defer  \nto future work the incorporation of the e ff ects associated with the fluctuating SFE.', '5.1. Steady wind': 'We consider the outflow generated by supernovae from a stellar system within a sphere of radius R , assuming here for simplicity a constant SFR during the ∼ 100 Myr formation time of the galaxy. The overall rates of mass and energy outflows ( ˙ M sn , ˙ E sn) increase over time as more stars explode in supernovae. The wind from each supernova is strong for about 20 Myr (see Figs. 1, 2 of D23). After such a period, the supernova outflow saturates in a steady wind with outflow rates that are determined by the SFR. Based on Leitherer et al. (1999), assuming Z ∼ 0 . 1 Z ⊙ , these outflows are  \n˙ M sn ≃ 0 . 2 SFR , ˙ E sn ≃ 7 × 10 41 erg s -1 SFR M ⊙ yr -1 . (20)  \nWe assume that the supernova ejecta loads additional interstellar gas into the wind, such that the mass outflow rate is related to SFR via a mass-loading factor η ,  \n˙ M out = η SFR . (21)  \nThe wind velocity is  \nV w ≃ 2 ˙ E sn ˙ M out ! 1 / 2 = 1 , 490 km s -1 η -1 / 2 . (22)  \nThe mass-loading factor η can be expressed in terms of the SFE ϵ , defined here to be  \nϵ = SFR ˙ M in , (23)  \nwhere ˙ M in is the gas inflow rate onto the galaxy. Assuming that the inflowing gas that fails to be converted to stars constitutes the di ff use medium to be loaded onto the supernova wind, we write  \n˙ M out = ˙ M sn + (1 -ϵ ) ˙ M in . (24)  \nUsing eq. (20), eq. (21), eq. (23) and eq. (24), we obtain  \nη = 0 . 2 + 1 -ϵ = ϵ -1 -0 . 8 . (25)  \nϵ  \nBased on eq. (22) with eq. (25), the wind velocity is predicted to be high and growing with the SFE, reflecting the associated dilution of the ISM,  \nV w = ( 3 , 333 km s -1 ϵ = 1 727 km s -1 ϵ = 0 . 2 . (26)  \nThe observed line width at the wings is expected to be roughly FWHM ≃ 2 V w, and the line may be split to two peaks symmetric about the zero velocity, depending on the actual geometry of the inflow and outflow. For high SFE in the FFB phase, where the ISM is rather dilute, these high supernova-driven velocities may be comparable to the velocities expected from AGN-driven winds (e.g., Kokorev et al. 2023; Harikane et al. 2023b). The FFB supernova-driven winds may be distinguishable by the associated lower metallicity (§7.2 below) and lower ionization, in terms of line ratios in the BPT diagram or electron densities (e.g., Kocevski et al. 2023; Maiolino et al. 2024). We note, however, that the line ratios may be quite di ff erent than usual in nonFFB  \nstar-forming galaxies or AGN due to excess emission from the supernova wind at > 70 eV (Sarkar et al. 2022). Being only a function of ϵ , the observed line width at the wings can be used to constrain the maximum SFE in the FFB phase with ϵ ≳ 0 . 2.  \nThe potential e ff ect of radiative losses on the wind is discussed in §8.4, where we find that they are expected to be ine ff ective at high SFE though they may slow down the wind for ϵ ≲ 0 . 2. The possible e ff ect of galactic radiative-driven winds is discussed in §9.3.', '5.2. Gas density profile': 'Given the outflow rates of mass and energy and the galaxy radius R , the steady wind has a mass density profile with di ff erent asymptotic forms interior and exterior to R (Chevalier & Clegg 1985),  \nρ ( r ) =          0 . 3 GLYPH<16> ˙ M out ˙ E GLYPH<17> 1 / 2 ˙ M out R 2 r < R 0 . 05 GLYPH<16> ˙ M out ˙ E GLYPH<17> 1 / 2 ˙ M out r 2 r > R . (27)  \nThe gas density within the galaxy radius R is approximately constant. Inserting the outflow rates from eq. (21), we obtain  \nρ ( r ) ≃        1 . 1 × 10 -2 m p cm -3 η 3 / 2 SFR M ⊙ yr -1 R -2 kpc r < R 1 . 9 × 10 -3 m p cm -3 η 3 / 2 SFR M ⊙ yr -1 r -2 kpc r > R . (28)  \nThe discontinuity at r = R is due to a sharp transition from a subsonic wind within R to a supersonic wind outside, which in practice is materialized by a continuous steep drop near R .', '6. Galaxy radius': 'The typical radii of FFB galaxies can be predicted as a function of mass and redshift in the two extreme cases of disks and shells, based on the estimate of D23 for disks and a slightly revised new estimate based on a steady wind for shells. Realistic galaxies might be mixtures of these two idealized cases.', '6.1. Spherical shell model': 'The outflow comes to a halt when it encounters the inflow and creates a strong shock at a stalling radius R sh, which can be associated with the galaxy boundary. The balance between the ram pressures of the outflow and inflow leads to  \nR 2 sh R 2 str = V w ˙ M out 4 V in ˙ M in = 1 . 13 f ( ϵ ) M -1 / 3 v , 10 . 8 (1 + z ) -1 / 2 10 , (29)  \nwhere the first equality is based on eq. 69 of D23 and the second equality is derived here based on the steady wind model, §5, replacing eq. 70 of D23. Here R str is the e ff ective radius of the streams, ˙ M in = SFR /ϵ and ˙ M out = η SFR are the inflow and outflow rates, V w is the wind speed (eq. 22), V in is the inflow stream speed, assumed to equal the virial velocity, eq. 24 in D23. The factor f ( ϵ ) ≡ ϵ p 5 η = √ ϵ (5 -4 ϵ ) absorbs the dependence on ϵ , with f ( ϵ ) ≃ 1 for ϵ ≳ 0 . 2. Using the e ff ective stream radius from eq. 35 of D23, the shell radius becomes  \nR sh = 0 . 56 kpc × f ( ϵ ) 1 / 2 λ s ,. 025 M 1 / 6 v , 10 . 8 (1 + z ) -3 / 4 10 , (30)  \n0.5  \n<!-- image -->  \nFig. 10. Galaxy size. Left: Galaxy size as a function of halo mass and redshift as predicted by FFB disk model assuming an interim value of ϵ max = 0 . 5. Right: Predicted galaxy half-mass radius as a function of redshift for galaxies at the FFB threshold. The estimates by the two limiting FFB models of shell (blue) and disk (orange), from eq. (31) and eq. (33), are pretty similar, predicting roughly R e ∼ 0 . 3 kpc (1 + z ) -2 . 5 10 . The color bands crudely indicate the uncertainty considering the scatter in the spin parameter. The dependence on ϵ max is weak. Shown for comparison are tentative JWST observations (Baggen et al. 2023; Robertson et al. 2023b; Tacchella et al. 2023a; Finkelstein et al. 2023b; Casey et al. 2024; Sun et al. 2024, McGrath et al. in prep.). The model provides a rough but reasonable estimate of the radius and e ff ectively captures the observed decline of radius with redshift.  \n<!-- image -->  \nwhere λ s ,. 025 ≡ λ s / 0 . 025 is the contraction factor of the stream width with respect to the width of the underlying dark-matter filament (Mandelker et al. 2018).  \nAt the FFB threshold M v , 10 . 8 (1 + z ) 6 . 2 10 = 1 from eq. (2), we obtain  \nR sh , ff b = 0 . 56 kpc × f ( ϵ ) 1 / 2 λ s ,. 025 (1 + z ) -1 . 78 10 . (31)  \nThe e ff ective stellar radius may be very crudely assumed to be R e ≃ 0 . 5 R sh, allowing a contraction to virialization by a factor of two.  \nWe note that the shell radii obtained here using the steady wind model are only slightly di ff erent from those obtained based on eq. 70 of D23. Both are determined specifically for the shell version of the FFB scenario.', '6.2. Disk model': 'In the case of a disk, the typical gas disk radius is derived from the halo virial radius R v and the contraction factor λ to be  \nR d , ff b = λ R v = 0 . 31 kpc λ. 025 M 1 / 3 v , 10 . 8 (1 + z ) -1 10 , (32)  \nwhere R v = 12 . 3 kpc M 1 / 3 v , 10 . 8 (1 + z ) -1 10 from eq. 23 of D23. We note that this disk radius is generic, not special for the FFB scenario.  \nAt the FFB threshold M v , 10 . 8 (1 + z ) 6 . 2 10 = 1, we obtain  \nR disc = 0 . 31 kpc λ. 025 (1 + z ) -3 . 07 10 . (33)  \nThis radius may reflect the stellar e ff ective radius of the disk.  \nOne may crudely assume that the scatter in the size is dominated by the scatter in λ . If λ reflects the halo spin parameter, it is well fitted in Λ CDMcosmological simulations by a lognormal distribution with a dispersion σ ln R e ∼ 0 . 5 in the natural logarithm (Bullock et al. 2001; Danovich et al. 2015; Yung et al. 2024b), though this has not been reliably examined for ∼ 10 11 M ⊙ halos at z ∼ 10.  \nIn summary, based on eq. (30) and eq. (32), we adopt for the galaxy radius (corresponding to ∼ 2 R e) in the two toy models of shell and disk  \nR ≃        0 . 56 kpc f ( ϵ ) 1 / 2 λ s ,. 025 M 1 / 6 v , 10 . 8 (1 + z ) -3 / 4 10 (shell) 0 . 62 kpc λ. 025 M 1 / 3 v , 10 . 8 (1 + z ) -1 10 (disc) . (34)  \nAt the FFB threshold of eq. (2), the associated galaxy radius as a function of redshift is  \nR ≃ ( 0 . 56 kpc f ( ϵ ) 1 / 2 λ s ,. 025 (1 + z ) -1 . 78 10 (shell) 0 . 62 kpc λ. 025 (1 + z ) -3 . 07 10 (disc) . (35)  \nWhile the radii in the two cases are comparable, the estimate for a shell is specific to the FFB model, while the estimate for a disk is rather generic to disk formation.  \nThe predicted galaxy half-mass radius, R e (assumed to be = 0 . 5 R of eq. 34), is shown as a function of halo mass and redshift in Fig. 10 for the toy models of shell and disk with ϵ ( M v , z ) computed as explained in §2.1. Also shown is the half-mass radius of typical galaxies at the FFB threshold of eq. (2) as a function of z (eq. 35), in comparison with a collection of tentative JWST observational measurements. In order to focus on potential FFB galaxies, we display only the galaxies with an estimated stellar mass that is larger than 10% of the FFB threshold halo mass at the corresponding redshift.', '7.1. Gas fraction': 'Integrating over the density profile from eq. (28), the total gas mass in the galaxy is  \nM g = 1 . 16 × 10 6 M ⊙ η 3 / 2 SFR M ⊙ yr -1 R 1 kpc . (36)  \nInserting η ( ϵ ) from eq. (25), SFR from eq. (16), and R from eq. (34), we obtain the gas mass M g for the disk and shell models as a function of ϵ , M v and z . Dividing M g by the stellar mass  \nFig. 11. Gas fraction. Shown is f gas = M g / ( M g + M s) within R ≃ 2 R e as a function of redshift at the FFB threshold, based on eq. (38), for ϵ max = 0 . 5 in the shell and disk models. The shaded bands indicate the uncertainty due to ϵ max ranging from 0.2 (top) to 1 (bottom). The gas fraction is rather low in the whole FFB regime, largely due to the low mass-loading factor associated with the high SFE.  \n<!-- image -->  \nM s ≃ ϵ f b M v, we obtain the mass ratio  \nM g M s ≃        0 . 61 × 10 -2 η 7 / 4 1 . 2 ϵ 1 / 2 0 . 5 λ s ,. 025 M 0 . 31 v , 10 . 8 (1 + z ) 1 . 75 10 (shell) 0 . 60 × 10 -2 λ. 025 η 3 / 2 1 . 2 M 0 . 47 v , 10 . 8 (1 + z ) 3 / 2 10 (disc) . (37)  \nHere ϵ 0 . 5 = ϵ/ 0 . 5, and η 1 . 2 ≡ η/ 1 . 2 refers to ϵ = 0 . 5. For ϵ = 1 or 0 . 2, η = 0 . 2 or 4 . 2 respectively, so the ratio should be multiplied by ∼ 0 . 06 or 6 respectively. At the FFB threshold given in eq. (2), we obtain as a function of redshift  \nM g M s ≃ ( 0 . 61 × 10 -2 η 7 / 4 1 . 2 ϵ 1 / 2 0 . 5 λ s ,. 025 (1 + z ) -0 . 15 10 (shell) 0 . 60 × 10 -2 λ. 025 η 3 / 2 1 . 2 (1 + z ) -1 . 43 10 (disc) . (38)  \nWhen this ratio is small it is a proxy for the gas fraction, f gas = M g / ( M g + M s).  \nFig. 11 shows the gas fraction as a function of redshift at the FFB threshold for the shell and disk models. It is predicted to be rather low, f g < 0 . 1 for ϵ > 0 . 2. A higher SFE leads to a smaller f g due to its dependence on η in eq. (38), reflecting the smaller amount of residual gas left after star formation.', '7.2. Metallicity': "Based on their Fig. 1, showing the Starburst99 mechanical energy injection rate of stellar winds from a 10 6 M ⊙ starburst, D23 expressed the worry that the winds may suppress the SFR if the metallicity in the star-forming clouds is Z ≥ 0 . 2 Z ⊙ . On the other hand, the simulations of Lancaster et al. (2021) indicated that mixing and cooling at the interface of the hot bubble and the cold ISM are likely to make the stellar winds ine ff ective for a few free-fall times even when the metallicity is as high as solar. Even if this conclusion is not valid under all circumstances, one can argue that the stellar winds should not be disruptive because, as we show here, the metallicity in the FFB star-forming clouds is expected to be low. The worry is that metals carried by supernova winds from the earlier generations of stars in the galaxy may contaminate the clouds in the outer shell or ring where the FFBs are expected to occur. This is unlikely for two reasons. First, the mixing of the metals in the clouds is expected to occur on a Kelvin-Helmholtz timescale. As evaluated in §4 of D23 in the context of the cloud shielding against the supernova-driven  \nFig. 12. Metallicity. Shown is the upper limit for the metalicity in FFB star-forming clouds, assuming full mixing between the outflowing metals Z sn = 1 Z ⊙ and inflowing gas Z in = 0 . 1 Z ⊙ (blue) or 0 . 02 Z ⊙ (orange) in eq. (42). The metalicity is computed for halos at the FFB threshold mass. The bottom and top bounds of the shaded areas represent ϵ max = 0 . 2 and unity, respectively. Despite the mixing with supernovaenriched gas from earlier generations, the metallicity remains low and close to Z in, at Z ∼ 0 . 1 Z ⊙ or below. Shown for reference are tentative observations (Curti et al. 2023, 2024; Jones et al. 2023; Hu et al. 2024; also cf Arellano-Córdova et al. 2022), which suggest low but not negligible metalicities in the same ball park, though these galaxies are not necessarily in the FFB regime.  \n<!-- image -->  \nwind, this timescale is expected to be longer than the cloud freefall time once the clouds are more massive than ∼ 10 4 M ⊙ . Second, we show below that even in the event of full mixing, the metallicity in the clouds is expected to remain relatively low, at the level of Z ∼ 0 . 1 Z ⊙ .  \nAssuming that the FFB clouds in the shell are limited to the areas of the shell that come into direct contact with the inflowing streams, only a fraction of the outflowing metals participates in the mixing, f Ω ≡ Ω str / 4 π , where Ω str is the solid angle corresponding to the streams. According to eq. (29), we have  \nf Ω ≡ Ω str 4 π = π R 2 str 4 π R 2 sh = 0 . 22 f ( ϵ ) -1 M 1 / 3 v , 10 . 8 (1 + z ) 1 / 2 10 , (39)  \nwhere f ( ϵ ) ≡ √ ϵ (5 -4 ϵ ) ≃ 1 for ϵ ≳ 0 . 2. This gives in the FFB regime f Ω ∼ 0 . 2 rather insensitively to the value of ϵ .  \nThe mass of gas that outflows into the stream-fed regions of the shell (denoted by ' ) is  \n˙ M ' out = f Ω ˙ M out = f Ω [ ˙ M sn + (1 -ϵ ) ˙ M in] , (40)  \n(cf eq. 24). Similarly, the outflow rate of metal mass into these areas is  \n˙ M ' out , Z = f Ω [ ˙ M sn Z sn + (1 -ϵ ) ˙ M in Z in] , (41)  \nThe first term represents the metals ejected by the supernovae, with ˙ M sn = 0 . 2 SFR = 0 . 2 ϵ ˙ M in (see eq. 20) and metalicity Z sn. The second term represents the gas that came in by inflow, with metallicity Z in, and has not turned into stars.  \nThe outflowing metals and mass are assumed here to mix with the inflows, ˙ M in , Z = ˙ M in Z in. The mixed metallicity is then  \nZ mix = ˙ M ' out , Z + ˙ M in , Z ˙ M ' out + ˙ M in = Z in + 0 . 2 ϵ f Ω 1 + (1 -0 . 8 ϵ ) f Ω ( Z sn -Z in) < Z in + 0 . 2 ϵ f Ω Z sn . (42)  \nFig. 12 shows the mixed metallicity as a function of redshift at the FFB threshold of eq. (2), assuming Z sn ∼ 1 Z ⊙ and Z in ∼ 0 . 1 Z ⊙ for ϵ max = 0 . 2 -1. The mixed metallicity is in the ball park of Z mix ∼ 0 . 1 Z ⊙ at all redshifts, indicating that the enrichment by outflow is rather weak (because of the low f Ω of penetrating cold streams) and quite insensitive to redshift or ϵ max. This is a conservative upper limit for the metallicity in the FFB clouds. The metallicity becomes even lower when assuming lower metallicity in the accreted gas, e.g., Z in = 0 . 02 Z ⊙ . Shown for comparison are tentative JWST observational measurements, which indicate low but non-negligible metallicities in the ball park of the model predictions.  \nWe conclude that as long as the observed gas metallicity is dominated by the gas in the FFB star-forming clouds, from which most of the UV luminosity is emitted, the predicted gas metallicity to be observed in a galaxy during its FFB phase is in the ball park of Z ∼ 0 . 1 Z ⊙ . This is high enough for the cooling time to be shorter than the free-fall time, as required for a burst at the FFB densities, and, in turn, it is low enough to guarantee that the stellar winds are not suppressing the FFB.", '8.1. Attenuation at a given SFR': 'Here, we use the global steady-wind model of an FFB galaxy to estimate its expected dust attenuation as a function of halo mass and redshift. The results are to be used to correct the predicted observable luminosities in the FFB regime, e.g., in the luminosity function.  \nIntegrating over ρ ( r ) in eq. (28) along a radial line of sight in the two regimes inside and outside the galaxy, we obtain the column density of Hydrogen atoms,  \nNH = Z R 0 m -1 p ρ ( r )d r + Z ∞ R m -1 p ρ ( r )d r = (3 . 44 + 0 . 58) × 10 19 cm -2 η 3 / 2 R -1 kpc SFR M ⊙ yr -1 , (43)  \nwhere the first and second terms refer to integration over r < R and r > R respectively.  \nGiven NH , the optical depth at 1500Å is  \nτ 1500 = NH σ ext , mw f d , (44)  \nwhere σ ext , mw ≃ 10 -21 cm 2 / H is the corresponding extinction cross section of the dust per H nucleon in the Milky Way (Weingartner & Draine 2001) and f d is the dust-to-gas ratio in the outflow with respect to its value in the Milky Way. The factor f d in principle has contributions both from the supernovae ( f d , sn) and from the inflowing gas. For low metallicity in the inflow, Z in ≪ Z ⊙ , relevant for FFB, the contribution from the inflow turns out to be negligible (see §8.3), such that  \nf d ≃ 0 . 2 η -1 f d , sn , (45)  \nwith f d , sn = 5 -8 ≃ 6 . 5 for Z ∼ 0 . 1 Z ⊙ in stars (Marassi et al. 2019). Inserting NH from eq. (43) in eq. (44) we obtain  \nτ 1500 ≃ (6 . 89 + 1 . 16) × 10 -3 f d , sn η 1 / 2 R -1 kpc SFR M ⊙ yr -1 . (46)  \nThe corresponding obscured fraction is determined from the optical depth as f obsc = 1 -e -τ . The attenuation in UV is  \nA UV = 2 . 5 τ/ ln(10) = -2 . 5 log 10 (1 -f obsc) . (47)  \n1  \nFig. 13. Obscured fraction f obsc = 1 -e -τ as a function of redshift z along the threshold line for FFB, according to eq. (49). The corresponding extinction A UV is shown on the right y -axis. The toy models of the shell (blue) and disk (orange) are shown, with similar results. The color bands show the uncertainty corresponding to ϵ max = 0 . 2 to 1 about ϵ max = 0 . 5. We see a significant decline in obscuration with redshift. The solid lines are for Z in = 0 . 1, typical in the FFB regime. Shown for comparison (dashed line) are the results for Z in = 1, valid at later times, where the obscuration is stronger. Shown for reference are tentative observations (Curti et al. 2023; Jones et al. 2023; Hu et al. 2024; assuming a Calzetti et al. 2000 attenuation law), which lie between the FFB predictions for shells and disks.  \n<!-- image -->  \nAssuming that the new FFB star formation occurs mostly in a shell or a ring near R , such that the luminosity originates mostly there, the more relevant obscuration is by dust at r > R . In the shell toy model, one could very crudely assume that half the shell closer to the observer is obscured only by the dust in r > R , while the other half is obscured by the dust in r > 0, such that the e ff ective obscured fraction is the average of f obsc , r > R and f obsc , r > 0. In the case of a ring in the outskirts of a disk, for a line-of-sight that does not lie near the disk plane, the e ff ective obscuration may be very crudely estimated by f obsc , r > R alone. For example, for ϵ = 0 . 5, SFR = 65 M ⊙ yr -1 and R = 0 . 5 kpc, we estimate f obsc in the ball park of 0 . 6 -0 . 8.', '8.2. Attenuation in the FFB regime': 'We can further express τ as a function of ϵ , M v and z . The average SFR, based on eq. 31 of D23, is given by eq. (16). Using the galaxy radius from eq. (34) for the two toy models of shell and disk adopting a contraction factor λ. 025 = 1, and inserting these SFR and R in eq. (46), we obtain  \nτ 1500 ≃      (2 . 32 + 0 . 39) f ( ϵ ) 0 . 5 M 0 . 97 v , 10 . 8 (1 + z ) 3 . 25 10 (shell) (2 . 12 + 0 . 36) f ( ϵ ) M 0 . 81 v , 10 . 8 (1 + z ) 3 . 5 10 (disc) , (48)  \nwhere f ( ϵ ) ≡ ϵ p 5 η = √ ϵ (5 -4 ϵ ). For the disk model we only consider r > R , while for the shell model we adopt the average f obsc in the ranges r > 0 and r > R , resulting in an e ff ective prefactor of 0.36 for disks and 0.99 for shells in eq. (48). Interestingly, the dependence of τ on the SFE ϵ is rather weak when ϵ ≳ 0 . 2, with f ( ϵ ) = 1 , 1 . 22 , 0 . 92 for ϵ = 1 , 0 . 5 , 0 . 2 respectively. This is because the increase of SFR with ϵ is roughly balanced by the decrease of η with ϵ .  \nAt the FFB threshold as estimated in eq. (2), we obtain in eq. (48)  \nτ 1500 ≃ ( 0 . 99 f ( ϵ ) 0 . 5 (1 + z ) -2 . 78 10 (shell) 0 . 36 f ( ϵ ) (1 + z ) -1 . 50 10 (disc) . (49)  \nFig. 13 shows the obscured fraction as a function of redshift, based on eq. (49) with ϵ max = 0 . 5, valid for Z in ≪ 1 (solid lines). The shell and disk toy models are compared. We see a significant decline of f obsc with redshift for typical galaxies at the FFB threshold. The tentative JWST observational measurements seem to lie between the FFB predictions for shells and disks, though they do not necessarily correspond to FFB galaxies. The model predictions are also consistent with the low dust suggested by the very blue UV slopes for highz galaxies (Topping et al. 2024; Morales et al. 2024).', '8.3. Attenuation at low redshift and high metallicity': 'For completeness, we comment that at lower redshifts, away from the FFB regime, when Z in (in solar units) is not much smaller than unity, the inflowing gas can make f d grow significantly, and thus increase τ 1500. To see this, define δ to be the dust-to-gas ratio in the di ff erent components. Based on RémyRuyer et al. (2014), δ in the inflowing gas is found empirically to be  \nδ in = Z 1 . 6 in δ mw , (50)  \nwhere δ mw is the ratio at solar metallicity, assumed to be valid in the ISM of the Milky Way. By definition, f d = δ out /δ mw, where δ out = ˙ M dust / ˙ M out. For the dust we write  \n˙ M dust = ˙ M sn δ sn + (1 -ϵ ) ˙ M in δ in . (51)  \nUsing δ sn = f d , sn δ mw based on the definition of f d , sn, as well as eq. (23) and eq. (50), we get  \n˙ M dust = ˙ M sn f d , sn δ mw + ϵ -1 (1 -ϵ ) SFR δ mw Z 1 . 6 in . (52)  \nInserting ˙ M dust in δ out = ˙ M dust / ˙ M out, and using eq. (20) and eq. (21), we finally obtain  \nf d = ˙ M dust ˙ M out δ mw = η -1 [0 . 2 f d , sn + ( ϵ -1 -1) Z 1 . 6 in ] . (53)  \nThe second term is indeed negligible for Z in ∼ 0 . 1 Z ⊙ , namely in FFB conditions, but it dominates if Z in ∼ Z ⊙ and ϵ ≪ 1, valid at low redshifts. The obscured fraction for Z in = Z ⊙ is also shown in Fig. 13 for comparison.', '8.4. The effect of cooling in the wind': 'The analysis so far neglected the possible e ff ect of radiative cooling of the wind within the galaxy, which could, in principle, weaken the outflow and enhance the dust attenuation if the e ff ective cooling radius r cool is smaller than the galaxy radius R . Eq. 6 of Thompson et al. (2016) provides an estimate of the cooling radius in terms of the galaxy radius,  \nr cool ≃ 4 kpc η -2 . 92 R 1 . 79 0 . 3 (0 . 3 SFR10) -0 . 79 , (54)  \nwhere η = ϵ -1 -0 . 8 is the mass loading factor (eq. 25), R 0 . 3 = R / 0 . 3 kpc and SFR10 = SFR / (10 M ⊙ yr -1 ). For a crude estimate,  \nwe insert the galaxy radius in the shell model from eq. 30 and the SFR from eq. (16) in eq. (54) to obtain  \nr cool R ≃ 1 , 418 × (5 ϵ -1 -4) -2 . 72 ϵ -0 . 4 M -0 . 77 v , 10 . 8 (1 + z ) -2 . 57 10 . (55)  \nThis ratio also scales with λ 0 . 79 s ,. 025 , which we assume here to equal unity. For the halos of the FFB threshold mass M h , 10 . 8 = (1 + z ) -6 . 2 10 , this ratio becomes  \nr cool R ≃ 1 , 418 × (5 ϵ -1 -4) -2 . 72 ε -0 . 4 (1 + z ) 2 . 2 10 . (56)  \nWe learn that for the high values of SFE and z that are relevant for most of the FFB regime, the cooling radius is larger than the galaxy radius, so the correction for cooling is small. Cooling may become marginally e ff ective at low ϵ or low z . For example, according to eq. (56), at (1 + z ) = 10, r cool ∼ R for ϵ < 0 . 23, namely cooling may be important at the low bound of the FFB regime. At z ≳ 11, on the other hand, r cool > R even for ϵ = 0 . 2, so cooling is unimportant in the whole range of SFE that is relevant for FFB. In the nonFFB regime, where ϵ ≪ 1, r cool drops below R due to the strong ϵ dependence in eq. (56). When r cool ≤ R , the wind terminates at r < R inside the galaxy, so the gas and dust are confined to radii smaller than R . In this case, if the star formation is mostly in an outer shell or a ring of radius R , the obscured fraction is expected to be in the range 0 -0 . 5 for ring or shell respectively, implying that less than half of the light from the star-forming region is obscured. This is actually comparable to the result for a wind with no cooling.  \nThere are several additional factors that can suppress the radiative cooling in the wind beyond the cooling rate that led to eq. (54). First, Sarkar et al. (2022) showed that the cooling can be suppressed by nonequilibrium ionization in the wind and by excess radiation produced by the hot gas in the star-forming region. Second, entropy can be regenerated in the hot wind by bow shocks created at the interface between the wind and cold clouds (Conroy et al. 2015; Nguyen & Thompson 2021; Fielding & Bryan 2022), though it is possible that mixing of cold and hot gas near these clouds may overwhelm this increase in entropy. Based on all the above, we can conclude that the estimates of gas fraction, metallicity, and dust attenuation that neglected cooling in the wind are viable approximations in the FFB regime.', '9.1. Additional mechanisms for bright galaxy excess': "Beyond the enhanced SFE due to FFBs presented in D23 and discussed here, there are several additional mechanisms that can potentially contribute to the excess of bright galaxies at cosmic dawn compared to extrapolations of the standard model and current simulations. Four di ff erent general types of additional mechanisms come to mind as follows.  \n- · Scatter in the luminosity at a given halo mass, e.g., due to burstiness in the star-formation rate, which preferentially scatters low-mass halos to higher masses at the massive end, where the halo mass function is declining steeply (Sun et al. 2023b; Yung et al. 2024a).\n- · Enhanced ratio of UV luminosity to stellar mass, e.g., by a topheavy IMF, or by a nonstellar UV source such as AGN (Yung et al. 2024a; Tacchella et al. 2023b; Grudi'c & Hopkins 2023).  \n- · Reduction in dust attenuation at higher redshifts (Ferrara et al. 2023; Ferrara 2024).\n- · Modifications of the standard Λ CDM cosmology that are associated with an enhancement in the abundance of massive halos, such as the introduction of Early Dark Energy (Klypin et al. 2021) or enhanced primordial power spectrum (Parashari &Laha 2023).  \nThe enhanced SFE is robustly predicted by the physical model of FFB, which is naturally expected preferentially at the high densities and low metallicities present in high-mass halos at cosmic dawn. This mechanism alone can enhance the observed luminosities by an order of magnitude. As shown in §4, the FFB scenario naturally involves a characteristic bursty SFH at the bright end, which could bias the luminosity further to higher values. Furthermore, as estimated in §8, the dust attenuation in the FFB regime is expected to be low and to decrease further with redshift, which also helps to make the galaxies brighter.  \nThe consideration of a nonstandard IMF is currently under both theoretical and observational uncertainties, potentially either suppressing or enhancing the FFB mechanism and the abundance of early bright galaxies (see a more detailed discussion in §9.2). The involvement of AGN in enhancing the UV luminosity is also conceivable but remains rather uncertain and yet to be explored; tentative indications suggest that the contribution of AGN to the UV at cosmic dawn is subdominant (Scholtz et al. 2023). Finally, turning to modifications of Λ CDM would necessarily involve ad hoc free parameters, potentially conflicting with the spirit of the simplicity principle of Occam's razor.  \nOverall, it is important to come up with observable features that could constrain the contribution of each mechanism to the excess of bright galaxies. The main observable distinguishing feature of the FFB scenario is the preferred appearance of high SFE at high redshifts and possibly at high halo masses. Such a mass dependence does not seem to emerge as a natural intrinsic feature if a top-heavy IMF is responsible for the bright excess, or if variations in dust attenuation are responsible for the redshift dependence of the excess (see Finkelstein et al. 2023a).", '9.2. Initial mass function': "The IMF at cosmic dawn is unknown. Given both the theoretical and observational uncertainties, it could deviate from the standard low-redshift IMF either way, towards a top-heavy or a bottom-heavy distribution. Such deviations could either suppress or enhance the FFB nature of galaxy formation and the excess of early bright galaxies.  \nOn one hand, there are observational indications that the stellar populations in the cores of massive elliptical galaxies at low redshifts have a bottom-heavy IMF (Conroy & van Dokkum 2012; van Dokkum et al. 2017; Hallakoun & Maoz 2021). If these are the descendants of the early FFB galaxies, in agreement with the average halo mass growth of three orders of magnitude from z ∼ 10 to the present, it would indicate that the FFB population had a bottom-heavy IMF. In agreement, simulations that include radiative and proto-stellar outflow feedback find that these processes tend to induce an IMF with a lower characteristic stellar mass and an overall narrower mass range in high surface density environments (Tanvir et al. 2022), typical of the starforming clouds at cosmic dawn. When the IMF is bottom-heavy, the deficiency in OB stars would lead to suppressed supernova and stellar feedback, enhancing the FFB nature of galaxy forma-  \non at cosmic dawn even beyond the direct primary e ff ects of high densities.  \nOn the other hand, a top-heavy IMF may be expected at cosmic dawn due to the high cosmic microwave background (CMB) temperature or the high ISM pressure and density, as indicated in other simulations (Grudi'c & Hopkins 2023; Bate 2023; Yung et al. 2024a). A top-heavy IMF is also indicated under certain assumptions from observations, e.g., in two nebular-dominated galaxies at z ∼ 6 -8 (Cameron et al. 2023) or by fitting the photometry of high-redshift galaxies with revised templates (Steinhardt et al. 2023). By enhancing the UV luminosity-to-mass ratio, a top-heavy IMF could contribute to an increased abundance of bright galaxies. However, typical simulations indicate only a limited factor of 2-3 in luminosity enhancement (Yung et al. 2024a), which may not be su ffi cient by itself to explain the full excess of bright galaxies. Furthermore, the excessive radiation pressure associated with a top-heavy IMF will tend to suppress the SFE and thus act to reduce the excess of bright galaxies, leaving the net e ff ect uncertain. Simulations by Menon et al. (in prep.) find that a top-heavy IMF with a luminosity-to-mass ratio that is larger by an order of magnitude than the standard could reduce the SFE by a factor of two, but the suppression is smaller at the low metallicities and dust content expected at cosmic dawn, where the radiative feedback in high-surface-density clouds is largely by IR radiation pressure on dust. This indicates that a top-heavy IMF and a high SFE can co-exist.  \nIn parallel, a top-heavy IMF can generate FFB in an alternative way, by accelerating the core collapse to black holes at the centers of the star-forming clusters (Dekel et al., in prep.). This process is driven by the inward migration of the massive stars due to dynamical friction, which would be naturally enhanced when more massive stars are present (Portegies Zwart & McMillan 2002; Rizzuto et al. 2023). The core collapse of clusters is further sped up by the rotation and flattening of the clusters that form in galactic disks (Ceverino et al. 2012), via a process of gravo-gyro instability (Hachisu 1979; Kamlah et al. 2022). In clouds where this core collapse occurs within a few Myr, it is expected to eliminate the massive stars before they can generate e ff ective feedback at the end of their lifetimes, thereby leading to FFB. From the above discussion, we conclude that the IMF at cosmic dawn and its implications on FFB and the excess of bright galaxies remain largely uncertain, necessitating more thorough exploration.", '9.3. Radiatively driven versus supernova-driven wind': "The analysis of gas fraction, metallicity and dust in §7 and §8 is based on an assumed supernova-driven wind in steady state, as presented in §5. This wind is assumed to be driven by supernovae from earlier generations of stars, in their active supernovae phase, during most of the ∼ 100 Myr duration of the FFB galaxy formation phase. Here we comment on the validity of this supernova-driven wind in view of a possible global gas (and dust) ejection from the galaxy by super-Eddington radiativedriven winds, as proposed by Ferrara (2024). Based on eq. 2 of Ferrara (2024), super-Eddington ejection is expected to occur once the specific SFR in the galaxy is above a threshold,  \nsSFR > sSFRcrit ≃ 25 Gyr -1 100 A ! 2 f bol ! f ts . (57)  \nHere A is the ratio of dust to Thomson scattering cross-section and f bol is the bolometric correction with respect to the UV luminosity at 1500Å, both set to their fiducial values. We introduce  \nthe correction factor f ts, the ratio of total to stellar mass within the galaxy, to account for the neglected contribution of the nonstellar mass to the Eddington luminosity. With a gas-to-stellar mass ratio of ∼ 0 . 5 and a moderate additional contribution of dark matter one can assume f ts ∼ 2. However, this factor may be balanced by a larger value of A , so we conservatively ignore it.  \nFerrara (2024) assumes that the radiation dominates over the supernovae in driving the wind. Tentatively adopting this hypothesis, we evaluate the redshift range of validity of the condition in eq. (57). Based on Dekel et al. (2013), the analytic estimate in their eq. 7 and their simulation results, we approximate the sSFR by the specific accretion rate onto the halo,  \nsSFR ≃ 6 . 9 Gyr -1 (1 + z ) 5 / 2 10 M 0 . 14 v , 11 . (58)  \nThis redshift dependence is consistent with the findings of Zhao et al. (2009) and it has been shown to be a better and more physical approximation than the crude estimate that is used by Ferrara (2024) based on the ratio of halo mass to halo crossing time (see the di ff erence between eqs. 31 and 32 of D23). Inserting the sSFR from eq. (58) in eq. (57), one expects super-Eddington ejection at  \n1 + z > 1 + z crit ≃ 16 . 8 M -0 . 056 v , 11 f 0 . 4 st ( A / 100) -0 . 4 ( f bol / 2) -0 . 4 . (59)  \nAssuming the fiducial values of the parameters, this implies a conservative lower limit of z crit ≃ 16 for halos of M v ∼ 10 11 M ⊙ . If we adopt the weak mass dependence in eq. (58), and insert the declining M v , ff b( z ) at the FFB threshold from eq. (2), we can solve eq. (59) for z to obtain z crit ≃ 22. The value of z crit = 16, and certainly z crit = 22, are beyond the redshift range of the current JWST observations in which the excess of bright galaxies is detected, making the radiative-driven global winds less relevant for the main FFB phase.  \nAt these extremely high redshifts when super-Eddington winds prevail, we return to the question of whether it is justified to neglect the SN-driven winds in the first place. A quick comparison with the equivalent pressure produced by the SNe and radiation shows that these pressures are almost equivalent (Thompson et al. 2005; Murray et al. 2011). However, the SNedriven pressure may be reduced due to radiative cooling such that the SNe pressure becomes much lower, and an energy-driven solution may not be viable for individual SNe. In such cases, the common wisdom is that a momentum-driven wind dominates. Even in this case, it is found that the wind-driving forces exerted by supernovae and by radiation from stars are comparable, at ˙ p / m ∗ ∼ 6 × 10 -8 cms -2 (Grudi'c et al. 2020), and both are comparable to the additional contribution from stellar winds. We recall that the momenta from clustered sources, as expected in particular in the FFB compact clouds in compact galaxies at z ∼ 10, add up super-linearly for supernovae (Gentry et al. 2017), but linearly for radiative pressure, making the supernovae more important. At significantly clustered star formation regions (as is the case for the high redshift galaxies under consideration), the wind becomes energy-driven again. Such highly clustered SNe can retain 20-100% of the SNe energy (Fielding et al. 2018), as estimated observationally in M82 (Strickland & Heckman 2009).  \nFurthermore, the pre-blast wave energy losses of SN bubbles due to their interaction with the ISM, highlighted in Ferrara (2024), are relevant only when the ISM is dense ( ≳ 10 5 m p cm -3 ; Terlevich et al. 1992), as it is in the first free-fall times of the first starbursts, prior to the onset of e ff ective SN feedback. Soon after, the ISM becomes dilute (in gas density, metals and dust) due to both the high SFE of conversion of gas into stars and the  \nsuper-Eddington winds. The SN-driven steady wind (§5) keeps the ISM dilute for the rest of the ∼ 100 Myr global FFB phase (§7). In such a dilute ISM, the radiation-driven wind is expected to become negligible, as there is only little dust to respond to the radiation pressure. In turn, the dilute ISM allows the SN-driven wind stream freely and be as e ff ective as it could. In more detail, there could be repeating episodes of a dense ISM (dusty and red) followed by a dilute ISM (dust-free and blue) in a characteristic duty cycle, as eluded to in Ferrara (2024). This is as predicted in our FFB scenario, §4 and D23. During the dense episodes, one may expect bursty super-Eddington ejection on top of the steady SN-driven wind, but typically only at the very high redshifts indicated by eq. (59).  \nWe conclude that super-Eddington ejections are potentially valid only at very high redshifts beyond the current main observed range, and even when they are, they are not likely to overwhelm the steady supernova winds that we use as a basis for our estimates of gas fraction, metallicity and dust attenuation. Nevertheless, based on the SN-driven winds, we do predict low dust attenuation at high redshifts, which becomes negligible at very high redshifts, where it is also predicted to be removed in episodes of super-Eddington winds.  \nWhile Ferrara (2024) argues that being dust-free is enough to bring the observed luminosities and SFRs to the level of the JWST observations, this is not the case in massive halos at cosmic dawn based on the tentative comparisons shown in figures above between JWST observations and the standard nonFFB models; for example, as represented by the UM empirical model (B19) naively extrapolated to high redshifts and the standard SAM (Yung et al. 2024a), both with no dust correction. Being dust-poor may help, especially at very high redshifts, but the main contribution to the excess of bright galaxies and high SFR at cosmic dawn is likely to be the enhanced SFE because of FFB, helped by the associated bursty SFE, maybe by a somewhat topheavy IMF, and possibly by UV from AGN.", '9.4. Choice of \'standard" UM model': "A word of caution is called for concerning the use of the empirical UM model (B19). We used it to normalize the SFE at low redshifts and low masses away from the FFB regime. Then, in certain occasions, we compared our predictions for FFB with extrapolations of the UM model to the relevant high redshifts. One should be aware of the fact that the current version of the UMmodel has been derived based on observations at lower redshifts, and the results with which we were comparing are based on uncertain extrapolations to z ∼ 10. These 'predictions' of the UMat high redshifts should therefore be considered with a grain of salt. This uncertainty in the UM model does not a ff ect our predictions in the FFB regime, in which the UM predictions are not involved. The UM can eventually be revised to refer to z ∼ 10 by fitting the model to new JWST observations and thus recover the excess of bright galaxies at high redshifts.", '9.5. FFB at lower redshifts': 'In principle, the FFB phenomenon may occur and be observable at redshifts and / or masses below the FFB threshold of eq. (2), with decreasing abundance further below the threshold line. These FFBs can naturally arise in our simple FFB model from the scatter in the parameters that control the density in the star-forming regions, which, once in the tails of their distribu-  \ntions, may bring the density to above the required threshold for FFB. For example, this may occur in disks of particularly small sizes that derive from small values of λ , or where the gas inflows are denser and / or more supersonic than the average for the given halo mass. Alternatively, di ff erent mechanisms, such as wet galaxy mergers, intense instreaming, or other compaction events, may drive the required high densities and generate starbursts with high SFEs even at later times (Krumholz et al. 2012; Zolotov et al. 2015; Turner et al. 2015; Lee et al. 2016; Usero et al. 2015; Fisher et al. 2022; Querejeta et al. 2019; Lapiner et al. 2023). These later high-SFE starbursts, however, may not necessarily follow the \'standard" FFB scenario worked out above and may therefore not share all the properties that were predicted for highz FFB galaxies.  \nDensities above the FFB threshold of n crit ∼ 3 × 10 3 cm -3 apparently appear in dense regions within giant molecular clouds (GMCs) in low-redshift disk galaxies. However, contrary to the FFB clouds, the GMCs are believed to have low overall SFEs that are likely suppressed by feedback, which actually tends to disrupt the clouds in a few dynamical times (e.g., Krumholz 2017). Two obvious distinguishing features between GMCs and highz FFBs are the metallicity and the background disk density. The metallicity is near solar in GMCs while it is well below solar in FFBs, and the background disk density is much lower than the FFB density in GMCs while it is not far below the FFB density in FFBs. In GMCs, the low disk density makes the cloud contraction (either at the original cloud formation or at its reformation after disruption) be slow, on the disk free-fall timescale, which is longer than the ∼ 1 Myr time window for FFB. During this contraction, the high metallicity makes t cool < t ff already in the early phases of cloud contraction, allowing star formation to occur before the cloud develops the dense core, thus leaving no time window free of supernova feedback. In FFBs, on the other hand, the high disk density makes the collapse rapid, on a cloud free-fall timescale of ∼ 1 Myr, and the low metallicity makes the cooling time drop below the free fall time and allow star formation only once the density exceeds the critical FFB density. This permits a burst free of feedback.  \nFurthermore, based on measurements of HCN / CO ratio, the regions in GMCs where the density is above n crit, which tend to be filamentary structures, are estimated to occupy only ∼ 10 -3 of the GMC volume in local disks, with a clear trend upwards with larger galactic-averaged density (Mark Krumholz, private communication; based on, e.g., Jiménez-Donaire et al. 2019). This makes the surface density in local GMCs of size < 100 pc lower than the critical surface density of Σ crit ∼ 3000 M ⊙ pc -2 required for preventing significant gas ejection by stellar radiation pressure (Menon et al. 2023), but probably not so in high-redshift FFB clouds that reside in much denser background disks. Therefore, the SFR in local GMCs is not expected to be free of radiative feedback the way the highz clouds are.', '10. Conclusion': "As proposed analytically by D23, feedback-free starbursts (FFB) are expected in gas clouds where the gas density exceeds a critical value of n ∼ 3 × 10 3 cm -3 . This ensures that the free-fall time is ≤ 1 Myr, which is shorter than the timescale for feedback from supernovae and stellar winds, and also ensures that the cooling time is shorter than the free-fall time, allowing a starburst to occur on a free-fall timescale. Moreover, the gas surface density in the FFB clouds should be above Σ ∼ 2 × 10 3 M ⊙ pc -2 to guarantee that radiative feedback is ine ff ective.  \nConsequently, enhanced SFE, of namely SFE ∼ 0 . 2 -1, is expected in halos above a threshold in mass and redshift, roughly defined by the curve M v , 10 . 8 ∼ (1 + z ) -6 . 2 10 , leading to distinguishable observable predictions, as follows.  \n- · The average cosmological SFR density should be enhanced compared to the 'standard' low-SFE scenario -a small overdensity at z ∼ 8 and a large overdensity at z ∼ 12.\n- · The galaxy mass and luminosity functions are expected to be correspondingly enhanced at z ≥ 8, and gradually more so at higher masses and brighter magnitudes. For example, at z ∼ 12, galaxies of UV magnitudes of -20 are expected to be overabundant by more than an order of magnitude compared to the extrapolated standard empirical model of low SFE (B19).\n- · The SFH, during the global galaxy formation time of ∼ 100 Myr, is expected to be bursty on a timescale of ∼ 10 Myr. In each such generation, a peak in SFR is expected to be followed by a low-SFR phase. This duty cycle should be translated to a global average SFE that is smaller than unity, possibly as small as 0.2. The predicted bursty SFH should be detectable by the distribution of the ratio of H α (or other Balmer lines) to UV SFR indicators. This burstiness should also contribute to the enhanced observed luminosities at the bright end.\n- · The FFB galaxies are expected to be compact, with typical stellar radii of R e ∼ 0 . 3 kpc (1 + z ) -2 . 5 10 .\n- · The gas in a galaxy during its FFB phase is expected to form a steady wind driven by supernovae from post-FFB star clusters of earlier generations. The wind velocity is expected to be V w ∼ 3 , 300 km s -1 for ϵ = 1 and ∼ 700 km s -1 for ϵ = 0 . 2, observable through broad line wings of FWHM ∼ 2 V w, which could be confused with AGN-driven winds. The supernova-driven winds could possibly show two peaks centered on zero velocity, as well as low metallicity and low ionization in BPT-like diagrams.\n- · Using the steady density profile associated with the wind, we predict (a) a low gas fraction of approximately f g ≤ 0 . 05 (1 + z ) -1 10 , (b) low metallicity in the FFB star-forming clouds of Z ∼ 0 . 1 Z ⊙ , and (c) low attenuation by dust with obscured fractions of ∼ 0 . 4 and 0 . 2 at z = 9 and 15 , respectively. This diminished attenuation is also expected to contribute to the enhanced observed luminosity at high redshifts.  \nIn addition, it is worth reiterating several other observable predictions made by D23.  \n- · A galaxy at the end of the FFB phase is predicted to be a compact assembly of thousands of star clusters, typically of ∼ 10 6 M ⊙ each. These clusters form in gas clouds of ≥ 7 pc in radius.\n- · The feeding of the galaxy is expected to be through inflowing cold streams from the cosmic web, with no shock-heated gas in the halo, and minor outflows.  \nSeveral galaxy properties are not significantly constrained by the FFB model. For example, the global morphology could be of clumpy disks or clumpy spherical shells or an irregular morphology in between. As another example, an ad hoc top-heavy IMF and / or UV from AGN may also contribute to the enhanced luminosities, but these factors are not necessary to explain the bright excess.  \nOngoing and future studies of the FFB scenario and its implications include the following items:  \n- · In Libanore et al. (2023), we show that the FFB phase is expected to be detectable in the power spectrum of the 21cm signal, which will be observed by the Hydrogen Epoch of Reionization  \nArray (HERA, DeBoer et al. 2017) at z ≥ 15, but the net e ff ect on the reionization history is expected to be negligible.  \n- · The FFB star clusters could be the natural sites of massive black hole (BH) seeds due to core collapse that occurs in a few free-fall times, and is sped up by the broad stellar mass spectrum spectrum in the young clusters and the possible cluster rotation and flattening. While the stellar clusters disrupt each other rather quickly, these seeds, which populate compact galactic disks, merge into supermassive black holes with high BH-tostellar mass ratios (Dekel et al. in prep.). This post-FFB phase is yet to be properly investigated and simulated.\n- · The FFB scenario and its observable predictions are to be studied using cosmological simulations in order to deepen, sharpen, and quantify our theoretical understanding of galaxy formation at cosmic dawn. This can be done via semi-analytic simulations in which FFB is incorporated by subgrid models. Ultimately, it should be done using full hydro-gravitational cosmological simulations, which resolve the processes of star formation and feedback that should enable FFB.  \nIn summary, we obtained a large body of observable predictions, which should help us to distinguish the robust FFB scenario from other potential mechanisms that could contribute to the overabundance of bright galaxies at high redshifts. Tentative, qualitative comparisons to early observational results from JWST seem to indicate qualitative agreement with the FFB predictions. These predictions should be quantitatively compared to upcoming observations by JWST.  \nTo help others perform comparisons to other models, simulations, and observations, we provide the digital tables and code to produce the predictions online. 8 Data and results underlying this article will be shared upon reasonable request to the corresponding author.  \nAcknowledgements. We are grateful for stimulating discussions with Josephine Baggen, Guilermo Barro, Peter Behroozi, Caitlin Casey, Katherine Chworowsky, Andrea Ferrara, Steven Finkelstein, Yuichi Harikane, Thomas Harvey, Mark Krumholz, Elizabeth McGarth, Shyam H. Menon, Yongzhong Qian, Brant Robertson, Xuejian Shen, Rachel Somerville, Romain Teyssier, Andrea Weibel, Haojing Yan, Zhiyuan Yao, and Aaron Yung. ZL thanks Hengxiao Guo, Peng Wang, Jiaxin Han, and Rachel Somerville for hosting his visit at SHAO, SJTU, and CCA / Flatiron during the preparation of this work. This work was supported by the Israel Science Foundation Grants ISF 861 / 20 (ZL, AD) and 2190 / 20 (KCS). ZL acknowledges the funding from the European Union's Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant 101109759 ('CuspCore'). HA acknowledges support from the Zuckerman Postdoctoral Scholars Program.", 'References': '```\nAdams N. J., et al., 2023, MNRAS, 518, 4755 Adams N. J., et al., 2024, ApJ, 965, 169 Arellano-Córdova K. Z., et al., 2022, ApJ, 940, L23 Baggen J. F. W., et al., 2023, ApJ, 955, L12 Bate M. R., 2023, MNRAS, 519, 688 Behroozi P. S., Wechsler R. H., Conroy C., 2013, ApJ, 762, L31 Behroozi P., Wechsler R. H., Hearin A. P., Conroy C., 2019, MNRAS, 488, 3143 Bouwens R., et al., 2020, ApJ, 902, 112 Bouwens R. J., et al., 2021, AJ, 162, 47 Bouwens R., Illingworth G., Oesch P., Stefanon M., Naidu R., van Leeuwen I., Magee D., 2023a, MNRAS, 523, 1009 Bouwens R. J., et al., 2023b, MNRAS, 523, 1036 Bowler R. A. A., Jarvis M. J., Dunlop J. S., McLure R. J., McLeod D. J., Adams N. J., Milvang-Jensen B., McCracken H. J., 2020, MNRAS, 493, 2059\n```  \nMandelker N., van Dokkum P. G., Brodie J. P., van den Bosch F. C., Ceverino'}