Imagine you hire ten interns to rewrite your company's Excel macros in Python. Some produce clean, runnable code on the first try. Others need you standing over their shoulder. A couple surprise you with genuinely clever refactors. But here's the kicker: when you go back and just retune the original Excel formulas with fresh data, they work almost as well — and they're faster, cheaper, and you actually understand what they're doing. That's the core result of this paper, applied to the physics equations inside Earth system models. The committed claim: LLM agents can synthesize scientific literature, reimplement published bulk parameterizations for warm-rain collision-coalescence as differentiable JAX code, and — in rare cases — propose new functional forms that slightly outperform retuned classical schemes. Specifically, 2 out of 100 unattended agent-generated closures achieved lower rain-mass error than the best retuned physical scheme on held-out superdroplet simulations. That's the headline number, and it's honest about both the promise and the strike rate. The ladder here is revealing. The paper compares three families: conventional bulk parameterizations (the textbook equations climate models actually use), data-driven emulators (neural nets trained directly on simulation data), and agent-generated closures (LLM-written code). The reference standard is high-fidelity superdroplet simulations — essentially a particle-resolved ground truth for raindrop formation. The surprise is that retuned physical schemes — classical equations with fresh coefficients — matched or outperformed many learned parameterizations at substantially lower computational cost. The agents' advantage is narrow: only 2 of 100 unattended runs beat the best retuned scheme, and only on rain-mass error. Architecturally, this is an agentic-AI-for-scientific-code pipeline. Ten language models (not named individually in the abstract, but spanning the frontier) are given the task of reading literature, reimplementing known schemes, and proposing new differentiable closures in JAX. The key structural choice is differentiability and mass conservation as hard constraints — the code must be autodiff-compatible and physically consistent. Iterative feedback (expert-in-the-loop or automated rollout error signals) lowered median rollout error for 8 of 10 models. Unattended performance varied strongly by scheme and model, which tells you the pipeline is brittle without guidance. The integrity regime is solid for a workshop paper. The ground truth is superdroplet simulations — a well-established community reference for cloud microphysics. Held-out simulations provide a real test set. The comparison spans conventional, learned, and agent-generated approaches on the same benchmark. The 100-run unattended experiment across ten models is a genuine attempt at statistical coverage, not a cherry-picked best run. The main weakness: no independent replication, no pre-registration, and the warm-rain test case is deliberately simple — collision-coalescence is one of the easier microphysics problems. The milestone question is where this gets interesting. The paper's implicit trajectory is: if agents can write valid, mass-conserving, differentiable parameterizations for warm rain, can they handle ice microphysics, convection, or turbulence closures? Those are the parameterizations that actually dominate climate model uncertainty. The gap between warm-rain collision-coalescence and, say, deep convective parameterization is enormous — not just in physics complexity but in the number of interacting variables and the sensitivity of the coupled system. A concrete next milestone: agent-generated ice-phase microphysics closures that run stably inside a single-column model for 1,000 simulated days. The obvious experiment not run: coupling an agent-generated closure into an actual ESM column or global run and measuring whether it stays stable and physically consistent over long integrations. The honest read is (a) — compute and engineering scope. Plugging a new parameterization into an ESM is a months-long integration task, not something you do for a workshop paper. But it's the experiment that would convert this from a proof of concept to a result that climate modelers would actually use.
