REVIEW 2 major objections 3 minor 47 references
Bayesian Doppler Imaging: Simultaneous Inference of Surface Maps and Geometric Parameters
T0 review · 2 major / 3 minor · reviewed 2026-05-12 · grok-4.3
Pith's one-line read A Bayesian pixel-based framework simultaneously infers surface brightness maps and geometric parameters such as inclination and equatorial rotation velocity from high-resolution spectral time series.
desk verdict The paper's main contribution is a Bayesian framework that jointly infers surface maps and geometric parameters in Doppler imaging via GP priors and HMC, with synthetic recovery of longitudes and a Luhman 16B application that adds uncertainties to prior results. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
Gaussian Process prior over pixel intensities in a Bayesian linear inverse problem, which permits analytical marginalization of map coefficients and Hamiltonian Monte Carlo sampling of nonlinear geometric parameters.
What would settle it
Independent measurements of Luhman 16B's inclination or equatorial rotation velocity falling outside the reported ranges of 48.7 to 75.3 degrees and 28.1 to 36.5 km/s would challenge the joint inference.
Extended reading notes
Core claim
We present a fully Bayesian, pixel-based Doppler imaging framework that enables the simultaneous inference of surface brightness maps and geometric parameters, including the inclination i and equatorial rotation velocity v_rot, from high-resolution spectral time series. We treat the inference as a Bayesian linear inverse problem conditioned on nonlinear geometric parameters. The surface map is modeled as a Gaussian Process prior over pixel intensities, introducing a characteristic spatial scale that sets the map resolution. This allows analytical marginalization of the linear coefficients and efficient sampling of the nonlinear parameters with Hamiltonian Monte Carlo. Validation with synetht
Load-bearing premise
The surface map is modeled as a Gaussian Process prior over pixel intensities that introduces a characteristic spatial scale setting the map resolution, and the method recovers features under the adopted model assumptions including the limited latitudinal sensitivity intrinsic to Doppler imaging.
Editorial extensions
Reading between the lines
- The joint-inference approach could be applied to other rapidly rotating objects to reduce biases that arise when geometry is fixed in advance.
- The noted limited latitudinal sensitivity implies that Doppler imaging alone may always require supplementary data types to resolve features near the poles.
- Public release of the code allows direct tests on new spectral datasets to check whether recovered dark regions persist under varied GP length scales.
- If the derived radius from v_rot and i is combined with evolutionary models for other brown dwarfs, it could tighten constraints on their internal structure.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper presents a fully Bayesian pixel-based Doppler imaging framework for jointly inferring surface brightness maps (via Gaussian Process priors over pixels) and nonlinear geometric parameters (inclination i and equatorial rotation velocity v_rot) from high-resolution spectral time series. Linear map coefficients are analytically marginalized for fixed geometry, with the nonlinear parameters sampled via Hamiltonian Monte Carlo. Synthetic data validation recovers longitudes of inhomogeneities and constrains geometry under the adopted model assumptions (including a single characteristic GP spatial scale and limited latitudinal sensitivity), while the application to VLT/CRIRES observations of Luhman 16B yields a mid-latitude dark region with spatially resolved uncertainties plus i = 61.0_{-12.3}^{+14.3} deg and v_rot = 31.2_{-3.1}^{+5.3} km s^{-1}. The code is made publicly available.
Significance. If the central results hold, the work provides a meaningful advance in Doppler imaging by enabling simultaneous posterior inference of maps and geometry without fixing v sin i or i a priori, along with spatially resolved uncertainty estimates. The synthetic validation and public code are strengths that support reproducibility and falsifiability. The Luhman 16B application demonstrates consistency with prior studies while adding quantitative uncertainties, with potential impact for atmospheric studies of brown dwarfs and directly imaged exoplanets.
major comments (2)
- [§4] §4 (Synthetic Validation): Recovery of v_rot and i is shown only for surface maps generated with the same fixed GP characteristic spatial scale used in the inference; no mismatch tests (e.g., injected maps with different correlation lengths or non-GP structure) are presented. This is load-bearing for the claim that geometry is robustly constrained independently of the prior scale choice, as the forward operator and marginalized likelihood could shift under realistic model mismatch.
- [§5.3] §5.3 (Luhman 16B results): The reported posterior for i (with ~13° uncertainties) is presented as a data-driven constraint, but the manuscript notes intrinsic latitudinal insensitivity of Doppler imaging. Without a prior-only comparison, information-gain metric, or explicit sensitivity test to the GP scale, it remains unclear whether the i posterior is meaningfully informed by the data or largely prior-dominated.
minor comments (3)
- [Abstract and §3] The abstract and §3 should more explicitly state whether the GP characteristic spatial scale is fixed a priori or optimized/marginalized, and how its value was chosen for the Luhman 16B analysis.
- [Figure 4] Figure 4 or equivalent (posterior maps): The uncertainty visualization would be clearer with an additional panel or colorbar showing the ratio of posterior standard deviation to prior standard deviation to highlight data-informed regions.
- [§2] Notation in the forward model (likely Eq. 5-7): The projection and line-of-sight velocity operators could be defined with explicit symbols for the pixel grid and Doppler shift to improve readability for readers outside the subfield.
Simulated Author's Rebuttal
We thank the referee for the constructive review and the recommendation for minor revision. We address each major comment below with proposed changes to the manuscript.
read point-by-point responses
-
Referee: §4 (Synthetic Validation): Recovery of v_rot and i is shown only for surface maps generated with the same fixed GP characteristic spatial scale used in the inference; no mismatch tests (e.g., injected maps with different correlation lengths or non-GP structure) are presented. This is load-bearing for the claim that geometry is robustly constrained independently of the prior scale choice, as the forward operator and marginalized likelihood could shift under realistic model mismatch.
Authors: We agree that mismatch tests would strengthen the validation. However, the manuscript already qualifies all synthetic results as holding 'under the adopted model assumptions,' which explicitly includes the single fixed GP spatial scale. We do not claim that geometry is constrained independently of the prior scale. Performing full mismatch tests with non-GP structures or varied correlation lengths would require substantial additional modeling and computation beyond the scope of this paper. In revision we will expand §4 with a discussion of how the marginalized likelihood depends on the GP hyperparameter and note that exploring model mismatch is an important direction for future work. revision: partial
-
Referee: §5.3 (Luhman 16B results): The reported posterior for i (with ~13° uncertainties) is presented as a data-driven constraint, but the manuscript notes intrinsic latitudinal insensitivity of Doppler imaging. Without a prior-only comparison, information-gain metric, or explicit sensitivity test to the GP scale, it remains unclear whether the i posterior is meaningfully informed by the data or largely prior-dominated.
Authors: The broad uncertainties on i (~13°) already signal the limited latitudinal sensitivity of Doppler imaging, as stated in the text. To clarify the data contribution, we will add to the revised §5.3 a direct overlay of the prior and posterior distributions for both i and v_rot, providing a visual information-gain assessment. We will also include a short sensitivity test showing how the geometric posteriors for Luhman 16B change when the GP characteristic scale is varied around the adopted value. revision: yes
Circularity Check
No significant circularity; derivation self-contained
full rationale
The framework models the surface map via a GP prior with a fixed characteristic spatial scale (a hyperparameter) to enable analytical marginalization of linear pixel coefficients for fixed nonlinear geometry (i, v_rot), followed by HMC sampling. Synthetic validation recovers injected features only under the stated model assumptions, which is the expected behavior for a correctly specified forward model rather than a tautology. Real-data posteriors on i and v_rot are presented as independent outputs from the marginalized likelihood applied to VLT/CRIRES spectra, without any reduction of the target quantities to quantities defined solely in terms of the fitted parameters or prior self-citations. No self-definitional, fitted-input-as-prediction, or load-bearing self-citation steps appear in the derivation chain.
Assumptions & free parameters
free parameters (1)
- characteristic spatial scale
assumptions (1)
- domain assumption The surface map can be modeled as a Gaussian Process prior over pixel intensities
Cite this review
Pith. "Pith review of Bayesian Doppler Imaging: Simultaneous Inference of Surface Maps and Geometric Parameters." pith.science (2026). https://pith.science/paper/2605.08544
@misc{pith2026260508544,
author = {Pith},
title = {Pith review of: Bayesian Doppler Imaging: Simultaneous Inference of Surface Maps and Geometric Parameters},
year = {2026},
howpublished = {\url{https://pith.science/paper/2605.08544}},
note = {Machine review of arXiv:2605.08544}
}
abstract
We present a fully Bayesian, pixel-based Doppler imaging framework that enables the simultaneous inference of surface brightness maps and geometric parameters, including the inclination $i$ and equatorial rotation velocity $v_{\mathrm{rot}}$, from high-resolution spectral time series. We treat the inference as a Bayesian linear inverse problem conditioned on nonlinear geometric parameters. The surface map is modeled as a Gaussian Process prior over pixel intensities, introducing a characteristic spatial scale that sets the map resolution. This allows analytical marginalization of the linear coefficients and efficient sampling of the nonlinear parameters with Hamiltonian Monte Carlo. {Validation with synthetic data demonstrates that our method recovers the longitudes of large-scale surface inhomogeneities and constrains $v_{\mathrm{rot}}$ and $i$ under the adopted model assumptions, while also revealing the limited latitudinal sensitivity intrinsic to Doppler imaging.} We applied this framework to high-resolution VLT/CRIRES observations of the brown dwarf Luhman 16B. Our analysis reveals a large-scale dark region at mid-latitudes, consistent with previous studies but now with spatially resolved uncertainty estimates. Furthermore, we successfully constrained the geometric parameters without fixing \(v_{\mathrm{rot}}\sin i\) or $i$ to literature values, deriving an inclination of $i = 61.0_{-12.3}^{+14.3}$ degrees and an equatorial rotation velocity of $v_{\mathrm{rot}} = 31.2_{-3.1}^{+5.3}~\mathrm{km\,s^{-1}}$. These results indicate a radius broadly consistent with evolutionary models and suggest a possible spin-axis misalignment under the assumption of comparable equatorial rotation velocities for the two components. Our code is publicly available under the MIT license.
Figures
Figures from the paper (10 more)
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.lean (J-cost uniqueness, Aczél classification)washburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
The surface map is modeled as a Gaussian Process prior over pixel intensities, introducing a characteristic spatial scale that sets the map resolution. This allows analytical marginalization of the linear coefficients and efficient sampling of the nonlinear parameters with Hamiltonian Monte Carlo.
-
IndisputableMonolith/Foundation/AlphaCoordinateFixation.lean (higher-derivative calibration of CostAlphaLog)costAlphaLog_high_calibrated_iff unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
Σ_a[j,j′] = σ_a² exp(−d_jj′²/(2ℓ²))
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
Agol, E., Cowan, N. B., Knutson, H. A., et al. 2010, ApJ, 721, 1861, doi: 10.1088/0004-637X/721/2/1861
-
[2]
2024, MNRAS, 532, 1361, doi: 10.1093/mnras/stae1549
Aizawa, M., Muto, T., & Momose, M. 2024, MNRAS, 532, 1361, doi: 10.1093/mnras/stae1549
-
[3]
Apai, D., Nardiello, D., & Bedin, L. R. 2021, ApJ, 906, 64, doi: 10.3847/1538-4357/abcb97 Artigau, ´E. 2018, in Handbook of Exoplanets, ed. H. J. Deeg & J. A. Belmonte (Springer International Publishing AG), 94, doi: 10.1007/978-3-319-55333-7 94 Asensio Ramos, A., D´ ıaz Baso, C. J., & Kochukhov, O. 2022, A&A, 658, A162, doi: 10.1051/0004-6361/202142027
-
[4]
Hauschildt, P. H. 2003, A&A, 402, 701, doi: 10.1051/0004-6361:20030252
-
[5]
Betancourt, M. 2017, arXiv e-prints, arXiv:1701.02434. https://arxiv.org/abs/1701.02434
work page Pith review arXiv 2017
-
[6]
2017, The Astronomical Review, 13, 1, doi: 10.1080/21672857.2017.1303105
Biller, B. 2017, The Astronomical Review, 13, 1, doi: 10.1080/21672857.2017.1303105
-
[7]
2018, JAX: composable transformations of Python+NumPy programs, 0.2.5
Bradbury, J., Frostig, R., Hawkins, P., et al. 2018, JAX: composable transformations of Python+NumPy programs, 0.2.5. http://github.com/google/jax
work page 2018
-
[8]
C., Weiner Mansfield, M., Cubillos, P
Challener, R. C., Weiner Mansfield, M., Cubillos, P. E., et al. 2025, Nature Astronomy, 9, 1821, doi: 10.1038/s41550-025-02666-9 22Ureshino et al
Show all 47 references
-
[9]
A., Vos, J
Chen, X., Biller, B. A., Vos, J. M., et al. 2024, Monthly Notices of the Royal Astronomical Society, 533, 3114–3143, doi: 10.1093/mnras/stae1995
2024 doi
-
[10]
Collin, G. H. 2024, healjax: Port of healpix routines to JAX, 0.2.2. https://github.com/ghcollin/healjax
2024
- [11]
-
[12]
Crossfield, I. J. M. 2014, A&A, 566, A130, doi: 10.1051/0004-6361/201423750
2014 doi
-
[13]
Crossfield, I. J. M., Biller, B., Schlieder, J. E., et al. 2014, Nature, 505, 654, doi: 10.1038/nature12955 de Wit, J., Gillon, M., Demory, B.-O., & Seager, S. 2012, A&A, 548, A128, doi: 10.1051/0004-6361/201219060
2014 doi
-
[14]
Deutsch, A. J. 1958, in IAU Symposium, Vol. 6, Electromagnetic Phenomena in Cosmical Physics, ed. B. Lehnert, 209
1958
-
[15]
D., Pendleton, B
Duane, S., Kennedy, A. D., Pendleton, B. J., & Roweth, D. 1987, Physics Letters B, 195, 216, doi: 10.1016/0370-2693(87)91197-X
1987 doi
- [16]
-
[17]
2018, AJ, 156, 146, doi: 10.3847/1538-3881/aad775
Robinson, T. 2018, AJ, 156, 146, doi: 10.3847/1538-3881/aad775
2018 doi
-
[18]
2016, Journal of Open Source Software, 1, 24, doi: 10.21105/joss.00024
Foreman-Mackey, D. 2016, Journal of Open Source Software, 1, 24, doi: 10.21105/joss.00024
2016 doi
-
[19]
2012, ApJ, 755, 101, doi: 10.1088/0004-637X/755/2/101
Fujii, Y., & Kawahara, H. 2012, ApJ, 755, 101, doi: 10.1088/0004-637X/755/2/101
2012 doi
-
[20]
Gillon, M., Triaud, A. H. M. J., Jehin, E., et al. 2013, A&A, 555, L5, doi: 10.1051/0004-6361/201321620
2013 doi
-
[21]
Yagola, A. G. 1982, Soviet Ast., 26, 690
1982
-
[22]
Yagola, A. G. 1977, Soviet Astronomy Letters, 3, 147 G´ orski, K. M., Hivon, E., Banday, A. J., et al. 2005, ApJ, 622, 759, doi: 10.1086/427976
1977 doi
-
[23]
J., Challener, R
Hammond, M., Bell, T. J., Challener, R. C., et al. 2024, AJ, 168, 4, doi: 10.3847/1538-3881/ad434d
2024 doi
-
[24]
D., & Gelman, A
Hoffman, M. D., & Gelman, A. 2011, arXiv e-prints, arXiv:1111.4246. https://arxiv.org/abs/1111.4246
2011
-
[25]
2010, ApJ, 720, 1333, doi: 10.1088/0004-637X/720/2/1333 —
Kawahara, H., & Fujii, Y. 2010, ApJ, 720, 1333, doi: 10.1088/0004-637X/720/2/1333 —. 2011, ApJL, 739, L62, doi: 10.1088/2041-8205/739/2/L62
2010 doi
-
[26]
2022, The Astrophysical Journal Supplement Series, 258, 31, doi: 10.3847/1538-4365/ac3b4d
Kawahara, H., Kawashima, Y., Masuda, K., et al. 2022, The Astrophysical Journal Supplement Series, 258, 31, doi: 10.3847/1538-4365/ac3b4d
2022 doi
-
[27]
2020, ApJ, 900, 48, doi: 10.3847/1538-4357/aba95e
Kawahara, H., & Masuda, K. 2020, ApJ, 900, 48, doi: 10.3847/1538-4357/aba95e
2020 doi
-
[28]
2025, The Astrophysical Journal, 985, 263, doi: 10.3847/1538-4357/adcba2
Kawahara, H., Kawashima, Y., Tada, S., et al. 2025, The Astrophysical Journal, 985, 263, doi: 10.3847/1538-4357/adcba2
2025 doi
-
[29]
2025, ApJ, 988, 53, doi: 10.3847/1538-4357/adddbd
Kawashima, Y., Kawahara, H., Kasagi, Y., et al. 2025, ApJ, 988, 53, doi: 10.3847/1538-4357/adddbd
2025 doi
-
[30]
Khokhlova, V. L. 1976, Astronomische Nachrichten, 297, 203, doi: 10.1002/asna.19762970408
1976 doi
-
[31]
Kipping, D. M. 2013, MNRAS, 435, 2152, doi: 10.1093/mnras/stt1435
2013 doi
-
[32]
A., Charbonneau, D., Allen, L
Knutson, H. A., Charbonneau, D., Allen, L. E., et al. 2007, Nature, 447, 183, doi: 10.1038/nature05782
2007 doi
-
[33]
R., Parmentier, V., et al
Kreidberg, L., Line, M. R., Parmentier, V., et al. 2018, AJ, 156, 17, doi: 10.3847/1538-3881/aac3df
2018 doi
-
[34]
2021, arXiv e-prints, arXiv:2110.06271, doi: 10.48550/arXiv.2110.06271
Luger, R., Bedell, M., Foreman-Mackey, D., et al. 2021, arXiv e-prints, arXiv:2110.06271, doi: 10.48550/arXiv.2110.06271
2021 doi
-
[35]
Majeau, C., Agol, E., & Cowan, N. B. 2012, ApJL, 747, L20, doi: 10.1088/2041-8205/747/2/L20 Pall´ e, E., Biazzo, K., Bolmont, E., et al. 2025, Experimental Astronomy, 59, 29, doi: 10.1007/s10686-025-10000-4
2012 doi
-
[37]
2019, arXiv preprint arXiv:1912.11554
Phan, D., Pradhan, N., & Jankowiak, M. 2019, arXiv preprint arXiv:1912.11554
2019
-
[38]
E., Tuominen, I., & Vilhu, O
Piskunov, N. E., Tuominen, I., & Vilhu, O. 1990, A&A, 230, 363
1990
-
[39]
K., & Wang, J
Plummer, M. K., & Wang, J. 2022, The Astrophysical Journal, 933, 163, doi: 10.3847/1538-4357/ac75b9 —. 2023, The Astrophysical Journal, 951, 101, doi: 10.3847/1538-4357/accd5d
2022 doi
-
[40]
Snellen, I. A. G., Brandl, B. R., de Kok, R. J., et al. 2014, Nature, 509, 63, doi: 10.1038/nature13253
2014 doi
-
[41]
S., et al
Tan, X., Zhang, X., Marley, M. S., et al. 2025, Science Advances, 11, doi: 10.1126/sciadv.adv3324
2025 doi
-
[42]
2005, Inverse Problem Theory and Methods for Model Parameter Estimation (the Society for Industrial and Applied Mathematics)
Tarantola, A. 2005, Inverse Problem Theory and Methods for Model Parameter Estimation (the Society for Industrial and Applied Mathematics)
2005
-
[43]
2026, BayesianDI, 1.0.0, Zenodo, doi: 10.5281/zenodo.20024778
Ureshino, Y. 2026, BayesianDI, 1.0.0, Zenodo, doi: 10.5281/zenodo.20024778
2026 doi
-
[44]
S., & Penrod, G
Vogt, S. S., & Penrod, G. D. 1983, Publications of the Astronomical Society of the Pacific, 95, 565, doi: 10.1086/131208
1983 doi
-
[45]
S., Penrod, G
Vogt, S. S., Penrod, G. D., & Hatzes, A. P. 1987, ApJ, 321, 496, doi: 10.1086/165647
1987 doi
-
[46]
W., M´ erand, A., Thompson, W., et al
Xuan, J. W., M´ erand, A., Thompson, W., et al. 2024, Nature, 634, 1070, doi: 10.1038/s41586-024-08064-x
2024 doi
-
[47]
2025, arXiv e-prints, arXiv:2511.23018, doi: 10.48550/arXiv.2511.23018 Bayesian Doppler Imaging23
Yama, H., Masuda, K., Kawashima, Y., & Kawahara, H. 2025, arXiv e-prints, arXiv:2511.23018, doi: 10.48550/arXiv.2511.23018 Bayesian Doppler Imaging23
2025 doi
-
[48]
2019, Journal of Open Source Software, 4, 1298, doi: 10.21105/joss.01298
Zonca, A., Singer, L., Lenz, D., et al. 2019, Journal of Open Source Software, 4, 1298, doi: 10.21105/joss.01298
2019 doi
Reviewed May 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.