{"id":"88ab5cd3-6a9c-4746-a27f-5ece259eda73","arxiv_id":"2411.14225","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"PICASO can now compute reflected-light phase curves from 3D atmospheres, and for Kepler-7b low-sedimentation multi-species clouds give the closest, though still three-times-too-dim, match to Kepler data.","lead":"Scientists added a new capability to the open-source PICASO model that computes how much starlight a hot Jupiter reflects at each point in its orbit. When applied to Kepler-7b, the model reproduces the observed phase-curve shape only after assuming a specific mix of cloud species, and it still predicts about one-third of the observed brightness.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No validation of the new reflected-light phase-curve mode against an analytic or 1D standard; the factor-of-three Kepler amplitude gap could be a normalization/geometry bug rather than missing cloud physics.","rationale":"The paper honestly reports the factor-of-three amplitude deficit and lists plausible physical causes. The reader's CONDITIONAL verdict centers on the GCM lacking cloud radiative feedback and on post-hoc fsed selection; I agree those are real limitations, but the more load-bearing issue for the paper's stated central claim (a new PICASO capability) is that the new geometric integration is never tested against a known solution. The 3D phase-curve mode adapts the thermal phase-curve routine of Robbins-Blanch et al. (2022), but reflected light has different visibility and geometry, and no equivalent validation is presented. The Appendix only demonstrates grid convergence, not absolute accuracy. Since the benchmark itself shows a factor-of-three offset, one cannot distinguish physical under-prediction from an implementation error in absolute flux normalization or phase-angle weighting. This is a correctness risk, not a disagreement with consensus. The proposed Lambertian-sphere and 1D-comparison tests would settle it. The Kepler conclusions about best fsed and cloud species are also conditional because they depend on the GCM's missing cloud radiative feedback, matching the reader's concern. The verdict therefore remains CONDITIONAL: accept the open-source engineering contribution only with the condition that absolute validation of the new phase-curve mode be added.","tokens_in":21375,"tokens_out":5850,"duration_ms":60796,"concrete_test":"Apply the new routine to a uniform Lambertian sphere: every facet with the same reflectance and no gas opacity. The computed Fp/Fs must follow [sin α + (π - α) cos α]/π to within 1% at all 15 phase angles, and the full-phase geometric albedo must equal 2/3. Independently, at phase 0, run the 3D mode with a horizontally uniform T-p profile and compare its absolute Fp/Fs against the existing 1D PICASO reflected-light mode for the same profile; the two must agree within numerical tolerance. If either check fails, the factor-of-three Kepler amplitude gap is at least partly a code artifact rather than missing cloud physics, and the central capability claim is unverified.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The paper's central deliverable is the new reflected-light phase-curve mode in PICASO (Section 2.2.1), yet nothing in the manuscript validates the absolute photometry or geometry of that mode. The only benchmark is the Kepler-7b comparison (Section 3.2.1), which misses the observed amplitude by roughly a factor of three: the three-cloud, fsed=0.03 model peaks at ~11-14 ppm versus the observed 50±2 ppm. The paper attributes this gap to missing zonal transport, cloud radiative feedback, and nonspherical particles, but a normalization or integration error in the newly written geometry would produce the same symptom. The Appendix only tests spatial resolution (10x10 versus 20x20 versus 30x30); it demonstrates grid convergence, not absolute correctness. Because this factor-of-three deficit is the largest quantitative discrepancy in the paper, the 'enhanced capability' claim is not yet supported: one cannot distinguish a physical under-prediction from an implementation error in the reflected-light flux normalization or in the Chebyshev-Gauss phase-angle weighting of Equation (6). The issue is not that the Kepler model disagrees with consensus; it is that the new code itself has no known-answer test.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper introduces a new capability in the open-source radiative transfer code PICASO for computing reflected-light phase curves from three-dimensional GCM atmospheres. The routine rotates and regrids GCM temperature, chemistry, and Virga cloud fields onto the visible dayside hemisphere using Chebyshev–Gauss integration. The authors apply it to Kepler-7b using the SPARC/MITgcm thermal structure from Adams et al. (2022), modeling three cloud compositions (MgSiO3-only, Mg2SiO4-only, and Mg2SiO4+Al2O3+TiO2) at fsed = 0.03, 0.1, 0.3, and 1.0, and compare KRF-weighted phase curves to the Demory et al. (2013) Kepler data. The models reproduce the general westward cloud asymmetry but under-predict the observed 50 ± 2 ppm amplitude by roughly a factor of three; the lowest reduced chi-square is 1.94 for the three-cloud, fsed = 0.03 model. The paper concludes that Al2O3 and TiO2 contribute substantially to the albedo and that better treatments of zonal transport, cloud radiative feedback, and particle shape are needed.","tokens_in":21601,"tokens_out":8077,"duration_ms":72506,"significance":"If the code is correct, this is a useful contribution: PICASO and Virga are open source, the authors provide a public tutorial, and the comparison to Kepler-7b is transparent, including the factor-of-three amplitude deficit. The use of a standard Toon et al. (1989) radiative-transfer core and the explicit resolution-convergence test in the Appendix are strengths. The paper also makes falsifiable predictions for future reflected-light observations of hot Jupiters. However, the central new routine is not tested against any known solution, and the only benchmark under-predicts the observed amplitude by a factor of three. The scientific conclusions about Kepler-7b's clouds are therefore conditional; the main value of the paper is the tool description and the transparent case study.","major_comments":[{"comment":"The new reflected-light phase-curve mode is the central deliverable, yet no known-answer test of its absolute normalization or geometry is provided. The Appendix compares 10×10, 20×20, and 30×30 grids and demonstrates self-convergence, but a grid-convergence test cannot catch a global factor error in the flux normalization of Eq. (5) or in the Chebyshev–Gauss weights of Eq. (6). The factor-of-three amplitude deficit against Kepler-7b (Section 3.2.1) is exactly the symptom such an error would produce, so the current paper cannot distinguish an implementation error from the physical explanations offered. Please add a benchmark against an analytic Lambert-sphere phase curve (with a known full-phase flux) or a fixed-phase comparison against the existing 1D reflected-light mode of PICASO, and verify Eq. (5) with the time-dependent planet–star distance and phase-angle convention.","section":"Section 2.2.1 and Appendix"},{"comment":"The claim that the three-cloud, low-fsed models “match best” is based on reduced chi-square values of 1.94–2.75, all well above unity, with every model under-predicting the observed amplitude by roughly a factor of three. Since no model is an acceptable absolute fit, the ranking may reflect the phase-curve shape and the chosen exclusion of 0° ± 9° rather than a robust amplitude constraint. Please report the number of independent data points used, the number of fitted parameters, and ideally confidence intervals on fsed, and state more explicitly that the preference for fsed = 0.03–0.1 is conditional on the assumed GCM and cloud physics.","section":"Section 3.2.1"},{"comment":"The two dominant explanations for the amplitude deficit—missing zonal transport and missing cloud radiative feedback—are not tested in this paper. In particular, the SPARC/MITgcm thermal structure (Section 2.2, Figure 1) does not include cloud radiative feedback, and Virga cloud locations and optical depths depend on that temperature structure. A GCM with radiative feedback could change the dayside temperature field, the cloud map, and the phase curve. Please either quantify the sensitivity of the phase curve to the temperature structure or clearly label the Kepler-7b conclusions as dependent on this no-feedback assumption.","section":"Section 3.2.1"}],"minor_comments":[{"comment":"“RP is the stellar radius” should read “planetary radius”; the equation and Table 1 otherwise imply R_P is the planet radius.","section":"Equation (5) and surrounding text"},{"comment":"“N. J. Kasdin (Kasdin et al. 2020; S. R. Vaughan et al. 2023)” should be “N. J. Kasdin et al. (2020)”.","section":"Section 4"},{"comment":"Clarify how negative Kzz values are treated; the text says the minimum is set to 10^5 cm2 s−1, and the caption says “Low/negative Kzz values at pressures >0.1 bar are set to the Kzz minimum”, which should be stated in the main text.","section":"Section 2.1 and Figure 1 caption"},{"comment":"The sentence describing the 20×20 grid (“grid resolution fluctuates... spans around 167° of longitude for a phase of 0°”) is confusing; it would help to state explicitly that the number of grid points is fixed while the longitudinal extent of the visible dayside shrinks with phase.","section":"Section 2.2.1"},{"comment":"“All our Virga models reproduce the cloudy region west of the substellar point expected from previous studies” should be “produce” or “recover”, since the models are not fitting those previous results.","section":"Abstract and Section 3.2.1"}],"recommendation":"major_revision","confidential_remarks":"The paper is a genuine methods contribution with a transparent case study, and the open-source code release is a strength. The missing known-answer test for the new reflected-light mode is the main obstacle; the factor-of-three amplitude deficit could be physical, but the manuscript does not currently rule out a normalization bug. This is fixable within the scope of the paper by adding an analytic or 1D benchmark. No concerns about citation ethics or novelty disclosure."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The useful thing here is a real new capability: PICASO can now turn GCM output into reflected-light phase curves, tracking the waxing and waning dayside as the planet orbits. The paper is open about what it does, ships code and a tutorial, and the Kepler-7b comparison is transparent, including the honest admission that the models peak at ~11–14 ppm against 50±2 ppm observed. The new cloud maps, with Al2O3 and TiO2 opacity near the substellar point and magnesium silicates at the limbs, are also genuinely new and worth having.\n\nThe soft spot is the one the stress-test flags, and it is real. The new geometry in Section 2.2.1 is the whole point of the paper, yet nothing validates its absolute normalization. The only benchmark is Kepler-7b, which misses by a factor of three. The authors attribute that gap to zonal transport, cloud radiative feedback, and nonspherical particles, and those are plausible. But a normalization error or an integration bug in the Chebyshev–Gauss weighting of Equation (6) would produce exactly the same symptom. The appendix shows grid convergence, not absolute correctness. I would want a known-answer test: a Lambertian sphere phase curve, or a simple 1D radiative-transfer comparison at a few phase angles, before trusting the code's absolute fluxes. Without that, the 'enhanced capability' claim is not fully supported.\n\nTwo smaller issues. First, the fsed conclusion is a post-hoc chi-square selection on the same dataset, so the 'best fit' is a fit result, not a prediction. The paper mostly frames it that way, but the abstract overstates the preference for low fsed. Second, the GCM without cloud radiative feedback and Virga without zonal transport are acknowledged limitations, and they likely explain part of the amplitude gap. Those are fine as stated, but they make the Kepler-7b comparison a demonstration, not a validation.\n\nThe paper deserves a serious referee. The community needs exactly this kind of open-source tool for Roman and HWO planning, and the core geometry change is worth scrutiny. But I would send it back demanding the validation test, and I would not cite the absolute phase-curve amplitudes in my own work until that test is in the literature.","headline":"A genuinely new reflected-light phase-curve mode for PICASO, but the missing known-answer test for the new geometry means the factor-of-three Kepler under-prediction could be a code bug rather than missing physics.","tokens_in":22193,"tokens_out":1476,"would_cite":false,"duration_ms":16918,"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":"This paper claims that PICASO can now model reflected-light phase curves from GCM atmospheres, and that for Kepler-7b a three-cloud, low-sedimentation scenario best matches the observed phase curve but under-predicts its brightness by…","keywords":["reflected light phase curves","PICASO","Virga","Kepler-7b","hot Jupiter","exoplanet clouds","radiative transfer","general circulation models"],"falsifier":"Measure a spatially resolved dayside brightness map of Kepler-7b at near-infrared wavelengths (e.g., with JWST) and compare the longitude and contrast of the bright western region with the model's cloud map; if the observed bright region is substantially broader or brighter than the Virga predictions, the paper's factor-of-three gap is explained, and if not, the missing brightness must come from something else such as zonal transport or non-spherical particles.","tokens_in":21164,"feed_emoji":"☁️","tokens_out":5496,"duration_ms":46013,"temperature":0.7,"pith_summary":"This paper introduces a new routine in the open-source radiative transfer code PICASO that computes reflected-light phase curves directly from three-dimensional GCM atmospheres, with clouds supplied post hoc by the Virga cloud model. The authors apply it to Kepler-7b, a hot Jupiter whose optical phase curve is dominated by starlight reflected off clouds. They find that the best match to the Kepler data comes from a three-condensate scenario (Mg2SiO4, Al2O3, TiO2) with low sedimentation efficiency (fsed = 0.03, reduced chi-square 1.94), but every model under-predicts the observed brightness by about a factor of three. The work matters because reflected-light phase curves will be a primary observable for the next generation of direct-imaging and phase-curve telescopes, and this tool makes such predictions possible from published GCM output.","feed_headline":"PICASO models reflected-light phase curves from 3D climates","feed_subtitle":"Kepler-7b runs match the observed shape with three cloud species but under-predict brightness threefold.","key_machinery":"The load-bearing machinery is the new reflected-phase-curve routine in PICASO, which adapts the thermal phase-curve method by tracking how much of the illuminated dayside is visible at each orbital phase and integrating facet intensities with Chebyshev-Gauss weights. Cloud optical properties (single-scattering albedo, asymmetry parameter, optical depth per pressure level) are computed by Virga at every GCM grid point and rotated/regridded onto a 20 x 20 dayside facet grid before radiative transfer; the choice of 20 x 20 keeps errors below two percent versus 30 x 30.","core_discovery":"The central claim is that PICASO can now take a GCM's temperature-pressure and eddy-diffusion fields, run Virga to compute equilibrium condensate clouds, and produce a full-orbit reflected-light phase curve whose shape and offset are set by the three-dimensional cloud distribution. For Kepler-7b, the model reproduces the expected cloud bank west of the substellar point, high-latitude and eastern-limb magnesium silicate clouds, and Al2O3/TiO2 clouds near the substellar point; the low-fsed three-cloud models match the observed phase-curve shape best, yet yield intensities roughly one-third of the observed values. The paper's proposed explanations for the shortfall are an under-bright or undersized western cloud region, missing zonal transport of cloud particles, neglected cloud radiative feedback, and the assumption of spherical particles.","pith_inferences":["Including cloud radiative feedback in the GCM would likely brighten the dayside and may close part of the factor-of-three gap, since the paper notes feedback is neglected.","Zonal transport of cloud particles from the nightside toward the substellar point is a testable extension: a GCM coupled to a transport-aware cloud microphysics scheme could be compared against the Virga static-cloud prediction.","The factor-of-three shortfall suggests assuming spherical, compact particles underestimates scattering; computing scattering with porous or aggregate particles would give a quantitative check.","The paper's resolution tests show 20 x 20 is a reasonable default; a 10 x 10 grid introduces more than five percent errors at some phases, so future users should not drop below 20 x 20."],"forward_implications":["The community gets a validated, open-source path from GCM output to reflected-light phase curves, applicable to any exoplanet with a 3D climate model.","Kepler-7b's albedo appears to require multiple cloud species: Al2O3 and TiO2 can contribute up to half of the reflected light, despite magnesium silicates dominating the western cloud bank.","Phase-curve offsets are directly readable as cloud-location diagnostics: low-fsed, multi-cloud models place the offset at about 30 degrees westward, matching observations.","Future Roman and Habitable Worlds Observatory reflected-light observations can be interpreted with physical cloud models rather than free-parameter albedo fits.","Thermal phase-curve and reflected-light capabilities now share one code, enabling joint visible and infrared modeling of the same planet."],"supporting_citations":[{"why":"Supplies the SPARC/MITgcm temperature, pressure, and Kzz fields for Kepler-7b that drive all cloud and phase-curve calculations.","marker":"D. J. Adams et al. 2022"},{"why":"Provides the thermal-emission phase-curve routine and Chebyshev-Gauss integration that the reflected-light update modifies.","marker":"N. Robbins-Blanch et al. 2022"},{"why":"Presents the PICASO radiative transfer model and its reflected-light formalism that the new routine builds on.","marker":"N. E. Batalha et al. 2019"},{"why":"Defines the EddySed sedimentation and mixing balance that Virga implements to set cloud vertical extent and particle sizes.","marker":"A. S. Ackerman & M. S. Marley 2001"},{"why":"Provides the Kepler-7b phase-curve observations and best-fit model used as the benchmark for comparison.","marker":"B.-O. Demory et al. 2013"},{"why":"Predicted westward-shifted reflected phase curves from cloud formation west of the substellar point, which this paper's cloud maps reproduce.","marker":"V. Parmentier et al. 2016"},{"why":"Showed cloud radiative feedback affects dayside temperature structure; the paper cites it as a missing process in its GCM.","marker":"M. Roman & E. Rauscher 2017"},{"why":"Supplies the Resampled Opacity Database for PICASO, providing the gas opacities used in the Kepler-band calculations.","marker":"N. Batalha et al. 2020a"},{"why":"Earlier Virga-based cloud models of Kepler-7b with which the new 3D results are compared for consistency.","marker":"M. W. Webber et al. 2015"},{"why":"Microphysical nucleation model predicting which condensates form in significant amounts, motivating the three-cloud scenario.","marker":"D. Powell et al. 2019"}],"fun_headline_variants":["PICASO 3D clouds match Kepler-7b phase shape, but dim","New PICASO phase curves: right clouds, one-third too faint","Kepler-7b: 3D cloud models catch shape, miss brightness","PICASO upgrade adds reflected-light curves from 3D atmospheres","Reflected light phase curves: PICASO with Virga for hot Jupiters"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The GCM thermal structure from Adams et al. (2022), which does not include cloud radiative feedback, is a faithful representation of Kepler-7b's atmosphere, because that temperature field determines where Virga places clouds and therefore the shape and amplitude of the phase curve.","fun_headline_variants_meta":{"raw":{"variants":["PICASO 3D clouds match Kepler-7b phase shape, but dim","New PICASO phase curves: right clouds, one-third too faint","Kepler-7b: 3D cloud models catch shape, miss brightness","PICASO upgrade adds reflected-light curves from 3D atmospheres","Reflected light phase curves: PICASO with Virga for hot Jupiters"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000323,"raw_usage":{"total_tokens":1854,"prompt_tokens":1022,"completion_tokens":832,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":638,"completion_tokens_details":{"reasoning_tokens":726}},"tokens_in":638,"tokens_out":832,"duration_ms":7385,"temperature":1.0,"reasoning_tokens":726,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T15:24:19.764845+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure a spatially resolved dayside brightness map of Kepler-7b at near-infrared wavelengths (e.g., with JWST) and compare the longitude and contrast of the bright western region with the model's cloud map; if the observed bright region is substantially broader or brighter than the Virga predictions, the paper's factor-of-three gap is explained, and if not, the missing brightness must come from something else such as zonal transport or non-spherical particles.","supporting_citations":[],"review_version":1}