{"id":"97174c94-4e30-464e-8e13-535a9d4d5b30","arxiv_id":"2605.08544","paper_version":1,"verdict":"CONDITIONAL","confidence":"LOW","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"A fully Bayesian pixel-based Doppler imaging framework uses Gaussian Process priors and Hamiltonian Monte Carlo to simultaneously infer surface maps and geometric parameters from spectral data.","lead":"This paper introduces a Bayesian Doppler imaging method that jointly infers surface brightness maps and geometric parameters like inclination and rotation velocity from high-resolution spectral time series. A smart generalist might read it because the approach provides uncertainty estimates on both maps and geometry without fixing parameters to external values, which could improve studies of brown dwarfs and similar rotating objects.","discovery_kind":"new_method","skeptic_critique":{"model":"grok-4.3","headline":"Single characteristic scale in GP surface prior may bias joint inference of i and v_rot under model mismatch","rationale":"This sharpens the reader's identified weakest assumption (GP prior and model assumptions) into a concrete, testable risk for the joint geometry inference that is central to the strongest claim. It does not warrant REJECT because the method remains internally consistent under its stated assumptions and the real-data results include uncertainties; the low reader confidence is appropriate pending such a robustness check.","tokens_in":1895,"tokens_out":351,"duration_ms":73374,"concrete_test":"Generate synthetic time-series spectra from a surface map containing both a broad mid-latitude feature and smaller-scale spots (different from the single GP lengthscale used at inference time), forward-model at known true i and v_rot with realistic noise; re-run the full Bayesian pipeline across a grid of GP lengthscales and check whether the recovered posterior medians for i or v_rot deviate from truth by more than the reported 1-sigma uncertainties.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The framework treats the surface map as a GP prior over pixel intensities with one characteristic spatial scale that sets resolution and enables analytical marginalization of the linear map coefficients for fixed nonlinear geometry (i, v_rot). Synthetic validation recovers longitudes and constrains geometry only under the adopted model assumptions, including this scale choice and the intrinsic latitudinal insensitivity of Doppler imaging. For the Luhman 16B application, if the observed mid-latitude dark region contains structure at scales differing from the fixed hyperparameter, the effective forward operator and marginalized likelihood could systematically shift the HMC posterior for the geometric parameters away from truth, even while producing plausible uncertainties.","agreement_with_reader":"partial"},"referee_report":{"model":"grok-4.3","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.","tokens_in":2031,"tokens_out":731,"duration_ms":42769,"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":[{"comment":"§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.","section":"§4"},{"comment":"§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.","section":"§5.3"}],"minor_comments":[{"comment":"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.","section":"Abstract and §3"},{"comment":"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.","section":"Figure 4"},{"comment":"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.","section":"§2"}],"recommendation":"minor_revision","confidential_remarks":"The manuscript fits well within the scope of an astronomy journal such as ApJ or MNRAS. No obvious citation or novelty issues; the public code is a positive for the field."},"author_rebuttal":{"model":"grok-4.3","summary":"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.","responses":[{"response":"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_made":"partial","referee_comment":"§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."},{"response":"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_made":"yes","referee_comment":"§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."}],"tokens_in":1590,"tokens_out":483,"duration_ms":52796,"standing_objections":[]},"desk_editor":{"model":"grok-4.3","letter":"The paper introduces a fully Bayesian pixel-based Doppler imaging setup that treats the surface map as a linear inverse problem conditioned on nonlinear geometry. A Gaussian Process prior on pixel intensities allows analytical marginalization of the map coefficients, after which Hamiltonian Monte Carlo samples inclination and equatorial rotation velocity. Synthetic tests recover the longitudes of large-scale features and place constraints on i and v_rot under the model assumptions, while the Luhman 16B analysis yields a mid-latitude dark region with spatially resolved uncertainties plus i = 61.0 with asymmetric errors and v_rot = 31.2 km/s without fixing v sin i to literature values. The code is released publicly under MIT license, which is a clear plus for reproducibility. The single characteristic spatial scale in the GP prior is the clearest soft spot. If the true surface structure has power at scales different from the chosen hyperparameter, the effective forward model and marginalized likelihood can shift the posterior on the geometric parameters even while producing plausible error bars. The abstract already flags the limited latitudinal sensitivity, which is intrinsic to Doppler imaging and means the map is better constrained in longitude than latitude. Validation remains conditional on the adopted assumptions, and full details on data reduction would be needed to judge robustness. This work is aimed at researchers doing high-resolution spectroscopy of brown dwarfs and directly imaged planets who want to avoid fixing geometry when mapping surfaces. Readers interested in Bayesian inverse problems or Gaussian Process methods in astronomy will get the most out of it. It deserves peer review because the joint-inference framework is new, the synthetic tests are informative, and the real-data example shows the method can be applied without external priors on i or v_rot.","headline":"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.","tokens_in":2543,"tokens_out":425,"would_cite":true,"duration_ms":23955,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.3","headline":"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.","keywords":["Bayesian Doppler imaging","surface mapping","geometric parameters","brown dwarf","Luhman 16B","Gaussian process","spectral time series","Hamiltonian Monte Carlo"],"falsifier":"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.","tokens_in":2784,"feed_emoji":"🔭","tokens_out":686,"duration_ms":32499,"temperature":0.7,"pith_summary":"The paper develops a fully Bayesian Doppler imaging method that models the surface as a pixel grid with a Gaussian process prior on intensities. Geometric parameters including inclination and rotation speed are treated as nonlinear variables sampled by Hamiltonian Monte Carlo, while linear map coefficients are marginalized analytically. Validation on synthetic data shows the approach recovers longitudes of large-scale features and constrains the geometry parameters, and the method is then applied to VLT/CRIRES observations of the brown dwarf Luhman 16B. A reader would care because the joint inference supplies uncertainty estimates on both the map and the geometry without relying on fixed literature values for rotation or inclination.","feed_headline":"Bayesian method jointly infers brown dwarf map and geometry","feed_subtitle":"Applied to Luhman 16B it recovers a mid-latitude dark region while constraining inclination to 61 degrees and rotation to 31 km/s with full-","key_machinery":"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.","core_discovery":"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","pith_inferences":["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."],"forward_implications":[],"fun_headline_variants":["Bayesian model jointly infers brown dwarf surface map and geometry","Doppler imaging with Bayesian priors constrains rotation and inclination","Bayesian framework maps dark region on Luhman 16B and parameters","Pixel-based Bayesian method infers inclination and equatorial velocity"],"cache_read_input_tokens":2112,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Bayesian model jointly infers brown dwarf surface map and geometry","Doppler imaging with Bayesian priors constrains rotation and inclination","Bayesian framework maps dark region on Luhman 16B and parameters","Pixel-based Bayesian method infers inclination and equatorial velocity"]},"model":"grok-4.3","cost_usd":0.00948,"raw_usage":{"total_tokens":4300,"prompt_tokens":801,"num_sources_used":0,"completion_tokens":60,"cost_in_usd_ticks":94799500,"prompt_tokens_details":{"text_tokens":801,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":3439,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":801,"tokens_out":60,"duration_ms":40707,"temperature":1.0,"reasoning_tokens":3439,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-05-12T01:16:17.962140+00:00","model_set":{"reader":"grok-4.3"},"falsifier":"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.","supporting_citations":[],"review_version":1}