An Extreme Scattering Event Toward PSR J2313+4253

Technical Summary: An Extreme Scattering Event Toward PSR J2313+4253

Problem and Context
Extreme Scattering Events (ESEs) are rare phenomena characterized by frequency-dependent variations in the flux density of radio sources, theorized to arise from small-scale density variations in the ionized Interstellar Medium (ISM). While initially observed in quasars, ESEs have been detected in pulsars, manifesting as correlated changes in flux, Dispersion Measure (DM), scintillation bandwidth, and timescale. The physical origin of these events remains debated, with proposed models ranging from Gaussian plasma lenses to catastrophe theory involving $A_3$ cusp singularities. A critical challenge in studying ESEs is the lack of high-cadence monitoring, which limits the ability to track the evolution of these events and resolve fine-scale structural details of the scattering screens.

Methodology
This study utilizes high-cadence observations of the pulsar PSR J2313+4253 (B2310+42) obtained with the Green Bank Observatory 20m telescope as part of the Pulsar Science Collaboratory (PSC) Scintillation Survey. The dataset spans 777 days (MJD 60006–60783) with a near-biweekly cadence, increasing to near-weekly when a potential ESE was detected. Observations covered a frequency range of 1300–1800 MHz, analyzed within a 1350–1525 MHz band to minimize Radio Frequency Interference (RFI).

The analysis pipeline involved:

1. Dynamic Spectrum Generation: Using pypulse to extract dynamic spectra ($S(\nu, t)$) by subtracting off-pulse power from on-pulse power.
2. Scintillation Parameter Extraction: Calculating scintillation bandwidth ($\Delta\nu_d$) and timescale ($\Delta t_d$) by fitting 2D autocorrelation functions (ACF) with Lorentzian profiles, consistent with the ISM's pulse broadening function.
3. Secondary Spectrum Analysis: Generating 2D Fourier transforms of dynamic spectra to identify scintillation arcs and calculate arc curvature ($\eta$) using the Hough transform technique (scintools).
4. Statistical Validation: Employing sliding window t-tests to identify statistically significant deviations in scintillation parameters and refractive timescales during the event period.
5. Distance and Size Estimation: Deriving scattering screen distances ($D_s$) using both scintillation parameters (assuming a Kolmogorov spectrum and transverse velocity dominance) and arc curvature. Transverse sizes were estimated using the event duration and effective velocity.

Key Results
The authors present evidence of an ESE toward PSR J2313+4253 with the following characteristics:

• Temporal Evolution: The event lasted approximately 220 days, beginning around MJD 60306. It is characterized by a sharp drop in scintillation bandwidth, followed by a pair of spikes and a subsequent decline, accompanied by a simultaneous increase in scintillation timescale.
• Statistical Significance: A t-test on the scintillation bandwidth during the event yielded a p-value of 0.012, indicating the values are drawn from a different distribution than pre- and post-event epochs. The refractive timescale showed a statistically significant dip (p-value 0.02) at MJD 60442 and 60463.
• Scattering Screen Distance: The estimated distance to the dominant scattering screen shifted during the event. Pre- and post-event measurements yielded a weighted average of $1.00 \pm 0.04$ kpc. During the ESE, the distance was more tightly constrained at $1.04 \pm 0.01$ kpc. This shift was statistically significant (p-value 0.0006).
• Double-Lensing Signature: Secondary spectra processed during the event revealed a "detached feature" at a delay of 0.6 $\mu$s, distinct from the main scintillation arc (which terminated at 0.2 $\mu$s). This feature, observed only during the event, implies a double-lensing geometry where the ESE structure acts as a secondary lens in front of the primary scattering screen.
• Physical Dimensions: Based on the event duration and the pulsar's transverse velocity, the transverse size of the scattering structure is estimated to be $15.6 \pm 1.6$ AU.

Significance and Claims
The paper claims that these results provide a detailed case study of an ESE, demonstrating the utility of high-cadence monitoring in resolving the temporal evolution of scattering structures. The detection of a detached feature in the secondary spectrum offers direct evidence of a double-lensing event, supporting models where ESEs arise from discrete, lens-like structures within the ISM.

The authors assert that such rare events offer unique opportunities to probe AU-scale structures in the ISM that are otherwise unresolvable. Furthermore, they note that sudden changes in scattering along a line of sight can introduce timing offsets and noise, which are relevant for Pulsar Timing Arrays (PTAs). The study concludes that continuous monitoring of scintillation parameters is a viable method for identifying ESEs, potentially enabling triggered follow-up observations to further constrain theories regarding the origins of these events.