{"id":"0e55cf5f-5e6a-452d-af1a-200c04610e26","arxiv_id":"2607.22754","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Adding phase-dependent reflected-light photometry to astrometry-based orbit fits improves inclination precision by 12% (SNR=3) to 33% (SNR=10) in simulated Roman CGI data.","lead":"This paper modifies the orbitize! exoplanet-orbit fitting code so it can combine astrometry with reflected-light photometry, then runs simulated Roman Coronagraph observations. The result: including photometry tightens orbital constraints, improving inclination precision by 12% at SNR=3 and 33% at SNR=10.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Photometric precision gains are demonstrated only under a self-consistent Lambertian model; a realistic non-Lambertian phase function could bias or eliminate the 12–33% improvement.","rationale":"The reader's weakest_assumption—that the Lambertian model is used both to generate and fit the data—is the most load-bearing soft spot. I agree with that assessment. The central claim is specifically about posterior precision, which can improve even with a misspecified likelihood, so the paper's own stated goal ('improves orbital parameter estimates') requires that the photometric model not systematically misrepresent the phase-dependent signal. The authors' self-acknowledged 'simplistic' model is a limitation, not fraud, and the paper is upfront about future work. The missing piece is a misspecification test: generate data with a different, more realistic phase function and fit with the current model. Until that is done, the quantitative 12–33% gains should be treated as conditional on the Lambertian assumption. This is consistent with the reader's CONDITIONAL verdict, so no change to the verdict is needed; the concern strengthens the case for conditional acceptance but does not overturn it.","tokens_in":4087,"tokens_out":6400,"duration_ms":60375,"concrete_test":"Regenerate mock photometry using a non-Lambertian phase function (e.g., a Henyey-Greenstein phase function with g=0.3–0.5, or an independently calibrated cloudy-atmosphere model) for the same orbital parameters and SNR=3, 5, 10, while keeping astrometry generated as in the paper. Fit these data with the current Lambertian model in orbitize!, and compare the joint-fit posterior to the astrometry-only posterior across multiple noise realizations. If the joint fit's posterior width reduction is preserved and the true values remain within the credible intervals, the claim is robust; if the posterior becomes biased or the precision gain shrinks below the quoted 33%/12%, the central claim is an artifact of model self-consistency.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central quantitative claim—joint astrometry+photometry improves orbital posterior precision, with 33% at SNR=10 versus 12% at SNR=3—rests entirely on the photometric likelihood built from the Lambertian reflection model (Appendix Eq. 5). The mock data are generated with the same model used for fitting (§2.1), so the experiment is a self-consistency check, not a test against real reflected-light behavior. The authors explicitly acknowledge this limitation in §3.1: the model is 'simplistic' and will be replaced by more realistic atmosphere models. This matters because if the true phase function differs (e.g., due to clouds, Rayleigh scattering, or specular reflection), the photometric likelihood is misspecified. In a joint fit, a misspecified likelihood can still narrow the posterior while shifting it away from the true orbit, so the stated 'improves orbital parameter estimates' is not established for real Roman data. The precision metric itself is undefined (e.g., standard deviation vs. credible interval width), and results come from a single noise realization, but the model-misspecification issue is more load-bearing because it threatens the external validity of the entire effect, not just its uncertainty.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript describes an extension to the open-source orbitize! package for jointly fitting astrometric and reflected-light photometric observations of directly imaged exoplanets, using a Lambertian disk phase function. The authors generate mock Roman Coronagraph Instrument (CGI) astrometric and photometric measurements at SNRs of 3, 5, and 10, then run parallel-tempered MCMC orbit fits with and without the photometric likelihood. They report that including photometry improves posterior precision, with the effect increasing with SNR: 12% improvement in inclination precision at SNR=3 and 33% at SNR=10. The Appendix derives the Lambertian brightness model in the orbitize! coordinate system. The paper is framed as a technical Note for the orbitize! community, with the main deliverable being the new joint-modeling capability and a demonstration of its potential value for Roman observations.","tokens_in":4314,"tokens_out":4120,"duration_ms":45702,"significance":"If the quantitative result holds, the paper provides a useful, forward-looking capability for planning Roman CGI observations and for orbit fitting of reflected-light exoplanets. The open-source implementation in orbitize! and the explicit derivation of the Lambertian likelihood are strengths; the paper is reproducible in principle and directly extends a widely used community tool. However, the central quantitative claim—that joint fitting improves posterior precision by 12–33%—rests on a single noise realization per SNR, an undefined precision metric, and the assumption that the Lambertian model used to generate the data also describes the true phase function. These issues do not invalidate the tool itself, but they currently limit the strength of the conclusions about real Roman data. With additional robustness tests and clarified metrics, the paper could be a valuable contribution to the direct-imaging and orbit-fitting literature.","major_comments":[{"comment":"The central quantitative claim, e.g., the '33% improvement in inclination precision at SNR=10', is not defined. It is unclear whether 'improvement' refers to the width of the 68% credible interval, the standard deviation, the variance, or some other measure. The abstract and Section 3 present these numbers as headline results, so the metric must be specified. I request a precise definition and a statement of how the improvement is computed from the posterior samples.","section":"Abstract and Section 3"},{"comment":"Each SNR case appears to use a single random noise realization (one mock dataset per SNR). The reported 12% and 33% values are therefore point estimates from one realization and carry no uncertainty. Precision gains could easily vary by several percent across realizations. The authors should repeat the analysis with multiple mock datasets per SNR (e.g., 10–50 random seeds) and report the distribution of improvement, or at least the median and spread. This is load-bearing for the abstract's quantitative claim.","section":"Sections 2.1 and 3"},{"comment":"The mock data are generated using the same Lambertian disk model that is subsequently used in the fit. This is a standard self-consistency test, but it means the reported precision gains are conditional on the Lambertian model being exactly correct. The manuscript itself acknowledges in Section 3.1 that the model is 'simplistic' and will be replaced by more realistic atmosphere models. Given that the abstract states the joint fit 'improves orbital parameter estimates' in general terms, I ask for either (a) a sensitivity test with an alternative non-Lambertian phase function (e.g., a Henyey-Greenstein or empirical phase function) to see whether the precision gain persists, or (b) a clear limitation statement in the abstract and conclusions that the quoted gains are valid only under the Lambertian assumption. Without this, the external validity of the central claim is not established.","section":"Section 2.1 and Appendix Eq. (5)"},{"comment":"MCMC convergence is assessed 'by eye' (Section 2.2). While the chain lengths are generous, the posterior precision numbers in Section 3 are the main result, and visual convergence checks alone are not sufficient for a quantitative claim. Please provide convergence diagnostics (e.g., Gelman-Rubin statistics, integrated autocorrelation time, or acceptance rates) or, at minimum, overplot traces for a few representative parameters. This concern is secondary to the single-realization issue but still affects the reliability of the quoted improvements.","section":"Section 2.2"}],"minor_comments":[{"comment":"The symbol R is reused: Eq. (1) defines R as the instantaneous star–planet distance, while the text below Eq. (5) says R is the planetary radius. This is confusing; please use distinct symbols (e.g., r for separation and R_p for planetary radius) and clarify how the relative brightness amplitude is parameterized in the orbitize! implementation.","section":"Appendix Eq. (5)"},{"comment":"The parameter 'τ 58849 = 0' is not defined. It appears to be the time of periastron passage in MJD, but as written it is unclear. Please define all orbital parameters in the simulation setup.","section":"Section 2.1"},{"comment":"The term 'Lambertain' is misspelled; it should be 'Lambertian' (Sections 2.1, 3.1, and the Appendix). Also, the abstract says 'Nancy Roman Grace Space Telescope'; the correct name is 'Nancy Grace Roman Space Telescope'.","section":"Throughout"},{"comment":"In the text, both 'Blunt et al. 2020a' and 'Blunt et al. 2020b' are cited, but the reference list appears to give the same journal article twice (2020a and 2020b). If these are the same paper, consolidate; if they are distinct, provide the appropriate bibliographic entries.","section":"References"},{"comment":"The bottom six panels in Figure 1 are described as showing posterior constraints with and without photometry, but the caption does not indicate which two parameters are plotted (presumably eccentricity and inclination). Please label the panels or expand the caption so the reader can interpret the figure without referring to the text.","section":"Figure 1"},{"comment":"The derivation would benefit from a sentence defining the physical angles. In particular, α is called 'the angle subtended by the path between the star, the planet, and the observer'—this is the phase angle—and Eq. (4) is the Lambertian phase function. A brief remark connecting β and α to the standard phase-angle notation would make the Appendix clearer.","section":"Appendix Eq. (3) and (4)"}],"recommendation":"major_revision","confidential_remarks":"The paper's open-source contribution and forward-looking simulation are worthwhile. The main issue is that the headline quantitative results are not yet robustly established: the precision metric is undefined, each SNR uses one mock realization, and the generative and fitted photometric models are identical. These are fixable within the scope of the manuscript through additional simulations and clarification. I would not go directly to rejection, but the current version is not ready for publication without these revisions."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know: this is a short software note, not a full paper, and the central effect is real but the quantitative claim is thinner than the abstract makes it sound. The authors add a reflected-light photometric likelihood (Lambertian disk) to orbitize! and run posterior recovery on mock Roman CGI data. The qualitative direction—more data narrows posteriors—is unsurprising; the useful part is the SNR-dependent estimate of how much.\n\nWhat's genuinely new: the specific simulation of Roman CGI partial orbits with SNR-dependent astrometric and photometric errors, and the implementation in an open-source package. The geometry in the appendix is standard. The paper is honest about the Lambertian limitation.\n\nSoft spots: the 12% and 33% numbers come from a single noise realization per SNR. No error bars on those improvements, no definition of which interval width (68%? 95%?) was measured, and convergence was assessed by eye. The mock data are generated from the same model used in the fit, so this is a self-consistency check. The stress-test note is right: if real reflected-light phase functions differ, the likelihood is misspecified and the precision gain could turn into bias. The authors acknowledge this in §3.1, but it means the external validity is not established. Also, no code or data were shipped, and there's a small version inconsistency (Sec. 2.2 says 3.3.0, software list says 3.4.0).\n\nNone of these are fatal. The central argument holds: adding a correct likelihood cannot hurt, and the improvement increasing with SNR is expected. The paper just needs to be more careful about the strength of its claims.\n\nWho it's for: people working on direct imaging orbit fitting, especially with Roman CGI. It's a useful prompt for the community to think about joint fitting. It deserves a serious referee—the methodology is easy to fix and the results, once robust, would be a helpful reference.\n\nRecommendation: send it to peer review, but ask for repeated simulations, a defined metric, better convergence diagnostics, and ideally a release of the code and data.","headline":"A plausible, well-scoped simulation shows joint photometry+astrometry tightens orbital posteriors, but the headline 12–33% numbers rest on one noise draw and a self-consistent Lambertian model.","tokens_in":4832,"tokens_out":1821,"would_cite":false,"duration_ms":18445,"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":"Jointly fitting photometry with astrometry tightens exoplanet orbit posteriors by up to 33%.","keywords":["exoplanet direct imaging","Roman Coronagraph Instrument","reflected-light photometry","orbital parameter estimation","Bayesian orbit fitting","phase function","Lambertian scattering","astrometry"],"falsifier":"Observe a directly imaged exoplanet with Roman at SNR~10 over a partial orbit and compare the inclination posterior from astrometry alone against the joint fit; if the true phase curve deviates from Lambertian in a way that the model cannot absorb, the recovered inclination would be biased or the precision gain would not materialize. A controlled version: generate mock data with a non-Lambertian phase function (e.g., Rayleigh or cloudy scattering) and see whether the Lambertian-based joint fit improves or degrades accuracy.","tokens_in":3912,"feed_emoji":"🔭","tokens_out":4364,"duration_ms":40262,"temperature":0.7,"pith_summary":"The paper tests whether adding reflected-light photometry to astrometric orbit fits improves the recovered orbital parameters for exoplanets imaged by the Roman Coronagraph Instrument. Simulating realistic partial-orbit observations at signal-to-noise ratios of 3, 5, and 10, the authors find that joint fits always tighten the posterior, with inclination precision improving 33% at SNR=10 versus 12% at SNR=3. The gain comes from the phase dependence of reflected light, which breaks degeneracies that astrometry alone leaves. A sympathetic reader would care because Roman is expected to produce the first reflected-light direct imaging of exoplanets, so this is a practical route to better orbits before the Habitable Worlds Observatory.","feed_headline":"Photometry tightens exoplanet orbits by a third","feed_subtitle":"Roman-era joint fits of light and position yield 33% sharper inclinations at high signal-to-noise.","key_machinery":"The load-bearing mechanism is the Lambertian disk reflection model, which gives the planet's brightness as proportional to A α / R², where A is albedo, R is the star-planet distance, and α is a phase-angle function computed from the geometry of the star, planet, and observer. This photometric model is added to an open-source Bayesian orbit-fitting package, so that each orbit hypothesis predicts both an astrometric position and a brightness; the fit then updates the posterior from both data streams. The phase-dependent term is what carries the extra information.","core_discovery":"The central claim is that a joint model of astrometric positions and phase-dependent reflected-light brightness, both generated from the same orbital geometry and fit simultaneously, yields narrower posterior distributions for orbital elements than astrometry alone. The effect is SNR-dependent: at SNR=3 photometry adds 12% improvement in inclination precision, while at SNR=10 it adds 33%. The paper further shows that the improvement appears in eccentricity as well as inclination, and that the added constraint grows as data quality improves.","pith_inferences":["If real directly imaged planets scatter light non-Lambertianly (clouds, Rayleigh scattering, specular highlights), the actual precision gain may differ; a mismatch between model and truth could bias orbital parameters rather than just broaden posteriors.","The 33% figure is for one fiducial system (a=50 AU, e=0.3, i=30°); the gain likely depends on orbital phase coverage and viewing geometry, and future work could map where photometry helps most.","Joint fitting may help break the known degeneracy between inclination and other orbital elements in partial orbits, which is why inclination improves most.","The method assumes Gaussian photon-noise-limited uncertainties; real coronagraphic data include systematics such as speckle and calibration errors that may dilute the gain."],"forward_implications":["For Roman Coronagraph Instrument observations, adding photometry to orbit fits should be standard, yielding noticeably tighter inclinations and eccentricities at high SNR.","Higher-SNR imaging is especially valuable: the fractional gain from photometry grows from 12% to 33% as SNR goes from 3 to 10, so investments in achieving high SNR pay off doubly.","The same joint-fitting approach can be applied to future Habitable Worlds Observatory data, where reflected-light phase curves will be a primary observable."],"fun_headline_variants":["Joint light+position fits sharpen exoplanet orbits","Roman-era joint fits boost orbit precision by 33%","Add photometry to astrometry: 33% sharper inclinations","Combined astrometry and photometry tighten orbital estimates","Modeling light and position together refines exoplanet orbits"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The photometric model assumes a Lambertian disk (matte, uniformly scattering surface) for reflected light, and the same model is used to generate the mock data and to fit them; if real planets reflect light differently, the size of the improvement could change.","fun_headline_variants_meta":{"raw":{"variants":["Joint light+position fits sharpen exoplanet orbits","Roman-era joint fits boost orbit precision by 33%","Add photometry to astrometry: 33% sharper inclinations","Combined astrometry and photometry tighten orbital estimates","Modeling light and position together refines exoplanet orbits"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000344,"raw_usage":{"total_tokens":1683,"prompt_tokens":660,"completion_tokens":1023,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":404,"completion_tokens_details":{"reasoning_tokens":951}},"tokens_in":404,"tokens_out":1023,"duration_ms":7023,"temperature":1.0,"reasoning_tokens":951,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T06:51:29.497056+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Observe a directly imaged exoplanet with Roman at SNR~10 over a partial orbit and compare the inclination posterior from astrometry alone against the joint fit; if the true phase curve deviates from Lambertian in a way that the model cannot absorb, the recovered inclination would be biased or the precision gain would not materialize. A controlled version: generate mock data with a non-Lambertian phase function (e.g., Rayleigh or cloudy scattering) and see whether the Lambertian-based joint fit improves or degrades accuracy.","supporting_citations":[],"review_version":1}