{"id":"a60ed243-b2ff-4c98-bfde-5b0afa95805b","arxiv_id":"2411.10305","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":9,"one_line_summary":"A GCM-coupled 2-moment microphysical cloud scheme produces long-lived patchy KCl clouds and 0.5-1% rotational variability in a WISE 0359-54 Y-dwarf model.","lead":"The authors build a fast two-moment cloud microphysics model, tracking particle number and mass density, and couple it to a 3D general circulation model of the Y-dwarf WISE 0359-54. The simulation produces patchy KCl clouds and 0.5-1% rotational brightness variations, which can be compared with JWST and Spitzer observations.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The monodisperse delta-function closure (Appendix A.1) is the load-bearing assumption: cloud dynamics, settling, collisions, and opacity feedback all use a single radius, while the synthetic spectra use a broad exponential distribution (Section 5/Appendix B).","rationale":"The paper is a clear, self-aware modeling study with standard moment equations and a public code, and the forward JWST comparison is a legitimate test rather than a fit. I agree with the reader that the monodisperse closure is the load-bearing assumption. The GCM's dynamical cloud feedback is computed with delta-function particles, but the spectra are post-processed with an exponential distribution, so the observable comparison does not validate the feedback assumption used inside the simulation. Appendix B explicitly shows that the choice of distribution shape matters for opacity, and the paper acknowledges this limitation. A sensitivity test with a fixed-width gamma closure is the direct way to decide whether the central suitability claim survives. The Kzz = 10^-5 vs 10^5 inconsistency in Section 4 should be corrected, but it is secondary and does not by itself overturn the method claim. If the gamma-closure test shows negligible changes in cloud mass, rc, and variability, the delta-function closure is adequate and the paper's claim stands; if it shows large changes, the claim must be qualified. Since the reader's conditional verdict already captures this risk, my read does not change the verdict.","tokens_in":24449,"tokens_out":6092,"duration_ms":60278,"concrete_test":"Run a controlled sensitivity test: in the same GCM framework, replace the delta-function closure with a gamma-distribution 2-moment closure of fixed width (e.g., alpha = 2, as in Ziegler 1985) while keeping all other parameters identical, and compare globally averaged cloud mass, rc, cloud optical depth, and Spitzer [3.6]/[4.5] variability. If these quantities differ by more than the model's own 0.5-1% variability signal (or by more than the JWST fitting residuals), the monodisperse closure is load-bearing and the central claim must be softened. A 1D column version of the same comparison would be a faster first step.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Central claim: the 2-moment scheme is 'highly suitable' for investigating cloud characteristics and feedback in GCMs. The scheme closes with f(m) = Nc delta(m - mc), so Nc and rho_c only evolve a single radius rc. Equations A.8 and A.9 evaluate the settling velocity, condensation rate, and collision kernel at mc, and Section 3.1 computes opacity from rc. Appendix B explicitly demonstrates that distributions with the same E[m] and V[m] (log-normal, inverse gamma, exponential, Rayleigh) yield materially different radius distributions and distribution-integrated opacities. The paper itself chooses an exponential distribution for the synthetic spectra in Section 5, so the observable comparison is not generated by the same size-distribution assumption that supplied the cloud opacity feedback inside the GCM. If the true KCl size distribution is broad rather than delta-like, the simulated cloud mass, vertical extent, and radiative feedback, and hence the 0.5-1% variability and the JWST comparison, could shift. This is not an internal contradiction, but the central 'suitability' claim is not established until the sensitivity to this closure is quantified.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a new two-moment bulk microphysical cloud scheme for sub-stellar atmospheres, derived from the mass moment equations of the particle size distribution. The scheme evolves the zeroth and first moments (number density and mass density of cloud particles) with source/sink terms for homogeneous nucleation, condensation/evaporation, collisional growth (Brownian coagulation and gravitational coalescence), seed-particle evaporation, and deep vapour replenishment. The scheme is coupled to the Exo-FMS GCM and applied to a WISE 0359-54 Y-dwarf parameter regime with KCl clouds. The authors report a sluggish atmospheric circulation, patchy long-lived cloud structures, synthetic spectra that reproduce the JWST observations shortward of about 7 microns only if additional cloud opacity is added, and 0.5-1% Spitzer-band variability. The central claim, stated in the Conclusions, is that the two-moment scheme is 'highly suitable' for investigating cloud characteristics and feedback in GCMs.","tokens_in":24664,"tokens_out":5051,"duration_ms":50439,"significance":"If the central claim is robust, the scheme would be a valuable contribution: it adds time-dependent cloud microphysics with collisions to 3D GCMs at moderate computational cost while retaining the main dynamical and radiative feedbacks. The derivation of the moment equations from the Smoluchowski equation (Section 2 and Appendix A) is clear and standard, the code is publicly available, and the JWST comparison is a genuine forward test in which the authors explicitly report the model's shortcomings rather than tuning to the data. The main caveat is that the predicted cloud structures, spectra, and variability rest on the monodisperse delta-function closure, and the paper itself shows in Appendix B that different closures give materially different opacities. This limitation is acknowledged but not quantified, so the central 'suitability' claim is not yet fully established.","major_comments":[{"comment":"There is a direct internal inconsistency in the quoted minimum vertical mixing coefficient: Section 3 states 'We assume a minimum Kzz = 10^5 cm2 s−1' while Section 4 states 'Vertical mixing is therefore dominated by the prescribed minimal Kzz = 10−5 cm2 s−1'. These differ by ten orders of magnitude and the correct value is load-bearing for the discussion of chemical quenching and vertical transport. The authors must correct the typo and, more importantly, verify that all conclusions (e.g., the NH3 fit in Section 5) correspond to the actual value used in the simulation.","section":"Sections 3 and 4"},{"comment":"The GCM cloud opacity feedback and all microphysical rates assume a monodisperse size distribution f(m) = Nc δ(m − mc) (Appendix A.1), but the synthetic spectra in Section 5 are post-processed with an exponential size distribution reconstructed from the moments (Appendix B). Because Appendix B demonstrates that distributions with identical moments (log-normal, inverse gamma, exponential, Rayleigh) produce substantially different radius distributions and distribution-integrated opacities, the observable comparison is not generated under the same closure that supplied the cloud opacity feedback inside the GCM. The authors should quantify the sensitivity of the simulated cloud structure, spectra, and variability to the assumed closure, for example by repeating the post-processing with the monodisperse distribution or by rerunning a short simulation with a gamma distribution with fixed shape parameter.","section":"Sections 3.1, 5, and Appendices A.1 and B"},{"comment":"The gravitational coalescence term is problematic under the monodisperse closure. With a delta-function size distribution all particles have the same radius and settling velocity, so the differential settling velocity in the kernel (Eq. 33) is identically zero. The model bypasses this by introducing the ad hoc parameter ϵ = 0.5 to estimate the relative velocity (Eq. 34), a value taken from protoplanetary disk grain-growth simulations (Sato et al. 2016). Since collisional growth is a new component of the scheme and directly affects particle sizes and cloud opacity, the sensitivity of the results to ϵ should be quantified, or the approximation should be justified for sub-stellar conditions.","section":"Section 2.1.2"},{"comment":"Equation (B.5) defines σm = exp(sqrt(ln(1 + V[m]^2/E[m]^2))). The argument V[m]^2/E[m]^2 is dimensionless only if the variance is dimensionless, but V[m] has units of g^2, so the expression is dimensionally inconsistent. The correct dimensionless ratio is V[m]/E[m]^2. As written, this formula would produce incorrect distribution widths in the reconstructed log-normal distribution, which feeds directly into the spectral post-processing in Section 5.","section":"Appendix B.1, Eq. (B.5)"}],"minor_comments":[{"comment":"There are numerous OCR-style typos and spacing errors in the manuscript text (e.g., 'o ffers', 'di fferent', 'with it’s', 'we not that'), which should be corrected in the final version.","section":"Throughout"},{"comment":"The text says 'This patchiness is long-lived in the GCM simulations, with the features in the simulations lasting several rotation rates.' This should likely read 'several rotation periods', and the concept of a 'rotation rate' as a duration is confusing.","section":"Section 4"},{"comment":"The phrase 'sub-rotational variability' is used to describe the light curves; it is unclear whether this means variability at timescales shorter than the rotation period, and the term should be defined.","section":"Section 5"},{"comment":"The thermal timescale estimate uses Teff = 548 K while the model setup uses Tint = 458 K; the relationship between these values and their definitions should be clarified to avoid apparent inconsistency.","section":"Section 6"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the journal's scope and the core derivation is sound. The main barrier to acceptance is the lack of a quantitative sensitivity analysis for the monodisperse closure and the internal inconsistency in the Kzz value. The closure issue is particularly important because the paper's own Appendix B shows that the choice of size distribution changes the opacity substantially, so the 'highly suitable' claim needs support. I would encourage the editor to request a revision that addresses these points rather than a rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper is worth a serious referee. What is new is not the moment equations — the authors openly say those follow Ohno & Okuzumi (2018) — but the coupling of a collision-inclusive 2-moment scheme into the Exo-FMS GCM and the first 3D cloud simulations for a Y-dwarf, including patchy KCl clouds and rotational variability of 0.5–1%. That is a tool the field will use, especially for JWST interpretation.\n\nThe paper does several things well. The derivation in Appendix A is clear and standard. The code is on GitHub. The authors are candid about limitations: the monodisperse closure, missing opacity below 7 μm, the need for extra cloud species. The JWST comparison is a forward test, not a fit, so there is no circularity. Self-citations to Lee (2023) and Ohno & Okuzumi (2018) are appropriate because the framework is explicitly built on those.\n\nThe soft spots are real but not fatal. The load-bearing assumption is the delta-function size distribution in Appendix A.1: settling, collisions, condensation, and opacity feedback all use a single radius. The paper itself shows in Appendix B that other distributions with the same moments give noticeably different radius distributions and opacities, and the synthetic spectra in Section 5 use an exponential distribution rather than the delta function used inside the GCM. So the observable comparison is not generated under the same closure that produced the cloud feedback. That means the central claim that the scheme is “highly suitable” is not fully established until the sensitivity to this closure is quantified. There is also an internal inconsistency in the quoted minimum Kzz: Section 3 says 10^5 cm2/s, Section 4 says 10^-5 cm2/s. A typo, probably, but it must be fixed. And τ_evap = 1 s is arbitrary and untested.\n\nThese issues do not undermine the basic contribution. The model is computationally cheap, physically motivated, and the authors have done an honest job reporting where it fails against data. The paper is aimed at brown dwarf and exoplanet GCM modelers and anyone trying to connect cloud microphysics to JWST spectra and variability. I would bring it to my reading group and would cite it if I worked in this area. It deserves peer review, with the expectation of a substantive revision: fix the Kzz typo, add a sensitivity test to the closure (even one 1D case), and discuss the distribution mismatch between the GCM and the post-processing.\n\nSend it out. The flaws are fixable and the tool is useful.","headline":"Useful and honest 3D GCM cloud microphysics paper; the delta-function closure is a real limitation, but the forward JWST test and available code make it worth refereeing.","tokens_in":25300,"tokens_out":2505,"would_cite":true,"duration_ms":24615,"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":"The paper demonstrates that a two-moment cloud microphysics scheme—just particle number plus condensed mass per condensate—can capture time-dependent cloud formation, patchiness, and spectra in a 3D model of a Y-dwarf atmosphere at…","keywords":["brown dwarfs","Y dwarfs","cloud microphysics","2-moment bulk scheme","general circulation models","atmospheric variability","KCl clouds","JWST spectroscopy"],"falsifier":"Run the identical WISE 0359-54 GCM setup with a bin-resolving scheme that tracks dozens of particle sizes (the CARMA-style approach of Gao et al. 2018) and compare cloud structure, spectra, and variability amplitudes; substantial differences would show that the single-size closure is the limiting assumption. On the observational side, a JWST spectrum of a Y-dwarf whose retrieved cloud particle size distribution is clearly broad rather than narrow around one mean radius would test that same assumption directly.","tokens_in":24087,"feed_emoji":"☁️","tokens_out":17942,"duration_ms":148239,"temperature":0.7,"pith_summary":"The paper argues that a bulk cloud model tracking only two moments of the particle size distribution—total particle number and total condensed mass—can carry time-dependent cloud microphysics into three-dimensional circulation models of brown dwarfs and exoplanets at acceptable cost. The authors derive moment equations for homogeneous nucleation, condensation and evaporation, and collisional growth, close the distribution with a single mean particle size, and couple the scheme to a GCM of the Y-dwarf WISE 0359-54 with KCl clouds. The simulation produces a global but patchy cloud field, an equatorial belt of larger particles, emission spectra that broadly match the JWST observations, and roughly 0.5-1% rotational variability in the Spitzer bands. If this demonstration holds, cloud microphysics stops being a bottleneck for 3D modelling and becomes a standard module for interpreting the JWST spectra of cold sub-stellar atmospheres.","feed_headline":"Two cloud numbers now carry cloud physics into 3D brown-dwarf models","feed_subtitle":"The 2-moment scheme matches Y-dwarf spectra beyond 7 µm and reproduces rotational variability at a cost 3D models can afford.","key_machinery":"The load-bearing object is the two-moment mass-moment system. The zeroth and first moments of the particle mass distribution, $M^{(0)} = N_c$ (number density) and $M^{(1)} = \\rho_c$ (mass density), are advanced by coupled ODEs that sum the rates of homogeneous nucleation (with seed mass $m_{\\rm seed}$), condensation and evaporation evaluated at the mean particle mass, Brownian coagulation and gravitational coalescence reduced to same-size collisions at the mean radius, seed-particle evaporation, and a deep vapour-replenishment relaxation. Closure is the monodisperse delta-function distribution $f(m) = N_c\\,\\delta(m - m_c)$, so every process—settling velocity, collision rates, and cloud opacity—is evaluated at the single mass-weighted mean radius $r_c = (3m_c/4\\pi\\rho_d)^{1/3}$. The practical payoff is that only two tracers per condensate are advected in the GCM, which is what makes time-dependent microphysics affordable inside 3D simulations.","core_discovery":"The central result, in the authors' own framing, is a capability demonstration: a time-dependent two-moment bulk microphysical scheme—evolving the zeroth and first moments of the particle mass distribution, $N_c$ and $\\rho_c$, together with the condensable vapour density—can be coupled directly to a 3D general circulation model and still reproduce the main observable consequences of clouds in a sub-stellar atmosphere. The scheme unifies homogeneous nucleation, condensation and evaporation, Brownian coagulation, and gravitational coalescence in one moment framework, closes the size distribution with a monodisperse delta function at the mass-weighted mean radius, and is integrated on the same timestep as the GCM's chemistry. Applied to the Y-dwarf WISE 0359-54, it yields a global KCl cloud layer with long-lived patchiness at higher latitudes, an equatorial belt of larger particles, synthetic spectra consistent with the JWST data beyond about 7 $\\mu$m, and 0.5-1% periodic and sub-rotational variability in the Spitzer [3.6] and [4.5] bands. The authors conclude that the 2-moment scheme is a highly suitable method for investigating cloud characteristics and feedback in GCMs and other large-scale simulations.","pith_inferences":["The paper's own Appendix B shows that log-normal, inverse-gamma, exponential, and Rayleigh distributions with the same moments give noticeably different radius distributions and opacities; a consequence the authors leave implicit is that the specific patchiness patterns and variability amplitudes should be treated as contingent on the monodisperse closure rather than as predictions that would surv","The missing opacity below about 7 $\\mu$m and the low variability amplitude compared with objects like WISE 1405 (Cushing et al. 2016) point to a concrete next test: letting H2O condense on KCl seeds, which the paper flags as future work, could add the high-altitude opacity that the current single-species scheme cannot provide.","Re-closing the same two moments with a gamma or log-normal distribution instead of a delta function—the route Earth cloud models took with fixed shape parameters—would cost almost nothing and would let JWST retrievals of size distribution width feed directly back into the GCM.","Because the equatorial belt of larger particles is sustained by the deep vapour-replenishment boundary condition, a parameter sweep over the deep mixing strength ($K_{zz}$ floor) or the deep KCl reservoir abundance would quantify how much of the predicted 0.5-1% variability is a boundary-condition choice rather than atmospheric dynamics."],"forward_implications":["Time-dependent microphysics—nucleation, condensation, and collisional growth—can be run inside a 3D GCM for the cost of two advected tracers per condensate, avoiding the 40-70 bins that bin schemes require.","For the Y-dwarf regime, KCl clouds form a global layer that is patchy at higher latitudes and shows longitudinally varying particle sizes at the equator, with the patchy features lasting over many rotations.","The synthetic spectra match the JWST spectrum of WISE 0359-54 longward of about 7 $\\mu$m but over-predict flux at shorter wavelengths, so additional cloud species (Na2S, ZnS, or NaCl) or mixed-composition grains are needed to add opacity.","The model yields roughly 0.5-1% periodic and sub-rotational variability in the Spitzer [3.6] and [4.5] bands with the right rotational periodicity, but a lower amplitude than the variability observed on other Y dwarfs.","Cloud opacity feedback in this regime is enough to help force weak equatorial jets even in an otherwise sluggish, stagnant atmosphere, tying cloud radiative forcing to the dynamical regime."],"supporting_citations":[{"why":"Supplies the general method, cited at the outset, for deriving the time derivatives of the mass moments that the whole 2-moment scheme is built on.","marker":"Gail & Sedlmayr (2013)"},{"why":"Provides the mass-moment generative function of the Smoluchowski equation from which the collisional gain-and-loss rates for $N_c$ are derived.","marker":"Drake (1972)"},{"why":"The 2-moment framework and the settling-velocity expression that the paper says its master equations are equivalent to and that it extends with nucleation and collisions.","marker":"Ohno & Okuzumi (2018)"},{"why":"The mini-cloud code package and cloud-opacity calculation scheme that this work expands into a 2-moment, collision-inclusive GCM-coupled model.","marker":"Lee (2023)"},{"why":"Supplies the condensation and evaporation growth rates in both the diffusion-limited and free-molecular regimes.","marker":"Woitke & Helling (2003)"},{"why":"Gives the modified classical nucleation theory rate used for homogeneous KCl seed formation.","marker":"Helling & Fomins (2013)"},{"why":"Supplies the interpolation-function coagulation kernel actually used in the scheme, and contrasts with the more expensive bin-resolving approach the paper compares bulk schemes against.","marker":"Gao et al. (2018)"},{"why":"Supplies the deep vapour-replenishment relaxation formula, and its saturation-adjustment cloud model is the baseline replaced in the GCM test.","marker":"Komacek et al. (2022)"},{"why":"The JWST spectrum of WISE 0359-54 against which the GCM emission spectra are compared.","marker":"Beiler et al. (2023)"},{"why":"Supplies the adopted WISE 0359-54 parameters (internal temperature, metallicity, gravity) and the retrieved quench abundances used to validate the GCM chemistry.","marker":"Kothari et al. (2024)"}],"fun_headline_variants":["Two-moment cloud microphysics now runs inside 3D Y-dwarf models","Y-dwarf cloud patchiness and variability from two-moment microphysics","3D cloud microphysics matches Y-dwarf spectra beyond 7 microns","Two cloud moments now describe brown-dwarf cloud evolution in 3D","Cloud microphysics at GCM cost: two moments capture Y-dwarf clouds"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The entire calculation is anchored to a single particle size: the model assumes every cloud at a given location has the same radius, and it computes all growth rates, settling speeds, and the amount of light blocked using that one radius, so if real Y-dwarf clouds contain a wide spread of particle sizes, the simulated cloud patterns, spectra, and variability would change.","fun_headline_variants_meta":{"raw":{"variants":["Two-moment cloud microphysics now runs inside 3D Y-dwarf models","Y-dwarf cloud patchiness and variability from two-moment microphysics","3D cloud microphysics matches Y-dwarf spectra beyond 7 microns","Two cloud moments now describe brown-dwarf cloud evolution in 3D","Cloud microphysics at GCM cost: two moments capture Y-dwarf clouds"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00074,"raw_usage":{"total_tokens":3407,"prompt_tokens":1153,"completion_tokens":2254,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":769,"completion_tokens_details":{"reasoning_tokens":2155}},"tokens_in":769,"tokens_out":2254,"duration_ms":14835,"temperature":1.0,"reasoning_tokens":2155,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T19:45:49.047323+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the identical WISE 0359-54 GCM setup with a bin-resolving scheme that tracks dozens of particle sizes (the CARMA-style approach of Gao et al. 2018) and compare cloud structure, spectra, and variability amplitudes; substantial differences would show that the single-size closure is the limiting assumption. On the observational side, a JWST spectrum of a Y-dwarf whose retrieved cloud particle size distribution is clearly broad rather than narrow around one mean radius would test that same assumption directly.","supporting_citations":[{"cited_title":"& Okuzumi , S","cited_arxiv_id":null,"evidence_quote":"The 2-moment framework and the settling-velocity expression that the paper says its master equations are equivalent to and that it extends with nucleation and collisions."},{"cited_title":"C., Burningham , B., et al","cited_arxiv_id":null,"evidence_quote":"Supplies the adopted WISE 0359-54 parameters (internal temperature, metallicity, gravity) and the retrieved quench abundances used to validate the GCM chemistry."}],"review_version":1}