{"id":"72f491b1-62b4-4098-8437-6053b4d5048d","arxiv_id":"2506.12272","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Position-angle shifts in polarized radio emission from three black widow pulsars reveal milligauss-strength line-of-sight magnetic fields in their eclipse media, with evidence of field reversals in two systems.","lead":"Three black widow pulsars were observed with the FAST radio telescope, and the polarization angle of their radio light rotated near each eclipse, implying magnetic fields of a few to tens of milligauss in the eclipsing material. The observations also suggest the magnetic field direction reverses in two of the three systems, a new clue to how these pulsars strip and interact with their companion stars.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The magnetic-field inference depends entirely on the untested assumption that the observed position-angle shifts are pure Faraday rotation; no frequency-dependence or profile-constancy check is presented.","rationale":"The reader's weakest assumption exactly matches the key unvalidated step: subtracting the out-of-eclipse PA swing from each sub-integration and attributing the residual entirely to Faraday rotation. The paper does not demonstrate the λ² scaling that would distinguish propagation effects from intrinsic PA variability, nor does it show that the PA residual is a constant offset across the pulse profile. These omissions leave the central claim conditional on a plausible but unverified assumption. The concern does not by itself warrant rejection, because the method is standard in the field and the results may well be correct; it does, however, justify the CONDITIONAL verdict already given. The proposed sub-band test is decisive and can be performed with the existing data, since the 400-MHz bandwidth intrinsically permits a λ² consistency check. The reader's other points (the sigma_RM values being derived from an arbitrary 5% depolarization assumption and the numerical inconsistency in Section 4.2) are secondary but reinforce the need for caution. Thus no verdict adjustment is recommended beyond what the reader already proposed.","tokens_in":11725,"tokens_out":7133,"duration_ms":86930,"concrete_test":"Split the calibrated FAST observations into four 100-MHz sub-bands (1.05–1.15, 1.15–1.25, 1.25–1.35, 1.35–1.45 GHz). For each sub-integration flagged as having a significant PA shift in Figure 3, compute the PA shift relative to the same out-of-eclipse reference in each sub-band. Fit ΔPA(ν) = ΔRM (c/ν)² + c₀. If the Faraday interpretation is correct, a single ΔRM must fit all sub-bands with c₀ ≈ 0 (mod π). Report the best-fit ΔRM and reduced chi-square; if the fit is poor or ΔRM is inconsistent across sub-bands, the PA shifts are not purely propagation-induced and the inferred magnetic fields are not supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 3.2.1 introduces the method: the PA swing of the average out-of-eclipse profile is used as the reference, and the PA shift (ΔPA) of each sub-integration is attributed to Faraday rotation in the eclipse medium via Equation (1), ΔPA = ΔRM λ². This is the sole channel by which the magnetic-field strengths and reversals are inferred. The paper does not test two falsifiable predictions of this model. First, a Faraday shift must scale as λ² across the 400-MHz band, yet the analysis uses only the full-band PA and a single effective wavelength (24 cm). Second, a pure rotation leaves the shape of the PA swing unchanged, so the per-pulse-phase difference between any sub-integration and the reference should be a constant offset modulo π. If the difference varies across pulse phase, intrinsic changes in emission geometry or orthogonal-mode ratio are present, and the single-number ΔPA is not a valid measure of Faraday rotation. These possibilities are not addressed for any of the three events; for J2055+3829 the detection rests on a single sub-integration, as explicitly acknowledged in Section 3.2.2. If the PA shifts are not purely Faraday, the claimed lower-limit magnetic fields and the evidence for reversals in B1957+20 and J1544+4937 would be unsupported. The central claim therefore stands or falls on this untested assumption.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents FAST 1250 MHz polarimetric observations of three black widow pulsars, PSRs B1957+20, J2055+3829, and J1544+4937, during and around their radio eclipses. For each pulsar, the authors measure flux density, dispersion measure, rotation measure, and linear/circular polarization as functions of orbital phase. They report position-angle (PA) shifts at the eclipse boundaries in all three systems, convert these shifts to rotation-measure changes via ΔPA = ΔRM λ², and then estimate lower limits on the line-of-sight magnetic field strength in the eclipse medium using ΔB_∥ = 1.23 μG ΔRM / ΔDM. They claim lower-limit field strengths of a few mG to tens of mG, and interpret sign changes of ΔRM in B1957+20 and J1544+4937 as evidence for magnetic field reversals. The paper also discusses depolarization at eclipse boundaries and compares the inferred σ_RM values with those seen in repeating fast radio bursts.","tokens_in":11934,"tokens_out":4877,"duration_ms":59655,"significance":"If the central claim is correct, this would provide the first direct measurement of the line-of-sight magnetic field strength in the eclipse medium of PSR B1957+20, and would add two more systems with evidence for ordered, sign-changing milligauss-level magnetic fields in black widow eclipses. Such measurements are important for understanding eclipse mechanisms and the magnetoionic environment of spider pulsars. The paper has several strengths: the PA shifts for B1957+20 and J1544+4937 are large compared with the quoted uncertainties and are directly visible in the figures; the observations are polarization-calibrated with standard tools; and the authors explicitly acknowledge the limitation of the J2055+3829 detection resting on a single sub-integration. However, the central inference depends on an untested assumption that the observed PA shifts are entirely due to Faraday rotation, and the lack of a frequency-dependence or profile-shape check is a significant gap.","major_comments":[{"comment":"The central inference that the observed ΔPA is purely Faraday rotation is not tested. A Faraday rotation measure would produce a PA shift that scales as λ² across the 400 MHz band, and a pure rotation would leave the shape of the PA swing unchanged (i.e., the difference between a sub-integration and the reference would be a constant offset modulo π at every pulse phase). The analysis uses only the full-band PA and a single effective wavelength of 24 cm, and does not check either prediction. If the intrinsic PA swing changes at the eclipse boundary (for example, because different emission regions or modes become visible), the inferred ΔRM and ΔB_∥ would be spurious. I request a frequency-resolved analysis (e.g., sub-band ΔPA values) and a per-pulse-phase residual analysis to support the Faraday interpretation before the magnetic field and reversal claims can be accepted.","section":"Section 3.2.1, Eq. (1)"},{"comment":"The detection of a PA shift for PSR J2055+3829 rests on a single sub-integration at orbital phase 0.402, with ΔPA = -23±4 deg and a 6% fractional linear polarization outside eclipse. The authors acknowledge this, but a single measurement is not sufficient to claim a PA shift at the eclipse boundary for this pulsar, especially given the low linear polarization fraction and the possibility of a noise fluctuation. The abstract and conclusions state that PA shifts are observed in all three pulsars; the J2055+3829 result should either be supported by additional independent sub-integrations or explicitly downgraded to a tentative detection.","section":"Section 3.2.2"}],"minor_comments":[{"comment":"The σ_RM values quoted in the discussion (13, 16, and 20 rad/m² for PSRs B1957+20, J2055+3829, and J1544+4937) are inconsistent with the values derived in Section 3.2 (15, 5, and 13 rad/m², respectively). Please reconcile these numbers and clarify the order in which the pulsars are listed.","section":"Section 4.2"},{"comment":"The statement that linear polarization depolarizes at orbital phases 0.120-0.158 and 0.195-0.348 while it remains unchanged at 0.158-0.220 is confusing because 0.158-0.220 overlaps with the second depolarized interval. Please clarify whether 0.158-0.220 is a region where the linear polarization is measured to remain unchanged, and if so, how this is consistent with the claimed depolarization in the overlapping range.","section":"Section 3.2.3"},{"comment":"The caption refers to 'vertical filled red areas in the fourth panel', but the panels are not numbered in the figure. Please add explicit panel labels (e.g., a, b, c) or describe the panels by their content to avoid ambiguity.","section":"Figure 3 caption"},{"comment":"The sentence 'For PSRs J2055+3829, and J1544+4937, our observations cover the entire binary orbital phases' contains a grammatical error and should read 'For PSRs J2055+3829 and J1544+4937, our observations cover the entire binary orbit'.","section":"Section 2"}],"recommendation":"major_revision","confidential_remarks":"The manuscript reports potentially important observations, but the central magnetic-field inference hinges on the Faraday-only assumption for the PA shifts. The data, which span 400 MHz of bandwidth, should permit a sub-band analysis to test the λ² dependence; the authors should be asked to perform this check. The J2055+3829 detection is too weak to carry the weight of the 'all three pulsars' claim in its current form. I recommend major revision rather than rejection because the requested tests appear feasible with the existing dataset."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This paper is worth a careful read. It reports the first direct line-of-sight magnetic field estimate for the eclipse medium of PSR B1957+20, the prototype black widow, and presents evidence for magnetic field reversals in B1957+20 and J1544+4937. The FAST data are good: new DM, RM, and eclipse-duration measurements for three spiders, with PA shifts shown directly in the sub-integration plots. The method is borrowed from Crowter et al. (2020), but the application to these pulsars is new, and the B1957+20 measurement is a legitimate first.\n\nThe soft spot is the central interpretation. The PA shifts are assumed to be pure Faraday rotation (ΔPA = ΔRM λ²), but the analysis uses only a single effective wavelength, so there is no λ² check. And the reference PA swing subtraction assumes the intrinsic PA swing does not change across orbital phase; if the emission geometry shifts at those phases, the inferred ΔRM values would be spurious. The stress-test note is right that neither falsifiable prediction is tested. For J2055+3829, the detection rests on one sub-integration, which the authors acknowledge. These are not fatal flaws—the shifts occur at eclipse boundaries and a propagation effect is plausible—but they are gaps that a referee should push on.\n\nThere is also an internal inconsistency: Section 3 derives σ_RM values of 15, 5, and 13 rad m⁻² for the three pulsars, but Section 4.2 quotes 13, 16, and 20, calling them \"measured\" when they come from an arbitrary 5% depolarization assumption. That needs correcting. The quoted uncertainties on ΔPA are statistical only, and the lack of a frequency-resolved analysis is not mentioned in the limitations.\n\nIf the Faraday assumption holds, the sign changes in ΔRM are large compared to the errors, so the reversals are likely real. But the paper does not yet close the sale. It deserves a serious referee because the systems are important and the data are new, but it needs major revision before the B-field and reversal claims can be relied on. My recommendation: send it to peer review, with explicit requests for a λ² test or at least a demonstration that the PA swing shape is constant across the eclipse boundary, and a fix to the σ_RM reporting.","headline":"First direct B-field estimate for B1957+20's eclipse medium and reversal evidence in two spiders, but the core inference rests on untested Faraday assumptions and a fixable internal inconsistency.","tokens_in":12583,"tokens_out":2276,"would_cite":true,"duration_ms":28346,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Using FAST observations, this paper shows that the linear polarization position angle shifts at the eclipse boundary of three black widow pulsars, implying lower-limit line-of-sight magnetic fields of a few to tens of milligauss, with…","keywords":["black widow pulsars","eclipse medium","polarization position angle","Faraday rotation","rotation measure","magnetic field reversal","millisecond pulsars","FAST observations"],"falsifier":"Observe the three pulsars at two or more widely spaced frequencies, such as 600 MHz and 2 GHz, and test whether the eclipse-boundary PA shift scales as $\\lambda^2$ with the same $\\Delta\\mathrm{RM}$; a deviation from the $\\lambda^2$ law would mean part of the shift is intrinsic and the inferred fields are spurious. The single-sub-integration detection for J2055+3829 could also be checked with a longer, higher-signal-to-noise observation.","tokens_in":11458,"feed_emoji":"🕷️","tokens_out":11768,"duration_ms":122180,"temperature":0.7,"pith_summary":"This paper reports FAST polarimetry of three black widow pulsars — PSRs B1957+20, J2055+3829, and J1544+4937 — at 1250 MHz and argues that the eclipse medium around their companions is magnetized at the milligauss level. In all three systems the linear polarization position angle shifts at the eclipse boundary; interpreting the shift as Faraday rotation gives lower-limit line-of-sight magnetic field strengths of a few to tens of milligauss. For PSRs B1957+20 and J1544+4937 the inferred rotation measure changes sign across the eclipse, which the paper takes as evidence of magnetic field reversal in the eclipse medium. These would be the first direct line-of-sight field estimates for B1957+20 and would show that black widow eclipse media contain ordered, spatially structured fields rather than only dense turbulent plasma. The paper also links the observed depolarization to rapid rotation-measure fluctuations, with $\\sigma_{\\rm RM}$ values of 13–20 rad m$^{-2}$, comparable to values inferred for repeating fast radio bursts.","feed_headline":"Black widow pulsar eclipses hide milligauss magnetic fields","feed_subtitle":"Position-angle shifts in three pulsars reveal ordered, sign-reversing fields in the plasma around their companions.","key_machinery":"The mechanism carrying the argument is Faraday rotation of the pulsar signal as it traverses the ionized eclipse medium. The paper measures, for each 30-second sub-integration, the shift $\\Delta\\mathrm{PA}$ in the linear polarization position angle relative to the out-of-eclipse average profile, and converts it to a rotation-measure change using $\\Delta\\mathrm{PA} = \\Delta\\mathrm{RM}\\,\\lambda^2$, where $\\lambda$ is the observing wavelength; the line-of-sight magnetic field then follows from $\\Delta B_{\\parallel} = 1.23\\,\\mu\\mathrm{G} \\times \\Delta\\mathrm{RM}/\\Delta\\mathrm{DM}$. A second element is the depolarization model $L/L_0 = \\exp(-2\\lambda^4 \\sigma_{\\rm RM}^2)$, used to attribute the loss of linear polarization at eclipse boundaries to rapid rotation-measure fluctuations within a sub-integration rather than to changes in the pulsar's intrinsic emission. The eclipse edges are located by Fermi-Dirac fits to the flux density, following the approach of Polzin et al. (2019).","core_discovery":"On the paper's own terms, the central discovery is that the position angle of linearly polarized emission shifts at the eclipse boundary of all three black widow pulsars, and that this shift is best explained as Faraday rotation produced in the eclipse medium. From $\\Delta\\mathrm{PA} = \\Delta\\mathrm{RM}\\,\\lambda^2$, the authors obtain rotation-measure changes of about $-9$ to $+30$ rad m$^{-2}$ for PSR B1957+20 and about $+25$ to $-10$ rad m$^{-2}$ for PSR J1544+4937; combined with a conservative choice of $\\Delta\\mathrm{DM}$ (taken to be the DM measurement uncertainty), these imply lower-limit line-of-sight fields of roughly $-8$ to $+23$ mG for B1957+20 and $+7.3$ to $-4.8$ mG for J1544+4937. For PSR J2055+3829 a single sub-integration gives $\\Delta\\mathrm{RM} = -7 \\pm 1$ rad m$^{-2}$ and a lower-limit field of $-2.2 \\pm 0.3$ mG. The sign changes in $\\Delta\\mathrm{RM}$ for B1957+20 and J1544+4937 are presented as evidence that the magnetic field reverses inside the eclipse medium. The paper explicitly notes that the field-strength values are lower limits because the DM change is not independently resolved, and that the J2055+3829 detection rests on one sub-integration.","pith_inferences":["A direct test of the Faraday-rotation interpretation would be to observe these pulsars at two widely separated frequencies and check that the same $\\Delta\\mathrm{RM}$ reproduces the PA shifts at both; if $\\Delta\\mathrm{PA}$ does not scale as $\\lambda^2$, part of the shift is intrinsic to the pulsar.","The inferred field reversals, if confirmed, suggest that the companion wind and the pulsar wind interact to create current sheets or magnetic loops; time-resolved mapping of $\\Delta\\mathrm{RM}$ across the orbit could constrain the geometry and scale of these structures.","Because the depolarization scatter $\\sigma_{\\rm RM}$ values (13–20 rad m$^{-2}$) overlap those inferred for repeating fast radio bursts, the black widow eclipse medium may serve as a local laboratory for the magnetoionic environments that depolarize FRBs; bright lensed pulses near eclipse, like those seen in PSR B1744-24A, could test this connection."],"forward_implications":["If the field estimates hold, the eclipse media of black widow pulsars are threaded by ordered magnetic fields of order a milligauss or more, not just by dense, unstructured plasma.","The sign reversals seen in PSRs B1957+20 and J1544+4937 imply that the magnetic field structure in the eclipse medium is coherent over at least a fraction of the orbit and changes direction within it.","Depolarization setting in before measurable DM changes points to magnetic-field fluctuations, rather than density fluctuations, as the dominant cause of the loss of linear polarization at eclipse boundaries.","PSR J1544+4937, whose radio emission survives throughout the eclipse at 1250 MHz, offers a way to map the magnetic field across the full eclipse, making it a useful laboratory for eclipse-medium structure."],"supporting_citations":[{"why":"Supplies the method of subtracting the out-of-eclipse PA swing as the reference for measuring ΔPA at eclipse boundaries.","marker":"Crowter et al. (2020)"},{"why":"Gives the depolarization model used to interpret the loss of linear polarization as rapid RM fluctuations.","marker":"You et al. (2018)"},{"why":"Provides the Fermi-Dirac eclipse-fitting procedure and the upper-limit relation for the perpendicular magnetic field.","marker":"Polzin et al. (2019)"},{"why":"Discovered PSR B1957+20 and established the black widow phenomenon of eclipsing millisecond pulsars.","marker":"Fruchter et al. (1988)"},{"why":"Reported the earlier 1400 MHz eclipse duration for B1957+20 that the new measurement is compared with.","marker":"Ryba & Taylor (1991)"},{"why":"Discovered PSR J1544+4937 and measured its eclipses at 322 and 607 MHz, defining the expected eclipse window used here.","marker":"Bhattacharyya et al. (2013)"},{"why":"Discovered PSR J2055+3829 and measured its 10 per cent eclipse duration at 1.4 GHz, the baseline for the new duration.","marker":"Guillemot et al. (2019)"},{"why":"Provided the plasma-lensing constraint on B1957+20's eclipse medium that argues against cyclotron damping at L band.","marker":"Li et al. (2019)"}],"fun_headline_variants":["Pulsar eclipses expose milligauss fields and reversals","Black widow pulsing reveals magnetic reversals in eclipse","Eclipse position-angle shifts reveal pulsing magnetic fields","Field reversals found in black widow pulsar eclipse plasma","Milligauss fields detected in black widow pulsar eclipses"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the position-angle shifts come entirely from Faraday rotation in the eclipse medium, with the pulsar's own position-angle swing unchanged at those orbital phases; if that premise fails, the rotation measures and field strengths are not real.","fun_headline_variants_meta":{"raw":{"variants":["Pulsar eclipses expose milligauss fields and reversals","Black widow pulsing reveals magnetic reversals in eclipse","Eclipse position-angle shifts reveal pulsing magnetic fields","Field reversals found in black widow pulsar eclipse plasma","Milligauss fields detected in black widow pulsar eclipses"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000244,"raw_usage":{"total_tokens":1585,"prompt_tokens":1051,"completion_tokens":534,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":667,"completion_tokens_details":{"reasoning_tokens":452}},"tokens_in":667,"tokens_out":534,"duration_ms":63983,"temperature":1.0,"reasoning_tokens":452,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T00:55:02.463294+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Observe the three pulsars at two or more widely spaced frequencies, such as 600 MHz and 2 GHz, and test whether the eclipse-boundary PA shift scales as $\\lambda^2$ with the same $\\Delta\\mathrm{RM}$; a deviation from the $\\lambda^2$ law would mean part of the shift is intrinsic and the inferred fields are spurious. The single-sub-integration detection for J2055+3829 could also be checked with a longer, higher-signal-to-noise observation.","supporting_citations":[{"cited_title":"F., Taylor J","cited_arxiv_id":null,"evidence_quote":"Reported the earlier 1400 MHz eclipse duration for B1957+20 that the new measurement is compared with."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Discovered PSR J2055+3829 and measured its 10 per cent eclipse duration at 1.4 GHz, the baseline for the new duration."}],"review_version":1}