Pith. sign in

REVIEW 3 major objections 4 minor 1 cited by

The Effects of Cuboid Particle Scattering on Reflected Light Phase Curves: Insights from Laboratory Data and Theory

T0 review · 3 major / 4 minor · reviewed 2026-08-06 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read Simulated reflected-light phase curves of GJ 1214b stay within 3 ppm when cuboid KCl particle scattering replaces the TTHG approximation, with only perfectly regular cuboids exceeding that.

desk verdict A carefully scoped comparison showing TTHG approximates cuboid KCl scattering well, but the headline 'less than 3 ppm' rests on an unmeasured backscatter extrapolation and missing error propagation. read the letter →

arxiv 2507.05485 v1 pith:IUR6Y22N submitted 2025-07-07 astro-ph.EP astro-ph.IM

Colin D. Hamill , Alexandria V. Johnson , Matt Lodge , Peter Gao , Rowan Nag , Natasha Batalha , Duncan A. Christie , Hannah R. Wakeford This is my paper · ORCID
classification astro-ph.EPastro-ph.IM
keywords exoplanetatmospherescloudparticlesKClcuboidsphasefunctionstwo-termHenyey-GreensteindiscretedipoleapproximationreflectedlightcurvesGJ1214b
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper asks whether the simplified two-term Henyey–Greenstein (TTHG) scattering functions built into the PICASO radiative-transfer code are good enough for cloudy exoplanets whose condensates are not spheres, focusing on KCl, which crystallizes as cubes and irregular cuboids. It compares TTHG phase functions against laboratory measurements and discrete dipole approximation (DDA) calculations for three KCl particle size distributions at 532 nm, then inserts each phase function into models of the sub-Neptune GJ 1214b. The simulated reflected-light phase curves agree to within 3 ppm for cubic and irregular cuboid particles, with only perfectly regular cuboids producing larger deviations of 4–10 ppm near secondary eclipse. The conclusion is that TTHG, although a poor stand-in for Mie scattering from spheres, happens to reproduce the weak backscattering of cuboid particles well enough for current reflected-light observations.

What carries the argument

The two-term Henyey–Greenstein (TTHG) phase function is the object under test: an analytic scattering profile with a forward lobe, a backscattering lobe, and a weight $f$, parameterized by $g_f$, $g_b$, and $f$ as in PICASO. Its suppressed backscattering relative to Mie spheres is what makes it a closer match to cubes and irregular cuboids. The comparison machinery is a pipeline of four inputs: laboratory phase functions measured between 20° and 169° and extrapolated to all angles using Mie theory forward and a quadratic fit with cubic splines backward; DDA phase functions for cubes, regular cuboids, and irregular cuboids computed with equal-volume dipole shapes; the cloud model Virga supplying vertical cloud structure; and the PICASO radiative-transfer model, which turns each phase function into a reflected-light phase curve for GJ 1214b.

What would settle it

A decisive test would be a laboratory measurement of KCl scattering intensity at scattering angles from 169° to 180° for the same three size distributions, using a detector geometry that does not block the laser. If the measured near-180° intensity deviates substantially from the quadratic extrapolation, the GJ 1214b phase-curve differences should be recomputed; alternatively, a reflected-light phase curve of GJ 1214b at 532 nm with per-point precision better than about 1 ppm would distinguish the regular-cuboid DDA model from the TTHG model at secondary eclipse.

Watch

Extended reading notes

Core claim

The central claim is that reflected-light phase curves of a cloudy exoplanet are largely insensitive to the choice of scattering phase function, provided the particles are not perfectly regular cuboids. Using the same cloud distributions, single-scattering albedo, and extinction from Virga, the authors find that replacing PICASO's native TTHG phase functions with laboratory-measured or DDA-computed phase functions for cubic and irregular cuboid KCl changes the 532 nm phase curve of GJ 1214b by less than 3 ppm; regular cuboid shapes change it by 4–10 ppm. Because the differences stay below the typical precision of current exoplanet phase-curve observations, the paper argues that fast TTHG approximations are adequate for modeling reflected light from these nonspherical cloud particles, and that the low backscattering of TTHG accidentally matches the measured dim backscattering of cuboids.

Load-bearing premise

The load-bearing premise is that the unmeasured backscattering intensity at 180 degrees, obtained by a quadratic fit to laboratory data from 135 to 169 degrees and then splined to the endpoint, is close to the true value; if that estimate is wrong, the claimed 2–3 ppm agreement between laboratory and TTHG phase curves near secondary eclipse could change.

Editorial extensions

If this is right

  • If the central claim is correct, PICASO's existing TTHG treatment is a valid shortcut for reflected-light modeling of cubic and irregular cuboid KCl clouds at 532 nm, with errors under 3 ppm.
  • Order-of-magnitude errors in backscattering phase function intensity translate to only about 1 ppm in reflected albedo at most orbital phases, so phase curves are robust to large errors in scattering-angle detail.
  • Regular cuboid shapes are the exception: they can brighten secondary eclipse by 4–10 ppm, so idealized perfectly rectangular particles are the main shape-related risk in reflected-light models.
  • The conclusion is stable across 1× and 100× solar metallicity and across sedimentation efficiencies 0.1 and 1.0, as shown in the appendix.
  • Current observational precision is too coarse to distinguish cubic from irregular cuboid scattering in reflected light at this wavelength.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If this accidental match between TTHG and cuboid backscattering holds for other shapes and wavelengths, reflected-light retrievals may not need shape-resolved phase functions until per-point precision improves below roughly 1 ppm; the same would not necessarily hold for transmission or emission geometry.
  • The quadratic extrapolation of laboratory data from 169° to 180° is the main unmeasured hinge: a laboratory measurement that directly samples near-180° backscattering for the same KCl size distributions would test whether the <3 ppm agreement survives.
  • Because the study fixes single-scattering albedo and extinction at the spherical Virga values, a self-consistent treatment where shape changes those quantities as well could shift albedos beyond 3 ppm even if the phase function alone does not.
  • Extending the comparison to longer wavelengths and larger size parameters would show whether the <3 ppm result is specific to 532 nm and roughly 1 μm KCl particles or generalizes.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 4 minor

Summary. The paper compares phase functions for KCl particles at 532 nm computed with Mie theory, two-term Henyey–Greenstein (TTHG) approximations as implemented in PICASO, laboratory measurements from the ExCESS system, and discrete dipole approximation (DDA) simulations. Three particle size distributions are considered: a narrow small distribution of cubic particles (mean radius ~0.6 μm) and two broader medium and large distributions of irregular cuboids (~0.6 and 1.2 μm). The phase functions are inserted into PICASO with cloud distributions from Virga for a clear-sky GCM of the sub-Neptune GJ 1214b, and reflected light phase curves at 532 nm are compared. The main finding is that phase curves produced with TTHG, laboratory, and DDA phase functions for cubic and irregular cuboid shapes differ by less than 3 ppm, with the exception of regular cuboid DDA cases that differ by 4–10 ppm. The authors conclude that TTHG approximations are adequate for these nonspherical KCl shapes in reflected light at this wavelength, within current observational uncertainties.

Significance. If the result is robust, it is practically important: it would validate the continued use of fast TTHG phase functions in retrieval and forward models for cloudy exoplanets, at least for cubic and irregular cuboid KCl at 532 nm, and it identifies secondary eclipse as the phase angle where shape effects are largest. The study is notable for combining independent laboratory measurements and DDA benchmarks, avoiding circular use of TTHG. The authors are transparent about the narrow parameter space (single wavelength, single condensate, fixed single-scattering albedo and extinction) and provide public data and code, which supports reproducibility. The main quantitative claim, however, depends on an extrapolation of laboratory data beyond the measured angle range, which is not currently propagated into the reported uncertainties.

major comments (3)
  1. [Section 2.1.3 and Section 3.2 (Figure 7)] The laboratory phase functions are only measured to 169°, and the 169°–180° interval is filled by a quadratic fit to the 135°–169° data followed by a cubic spline forced to zero slope at 180°; no uncertainty from this extrapolation is propagated into the phase-curve differences. Since the backscattering direction controls the reflected flux at secondary eclipse (as stated in Section 3.2), the quantitative '<3 ppm' claim in the abstract is not yet secured. Please add a sensitivity study that varies the extrapolation method (e.g., linear extrapolation, different fit ranges, or different spline boundary conditions) and shows how the laboratory-minus-TTHG difference at secondary eclipse changes. If the difference can exceed 3 ppm under plausible alternatives, the abstract and conclusions should be rephrased as a qualitative agreement rather than a quantitative bound.
  2. [Abstract and Section 5] The abstract states that phase curves 'produced using the TTHG, laboratory, and DDA phase curves differ by less than 3ppm,' but Section 5 reports 4–10 ppm differences for the regular cuboid DDA cases. This inconsistency could mislead readers; the abstract should either exclude the regular cuboid cases or explicitly state that these shapes are not considered representative.
  3. [Section 2.1.4 and Appendix A.1] The DDA calculations assume a single identical shape for all particles in each size distribution, and the irregularity level used in the main study is chosen by eye to match SEM images. The appendix shows that phase functions vary with irregularity level. While the authors acknowledge this limitation, the robustness of the reported 2–3 ppm differences for irregular cuboids across plausible shape ensembles is not demonstrated. A brief test with an alternative irregular shape or a short discussion of the expected variation would strengthen the DDA leg of the conclusion.
minor comments (4)
  1. [Section 2.1.3] There are typos in the particle size units: 'mean particle radius of 0.6 ∝m' should be '0.6 μm', and similarly '1.2 ∝m' should be '1.2 μm'.
  2. [Figure 7] The difference panels (d–f) would be easier to interpret if a horizontal line at ±3 ppm were drawn to calibrate the eye against the paper's central claim.
  3. [Abstract] The planet name is written as 'GJ1214b' in the abstract but as 'GJ 1214b' elsewhere; please use the latter consistently.
  4. [Section 3.2] The sentence 'The peak in reflected light at secondary eclipse is due to the increased backscattering intensity ... at a scattering angle of 180°' is slightly imprecise because the exact 180° angle has zero weight in the normalized phase-function integral; the peak is driven by the nearby near-backscattering angles. Consider rewording to 'near 180°'.

Circularity Check

0 steps flagged · score 1.0 of 10

No significant circularity: TTHG functions are constructed from independent Mie asymmetry parameters, and laboratory/DDA phase functions serve as external benchmarks; the <3 ppm agreement is an emergent comparison, not an input.

full rationale

The paper's central comparison is not circular. The TTHG phase functions are generated from the PICASO-default relation (g_f = g, g_b = -g/2, f = 1 - g_b^2) with the asymmetry parameter g computed from Mie phase functions via Equation (8) in Sections 2.1.2 and 2.2.2; none of the TTHG parameters are fitted to the laboratory or DDA phase functions. The laboratory phase functions were measured independently in C. D. Hamill et al. (2024b) with ExCESS, and the DDA calculations use ADDA on shapes chosen from SEM images (by eye) rather than optimized to reproduce the phase curves. Because the PICASO runs hold cloud location, single-scattering albedo, and extinction fixed, the <3 ppm phase-curve differences are an emergent result of comparing independent phase functions, not an input. The acknowledged backscattering extrapolation (quadratic fit to 135-169 degrees followed by cubic splines forced to zero slope at 180 degrees) is a measurement limitation and a source of uncertainty, but it is not a parameter fitted to TTHG, so it does not constitute circularity. Self-citations to Hamill et al. (2024a,b), Lodge et al. (2023), and Lodge (2024) provide the lab data, Mie curves, and shape-generation code; these are independent, reproducible inputs rather than self-supporting uniqueness claims. No load-bearing step reduces to its own inputs by definition.

Assumptions & free parameters 4 free parameters · 6 assumptions · 0 invented entities

The central claim rests on measured and modeled inputs from prior work plus two subjective choices (DDA shape selected by eye; quadratic backscatter extrapolation). No new physical entities are introduced. The robustness check across metallicity and sedimentation values in the Appendix reduces the dependence on any single atmospheric parameter.

free parameters (4)
  • sedimentation efficiency f_sed = 0.1 (main case), 1.0 (appendix)
    Virga sedimentation efficiency; set from prior GJ 1214b studies (Morley et al. 2013, 2015; Christie et al. 2022). The headline <3 ppm differences are unchanged at f_sed = 1.0 in the Appendix.
  • DDA cuboid irregularity level = medium irregularity
    The main irregular cuboid shape was chosen by eye to match SEM images; the Appendix shows slight, medium, and extreme options with modest phase function changes.
  • backscatter quadratic extrapolation coefficients = not reported
    A quadratic fit to lab intensities at 135 to 169 degrees fixes the 180 degree value; used to complete the phase function for PICASO and affects secondary eclipse comparison.
  • K_ZZ,0 eddy diffusion coefficient = 7e2 and 3e3 m2/s (1x and 100x solar)
    From Charnay et al. 2015a; controls vertical cloud distribution in Virga. Input from prior literature, not fitted here.
assumptions (6)
  • domain assumption KCl forms cubic and cuboid crystals during condensation and grinding.
    Premise for wet-generation cubic particles and dry-generation irregular cuboids; supported by Walker et al. 2004 and SEM images from Hamill et al. 2024b (Section 2.1.3).
  • domain assumption Mie theory can extrapolate forward scattering (0 to 20 degrees) of cuboids because forward scattering is shape-insensitive for moderate aspect ratios.
    Section 2.1.3; based on Mishchenko et al. 1996, 1997 and Liu et al. 2003.
  • domain assumption PICASO's TTHG parameterization (gf = g, gb = -g/2, f = 1 - gb^2) represents the default two-stream scattering approximation.
    Section 2.2.2; from Cahoy et al. 2010; the comparison is framed against this specific parameterization.
  • ad hoc to paper DDA with a single identical shape per size distribution adequately represents an ensemble of real particles.
    Section 2.1.4; authors acknowledge all particles share one geometry, which is not realistic but a first step; resonance features in the DDA phase functions show this simplification.
  • domain assumption Virga spherical-particle single-scattering albedo and extinction can be held fixed while varying only phase functions.
    Section 2.2.3 and Discussion; isolates phase function effects but is not self-consistent with nonspherical particles.
  • domain assumption GJ 1214b clear-sky 100x solar GCM from Christie et al. 2022 is adequate for comparing scattering phase functions.
    Section 2.2.1; cloud radiative feedback is excluded because this study only treats KCl scattering, following common practice.

how reviews work

0 comments
Cite this review

Pith. "Pith review of The Effects of Cuboid Particle Scattering on Reflected Light Phase Curves: Insights from Laboratory Data and Theory." pith.science (2026). https://pith.science/paper/IUR6Y22N

@misc{pith2026250705485,
  author       = {Pith},
  title        = {Pith review of: The Effects of Cuboid Particle Scattering on Reflected Light Phase Curves: Insights from Laboratory Data and Theory},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/IUR6Y22N}},
  note         = {Machine review of arXiv:2507.05485}
}
read the original abstract

Understanding the optical properties of exoplanet cloud particles is a top priority. Many cloud condensates form as nonspherical particles and their optical properties can be very different from those of spheres. In this study, we focus on KCl particles, which likely form as cuboids in warm (T=500-1000K) exoplanet atmospheres. We compare the phase functions (at 532 nm wavelength) of KCl particles computed with Mie theory, the two-term Henyey-Greenstein (TTHG) approximation, laboratory data, and the discrete dipole approximation (DDA). Mie theory assumes scattering from spheres, while TTHG functions are used to approximate cloud scattering in two-stream radiative transfer models like PICASO. Laboratory measurements and DDA allow for a robust understanding of scattering from cuboid and deformed cuboid particle shapes. We input these phase functions into PICASO using cloud distributions from the cloud model Virga, to determine how different phase functions can impact the reflected-light intensities of the benchmark sub-Neptune, GJ 1214b. Simulated reflected light phase curves of GJ1214b produced using the TTHG, laboratory, and DDA phase curves differ by less than 3ppm. Our findings suggest that TTHG phase functions may be useful for approximating the scattering intensity of certain cuboid and irregular particle shapes. Future work should expand upon the wavelengths and particles considered to better determine when scattering approximations, like TTHG, may be useful in lieu of more accurate, but time-consuming laboratory measurements and/or nonspherical scattering theory.

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. Circular polarization as a probe of cloud properties and asymmetries in giant exoplanet atmospheres

    astro-ph.EP 2026-08 conditional novelty 6.0 of 10

    Circular polarization of starlight reflected by giant exoplanets is sensitive to cloud particle composition and atmospheric asymmetries, but the predicted signal is at most a few times 1e-4 of the planetary flux and b...

Reference graph

Works this paper leans on

2 extracted references · cited by 1 Pith paper

  1. [1]

    S., & Marley, M

    Ackerman, A. S., & Marley, M. S. 2001, ApJ, 556, 872 Adams, D., Gao, P., Pater, I. de, & Morley, C. V. 2019, ApJ, 874, 61 Adams, D., Kataria, T., Batalha, N., Gao, P., & Knutson, H. 2022, BAAS, 926, 157 Amundsen, D. S., Baraffe, I., Tremblin, P., et al. 2014, A&A, 564, A59 Amundsen, D. S., Tremblin, P., Manners, J., Baraffe, I., & Mayne, N. J. 2017, A&A, ...

    2001

  2. [9]

    Four different combinations of atmospheric metallicity and sedimentation efficiency are shown

    (Top) Reflected light phase curves computed from our GJ 1214b GCM at 532 nm for the small particle size distribution with the TTHG (black), laboratory (orange), or cubic DDA (sky blue) phase functions. Four different combinations of atmospheric metallicity and sedimentation efficiency are shown. Line opacities are used to differentiate the cases further, ...

    arXiv 2025

Pith tools

Reviewed August 6, 2026 · model on record in the stance chip above.