---
title: 'PSR J2339-0533: Redback Millisecond Pulsar Binary'
url: https://www.emergentmind.com/topics/psr-j2339-0533
type: topic
---

# PSR J2339-0533: Redback Millisecond Pulsar Binary

PSR J2339−0533 is a spider millisecond-pulsar binary that is now generally modeled as a redback: a pulsar with spin period \(P \approx 2.9\) ms in a \(P_B \approx 4.6\) hr orbit with a strongly irradiated companion of mass \(\approx 0.3\,M_\odot\). The system shows nearly continuous radio eclipses, large orbital-period instabilities, hard orbitally modulated X-rays, and low-energy pulsed \(\gamma\)-ray orbital modulation. These properties have made it a central case for studies of intrabinary shocks, anisotropic pulsar winds, companion-surface heating, and time-variable orbital dynamics [1905.12591] [2005.07888] [1504.07466] [2009.04105].

## 1. Identification and classification history

PSR J2339−0533 entered the literature as the optical and X-ray counterpart of the bright, steady Fermi-LAT source 0FGL J2339.8−0530, later cataloged as 2FGL J2339.6−0532. Early optical photometry showed a large orbital modulation exceeding \(2.5\) mag, and spectroscopy revealed strong phase-dependent heating. On that basis, the source was initially interpreted as a likely “black widow” MSP binary: the 2011 combined light-curve and radial-velocity fit gave \(M_c = 0.075 \pm 0.007\,M_\odot\), \(M_p = 1.40 \pm 0.04\,M_\odot\), and \(i = 57.4^\circ \pm 0.5^\circ\) [1111.3074].

Later work revised this picture substantially. Gamma-ray timing, optical spectroscopy, and X-ray intrabinary-shock modeling instead treat J2339−0533 as a redback with a relatively massive companion of \(\approx 0.3\,M_\odot\) [1504.07466] [1905.12591] [2009.04105]. The redback interpretation is supported by the hard X-ray spectrum, the characteristic double-peaked X-ray orbital light curve bracketing \(\phi_B \approx 0.75\), strong day-side optical heating, and the inference that the companion wind dominates the pulsar-wind interaction region [2005.07888].

The classification history is itself astrophysically informative. The early “black widow” solution emerged from an irradiation model with a very low companion mass, whereas later optical analyses introduced asymmetric heating and magnetic-pole hot spots, and X-ray analyses modeled a redback-type intrabinary shock. This suggests that the inferred binary parameters are sensitive to how irradiation, center-of-light radial velocities, and the companion’s surface-brightness asymmetry are treated [1111.3074] [2009.04105].

## 2. Spin, orbit, and binary architecture

The pulsar spin frequency is reported as \(\nu = 346.71337922047(2)\ {\rm s}^{-1}\), with spin-down rate \(\dot{\nu} = -1.6945(2)\times 10^{-15}\ {\rm s}^{-2}\), corresponding to \(P \approx 2.9\) ms [2005.07888]. A phase-coherent LAT timing solution over more than six years gave \(f = 346.71337922051(2)\) Hz, \(\dot f = -1.6952(8)\times 10^{-15}\ {\rm Hz\ s}^{-1}\), \(P = 2.8842267415473(2)\) ms, and \(\dot P = 1.4102(6)\times 10^{-20}\ {\rm s\ s}^{-1}\) [1504.07466].

The orbit is compact: \(P_{\rm orb} = 0.1930984018(3)\) d in the timing solution, or \(P_B = 0.19\) d \(\approx 4.6\) hr in later high-energy analyses [1504.07466] [2005.07888]. The projected semi-major axis is measured as \(x = 0.611656(4)\) lt-s in the timing work and \(A1 = 0.611668(3)\) lt-s in the 11-year LAT analysis [1504.07466] [2005.07888]. The timing model uses an ELL1-like parametrization together with high-order orbital-frequency derivatives to follow the strong time variability of the orbit [1504.07466].

The orbital eccentricity is small but measured: \(e = 2.1(1)\times 10^{-4}\), with \(\epsilon_1 = +1(1)\times 10^{-5}\) and \(\epsilon_2 = -21(1)\times 10^{-5}\) [1504.07466]. In contrast, the later optical hot-spot analysis adopts a circular orbit, \(e \approx 0\), for light-curve and radial-velocity modeling [2009.04105]. This difference reflects methodology rather than a contradiction in the existence of a very low eccentricity.

Representative masses and dimensions also depend on the analysis. The 2020 hot-spot optical fit gives \(M_{\rm NS} = 1.47 \pm 0.09\,M_\odot\), \(M_C = 0.30 \pm 0.02\,M_\odot\), \(q \equiv M_C/M_{\rm NS} = 0.204 \pm 0.019\), \(a = 1.71 \pm 0.03\,R_\odot\), and companion radius \(R_C \approx 0.42 \pm 0.02\,R_\odot\) [2009.04105]. By contrast, the anisotropic IBS X-ray fit adopts \(M_{\rm NS} \approx 2.16\,M_\odot\) and \(M_c \approx 0.3\,M_\odot\) for representative geometry, giving \(a \approx 1.3\times 10^{11}\) cm \(\approx 1.9\,R_\odot\) [1905.12591].

Distance estimates are similarly model-dependent. Reported values include a radio dispersion-measure distance of \(\sim 450\) pc, \(d = 1.25\) kpc in the J2339_X X-ray fit, \(d = 1.87 \pm 0.01\) kpc in the GROND hot-spot solution, and \(d = 1.1\) kpc in the broadband X-ray/\(\gamma\)-ray IBS modeling [1504.07466] [1905.12591] [2009.04105] [2402.02674].

## 3. Companion heating, surface asymmetry, and optical constraints

The companion is strongly irradiated and shows large-amplitude, color-dependent orbital modulation. Early spectroscopy found a phase-dependent effective temperature ranging from \(\approx 6900\) K near superior conjunction to \(\approx 2900\) K near inferior conjunction, with the spectrum evolving from approximately F3 at maximum to approximately M5 at minimum [1111.3074]. Later multi-epoch photometry and spectroscopy describe day–night temperature contrasts from \(\approx 7000\) K on the irradiated side to \(\approx 3500\) K near minimum light [2009.04105].

A key development was the move from symmetric heating prescriptions to explicitly asymmetric ones. Using GROND, Keck, SOAR, WIYN, OISTER, and HET data aligned with the updated Fermi-LAT ephemeris, three heating prescriptions were compared: direct heating (DH), wind heating (WH), and hot-spot heating (HS). DH gave \(\chi^2/{\rm DoF} = 1388/553\), WH improved this to \(877/552\), and HS yielded the best fit at \(671/550\) [2009.04105]. In the HS model, a Gaussian excess on the companion surface represents a magnetic-cap hot spot produced by pulsar-wind or IBS particles ducted by the companion field.

The preferred GROND HS solution gives \(i = 69.3^\circ \pm 2.3^\circ\), \(K_C = 347.0 \pm 3.7\ {\rm km\ s}^{-1}\), systemic velocity \(\gamma = 17.7 \pm 2.1\ {\rm km\ s}^{-1}\), Roche-lobe fill factor \(f_c = 0.97 \pm 0.02\), base temperature \(T_N = 3307^{+64}_{-90}\) K, and irradiation parameter \(L_P = (1.48 \pm 0.03)\times 10^{34}\ {\rm erg\ s}^{-1}\) [2009.04105]. The maxima lag superior conjunction by \(\Delta\phi \approx -0.03\), and the leading day side is systematically brighter and bluer than the trailing side [2009.04105].

The hot spots are large, lie in the southern hemisphere, and move with epoch. Reported best-fit spot parameters include \(\theta_{\rm hs} = 65.2 \pm 2.4^\circ\), \(\phi_{\rm hs} = -79.3 \pm 2.2^\circ\) for WIYN+OISTER; \(\theta_{\rm hs} = 85.0 \pm 8.0^\circ\), \(\phi_{\rm hs} = -80.1 \pm 5.7^\circ\) for SOAR; \(\theta_{\rm hs} = 70.1 \pm 1.1^\circ\), \(\phi_{\rm hs} = -53.3^{+5.1}_{-4.6}{}^\circ\) for GROND; and \(\theta_{\rm hs} = 124.5 \pm 17.2^\circ\), \(\phi_{\rm hs} = -59.0^{+26.0}_{-13.8}{}^\circ\) for Keck [2009.04105]. A two-spot dipolar test for the GROND epoch gave \(A_1/A_2 \approx 10\), consistent with more efficient IBS particle capture on field lines directed toward the pulsar nose [2009.04105].

These optical results matter beyond companion phenomenology. The paper explicitly links the hot spots to magnetic poles and to ducting of pulsar-wind or IBS particles, thereby connecting companion heating geometry to intrabinary-shock asymmetry [2009.04105]. A plausible implication is that some of the phase shifts seen in X-rays and optical bands are manifestations of the same magnetically mediated redistribution process.

## 4. Intrabinary shock structure and X-ray emission

The X-ray phenomenology is interpreted primarily in terms of synchrotron emission from an intrabinary shock formed where the pulsar wind collides with a strong companion outflow. In the anisotropic IBS model, the pulsar-wind energy or momentum flux is taken to be equatorially concentrated,
\[
F_p(\theta_\ast)\propto \sin^{2n}\theta_\ast,
\]
with \(\theta_\ast\) measured from the pulsar spin axis. For PSR J2339−0533 the modeling assumes \(n=1\), so \(F_p\propto \sin^2\theta_\ast\) [1905.12591]. The thin-shell contact discontinuity obeys a generalized Cantó-type solution,
\[
r(\theta,\phi)=D\,\sin\theta_1\,\csc(\theta+\theta_1),
\]
with momentum-flux ratio
\[
\beta = \dot{M}_c v_c c / \dot{E}.
\]
In the redback regime, \(\beta \gtrsim 1\), the shock wraps around the pulsar, and increasing \(n\) flattens the shock near the poles, producing an hourglass cross-section [1905.12591].

For J2339−0533, the X-ray IBS fit gives \(i = 54.05^\circ \pm 6.38^\circ\), \(\beta = 5.19 \pm 0.23\), \(d = 1.25\) kpc, and \(\dot E_{\rm IBS} = (4.70 \pm 0.30)\times 10^{34}\ {\rm erg\ s}^{-1}\) [1905.12591]. Because \(\beta>1\), the companion wind dominates and wraps the shock around the pulsar. The paper notes that the isotropic stand-off proxy would be \(R_0 \sim 0.3\,a \approx 4\times 10^{10}\) cm, but emphasizes that the full anisotropic geometry is solved numerically rather than by a single closed-form \(R_0\) [1905.12591].

The shocked flow is mildly relativistic and Doppler boosted. The flow Lorentz factor is parameterized as
\[
\Gamma(s)=\Gamma_0\big(1+k\,s/r_0\big),
\]
with redback values \(\Gamma_0 \approx 1.1\) and \(k \approx 0.2\), and the observed emission is transformed with \(\delta = [\Gamma(1-\beta_{\rm flow}\cos\psi)]^{-1}\) [1905.12591]. Synchrotron emission is computed using
\[
\nu_c = \frac{3eB}{4\pi m_e c}\gamma^2,\qquad
t_{\rm syn} = \frac{6\pi m_e c}{\sigma_T B^2\gamma},
\]
with particle injection
\[
Q(\gamma)=N_0\,\gamma^{-p}\,d\gamma.
\]
The very hard X-ray photon indices imply \(p\sim 1\) for uncooled synchrotron, which the paper interprets as likely evidence that reconnection dominates acceleration in a high-\(\sigma_w\) wind [1905.12591].

The observational basis is extensive. The anisotropic IBS study used archival Chandra ACIS-S (20 ks; ObsID 11791), Swift (49.4 ks cumulative), XMM-Newton (182 ks; ObsIDs 721130101 and 790800101), Suzaku (104 ks; ObsID 406007010), and NuSTAR (163 ks; ObsID 30202020002), with phase-resolved light curves and spectra fitted simultaneously in XSPEC [1905.12591]. The \(3\)–\(15\) keV orbital light curve shows a prominent double peak bracketing pulsar inferior conjunction at \(\phi_B = 0.75\), with separation \(\Delta\phi_B \approx 0.35\) and a strong bridge, especially at low energies [1905.12591]. The phase-resolved absorbed power-law fits are consistently hard: Off \((\phi_B = 0.05\text{–}0.45)\), \(\Gamma = 1.19 \pm 0.11\); P1 \((0.45\text{–}0.65)\), \(\Gamma = 0.96 \pm 0.11\); Bridge \((0.65\text{–}0.85)\), \(\Gamma = 0.80 \pm 0.13\); and P2 \((0.85\text{–}0.05)\), \(\Gamma = 0.77 \pm 0.12\) [1905.12591].

The hard-X-ray behavior extends the picture. Orbital modulation persists above \(40\) keV with chance probability \(p \approx 10^{-4}\), but the double peaks weaken or disappear above \(\approx 15\) keV and are replaced by a single or tight-double component centered at \(\phi_B = 0.75\) [1905.12591]. The paper argues that this evolution appears steeper than a simple cooling break and may instead probe the electron \(\gamma_{\max}\) in the peak zones [1905.12591].

A 2025 joint XMM-Newton+NuSTAR analysis sharpened the spectral constraints. Using XMM \(0.5\)–\(10\) keV and NuSTAR \(3\)–\(78\) keV data, with an off-phase baseline power law plus an additional IBS power law in the peak phase, the Bayesian Information Criterion prefers an unbroken IBS power law with \(N_H < 0.38\times 10^{21}\ {\rm cm}^{-2}\), \(\Gamma_{\rm IBS}=0.94^{+0.18}_{-0.20}\), \(\Gamma_0=1.71^{+0.42}_{-0.15}\), and \(\chi^2/{\rm DoF}=1.11\) [2503.24384]. Cooling-break and cutoff models are not BIC-preferred, but they place lower limits \(E_{b,\rm IBS}>7.1\) keV and \(E_{c,\rm IBS}>26\) keV, implying \(B_{\rm IBS}<63\) G and \(\gamma_{\max}>1.9\times 10^5\) [2503.24384]. The same paper interprets the source as likely viewed at high inclination, \(i \approx 70^\circ\), such that striped-wind field annihilation across the shock leaves a relatively modest post-shock field in the visible flow [2503.24384].

## 5. Orbital \(\gamma\)-ray modulation and competing high-energy scenarios

PSR J2339−0533 also shows orbital modulation in Fermi-LAT data, but with a phenomenology distinct from the X-ray intrabinary shock. Using approximately 11 years of Pass 8 LAT observations from 2008 Aug 04 to 2019 Jul 28, the modulation is detected only in the low-energy \(100\)–\(600\) MeV band. The signal is approximately sinusoidal, peaks near pulsar superior conjunction at \(\phi_B \approx 0.25\), and is confined to the on-pulse interval \(\phi = 0.20\text{–}0.78\) rather than the off-pulse interval \(\phi = 0.78\text{–}1.20\) [2005.07888]. The weighted H-test gives \(H \approx 40\) and chance probability \(p \approx 10^{-7}\) for the on-pulse \(100\)–\(600\) MeV light curve, whereas the off-pulse band gives \(H \approx 1\) [2005.07888].

This behavior sharply differs from standard IBS expectations. In redbacks, IBS beaming would normally produce an unpulsed GeV modulation strongest off-pulse and often aligned with the X-ray IBS geometry near inferior conjunction. J2339−0533 instead shows a pulsed modulation peaking near superior conjunction, opposite to the X-ray phase [2005.07888]. The phase-resolved spectral fits reinforce that distinction: the phase-averaged spectrum follows
\[
\frac{dN}{dE}=N_0(E/E_0)^{-\Gamma_1}\exp[-aE^{\Gamma_2}],
\]
with \(E_0=1.1\) GeV, \(\Gamma_2=0.67\), \(N_0=(9.5\pm0.8)\times10^{-12}\ {\rm ph\ MeV^{-1}\ cm^{-2}\ s^{-1}}\), \(\Gamma_1=1.06\pm0.07\), \(a=(7.8\pm0.6)\times10^{-3}\), and \(F_{0.1\text{–}300\,{\rm GeV}}=(2.2\pm0.1)\times10^{-8}\ {\rm ph\ cm^{-2}\ s^{-1}}\) [2005.07888]. At orbital maximum \(\Gamma_1 = 1.39 \pm 0.14\), whereas at minimum \(\Gamma_1 = 0.6 \pm 0.3\), so the variable component is inferred to be soft relative to a hard baseline [2005.07888].

The 2020 LAT study considers several mechanisms. Absorption of pulsed emission by the shocked pulsar wind is modeled with a Klein–Nishina cross section,
\[
\sigma_{\rm KN} = \frac{3\sigma_T}{4}\left[\cdots\right],\qquad x = \frac{E_\gamma}{\Gamma_w m_e c^2},
\]
and an attenuation factor \(e^{-\tau(E)}\) with \(\tau(E)=C(\phi_B)\sigma_{\rm KN}(E)\). However, the required optical depths are far too large: the Thomson depth of a cold pair wind is estimated as
\[
\tau_T \approx 7\times10^{-7}\frac{\dot E_{34}}{a_{11}\Gamma_w},
\]
with further suppression in the KN regime, so absorption is ruled out quantitatively [2005.07888]. The paper instead argues that inverse-Compton upscattering of companion photons by a highly coherent striped pulsar wind remains viable. In the Thomson limit,
\[
E_\gamma \approx \frac{4}{3}\gamma^2 E_{\rm seed},
\]
and the estimated luminosity can match the soft excess if \(\Gamma_w\delta r_\odot \approx 2\times10^4\) [2005.07888].

A later broadband study reached a different modeling conclusion. Using Chandra, XMM-Newton, NuSTAR, and Fermi-LAT data, it found that the X-ray emission is well explained by IBS synchrotron, but the modulated \(\gamma\)-ray component is difficult to reproduce with IBS components alone [2402.02674]. Two alternatives were examined. Scenario (1), inverse-Compton emission in the upstream unshocked wind, requires the wind to decelerate to \(\sim 1000\ {\rm km\ s^{-1}}\) before reaching the shock, which the paper judges inconsistent with strong-shock X-ray phenomenology [2402.02674]. Scenario (2), synchrotron radiation from shock-penetrating primaries in a companion magnetic field, can fit the modulation with \(B_c \approx 1.39\) kG, \(\gamma_p \approx 4\times10^7\), \(\zeta \approx 0.11\), and \(\eta_p \approx 0.82\), while the X-ray IBS parameters remain \(p \approx 1.37\), \(\gamma_{\max} \approx 6\times10^6\), \(B_s \approx 4.10\) G, \(\Gamma_D \approx 1.45\), and \(\beta \approx 0.20\) [2402.02674].

Taken together, these studies establish a point that is easy to misunderstand: the orbitally modulated \(\gamma\)-ray signal in J2339−0533 is not a straightforward GeV analog of the X-ray intrabinary shock. The 2020 analysis argues for a pulsed IC excess from a coherent striped wind, whereas the 2024 modeling prefers synchrotron emission in a kilogauss companion field after subtraction of a constant pulsed component [2005.07888] [2402.02674].

## 6. Orbital-period variability, companion magnetic activity, and system significance

One of the most distinctive dynamical features of PSR J2339−0533 is its extreme orbital-period variability. Precision gamma-ray timing over more than six years revealed alternating epochs of orbital-period decrease and increase, requiring a polynomial expansion of the orbital frequency with derivatives up to \(f_{\rm orb}^{(6)}\) [1504.07466]. The modulation is well described by
\[
\Delta f_{\rm orb}(t)=\dot f_{{\rm orb},\ast} t + A\sin\!\left[2\pi \nu_{\rm mod} t + \pi \dot\nu_{\rm mod} t^2 + \psi\right],
\]
with best-fit amplitude \(A = 8.48(2)\times 10^{-11}\ {\rm s}^{-1}\), modulation period \(P_{\rm mod}=4.15(1)\) yr, residual secular orbital-period derivative \(\dot P_{{\rm orb},\ast}=-1.57(2)\times10^{-11}\ {\rm s\ s}^{-1}\), and fractional amplitude \(\Delta P_{\rm orb}/P_{\rm orb} \approx 2.3\times10^{-7}\) [1504.07466]. Extremes of the instantaneous orbital-period derivative are \(\dot P_{\rm orb,min}\approx -5.8\times10^{-10}\ {\rm s\ s}^{-1}\) and \(\dot P_{\rm orb,max}\approx +2.7\times10^{-10}\ {\rm s\ s}^{-1}\) [1504.07466].

The favored explanation is gravitational quadrupole coupling in the companion, i.e. Applegate-type modulation driven by magnetic activity in the convective envelope. The key relation used is
\[
\frac{\Delta P_{\rm orb}}{P_{\rm orb}} = -\frac{9\,\Delta Q}{M_c a^2}.
\]
For the observed modulation, the required fractional outer-layer spin change is small, \(\Delta\Omega/\Omega \approx 10^{-5}\) for a shell mass \(M_s \approx 0.1 M_c\), and the implied luminosity variation is \(\Delta L \approx 4.1\times10^{28}\ {\rm erg\ s}^{-1}\), only \(\approx 0.03\%\) of the companion’s intrinsic luminosity \(L \approx 1.3\times10^{32}\ {\rm erg\ s}^{-1}\) [1504.07466].

Alternative mechanisms were quantified and disfavored. Gravitational-wave emission gives \(\dot P_{\rm orb}^{\rm GW}\approx -1.3\times10^{-13}\ {\rm s\ s}^{-1}\), Shklovskii and Galactic acceleration terms are of order \(10^{-14}\ {\rm s\ s}^{-1}\), and even extreme mass-loss assumptions give \(\dot P_{\rm orb}^{M}\approx -1.2\times10^{-11}\ {\rm s\ s}^{-1}\), all far too small and, crucially, monotonic rather than alternating [1504.07466]. The residual eccentricity \(e=2.1(1)\times10^{-4}\) is likewise consistent with ongoing convective or magnetic activity in the companion [1504.07466].

This dynamical picture connects naturally to the broader phenomenology. The system’s nearly continuous radio eclipses, hard X-ray intrabinary-shock emission, and strong companion heating all point to a pervasive companion wind and a magnetically active, near-Roche-lobe-filling star [1905.12591] [2009.04105]. Several observational tests have been proposed: deeper phase-resolved X-ray spectroscopy to detect spectral gradients across the IBS peaks, improved radio eclipse mapping to constrain the contact discontinuity, continued LAT monitoring to refine the low-energy pulsed modulation, and sub-GeV measurements to distinguish inverse-Compton from synchrotron scenarios [1905.12591] [2005.07888].

In that sense, PSR J2339−0533 is best understood not as a single-purpose laboratory but as an unusually coupled system: orbital dynamics track the companion’s internal magnetic activity, optical light curves encode asymmetric heating and magnetic poles, X-rays map a hard synchrotron-emitting intrabinary shock, and low-energy \(\gamma\) rays probe wind coherence or companion-field interactions on scales beyond the standard IBS picture [1504.07466] [1905.12591] [2005.07888] [2402.02674].

Source: https://www.emergentmind.com/topics/psr-j2339-0533