Multiband Optical Light Curves of Black-Widow Pulsars

Alexei V. Filippenko; Fernando Camilo; Jules P. Halpern; Paul Draghis; Roger W. Romani; Thomas G. Brink; Weikang Zheng

REVIEW 2 major objections 5 minor 1 cited by

Multiband Optical Light Curves of Black-Widow Pulsars

T0 review · 2 major / 5 minor · reviewed 2026-08-14 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Direct pulsar heating can explain the optical light curves of nine black-widow binaries, though only with imperfect beaming corrections.

desk verdict 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. read the letter →

arxiv 1908.00992 v2 pith:ZSSAKPQJ submitted 2019-08-02 astro-ph.HE

Paul Draghis , Roger W. Romani , Alexei V. Filippenko , Thomas G. Brink , WeiKang Zheng , Jules P. Halpern , Fernando Camilo This is my paper · ORCID

classification astro-ph.HE

keywords black-widowpulsarsmillisecondopticallightcurvesdirectheatingmodelICARUSgamma-raybeamingRoche-lobefillingcompanionevolution

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

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.

What carries the argument

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.

What would settle it

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.

Watch

Extended reading notes

Core claim

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.

Load-bearing premise

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.

Editorial extensions

If this is right

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.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

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.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the 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.

What 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.

The 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."

Should 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.

Referee Report

2 major / 5 minor

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.

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 (2)

[Sec. 4.1, Sec. 4.3, Fig. 8] 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.
[Sec. 4, Table 3] 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.

minor comments (5)

[References, Sec. 2] 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.
[Sec. 2, Table 1] The text says J1124 was found in the 'GBT 350 GHz drift-scan survey'; this should presumably be 350 MHz, not 350 GHz.
[Sec. 3, Table 2] 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.
[Fig. 8] 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.
[Sec. 4.1] 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.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the core comparisons rely on external Fermi-LAT fluxes, DM distances, and an independent timing parallax, not on fitted quantities renamed as predictions.

full rationale

Walking the claimed derivation chain, the paper fits companion light curves with the ICARUS direct-heating model to obtain distance, fill factor, inclination, and heating power (Section 4), then compares the fitted heating power with observed Fermi-LAT GeV fluxes and with DM/timing distances (Sections 4.1-4.3). The Fermi-LAT fluxes are external observables and are not constructed from the fitted parameters, so the heating-versus-GeV comparison is not circular. The outer-gap beaming factors are taken from Romani & Watters (2010), a self-citation, but the paper explicitly presents them as a 'concrete example' with known limitations, averages over the allowed angle ranges, and concludes that residual eta_gamma versus eta_H disagreements indicate imperfect beaming corrections, distance errors, or unmodeled heating. That is a falsifiable comparison, not a self-citation used to force the conclusion. The limitations emphasized by the skeptic view, such as J0952's distance-fill-factor-metallicity degeneracy and J1124's photometric distance being 2.5 times the DM estimate, are explicitly acknowledged fit-degeneracy concerns rather than inputs relabeled as predictions. The J0636 timing parallax provides an independent distance check, and radial velocities for J1301 and J1959 are external constraints. No equation in the paper reduces to its own input by construction, and the central quantitative claims are conditional on model assumptions that the authors themselves flag; that is a correctness or model-validation risk, not circularity. Score 0.

Assumptions & free parameters 7 free parameters · 5 assumptions · 0 invented entities

The central results depend on five model parameters fitted per target, a metallicity grid choice, and several external assumptions about heating geometry, model atmospheres, extinction, and neutron star mass. No new particles, forces, or entities are introduced.

free parameters (7)

LX (bolometric heating luminosity) = 0.97 to 30.1 x 10^33 erg/s for the nine targets (Table 3)
Central fitted quantity in the ICARUS direct-heating model; it determines the companion heating and is compared with gamma-ray luminosity.
d (source distance) = 1.05 to 3.97 kpc for the nine targets (Table 3)
Fitted distance from the light curve, sometimes inconsistent with DM distances; used in efficiency calculations.
fc (Roche-lobe fill factor) = 0.48 to 0.99 depending on target (Table 3)
Ratio of companion radius to L1 radius; fitted.
TN (night-side temperature) = 1140 to 3340 K (Table 3)
Uniform underlying temperature of the unheated side; fitted.
i (binary inclination) = 23.3 to 77 degrees (Table 3)
Orbital inclination fitted from the light curve shape; critical for mass estimates.
log[Z] (metallicity) = Chosen among -1, -0.5, 0, +0.3, +0.5 BTSettl models; best fit per target (Figure 7)
Metallicity of the companion atmosphere is not fitted continuously but selected from a grid.
AV (extinction) = Adopted from dust maps; for J2052, free fit stays near 0.3 mag
For most targets AV is adopted from Bayestar2019; for J2052 the model allows AV to vary.

assumptions (5)

domain assumption The companion is heated directly by a point-source bolometric luminosity propagating linearly from the pulsar (the ICARUS direct-heating model).
This is the physical model stated in Section 4 and used for all fits; the paper notes that wind and shock heating are ignored.
domain assumption The BTSettl model atmospheres of Allard et al. (2012) describe the surface flux of the companions across the full temperature range.
Used to convert temperature distributions to filter magnitudes (Section 4); errors in these atmospheres would propagate into all fitted parameters.
ad hoc to paper The neutron star mass is assumed to be 1.5 Msun for objects without radial-velocity measurements (J0023, J0251, J0636, J0952, J1124, J2052, J2241).
Stated in Section 4; companion masses scale as MNS^(2/3), so this assumption affects the mass and density conclusions.
domain assumption The outer-gap beaming model of Romani and Watters (2010), computed for young pulsars, provides valid flux correction factors for these millisecond pulsars.
Used in Section 4.1 and 4.3 to compare heating and gamma-ray efficiencies; the paper itself notes this model may not apply to BWs.
domain assumption The extinction values from the Bayestar2019 dust map and from Schlafly and Finkbeiner (2011) are correct.
Adopted in Section 4 for the observed colors; wrong extinction would change the fitted temperature and distance.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Multiband Optical Light Curves of Black-Widow Pulsars." pith.science (2026). https://pith.science/paper/ZSSAKPQJ

@misc{pith2026190800992,
  author       = {Pith},
  title        = {Pith review of: Multiband Optical Light Curves of Black-Widow Pulsars},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/ZSSAKPQJ}},
  note         = {Machine review of arXiv:1908.00992}
}

abstract

We collect new and archival optical observations of nine "black-widow" millisecond pulsar binaries. New measurements include direct imaging with the Keck, Gemini-S, MDM, and LCO 2~m telescopes. This is supplemented by synthesized colors from Keck long-slit spectra. Four black-widow optical companions are presented here for the first time. Together these data provide multicolor photometry covering a large fraction of the orbital phase. We fit these light curves with a direct (photon) heating model using a version of the ICARUS light-curve modeling code. The fits provide distance and fill-factor estimates, inclinations, and heating powers. We compare the heating powers with the observed GeV luminosities, noting that the ratio is sensitive to pulsar distance and to the gamma-ray beaming. We make a specific correction for "outer-gap" model beams, but even then some sources are substantially discrepant, suggesting imperfect beaming corrections and/or errors in the fit distance. The fits prefer large metal abundance for half of the targets, a reasonable result for these wind-stripped secondaries. The companion radii indicate substantial Roche-lobe filling, $f_c \approx 0.7-1$ except for PSR J0952$-$0607, which with $f_c< 0.5$ has a companion density $\rho \approx 10\,{\rm g\,cm^{-3}}$, suggesting unusual evolution. We note that the direct-heating fits imply large heating powers and rather small inclinations, and we speculate that unmodeled effects can introduce such bias.

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

Optical detection of the black widow binary PSR J2052+1219
astro-ph.SR 2019-09 conditional novelty 5.0 of 10

The optical companion of PSR J2052+1219 is detected, and light curve modeling yields a distance of 3.94±0.16 kpc, a 0.034 M_sun companion close to filling its Roche lobe.

Reference graph

Works this paper leans on

24 extracted references · 19 canonical work pages · cited by 1 Pith paper

[1]

2012, RSPTA, 370, 2765

Allard, F., Homeier, D., & Freytag, B. 2012, RSPTA, 370, 2765

2012
[2]

M., Oezel, F., Barr, E., Champion , D

Antionadis, J., Tauris, T. M., Oezel, F., Barr, E., Champion , D. J. et al. 2016, ArXiv:1605.01665

arXiv 2016
[3]

G., Pleunis, Z., Hessels, J

Bassa, C. G., Pleunis, Z., Hessels, J. W. T., et al. 2017, ApJ, 846, L20

2017
[4]

S., Ransom, S

Boyles, J., Lynch, R. S., Ransom, S. M., et al. 2013, ApJ, 763, 80

work page 2013
[5]

P ., Rappaport, S

Breton, R. P ., Rappaport, S. A., van Kerkwijk, M. H., & Carter , J. A. 2012, ApJ, 748, 115

work page 2012
[6]

T., Camilo, F., Kerr, M., et al

Cromartie, H. T., Camilo, F., Kerr, M., et al. 2016, ApJ, 819, 34

2016
[7]

Draghis, P ., & Romani, R. W. 2018, ApJ, 862, L6

work page 2018
[8]

M., Schlaﬂy, E

Green, G. M., Schlaﬂy, E. F., Finkbeiner, D. P ., et al. 2015, A pJ, 810, 25

work page 2015

Show all 24 references

[9]

K., & Kazanas, D

Kalapotharakos, C., Brambilla, G., Timokhin, A., Harding, A. K., & Kazanas, D. 2018, ApJ, 857, 44

2018
[10]

L., Stovall, K., van Kerkwijk, M

Kaplan, D. L., Stovall, K., van Kerkwijk, M. H., Fremling, C., & Istrate, A. G. 2018, ApJ, 864, 15

2018
[11]

J., Johnston, S., Ray, P

Kieth, M. J., Johnston, S., Ray, P . S., et al. 2011, MNRAS, 414 , 1292

2011
[12]

P ., & Thorstensen, J

Li, M., Halpern, J. P ., & Thorstensen, J. R. 2014, ApJ, 795, L1 15

2014
[13]

J., Bassa, C

Nieder, L., Clark, C. J., Bassa, C. G., et al. 2019, ArXiv1905 11352

2019
[14]

Perley, D. A. 2019, ArXiv190307629

2019
[15]

T., Callanan, P

Reynolds, M. T., Callanan, P . J., Fruchter, A. S., et al. 2007 , MNRAS, 379, 1117

2007
[16]

W., Graham, M

Romani, R. W., Graham, M. L., Filippenko, A. V ., et al. 2015, A pJ, 809, L10

2015
[17]

W., Graham, M

Romani, R. W., Graham, M. L., Filippenko, A. V ., & Zheng, W.K.2016, ApJ, 833, 138

2016
[18]

W., & Sanchez, N

Romani, R. W., & Sanchez, N. 2016, ApJ, 828, 7

2016
[19]

W., & Watters, K

Romani, R. W., & Watters, K. P . 2010, ApJ, 714, 810

2010
[20]

Sanchez, N., & Romani, R. W. 2017, ApJ, 845, 42 Schlaﬂy, E. F.., & Finkbeiner, D. P . 2011, ApJ, 737, 103

2017
[21]

A., Bruel, P ., Cognard, I., et al

Smith, D. A., Bruel, P ., Cognard, I., et al. 2019, ApJ, 871, 78

2019
[22]

M., Archibald, A

Stovall, K., Ransom, S. M., Archibald, A. M., et al. 2014, ApJ , 791, 67

2014
[23]

2019, ApJ, 872, 42

Strader, J., Swihart, S., Chomiuk, L., Bahramian, A., Britt , C., et al. 2019, ApJ, 872, 42

2019
[24]

2016, MNRAS, 457, 3384 van Kerkwijk, M

Tchekhovskoy, A., Philippov, A., & Spitkovsky, A. 2016, MNRAS, 457, 3384 van Kerkwijk, M. H., Breton, R. P ., & Kulkarni, S. R. 2011, ApJ , 728, 95 Y ao, J. M., Manchester, R. N., & Wang, N. 2017, ApJ, 835, 29

2016

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.

[1] [1]

2012, RSPTA, 370, 2765

Allard, F., Homeier, D., & Freytag, B. 2012, RSPTA, 370, 2765

2012

[2] [2]

M., Oezel, F., Barr, E., Champion , D

Antionadis, J., Tauris, T. M., Oezel, F., Barr, E., Champion , D. J. et al. 2016, ArXiv:1605.01665

arXiv 2016

[3] [3]

G., Pleunis, Z., Hessels, J

Bassa, C. G., Pleunis, Z., Hessels, J. W. T., et al. 2017, ApJ, 846, L20

2017

[4] [4]

S., Ransom, S

Boyles, J., Lynch, R. S., Ransom, S. M., et al. 2013, ApJ, 763, 80

work page 2013

[5] [5]

P ., Rappaport, S

Breton, R. P ., Rappaport, S. A., van Kerkwijk, M. H., & Carter , J. A. 2012, ApJ, 748, 115

work page 2012

[6] [6]

T., Camilo, F., Kerr, M., et al

Cromartie, H. T., Camilo, F., Kerr, M., et al. 2016, ApJ, 819, 34

2016

[7] [7]

Draghis, P ., & Romani, R. W. 2018, ApJ, 862, L6

work page 2018

[8] [8]

M., Schlaﬂy, E

Green, G. M., Schlaﬂy, E. F., Finkbeiner, D. P ., et al. 2015, A pJ, 810, 25

work page 2015