{"id":"5d4768a0-7e57-4ab4-a985-1c071f06b9fb","arxiv_id":"2505.00184","paper_version":1,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A stacking analysis of 57 millisecond pulsars in HAWC data finds no TeV halo emission, setting upper limits below the inferred efficiency of known TeV halos.","lead":"HAWC searched 2565 days of data for TeV gamma-ray halos around 57 millisecond pulsars and found no significant emission. The null result argues that millisecond pulsars are less efficient than isolated pulsars at producing these halos, affecting how the Galactic Center gamma-ray excess is interpreted.","discovery_kind":"replication","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Efficiency limits assume MSP halos share Geminga's physical radius; if MSP halos are more extended, the Gaussian template captures only a fraction of the flux and the reported limits may be artificially low.","rationale":"The reader identified Eq. (1) as the weakest assumption; this stress-test agrees and sharpens it. The paper's Table IV extension systematic uses a single 0.468 deg template for all sources. That is not a test of the scaling of physical halo radius. For sources with distance greater than ~0.9 kpc, the fixed angular extension corresponds to a larger physical radius than Eq. (1), but for the closest sources (e.g., PSR J0030+0451 at 0.32 kpc, Eq. (1) gives 1.54 deg) the fixed template is much smaller, so the test does not coherently probe the 'larger halos' hypothesis. The point-source stacking in the Appendix provides the compact limit and shows similar flux limits, but a larger-than-Geminga halo is precisely the case where both the extended and point-source templates are mismatched. Quantitatively, for a true Gaussian source of total flux F_true and width sigma_true, fitting a Gaussian template of width sigma_model yields a best-fit total flux K ~ F_true * (sigma_model/sigma_true)^2. Thus the efficiency limits in Figure 2 are limits on the flux within the assumed Geminga-size aperture, not on the total halo luminosity, and can be artificially low by the area ratio if the true halos are larger. This is load-bearing for the central claim, because the paper's own Discussion offers the 'no slow diffusion region around MSPs' scenario, which would naturally produce larger, fainter halos. The proposed test directly settles whether the claim survives a coherent factor-2 or factor-4 rescaling of the physical radius. If the limits are stable, the conclusion stands; if they cross the isolated-pulsar efficiency bands, the conclusion should be conditional on the size assumption. The paper is otherwise internally consistent, with a clear null result and appropriate caveats about spectral index and detector response, so a conditional acceptance with a requested robustness check, rather than a rejection, is the appropriate outcome.","tokens_in":16274,"tokens_out":11910,"duration_ms":126699,"concrete_test":"Repeat the stacking analysis of the Methods section with the Gaussian template radius multiplied by a global factor f = 2 and f = 4 (physical radii ~17 and ~35 pc, respectively) for all 57 MSPs, keeping the same source weights, spectral index, and energy bins. Recompute the 95% confidence upper limits on the common efficiency in the full band (0.32-100 TeV) and in bins 2 (3.16-10 TeV) and 3 (10-31.6 TeV), and compare them to the isolated-pulsar halo efficiency bands in Figure 2. If the rescaled limits lie above those bands, the paper's headline conclusion must be qualified.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim that MSPs are less efficient than isolated pulsars in producing TeV halos rests on the flux limits obtained with Eq. (1), which assigns every MSP halo the same physical radius as Geminga (sigma_source = d_Geminga/d_source * 2.0 deg). For a Gaussian template normalized to unit integral, the fitted normalization K equals the total flux only if the template size matches the true source size. If MSP halos are physically larger than Geminga's (plausible if, as the paper itself argues, MSPs lack the intrabinary activity to generate the slow-diffusion region that confines isolated-pulsar halos), the fitted K for the too-small template is suppressed by approximately (sigma_template/sigma_true)^2 relative to the true total flux. The upper limits in Table I and Figure 2 would then underestimate the true TeV luminosity, and the conclusion 'MSPs are not as efficient as isolated pulsars' would not follow. The paper's extension systematic in Table IV replaces the per-source Eq. (1) values with a single angular extension of 0.468 deg; this is not a coherent rescaling of the physical radius, and for nearby sources it actually shrinks the template, testing the opposite direction. The point-source analysis brackets the compact limit but not the extended limit. Thus the key assumption of the efficiency comparison is untested in the direction that matters.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a search for very-high-energy (TeV) gamma-ray halos around 57 millisecond pulsars (MSPs) using 2565 days of HAWC data. The authors find no significant emission from individual MSPs and, combining likelihood profiles under two weighting schemes (proportional to spin-down flux and to Fermi-LAT GeV flux), find that the stacked signal is consistent with background. They derive upper limits on a common halo efficiency and conclude that MSPs are less efficient than isolated pulsars in producing TeV halos, with implications for the Galactic Center GeV excess and the galactic diffuse emission. The analysis is careful in its treatment of background trials and systematic uncertainties, and it re-examines earlier tentative detections with more data.","tokens_in":16546,"tokens_out":12621,"duration_ms":127521,"significance":"If upheld, this is the first population-level constraint on TeV halos around MSPs and would weaken the case that unresolved MSPs explain the Galactic Center GeV excess via TeV counterparts. The paper is valuable for its null result, its explicit refutation of previous tentative detections, and its discussion of implications. The likelihood methodology is standard and the authors include Monte Carlo background distributions and systematic checks on spectral index, detector response, and source extension. However, the central comparative claim about halo efficiency rests on assumptions about the angular size of MSP halos that are not adequately tested.","major_comments":[{"comment":"The upper limits on the halo efficiency in Table I and Figure 2 are derived from a Gaussian template whose angular size is fixed by Eq. (1) to the same physical radius as the Geminga halo. If MSP halos are physically more extended than Geminga's, the fitted normalization K is biased low by approximately (sigma_template/sigma_true)^2, so the reported limits underestimate the true total flux and efficiency. The 'Extension' systematic in Table IV replaces the per-source angular sizes with the sample mean (0.468 deg), which is not a coherent rescaling of the assumed physical radius and does not test the direction that matters; for nearby sources it actually shrinks the template. The Discussion explicitly suggests that MSPs may lack the slow-diffusion regions that confine isolated-pulsar halos, which would make MSP halos larger. The authors should present limits for a range of assumed physical radii (e.g., 0.5, 1, 2, and 4 times Geminga's radius) and assess whether the conclusion that MSPs are less efficient than isolated pulsars is robust, or alternatively qualify the abstract's unconditional claim.","section":"Methods, Eq. (1); Table IV"},{"comment":"The stacking analysis assumes a single common efficiency eta for all MSPs, with source weights C_i fixed to the spin-down or GeV flux scaling. If only a subset of MSPs (e.g., those with favorable conditions) produce halos with efficiencies comparable to isolated pulsars, the stacked signal could be diluted and the upper limit on the common efficiency would not exclude such a population. The individual-source upper limits partially address this, but they are not combined into a statement about the fraction of MSPs that could be efficient. The authors should either quantify the sensitivity to a fraction f of MSPs contributing with efficiency eta, or restrict the conclusion accordingly.","section":"Results, Table I; Eq. (3)"}],"minor_comments":[{"comment":"The text says the background simulation is repeated 3093 times, while the caption of Figure 5 states the background distribution is composed of 11068 trials; these numbers are inconsistent and should be reconciled.","section":"Appendix: Distribution of the Stacked TS"},{"comment":"The source-count statement reports 57 sources selected from ATNF (53), 3PC (37, 34 overlap with ATNF), WVU (3), and LOTAAS (1); a direct sum gives 60 unique sources if the WVU and LOTAAS entries do not overlap with the others, so the overlapping membership of these catalogs should be clarified.","section":"Methods: source selection"},{"comment":"The paper states the spectral index is set to p=2.7 and that its uncertainty is studied as a systematic, but it does not report the range of p values tested; please specify the tested values in the systematic uncertainty section.","section":"Methods: spectral index"},{"comment":"The sentence beginning 'As the best fit TS of the two weight models over the full energy range is 2.43 and 5.04, which indicates that' is grammatically incomplete and should be rewritten.","section":"Results: point-source stacking"},{"comment":"The entry PSRJ1939+2134 lists RA=19.66 degrees, which is far from the pulsar's actual right ascension; this appears to be a typo and should be corrected.","section":"Table II"}],"recommendation":"major_revision","confidential_remarks":"The paper is a well-executed null-result analysis from the HAWC collaboration, and the refutation of earlier tentative detections is a solid contribution. The main concern is the untested assumption on the halo angular size, which could weaken the central claim about MSP efficiency. I recommend requesting a dedicated systematic study with larger template sizes before publication, and possibly a rephrasing of the abstract to state the conclusion conditionally on the assumed halo morphology."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Hi — quick take on 2505.00184. The HAWC collaboration reports a non-detection of TeV halos around 57 MSPs in 2565 days of data, and their re-analysis of the four sources Hooper & Linden claimed shows all TS < 2. That part is convincing: with roughly 70% more data and better reconstruction, the earlier signals look like fluctuations. The stacking analysis is also clean, with Monte Carlo background trials and systematics on spectral index and detector response. This is a real result worth taking seriously.\n\nThe soft spot is the halo template. Eq. (1) sets each MSP halo to the same physical radius as Geminga (sigma = 8.7 pc/d). The paper’s own argument — that MSPs may lack the slow-diffusion region that confines isolated-pulsar halos — points in the opposite direction: if anything, MSP halos would be more extended, not the same size. For a Gaussian template whose size is underestimated, the fitted normalization K is suppressed by roughly (sigma_template/sigma_true)^2 relative to the total flux, so the reported upper limits on efficiency would be too low. The extension systematic in Table IV replaces per-source sizes with a single 0.468° template; for nearby sources this is smaller than the Geminga-scaled size, so it tests the compact direction, not the extended direction. The point-source analysis also brackets only the compact limit. So the key assumption is not really stress-tested in the direction that matters.\n\nHow much does this hurt? The direct refutation of Hooper & Linden does not depend on the template size, and even a factor-several weakening of the stacking limits would still leave the MSP efficiency below the Geminga-like band in the few-TeV bins. The claim that MSPs are not as efficient as isolated pulsars is stated as a suggestion, and it is probably true in broad strokes, but it is not as airtight as the abstract implies. A systematic that inflates the template by, say, a factor of 2–5 would make the limit robust.\n\nBottom line: this paper deserves a serious referee. The referee should push on the extension assumption, but I would not desk-reject over it. My inclination is to accept after a moderate revision that either tests larger templates or softens the efficiency conclusion.","headline":"Solid HAWC non-detection that convincingly refutes earlier tentative MSP-halo claims; main caveat is the Geminga-sized template, which is not stress-tested in the extended direction.","tokens_in":17615,"tokens_out":5649,"would_cite":true,"duration_ms":61608,"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":"Millisecond pulsars show no TeV halos in a stacked search of 2565 days of HAWC data.","keywords":["TeV halos","millisecond pulsars","HAWC Observatory","gamma-ray stacking analysis","pulsar wind nebulae","Galactic Center GeV excess","inverse Compton emission","very-high-energy gamma rays"],"falsifier":"Observing one millisecond pulsar with a TeV halo whose efficiency matches Geminga's, or a stacked test statistic above the paper's $2\\sigma$ background threshold of ${\\rm TS}=3.99$ in a reanalysis with more exposure, would directly contradict the central claim.","tokens_in":16120,"feed_emoji":"🌌","tokens_out":6726,"duration_ms":67212,"temperature":0.7,"pith_summary":"The paper tries to establish that millisecond pulsars (MSPs) do not power TeV halos as efficiently as isolated, middle-aged pulsars do. Using 2565 days of HAWC Observatory data, it stacks the gamma-ray likelihoods of 57 MSPs under two assumptions about how halo brightness scales (with spin-down power or with GeV brightness) and finds the combined signal consistent with background fluctuations. The 95% confidence upper limits place MSP halo efficiency below the measured efficiency of halos around pulsars like Geminga, especially above 10 TeV. If right, the result also means earlier tentative detections around a few MSPs were statistical fluctuations, and it removes a major astrophysical objection to unresolved MSPs as an explanation of the Galactic Center GeV excess.","feed_headline":"Stacked search finds no TeV halos around millisecond pulsars","feed_subtitle":"Combining 57 pulsars and 2565 days of data caps halo efficiency below Geminga-like levels.","key_machinery":"The load-bearing machinery is a stacking maximum-likelihood search. Each MSP is modeled as a bidimensional Gaussian halo with extension $\\sigma_{\\rm source} = (d_{\\rm Geminga}/d_{\\rm source})\\,2.0^\\circ$, a power-law spectrum of index 2.7, and a common normalization $K$ shared by all sources; source-dependent weights $C_i$ set the relative flux either from spin-down flux or GeV flux. The log-likelihood profiles of all sources are summed to constrain the common efficiency $\\eta$, and the resulting stacked test statistic is compared against a Monte Carlo background built from random sky positions that satisfy the same selection criteria. This machinery lets a population of individually undetectable halos be tested jointly, and converts a null result into upper limits on halo efficiency.","core_discovery":"The central claim is that TeV halos around millisecond pulsars are absent, or at least far less efficient than those around isolated pulsars. The paper reaches this by fitting each of 57 MSPs with a Gaussian halo template whose angular size is scaled from Geminga's $2.0^\\circ$ halo by distance, then adding the individual likelihood profiles with a common normalization for halo flux. In both weighting schemes (flux proportional to spin-down flux $L_{\\rm sd}/(4\\pi d^2)$ and flux proportional to Fermi 0.1-100 GeV flux $G_{100}$), the stacked test statistic matches background expectations in every energy bin from 0.32 to 100 TeV. The resulting 95% confidence upper limits on the efficiency $\\eta = E^2\\Phi/(L_{\\rm sd}/4\\pi d^2)$ or $\\eta_{\\rm GeV} = E^2\\Phi/G_{100}$ lie below the efficiency band of detected TeV halos, particularly above 10 TeV. The paper also re-fits the four MSPs previously reported as tentative detections and finds all have $\\sqrt{\\rm TS}<2$ with the updated data, attributing the earlier excesses to stray Galactic plane emission.","pith_inferences":["If MSPs genuinely lack the slow-diffusion regions that confine electrons around isolated pulsars, the confinement mechanism may require large-scale pulsar wind nebulae or long-lasting particle injection; a targeted search for small halos around the most energetic MSP binaries would test this.","The null result sharpens the puzzle of what accelerates particles around middle-aged pulsars: MSPs have comparable GeV luminosity but apparently do not make TeV halos, so the efficiency difference may point to the wind nebula rather than the magnetosphere as the essential ingredient.","A southern-sky survey with next-generation instruments could test whether the absence is intrinsic to MSPs or an artifact of the northern-sky HAWC sample, and could also probe MSP-rich globular clusters for faint TeV emission.","The stacking technique of scaling a single halo template across many sources and summing likelihoods could be transferred to other faint source populations, such as radio-quiet pulsars or globular clusters, to constrain their TeV emission."],"forward_implications":["MSPs do not contribute significantly to TeV and higher-energy Galactic diffuse emission, leaving unresolved halos of isolated middle-aged pulsars as the viable leptonic explanation at a few TeV.","The MSP interpretation of the Galactic Center GeV excess is strengthened, because MSPs can produce GeV emission without accompanying TeV halos, so previous TeV-halo constraints on that scenario are weakened.","The four tentative MSP halo detections reported in earlier work are consistent with background fluctuations once more data and better event reconstruction are used.","For individual bright MSPs, the current upper limits are approaching but still above the magnetospheric TeV emission models of PSR B1957+20 and PSR J2339-0533, implying that a 25-60% sensitivity improvement could lead to a detection."],"supporting_citations":[{"why":"It supplies the Geminga halo detection and defines the efficiency ratio $\\eta$ used throughout the paper.","marker":"[1]"},{"why":"It provides Geminga's $2.0^\\circ$ Gaussian angular extension and the slow-diffusion interpretation that motivates the halo template.","marker":"[5]"},{"why":"It supports the claim that TeV halos have similar physical sizes, which justifies scaling the angular template by distance.","marker":"[6]"},{"why":"It supplies the Fermi-LAT pulsar catalog from which the MSP sample and the $G_{100}$ GeV fluxes are drawn.","marker":"[9]"},{"why":"It is the earlier work reporting tentative TeV emission around MSPs that this paper re-examines and finds consistent with fluctuations.","marker":"[12]"},{"why":"It is the earlier stacked analysis claiming MSP TeV halos with Geminga-like efficiency, which this paper directly contradicts.","marker":"[13]"},{"why":"It provides the third HAWC catalog used to exclude MSPs near known TeV sources and thereby reduce source confusion.","marker":"[2]"},{"why":"It supplies the 3ML maximum-likelihood framework used for the individual source fits and the stacking analysis.","marker":"[34]"}],"fun_headline_variants":["No TeV halos around millisecond pulsars","Stacked 57 MSPs: no TeV halo excess","Millisecond pulsars produce no TeV halos","MSP TeV halos absent in HAWC stack","TeV halo efficiency of MSPs capped low"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The analysis assumes every millisecond pulsar halo is a Gaussian of the same physical radius as Geminga's halo, with flux proportional either to spin-down power or to GeV brightness; if real MSP halos are smaller, larger, or differently shaped, the stacked search could miss them.","fun_headline_variants_meta":{"raw":{"variants":["No TeV halos around millisecond pulsars","Stacked 57 MSPs: no TeV halo excess","Millisecond pulsars produce no TeV halos","MSP TeV halos absent in HAWC stack","TeV halo efficiency of MSPs capped low"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000225,"raw_usage":{"total_tokens":1472,"prompt_tokens":961,"completion_tokens":511,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":577,"completion_tokens_details":{"reasoning_tokens":433}},"tokens_in":577,"tokens_out":511,"duration_ms":5228,"temperature":1.0,"reasoning_tokens":433,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T04:48:37.117676+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Observing one millisecond pulsar with a TeV halo whose efficiency matches Geminga's, or a stacked test statistic above the paper's $2\\sigma$ background threshold of ${\\rm TS}=3.99$ in a reanalysis with more exposure, would directly contradict the central claim.","supporting_citations":[{"cited_title":"Hooper and T","cited_arxiv_id":null,"evidence_quote":"It is the earlier work reporting tentative TeV emission around MSPs that this paper re-examines and finds consistent with fluctuations."},{"cited_title":"Very-High-Energy Emission From Pulsars","cited_arxiv_id":"2110.09412","evidence_quote":"It is the earlier stacked analysis claiming MSP TeV halos with Geminga-like efficiency, which this paper directly contradicts."},{"cited_title":"Galactic millisecond pul- sars,","cited_arxiv_id":null,"evidence_quote":"It supplies the 3ML maximum-likelihood framework used for the individual source fits and the stacking analysis."}],"review_version":1}