{"id":"1ba89f0e-ea91-4eb4-ae88-f900976c66ec","arxiv_id":"2509.09605","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":9,"one_line_summary":"PSR J1544-2555 is a new 2.39 ms black-widow pulsar in a 2.7-hour orbit, identified via optical variability and confirmed by radio and gamma-ray pulsations.","lead":"Astronomers discovered a new black-widow millisecond pulsar, PSR J1544-2555, by spotting its 2.7-hour orbital signal in optical data and then confirming it with radio and gamma-ray observations. The find expands the known population of these extreme binary systems and demonstrates a promising pathway for discovering more.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Spot-model optical fit drives distance and derived masses into tension with the DM distance; this is the main condition on the interpretation but does not threaten the pulsar discovery.","rationale":"I agree with the reader's weakest assumption. The central claim—that PSR J1544-2555 is a newly discovered black-widow MSP—is robust: the radio pulsations establish the 2.39 ms spin period, the optical period and radio timing jointly fix the 0.1135 d orbit, and the gamma-ray detection is independently significant (H=76.3 even before the OPV fit). The spot-model distance tension does not threaten the pulsar's existence or its classification as a compact spider binary, but it does affect the inferred companion mass, inclination, and X-ray luminosity, which are secondary derived quantities. The paper itself notes the absence of radial velocities, the only external anchor for the heating models. The missing OPV code release is a secondary reproducibility concern, but it does not change the verdict because the discovery does not hinge on the final H=627.8 value. The CONDITIONAL verdict is therefore appropriate, and my stress-test does not alter it.","tokens_in":21782,"tokens_out":10183,"duration_ms":120354,"concrete_test":"Perform high-S/N phase-resolved optical spectroscopy of the companion to measure its radial-velocity semi-amplitude K_c. The spot model predicts K_c=470±20 km s^-1 (§6.2, Table 4), while the direct-heating/convection models predict 300–400 km s^-1. If K_c is measured outside ~450–490 km s^-1, the spot-model inclination/distance solution is excluded and the distance conflict resolves in favor of the DM estimate; if K_c=470 km s^-1, the spot model and its d≈2 kpc are independently supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The discovery of PSR J1544-2555 as a black-widow MSP is secure: radio pulsations at 2.39 ms, a 2.7-hr optical orbital modulation, and gamma-ray pulsations independently support the binary nature. The load-bearing conditional is the Icarus spot-model interpretation (§6.2). The spot model adds four free parameters (longitude, colatitude, radius, temperature) to direct heating, and although it is claimed to improve the fit, no Δln Z or equivalent model-comparison statistic is reported; the improvement could be overfitting of the light-curve asymmetry. Under this model, the distance is 2.0^{+0.1}_{-0.3} kpc, in tension with the DM-derived 1.02–1.06 kpc (Table 3). It also changes inclination from the 46–50° range of the other models to 65^{+6}_{-5}°, companion mass to 0.095 M⊙, and thus the X-ray luminosity (§5). Since no radial-velocity measurement anchors K_c, the mass/inclination/distance solution is not externally verified. The paper admits this in §7. This is a real conditional, but not a fatal one: even if the spot model is wrong, the pulsar and its orbital period stand.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports the discovery of PSR J1544−2555, a black-widow millisecond pulsar associated with the Fermi-LAT source 4FGL J1544.2−2554. ULTRACAM photometry revealed a ~2.7-hour optical modulation; follow-up MeerKAT observations detected 2.39-ms radio pulsations. A joint radio and gamma-ray timing solution, spanning 16 years of Fermi-LAT data, yields a precise orbital ephemeris with H=627.8 and shows orbital-period variations typical of spider pulsars. eROSITA X-ray data indicate non-thermal emission, and Icarus light-curve modeling favors a spot model over direct-heating and heat-redistribution models, giving a distance of ~2.0 kpc, an inclination of ~65°, and a companion mass of ~0.095 M_sun. The authors acknowledge that the lack of radial-velocity measurements leaves the neutron-star mass unconstrained and that the spot model is one of several possible interpretations.","tokens_in":22185,"tokens_out":2992,"duration_ms":35767,"significance":"If the discovery claims hold, this is a valuable new addition to the small population of black-widow pulsars and a demonstration of an effective discovery path from Fermi unIDs through optical variability to radio confirmation. The core detection is robust: the 2.4-ms radio pulsations, the 2.7-hour optical modulation, and the gamma-ray pulsations are mutually independent and jointly establish the binary nature. The timing solution spans 16 years with high significance and is a clear strength. The optical light-curve modeling is more speculative, but the paper's main discovery does not depend on it. The paper also makes useful comparisons with the independent K24 study and places the source in the context of spider-pulsar populations.","major_comments":[{"comment":"The claim that the spot model 'exhibits a significantly improved goodness of fit' and is favored over the DH, C, and D+C models is not supported by any reported model-comparison statistic. The text mentions 'slightly better Bayesian evidence' for the heat-redistribution models but gives no Δln Z, BIC, or equivalent. Since the spot model changes the derived distance from ~1 kpc to ~2 kpc, the inclination from 46–50° to 65°, and the companion mass, the evidence for this extra four-parameter model is load-bearing for Table 4 and §5. Please report the Bayesian evidence (or AIC/BIC) for all four models, along with posterior predictive checks or residuals, so the reader can judge whether the improvement is not just overfitting the light-curve asymmetry.","section":"§6.2"},{"comment":"The X-ray luminosity L_X ≈ 4.0×10^31 erg/s is computed using the spot-model distance of 2.0 kpc, but the timing-derived DM distance is 1.02–1.06 kpc (Table 3). At 1.04 kpc the same flux would give L_X ≈ 1.0×10^31 erg/s, changing the stated X-ray efficiency by a factor of four. This is a systematic uncertainty that should be stated explicitly. Please give L_X for both the spot-model distance and the DM-based distance, or at least note that the quoted efficiency is conditional on the spot-model distance.","section":"§5"},{"comment":"The distance prior is built from the YMW16 DM distance of 1.02 kpc and the Galactic MSP distribution, yet the spot-model posterior is 1.7–2.0 kpc (Table 4). The paper notes this exceeds the DM distances but does not discuss whether this indicates a problem with the spot model, the DM distance, or the Galactic electron-density model. Since the distance enters the X-ray luminosity and other derived quantities, a short quantitative discussion (e.g., a posterior predictive check or a comparison of the DM distance with the modeled distance under different assumptions) would strengthen the interpretation.","section":"§6.1–6.2"}],"minor_comments":[{"comment":"The caption says '0.25 corresponds to the pulsar's superior conjunction'; please define the phase convention explicitly in the text (phase 0 = ascending node) so the reader can connect the optical and radio ephemerides.","section":"§2, Figure 1"},{"comment":"The abbreviation 'K24' is used for Karpova et al. (2024); please define it at first use in the text or in the comparison section, since it is not a standard abbreviation.","section":"§6.2"},{"comment":"Typo: 'disentagle' should be 'disentangle'.","section":"§7"},{"comment":"The 'detection likelihood around 20' is not defined; please state whether this is a binned Poisson likelihood, a maximum-likelihood ratio, or another statistic, and give the corresponding significance in Gaussian sigma if possible.","section":"§5"},{"comment":"The tidal-locking assumption (ω=1) is stated but not tested. Since no radial velocities are available, the derived masses and inclinations depend on this assumption. A brief comment on how non-synchronous rotation would affect the results would be useful.","section":"§6.1"}],"recommendation":"minor_revision","confidential_remarks":"The discovery of the pulsar is secure and well presented. The main revision request is to add quantitative model-comparison evidence for the spot model and to make the distance dependence of the X-ray luminosity explicit. Both are local fixes; the paper does not require new observations to address them. The manuscript fits the scope of MNRAS and will be of interest to the spider-pulsar and Fermi-LAT communities."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth your time. PSR J1544-2555 is a real new black-widow MSP: 2.39-ms radio pulsations from MeerKAT, a clean joint radio+gamma-ray timing solution spanning 16 years with H=627.8, and optical orbital modulation at 2.7 hours. The discovery path — optical period search first, then targeted radio follow-up — is well executed and the first time a black-widow has been confirmed this way. The OPV measurement with a Matérn Gaussian process is also a nice piece of work, and the paper is honest about its own limitations.\n\nThe soft spots are concentrated in the optical light-curve modeling, and the stress-test note gets this exactly right. The spot model improves the fit, but no model-comparison statistic (Δln Z or similar) is reported, so \"significantly improved goodness of fit\" is asserted rather than demonstrated. More importantly, the spot model drives the distance to 2.0 kpc, in tension with the DM-derived 1.0–1.1 kpc, and shifts the inclination from 46–50° to 65° and the companion mass to 0.095 M⊙. Without a radial-velocity anchor for Kc, these derived values are conditional on the spot model being correct. The X-ray luminosity in Section 5 inherits that condition because it uses the spot-model distance. The paper admits this in Section 7, which is good, but the abstract and conclusions lean harder on the spot interpretation than the evidence supports. Also, the OPV timing code is not yet released, which limits independent verification — though the authors say it will come.\n\nNone of this threatens the discovery. The pulsar exists, the orbital period is secure, and the system is clearly a black-widow. The modeling questions are open, but they are exactly the kind of thing peer review should probe.\n\nWho should read this: anyone working on spider pulsars, especially the discovery/classification pipeline and OPV methodology. It deserves a serious referee, not a desk reject. I'd send it to review with one request: quantify the evidence for the spot model over the simpler alternatives, and discuss the DM-distance tension explicitly, even if no resolution is available. Code release would be a plus but I wouldn't hold acceptance on it.\n\nI'd cite this for the timing solution and the discovery, and I'd bring it to reading group.","headline":"A genuinely new, securely confirmed black-widow MSP with a solid timing package; the spot-model interpretation is the main conditional and should be flagged, but the discovery itself is not in doubt.","tokens_in":22720,"tokens_out":1219,"would_cite":true,"duration_ms":15602,"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 the gamma-ray source 4FGL J1544.2−2554 is a new black-widow millisecond pulsar, PSR J1544−2555, a 2.39-ms radio pulsar in a 2.7-hour orbit with a low-mass, strongly heated companion.","keywords":["millisecond pulsar","black-widow pulsar","spider binary","gamma-ray pulsation","orbital period variation","optical light curve","X-ray emission","radio timing"],"falsifier":"Measure the companion's radial velocity curve with high-resolution spectroscopy: the spot model predicts a semi-amplitude of roughly 450–500 km/s, while the simpler heating models predict 300–400 km/s; a measured value near 400 km/s with small uncertainty would rule out the spot model's inclination and mass, and a value far outside 300–500 km/s would call the entire light-curve interpretation into question.","tokens_in":21711,"feed_emoji":"🕷️","tokens_out":3451,"duration_ms":41274,"temperature":0.7,"pith_summary":"The paper reports the discovery and multiwavelength characterization of a new black-widow millisecond pulsar, PSR J1544−2555, identified with an unassociated gamma-ray source. Optical monitoring revealed a 2.7-hour periodic variation, radio observations then detected 2.39-ms pulsations, and a joint 16-year gamma-ray and radio timing solution confirmed the binary and tracked orbital period variations. Light-curve modeling shows a tidally distorted companion with a localized hot spot, and X-ray data indicate non-thermal emission. The discovery validates a search strategy that uses optical variability to guide targeted radio and gamma-ray follow-up of unidentified gamma-ray sources.","feed_headline":"A new 2.39-ms pulsar hides in a 2.7-hour binary","feed_subtitle":"Optical flickers led to radio and 16-year gamma-ray detections, adding a fresh black widow to the spider pulsar family","key_machinery":"The argument is carried by a joint radio and gamma-ray timing solution that spans 16 years, where orbital period variations are modeled as a Gaussian process with a Matérn covariance function, plus a binary light-curve synthesis model that treats the companion as a tidally locked, irradiated star with an optional localized Gaussian hot spot. The timing solution fixes the spin, orbital, and astrometric parameters, while the light-curve model converts the observed optical modulation into inclination, companion mass, temperatures, and distance. The spot model is the key new element that captures the light-curve asymmetry and shifts the inferred inclination from about 48 to 65 degrees.","core_discovery":"PSR J1544−2555 is a genuine black-widow millisecond pulsar: a 2.390-ms pulsar in a 0.1135-day (about 2.7-hour) circular orbit with a companion of roughly 0.06–0.12 solar masses. The identification rests on three independently consistent pillars: radio pulsations detected in a targeted search, a 16-year gamma-ray timing solution with high significance that includes orbital period variations characteristic of spider pulsars, and optical light curves showing a strongly irradiated, tidally locked companion. The optical asymmetry is best matched by a model with a localized hot spot near the companion's pole, and a slightly bluer color at the minimum suggests possible non-thermal emission from an","pith_inferences":["If the spot model is correct, the light-curve distance of about 2 kpc conflicts with the dispersion-measure distance of about 1 kpc; this discrepancy, if not a modeling artifact, implies that electron-density models along this line of sight underestimate the distance or that the system has unusual properties, which would affect the derived spin-down luminosity and gamma-ray efficiency.","A single-epoch radial velocity measurement of the companion would immediately break the degeneracy between the spot model and the simpler heating models, providing a direct test of the paper's favored geometry and a way to measure the neutron star mass.","The bluer color at the orbital minimum, if confirmed as synchrotron emission from an intra-binary shock, would make this system a useful laboratory for studying how pulsar winds interact with low-mass companions; a deeper X-ray spectrum could verify this interpretation.","The successful optical-first discovery path could be applied to the hundreds of remaining unidentified gamma-ray sources in the catalog, potentially yielding more black widows and redbacks without the need for wide-field radio pulsation surveys."],"forward_implications":["PSR J1544−2555 becomes a new member of the black-widow pulsar population, adding to the census of compact spider binaries that can constrain neutron-star masses and the equation of state of dense matter.","The discovery demonstrates that optical periodicity searches of unidentified gamma-ray sources can efficiently select black-widow candidates that are then confirmed by short, targeted radio observations.","The 16-year timing solution, including orbital period variations, provides a new data point for studying the orbital dynamics and companion-driven mechanisms common to spider pulsars.","The inferred non-thermal X-ray emission, with an X-ray to spin-down luminosity ratio of about 1.4 times 10 to the minus 3, is consistent with the range seen in other millisecond pulsars and supports an intra-binary shock origin."],"fun_headline_variants":["New black-widow pulsar: 2.39-ms spin, 2.7-hour orbit","Optical flickers reveal 2.39-ms pulsar in a 2.7-hour binary","16-year gamma-ray timing exposes new spider pulsar","Spot model explains black-widow pulsar's asymmetric light curve"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The derived distance, companion mass, inclination, and X-ray luminosity all depend on the light-curve modeling assumption that the companion is tidally locked and that the asymmetric light curve is produced by a single localized Gaussian hot spot; if that geometry is wrong, those numbers shift substantially, although the pulsar's existence and binary nature do not.","fun_headline_variants_meta":{"raw":{"variants":["New black-widow pulsar: 2.39-ms spin, 2.7-hour orbit","Optical flickers reveal 2.39-ms pulsar in a 2.7-hour binary","16-year gamma-ray timing exposes new spider pulsar","Spot model explains black-widow pulsar's asymmetric light curve"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000419,"raw_usage":{"total_tokens":2065,"prompt_tokens":887,"completion_tokens":1178,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":631,"completion_tokens_details":{"reasoning_tokens":1089}},"tokens_in":631,"tokens_out":1178,"duration_ms":12653,"temperature":1.0,"reasoning_tokens":1089,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T18:44:08.137908+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the companion's radial velocity curve with high-resolution spectroscopy: the spot model predicts a semi-amplitude of roughly 450–500 km/s, while the simpler heating models predict 300–400 km/s; a measured value near 400 km/s with small uncertainty would rule out the spot model's inclination and mass, and a value far outside 300–500 km/s would call the entire light-curve interpretation into question.","supporting_citations":[],"review_version":1}