{"id":"37f3b223-a044-4e7f-8be9-29d03ebc54c8","arxiv_id":"2607.19585","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"International LOFAR stations are 20–45% more sensitive to sources on the rising side of the sky than the setting side, an asymmetry observed in all 11 tracked pulsars and across three stations.","lead":"Radio astronomers mapped how the Irish LOFAR telescope's sensitivity varies across the sky by tracking 11 bright pulsars over 170 hours. They found the station sees sources 20–45% better when they are rising rather than setting — an asymmetry the standard beam model does not predict.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Unmodelled LST-dependent sky noise (Fig. 10) may explain the rise/set S/N asymmetry, so the claimed excess over the Hamaker beam model is not yet established.","rationale":"The reader's verdict CONDITIONAL already flags the unmodelled Tsys(LST) issue as a secondary premise. My stress-test elevates this to the single most load-bearing concern because the paper's own Fig. 10 provides direct evidence that Tsys varies with LST at a fixed elevation, and the claimed beam-model discrepancy is quantified only after asserting Tsys is 'more or less constant' (Sect. 5). The observed consistency across pulsars and stations argues against source variability as the primary cause, but it does not separate Aeff from Tsys. A joint fit of beam model plus Tsys(LST) is needed to establish the magnitude of any beam-model error. The paper otherwise provides strong empirical support for a real rise/set asymmetry, so the verdict should remain CONDITIONAL pending this test rather than moving to ACCEPT or REJECT.","tokens_in":25213,"tokens_out":4781,"duration_ms":47502,"concrete_test":"Using the SE607 Fig. 10 data (or PyGDSM sky models for IE613), construct Tsys(LST) at the relevant elevation. For each pulsar in Table 3, predict S/N(t) = G_DreamBeam(el,az)/Tsys(LST(t)), with gain from DreamBeam and Tsys from the measured noise curve; then compute the rise/set S/N ratio at matched elevations. If the predicted ratio reproduces the observed 20–45% asymmetry, the case for an unmodelled beam-model deficit collapses. If it leaves a residual >10%, the X/Y-weighting explanation is supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central attribution — that the observed rise/set S/N asymmetry exceeds the Hamaker/DreamBeam beam model and may reflect imbalanced X/Y weighting — rests on treating Tsys as 'more or less constant' over each track (Sect. 5). But Fig. 10 (SE607 fixed-elevation, two azimuths) shows LST-dependent sky noise variation of a scale comparable to the reported 20–45% asymmetry. Since a pulsar's rising and setting phases occur at different LSTs (separated by ~12h for circumpolar sources), a time-varying Tsys(LST) can masquerade as an azimuthal beam asymmetry. The paper never jointly models Tsys(LST) with the beam response; it only notes the noise variation qualitatively. Without this joint model, the quantitative claim that the beam model under-predicts the asymmetry is not closed. The qualitative asymmetry is likely robust, but its attribution to Aeff rather than Tsys is not.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a characterization of the high-band antenna (HBA) response of the Irish LOFAR station (IE613) using long-track observations of 11 bright pulsars, totaling ~170 hours. The authors measure the S/N of folded pulse profiles as a function of elevation and azimuth, and compare the resulting sensitivity pattern with the DreamBeam implementation of the Hamaker/Hamaker beam model. Their main empirical finding is a pronounced azimuthal asymmetry: for every pulsar observed, the S/N is higher on the rising side of the track than on the setting side, with the difference often reaching ~20% and up to 45%. The same trend is reported for the Swedish station (SE607) using PSR B0329+54 and is claimed to be consistent with earlier Polish station (PL612) results. The paper argues that the standard beam model predicts a near-symmetric response and proposes, as a candidate explanation, an imbalanced weighting of the X and Y polarization channels. The work concludes with operational recommendations for LOFAR 2.0 and SKA-Low.","tokens_in":25328,"tokens_out":6608,"duration_ms":58747,"significance":"If the quantitative claim holds, this is an important result for low-frequency aperture-array instruments: international LOFAR stations would carry a position-dependent sensitivity asymmetry at the tens-of-percent level that is not captured by the standard Hamaker beam model, with direct implications for flux calibration, polarimetry, scheduling, and the design of LOFAR 2.0 and SKA-Low. The paper has substantial strengths: a large, homogeneous dataset (170 h over 11 pulsars), an independent check at a second station (SE607), a qualitative comparison with previously published PL612 data, and a transparent acknowledgment of limitations (e.g., the scintillating pulsar B0809+74 and the unclear physical origin of the effect). The dataset in Appendix B is a useful community resource. However, the central quantitative claim is not yet error-budgeted, and the attribution of the asymmetry to Aeff rather than to time-variable sky noise is not closed. These issues are fixable and do not negate the qualitative finding, but they are load-bearing for the paper's main conclusion.","major_comments":[{"comment":"The attribution of the rise/set S/N asymmetry to the beam response (through an imbalanced X/Y weighting) rests on the statement that 'except near the horizon ... the Tsys factor is more or less constant for the majority of the trajectory of each source.' However, the SE607 fixed-elevation test in Fig. 10 shows that sky noise varies with LST at a scale that is comparable to the reported 20–45% asymmetry. For a circumpolar source, the rising and setting phases occur at LSTs separated by roughly 12 hours, so a time-varying Tsys(LST) can masquerade as an azimuth-dependent Aeff. The paper notes this variation qualitatively but never jointly models Tsys(LST) with the DreamBeam Aeff model. To establish that the standard beam model under-predicts the observed asymmetry, the authors must either subtract a data-derived Tsys(LST) model from the S/N or quantitatively demonstrate that the observed LS","section":"Section 5 (Discussion), Fig. 10"},{"comment":"The main quantitative claims — 'responses often being ~20% better on the rising side' (Sect. 4) and 'a discrepancy that reaches 45%' (Sect. 6) — are not supported by error estimates. Table 3 lists S/N values to two decimal places but gives no uncertainties, no off-pulse noise estimates, and no measure of profile-fitting or scrunching errors. Moreover, the statement that the effect is seen 'for every pulsar' is presented without a statistical test. The authors should provide per-point S/N errors, propagate them into rise/set ratios, and report a significance for the sign of the asymmetry for each pulsar and for the sample as a whole (for example, a sign test). Without this, the magnitude and universality of the effect are not quantitatively established.","section":"Section 4 (Results), Table 3 / Appendix B"},{"comment":"The paper explicitly flags PSR B0809+74 as showing interstellar scintillation and advises 'extra caution ... when interpreting the measured values' (Sect. 4 and the note to Fig. 11). Nevertheless, the same section states that 'every pulsar we observe is seen more significantly as it rises,' and this pulsar is included in the aggregate claim. The Table 3 S/N values for B0809+74 show large, non-monotonic scatter (e.g., 872, 1487, 1419, 1139 at adjacent elevations of 38.6, 40.3, 42.7, 45.7 degrees), which is consistent with scintillation rather than a smooth antenna response. This is an internal inconsistency: a source the authors themselves flag should not be included in an 'all pulsars' claim without a demonstration that the asymmetry survives its exclusion. Please re-run the aggregate statistics with and without B0809+74.","section":"Section 4 (Results), Fig. 11 and Table 3 (B0809+74)"},{"comment":"The assertion that the observed asymmetry 'is more than one can explain/model from' the Hamaker beam model (Sect. 6) is not quantitatively demonstrated. Fig. 8 provides a single-frequency (150 MHz) DreamBeam calculation for one pulsar (PSR B1508+55) and states qualitatively that the equal-weighting response is 'highly symmetric.' Fig. 5 compares the data to an ad hoc cos^2(ZA) fit, not to the DreamBeam model. Since the observed asymmetry is claimed to be present across the frequency band (Fig. 6), the model comparison should be made per frequency sub-band, and the residual (data minus model Aeff/Tsys) should be presented as a quantitative function of elevation, azimuth, and LST. As written, the 'exceeds the model' claim is a visual impression rather than a measured discrepancy.","section":"Section 5 (Discussion), Fig. 8 and Fig. 5"}],"minor_comments":[{"comment":"The expression for gain, G = 4π^2 Aeff/λ^2, should be G = 4π Aeff/λ^2 (standard effective-aperture relation). The extra factor of π does not affect the rest of the analysis, but the formula should be corrected.","section":"Eq. (1)"},{"comment":"The caption says 'The dotted orange line shows the combined Stokes I response. A representative unequal weighting is shown in dotted orange.' The color/line-style designation appears duplicated; clarify which line corresponds to equal weighting and which to the illustrative unequal (a = 0.5) weighting.","section":"Fig. 8 caption"},{"comment":"'Manch-ester et al.' is a typo for 'Manchester et al.' (the ATNF pulsar catalogue reference).","section":"Section 3"},{"comment":"'calbrating' should be 'calibrating' in the opening paragraph.","section":"Section 1"},{"comment":"The phrase 'responses often being ~20% better on the rising side' would benefit from a precise definition of 'often' — for example, the fraction of pulsars or the fraction of elevation-matched track points for which the rising S/N exceeds the setting S/N.","section":"Section 4, summary paragraph"},{"comment":"The column header says 'relative S/N' while the text (Sect. 3) says the reported values are the S/N from pdmp. Clarify whether the values are normalized in any way or are absolute S/N values.","section":"Table 3"}],"recommendation":"major_revision","confidential_remarks":"The manuscript fits the journal's scope well and the dataset is valuable. The main risk is the Tsys(LST) confound: the SE607 fixed-elevation test (Fig. 10) shows noise variations of a scale comparable to the claimed asymmetry, so the Aeff-vs-Tsys attribution needs to be made quantitatively. I would encourage the editor to send the revised manuscript to a referee familiar with LOFAR station calibration and the RIME formalism. The qualitative asymmetry is credible, but the current version does not close the quantitative case."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know: this is an empirical paper with a genuinely new observational result—an azimuthal (rise/set) asymmetry in LOFAR HBA station sensitivity, seen in 11 pulsars at IE613, reproduced at SE607, and back-read into PL612 data. The data are extensive: 170 hours, full tracks, careful RFI cleaning. That much is solid. The weak point is that the quantitative size of the effect (20–45%) is quoted without propagated uncertainties, and the argument that it is a beam (Aeff) effect rather than a time-varying sky-noise (Tsys) effect is not closed.\n\nWhat is new: prior elevation-only fits (Noutsos, Błaszkiewicz) did not report azimuthal asymmetry. The paper extends to full azimuthal characterization, uses DreamBeam/Hamaker comparison, and proposes a plausible X/Y weighting imbalance. It is also transparent about limits: it flags B0809+74 scintillation, says the origin is unclear, and the three-station concordance is good evidence.\n\nThe main soft spot, as the stress-test notes, is Tsys. The paper asserts Tsys is “more or less constant” over most of the track (Sect. 5), but Fig. 10 shows LST-dependent sky noise at fixed elevation at SE607 at a scale that could be comparable to the claimed asymmetry. Since rising and setting phases occur ~12 hours apart in LST, a time-varying Tsys can masquerade as an azimuthal beam asymmetry. The paper never jointly models Tsys(LST) with the beam response; it only mentions it qualitatively. So the claim that the beam model under-predicts the asymmetry is not quantitatively established. Also Table 3 has no error bars, and the 20% vs 45% statements are not reconciled. The illustrative X/Y ratio a=0.5 in Fig. 8 is not fitted.\n\nThat said, the qualitative asymmetry is likely robust: it appears for all pulsars (with the B0809 caveat), across frequency sub-bands, and at multiple stations. The stress-test concern does not kill the paper; it just pushes the quantitative claim back to “needs a joint model.” This is a paper for pulsar astronomers, LOFAR calibration people, and the SKA-Low beam-modelling community. I’d bring it to reading group and would cite it—the empirical asymmetry is likely real and the paper documents the effect.\n\nRecommendation: send it to peer review. It deserves a serious referee. The referee should ask for a joint Tsys(LST)+beam fit and uncertainties on the asymmetry parameters. Don’t desk-reject.","headline":"A credible, well-resourced measurement of a rise/set sensitivity asymmetry at LOFAR stations, but the quantitative amplitude and the attribution to the beam rather than sky noise need more work before it becomes a calibration-grade result.","tokens_in":26025,"tokens_out":2772,"would_cite":true,"duration_ms":25455,"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":"A LOFAR station detects pulsars with roughly 20% higher signal-to-noise as they rise than as they set, a bias standard beam models do not predict.","keywords":["LOFAR","international station","beam model","azimuthal asymmetry","pulsars","signal-to-noise ratio","aperture array calibration","Hamaker formalism"],"falsifier":"Track a bright continuum calibrator with known stable flux over a full rise/set pass at the same station on two consecutive days, repeating the same LST window, and compare S/N at identical elevations on the rising and setting sides: if the S/N difference flips sign or disappears when the source's azimuth is mirrored, the asymmetry is caused by sky noise, not by a fixed polarisation imbalance. Alternatively, re-beamform recorded data with digitally equalised X/Y weights; if the hysteresis persists, the weighting is not the cause.","tokens_in":24998,"feed_emoji":"📡","tokens_out":4241,"duration_ms":39653,"temperature":0.7,"pith_summary":"This paper sets out to characterise the full-sky response of the high-band antennas of an international LOFAR station, using 170 hours of track observations of 11 bright, supposedly flux-stable pulsars. The central claim is that the station's sensitivity is not symmetric about the zenith: every measured pulsar shows a higher signal-to-noise ratio as it rises than when it sets, with the difference usually around 20% and reaching 45%. The standard Hamaker beam model, as implemented in the DreamBeam software, predicts an essentially symmetric response, so the observed asymmetry is not explained by current beam modelling. If the claim is right, calibration, flux-density measurement, and scheduling for LOFAR 2.0 and SKA-Low must account for a station-level, direction-dependent sensitivity bias of tens of percent. The authors suggest an imbalanced weighting of the two polarisation channels as a plausible cause.","feed_headline":"Rising pulsars show 20% higher signal at LOFAR","feed_subtitle":"International stations' sensitivity is asymmetric in azimuth, a bias standard beam models do not predict.","key_machinery":"The carrier of the argument is the station beam model built on the Radio Interferometer Measurement Equation and the Hamaker formalism, implemented in the DreamBeam software as Jones matrices; the relevant quantity is the Mueller element m_II, which maps incident Stokes I to measured power. The paper models the measured S/N as (a·X + b·Y)/Tsys, where X and Y are the two polarisation channel responses and a and b are their relative weights; an equal-weight model (a=b) gives the symmetric response, while a≠b produces the observed rise/set asymmetry.","core_discovery":"The discovery is an azimuthal asymmetry in the measured beam response of the station: the effective sensitivity, encoded in the signal-to-noise ratio (effectively Aeff/Tsys), is higher on the rising side of a source's trajectory than on the setting side at the same elevation. The effect is broadband, consistent across all 11 pulsars in the sample, and reproduced at the Swedish and Polish international stations, indicating a systematic, not local, property. The paper argues this cannot come from the Hamaker beam model alone, and that an unequal weighting of the X and Y polarisation contributions in the station beamforming can reproduce the hysteresis-like pattern in elevation.","pith_inferences":["If the asymmetry is a fixed X/Y weighting in station firmware, a testable fix is to equalise the two channels digitally and re-measure the S/N curve: the hysteresis should flatten, providing a cheap correction and confirming the explanation.","The paper does not jointly model the LST-dependent sky noise measured in the fixed-elevation SE607 test with the polarisation weighting; a joint fit of (aX+bY)/Tsys(LST) could separate an instrumental cause from a sky-noise cause.","Tracking a strong continuum source (e.g. Cassiopeia A) through the same rise/set geometry would test whether the effect is intrinsic to the beam or specific to pulsar processing such as folding and dedispersion.","Because the dipole orientation varies between international stations, comparing stations with very different array rotations could isolate the geometric mechanism behind the imbalanced weighting."],"forward_implications":["Observations scheduled on the rising side gain tens of percent in sensitivity, translating into large time savings since S/N scales only with the square root of observing time.","Flux densities and spectral indices measured with international LOFAR stations will be biased by this asymmetric response unless calibration explicitly includes per-polarisation weighting.","The same pulsar-track test can serve as an early validation diagnostic for LOFAR 2.0 and SKA-Low, before those systems are fully complete.","Unmodelled asymmetries of this size can contaminate pulsar timing, polarimetry, and transient searches, not just imaging.","The effect, if confirmed at other stations, implies that aperture-array calibration software must treat the beam as a function of azimuth as well as elevation and frequency."],"fun_headline_variants":["LOFAR sensitivity biased by azimuth: rise beats set by 20%","Azimuthal asymmetry in LOFAR beam challenges standard models","Pulsar tracks reveal LOFAR sees rising sources better","Unequal polarisation weighting explains LOFAR's directional bias","LOFAR beam model fails to predict 20% rise-set asymmetry"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The claim rests on the pulsars being flux-stable across each track and on Tsys being essentially constant with azimuth, so that the measured S/N differences trace the instrument rather than the source or the sky; the paper itself flags that PSR B0809+74 shows interstellar scintillation yet still includes it in the all-pulsars result, and the SE607 fixed-elevation test shows sky-noise variations of comparable size.","fun_headline_variants_meta":{"raw":{"variants":["LOFAR sensitivity biased by azimuth: rise beats set by 20%","Azimuthal asymmetry in LOFAR beam challenges standard models","Pulsar tracks reveal LOFAR sees rising sources better","Unequal polarisation weighting explains LOFAR's directional bias","LOFAR beam model fails to predict 20% rise-set asymmetry"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000242,"raw_usage":{"total_tokens":1403,"prompt_tokens":825,"completion_tokens":578,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":569,"completion_tokens_details":{"reasoning_tokens":487}},"tokens_in":569,"tokens_out":578,"duration_ms":5887,"temperature":1.0,"reasoning_tokens":487,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T12:19:12.216165+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Track a bright continuum calibrator with known stable flux over a full rise/set pass at the same station on two consecutive days, repeating the same LST window, and compare S/N at identical elevations on the rising and setting sides: if the S/N difference flips sign or disappears when the source's azimuth is mirrored, the asymmetry is caused by sky noise, not by a fixed polarisation imbalance. Alternatively, re-beamform recorded data with digitally equalised X/Y weights; if the hysteresis persists, the weighting is not the cause.","supporting_citations":[],"review_version":1}