{"id":"90a8c97c-796a-482d-bb5e-fa7c75e312f2","arxiv_id":"2509.19427","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Early dark energy in large hydrodynamic simulations reproduces JWST's excess of bright, massive, high-redshift galaxies while preserving low-redshift Lambda-CDM behavior.","lead":"This paper uses the first large cosmological simulations of an early-dark-energy universe to show that it naturally produces more bright, massive galaxies at early times, matching JWST observations without retuning the galaxy-formation model. The same model converges back to standard cosmology by redshift 3, so it claims to preserve everything Lambda-CDM already gets right.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The EDE-vs-ΛCDM conclusion is only as secure as the IllustrisTNG star-formation efficiency at z=4–14; Fig. 4 shows alternative high-z SFE models with 0.5–1 dex higher values that could mimic EDE.","rationale":"The reader's weakest assumption identifies the same load-bearing issue: the TNG subgrid model, calibrated at low redshift, is assumed to be accurate at z=4–14 without recalibration. My stress-test sharpens this into a concrete degeneracy: EDE primarily boosts the halo mass function, while alternative high-z star-formation models can boost the galaxy-halo connection by a similar amount. The paper's own Fig. 4 displays this range of SFE, so the concern is grounded in the manuscript's evidence, not in outside speculation. I am not claiming the model is wrong; I am claiming the central inference is conditional on the SFE assumption being correct. The proposed empirical-model test would settle whether the EDE signal is distinguishable from a baryonic SFE enhancement. Because the reader already assigned CONDITIONAL and my concern reinforces that condition rather than overturning the paper, I recommend UNCHANGED.","tokens_in":32078,"tokens_out":6538,"duration_ms":48794,"concrete_test":"Use the paper's own empirical model (§2.4) to compute the z=4–14 UV LF and z=3–9 SMF in ΛCDM with a mass- and redshift-dependent SFE boost equal to the EDE/ΛCDM halo mass function ratio (from Fig. A1 or the simulation subhalo mass functions). If the boosted-ΛCDM curves agree with the EDE curves within ~0.1–0.2 dex, then the EDE advantage is fully degenerate with a high-z SFE increase and the 'no recalibration' claim is unsupported. If the boosted curves still fall below EDE (e.g., because of age or morphological differences), the degeneracy is broken and the EDE interpretation is strengthened.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's headline inference—that EDE is required to match JWST—depends on the IllustrisTNG subgrid model's high-redshift star-formation efficiency being correct. Both runs share the same TNG physics (§2.2), and Fig. 4 shows that TNG's SFE is nearly redshift-independent and sits at the lower envelope of published high-z models. Alternative prescriptions (FIREbox at z=6, thesan-zoom at z=6, FFB at z=10) give SFE up to ~0.5–1 dex higher at M_halo ~ 10^10–10^11 M_sun (Fig. 4; Li+23; Feldmann+25; Shen+25). Since EDE's effect is to raise the halo mass function by ~0.2–0.8 dex over the relevant range (§3.1, Fig. A1), a modest upward revision of TNG's high-z SFE in ΛCDM could produce the same UV LF and SMF as EDE+TNG. The 'essentially no additional calibration' claim is therefore relative to TNG, not to the broader space of plausible baryonic models. The argument is not internally inconsistent, but the discriminating power of the JWST match is not quantified. The quenched-galaxy underprediction (§3.4) is an acknowledged incompleteness; the more central issue is that the success channel itself is degenerate with high-z star-formation efficiency.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript presents the first large-scale cosmological hydrodynamic simulations of an Early Dark Energy (EDE) cosmology using the IllustrisTNG galaxy formation model, alongside a matched ΛCDM simulation. EDE parameters are taken from a CMB fit (Smith et al. 2022) and are not adjusted to galaxy data. The central claim is that EDE, through an enhanced small-scale matter power spectrum, boosts the abundance of UV-bright and massive galaxies at z ≃ 4–14, bringing simulated UV luminosity functions and stellar mass functions into 'excellent agreement' with JWST measurements, while predictions converge to ΛCDM by z ≲ 3. The paper also reports earlier emergence of disky galaxies in EDE and a surprising delay in quenched galaxies due to the younger universe. The comparison is supported by an empirical model and by an appendix isolating the power-spectrum change from the H0 change.","tokens_in":32345,"tokens_out":5256,"duration_ms":44263,"significance":"If the headline claim holds, this is an important result: it demonstrates that modifications to pre-recombination cosmology can simultaneously ease several JWST 'too early, too massive, too disk-like' tensions without retuning the IllustrisTNG baryonic model. The controlled EDE-versus-ΛCDM comparison with identical subgrid physics is a genuine strength, as is the authors' use of an independent empirical model and the clear decomposition of the effect into power-spectrum versus H0 changes. The paper is also commendably transparent about the quenched-galaxy underprediction. However, the absolute 'excellent agreement with JWST' is less secure than the differential EDE-ΛCDM signal: it depends on post-processing choices (zero dust attenuation at z ≥ 10, empirically calibrated dust at z < 10) and on the IllustrisTNG star-formation efficiency at high redshift, which Fig. 4 shows is at the lower envelope of current models. These degeneracies do not invalidate the differential statement, but they weaken the inference that EDE is required by the JWST abundance data.","major_comments":[{"comment":"The claimed 'excellent agreement' at z ≥ 10 relies on setting dust attenuation to zero at z ≥ 10 while using empirically calibrated dust at lower redshifts. This is an upper-limit choice in brightness: simulated galaxies do contain dust, and the motivation 'observed blue UV slopes' is partly circular when the goal is to match observed UV luminosity functions. The EDE-ΛCDM differential remains robust, but the absolute normalization of the z = 12–14 UV LFs, which is a headline result, is not. Please show the intrinsic A_UV distribution of the simulated galaxies at z ≥ 10, or test the sensitivity of Fig. 5 to a small nonzero attenuation (e.g., A_UV = 0.2–0.5 mag).","section":"§3.1, Fig. 5; §2.3–2.4"},{"comment":"The inference that EDE is needed to match JWST is degenerate with the high-redshift star-formation efficiency. The IllustrisTNG SFE in Fig. 4 is nearly redshift-independent and lies at the lower envelope of published models; FIREbox, thesan-zoom, and the FFB scenario show SFEs 0.5–1 dex higher at M_halo ~ 10^10–10^11 M_sun. Because the EDE boost to the halo mass function is roughly 0.2–0.8 dex over the relevant range, a modest upward revision of TNG's high-z SFE in ΛCDM could produce the same UV LF and SMF as EDE+TNG. The 'essentially no additional calibration' claim is therefore relative to TNG, not to the space of plausible baryonic models. Please quantify this degeneracy—e.g., show what a ΛCDM run with a 0.3–0.5 dex higher SFE would predict—or identify an observable that breaks it (clustering, stellar mass functions at fixed UV luminosity, or [C II] kinematics).","section":"§2.2, Fig. 4; §3.1"},{"comment":"The conclusions at z = 9–14 are based on a single (100 cMpc)^3 simulation volume per cosmology. In Fig. 5 the shaded 'one galaxy per mag' region indicates that the bright-end bins at z = 12 and 14 contain very few galaxies, and the Poisson errors do not include cosmic variance. The empirical model and Eq. (4) provide a useful volume correction, but the claimed 'almost perfect' agreement at z ≳ 12 could be affected substantially by sample variance. Please quantify cosmic variance—for example, by using the empirical model to estimate the expected field-to-field scatter in the z = 12–14 UV LF bins, or by running a second realization of at least one cosmology.","section":"§3.1, Fig. 5; Table 2"},{"comment":"The abstract's phrase 'simultaneously reconcile multiple high-redshift challenges' is too broad in light of the paper's own finding that quenched galaxies at M* ≳ 10^10 M_sun are underpredicted by an order of magnitude at z ≳ 3 in both cosmologies, and that EDE delays rather than accelerates quenching. The authors are appropriately candid in §3.4 and §4, but the abstract and conclusions should explicitly state that EDE does not address the early-quenching tension and may worsen it. This does not undermine the UV LF and SMF results, but it does limit the 'multiple challenges' claim.","section":"Abstract; §3.4, Fig. 10"}],"minor_comments":[{"comment":"The empirical-model parameter σ_UV = 0.75 mag is described as 'assumed' and 'purely for interpretation.' Since this scatter directly affects the LF shape and the stellar-mass-function derivation in Eq. (1), please state its source and the sensitivity of the empirical-model curves to it.","section":"§2.4"},{"comment":"The adopted EDE parameters from Smith et al. (2022) are now less favored by Planck PR4+BAO and are more aggressive than the ACT DR6/DESI-based fits cited in §2.1. The paper discusses this, but the abstract and conclusions should more clearly present the simulations as a proof-of-principle for a broader class of early-universe models, not as evidence for the specific best-fit model.","section":"§2.1, Table 1"},{"comment":"The disk thresholds (D/T)_* = 0.7 and (D/T)_gas = 0.8 are reasonable but somewhat arbitrary. A brief sensitivity test (e.g., varying the thresholds by ±0.1) would help confirm that the half-dex EDE-ΛCDM difference in disky-galaxy number density is not driven by the exact cut.","section":"§3.3"},{"comment":"The Schechter-function fits at z ≥ 5 fix the break mass to 10^11 M_sun. This is justified, but the choice could affect the integrated stellar-mass-density comparison; a short statement of the resulting systematic uncertainty would be useful.","section":"§3.2, Fig. 7"}],"recommendation":"major_revision","confidential_remarks":"This is a carefully executed simulation study with a robust differential EDE-ΛCDM comparison. My main concern is that the paper's headline claim moves from 'EDE can reconcile these observations under TNG baryonic physics' to 'EDE is required', and the current text does not quantify the SFE degeneracy or the no-dust assumption. These are fixable with additional analysis or a tempered interpretation, so I recommend major revision rather than rejection. The paper is a good fit for MNRAS if the authors address the absolute-comparison caveats."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThis is the first paper to put Early Dark Energy into large-volume cosmological hydrodynamic simulations, and that alone matters. The EDE-vs-LCDM comparison is clean: identical IllustrisTNG subgrid physics, both runs evolved to z=3, same box and resolution. The direction of the effect is robust — EDE boosts small-scale power and yields more UV-bright, massive, disky galaxies at z>6. The disk prediction is genuinely new and testable: about half a dex more gas disks at z~7, which ALMA should be able to check.\n\nWhere the paper gets soft is the absolute claim of “excellent agreement” with JWST. That agreement depends on post-processing choices and on the TNG star-formation efficiency. Figure 4 shows TNG’s SFE sits at the lower envelope of high-z models; FIREbox and thesan-zoom are 0.5–1 dex higher at the relevant halo masses. Since EDE shifts the halo mass function by a comparable amount, a plausible upward revision of TNG’s high-z SFE in LCDM could produce the same UV luminosity functions and stellar mass functions. The paper acknowledges the SFE model spread but does not quantify the degeneracy, and the “no additional calibration” framing can read as though the comparison is to all plausible baryonic models rather than to TNG itself. That is the central weakness. It is not fatal: the EDE-vs-LCDM difference is real and externally anchored by CMB-fixed parameters. But the JWST match alone does not single out EDE.\n\nThe dust treatment is also a bit circular at z>=10: the no-attenuation choice is justified by the observed blue UV slopes, which are the same data being matched. Minor, but worth stating. Also, one 100 Mpc box per cosmology means there is no cosmic variance error on the EDE-vs-LCDM ratio, and there is a residual ~1-sigma tension with Weibel et al. at z=9.\n\nThe quenched-galaxy underprediction is bigger than the paper’s tone suggests: both runs miss the observed abundance by an order of magnitude at z>=3. To their credit, the authors state this explicitly and tie it to SMBH feedback limitations, but it does mean the TNG model is not validated at high redshift in at least one important channel.\n\nOverall, the paper is honest, clearly written, and the EDE boost itself is secure. The overstatement is in the word “excellent” applied to the absolute match. The paper should go to peer review, with a request for a quantitative sensitivity analysis to SFE, dust, and cosmic variance.\n\nEnd of letter.","headline":"First hydro sims of EDE show a robust EDE-vs-LCDM offset in high-z galaxies, but the absolute match to JWST is degenerate with star-formation efficiency and dust choices.","tokens_in":32943,"tokens_out":2140,"would_cite":true,"duration_ms":20629,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Early dark energy, tuned only by CMB data, reproduces JWST's surplus of bright, massive, disky galaxies at cosmic dawn, then fades back to standard cosmology by z≈3.","keywords":["early dark energy","high-redshift galaxies","JWST","UV luminosity function","stellar mass function","disk galaxies","cosmological hydrodynamic simulations","Hubble tension"],"falsifier":"Measure the UV luminosity function at z≈12–14 with a larger, purely spectroscopic sample: if the bright-end number density turns out to be roughly 0.5 dex lower than current JWST photometric estimates, the EDE match disappears. Conversely, a robust measurement of the clustering bias of bright galaxies at z≈6–8—EDE predicts lower bias than ΛCDM at fixed luminosity—would discriminate the two.","tokens_in":31877,"feed_emoji":"🔭","tokens_out":4312,"duration_ms":685009,"temperature":0.7,"pith_summary":"The paper argues that a pre-recombination burst of early dark energy—an ingredient originally proposed to resolve the Hubble tension—also explains the surprisingly abundant, massive, star-forming disk galaxies JWST sees at z≈4–14. It does this with two large cosmological hydrodynamic simulations, one in standard ΛCDM and one with EDE parameters fixed by CMB fits; the galaxy formation model is the same out-of-the-box, low-redshift-calibrated model in both runs. In the EDE run, UV luminosity functions and stellar mass functions agree with JWST measurements with essentially no recalibration, and the number density of disky galaxies rises by about half a dex at z≈6–7. Predictions converge to ΛCDM at z≲3, preserving known low-redshift successes. The paper concludes that early-universe physics can simultaneously address several JWST anomalies and the Hubble tension.","feed_headline":"Early dark energy explains JWST's surplus of bright galaxies","feed_subtitle":"CMB-fixed early dark energy makes simulated galaxies at z=4–14 match JWST counts, then converges to ΛCDM by z≈3.","key_machinery":"The central object is early dark energy as a scalar field with an axion-like potential V(φ)≈[1−cos(φ/f)]³, with parameters taken from CMB fits that also yield H0≈74.8 km/s/Mpc. The field acts as a cosmological constant before recombination and then dilutes; its observable imprint is an enhanced small-scale linear power spectrum (higher n_s and A_s) that accelerates the collapse of the first halos. The simulations couple this cosmology to the IllustrisTNG galaxy formation model in the Arepo code, with identical subgrid physics in EDE and ΛCDM runs, so any galaxy differences come purely from the altered expansion history and power spectrum.","core_discovery":"The central claim is that EDE's accelerated early structure formation, driven mainly by an enhanced small-scale matter power spectrum with higher spectral index, primordial amplitude, and matter density, boosts the abundance of bright and massive galaxies enough to remove the 0.5–1 dex shortfall that standard ΛCDM simulations show against JWST at z≳10. Because EDE decays after recombination, the same run converges back to ΛCDM predictions at z≲3 with ≲0.2 dex differences. The faster assembly of massive halos also makes stellar and gaseous disks appear earlier: number densities of disky galaxies are about half a dex larger at z≈6–7, while the disky fraction at fixed stellar mass is unchanged,","pith_inferences":["Any beyond-ΛCDM model that increases the pre-recombination expansion rate generically tilts the small-scale power spectrum upward, so high-redshift galaxy abundance could serve as a broad probe of early-universe physics beyond EDE itself.","The paper's reduced-bias prediction—bright galaxies less clustered in EDE at fixed luminosity—offers a direct observational test: galaxy clustering at z≈6–10 from JWST fields could distinguish EDE from baryonic solutions like feedback-free starbursts or top-heavy initial mass functions.","If high-z star formation or feedback is later shown to be systematically different, the EDE match could be coincidental; the paper's own shortfall in quenched galaxies already signals that the subgrid black-hole model limits the comparison.","Lyman-α forest and 21-cm observations, usually used to rule out suppressed small-scale power, could be turned around to constrain enhanced small-scale power of the kind EDE produces."],"forward_implications":["JWST's excess of UV-bright galaxies at z≈10–14 does not require exotic star-formation physics; a CMB-consistent EDE cosmology with standard galaxy formation reproduces the observed counts.","Massive galaxy candidates in the COSMOS-Web field at z≈11–12, which appear as >3σ outliers in ΛCDM, fall within about 2–3σ of EDE predictions.","Stellar and gaseous disks form earlier in EDE, so ALMA and JWST disk detections at z≈6–8 become expected rather than surprising, with gaseous disk number densities roughly a factor of three higher at z≈6–7.","At z≲3 the EDE run converges to ΛCDM within ≲0.2 dex, preserving the established low-redshift successes of the standard model.","The model predicts fewer massive quenched galaxies than ΛCDM at z≳3.5, supporting the paper's view that the early-quenching puzzle is likely a baryonic or black-hole-physics limitation, not a cosmological one."],"fun_headline_variants":["Early dark energy fixes JWST's bright galaxy excess","EDE matches JWST counts, then yields to ΛCDM","EDE speeds up early galaxy and disk formation","Early dark energy explains JWST's luminous galaxies","EDE boosts early galaxies and disks to JWST levels"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The galaxy formation model, calibrated at low redshift, is assumed to stay accurate at z=4–14 with no recalibration; if real high-redshift star formation or feedback differs systematically, the EDE simulation's agreement with JWST could be coincidence rather than evidence for EDE.","fun_headline_variants_meta":{"raw":{"variants":["Early dark energy fixes JWST's bright galaxy excess","EDE matches JWST counts, then yields to ΛCDM","EDE speeds up early galaxy and disk formation","Early dark energy explains JWST's luminous galaxies","EDE boosts early galaxies and disks to JWST levels"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000855,"raw_usage":{"total_tokens":3587,"prompt_tokens":817,"completion_tokens":2770,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":561,"completion_tokens_details":{"reasoning_tokens":2692}},"tokens_in":561,"tokens_out":2770,"duration_ms":19389,"temperature":1.0,"reasoning_tokens":2692,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T15:21:53.805625+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the UV luminosity function at z≈12–14 with a larger, purely spectroscopic sample: if the bright-end number density turns out to be roughly 0.5 dex lower than current JWST photometric estimates, the EDE match disappears. Conversely, a robust measurement of the clustering bias of bright galaxies at z≈6–8—EDE predicts lower bias than ΛCDM at fixed luminosity—would discriminate the two.","supporting_citations":[],"review_version":1}