{"id":"a79abb06-02a9-4e4c-be07-9c01579a4cf2","arxiv_id":"2508.20175","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":8.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"Stellar radiation pressure markedly accelerates aerosol settling in hot exoplanets, reducing particle size and mass concentration, which steepens optical transmission slopes and increases near-infrared molecular feature amplitudes.","lead":"This paper shows that starlight can push aerosol particles in hot exoplanet atmospheres with a force comparable to the planet's gravity, often exceeding it by ten to twenty times. That extra push shrinks the particles and thins the haze, which changes how these planets look in transit and in reflected light.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Monodisperse, perfect-sticking coagulation model omits radiation-pressure-enhanced differential settling, so predicted size/mass-concentration reductions (Eqs. 10-11) are not yet robust.","rationale":"The reader's weakest assumption already identifies sticking efficiency, fragmentation, and broad size distributions as the key simplifications. I agree that this is the load-bearing area, but I sharpen it: the specific omission of differential-settling-driven coagulation is the most direct threat because radiation pressure scales the differential settling velocity by (1+β_p), so the direction and magnitude of the effect cannot be inferred from the monodisperse Brownian-only treatment. Other acknowledged simplifications (isothermal atmosphere, parameterized haze source, 2D equatorial band) are more likely to affect quantitative spectral retrievals than to reverse the basic force-balance argument. The beta_p estimates and the claim that radiation pressure can exceed gravity are on much firmer ground; the weakness is in the microphysical link from acceleration to particle size and concentration. The proposed resolved-bin test would settle whether the monodisperse assumption hides a compensating coagulation channel. Since the reader already issued a CONDITIONAL verdict conditioned on sensitivity to these simplifications, my concern does not move the verdict, but it does identify the specific test that should be run.","tokens_in":33266,"tokens_out":20126,"duration_ms":239880,"concrete_test":"Replace the monodisperse growth-settling model in Sections 3.3 and 4 with a resolved size distribution (e.g., 30-50 logarithmically spaced bins from 1 nm to 10 μm) using the same Ormel & Min (2019) relative-velocity kernels plus an explicit differential-settling term and a plausible sticking-efficiency range (0.1-1). Rerun the HAT-P-65 b tholin and soot cases at β_p = 0 and β_p = 10, with identical production rates and diffusion constants. Compare the mass-weighted mean particle size and haze mass concentration at the transmission photosphere, and the resulting optical slope and 1.4 μm water-feature amplitude. If the β_p = 10 case still yields the same qualitative reduction in size and concentration, the concern is resolved; if the differences shrink or reverse, Eqs. 10-11 and the spectral predictions in Section 5 need revision.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central size and mass-concentration scalings (Eqs. 10-11) are derived in Section 3.1 by equating the settling time to a collision time (Eq. 9) that includes only the Brownian relative velocity of equal-mass particles. In a real haze the size distribution is broad, and collisions between different-sized grains are driven substantially by differential settling. The settling velocity scales as (1+β_p) a, so radiation pressure amplifies these differential settling velocities. The 2D model (Section 3.3) instead evolves a single representative particle size per cell (Eq. 36) and treats the growth of freshly injected 10^-3 μm seeds only through a simplified mass-weighted average or sweep-up estimate (Section 3.4). This means the differential-settling coagulation channel is effectively absent from both the analytic scalings and the numerical simulations. If that channel is important, the net effect of radiation pressure on particle size and mass concentration is not necessarily the (1+β)^-4/9 and (1+β)^-5/9 scalings, and the quantitative transmission-spectrum predictions (Section 5) would change. The assertion in Section 6.1 that the basic result is unaffected by these simplifications is plausible but not demonstrated.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper proposes that stellar radiation pressure exerts a force on aerosol particles in hot exoplanet atmospheres that is often comparable to or larger than planetary gravity, with beta_p up to ~10-30 for low-gravity, highly irradiated planets. Assuming hazes with a monodisperse size and a growth-settling equilibrium, the authors derive analytic scalings in which radiation pressure reduces particle size and mass concentration; they then simulate haze dynamics in a 2D equatorial band with a single representative particle size, compute transmission spectra, and predict steeper optical slopes and larger near-IR molecular feature amplitudes at higher radiation pressure. They also discuss albedo effects, morning/evening terminator asymmetries, aerosol mass loss, and population-level correlations.","tokens_in":33622,"tokens_out":21197,"duration_ms":246404,"significance":"If the central scalings survive closer scrutiny, radiation pressure would be a genuinely new physical process for aerosol evolution on hot exoplanets, with direct consequences for the interpretation of HST and JWST transmission spectra and albedo measurements. The paper's analytic derivation is clean and does not fit observational data; the Mie opacities are computed from published laboratory optical constants; the simulation code is publicly available; and the main predictions (beta_p-driven steepening of optical slopes and stronger molecular features) are, in principle, testable with current observations. The paper is also appropriately cautious about implicit correlations between radiation pressure strength and other planetary/stellar parameters.","major_comments":[{"comment":"The analytic growth-settling equilibrium treats collisions as driven only by Brownian relative motion of equal-mass particles. This assumption is not justified at the pressures and sizes of interest. For the HAT-P-65 b analogue at P*=1 microbar with beta_p=10, Eq. (8) gives settling velocities of about 6e3 cm/s for 0.1 micron grains and 1.3e4 cm/s for 0.2 micron grains, so differential settling between different-sized grains is about 6e3 cm/s, while v_BM in Eq. (9) is about 20 cm/s. Since v_BM scales as a^{-3/2} and v_set as (1+beta_p) a / rho_gas, the hierarchy worsens for larger grains and lower gas densities. Radiation pressure multiplies differential settling velocities by (1+beta_p), so the scalings in Eqs. (10)-(11), and their use in Section 5, are not demonstrated for a realistic size distribution. Please provide a quantitative validity condition for the Brownian-only regime, or g","section":"Section 3.1, Eqs. (8)-(11)"},{"comment":"The numerical model evolves a single representative particle size per cell, yet the growth time is stated to come from Ormel & Min (2019), including differential particle drift. With only one size per cell it is unclear how the differential-drift relative velocity is evaluated; if it is evaluated between identical sizes it vanishes, and the model cannot capture the radiation-pressure-amplified differential-settling channel highlighted above. The separate treatment of freshly injected 10^-3 micron seeds (sweep-up vs. self-coagulation) addresses only the injection step, not the subsequent growth of the resident population. Please specify the assumed size distribution or effective relative velocity used in the Ormel & Min formula, and test the sensitivity of the size/concentration profiles (Figures 9, 11-12) and transmission spectra (Figures 14-16) to the single-size approximation. This is","section":"Section 3.3-3.4, Eq. (36)"},{"comment":"The collision model assumes perfect sticking ('all collisions result in growth') and no fragmentation. Radiation pressure can produce collision velocities of order km/s at the low pressures considered, where sticking efficiencies are plausibly below unity and fragmentation or erosion may become important. This could alter the growth-settling equilibrium and the predicted size scalings. The manuscript does not discuss this channel, even though it is directly relevant to the claim that the qualitative result is robust. At minimum, the authors should state the expected collision velocities and justify the perfect-sticking assumption for the parameter range simulated.","section":"Section 3.1, Eq. (9) and Section 6.1"}],"minor_comments":[{"comment":"The prefactor of 5 assumes F* = 4 sigma Teq^4 (full isotropic re-radiation) and zero albedo. State this explicitly when introducing Teq.","section":"Eq. (3)"},{"comment":"The caption refers to a 'HATP-67 b analogue' while the text uses HAT-P-65 b. Please correct the typo.","section":"Figure 9 caption"},{"comment":"The text says 'HAT-P-67 b having the largest value', while the square point in Figure 2 is HAT-P-65 b. Harmonize the notation.","section":"Section 2.1"},{"comment":"Please provide the fitted values of the eight alpha coefficients or a pointer to where they can be obtained; the functional form alone is not sufficient for reproduction.","section":"Eq. (40)"},{"comment":"Define the Rayleigh index b more precisely: Qpr is proportional to (a/lambda*)^{1+b} in the small-particle limit, and state the particle size range over which this approximation is used.","section":"Section 3.2"},{"comment":"Specify how tau_pro is computed (wavelength band, gas and particle opacity) and how the production shut-off exp(-tau_pro) is normalized relative to Eq. (39).","section":"Section 5.1"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is promising and likely of interest to the journal. My main concern is whether the Brownian/equal-mass analytic model and the single-size numerical model capture the dominant coagulation channel. I would be willing to accept after the authors either include a differential-settling term in the analytic scalings with a validity check, or explicitly test the single-size approximation against a two-moment or multi-size calculation. The robustness claim in Section 6.1 is asserted rather than demonstrated."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper's central claim is new and, I think, right: for highly irradiated, low-gravity planets, stellar radiation pressure on ~0.1-1 micron aerosols can exceed planetary gravity by an order of magnitude. That is a simple, parameter-free force comparison (Eq. 2-3, Figure 2), and it holds up. The consequence—faster settling, smaller particles, lower mass concentrations—is derived cleanly in the 1D growth-settling equilibrium, and the scalings (Eqs. 10-11, 13-16) are internally consistent and worth having. The 2D simulations are a sensible first step: they show the qualitative effects on optical slopes, terminator asymmetries, and feature amplitudes. The code is public, and the opacity fits are to Mie theory, not to the predicted observables, which is good practice. The paper is honest about its simplifications: isothermal atmosphere, parameterized production, single representative particle size, 2D band. For an initial exploration, that is acceptable.\n\nThe soft spots are where the stress-test note lands, and they are real but not fatal. The analytic scalings equate growth and settling time using only Brownian relative velocities of equal-mass particles. In a broad size distribution, differential settling drives coagulation, and radiation pressure amplifies differential settling velocities because settling velocity scales with (1+beta_p) times particle size. That channel is absent from both the analytic equilibrium and the single-size numerical model. So the specific exponents in Eqs. 10-11 and the quantitative spectral predictions in Section 5 should be treated as provisional. The Section 6.1 claim that the basic result is robust to these choices is plausible but not demonstrated; a proper sensitivity test with a moment or bin scheme would settle it. The isothermal assumption and the parameterized zonal wind are minor by comparison.\n\nAlso minor but worth noting: the paper's population-level correlation with the 1.4 micron water feature amplitude is speculative, as the authors themselves say, because they hold planetary parameters fixed while varying beta_p, and production rate probably correlates with radiation pressure in reality. That is a fair caveat, not a flaw.\n\nWho gets value: modelers working on haze microphysics and observers trying to interpret steep optical slopes or terminator asymmetries. The paper deserves a serious referee—it introduces a physical process that no one else has considered in this context, and it pushes the field toward a testable prediction. I would engage with it, but I would treat the size scalings as hypotheses to be stress-tested rather than established results.","headline":"Radiation pressure on aerosols is a genuinely new physical idea for hot exoplanet atmospheres, and the core force-balance argument is sound; the quantitative scalings and spectral predictions are real but rest on a simplified coagulation model, so the paper deserves a serious referee but also a clear sensitivity discussion.","tokens_in":34027,"tokens_out":654,"would_cite":true,"duration_ms":10662,"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":"Starlight exerts up to 30 times the acceleration of gravity on sub-micron haze particles in hot exoplanet atmospheres, and this extra push makes the particles smaller and more dilute, steepening optical transmission slopes and strengthening","keywords":["exoplanet atmospheres","aerosols","radiation pressure","transmission spectroscopy","haze microphysics","growth-settling equilibrium","hot Jupiters","sub-Saturns"],"falsifier":"Observe a sample of hot gas giants with homogeneous transmission spectra and measure the normalised 1.4 μm water amplitude and optical slope as a function of β_p; if planets with high β_p do not show larger water amplitudes and steeper slopes after accounting for scale height, the radiation-pressure mechanism is not the controlling factor. Alternatively, a laboratory measurement showing that soot/tholin aggregates bounce or fragment at the relevant collision speeds would falsify the perfect-sticking growth-settling equilibrium that the size scalings depend on.","tokens_in":33138,"feed_emoji":"☀️","tokens_out":8116,"duration_ms":74859,"temperature":0.7,"pith_summary":"This paper argues that stellar radiation pressure—the momentum of starlight absorbed and scattered by small particles—is a leading, previously neglected force in the atmospheres of hot exoplanets. For low-density, highly irradiated planets (sub-Saturns and inflated giants), the radiative acceleration on sub-micron haze particles can exceed the planet's gravity by a factor of 10–30, so haze particles settle much faster than they would under gravity alone. In the growth-settling equilibrium the authors derive, that faster settling leaves less time for coagulation, so radiation pressure produces smaller particles with lower mass concentrations. Using analytic scalings and 2D equatorial-band simulations, the paper shows these smaller, thinner hazes steepen optical transmission slopes, make near-infrared molecular features less muted, and create asymmetric morning/evening terminator spectra. If this is right, radiation pressure is needed to interpret transmission spectra and albedos of the very planets best suited for atmospheric characterisation.","feed_headline":"Starlight up to 30x gravity reshapes hot exoplanet hazes","feed_subtitle":"The extra acceleration makes haze particles shrink and thin out, steepening optical slopes and unmuting water features.","key_machinery":"The load-bearing mechanism is the radiation-pressure-to-gravity ratio β_p (Equation 2), which parametrises the extra downward acceleration on an aerosol particle. The argument then runs through the growth-settling equilibrium: the terminal settling velocity (Equation 8, Epstein drag with acceleration (1+β_p)g) and the coagulation timescale (Equation 9, ballistic Brownian collisions with perfect sticking) are equated, giving the particle size and mass concentration scalings (Equations 10–11). This equilibrium determines everything downstream: the aerosol opacity profile, the photospheric pressure, and finally the transmission spectrum. The 2D equatorial-band model (Section 3.3) carries this m","core_discovery":"The paper's central claim is that the ratio of radiative acceleration to planetary gravity, β_p, is typically larger than one for ~0.1–1 micron aerosols on highly irradiated giant planets, often reaching values of 10–30 (Equation 3, Figure 2). Because radiation pressure acts in the stellar direction, it adds to gravity in the settling of haze particles, raising their terminal velocity. Equating the coagulation timescale (Equation 9, assuming perfect sticking) with the settling timescale (Equation 8) gives the growth-settling equilibrium: particle size a ∝ (1+β_p)^(-2/9) and mass concentration X_p ∝ (1+β_p)^(-5/9) (Equations 10–11). Thus radiation pressure makes particles smaller and the aero","pith_inferences":["The paper's scalings suggest a direct observational test: a sample of hot giant planets with well-determined masses, radii and temperatures should show normalised water-feature amplitudes increasing with β_p even after controlling for scale-height effects; current observations of separate morning/evening terminator spectra could check the predicted asymmetry direction.","If the same acceleration applies to condensate clouds, as the paper speculates, radiation pressure would shrink cloud particles too, potentially explaining super-Rayleigh optical slopes without invoking extremely strong vertical mixing.","The perfect-sticking, single-size growth-settling equilibrium is the load-bearing premise; if laboratory studies of soot/tholin aggregates find low sticking efficiencies or fragmentation, the size scaling with β_p weakens and the spectral predictions would need revision.","The radiative feedback loop described in Section 6.1 (smaller particles → lower opacity → hotter atmosphere → longer settling time → larger particles) could create bistable or time-variable haze states, which would appear as transit-depth variations across epochs."],"forward_implications":["Microphysical models that omit radiation pressure will overestimate haze particle sizes and mass concentrations on highly irradiated planets, biasing retrievals of molecular abundances and cloud properties.","The paper predicts a population-level correlation: among hot gas giants, the normalised 1.4 μm water feature amplitude should increase with β_p (roughly with equilibrium temperature and inverse surface gravity), helping explain the scatter in measured feature amplitudes.","Morning and evening terminator spectra should differ systematically: the evening terminator is hazier, while the morning terminator has a steeper optical slope that can become super-Rayleigh under strong radiation pressure.","Radiation pressure can drive an outward flow of aerosols; for soot-like hazes the mass-loss rate is about 1.6% of the production rate, enough to deplete the upper-atmosphere haze precursor reservoir over the planet's lifetime unless replenished.","Radiation pressure narrows the haze production rates that produce high albedos, which could account for the lack of strong correlations between observed exoplanet albedos and standard planetary and stellar parameters."],"supporting_citations":[{"why":"Supplies the standard formulation of radiation pressure on small particles and the definition of the efficiency Q_pr.","marker":"Burns et al. 1979"},{"why":"Supplies the coagulation growth timescale, including ballistic and diffusive Brownian motion, used in both the analytic equilibrium and the 2D model.","marker":"Ormel & Min 2019"},{"why":"Provides soot optical constants used to compute radiative efficiencies for haze particles.","marker":"Lavvas & Koskinen 2017"},{"why":"Provides tholin optical constants used to represent organic haze.","marker":"Khare et al. 1984"},{"why":"Provides silicate (Mg2SiO4) optical constants for condensate cloud particles.","marker":"Jäger et al. 2003"},{"why":"Provides the Mie scattering code (bhmie) used to calculate absorption and scattering efficiencies and asymmetry parameters.","marker":"Bohren & Huffman 1998"},{"why":"Provides the standard no-radiation-pressure scaling limit that the paper extends and compares against.","marker":"Ohno & Kawashima 2020"},{"why":"Defines the normalised 1.4 micron water feature amplitude and the population sample the paper's correlation is framed against.","marker":"Dymont et al. 2022"},{"why":"Supplies the assumption that the 2D aerosol distribution can be treated as uniform on each terminator when computing transmission spectra.","marker":"Kempton et al. 2017"}],"fun_headline_variants":["Radiation pressure shrinks hot exoplanet haze particles","Starlight's push reshapes hazes on hot exoplanets","Radiative acceleration controls hot exoplanet aerosols","Hot planet haze sizes set by radiation pressure","Starlight pressure dominates hot exoplanet haze dynamics"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The predictions rest on a growth-settling equilibrium in which every collision sticks and the aerosol population is described by a single representative size; if sticking is inefficient, fragmentation occurs, or the size distribution is broad, the scalings and spectral signatures would change.","fun_headline_variants_meta":{"raw":{"variants":["Radiation pressure shrinks hot exoplanet haze particles","Starlight's push reshapes hazes on hot exoplanets","Radiative acceleration controls hot exoplanet aerosols","Hot planet haze sizes set by radiation pressure","Starlight pressure dominates hot exoplanet haze dynamics"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000445,"raw_usage":{"total_tokens":2116,"prompt_tokens":803,"completion_tokens":1313,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":547,"completion_tokens_details":{"reasoning_tokens":1233}},"tokens_in":547,"tokens_out":1313,"duration_ms":10137,"temperature":1.0,"reasoning_tokens":1233,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T15:13:28.055336+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Observe a sample of hot gas giants with homogeneous transmission spectra and measure the normalised 1.4 μm water amplitude and optical slope as a function of β_p; if planets with high β_p do not show larger water amplitudes and steeper slopes after accounting for scale height, the radiation-pressure mechanism is not the controlling factor. Alternatively, a laboratory measurement showing that soot/tholin aggregates bounce or fragment at the relevant collision speeds would falsify the perfect-sticking growth-settling equilibrium that the size scalings depend on.","supporting_citations":[],"review_version":1}