{"id":"0f870782-bf77-4705-a57f-9bf1e9661431","arxiv_id":"2607.14845","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"A new dimensionless ratio, ΨHALD, predicts the maximum haze particle radius that can reach a hot Jupiter's morning limb, and matches 3D climate simulations to within a factor of a few.","lead":"This paper develops a simple analytical formula that estimates the largest haze particles that can be carried from the evening to the morning limb of a hot Jupiter before falling out of the atmosphere. It gives observers a quick way to decide whether morning-evening haze asymmetries seen by JWST are consistent with particle transport and settling, before running expensive 3D climate models.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"ΨHALD=1 is a transport-survival threshold, not a limb-asymmetry threshold; the claimed morning-limb enhancement requires unmodeled gyre trapping, leaving the central claim conditional.","rationale":"The reader's weakest assumption — unmodeled gyre trapping — is also the single most load-bearing point in my reading. The 1D model (Eq. 26) is explicit: without trapping, χ_m/χ_e = exp(-1/ΨHALD) < 1, so ΨHALD=1 corresponds to a morning limb ~37% of the evening limb, not to a 'higher or comparable' morning concentration. The abstract and Table 2 connect ΨHALD=1 to the maximum radius for a morning-limb enhancement, but that connection is not a consequence of the transport equations; it is an imported expectation from prior 3D simulations (Steinrueck et al. 2021; M25; Lee et al. 2026). The framework therefore predicts a necessary condition (particles can reach the morning limb) but not a sufficient condition (they subsequently accumulate). The GCM validation gives order-of-magnitude support for the tested cases, and the paper is careful to call the result an estimate; hence the conditional verdict is right. A secondary concern is the arbitrary evaluation pressure (0.01 mbar) and its unproven mapping to the transmission photosphere, but the trapping issue is more fundamental because it changes the logical status of the headline claim. The proposed test uses existing GCM outputs to check whether the transport-only criterion actually coincides with the observed limb-asymmetry transition; if it does not, the discrepancy quantifies the missing trapping term.","tokens_in":22429,"tokens_out":13587,"duration_ms":113060,"concrete_test":"From the existing g=20 m/s^2 WASP-39b-like GCM runs (Sec. 3.1, radii 50, 60, 70, 100 nm), compute the paper's ΨHALD (Eq. 20) at p=0.005 mbar for each radius using the same parameter values as Fig. 3. Determine the GCM-diagnosed transition radius (where morning/evening concentration ratio crosses 1). If the transition occurs at ΨHALD≠1 by more than, say, a factor 2, the ΨHALD=1 contour is not an accurate proxy for the limb-asymmetry threshold, and the missing trapping/growth processes must be quantified before the framework can be used to predict JWST observables.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 3.1 defines the ΨHALD=1 contour as 'the estimated maximum radius for which haze particles can survive the hemispheric transport and reach the morning limb before vertical removal dominates' and then equates this with the maximum radius 'capable of producing a higher or comparable haze concentration over the morning limb' (Fig. 3 caption, Table 2). The framework's own 1D kinematic model, Eq. (26), gives χ_m/χ_e = exp(-1/ΨHALD), which is <1 for every finite ΨHALD; at ΨHALD=1 the morning limb is depleted by a factor e relative to the evening limb. The jump from 'survive transport' to 'higher or comparable concentration' is made entirely by the assumption that particles reaching the morning limb are trapped in nightside gyres (Secs. 3.1, 4.1). That trapping efficiency, residence time, and latitudinal extent are not modeled, so the diagnostic alone cannot predict the direction or magnitude of limb asymmetry. The GCM comparisons in Sec. 3.1 support the order-of-magnitude threshold for the three tested planets, but they do not validate the assumption across the parameter range the framework is intended to cover (Sec. 5). Thus the central claim is a conditional statement whose non-trivial condition (gyre trapping) is untested by the analytical framework.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper introduces an analytical framework, Ψ_HALD = τ_v/τ_a, intended to estimate the maximum photochemical haze particle radius that can be advected from the evening to the morning limb of a hot Jupiter before vertical removal dominates. The framework combines analytic scalings for jet strength (Zhang & Showman 2017), vertical velocity (Tan 2022), gravitational settling, and radiation pressure, and is compared with new UK Met Office UM GCM simulations for WASP-39b-like planets at two surface gravities, as well as with published GCM results for HD189733b, HD209458b, and WASP-39b. The authors state that the Ψ_HALD = 1 contour gives the maximum radius for which haze particles can survive transport to the morning limb, and, assuming subsequent trapping in nightside gyres, produce a higher or comparable morning-limb haze concentration.","tokens_in":22858,"tokens_out":8516,"duration_ms":70908,"significance":"If the framework is accepted as a transport-survival diagnostic, it is a genuinely useful, computationally inexpensive tool for planning JWST limb observations and for surveying which particle sizes may contribute to morning- versus evening-limb haze asymmetries. The paper's strengths include a transparent dimensionless criterion, original GCM experiments with passive-tracer particles as large as 1600 nm, an explicit smooth fit to the Socrates-computed Q_pr, and a candid discussion of limitations in Secs. 4.1–4.3. The framework is not fitted to the limb-asymmetry outputs: the Q_pr broken power law is fitted to precomputed optical efficiencies, not to the predicted transition radii, so the order-of-magnitude agreement with GCMs is genuine evidence for the survival threshold. The central caveat is that the framework as derived does not by itself produce a morning-limb enhancement; that step is imported through the unmodelled nightside-gyre trapping assumption.","major_comments":[{"comment":"The paper defines Ψ_HALD = 1 as 'the estimated maximum radius for which haze particles can survive the hemispheric transport' (Sec. 3.1) and then equates this with the maximum radius 'capable of producing a higher or comparable haze concentration over the morning limb' (Fig. 3 caption, Table 2). However, the 1D kinematic model yields χ_m/χ_e = exp(-1/Ψ_HALD) < 1 for every finite Ψ_HALD; at Ψ_HALD = 1 the prediction is χ_m = 0.37 χ_e. The morning-limb enhancement is therefore not a consequence of the framework; it is imported by the assumption that particles reaching the morning limb are trapped in nightside gyres (Sec. 4.1), a process that is not modelled. The central claim should be reframed as a transport-survival criterion, with the limb-asymmetry prediction stated as an additional, untested assumption, or the trapping efficiency should be explicitly parameterised.","section":"§3.1 and Eq. (26)"},{"comment":"The Ψ_HALD = 1 contours used to estimate r_max for HD189733b, HD209458b and WASP-39b are computed with a single generic parameter set (R_p = R_J, T = 1000 K, F* = 10^6 W m^-2, Ω = 2×10^-5 s^-1), not with the actual radius, equilibrium temperature, instellation, and rotation rate of each planet or of the GCM setups used for comparison. Since U_rms depends on these parameters (Eqs. 2–5) and W depends on U_rms (Eq. 10), the reported agreement in Table 2 may largely reflect the chosen generic values rather than the framework's predictive skill. The validation should be repeated with planet-specific parameters, or the comparison should be presented as a test of a 'generic hot-Jupiter' scaling only.","section":"§3.1, Fig. 3 and Table 2"},{"comment":"The quantitative predictions are evaluated at p = 0.01 mbar, a pressure chosen without explicit justification and well above the transmission photosphere introduced in Sec. 2.1.1 (p_{τ=1} ~ 10 mbar). The predicted r_max is strongly pressure-dependent (33.1 nm at 0.005 mbar vs 67.1 nm at 0.01 mbar for g = 20 m s^-2), so the headline numbers depend on a hand-picked level. The authors should either provide a rule for selecting p for a given observation, present the Ψ_HALD = 1 contour as a function of pressure, or restrict claims to the pressures at which GCM diagnostics are evaluated.","section":"§3.1 and Table 2"},{"comment":"For the WASP-39b-like low-gravity case (g = 4.3 m s^-2), the GCM simulations bracket the transition only as 50 nm < r_max(p = 0.01 mbar) < 1200 nm, and the framework's estimate of 1450 nm lies above the GCM upper bound. The statement that the framework is 'broadly consistent' to order of magnitude is fair, but the direction of the discrepancy (framework overpredicting r_max) is not discussed. Given the paper's claim that the framework performs best for higher-gravity planets, the manuscript should explicitly address why the survival threshold is overestimated in the low-gravity regime, for example by discussing the neglect of westward day-side flow or of vertical mixing (Secs. 4.1–4.2).","section":"§3.1, Figs. 4–5"}],"minor_comments":[{"comment":"Line 1: 'a analytical framework' should be 'an analytical framework'.","section":"Abstract"},{"comment":"The caption says '1D kinematic model based on Equation (6)', but the model is defined by Eqs. (21)–(26); Eq. (6) is the vertical advection timescale. The reference should be to Eq. (25) or to Sec. 2.1.3.","section":"Fig. 6 caption"},{"comment":"The broken power-law expression has a complex exponent; please ensure the notation in Eq. (A1) and the caption of Fig. A1 is consistent, and define Δ in the text. Also, the fit parameter r_b is quoted in metres (1.1994×10^-7 m), while the text uses nanometres; please state the unit explicitly.","section":"Appendix A, Eq. (A1)"},{"comment":"For HD209458b the GCM constraint 'r_max(p = 0.01 mbar) ≥ 1.5' is only a lower bound because only a 1.5 nm tracer was simulated; this weak data point should be flagged as a lower limit rather than a bracketing constraint.","section":"Table 2"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is likely within the scope of MNRAS and the underlying scaling framework is worth publishing, but the abstract and conclusions overstate the limb-asymmetry claim relative to what the analytical model actually derives. I would ask the authors to either add a simple parameterisation of gyre trapping or consistently reframe the paper as a transport-survival diagnostic, and to repeat the validation with planet-specific input parameters before acceptance. The choice of evaluation pressure should also be justified or made explicit as a sensitivity parameter."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper earns its keep. It gives you a dimensionless ratio, Psi_HALD = tau_v/tau_a, that collapses the competition between horizontal advection and vertical removal into a single contour, and it validates that contour against 3D GCMs across a range of surface gravities. The agreement is real order-of-magnitude stuff: 50-1200 nm bracketing your 1450 nm estimate at 0.01 mbar for WASP-39b, and 60-70 nm bracketing your 33.1 nm at 0.005 mbar for the g=20 case. No parameters were fitted to make the GCMs agree; the Q_pr fit is to an optical-efficiency curve, not to limb asymmetries. That is the right kind of validation for a cheap framework meant to triage more expensive simulations. The soft spot is the one the stress-test names, and it is real: Psi_HALD=1 is a transport-survival threshold, not a limb-asymmetry threshold. Your own 1D model, Eq. (26), gives chi_m/chi_e = exp(-1/Psi_HALD), which is less than 1 for every finite Psi. So the diagnostic alone predicts depletion, not enhancement; the morning-limb enhancement enters through the unmodeled nightside-gyre trapping step. You do not hide this - it is stated in the abstract, Sec. 3.1, and Sec. 4.1 - but it means the central claim is explicitly conditional. The referee should ask whether that condition is doing too much load-bearing work, especially when you want to use this to plan JWST observations. Two smaller issues. First, the evaluation pressure. You compute the headline numbers at 0.01 mbar while assuming p_tau=1 ~ 10 mbar for the wind scaling; there is no demonstrated mapping between that contour and the pressure probed by transmission spectroscopy. The maximum radius shifts by a factor of two between 0.01 and 0.005 mbar, so this is not a trivial detail. Second, there are no quantitative uncertainties on the input scalings or the predicted radii - given how many hand-chosen parameters go in (p_ref, Delta ln p, W prefactor, etc.), a sensitivity span would be worth adding. Those caveats do not sink the paper. It delivers what it promises: a cheap, physically motivated guess for which particle sizes can reach the morning limb, with honest about its limits. I would send this to a serious referee. The referee should push for a clearer statement that Psi_HALD diagnoses transport survival, not limb asymmetry per se, and for a sensitivity analysis on pressure and the gyre-trapping assumption. I would cite it in my own work on hot-Jupiter hazes.","headline":"A genuinely useful first-order diagnostic for predicting which haze sizes can reach the morning limb, but the jump from survival to a morning-limb enhancement rests on unmodeled gyre trapping - a limitation the paper states clearly.","tokens_in":751,"tokens_out":1275,"would_cite":true,"duration_ms":27767,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A single ratio of vertical-to-horizontal transport timescales predicts the largest haze particle that can accumulate over the morning limb of a hot Jupiter.","keywords":["hot Jupiters","photochemical haze","limb asymmetry","transmission spectroscopy","haze transport","gravitational settling","radiation pressure","analytical framework"],"falsifier":"A concrete test would be to run a 3D GCM with passive tracers of several particle sizes for a high-gravity hot Jupiter and measure whether the transitional radius between morning-limb and evening-limb enhancement falls within the range predicted by the ΨHALD = 1 contour at the relevant pressure, accounting for radiation pressure. Alternatively, observations of a hot Jupiter's transmission spectrum at two limbs that show no morning-limb enhancement for particle sizes the framework says should survive would falsify the trapping assumption or the timescale scalings.","tokens_in":22334,"feed_emoji":"🌫️","tokens_out":5966,"duration_ms":42276,"temperature":0.7,"pith_summary":"This paper argues that the relative concentration of photochemical haze between the morning and evening limbs of a hot Jupiter is set by a simple competition: how fast horizontal winds carry particles around the planet versus how fast they fall or are pushed downward. The authors encode this competition in a dimensionless ratio, ΨHALD, defined as the vertical removal timescale divided by the horizontal advection timescale. When ΨHALD exceeds 1, particles can traverse the nightside hemisphere before being removed, and the paper predicts that—assuming nightside gyres trap them—the morning limb should hold as much or more haze than the evening limb. The framework gives a maximum particle radius for morning-limb enhancement that depends on gravity, planetary radius, stellar flux, and particle size. If correct, it provides a quick, inexpensive way to plan transmission-spectroscopy limb observations and to decide when full 3D simulations are needed.","feed_headline":"A single ratio finds the haze size that survives to the morning limb","feed_subtitle":"A dimensionless ratio tells observers which haze particles can reach a hot Jupiter's morning limb.","key_machinery":"The key object is the dimensionless ratio ΨHALD, built from three timescales: horizontal advection τ_a = πR_p/U_rms, vertical advection τ_w = H/W, and gravitational settling τ_s = H/V_s, where V_s is a Stokes settling velocity corrected for the Cunningham slip factor and, optionally, enhanced by radiation pressure through a factor (1+β_rad)g. The vertical removal timescale τ_v is the inverse sum of 1/τ_w and 1/τ_s. ΨHALD = τ_v/τ_a is the ratio of vertical removal to horizontal transport; the ΨHALD = 1 contour gives the maximum particle radius that can make the hemispheric trip to the morning limb. Analytical scalings for U_rms and W are adopted from circulation theory so the framework can be","core_discovery":"The central discovery is the introduction of the Haze Asymmetric Limb Distribution ratio ΨHALD = τ_v/τ_a, where τ_v combines gravitational settling and large-scale downward advection (inverse-sum of their timescales) and τ_a is the hemispheric advection time πR_p/U_rms. The ΨHALD = 1 contour marks the estimated maximum particle radius that can be transported from the evening to the morning limb before vertical removal dominates. Comparing against 3D GCM simulations of WASP-39b and previously published models of HD189733b and HD209458b, the framework yields order-of-magnitude agreement for the particle radius separating morning-limb-enhancement from evening-limb-dominance. The paper further s","pith_inferences":["Because the framework reduces circulation to two timescales, its predictive power is strongest in the settling-dominated regime; a natural extension would be to replace the single hemispheric-average downward velocity with a distribution of vertical velocities to capture local upwelling that could counteract settling.","The assumption that particles reaching the morning limb are trapped in nightside gyres could be tested directly in existing GCM output by diagnosing gyre residence times; if residence times are short, the ΨHALD = 1 threshold would overestimate the likelihood of morning-limb enhancement.","The smooth broken-power-law fit for the radiation-pressure efficiency could be propagated to other haze compositions, allowing the framework to be applied to different aerosol types without recalibrating the radiative transfer."],"forward_implications":["For high-gravity hot Jupiters (e.g., HD189733b), only particles smaller than roughly 10–100 nm can survive transport to the morning limb; larger particles preferentially settle over the nightside, giving an evening-limb haze enhancement.","For low-gravity hot Jupiters (e.g., WASP-39b), even micrometer-sized particles can be advected to the morning limb, so limb asymmetries in transmission spectra can persist for large particle sizes.","Inclusion of radiation pressure reduces the maximum survivable radius by about 35–70% for the test cases, meaning radiation pressure can be the deciding factor in whether a given haze particle reaches the morning limb.","The framework's order-of-magnitude agreement with GCM simulations suggests it can be used to pre-screen targets and particle sizes for JWST limb-asymmetry observations, and to identify when expensive 3D simulations are necessary."],"fun_headline_variants":["Haze ratio flags which particles reach the morning limb","One dimensionless ratio estimates morning-limb haze size","Gravity sets if haze makes it to a hot Jupiter's morning limb","A simple ratio predicts haze arrival at a hot Jupiter's dawn","Quick ratio gives first-order size for morning-limb haze"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The entire morning-limb-enhancement prediction rests on the unmodeled assumption that haze particles which reach the morning limb are subsequently trapped and accumulate in nightside gyres; if gyre trapping is weak or particles recirculate, ΨHALD > 1 would not translate into a higher morning-limb haze concentration.","fun_headline_variants_meta":{"raw":{"variants":["Haze ratio flags which particles reach the morning limb","One dimensionless ratio estimates morning-limb haze size","Gravity sets if haze makes it to a hot Jupiter's morning limb","A simple ratio predicts haze arrival at a hot Jupiter's dawn","Quick ratio gives first-order size for morning-limb haze"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000613,"raw_usage":{"total_tokens":2726,"prompt_tokens":819,"completion_tokens":1907,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":563,"completion_tokens_details":{"reasoning_tokens":1825}},"tokens_in":563,"tokens_out":1907,"duration_ms":14744,"temperature":1.0,"reasoning_tokens":1825,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T00:53:22.552786+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A concrete test would be to run a 3D GCM with passive tracers of several particle sizes for a high-gravity hot Jupiter and measure whether the transitional radius between morning-limb and evening-limb enhancement falls within the range predicted by the ΨHALD = 1 contour at the relevant pressure, accounting for radiation pressure. Alternatively, observations of a hot Jupiter's transmission spectrum at two limbs that show no morning-limb enhancement for particle sizes the framework says should survive would falsify the trapping assumption or the timescale scalings.","supporting_citations":[],"review_version":1}