{"id":"e7c0c7b4-bcf9-4b48-b942-70e60f242b2f","arxiv_id":"2411.18763","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":6,"one_line_summary":"First full on-axis Mueller matrix calibration of the FAST 19-beam receiver, with time-variable parameters and usable average calibration for 2020 to 2022.","lead":"This paper calibrates the polarization response of the 19-beam receiver on China's Five-hundred-meter Aperture Spherical Telescope (FAST), measuring how the instrument distorts the polarization of incoming radio waves. It provides tables of calibration parameters and shows that the telescope's polarization distortion changes over months to years, so regular calibration is needed for reliable measurements.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Circular-polarization detection threshold is understated: I→V leakage error is 2σϵ, making 1.5% about 3σ, not the claimed 5σ.","rationale":"After checking the mathematics, the reader's α=0 concern is not the most load-bearing. Because M_F with χ=0 and M_SKY are both rotations, M_F M_SKY = R(θ+α); the α=0 fit simply rotates the fitted source Stokes parameters and leaves M_A M_IF (i.e., ΔG, ψ, ε, φ) unchanged. Thus the claimed time variability of ε and φ is not an artifact of a constant or epoch-varying rotation error. The real issue is the error propagation in Section 6.3: Table 2's σϵ is used as if it were the I-to-V leakage uncertainty, but the actual leakage is 2ε sin(φ+ψ), so the proper amplitude uncertainty is 2σϵ, plus φ/ψ covariance terms. This makes the 1.5% circular-polarization high-confidence threshold about 3σ rather than 5σ. The calibration itself remains valuable, but the abstract and conclusion overstate the reliability of weak circular-polarization detections. The verdict should remain CONDITIONAL pending a corrected threshold.","tokens_in":18912,"tokens_out":41833,"duration_ms":358560,"concrete_test":"Propagate the Table 2 uncertainties for the worst beam (e.g., M12) through the exact I-to-V Mueller entry 2ε sin(φ+ψ), including covariance terms from σφ and σψ, and compute the 5σ threshold on V/I. If that threshold exceeds 1.5%, the paper's headline circular-polarization detection limit is understated and must be revised.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The high-confidence detection thresholds in Section 6.3 are derived by equating the Mueller-matrix parameter uncertainties in Table 2 with the residual leakage uncertainties. This propagation is incorrect for the circular channel: in Eq. 10 the I-to-U and I-to-V leakages are 2ε cos(φ+ψ) and 2ε sin(φ+ψ) (α=0), so the leakage amplitude uncertainty is at least 2σϵ, not σϵ. For the worst beam, σϵ ≈ 0.30% (Table 2, M12), giving a 0.6% amplitude error before the substantial contributions from σφ (up to 30°) and σψ (up to 3°) are included. The abstract's claim that V/I ≳ 1.5% is a 5σ high-confidence detection therefore corresponds to roughly 2.5–3σ once the leakage is propagated properly. The linear-polarization threshold is conservative for the opposite reason (I-to-Q is ΔG/2, not ΔG), so the circular threshold is the load-bearing number that needs correction.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents the polarization calibration of the FAST L-band 19-beam receiver at 1420 MHz using spider and on-the-fly observations of polarized calibrators from 2018 to 2023. The authors fit a five-parameter Mueller matrix model (ΔG, ψ, ε, φ, and, for off-center beams, α) to derive on-axis polarization responses for the central and 18 off-center beams. They report that after correction, strong linearly polarized sources are measured to ~0.2% in fractional polarization and ~0.5° in angle, cross-checked against VLA measurements. They also provide weighted-mean 19-beam Mueller matrices for 2020–2022 and claim that, after applying these averages, fractional linear polarization ≳10% and fractional circular polarization ≳1.5% constitute high-confidence (5σ) detections. The paper additionally reports time variability of central-beam parameters and hints at an east-west asymmetry in the off-center beam response.","tokens_in":19186,"tokens_out":5708,"duration_ms":58978,"significance":"If the calibration is sound, this is a valuable long-term characterization of the FAST 19-beam receiver's polarimetric response, providing both a reference data set and a practical calibration product for the community. The paper includes error bars on all fitted parameters, external cross-checks against VLA polarization measurements, and makes the average Mueller matrices available for users. The main strengths are the multi-year baseline, the explicit treatment of systematic effects (noise diode, reflector surface), and the attempt to quantify detection thresholds after applying the calibration. The central calibration accuracy claim (0.2% and 0.5°) appears well supported by the VLA comparison. However, the derivation of the circular-polarization detection threshold contains a propagation error that weakens one of the headline claims.","major_comments":[{"comment":"The high-confidence detection threshold for circular polarization is understated. The authors use the tabulated σ_ε (up to 0.30% for M12 in Table 2) as the 1σ amplitude uncertainty of the I→U and I→V leakages. However, with α=0 the relevant matrix elements in Eq. (10) are 2ε cos(φ+ψ) and 2ε sin(φ+ψ), so the leakage amplitude uncertainty is at least 2σ_ε ≈ 0.6% before including the substantial contributions from σ_φ (up to ~30°–100° in Table 2) and σ_ψ (up to ~3°). For the worst beam the combined 1σ error in fractional circular polarization is ≳0.65%. The abstract's statement that V/I ≳ 1.5% is a 5σ high-confidence detection therefore corresponds to roughly 2.5–3σ. The linear-polarization threshold is conservative for the opposite reason (the I→Q leakage is ΔG/2, not ΔG), but the circular threshold is load-bearing and must be re-derived, and the abstract and Section 6.3 corrected accordingly.","section":"Section 4.2.2 and Section 5.1"},{"comment":"The degeneracy between α and the receiver rotation angle θ for the central beam is acknowledged, but its impact on the reported time variability is not assessed. The paper fixes α=0 for the central beam because θ may have a systematic error, yet Section 5.1 reports strong temporal trends in ε (from about −1% to +2%) and φ (from about 150° to −20°). If the rotation-angle error varies with time, these trends could be partly absorbed into the fitted ε and φ. Since the abstract claims that several Mueller-matrix parameters show time variability, the authors should quantify the sensitivity of ε and φ to plausible rotation-angle errors, for example by fitting α as a free parameter for a subset of epochs or by propagating a time-varying θ offset through the fit. Without such a test, the time-variability claim for ε and φ is not fully supported.","section":"Section 6.3"},{"comment":"The error-budget statement for linear polarization is also inaccurate: the text says that σ_ΔG of about 2% implies an I→Q leakage uncertainty of up to 2%, but the (1,2) element of the Mueller matrix in Eq. (10) is ΔG/2 (for α=0), so the leakage amplitude uncertainty is at most 1%, not 2%. This makes the quoted 10% threshold conservative (overestimated), but the description of the error budget should be corrected for consistency with the model equations.","section":"Section 6.3"}],"minor_comments":[{"comment":"The text spells 'Muller' where 'Mueller' is meant (e.g., 'solve the Muller matrices'); please correct throughout.","section":"§4.2.3"},{"comment":"The phrase 'high confident detections' should read 'high-confidence detections' for grammatical consistency.","section":"Section 6.3 and Conclusion"},{"comment":"The annotation 'PASRC (**UNCORRECTED FOR MASTRO**)' appears to contain a placeholder or undefined acronym; please clarify or remove it.","section":"Figure 5 caption"},{"comment":"The phrase 'full illumination of the Five-hundred-meter Aperture Spherical Telescope' is slightly awkward; consider 'full illumination of the Five-hundred-meter Aperture Spherical Telescope (FAST)' for readability.","section":"Abstract"}],"recommendation":"major_revision","confidential_remarks":"The paper is a solid instrumental calibration study with a long baseline and useful products. The main issue is the incorrect propagation of leakage uncertainties into the circular-polarization detection threshold, which affects a headline claim. The α=0 degeneracy also deserves a sensitivity check before the time-variability statement can be taken at face value. With these fixes the paper would be a strong contribution to the FAST calibration literature."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First, what's new: they give the first full on-axis Mueller matrix characterization of all 19 beams of the FAST L-band receiver, from spider observations for the central beam and OTF+spider for the off-center beams. The long 2018-2023 baseline shows real time variability in ΔG, ϵ, φ, and they see an east/west reflector asymmetry in off-center beams. They also provide an average 2020-2022 calibration table that observers can actually use. The cross-check against VLA polarimetry within ~0.2% in fractional polarization is solid evidence the method works.\n\nThe paper is honest and careful throughout. Errors are propagated in the fits, the degeneracy between α and θ is explicitly discussed (they set α=0 for the central beam), and they state the caveats for applying the average matrices clearly. The data it produces will be used by anyone doing 1420 MHz polarimetry with FAST.\n\nWhere it's soft:\n\n1. The circular-polarization detection threshold in §6.3 is understated. The I→U and I→V leakage amplitudes are 2ε cos(φ+ψ) and 2ε sin(φ+ψ) when α=0, so the amplitude uncertainty is at least 2σϵ, not σϵ. For the worst beam (M12, σϵ≈0.30%) that's 0.6% just from ϵ, before the sizable σφ (up to ~30°) and σψ contributions. So the abstract's claim that V/I ≳ 1.5% is a 5σ detection is more like ~3σ. The linear-polarization threshold (≳10%) is conservative for the opposite reason (I→Q leakage is ΔG/2), but the circular number needs a corrected propagation. This is a clear, fixable error.\n\n2. The α=0 degeneracy for the central beam means the reported time trends in ϵ and φ could partially absorb a slowly varying mechanical rotation error. The authors acknowledge the degeneracy but don't discuss this systematic. It doesn't invalidate the calibration—the correction still works—but the time-variability interpretation should be softened.\n\n3. Minor: the time-resolved per-epoch data are not released, only the average table. For re-analysis that would be useful, but it's not essential.\n\nOverall: this is a solid instrument-calibration paper with a load-bearing but addressable flaw in the detection-threshold propagation. It deserves a serious referee and should be published after the authors fix that number.","headline":"A genuinely useful, carefully documented FAST 19-beam polarization calibration whose main caveat is a factor-of-two understated circular-leakage threshold; worth refereeing with a requested fix.","tokens_in":19710,"tokens_out":3722,"would_cite":true,"duration_ms":118593,"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":"This paper establishes that a five-parameter Mueller matrix model, fitted to spider and on-the-fly observations, describes the FAST 19-beam receiver's on-axis polarization response well enough that strong sources are measured to about…","keywords":["polarization calibration","Mueller matrix","FAST 19-beam receiver","spider observations","on-axis leakage","Stokes parameters","radio polarimetry","fractional polarization"],"falsifier":"Measure the mechanical rotation angle with an independent encoder during a spider observation and compare the residuals of a fit that fixes $\\alpha=0$ with one that lets $\\alpha$ float; alternatively, illuminate the central feed with a laboratory source of known elliptical polarization to measure $\\alpha$ directly. If a nonzero $\\alpha$ or a rotation-encoder offset emerges, the central-beam parameters shown in the paper and the 2020–2022 average matrices would need re-derivation.","tokens_in":18738,"feed_emoji":"📡","tokens_out":5631,"duration_ms":45682,"temperature":0.7,"pith_summary":"This paper establishes that the on-axis polarization response of the Five-hundred-meter Aperture Spherical Telescope's 19-beam L-band receiver can be captured by a five-parameter Mueller matrix model, solved with spider drift scans for the central beam and on-the-fly maps for the other 18 beams. After the fitted matrices are applied, strong linearly polarized calibrators are recovered to about 0.2% in fractional polarization and 0.5 degrees in polarization angle. The data also show that several Mueller-matrix parameters drift on month-to-year timescales, so polarization calibration must be repeated fairly often. For users of the published 2020–2022 average matrices, the paper gives practical detection thresholds: on-axis fractional linear polarization above roughly 10% and on-axis fractional circular polarization above roughly 1.5% can be treated as high-confidence detections.","feed_headline":"Five-parameter fix calibrates FAST's 19 beams to 0.2 percent","feed_subtitle":"Average 2020–2022 matrices set detection thresholds of 10% linear and 1.5% circular polarization.","key_machinery":"The load-bearing object is the Mueller matrix formalism: a $4\\times4$ transfer matrix mapping intrinsic Stokes parameters $(I,Q,U,V)$ to observed ones, written as a product of a sky-rotation matrix and a receiving-system matrix. The receiving-system matrix factors into feed, feed-imperfection, and amplifier-chain matrices with parameters $\\Delta G$, $\\psi$, $\\alpha$, $\\epsilon$, and $\\phi$; the model keeps all orders in $\\alpha$, $\\phi$, and $\\psi$ and first order in $\\epsilon$ and $\\Delta G$. The observations that carry the argument are spider scans, in which a linearly polarized calibrator is drifted through the beam at receiver rotation angles $-60^\\circ$ to $+60^\\circ$, giving enough rotation-angle coverage to fit the parameters and the source Stokes parameters simultaneously. On-the-fly maps of several calibrators extend the same model to the 18 off-center beams using the central-beam source polarizations as known inputs.","core_discovery":"The central claim is that the on-axis polarization leakage of the FAST 19-beam receiver at 1420 MHz is fully described by the product structure $M_{\\mathrm{TOT}} M_{\\mathrm{SKY}}$, with $M_{\\mathrm{SKY}}$ a rotation by the receiver angle $\\theta$ and $M_{\\mathrm{TOT}}$ built from five parameters: the relative gain error $\\Delta G$, the residual electronic phase $\\psi$, the feed ellipticity angle $\\alpha$, and the cross-coupling amplitude $\\epsilon$ and phase $\\phi$. For the central beam, $\\alpha$ is fixed to zero because it is degenerate with a mechanical rotation error, leaving a four-parameter fit from spider observations; for the 18 off-center beams, $\\alpha$ is kept free and the source Stokes parameters are taken from the central-beam solution. The paper reports that after this calibration the recovered polarization percentage and angle of strong sources agree with independent measurements to about 0.2% and 0.5 degrees, that the calibrated fractional Stokes parameters are consistent across all 19 beams, and that the 2020–2022 average matrices make $\\gtrsim 10\\%$ on-axis fractional linear polarization and $\\gtrsim 1.5\\%$ on-axis fractional circular polarization high-confidence detections.","pith_inferences":["Beyond the paper, if the same average matrices are applied to archival FAST data from 2020–2022, sources with fractional linear polarization below the 10% threshold could have their polarization angles systematically affected by the uncalibrated $\\phi$ uncertainty, so the matrix should be used as a detection screen rather than a precision measurement for weak sources.","Beyond the paper, the time variability of $\\epsilon$ and $\\phi$ could be tested against independent noise-diode equivalent-temperature measurements across the same epochs; the paper already ties $\\Delta G$ to that fluctuation, so a similar check would separate receiver drift from calibration artifact.","Beyond the paper, the eastern/western reflector-surface difference suggests that a pointing- or surface-dependent correction may be needed for surveys that combine beams; a dedicated experiment rotating the receiver at a fixed source position across different surface sectors could quantify this.","Beyond the paper, because the same receiver feeds serve the pulsar backend, pulsar polarization observations using these matrices should be validated on a polarized pulsar of known rotation measure, which would also test applicability beyond continuum calibrators."],"forward_implications":["After calibration, strong linearly polarized sources can be measured to roughly 0.2% in fractional polarization and 0.5 degrees in polarization angle with the central beam.","The published 2020–2022 average Mueller matrices are usable for FAST spectral-line and pulsar observations at small zenith angle, with the caveat that the noise-diode phase must be calibrated first.","Users can treat on-axis fractional linear polarization $\\gtrsim 10\\%$ and fractional circular polarization $\\gtrsim 1.5\\%$ as high-confidence detections when using the average matrices; weaker signals need contemporaneous spider calibrations.","Because several Mueller-matrix parameters drift on month-to-year timescales, polarization calibration should be repeated frequently rather than assumed stable.","The central-beam parameters do not depend strongly on the reflector surface at small zenith angles, but the off-center beams show eastern/western differences, indicating surface-dependent leakage."],"supporting_citations":[{"why":"Supplies the Mueller-matrix treatment of on-axis leakage and the spider-observation technique used throughout the paper.","marker":"Heiles et al. 2001a"},{"why":"Formalizes single-dish polarization calibration with Mueller matrices, the framework the paper adopts and simplifies.","marker":"Heiles 2002"},{"why":"Provides the detailed Stokes-V convention and the matrix decomposition for the receiving system.","marker":"Robishaw & Heiles 2021"},{"why":"Gives the polarization standards (3C48, 3C138, 3C286) whose known fractional polarization and angle anchor the fits.","marker":"Perley & Butler 2013"},{"why":"Commissioning report supplying receiver parameters, pointing accuracy, and noise-diode stability context.","marker":"Jiang et al. 2020"},{"why":"The RHSTK package used for gain, phase, bandpass, and polarization calibration of the data.","marker":"Heiles & Robishaw 2022"},{"why":"Global ionospheric maps used to correct Faraday rotation in the measured polarization angles.","marker":"Meng et al. 2024"},{"why":"Earlier FAST HI Zeeman result that used the preliminary central-beam calibration, motivating the present work.","marker":"Ching et al. 2022"}],"fun_headline_variants":["FAST 19-beam polarimeter calibrated to 0.2% accuracy","Average Mueller matrices set high-confidence polarization limits","Five-parameter model fixes FAST beam polarization leakage","Time-variable Mueller matrices demand frequent recalibration","Off-center beams show east-west asymmetry in calibration"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument depends on setting the central beam's feed ellipticity $\\alpha$ to zero because $\\alpha$ and the receiver rotation angle $\\theta$ enter only as their sum; if the real feed has nonzero ellipticity, or if the mechanical rotation has a small time-dependent error, the fitted values of $\\Delta G$, $\\psi$, $\\epsilon$, and $\\phi$ could carry biases, and the reported time variability of $\\epsilon$ and $\\phi$ could be partly an artifact.","fun_headline_variants_meta":{"raw":{"variants":["FAST 19-beam polarimeter calibrated to 0.2% accuracy","Average Mueller matrices set high-confidence polarization limits","Five-parameter model fixes FAST beam polarization leakage","Time-variable Mueller matrices demand frequent recalibration","Off-center beams show east-west asymmetry in calibration"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000544,"raw_usage":{"total_tokens":2680,"prompt_tokens":1100,"completion_tokens":1580,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":716,"completion_tokens_details":{"reasoning_tokens":1504}},"tokens_in":716,"tokens_out":1580,"duration_ms":11329,"temperature":1.0,"reasoning_tokens":1504,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T10:54:20.294086+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the mechanical rotation angle with an independent encoder during a spider observation and compare the residuals of a fit that fixes $\\alpha=0$ with one that lets $\\alpha$ float; alternatively, illuminate the central feed with a laboratory source of known elliptical polarization to measure $\\alpha$ directly. If a nonzero $\\alpha$ or a rotation-encoder offset emerges, the central-beam parameters shown in the paper and the 2020–2022 average matrices would need re-derivation.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The RHSTK package used for gain, phase, bandpass, and polarization calibration of the data."}],"review_version":1}