---
title: 'LS I +61 303: Gamma-Ray Binary Overview'
url: https://www.emergentmind.com/topics/ls-i-61-303
type: topic
---

# LS I +61 303: Gamma-Ray Binary Overview

LS I +61 303 is a gamma-ray binary, and also a high-mass X-ray binary, composed of a B0 Ve star and a compact object in a $\sim 26.5$ day orbit, with emission that is variable and periodic across the electromagnetic spectrum from radio to very-high-energy gamma rays. A super-orbital modulation of $\sim 4.5$–$4.6$ years is also observed. For decades the nature of the compact object was debated between neutron-star/pulsar and black-hole or microquasar interpretations, but transient radio pulsations with $P=269.15508 \pm 0.00016$ ms provided the first evidence for pulsations from the source at any frequency and strongly argue for the existence of a rotating neutron star in LS I +61 303 [2203.09423].

## 1. System parameters and orbital descriptions

LS I +61 303 has been described as a system containing a massive B0 Ve star and a compact object at a distance of $\sim 2.0$ kpc, while another study gives $2.65 \pm 0.09$ kpc. The Be star mass is given as $10$–$15\,M_{\odot}$ in very-high-energy studies, and the source is identified as one of the brightest *Fermi* sources with orbitally modulated emission across the electromagnetic spectrum [2108.09235, 2203.09423, 1111.0042].

Published orbital solutions are not uniform. An earlier convention anchored the orbital phase $\phi({\rm TG})$ at zero at HJD $2{,}443{,}366.775$, with periastron at $\phi = 0.275$ and an eccentric orbit used in optical and TeV work [1111.0042]. By contrast, newer measurements cited in later very-high-energy and X-ray modeling papers suggest a relatively circular orbit with eccentricity $e < 0.15$ and periastron passage at phase $0.6$, while a precessing stellar disk model finds a preferred small eccentricity $e \simeq 0.06$ and periastron phase around $\Phi_{\rm p} \sim 0.6$ [2108.09235, 2409.04818]. The orbital motion on the sky is not well determined, and this leads to ambiguity in interpretation of the X-ray electric vector polarization angle measurement [2605.15972].

This multiplicity of orbital descriptions is central to the literature because phase assignments determine how radio morphology, H$\alpha$ spectroscopy, GeV/TeV light curves, and X-ray polarization are interpreted.

## 2. Orbital, precessional, and super-orbital periodicities

The fundamental orbital period is consistently near $26.5$ days. Long-term radio monitoring by the Owens Valley Radio Observatory found that the most powerful features in the periodogram are two peaks at $P_1 = 26.49 \pm 0.05$ d and $P_2 = 26.93 \pm 0.05$ d, with a long-term period $P_{\rm long} = 1698 \pm 196$ d. Dividing the OVRO data into three segments of equal length showed that the two periods, $P_1$ and $P_2$, are present in the periodogram of each of the consecutive long-term cycles, while the autocorrelation function of the full radio archive showed a regular pattern, proving the stability of the super-orbital modulation [2402.07719].

Within the radio timing framework, $P_1$ is interpreted as the orbital period, $P_2$ as a precession period, and $P_{\rm long}$ as a beat period satisfying
\[
\frac{1}{P_{\rm long}} = \left| \frac{1}{P_1} - \frac{1}{P_2} \right|.
\]
The same work reports a possible systematic modulation of the radio flux density with a timescale of approximately $40$ years that had so far remained unnoticed [2402.07719].

Very long baseline astrometry gives an independent measurement of the precession period. Multi-epoch VLBA phase-reference observations found that the observed elliptical trajectory of emission at $8.4$ GHz repeats after a $9$-year gap, and alignment of the 2006 and 2015 ellipses yielded a precession period of $26.926 \pm 0.005$ d. That study modeled the radio core as a precessing synchrotron-emitting jet tracing an elliptical trajectory on the sky [1711.07598]. Simultaneous long-term OVRO and *Fermi*-LAT monitoring further found an offset of the long-term modulation between radio and gamma-ray data, interpreted as evidence that the radio and GeV emission arise from different locations along a precessing jet [1804.07692].

A different super-orbital interpretation has been developed under the pulsar-wind scenario. A precessing stellar disk model attributes the $\sim 4.6$ year modulation to precession of the Be disk, reproduces asymmetric two-peak orbital X-ray light curves and sinusoidal-like long-term modulations, and argues that the neutron star is likely orbiting with a small eccentricity and periastron phase around $\Phi_{\rm p} \sim 0.6$ [2409.04818]. The literature therefore contains two distinct long-term geometric frameworks: precession of a relativistic jet and precession of the Be disk.

## 3. Multiwavelength phenomenology

LS I +61 303 shows emission from radio to TeV energies, with marked orbital and super-orbital structure. In broad terms, radio, X-ray, and TeV peaks occur at orbital phases $0.5$–$0.8$, especially near $0.6$–$0.7$, whereas the GeV emission peaks around orbital phase $\sim 0.25$ and is anti-correlated with X-ray and TeV emission. The GeV spectrum is described by a power law with an exponential cutoff at $3.9$ GeV [2108.09235].

Optical spectroscopy ties the high-energy phenomenology to the circumstellar disk. High-cadence H$\alpha$ spectra obtained over $35$ consecutive nights across a full orbit showed a dramatic burst of H$\alpha$ emission near $\phi({\rm TG}) \sim 0.6$, observed as a pronounced redshifted shoulder in the line profile. Difference spectra revealed an evolving S-shaped residual for $0.9 \leq \phi({\rm TG}) \leq 0.6$, indicative of a spiral density wave, while the equivalent width $W_{\rm H\alpha}$ decreases slightly before periastron and rises dramatically near $\phi \sim 0.6$. A correlated burst in radio, X-ray, and GeV emission is observed at the same orbital phase [1111.0042].

The H$\alpha$ event has been interpreted as the signature of interaction between a tidal mass stream drawn from the Be circumstellar disk and a relativistic pulsar wind, producing a compact pulsar wind nebula and a Balmer-dominated shock. In that picture, the redshifted shoulder traces high-velocity, receding post-shock material, while the equivalent-width changes trace transient reduction and later replenishment of the Be disk [1111.0042].

Multiwavelength campaigns have also examined cross-band correlations. During the 2014–2015 VERITAS, *Swift*-XRT, and *Fermi*-LAT campaign, no observable correlation was found between TeV and GeV fluxes, with a Pearson correlation coefficient of $0.07^{+0.36}_{-0.39}$, whereas a positive correlation was found between TeV and X-ray emission, with Pearson coefficient $0.80^{+0.14}_{-0.38}$ for contemporaneous observations [1508.06674]. On shorter timescales, simultaneous radio and X-ray monitoring found, for the first time, that the radio and X-ray emission are correlated up to $81$ per cent, with their few per cent variability correlated up to $40$ per cent; a radio quasi-periodic oscillation at $P = 1.75 \pm 0.10$ h was also detected [2011.13530].

These timing and correlation patterns are one of the main reasons the source has remained central to discussions of multi-zone emission, because they imply that different bands can track both common particle populations and physically distinct emission sites.

## 4. Compact object, pulsations, and the magnetar question

Before the pulsar detection, deep radio searches produced only upper limits. A phased-array search with the Giant Metrewave Radio Telescope at $1280$ MHz, centered at orbital phase $0.54$, found no pulses and set a minimum detectable mean flux density of $\sim 0.38$ mJy at the $8\sigma$ level for a putative pulsar with period $P > 2$ ms and duty cycle $D = 10\%$. That study concluded that the most favorable future searches should target frequencies $\sim 0.5$–$5$ GHz and orbital phases $0.6$–$0.7$, while noting that detection may still be impossible if the pulsar is not beamed at Earth or if binary absorption is strong [1206.0543].

The decisive observational change came from FAST. Four L-band observations totaling $\sim 10.2$ h yielded a significant periodic pulse train in one observation near orbital phase $\sim 0.59$, with period $P = 269.15508 \pm 0.00016$ ms, dispersion measure ${\rm DM} = 240.1~{\rm pc~cm^{-3}}$, pulse width $W = 33.30 \pm 0.96$ ms, and significance $>20\sigma$. Forty-two single, bright pulses were detected in the same dataset, while the other three observations showed no pulsations. The chance-coincidence probability for an unrelated pulsar within the FAST L-band beam was described as negligibly small, and the result was taken to settle the neutron star versus black hole or accreting microquasar debate in favor of a rotationally powered neutron star [2203.09423].

A distinct but related line of argument concerns magnetar-like behavior. A short soft-gamma burst detected by *Swift*-BAT from the direction of LS I +61 303 had $T_{90} \approx 0.24$ s and thermal spectrum with $kT \approx 7.5$ keV, resembling events generally labeled as magnetar-like. On that basis, a magnetar-containing flip-flop model was proposed, in which the system alternates from a rotationally powered regime at apastron to a propeller regime at periastron along the eccentric orbit [1109.5008]. Later analysis of transient radio pulses, spin evolution, accretion dynamics, age limits, and gamma-ray activity found that no conclusive evidence requires a “canonical” magnetar with $B_p > 10^{14}$ G. Instead, a strong propeller torque could account for much of the maximum allowed spin-down and weaken the inferred magnetic field by more than an order of magnitude, leaving the data compatible with $B_p \sim 10^{13}$–$10^{14}$ G [2210.09471].

Accordingly, current discussions of the compact object no longer center on neutron star versus black hole, but rather on the neutron star’s magnetic field strength, torque balance, and how rotationally powered and propeller-like states contribute to the orbital phenomenology.

## 5. High-energy and ultra-high-energy emission

LS I +61 303 is a persistent very-high-energy source with strong orbit-to-orbit variability. Earlier VERITAS observations from 2008 to 2010 accumulated $64.5$ h above $300$ GeV and found no strong evidence for TeV emission during the previously reliably active apastron phases, with upper limits less than $5\%$ of the Crab Nebula flux in that region. Instead, in late 2010 significant emission was detected close to superior conjunction, much closer to periastron passage, at $5.6\sigma$ post-trials significance [1105.0449]. This established that the TeV orbital light curve is not immutable.

The most dramatic TeV activity reported so far at that stage came in late 2014. VERITAS observed exceptionally bright flares around apastron with peak fluxes above $25\%$ and, in another summary, above $30\%$ of the Crab Nebula flux above $300$ GeV. The flares occurred near orbital phase $\phi \approx 0.60$, had rise and decay times shorter than one day, and produced photons with reconstructed energies up to $\sim 10$–$13$ TeV. Their spectra were well described by power laws, and no significant spectral variability was observed between the average and flare states [1601.01812, 1508.06800].

Longer-term VERITAS analysis using more than $150$ h, and specifically $\sim 164$ h of good livetime from October 2009 to January 2021, found a $33\sigma$ full-orbit detection above $260$ GeV. The nightly orbital-phase-binned light curve shows a clear outburst peaking at phase $\sim 0.65$, and the whole-dataset spectrum is a simple power law with normalization $N_0 = (1.35 \pm 0.05)\times 10^{-12}\ {\rm TeV}^{-1}{\rm cm}^{-2}{\rm s}^{-1}$ at $E_0 = 1$ TeV and spectral index $\Gamma = -2.55 \pm 0.05$. No strong spectral variation was found between orbital bins, although there is a possible hint of a cutoff above $10$ TeV in quiet states [2108.09235].

The source has now been reported at still higher energies. LHAASO detected LS I +61 303 with significances of $9.2\sigma$ in WCDA ($1.4$–$30.5$ TeV) and $6.2\sigma$ in KM2A ($25$–$267$ TeV). In KM2A alone, $16$ photon-like events above $100$ TeV were identified against an estimated $5.1$ background events, corresponding to a $3.8\sigma$ detection, and orbital modulation was observed at $4.0\sigma$ confidence between $25$ and $100$ TeV. The measured spectrum from $1.4$ to $267$ TeV is fit by a single power law with $\Gamma = 3.00 \pm 0.05$, and the results were interpreted in a composite scenario in which leptonic and hadronic processes jointly contribute [2510.23345].

The high-energy picture is therefore not limited to a single recurrent TeV outburst. It includes phase-dependent suppression, exceptionally bright day-scale flares, full-orbit TeV baseline emission, and now ultra-high-energy photons beyond $100$ TeV.

## 6. Radio morphology, astrometry, polarization, and geometric diagnostics

High-resolution radio imaging has long been central to the system’s interpretation. Reviews before the pulsation detection emphasized that the morphology of high-resolution radio images changes with orbital phase and had been interpreted either as a cometary tail in a pulsar-wind interaction picture or as a precessing jet in a microquasar picture [1710.00815]. The VLBI astrometric work strengthened the jet-precession interpretation by showing a stable elliptical trajectory of the $8.4$ GHz radio emission over a $9$-year baseline, with changes between one-sided and two-sided morphology and geometric parameters including precession cone angle $\psi = 21^\circ$, jet opening angle $\xi = 6^\circ$, orbital inclination $\zeta = 25^\circ$, and major-axis orientation $\Theta = -43.6^\circ$ [1711.07598].

That same astrometric study measured the absolute proper motion as $-0.150 \pm 0.006$ mas yr$^{-1}$ eastward and $-0.264 \pm 0.006$ mas yr$^{-1}$ northward. After removing Galactic rotation, it found a peculiar space velocity of $16~{\rm km~s^{-1}}$ and interpreted the small, $< 20~{\rm km~s^{-1}}$, non-circular motion as indicating a very low kick velocity when the black hole was formed [1711.07598]. In light of the later radio-pulsation detection, that specific compact-object identification is superseded, but the low peculiar motion remains an important kinematic constraint.

X-ray polarimetry now adds a direct probe of the magnetic field geometry at the particle-acceleration site. IXPE observations over orbital phases $0.74$ to $1.05$ detected polarization at $4.2\sigma$ significance with average polarization degree $13.1\% \pm 3.0\%$ in the $2$–$8$ keV band after background subtraction, and measured ${\rm EVPA} = -23.8^{\circ} \pm 6.6^{\circ}$. No statistically significant difference in polarization degree or EVPA was found between the two sampled orbits or between high and low flux intervals [2605.15972].

Interpretation of the EVPA depends on the adopted orbital elements. Using radial-velocity elements combined with radio imaging suggests an offset of approximately $26^\circ$–$28^\circ$ between the EVPA and the projected compact object–massive star axis; optical polarimetry gives either moderate alignment or an approximately $90^\circ$ offset depending on the adopted solution; and keV/TeV light-curve orbital elements give good alignment, with a difference $\lesssim 4^\circ$ [2605.15972]. This makes polarization not merely a magnetic-field diagnostic, but also a discriminator among competing orbital solutions.

Taken together, radio astrometry, long-term timing, optical spectroscopy, and X-ray polarimetry show that LS I +61 303 is not characterized by a single uncontested geometric model. Rather, it is a system in which precession, disk structure, shock geometry, and orbital uncertainty all remain active components of the interpretation, even after the compact object itself has been identified as a neutron star.

Source: https://www.emergentmind.com/topics/ls-i-61-303