Pole Trajectory Analysis in Pulsar Emission

• Pole Trajectory Analysis is a framework that characterizes the evolution of the polarization vector by tracking its geodesic path on the Poincaré sphere during mode transitions.
• It employs models such as coherent, partially coherent, nonorthogonal incoherent, and orthogonal with elliptical emission to explain abrupt or gradual PA and EA changes.
• The analysis links statistical mode intensity fluctuations with observable polarization fractions, offering diagnostics to distinguish between varied physical emission mechanisms.

Pole trajectory analysis is a theoretical and observational framework focused on understanding the evolution and transitions of the polarization vector in radio pulsar emission, particularly as it traverses the Poincaré sphere during polarization mode transitions. This analysis involves characterizing the behavior of the position angle (PA) and ellipticity angle (EA) associated with the polarization state, especially during abrupt or gradual changes that occur when dominance shifts between orthogonal polarization modes or when vector rotations arise from changes in phase lag. Four primary models—coherent, partially coherent, incoherent with nonorthogonal modes, and incoherent orthogonal modes with elliptically polarized emission—have been examined to quantify and interpret the trajectories these transitions follow, with a central finding that such pole trajectory transitions universally trace geodesics (great circles) connecting the polarization vectors of the involved modes on the Poincaré sphere (McKinnon, 2024).

1. Characterization of Polarization Mode Transitions

During polarization mode transitions in pulsar radio emission, the PA and EA can experience pronounced, sometimes asymmetric excursions. An abrupt or smooth switch in the dominant polarization mode often yields a PA shift of up to $90^\circ$, frequently accompanied by a significant and possibly asymmetric EA excursion. The details of this behavior depend on the relative intensities of the two modes (often parameterized by $m$, the normalized intensity difference) and, critically, on deviations from perfect orthogonality in the linear ($\delta_\ell$) and circular ($\delta_V$) polarization components. When orthogonality is not exact, the PA and EA transitions are displaced from symmetry; for instance, the PA transition may be offset from $m = 0$, and the EA excursion may be asymmetric about the transition point, properties that are especially evident in the incoherent nonorthogonal polarization mode (NPM) model.

A central physical insight is that the trajectory of the polarization vector during such a transition follows the geodesic (great circle) connecting the two mode orientation vectors on the Poincaré sphere. This geometrically determined path acts as a "mean trajectory" for the polarization state as it evolves through the mode transition.

2. Models for Polarization Transitions

Four main models have been systematically analyzed:

• Coherent Polarization Modes (COH): Assumes modes are coherent with a defined phase offset and a mode-mixing angle; the polarization vector rotates about a fixed axis, yielding a gradual, smooth transition which traces a geodesic.
• Partially Coherent Modes (PCOH): Generalizes COH to allow fractional coherence, parameterized by a coherence factor, a mode strength ratio, and a phase lag. As the coherence factor tunes from 1 to 0, the model transitions from the COH case to the incoherent regime, with corresponding changes in PA and EA excursion symmetry.
• Incoherent Modes with Nonorthogonal Polarization Vectors (NPM): An extension of the McKinnon–Stinebring statistical model, this includes explicit nonorthogonality, capturing both PA and EA transitions that are asymmetric with respect to $m$. For instance, the PA transition is:

$$\psi(m) = \frac{1}{2} \arctan{\left[\frac{\tan{\delta_\ell}(1 - m)}{\tan^2{\delta_\ell} + m}\right]}$$

with asymmetry governed by deviations from orthogonality.

• Incoherent Orthogonal Modes with Elliptically Polarized Emission (EPC): Combines two incoherent, orthogonally polarized modes (separated by $90^\circ$ in PA) with an additional independent elliptical component. The additional component, parameterized by $\epsilon$, generally reduces the total PA transition to less than $90^\circ$ and yields symmetric EA excursions about $m = 0$.

Despite the diversity of physical assumptions, each model predicts that the polarization mode transition populates a geodesic path on the Poincaré sphere.

3. Geodesic Trajectories on the Poincaré Sphere

The invariant feature across all four models is that mode transitions trace a great circle (geodesic) between the orientation vectors of the involved polarization modes. Mathematically, for a given transition, the relationship of EA to PA along the geodesic is

$$\chi(\psi) = \frac{1}{2}\arctan{\left[\frac{\sin(2\psi)\tan 2\delta_V}{\sin(2\delta_\ell)}\right]}$$

and the angular length (great circle distance) of the transition is

$$\zeta = \pi - \arccos\left(\cos(2\delta_V)\cos(2\delta_\ell)\right)$$

The parameters $\delta_\ell$ and $\delta_V$ fully specify the orientation and extent of the geodesic. In the case of exact orthogonality ($\delta_\ell = \delta_V = 0$), the transition extends over a full $90^\circ$ PA change with a symmetric EA excursion.

4. Effect of Mode Intensity Fluctuations

The observed polarization fraction and the character of EA excursions are strongly influenced by the statistical nature of mode intensity fluctuations. When the intensities are constant and deterministic, the mean polarization fraction is linear in the normalized intensity difference, $p(m) = |m|$. Introducing statistical fluctuations, such as exponential or Gaussian distributions for the mode intensities, modifies this relationship. For example:

• Exponential fluctuations: $p(m) = (1 + m^2)/2$, yielding a minimum polarization fraction $p_m = 0.5$ at $m = 0$.
• Gaussian fluctuations: $$p(m) = \beta\sqrt{{2}/{\pi} \exp(-{m^2}/2\beta^2)+ \operatorname{erf}(m/\beta\sqrt{2}) m}$$, with minimum $p_m$ dependent on the modulation index $\beta$.

Larger fluctuations tend to "fill in" the polarization, increasing the mean polarization fraction and suppressing large EA excursions. When mode intensities are quasi-stable, EA excursions are maximized; when intensities are highly stochastic, excursions are minimized and the mean polarization fraction is elevated.

5. Interpretation and Discrimination among Models

Although all models yield similar geodesic trajectories for the mean polarization vector, subtler features—including the symmetry or asymmetry of the EA excursion, the total PA change, and the polarization fraction as a function of $m$—permit discrimination between underlying physical mechanisms. For example:

• Asymmetric EA excursions, observed in certain pulsars such as PSR B0329+54, can only be reproduced by the nonorthogonal (NPM) or suitably enhanced EPC models.
• Both coherent vector rotation (resulting from a change in phase lag, as in generalized Faraday rotation scenarios) and genuine transitions in mode dominance yield geodesic trajectories, but vector rotations typically follow small circles, and produce distinctive, often smaller, EA changes.

These distinctions can, in principle, be exploited in high-sensitivity polarization measurements to probe physical models of pulsar emission and the effect of propagation through the magnetospheric plasma.

6. Observational Implications

The realization that all models predict pole trajectory transitions along geodesic paths places clear, model-independent constraints on the expected evolution of the PA and EA during mode changes. However, the specifics of the transition—including the total excursion, its symmetry, and the observed polarization fraction—are sensitive to details such as the degree of orthogonality and the statistical character of the mode intensity fluctuations.

Observations of abrupt PA transitions with large, possibly asymmetric EA excursions (near the poles of the Poincaré sphere) are interpreted as mode transitions or vector rotations, with the precise model assignment depending on secondary features. Accordingly, subtleties such as the asymmetry of the EA excursion or deviation from a $90^\circ$ PA swing may serve as diagnostics for mode nonorthogonality, propagation-induced mixing, or underlying emission anisotropy. Conversely, similarity in the mean trajectory across models implies that, in the absence of such details, interpretation of observed transitions within any single model is not uniquely determined.

7. Summary Table: Distinguishing Features of Polarization Models

Model	Geodesic Trajectory	EA Excursion Symmetry	Max PA Change
Coherent (COH)	Yes	Symmetric	90°
Partially Coherent (PCOH)	Yes	Symmetric/intermediate	Up to 90°
Nonorthogonal Incoherent	Yes	Asymmetric	90°
Orthogonal + Elliptical (EPC)	Yes	Symmetric	<90° (if ε≠0)

These results collectively provide a comprehensive quantitative and geometric framework for interpreting polarization mode transitions and their associated pole trajectories in pulsar radio emission, emphasizing both the universal geodesic structure on the Poincaré sphere and the importance of statistical and model-specific details for physical understanding (McKinnon, 2024).