REVIEW 4 major objections 6 minor 35 references
Jointly fitting photometry with astrometry tightens exoplanet orbit posteriors by up to 33%.
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
T0 review · deepseek-v4-flash
2026-08-01 06:51 UTC pith:33U62ZSG
load-bearing objection 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. the 4 major comments →
Jointly Modeling Roman Coronagraph Astrometry and Photometry Improves Orbital Parameter Estimates
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
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.
What carries the argument
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.
Load-bearing premise
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.
What would settle it
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.
If this is right
- 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.
Where Pith is reading between the lines
- 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.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
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.
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 (4)
- [Abstract and Section 3] 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.
- [Sections 2.1 and 3] 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 2.1 and Appendix Eq. (5)] 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 2.2] 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.
minor comments (6)
- [Appendix Eq. (5)] 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 2.1] 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.
- [Throughout] 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'.
- [References] 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.
- [Figure 1] 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.
- [Appendix Eq. (3) and (4)] 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.
Circularity Check
No circular derivation: the joint-fit improvement is demonstrated with injection-recovery simulations under a model the paper explicitly labels as simplistic; this limits external validity but does not make the result circular.
full rationale
The paper's central claim (joint astrometry+photometry improves orbital posterior precision) is established by simulation, not by definition. Mock data are generated with orbitize! using the Lambertian disk reflection model derived in the Appendix (Eqs. 1-5), and the same model is used in the fit. This is a standard self-consistency test: it validates the pipeline and quantifies the information content of photometry under the assumed phase function, but it does not reduce a prediction to its inputs. The photometric likelihood is a physical forward model (brightness ∝ Aα/R²), not a fitted parameter renamed as a prediction. The paper explicitly acknowledges the limitation in §3.1, calling the Lambertain disk reflection model 'simplistic' and stating that future work will replace it with more realistic atmosphere models. Self-citations to orbitize! (Blunt et al. 2020a, 2024) are normal software citations and are not load-bearing: the conclusion rests on the simulations run in this paper, not on an unverified prior result. No uniqueness theorem, ansatz-by-citation, or definitional equivalence is present. Concerns about non-Lambertian phase curves, single noise realization, and undefined precision metric are correctness/robustness issues, not circularity.
Axiom & Free-Parameter Ledger
free parameters (1)
- Relative brightness amplitude (A/R^2) =
not reported
axioms (4)
- domain assumption Lambertian disk reflection model describes the reflected-light phase curve
- domain assumption Noise is photon-noise limited: σ_ast = λ/D/SNR, σ_phot = flux/SNR
- standard math Keplerian orbit model and coordinate system from orbitize! are correct
- ad hoc to paper Mock data generated from the same model used for fitting
read the original abstract
Launching in 2027, the Nancy Roman Grace Space Telescope (Roman) has the potential to directly image exoplanets in reflected light for the first time. Roman imaging will introduce new constraints on exoplanet orbital parameters, since reflected-light intensity depends on orbital phase. In this Note, we discuss an addition to the open-source Python package orbitize!, which allows users to model exoplanet orbits using joint constraints from astrometry and photometric variations due to orbital phase. To investigate the impact of adding photometric data into the orbital model, we simulated realistic measurements of partial orbits, both including and excluding photometry in our model, and computed orbital posteriors. We found that fitting both astrometry and photometry improves posterior precision relative to fitting astrometry alone. This effect was more pronounced for higher-SNR images; for example, photometric data with SNR=10 yielded 33% improvement in inclination precision when including photometry, while SNR=3 data yielded only 12% improvement.
Figures
Reference graph
Works this paper leans on
-
[1]
Roman coronagraph instrument: engineering overview and status. Space Telescopes and Instrumentation 2022: Optical, Infrared, and Millimeter Wave , year = 2022, editor =. doi:10.1117/12.2630540 , adsurl =
-
[2]
Pathways to Discovery in Astronomy and Astrophysics for the 2020s. doi:10.17226/26141 , adsurl =
-
[3]
JWST/MIRI coronagraphic performances as measured on-sky , DOI= "10.1051/0004-6361/202244578", url= "https://doi.org/10.1051/0004-6361/202244578", journal =
-
[4]
2011 , eprint=
Keplerian Orbits and Dynamics of Exoplanets , author=. 2011 , eprint=
2011
-
[7]
Dynamic temperature selection for parallel tempering in Markov chain Monte Carlo simulations. , keywords =. doi:10.1093/mnras/stv2422 , archivePrefix =. 1501.05823 , primaryClass =
-
[8]
First VLTI/GRAVITY Observations of HIP 65426 b: Evidence for a Low or Moderate Orbital Eccentricity. , keywords =. doi:10.3847/1538-3881/ad06b7 , archivePrefix =. 2310.00148 , primaryClass =
-
[10]
American Astronomical Society Meeting Abstracts \#244 , year = 2024, series =
The Habitable Worlds Observatory: Status, Plans, and Opportunities. American Astronomical Society Meeting Abstracts \#244 , year = 2024, series =
2024
-
[11]
Hunter, J. D. , Title =. Computing in Science & Engineering , Volume =
-
[12]
The Nancy Grace Roman Space Telescope Coronagraph Instrument (CGI) technology demonstration. Space Telescopes and Instrumentation 2020: Optical, Infrared, and Millimeter Wave , year = 2020, editor =. doi:10.1117/12.2562997 , archivePrefix =. 2103.01980 , primaryClass =
Pith/arXiv arXiv 2020
-
[13]
Computing in Science and Engineering , keywords =
The NumPy Array: A Structure for Efficient Numerical Computation. Computing in Science and Engineering , keywords =. doi:10.1109/MCSE.2011.37 , archivePrefix =. 1102.1523 , primaryClass =
Pith/arXiv arXiv 2011
-
[14]
Nature Methods , keywords =
SciPy 1.0: fundamental algorithms for scientific computing in Python. Nature Methods , keywords =
-
[15]
2026 , note =
multiprocessing --- Process-based parallelism , howpublished =. 2026 , note =
2026
-
[16]
arXiv , author =:1307.6212 , journal =
doi:10.1051/0004-6361/201322068 , eid =. arXiv , author =:1307.6212 , journal =
-
[17]
The Astropy Project: Building an Open-science Project and Status of the v2.0 Core Package. , keywords =. doi:10.3847/1538-3881/aabc4f , archiveprefix =. 1801.02634 , primaryclass =
-
[18]
The Astropy Project: Sustaining and Growing a Community-oriented Open-source Project and the Latest Major Release (v5.0) of the Core Package. , keywords =. doi:10.3847/1538-4357/ac7c74 , archiveprefix =. 2206.14220 , primaryclass =
-
[19]
IPython: A System for Interactive Scientific Computing. Computing in Science and Engineering , year = "2007", month = "Jan", volume =. doi:10.1109/MCSE.2007.53 , adsurl =
-
[20]
Positioning and Power in Academic Publishing: Players, Agents and Agendas , editor =
Jupyter Notebooks -- a publishing format for reproducible computational workflows , author =. Positioning and Power in Academic Publishing: Players, Agents and Agendas , editor =
-
[21]
doi:10.5281/zenodo.21003741 , url =
The pandas development team , title =. doi:10.5281/zenodo.21003741 , url =
-
[22]
Ralf Gommers and Pauli Virtanen and Matt Haberland and Evgeni Burovski and Tyler Reddy and Warren Weckesser and Andrew Nelson and Travis E. Oliphant and David Cournapeau and Ilhan Polat and alexbrc and Pamphile Roy and Lucas Colley and Pearu Peterson and Josh Wilson and Jake Bowhay and endolith and Nikolay Mayorov and Albert Steppi and Stefan van der Walt...
-
[23]
Adam Ginsburg and Brigitta Sipőcz and C. E. Brasseur and Madhura Parikh and jcsegovia and Austen Groener and Henrik Norman and derdon and Manon Marchand and Michael Kelley and Sam Bianco and Thomas Robitaille and P. L. Lim and Eero Vaher and Christoph Deil and Michael Mommert and Jenny V Medina and Ricky Nilsson and Erik Tollerud and Matthieu Baumann and ...
-
[24]
, archiveprefix = "arXiv", eprint =
astroquery: An Astronomical Web-querying Package in Python. , archiveprefix = "arXiv", eprint =. doi:10.3847/1538-3881/aafc33 , adsurl =
-
[25]
Castillo Hair and Remy Prechelt and Marlon and Kyle Barbary and Johannes Buchner and James Tocknell and James Matthews and NeutralKaon and Ian Czekala , title =
Dan Foreman-Mackey and Adrian Price-Whelan and Will Vousden and Geoffrey Ryan and Matt Pitkin and Víctor Zabalza and jsheyl and Emily Rice and Michael Smith and Leo Singer and Gregory Ashton and Arfon Smith and Kelle Cruz and Zach Weiner and Wolfgang Kerzendorf and Thomas Vandal and Thomas A Caswell and Zhaozhou Li and Stephan Hoyer and Sebastian M. Casti...
-
[26]
The Journal of Open Source Software , year = 2016, month = jun, volume =
corner.py: Scatterplot matrices in Python. The Journal of Open Source Software , year = 2016, month = jun, volume =. doi:10.21105/joss.00024 , adsurl =
-
[27]
S. Behnel and R. Bradshaw and C. Citro and L. Dalcin and D. S. Seljebotn and K. Smith , journal =. Cython: The Best of Both Worlds , year =. doi:10.1109/MCSE.2010.118 , issn =
-
[28]
Dan Foreman-Mackey and Will Meierjurgen Farr and Anne Archibald and Erik Tollerud and David W. Hogg and Andrew Nelson and Nicholas Kern and Jeremy Sanders and Peter Williams and Dustin Lang and Manodeep Sinha and Bence Béky and joezuntz and Adrian Price-Whelan and Eric Dill and Miguel de Val-Borro and Will Vousden and Ilya and Oriol Abril-Pla and Marco Ta...
-
[29]
emcee: The MCMC Hammer. , keywords =. doi:10.1086/670067 , archiveprefix =. 1202.3665 , primaryclass =
-
[30]
Armando Solé and jialin and Daniel Hay Guest , title =
Andrew Collette and Thomas Kluyver and Thomas A Caswell and James Tocknell and Jerome Kieffer and Aleksandar Jelenak and Anthony Scopatz and Darren Dale and Chen and Thomas VINCENT and Matt Einhorn and payno and juliagarriga and Pierlauro Sciarelli and Valentin Valls and Satrajit Ghosh and Ulrik Kofoed Pedersen and Mark Kittisopikul and jakirkham and Mart...
-
[31]
2013 , publisher =
Python and HDF5 , author =. 2013 , publisher =
2013
-
[32]
Sarah Blunt and Jason Wang and Henry Ngo and isabelangelo and Devin Cody and TirthDS and vighneshnagpal and Sofia Covarrubias and Thea McKenna and Tomas Stolker and jorgellop and Robert De Rosa and Mireya Arora and Rodrigo Ferrer-Chávez and Malena Rice and William O Balmer and Lea Hirsch and Adam J. R. W. Smith and Dino Chih-Chun Hsu and eshelDror and Mic...
-
[33]
The Journal of Open Source Software , keywords =
orbitize! v3: Orbit fitting for the High-contrast Imaging Community. The Journal of Open Source Software , keywords =. doi:10.21105/joss.06756 , archiveprefix =. 2409.11573 , primaryclass =
-
[34]
orbitize!: A Comprehensive Orbit-fitting Software Package for the High-contrast Imaging Community. , keywords =. doi:10.3847/1538-3881/ab6663 , archiveprefix =. 1910.01756 , primaryclass =
Pith/arXiv arXiv 1910
-
[35]
Hanno Rein and Dan Tamayo and Shangfei Liu and Lukas Winkler and Pete Bartram and Ari Silburt and Garett Brown and Stanley A. Baronett and Tiger Lu and Robert Melikyan and Rosemary Dorsey and Francesco Biscani and Steven Rieder and Dave Spiegel and Luke Wriglesworth and Jose Antonio Castro and katvolk and Thomas Rometsch and Dang Pham and Albert Zhang and...
-
[36]
REBOUND: an open-source multi-purpose N-body code for collisional dynamics. , keywords =. doi:10.1051/0004-6361/201118085 , archiveprefix =. 1110.4876 , primaryclass =
-
[37]
arXiv e-prints , keywords =
Streamlining and standardizing software citations with The Software Citation Station. arXiv e-prints , keywords =
-
[38]
doi:10.5281/zenodo.17654855 , url =
Tom Wagg and Floor Broekgaarden and Phil Van-Lane and Kai Wu and Kayhan Gültekin , title =. doi:10.5281/zenodo.17654855 , url =
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.