You're a doctor estimating a patient's age from a photograph. You have a set of reference photos — people at known ages, all of whom aged gracefully with consistent diet and exercise. Your patient, though, spent a decade partying, then five years fasting, then bulked up recently. The face shows mostly the recent bulk. Your reference set has no template for that life history, so you match to the closest smooth-aging reference and call it. You'll get the current weight roughly right, but you'll overestimate their age (the weathering looks older) and underestimate how fast they're currently gaining weight (the recent burst gets averaged into a longer, slower trend). That's exactly what SED fitting does to Lyman break galaxies — and this paper quantifies how wrong it gets. The committed claim: single-component SED fitting with smoothly declining star formation histories produces large systematic biases when applied to high-redshift Lyman break galaxies — specifically, underestimated star formation rates and overestimated ages — because real galaxies have bursty, complex formation histories that smooth templates can't represent. Stellar masses are recovered reasonably well when redshift is known, but with substantial scatter. The category of result is a systematic calibration study: not a new method, but a quantitative audit of the standard method's failure modes. This is an Nth iteration of a known concern, but at exceptional depth — 85 pages, 34 figures, submitted to ApJS. The ladder: the baseline being tested IS the standard practice. SED fitting using stellar population synthesis templates with smoothly declining star formation rates (τ-models) is the workhorse method used by essentially every extragalactic survey to derive physical parameters of distant galaxies. Lee, Idzi, Ferguson, Somerville, Wiklind, and Giavalisco don't propose a replacement — they build a controlled experiment to measure how badly the standard method fails. They construct mock galaxy catalogues using ΛCDM hierarchical structure formation theory plus semi-analytic baryonic physics (the Somerville models), select Lyman break galaxies at z ~ 3.4, 4.0, and 5.0, generate broadband photometry, then run the standard SED-fitting pipeline on that photometry and compare recovered parameters to the known truth. This is the gold standard for method validation: you know the answer, you run the tool, you measure the error. The architecture is simulation-as-testbed. The 'ground truth' comes from semi-analytic models (SAMs) of galaxy formation embedded in ΛCDM N-body merger trees — these produce galaxies with complex, physically motivated star formation histories including mergers, bursts, and quenching episodes. The 'fitting method' is standard χ² SED fitting using Bruzual & Charlot-style stellar population synthesis templates with exponentially declining SFRs. The mismatch is structural: SAM galaxies have multi-component, bursty histories; the templates assume single-component, smooth decline. The key computational property being exploited is that SAMs are cheap enough to produce statistically meaningful mock catalogues while retaining physical realism. Integrity is the strongest feature of this paper. The validation is about as non-circular as you can get in this field: the ground truth (SAM) and the fitting method (SED templates) use different star formation histories by construction, so the test genuinely measures template mismatch rather than internal consistency. The weakness is that SAMs themselves are models — if the SAM star formation histories are wrong, the bias estimates are wrong. But the direction of the result is robust: real galaxies almost certainly have MORE complex histories than SAMs, not less, so the biases identified here are likely lower bounds. The 85-page length with 34 figures suggests exhaustive exploration rather than cherry-picking — this is a reference document, not a letter. The specific findings give you calibration numbers to carry around. Stellar masses: recovered reasonably well when redshift is known, with minor systematic underestimation and substantial scatter. Star formation rates: systematically underestimated by single-component fits. Ages: systematically overestimated. The mechanism is clean: young stars from the most recent burst dominate the broadband light, hiding older stellar populations underneath. The smooth-decline template interprets this composite light as a single, older population forming at a moderate rate — when the truth is a younger population forming rapidly on top of older stars it can't see. This is the 'outshining' problem, and this paper gives it quantitative teeth across three redshift bins. The obvious successor experiment: running the same calibration with multi-component SED fitting (two-burst models, non-parametric star formation histories, or Bayesian SED codes like BEAGLE or Prospector that allow more flexible SFHs) to quantify how much of the bias disappears when the template library better matches real formation histories. My read: (c) saving it for the next paper. The 85-page scope of this paper is already enormous; the natural follow-up is 'here's how to fix it,' and multi-component fitting is the obvious fix. The Somerville group and the GOODS/CANDELS collaboration (Ferguson, Giavalisco) have the infrastructure to do this.
