{"id":"454e99f4-6eaf-4e6a-a69e-69b242231c40","arxiv_id":"1908.00992","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"Multiband light curves of nine black-widow pulsar companions are fitted with a direct-heating model, giving distances, inclinations, and heating powers that often exceed the observed gamma-ray luminosities, suggesting beaming corrections, distance errors, or missing physics.","lead":"Astronomers measured how the optical brightness of nine 'black-widow' pulsar systems changes as their tiny companion stars are heated and evaporated. The brightness changes were fit with a simple heating model to estimate distances, orbital tilts, and heating power, revealing that several systems need more heating than the gamma-rays we detect.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The fitted distances are not robust against the model's d–f_c–metallicity degeneracy, and since ηγ scales as d², the claimed ηγ–ηH discrepancies may be a distance artifact rather than evidence of beaming or inclination bias.","rationale":"The reader's weakest_assumption correctly identifies direct-heating-only emission as the load-bearing premise. My stress-test finds a more specific and arguably more immediate vulnerability: even granting the direct-heating model, the fitted distances are demonstrably non-robust, and the paper's own examples (J0952's metallicity-driven distance collapse, J1124's ηγ = 2.5 at 2.5× DM, J0251's 2.8× DM distance) show that distance bias alone can produce the efficiency discrepancies the paper attributes to beaming or unmodeled heating. Since the GeV efficiency depends on d² while the heating efficiency is distance-independent, any systematic overestimate of d by the light-curve fit inflates ηγ and creates the observed scatter. The paper acknowledges this in part, but the central claim includes the fitted distances themselves, so the concern is load-bearing rather than peripheral. A refit with independent distance priors is the direct way to decide whether the distances and efficiencies survive. This does not change the reader's CONDITIONAL verdict: it strengthens the case for caution, but the original verdict already captures that the quantitative results should be used with caution.","tokens_in":19748,"tokens_out":2934,"duration_ms":31749,"concrete_test":"Re-fit the full sample with independent distance priors: fix d to the DM/timing-parallax values (e.g., J0636 from Arzoumanian et al. 2018) and marginalize over log[Z] with an uninformative prior, then compare the posterior distributions of i, f_c, LX, and ηH with Table 3. If the posteriors shift by more than the quoted 1σ errors, or if the ηγ–ηH scatter is substantially reduced, the headline distance/efficiency discrepancies are dominated by model degeneracy rather than astrophysical beaming.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central inference—that direct heating reproduces the light curves and that residual GeV/heating discrepancies trace beaming or unmodeled effects—depends on the ICARUS fits' distance, inclination, and heating power being trustworthy. That trust is not established. The fits are degenerate along d–f_c and sensitive to metallicity, as the paper itself shows: for J0952, log[Z] = −0.5 drives d to unphysically small values while log[Z] = 0 gives d ≈ 4 kpc (§4.2); for J1124, the best-fit d = 2.72 kpc is 2.5× the DM estimate and yields ηγ = 2.5, which the authors call unreasonably large; J0251 similarly fits at 2.8× its DM distance. Because ηγ ∝ d² while ηH is fixed by the model LX, any distance bias directly creates the efficiency mismatches that are attributed to imperfect beaming or unmodeled heating. The paper has only one independent distance check (J0636), and even that fit required an ad hoc Δφ = 0.02 phase shift to reach acceptable χ². Thus the central quantitative claims—distances, inclinations, heating efficiencies, and the resulting neutron-star-mass implications—are conditional on an untested degeneracy-breaking that the current data cannot provide.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents new and archival multiband optical photometry for nine black-widow millisecond pulsar binaries, including first-time optical counterparts for four systems, and fits the folded light curves with the ICARUS direct-heating model. The fits yield companion fill factors, inclinations, heating luminosities, night-side temperatures, distances, and preferred metallicities. The heating powers are then compared with Fermi-LAT GeV luminosities after applying outer-gap beaming corrections, expressed as efficiencies eta_gamma and eta_H. The paper finds that direct heating can broadly reproduce the light curves, but that reduced chi-square values are large for several objects, some fitted distances disagree strongly with dispersion-measure distances, and the eta_gamma-eta_H comparison shows only a mild correlation, which the authors attribute to distance errors, imperfect beaming corrections, and/or unmodeled heating effects.","tokens_in":20107,"tokens_out":3216,"duration_ms":36100,"significance":"If the reported fits are reliable, the paper provides a valuable homogeneous sample of black-widow companion parameters, including evidence for substantial Roche-lobe filling, low inclinations, and a possible high-metallicity preference, with implications for companion evolution and neutron-star mass measurements. The comparison with independently measured Fermi-LAT fluxes is a genuine external check and is not circular. The paper is also commendably transparent about its limitations, explicitly discussing degeneracies, phase shifts, flares, and the possibility of unmodeled effects. However, the central quantitative claims—distances, inclinations, heating efficiencies, and the resulting efficiency discrepancies—are currently conditional on trusting ICARUS fits that show high chi-square and poorly constrained distance-metallicity-fill-factor degeneracies. The significance of the paper will increase substantially if the distance-degeneracy issue is quantified and the efficiency comparison is shown to be robust to it.","major_comments":[{"comment":"The outer-gap beaming corrections are computed with a model designed for young pulsars, and Fig. 8 shows that the observed gamma-ray peak separations Delta_gamma would require i > 60 degrees for most sources, in tension with the photometric fit inclinations in Table 3 (mostly 40-60 degrees). The corrections fi and f* in Table 3 are averages over alpha and are consequently close to unity, which understates the model dependence of the comparison. The paper should test the sensitivity of the eta_gamma-eta_H conclusions to alternative beaming prescriptions, or at least state clearly that the efficiency comparison is contingent on the validity of the young-pulsar OG model applied to black widows.","section":"Sec. 4.1, Sec. 4.3, Fig. 8"},{"comment":"The reduced chi-square values are large for several fits (5.41 for J0023, 3.36 for J0636, 2.64 for J2241), and the fitting procedure includes post-hoc adjustments: J2052's MDM errors are doubled, J0952 flare points are excluded, J0636 requires an ad hoc phase shift of Delta_phi = 0.02 to reduce chi-square by about 200, and J1124 has an apparent broad maximum not matched by the model. The MCMC parameter uncertainties in Table 3 therefore do not capture the full systematic uncertainty in the model parameters. The central claims would be strengthened by an explicit systematic error budget or by demonstrating that the main conclusions persist when the most discrepant points are treated as upper limits or when the post-hoc adjustments are not applied.","section":"Sec. 4, Table 3"}],"minor_comments":[{"comment":"There are several typographical errors in the reference list: 'Antionadis' should be 'Antoniadis', 'Kieth' should be 'Keith', and 'Tiddaren' is used inconsistently in the text.","section":"References, Sec. 2"},{"comment":"The text says J1124 was found in the 'GBT 350 GHz drift-scan survey'; this should presumably be 350 MHz, not 350 GHz.","section":"Sec. 2, Table 1"},{"comment":"The J2052 u-band points are shown in Fig. 4 but are stated not to be fitted; the caption could be clearer that these are plotted for comparison only.","section":"Sec. 3, Table 2"},{"comment":"The color maps in Fig. 8 may be difficult to read in grayscale; consider adding contour labels or explicit numerical values for the key f_Omega levels.","section":"Fig. 8"},{"comment":"The Shklovskii-correction formula in Sec. 4.1 mixes symbol definitions (e.g., mu_mas/yr, d_kpc) without a full variable list; a short definition of each symbol would improve readability.","section":"Sec. 4.1"}],"recommendation":"major_revision","confidential_remarks":"The paper presents useful observational material and a thorough, honest exploration of the direct-heating model for black-widow companions. The main barrier to acceptance is the distance-degeneracy issue: because eta_gamma scales as d^2, the paper's central efficiency comparison is not yet load-bearing. I believe this can be fixed within the manuscript's scope by adding a quantitative treatment of the d-f_c-metallicity covariance and by softening or re-framing the conclusions that depend on the absolute distance scale. The paper is otherwise suitable for the journal and should not be rejected."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing you should know: this is the paper to cite for homogeneous multiband light curves of nine black-widow pulsars, and it presents four new companions. It is also more honest about its own problems than most papers in this area. The authors fix a flux-normalization error in their earlier code and explicitly flag where the fits are unreliable. The photometric dataset alone justifies reading it.\n\nWhat is genuinely new: four optical counterparts (J0251, J1124, J2052, J2241) appear here for the first time. The multicolor light curves are assembled on a consistent SDSS-based system, and the fits use a corrected version of ICARUS to give distances, inclinations, fill factors, and heating powers for the full sample. The comparison of heating power to Fermi-LAT GeV luminosity, with outer-gap beaming corrections, is a reasonable way to test the direct-heating picture, and it is not circular: the gamma-ray fluxes are external.\n\nThe soft spots are the usual ones for this kind of work, and the authors name most of them. The fits rely on the assumption that direct photon heating alone produces the optical modulation. The sharp step in J0023, the phase shift in J0636, and the flaring in J0952 all suggest unmodeled intrabinary-shock or nonthermal emission. With reduced chi-squares up to 5.4, the statistical case for the model is weak. More importantly, the fitted distances are degenerate with fill factor and metallicity; J0952's distance swings between 1.7 and 4 kpc depending on the adopted log[Z], and J1124's best fit is 2.5 times its DM distance, giving an implausible gamma-ray efficiency. Since eta_gamma scales as d^2, the claimed discrepancies between gamma-ray and heating efficiencies can be mostly a distance artifact. The stress-test note is right about that, though the paper itself already concedes it—it even uses the efficiency mismatch to argue that some distances are wrong. So the central conclusion is not \"the beaming model fails\"; it is \"the data cannot yet break the degeneracy between distance and model details.\"\n\nShould this be published? Yes. The observations are valuable, the analysis is transparent, and the qualitative result—that direct heating roughly works but unmodeled effects are present—is reasonable. A good referee should push for a clearer separation of fitted parameters from the gamma-ray comparison, for robustness checks against the distance degeneracy, and ideally for the release of the reduced light curves so others can refit. But desk rejection would be a mistake. For your own work, cite it for the data and the distances with caution.","headline":"Four new black-widow companions and a homogeneous multiband dataset, wrapped in an honest but degenerate direct-heating fit; cite for the data, not for the distances.","tokens_in":20609,"tokens_out":3164,"would_cite":true,"duration_ms":36255,"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":"Direct pulsar heating can explain the optical light curves of nine black-widow binaries, though only with imperfect beaming corrections.","keywords":["black-widow pulsars","millisecond pulsars","optical light curves","direct heating model","ICARUS","gamma-ray beaming","Roche-lobe filling","companion evolution"],"falsifier":"Measure a VLBI or timing parallax distance for J1124-3653 and J0251+2606; if the true distances match the dispersion-measure estimates rather than the ICARUS fit distances, the large fitted distances and the implied gamma-ray efficiencies above unity would be wrong, pointing to unmodeled heating or beaming rather than distance error.","tokens_in":19564,"feed_emoji":"🔭","tokens_out":6723,"duration_ms":61132,"temperature":0.7,"pith_summary":"This paper collects new and archival multiband optical photometry for nine black-widow millisecond pulsar binaries and fits each folded light curve with a direct photon-heating model. The fits return distances, Roche-lobe fill factors, inclinations, and heating powers for the companions. Comparing those heating powers with Fermi-LAT GeV luminosities, after applying outer-gap beaming corrections, leaves several systems substantially inconsistent, which the authors attribute to imperfect beaming corrections, distance errors, or heating effects not captured by the model. The fits also indicate that the companions mostly fill their Roche lobes, that several are metal-rich, and that J0952's companion is notably compact and dense.","feed_headline":"Direct heating fits nine black-widow pulsar light curves","feed_subtitle":"Multiband photometry gives distances, inclinations, and heating powers; comparing with Fermi-LAT gamma rays exposes beaming and distance…","key_machinery":"The load-bearing tool is the ICARUS light-curve synthesis code, which assumes the pulsar's bolometric luminosity propagates linearly from a point source and heats the companion's surface, with model-atmosphere spectra (BTSettl) folded through broadband filters. The model parameters are heating luminosity $L_X$, distance $d$, Roche-lobe fill factor $f_c$, underlying temperature $T_N$, and inclination $i$; the fit uses Markov Chain Monte Carlo to map the minimum. The comparison side uses outer-gap beaming factors $f_i$ and $f_*$ derived from Romani & Watters (2010) to convert Earth-directed Fermi-LAT fluxes and companion-directed heating fluxes to isotropic equivalents, with gamma-ray efficiency $\\eta_\\gamma$ and heating efficiency $\\eta_H$ defined relative to Shklovskii-corrected spindown power.","core_discovery":"The central claim is that a simple direct-heating model, in which a bolometric luminosity from the pulsar irradiates the companion's surface, can reproduce the multiband light curves of all nine black-widow systems well enough to extract physical parameters. The paper finds distances that often exceed dispersion-measure estimates, inclinations concentrated at intermediate values, and heating powers that are a substantial fraction of the spindown luminosity. When the heating power is compared with the Fermi-LAT GeV luminosity through outer-gap beaming corrections, some systems agree but others are discrepant by factors of several, and the paper argues that the discrepancy is a signature of incomplete beaming models, overestimated fit distances, or unmodeled heating channels such as intrabinary shocks. The authors state explicitly that the direct-heating fits imply large heating powers and rather small inclinations and speculate that unmodeled effects bias these quantities.","pith_inferences":["I infer that the $\\eta_\\gamma$ versus $\\eta_H$ comparison could be turned into an independent distance estimator for black-widow systems: fixing the beaming model, matching the two efficiencies would set $d$, and this could be checked against VLBI or timing parallaxes.","This suggests that if direct-heating fits systematically underestimate inclination, neutron-star masses derived from these systems could be overestimated, which would affect equation-of-state constraints.","A testable extension is simultaneous multiband optical and X-ray monitoring across the orbital phase, which should reveal whether the J0023 step and similar anomalies track an intrabinary shock rather than surface heating.","Comparing the same data with alternative gamma-ray beaming models, instead of only the outer-gap model, would show whether the discrepancies are model-specific."],"forward_implications":["If the fits are right, these nine systems have companions that mostly fill their Roche lobes with $f_c \\approx 0.7$–$1$, and J0952's companion is underfilled with density $\\rho \\approx 10\\,\\mathrm{g\\,cm^{-3}}$, implying unusual evolution.","The inferred heating powers imply gamma-ray efficiencies $\\eta_\\gamma$ and heating efficiencies $\\eta_H$ that disagree for several sources, so the comparison can act as a distance diagnostic when the beaming model is better constrained.","J0636's low inclination and very faint gamma-ray emission are consistent with outer-gap beaming away from Earth, providing a direct test of the beaming picture.","The preference for super-solar metallicity in about half the sample supports the idea that these companions are wind-stripped residues.","The J0023 light curve shows a sharp step near maximum that the thermal model cannot produce, suggesting a nonthermal or shock component contaminates some light curves."],"supporting_citations":[{"why":"Supplies the ICARUS light-curve code that the paper uses for all direct-heating fits.","marker":"Breton et al. 2012"},{"why":"Provided the original J0023 GMOS photometry that the paper remeasures and includes in its fits.","marker":"Breton et al. 2013"},{"why":"Provided the multiband photometry of J1959 that the paper re-fits and compares with earlier inclination estimates.","marker":"Reynolds et al. 2007"},{"why":"Supplies J1959's companion radial velocity and the high pulsar-mass argument that the revised inclination and distance bear on.","marker":"van Kerkwijk et al. 2011"},{"why":"Supplies the J1301 photometry and radial-velocity amplitude that anchor that system's fit.","marker":"Romani et al. 2016"},{"why":"Provides the outer-gap beaming factors used to correct Fermi-LAT and heating fluxes to isotropic equivalents.","marker":"Romani & Watters 2010"},{"why":"Describes direct-heating model extensions and previous distance estimates that this paper's corrected distances supersede.","marker":"Sanchez & Romani 2017"},{"why":"Gives the timing-parallax distance for J0636 that the photometric fit is checked against.","marker":"Arzoumanian et al. 2018"}],"fun_headline_variants":["Heating fits for nine black-widow pulsars expose beaming gaps","Nine black-widow pulsar light curves fit by direct heating model","Heating fits clash with gamma rays in black-widow pulsars","Dense companion spotlights odd black-widow evolution","Direct heating fits nine pulsars, but gamma-ray data disagree"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The companion's optical brightness is caused only by direct photon heating of the kind ICARUS computes, with intrabinary shocks, nonthermal emission, and heat transport too weak to shift the fitted inclination, distance, and heating power.","fun_headline_variants_meta":{"raw":{"variants":["Heating fits for nine black-widow pulsars expose beaming gaps","Nine black-widow pulsar light curves fit by direct heating model","Heating fits clash with gamma rays in black-widow pulsars","Dense companion spotlights odd black-widow evolution","Direct heating fits nine pulsars, but gamma-ray data disagree"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000756,"raw_usage":{"total_tokens":3402,"prompt_tokens":1031,"completion_tokens":2371,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":647,"completion_tokens_details":{"reasoning_tokens":2279}},"tokens_in":647,"tokens_out":2371,"duration_ms":16582,"temperature":1.0,"reasoning_tokens":2279,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T15:26:39.767517+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure a VLBI or timing parallax distance for J1124-3653 and J0251+2606; if the true distances match the dispersion-measure estimates rather than the ICARUS fit distances, the large fitted distances and the implied gamma-ray efficiencies above unity would be wrong, pointing to unmodeled heating or beaming rather than distance error.","supporting_citations":[{"cited_title":"P ., Rappaport, S","cited_arxiv_id":null,"evidence_quote":"Supplies the ICARUS light-curve code that the paper uses for all direct-heating fits."},{"cited_title":"T., Callanan, P","cited_arxiv_id":null,"evidence_quote":"Provided the multiband photometry of J1959 that the paper re-fits and compares with earlier inclination estimates."},{"cited_title":"W., & Watters, K","cited_arxiv_id":null,"evidence_quote":"Provides the outer-gap beaming factors used to correct Fermi-LAT and heating fluxes to isotropic equivalents."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Describes direct-heating model extensions and previous distance estimates that this paper's corrected distances supersede."}],"review_version":1}