{"id":"51410deb-94d7-46a6-9d16-7c2337642cff","arxiv_id":"1908.00991","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"Precise radial velocity measurements need a photon-weighted barycentric correction; using the photon-weighted midpoint time instead causes systematic errors up to about one meter per second for long exposures.","lead":"This paper shows that when astronomers measure a star's tiny wobble, they should correct for Earth's motion using the arrival time of every photon, not just the average time of the exposure. Ignoring this detail can create fake velocity shifts of up to one meter per second, enough to hide or mimic small exoplanets.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No significant objection identified","rationale":"The paper's central assertion is that applying the barycentric correction at the photon-weighted midpoint rather than photon-weighting the correction itself produces a second-order RV error of order 10 cm/s for 30-60 minute exposures and more than 1 m/s in realistic worst-case geometries. For this to hold, Eq. (9) must be a faithful estimate of the difference between v(⟨t⟩) and the photon-weighted average of v(t), and the omitted non-diurnal terms must be negligible over an exposure. Both conditions are met. The Taylor expansion is algebraically sound: expanding around ⟨t⟩ removes the linear term exactly, so only the second derivative and the variance of the arrival-time distribution matter. The uniform-flux variance gives Δt²/12, producing the 1.32 m/s per hour² coefficient, which matches the full barycorrpy simulations in Fig. 3 to the displayed precision. The analytical derivation also correctly predicts the 1.5× amplification for a centred V-shaped flux curve. The EXPRES comparison in Fig. 7 demonstrates that the uniform-flux approximation is realistic for short exposures, which supports the correction term proposed for archive data. The weakest point is that Fig. 7 does not compare against independently extracted RV residuals, and no exposure longer than 20 minutes is used; however, the effect under test is a mathematical property of how the barycentric correction enters the measured velocity, not an empirical effect requiring a separate detection. The known non-diurnal components of the barycentric velocity have negligible curvature on hourly timescales: the orbital acceleration is nearly constant over 1 hour, and its jerk contributes well below 1 cm/s. Thus the 60-minute predictions are a benign extrapolation of a validated model rather than a fragile assumption. I therefore agree with the reader's ACCEPT verdict and do not see a reason to condition or reject it.","tokens_in":11107,"tokens_out":15156,"duration_ms":169344,"concrete_test":"Re-run the Section 4 barycorrpy calculation for a 60-minute uniform exposure at Mauna Kea on a second epoch (e.g., 2019-01-15) and a non-equatorial target (e.g., δ=+45°), and compare the simulated v(⟨t⟩)−⟨v⟩ with Eq. (10). Agreement to within ~1 cm/s would close the single-epoch validation gap; separately compute the orbital jerk contribution to Eq. (9) to confirm it remains below 1 cm/s.","verdict_should_be":"UNCHANGED","load_bearing_attack":"After checking the derivation and simulations, I do not find a load-bearing flaw. The central claim rests on Eq. (9), a second-order Taylor expansion of the instantaneous barycentric velocity around the photon-weighted mean time; the first-order term vanishes by the definition of ⟨t⟩, leaving only the curvature term weighted by the variance of the arrival-time distribution. The numerical coefficient in Eq. (10) is reproduced by the barycorrpy simulations in Figs. 3-5 (the 60-min uniform case reaches ±1.00 m/s as stated), and the EXPRES data confirm that the uniform-flux approximation is adequate for real exposure-meter shapes at exposure times where the effect is measurable. The only caveat is that no real long-exposure (≥40 min) spectrum is used to verify that the predicted second-order term actually appears in extracted RVs, and Fig. 7 is a comparison of two calculations of the correction rather than an end-to-end RV test; this is a validation gap, not a reason to doubt the analytic argument. The diurnal-sinusoid model is the dominant curvature term; the orbital barycentric acceleration is nearly constant over an hour and its curvature (jerk) contribution is below the 1 cm/s level. The 60-minute extrapolation is therefore well supported.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript addresses a systematic error in barycentric corrections for radial-velocity (RV) exposures of finite duration. The authors show that the common practice of evaluating the barycentric correction at the photon-weighted midpoint time of an exposure is exact only if the correction varies linearly with time; because the dominant diurnal component is curved, a residual second-order offset remains between the correction evaluated at the photon-weighted mean time and the flux-weighted mean of the instantaneous corrections. They derive analytic estimates (Eqs. 8-10) from a Taylor expansion of a sinusoidal diurnal model, predict errors of 0.25 m/s for a 30 minute and 1.0 m/s for a 60 minute uniform exposure in a worst-case geometry, verify these scalings with barycorrpy simulations for different exposure times, declinations, flux shapes, and observatory latitudes, and use 1315 EXPRES exposures to show that a uniform-flux approximation captures the second-order error to within about 1 cm/s for exposures shorter than 20 minutes. The paper recommends that instruments record and archive exposure-meter flux curves and, for existing data without flux curves, apply a uniform-flux correction term.","tokens_in":11292,"tokens_out":13853,"duration_ms":148826,"significance":"If correct, this is an important and directly actionable result for the precision-RV community. The central derivation is a parameter-free analytic Taylor expansion with no fitted free parameters, and the resulting coefficient in Eq. (10) and the V-shape enhancement factor of 1.5 are concrete, falsifiable predictions. The analytic results are corroborated by simulations with barycorrpy, a widely used and externally benchmarked package, and by a real-data comparison with EXPRES exposure-meter data. The paper closes a gap left open in earlier barycentric-correction work and gives a practical recommendation for both current and future instruments, which is relevant for the 10 cm/s and 1 m/s RV precision regimes.","major_comments":[],"minor_comments":[{"comment":"The text reports 1315 EXPRES observations, while the Fig. 7 caption states 1316 observations; please reconcile the count.","section":"Section 5.2 / Fig. 7"},{"comment":"The statement that \"There is no indication that longer exposure times would change this picture\" is stronger than the data support, because Fig. 7 contains no exposure longer than 20 minutes; I suggest softening the wording or adding a simulation for longer exposures.","section":"Section 5.2"},{"comment":"Equation (10) is obtained from Eq. (9) via the spherical-trig identity cos(δ)sin(ψ) = -cos(alt)sin(az), but the derivation is relegated to a footnote; stating the identity explicitly in the text would make the derivation easier to verify and reproduce.","section":"Section 3"},{"comment":"The sentence \"If photons are concentrated towards one end of the exposure, the error second-order error decreases\" contains a typo; please also state explicitly the sign convention for azimuth so that the signs in Figs. 3 and 6 can be interpreted unambiguously.","section":"Section 3"},{"comment":"The claim that \"In our experience, a typical offset is of order 5% of the exposure time\" is anecdotal and unreferenced; providing a quantitative justification or a reference would strengthen this part of the mitigation discussion.","section":"Section 5.2"},{"comment":"Footnote 2 describes a step-function mitigation strategy for data without flux curves, but it does not specify exactly how to implement it; a brief formula would make the recommendation reproducible.","section":"Section 5.2"},{"comment":"The derivation in Eq. (4) assumes that non-diurnal contributions to the barycentric velocity have negligible curvature over an exposure, and although the barycorrpy simulations support this, the paper would be more self-contained with an explicit order-of-magnitude bound on the neglected terms.","section":"Section 3"}],"recommendation":"minor_revision","confidential_remarks":"The paper is well within the scope of an astronomical techniques journal, and I found no basis for concern about the use of barycorrpy, which is co-authored by two of the authors; it is a standard package with external validation, and the analytic derivation is independent of the code. The main limitation is the absence of direct long-exposure real-data validation, but this is a validation gap rather than a flaw in the central argument, and it can be addressed by softening the relevant wording and adding a brief discussion."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a solid methods paper. It formalizes something the field already knew existed—Fischer et al. (2016) and Blackman et al. (2017) flagged it—but it does the actual work: a clean Taylor-expansion derivation, explicit error formulas (Eqs. 8–10), barycorrpy simulations that reproduce the analytic curves, and a real EXPRES exposure-meter comparison. The 1.32 m/s coefficient in Eq. (10) and the 2.0 m/s-per-minute first-order coefficient are directly checkable from the stated sinusoid model, and they match.\n\nWhat's genuinely new: the analytic error estimates, including the V-shape amplification factor of 1.5, and the practical recommendation to store exposure meter flux curves so the correction can be photon-weighted rather than evaluated at the photon-weighted midpoint. The constant-flux archival correction is a reasonable, clearly-labeled approximation. The paper is honest about the regime where it matters: 30–60 minute exposures near the equator, east/west, where the error reaches 0.25–1.0 m/s.\n\nSoft spots, in proportion: they're minor. The real-data validation only covers exposures under 20 minutes, where the effect is a few cm/s, so the 60-minute magnitudes rest on the analytic model and simulation, not on an end-to-end long-exposure RV test. That's a validation gap, not a flaw—the model is well grounded. There's also a trivial 1315/1316 inconsistency in Fig. 7. No code or data is shipped, which is slightly annoying for a paper whose whole point is \"store your flux curves,\" but the equations are simple enough to reimplement in a few lines and the barycorrpy comparison anchors them. The citation pattern is fine; barycorrpy is co-authored by two of the authors but it's a standard, externally benchmarked package, and the prior work is properly credited.\n\nWho this is for: anyone building or operating a PRV instrument heading toward the 10 cm/s regime, and anyone reducing archival RV data from long exposures. It deserves a serious referee; it's a genuinely useful contribution with a correct central argument. I'd accept it. I'd also bring it to reading group—it's a nice example of a small systematic error being properly quantified.","headline":"A clean, useful derivation of the second-order barycentric correction error that deserves referee time; the central claim holds up.","tokens_in":723,"tokens_out":1060,"would_cite":true,"duration_ms":20292,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The standard barycentric correction for long exposures leaves a systematic radial-velocity error, sometimes over 1 m/s, because the correction curves with Earth's rotation rather than changing linearly.","keywords":["barycentric correction","photon-weighted midpoint time","exposure meter","radial velocity precision","second-order error","Earth rotation","exoplanet radial velocities"],"falsifier":"Take a 60-minute exposure at a low-latitude site on a target with declination near zero observed through the east, record the exposure-meter flux curve, and compare $\\langle v\\rangle$ computed from the full flux curve with $v(\\langle t\\rangle)$ computed at the photon-weighted midpoint; the predicted difference is about 1 m/s (scaled by the actual flux shape). A null result at the 10 cm/s level would show the diurnal-curvature model is incomplete.","tokens_in":10907,"feed_emoji":"🌍","tokens_out":8272,"duration_ms":75703,"temperature":0.7,"pith_summary":"The paper is trying to establish that the standard way of applying the barycentric correction to a long spectroscopic exposure—evaluating the correction once at the photon-weighted midpoint time—leaves a systematic second-order error because the correction curves over the exposure rather than changing linearly. That error is typically about 10 cm/s for 30–60 minute exposures but can pass 1 m/s for realistic worst-case targets observed low in the east or west near the celestial equator. The correct procedure, the paper argues, is to photon-weight the correction itself, using the exposure meter flux curve $f(t)$ as the weight in $\\langle v\\rangle = \\int v(t)f(t)\\,dt/\\int f(t)\\,dt$. The paper therefore recommends that precision radial velocity instruments record and archive the full exposure meter flux curve, and shows that archival data without the curve can be partly repaired with a constant-flux correction term.","feed_headline":"Photon-weight the barycentric correction, or lose 1 m/s","feed_subtitle":"A 60-minute exposure at a low-latitude site can carry a systematic second-order velocity error of a meter per second.","key_machinery":"The central object is the photon-weighted average in Equation (3), with the exposure meter flux curve $f(t)$ acting as the smoothing kernel applied to the instantaneous barycentric correction. Expanding $v(t)$ around $\\langle t\\rangle$ to second order, the missed term is proportional to the time variance of the photon arrival distribution, $\\langle t^2\\rangle-\\langle t\\rangle^2$, multiplied by the local curvature $2\\pi^2 V_0\\sin\\psi/(24\\,\\mathrm{h})^2$; for a uniform exposure the variance becomes $\\Delta t^2/12$, producing the compact scaling $\\propto \\cos(\\mathrm{lat})\\cos(\\delta)\\sin\\psi\\,(\\Delta t)^2$ in Equation (10). This machinery turns a correction that is normally a single number into a filter over the exposure, and it predicts that flux curves concentrated at one end shrink the error while curves with a mid-exposure dip grow it.","core_discovery":"At the center of the paper is the difference between two averages of the barycentric-corrected velocity over an exposure: the value at the photon-weighted midpoint time, $v(\\langle t\\rangle)$, and the true photon-weighted average $\\langle v\\rangle$ from Equation (3). Since the diurnal component of $v_B(t)$ is sinusoidal, evaluating at $\\langle t\\rangle$ always falls inside the curve's curvature, so the second-order error has a definite sign and scales with $\\cos(\\mathrm{lat})\\cos(\\delta)\\sin\\psi\\,(\\Delta t)^2$; for a uniform exposure it is about $1.32\\,\\mathrm{m\\,s^{-1}}\\cos(\\mathrm{lat})\\cos(\\delta)\\sin\\psi\\,(\\Delta t/\\mathrm{1\\,h})^2$. At a low-latitude observatory the worst case allowed by a 30-degree altitude limit gives 0.25 m/s for a 30-minute exposure and 1.0 m/s for 60 minutes, with a V-shaped flux dip amplifying the effect by about 1.5. Full numerical simulations of the barycentric correction reproduce the analytic scaling, and real exposure-meter data for exposures up to 20 minutes shows the constant-flux approximation recovers the correction to within about 1 cm/s. The paper's conclusion is that instruments aiming at 10 cm/s precision must store the flux curve and apply photon weights to the correction itself.","pith_inferences":["Beyond the paper's own tests, the same curvature argument should apply to any time-varying correction over an exposure, including the chromatic dependence of the barycentric correction, so instruments that archive only midpoint times may have hidden systematics in other correction terms as well.","The paper does not quantify survey-level impact, but because the error depends on hour angle, nightly-offset fits will partially absorb it; a testable prediction is that archival RV residuals should show a sinusoidal pattern in hour angle with amplitude growing as $\\Delta t^2$.","One could exploit the opposite sign of the error east and west of the meridian: pairs of observations taken at $\\pm$ hour angles would show a symmetric curl if this is the dominant systematic, providing a clean null test."],"forward_implications":["Storing the exposure meter flux curve in the raw data makes the barycentric correction reproducible at the 1 cm/s level even for hour-long exposures, independent of weather and guiding changes.","Without photon weighting, a 30-minute low-latitude exposure near the celestial equator can be wrong by 0.25 m/s, and a 60-minute exposure by 1.0 m/s, in the worst observing geometry.","For archived data that only has the photon-weighted midpoint time, adding the constant-flux correction term of Equation (10) removes most of the second-order error, and real data out to 20 minutes supports this approximation to about 1 cm/s.","Observers can reduce the error by scheduling near the meridian, shortening exposures, or choosing high-declination targets; the error is largest for rising or setting targets observed through the east or west."],"supporting_citations":[{"why":"Defines the 1 cm/s barycentric correction formalism and the multiplicative correction that this paper re-weights.","marker":"Wright & Eastman (2014)"},{"why":"Provides the numerical barycentric-correction implementation used for the full simulations in Section 4.","marker":"Kanodia & Wright (2018)"},{"why":"First described the multi-channel exposure meter and the chromatic dependence that motivates photon-weighted corrections.","marker":"Blackman et al. (2017)"},{"why":"Set the 10 cm/s precision context and first pointed out the second-order problem in an early draft of this work.","marker":"Fischer et al. (2016)"},{"why":"Supplies the real exposure-meter observations used to validate the constant-flux approximation on data up to 20 minutes.","marker":"Blackman et al. (2019)"},{"why":"Describes the exposure-meter light-splitting needed to measure the per-wavelength flux curve $f(t)$.","marker":"Landoni et al. (2014)"}],"fun_headline_variants":["Use photon weights for barycentric correction, or face >1 m/s errors","Barycentric correction: midpoint time is not enough for precise RVs","Exposure-time curvature in barycentric correction can cost >1 m/s","Store flux curves: photon-weighted barycentric correction matters","Second-order barycentric errors can reach >1 m/s in real observations"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The paper's error estimates rest on modeling the barycentric velocity over an exposure as a pure diurnal sinusoid with a constant hour-angle rate; everything else in the correction is assumed to curve negligibly over tens of minutes, and the real-data check only reaches 20-minute exposures.","fun_headline_variants_meta":{"raw":{"variants":["Use photon weights for barycentric correction, or face >1 m/s errors","Barycentric correction: midpoint time is not enough for precise RVs","Exposure-time curvature in barycentric correction can cost >1 m/s","Store flux curves: photon-weighted barycentric correction matters","Second-order barycentric errors can reach >1 m/s in real observations"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000926,"raw_usage":{"total_tokens":4006,"prompt_tokens":1020,"completion_tokens":2986,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":636,"completion_tokens_details":{"reasoning_tokens":2893}},"tokens_in":636,"tokens_out":2986,"duration_ms":21591,"temperature":1.0,"reasoning_tokens":2893,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T15:26:17.988783+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a 60-minute exposure at a low-latitude site on a target with declination near zero observed through the east, record the exposure-meter flux curve, and compare $\\langle v\\rangle$ computed from the full flux curve with $v(\\langle t\\rangle)$ computed at the photon-weighted midpoint; the predicted difference is about 1 m/s (scaled by the actual flux shape). A null result at the 10 cm/s level would show the diurnal-curvature model is incomplete.","supporting_citations":[{"cited_title":"T., Ong J","cited_arxiv_id":null,"evidence_quote":"Supplies the real exposure-meter observations used to validate the constant-flux approximation on data up to 20 minutes."},{"cited_title":"M., Cabral A., Cristiani S., Megevand D., 2014, in Ground-based and Airborne Instrumentation for Astronomy V","cited_arxiv_id":null,"evidence_quote":"Describes the exposure-meter light-splitting needed to measure the per-wavelength flux curve $f(t)$."}],"review_version":1}