Free Access
Issue
A&A
Volume 657, January 2022
Article Number A34
Number of page(s) 14
Section Stellar structure and evolution
DOI https://doi.org/10.1051/0004-6361/202140811
Published online 05 January 2022

© ESO 2022

1. Introduction

The physics of the pulsar magnetosphere, including the generation and propagation of the observed radio emission, remains poorly understood (e.g., Melrose & Yuen 2016; Melrose et al. 2021). A viable physical model should explain the origin of coherent emission over a broad range of frequencies, the high degree of polarisation, and a wide range of characteristic phenomena, such as nulling (e.g., Backer 1970), giant pulses (e.g., Romani & Johnston 2001; Hankins et al. 2003; Knight et al. 2006), the occurrence of orthogonally polarised modes (OPMs; e.g., Ekers & Moffet 1969; Manchester et al. 1975; Stinebring et al. 1984), and sub-pulse drifting (e.g., Drake & Craft 1968). The models proposed to date are able to explain some of the observed phenomena, but fail to describe others.

Drifting sub-pulses (Drake & Craft 1968) are observed in almost half of the pulsar population (e.g., Weltevrede et al. 2006, 2007; Basu et al. 2019). In addition to the spin period of the pulsar (P1), drifting is usually characterised by two periods: (1) the longitudinal interval or separation between successive drift bands (P2); and (2) the modulation interval, or drift cycle, expressed in pulsar rotations between drift bands (P3). Different methods exist to measure these periods. For example, Backer (1973) and Backer et al. (1975) introduce the longitude-resolved fluctuation spectrum (LRFS) method, based on one-dimensional Fourier transforms of single-pulse intensities as a function of pulse longitude. Deshpande & Rankin (2001) develops the harmonic-resolved fluctuation spectrum, and Edwards & Stappers (2002) elaborates an equivalent approach based on the two-dimensional Fourier transform. Edwards (2004, hereafter, E04) extends the LRFS to study the cross-spectral power between the four Stokes parameters, the polarisation LRFS (PLRFS). Drifting sub-pulses are also studied using the longitude-resolved cross-correlation function (e.g., Prószyński & Wolszczan 1986), two-dimensional auto-correlation function (e.g., Taylor et al. 1975), and P3 folding (e.g., Taylor et al. 1971; Edwards et al. 2003).

The origin of drifting sub-pulses is unclear. One scenario, known as the “rotating carousel” model (Ruderman & Sutherland 1975; Rankin 1993; Rankin et al. 2006), suggests that the observed radiation is produced by outflowing streams of charged particles that circulate around the magnetic pole due to E × B drift, where E and B are the electric and magnetic fields in the vicinity of the polar cap. Owing to the carousel motion, the streams drift across the primary beam of the pulsar, such that the associated sub-pulses appear to shift in pulse longitude as a function of time.

In relatively rare cases, sub-pulses are observed to periodically switch between orthogonally polarised states synchronously with the drift cycle of the modulated total intensity. OPM switching is evident as bi-modal distributions of single-pulse position angles (PAs), with peaks offset by 90° (e.g., Ekers & Moffet 1969; Manchester et al. 1975; Ramachandran et al. 2002). Many pulsars exhibit OPM switching (e.g., Stinebring et al. 1984); however, to date, only five pulsars are reported to exhibit periodic OPM switching that is synchronised with the drifting sub-pulse modulation: PSRs B0809+74 (Taylor et al. 1971; Ramachandran et al. 2002; E04), B1237+25 (Rankin & Ramachandran 2003), B0320+39, B0818−13 (E04), and B0031−07 (Ilie et al. 2020).

To explain the synchronicity of the OPM switching and drifting sub-pulses, Rankin & Ramachandran (2003) suggest a revision of the carousel model in which the emission is split into two rotating sub-beam structures, each corresponding to one of the observed OPMs. Owing to birefringence in the pulsar magnetosphere, the orthogonally polarised sub-beams are offset in both azimuthal and meridional directions with respect to the magnetic axis (Rankin et al. 2006), such that the observed switching between modes is synchronous with the modulation of total intensity caused by E × B drifting.

E04 presents examples of periodic modulations of the Stokes parameters in three pulsars with drifting sub-pulses that cannot be explained by cyclic switching between incoherently superposed OPMs (the result of which would have a linear distribution in Poincaré space). Rather, two of the three pulsars exhibit periodic fluctuations such that the tip of the Stokes polarisation vector draws an ellipse in Poincaré space. E04 discovers and studies these drift-modulated elliptical variations of the Stokes polarisation vector through eigenvalue analysis of the PLRFS. Summing the spectral density tensor of the PLRFS over all harmonics is equivalent to computing the longitude-resolved covariances between the Stokes parameters, which can be used to study stochastic fluctuations of the Stokes parameters (McKinnon 2004; Edwards & Stappers 2004). van Straten (2009) and van Straten & Tiburzi (2017) develop the mathematical foundation required to interpret the 4 × 4 matrix of covariances between the Stokes parameters by considering three different regimes of incoherent mode superposition (disjoint, superposed, and composite) and the effects of amplitude modulation.

Analysing the covariances between the Stokes parameters provides a way to extend the study of single-pulse variability, such as drifting and OPM switching, to faint pulsars. For the majority of pulsars, individual pulses cannot be detected and it is not possible to directly study single-pulse variability using scatter plots of the Stokes parameters (e.g., McKinnon 2004) or longitude-resolved histograms of polarisation state (e.g., Stinebring et al. 1984; Edwards & Stappers 2004). For the typical pulsar, there is also insufficient signal to detect the local sub-pulse drift period P3, which must be estimated in sufficiently short segments of time over which P3 can be considered stable in order to average the Stokes parameters synchronously with the drift period (e.g., Taylor et al. 1971; Weltevrede et al. 2006). Therefore, to study the majority of pulsars, it is necessary to develop statistical descriptions of sub-pulse variability that can be used to integrate over a large number of pulses. These can be applied to develop a more complete picture of single-pulse variability across the entire pulsar population and provide a deeper understanding of the physics of radio wave emission and propagation in the pulsar magnetosphere.

To prepare to study weak sources, we must first study bright pulsars, whose single-pulse variability can be analysed using methods that are applicable only at high signal-to-noise ratio (histograms, scatter plots, P3-folds, etc.) and connected with features derived from fourth-order moments (e.g., covariances between the Stokes parameters and the cross-correlations between them as a function of both longitudinal and temporal lags). These connections contribute toward the development of an interpretive framework from which we can infer the characteristics of single-pulse variability from statistical quantities that can be integrated over sufficiently long intervals.

In this article we present one further step in the development of such a framework through a case study of PSR B1919+21. As the first pulsar discovered by Jocelyn Bell (Hewish et al. 1968), PSR B1919+21 is among the most studied and, quite possibly as a consequence of such scrutiny, it has proven more exceptional than exemplary. For example, Rankin (1983) first regards its average profile as a “barely resolved conal-double”, in which the line of sight is neither tangential nor central to a hollow cone of emission. Rankin et al. (1989) later argue that PSR B1919+21 should be classified as a five-component multiple profile, based on evidence that the two main components are each composed of an unresolved pair of components produced as the line of sight crosses inner and outer cones, and that the saddle region between the two main components exhibits properties of core emission. With the caveat that there appears to be a weak trace of central core emission, Olszanski et al. (2019) classify PSR B1919+21 as having a conal quadruple profile produced by a line of sight that that cuts both inner and outer cones more centrally without passing near the core beam. With a spin-down energy loss rate of Ė ∼ 2.2 × 1031 erg s−1, the radio signal from PSR B1919+21 is expected to be dominated by conal emission (Rankin et al. 2020). Core emission that supports a five-component multiple profile classification appears to be more evident at radio frequencies below 80 MHz (Bondonneau et al. 2021).

Compounding the challenges associated with classifying its emission, PSR B1919+21 shows no evidence of the radius-to-frequency mapping thought to be associated with outer-cone component pairs (Mitra & Rankin 2002). Its main pulse is also preceded by weak emission (Weisberg et al. 1999) that appears more clearly at high radio frequencies as a separate leading component well beyond its putative outer cone (Olszanski et al. 2019). Furthermore, it is not possible to reliably constrain the magnetic field and viewing geometry of PSR B1919+21, owing to complex variations of its longitude-resolved position angle (P.A.) that depart significantly from the predictions of the rotating vector model (RVM; Radhakrishnan & Cooke 1969). For many pulsars, deviations from the RVM can be described by transitions in the dominance of orthogonally polarised modes of emission, most clearly evident as jumps of 90° in the P.A. curve (e.g., Everett & Weisberg 2001). However, at some radio frequencies, PSR B1919+21 exhibits a sharp 45° jump in P.A. near the profile peak, which can be interpreted as a sudden narrowing of the distribution of phase delay between the natural modes of wave propagation in the pulsar magnetosphere (Dyks 2019). In Fig. 18 of Mitra et al. (2015), some of the smaller distortions of the P.A. profile appear to be associated with absorption features in the total intensity that were first detected at frequencies below 200 MHz (Hankins 1973; Cordes 1975). These striking features of both the single-pulse and average profiles of PSR B1919+21 remain poorly understood.

Although the first-discovered pulsar clearly warrants further study, none of the above difficulties or exceptions motivate its selection for this work. Rather, as described in Sect. 2, PSR B1919+21 was objectively selected owing to the anomalous distribution of its polarisation state, as identified by applying novel statistical methods to publicly available data. In Sect. 3 we report on the results of applying a wide variety of methods to analyse the variability of its single-pulse polarisation state. In Sect. 4 we discuss and interpret the results, which are summarised in Sect. 5.

2. Data set

The observations used in this article are part of a larger single-pulse dataset collected at the Arecibo Radio Telescope, mainly from 1988 to 1992, at multiple frequencies, and presented and described in detail by Hankins & Rankin (2010). From the polarisation-calibrated full Stokes data summarised in Table 2 of Hankins & Rankin (2010), we selected 32 pulsars observed at one or two central frequencies, 1414 MHz and 430 MHz, as part of the P1260 program1.

After computing the eigenvalues of the longitude-resolved covariances between the Stokes parameters for each observation, we searched for behaviour that is inconsistent with incoherent OPM superposition (van Straten & Tiburzi 2017). That is, we searched for evidence of correlations between the Stokes parameters that are not consistent with a prolate spheroidal distribution of the polarisation vector. The observation of PSR B1919+21 at 1414 MHz clearly stands out as anomalous because, over a wide region of pulse longitude, the distribution of the Stokes polarisation vector is characterised by an oblate spheroid.

PSR B1919+21 produces complex patterns in the P.A. distributions of single pulse polarisation state (e.g., Fig. 18 of Mitra et al. 2015). It also exhibits sub-pulse drift with band separations P2 ≈ 10 − 17 ms and P3 ≈ 4 − 4.4 P1 (e.g., Backer 1973; Hankins 1973; Cordes 1975; Wolszczan 1978; Prószyński & Wolszczan 1986; Weltevrede et al. 2007). It has three distinct drifting components that are not well separated in pulse longitude; that is, the sub-pulses from different components drift through shared regions of pulse longitude. The definitions of the drift-band region boundaries depend on radio frequency; however, P3 is not frequency-dependent and is common to all three drift band regions (Prószyński & Wolszczan 1986). Although the single-pulse polarisation and sub-pulse drifting properties of PSR B1919+21 have been separately analysed, the connection between them has not been studied in detail2. Therefore, the remainder of our analysis focuses on PSR B1919+21.

3. Methods and results

The statistical methods applied in this work begin with an ensemble of sample mean Stokes 4-vectors

S ( ϕ , t ) = [ I ( ϕ , t ) , Q ( ϕ , t ) , U ( ϕ , t ) , V ( ϕ , t ) ] T $$ \begin{aligned} {S}(\phi ,t) = [{I}(\phi ,t),{Q}(\phi ,t), {U}(\phi ,t),{V}(\phi ,t)]^T \end{aligned} $$(1)

sampled as a function of pulse longitude ϕ, and pulsar rotation (or turn) t. Here, xT denotes the transpose of vector x; i.e., the Stokes 4-vector is treated as a column vector. The sample mean Stokes parameters are averaged over a single interval of pulse longitude, which for the PSR B1919+21 data presented in this study is ∼0.325°, or about 1.2 ms. The 4 × 4 matrix of cross-covariances between the Stokes parameters at two pulse longitudes, ϕ1 and ϕ2, as a function of time-lag in pulsar rotations τ > 0 is given by

K ( ϕ 1 , ϕ 2 ; τ ) = 1 N pulse t = τ + 1 N pulse S ¯ ( ϕ 1 , t τ ) S ¯ ( ϕ 2 , t ) T , $$ \begin{aligned} K(\phi _1, \phi _2; \tau ) = \frac{1}{N_\mathrm{pulse} } \sum _{t=\tau +1}^{N_\mathrm{pulse} } \bar{S}(\phi _1,t-\tau ) \otimes \bar{S}(\phi _2,t)^T, \end{aligned} $$(2)

where Npulse is the number of pulses, S ¯ $ \bar{S} $ represents the Stokes parameters after subtracting a local mean value, and ⊗ represents an outer product. If Kμν(ϕ1, ϕ2; τ) is the element in row μ and column ν, then the cross-covariance at negative lag is given by transposing both μ and ν, and ϕ1 and ϕ2; i.e.

K μ ν ( ϕ 1 , ϕ 2 ; τ ) = K ν μ ( ϕ 2 , ϕ 1 ; τ ) . $$ \begin{aligned} K_{\mu \nu }(\phi _1, \phi _2; -\tau ) = K_{\nu \mu }(\phi _2, \phi _1; \tau ). \end{aligned} $$(3)

3.1. Longitude-resolved covariances between the Stokes parameters

As in McKinnon (2004) and Edwards & Stappers (2004), we first studied the longitude-resolved 4 × 4 matrix of covariances between the Stokes parameters,

C ( ϕ ) = K ( ϕ , ϕ ; 0 ) $$ \begin{aligned} C(\phi ) = K(\phi , \phi ; 0) \end{aligned} $$(4)

by computing the eigendecomposition of the 3 × 3 partition of C corresponding to the three-dimensional polarisation vector, S = (Q, U, V). The eigenvectors and eigenvalues derived from this analysis describe the distribution of the polarisation vectors in Poincaré space. At a given pulse longitude ϕ, the eigenvectors define the principle axes of the distribution of the polarisation vectors; the eigenvalues λi, i ∈ {1, 2, 3}, are equal to the variances of the distribution in the directions of the corresponding eigenvectors. For example, when λ1 > λ2 = λ3, the distribution of the polarisation vectors has a prolate spheroidal shape, as expected for incoherent OPM superposition (van Straten & Tiburzi 2017).

The derived eigenvalues are shown in panel d of Fig. 1. Over the region of pulse longitude −5° ≲ ϕ ≲ 3°, the eigenvalues indicate an unexpected and uncommon oblate distribution of the polarisation state, which is the primary focus of the following section. The longitude-resolved standard deviation of the total intensity also highlights a component around longitude −12° that is most likely related to the precursor detected at 4.6 GHz (Olszanski et al. 2019). The phase-resolved modulation index of this component is ∼1.2 to 1.3, around two to three times the modulation index in the main profile, which varies between ∼0.4 and 0.7.

thumbnail Fig. 1.

Longitude-resolved mean and variability of the polarised emission from PSR B1919+21 at 1414 MHz. Panel a: both the ellipticity angle of the average Stokes parameters (points with error bars) and a histogram of the ellipticity angles of the individual pulses (greyscale intensity). Similarly, panel b: both the position angle of the average Stokes parameters and a histogram of the position angles of the individual pulses. Panel c: integrated profile in total intensity (black, upper solid), linear (red, dashed) and circular polarisation (blue, dot-dashed). These profiles were obtained by averaging together all the available data. In panel d, the black solid line shows the standard deviation of the total intensity, and the red dashed, green dot-dashed, and blue dotted lines display the square roots of the eigenvalues (from largest to smallest, respectively) of the 3 × 3 partition of the covariance matrix that describes the fluctuations of the polarisation vector S = (Q, U, V).

3.2. Distribution of the polarisation vectors

As in McKinnon (2004), we produced scatter plots of the Stokes polarisation vector using two-dimensional projections of the 3-dimensional Poincaré space. In Fig. 2, one such projected scatter plot is divided into a 25 × 25 grid and the number of points in each pixel is integrated. Rather than the expected concentration of points near the origin, the polarisation states are concentrated on the circumference of a ring that is offset from the origin.

thumbnail Fig. 2.

Integrated number of single-pulse polarisation states projected onto a two-dimensional plane that is normal to the minor axis of the oblate spheroidal distribution observed at pulse longitude 0°. The greyscale intensity represents the count of instances in each pixel after integrating 1100 single pulses; contours are spaced at intervals of 5 counts. The angle between the Stokes V axis and the normal to the Q′−U′ plane is 7.2°.

If the Stokes vectors are normalised by the total intensity, then the normalised polarisation vectors will lie within a sphere of unit radius known as the Poincaré sphere. Visualising the polarisation state in the Poincaré sphere offers a couple of advantages: first, the length of the normalised polarisation vector indicates the degree of polarisation; second, assuming that all four Stokes parameters are modulated by the same function, normalisation by Stokes I reduces additional scatter due to amplitude-modulation that is intrinsic to the pulsar and/or owing to scintillation in the interstellar medium.

Normalisation by Stokes I is possible only for very bright single-pulse observations, where the total intensity is significantly different from zero. Therefore, to verify that this normalisation is valid, we plot histograms of Stokes I in the on-pulse regions of interest. As an example, Fig. 3 shows that the total intensity is well above zero at ϕ = 0°, where the oblate distribution of the Stokes vectors is observed. The median of the off-pulse noise (i.e., the median value of the standard deviation estimated from the off-pulse longitudes of 200 consecutive pulses) varies as a function of time, most likely owing to the source moving away from the principal axis of the primary reflector. At the start of the observation, the median off-pulse noise is 75 mJy and at the end it was 134 mJy. The vertical dashed line in Fig. 3 depicts the maximum value of the standard deviation of the off-pulse noise (at the end of the observation). This figure shows that, at this pulse longitude, it is possible to visualise scatter plots of the polarisation state in either Poincaré space or within the Poincaré sphere. In Fig. 4, the projected scatter plots of the normalised polarisation vectors at ϕ = 0° show that the distribution is not an oblate spheroid, but a torus of revolution. A toroidal distribution is clearly detectable at −0.8° ≲ ϕ ≲ 0.8°, and a partial torus appears to be visible in adjacent pulse ranges −3.9° ≲ ϕ ≲ −1°, and 1° ≲ ϕ ≲ 1.6°. Scatter plots for four neighboring longitudes are presented in Fig. A.1; these show that, as a function of pulse longitude, the torus changes orientation and shifts its center. Moreover, the centers of the toroidal distributions are offset from the origin of Poincaré sphere axis at all longitudes. For comparison, we also present projected scatter plots of the normalised Stokes polarisation vector observed at ϕ ≃ −6.8° in Fig. A.2. The distribution shown in this scatter plot is consistent with Fig. 1, where the eigenvalues at this longitude indicate a prolate spheroid.

thumbnail Fig. 3.

Total intensity histogram at pulse longitude 0°, where the Stokes parameters are toroidally distributed. This example demonstrates the validity of Stokes parameters normalisation. Stokes I has a similar behaviour throughout all pulse longitudes. The vertical dashed line shows the maximum value of the median off-pulse noise.

thumbnail Fig. 4.

Distribution of the normalised Stokes parameters in the Poincaré sphere at ϕ 0.°. The two upper panels and the lower-left panel show the projected, 2D distribution of the Stokes parameters in the three planes normal to the basis vectors, and the 3D distribution is displayed in the lower-right panel.

3.3. Correlations of Stokes parameters

Together with the 3D distributions of the Stokes parameters it is important to know whether the torus formed in an orderly or a random way. To this end, we explored the spatial and temporal relationships between the polarisation vectors via

  1. The longitude-offset cross-correlation function of the Stokes parameters with respect to a reference longitude; and

  2. The longitude-resolved cross-correlation functions between two different Stokes parameters.

For each element of R(ϕ1, ϕ2; τ), we computed the Pearson correlation coefficient,

R μ ν ( ϕ 1 , ϕ 2 ; τ ) = K μ ν ( ϕ 1 , ϕ 2 ; τ ) σ μ ( ϕ 1 ) σ ν ( ϕ 2 ) $$ \begin{aligned} R_{\mu \nu }(\phi _1, \phi _2; \tau ) = \frac{K_{\mu \nu }(\phi _1, \phi _2; \tau )}{\sigma _\mu (\phi _1)\sigma _\nu (\phi _2)} \end{aligned} $$(5)

where μ, ν ∈ {I, Q, U, V}, and the longitude-resolved variance is given by

σ μ 2 ( ϕ ) = C μ μ ( ϕ ) = K μ μ ( ϕ , ϕ ; 0 ) $$ \begin{aligned} \sigma ^2_\mu (\phi ) = C_{\mu \mu }(\phi ) = K_{\mu \mu }(\phi , \phi ; 0) \end{aligned} $$(6)

Prior to computing R(ϕ1, ϕ2; τ), we subtracted the moving average of the Stokes parameters from each pulse longitude, computed with a boxcar size of 100 pulse periods, which is smaller than the scintillation time at 1.4 GHz (∼500 to 1500 s; Cordes 1986; Bhat et al. 1999).

3.3.1. Longitude-offset cross-correlations for each Stokes parameter

Figures 5 and 6 present the longitude-offset cross-correlation coefficients,

R μ ( Δ ϕ , τ ; ϕ 0 ) = R μ μ ( ϕ 0 , ϕ 0 + Δ ϕ ; τ ) , $$ R^\prime _\mu (\Delta \phi ,\tau ;\phi _0) = R_{\mu \mu }(\phi _0, \phi _0+\Delta \phi ; \tau ), $$

thumbnail Fig. 5.

Longitude-offset cross-correlations of each Stokes parameter as a function of lag τ and longitude offset Δϕ with respect to the reference pulse longitude ϕ0 ≃ −6.8°. In each panel, a horizontal black line marks zero lag and a vertical black line marks zero longitude offset.

thumbnail Fig. 6.

As for Fig. 5, but for reference longitude ϕ0 = 0°, in the middle of the region where the toroidal distribution of the polarisation vector in the Poincaré sphere is observed.

with μ ∈ {I, Q, U, V}, as a function of lag τ and longitude offset Δϕ with respect to the reference pulse longitude ϕ0.

The reference longitude of Fig. 5 is ϕ0 ≃ −6.8°, which is within a region of pulse longitude that is characterised by a prolate spheroidal distribution of the Stokes parameters that is approximately centered at the origin of the Poincaré sphere. Figure 5 shows three main drift regions in total intensity, as first noted by Prószyński & Wolszczan (1986). In this plot, the first drift-band lies between the beginning of the on-pulse region and a pulse longitude of approximately −5.4° (corresponding to a longitude offset of 1.4°), the second drift band starts at the end of the first band and spans up to a longitude of about −2.6° (longitude offset 4.5°), and the third drift band lies approximately between longitudes −0.6° and 3°.

The periodic modulation of total intensity at the reference longitude is remarkably correlated with the modulation in the rest of the on-pulse region, even as far away as longitude offset Δϕ ≃ 9.7°. However, this correlation is diminished over the saddle region 5° ≲ Δϕ ≲ 7° (corresponding to longitudes −1.8° ≲ ϕ ≲ 0.2°). This region of diminished correlation in total intensity partially overlaps the region (5.5° ≲ Δϕ ≲ 7.5°) where the periodic modulations of Stokes Q and U exhibit little or no drifting (i.e., the slope is flat). The periodic modulations of Stokes Q, U, and V remain strong in this region.

The reference longitude of Fig. 6 is ϕ0 = 0°, which is in the center of the longitude region where the toroidal distribution of Stokes (Q, U, V) in the Poincaré sphere is observed. Figure 6 reaffirms that the total intensity in this region is not highly correlated with intensity in the rest of the on-pulse region. In both Figs. 5 and 6, the longitudinal discontinuity in the modulation pattern of Stokes Q at ϕ ≃ −5.4° occurs at the transition between OPMs that is also evident in Fig. 1.

3.3.2. Longitude-resolved Stokes Q, U, V cross-correlations

Figure 7 presents the longitude-resolved cross-correlation coefficients,

R μ ν ( ϕ , τ ) = R μ ν ( ϕ , ϕ ; τ ) , $$ R^\prime _{\mu \nu } (\phi ,\tau ) = R_{\mu \nu }(\phi , \phi ; \tau ), $$

thumbnail Fig. 7.

Longitude-resolved Stokes Q, U, V cross-correlations as a function of lag τ and pulse longitude ϕ. From top to bottom are plotted the correlations between Stokes Q and U, Stokes U and V, and Stokes V and Q. In each panel, a horizontal black line marks zero lag.

with (μ, ν)∈{(Q, U),(U, V),(V, Q)}, as a function of lag τ and pulse longitude ϕ. Owing to symmetry, these three panels depict all of the information contained in the 6 off-diagonal components of the 4 × 4 matrix, R′(ϕ, τ); furthermore, taking the one-dimensional Fourier transforms of the components of R′(ϕ, τ) along the τ axis yields the PLRFS studied by E04. As in the PLRFS, the correlations between the Stokes parameters indicate their cyclic oscillation about an elliptical path in Poincaré space. They also indicate that the normal to this elliptical path precesses as a function of pulse longitude.

For example, over −1.5° ≲ ϕ ≲ 1°, R QU (ϕ,τ) $ R^\prime_{{QU}}(\phi,\tau) $ is antisymmetric about τ = 0, which indicates that the oscillations of Stokes Q and U are 90° out of phase. This is consistent with the periodic cycling of the polarisation vector around an elliptical path in a plane that is only slightly offset from the QU plane3 over this range of pulse longitudes.

At earlier pulse longitudes within this region, e.g., −1.5° ≲ ϕ ≲ 0°, R VQ (ϕ,τ) $ R^\prime_{{VQ}}(\phi,\tau) $ is negative and symmetric about τ = 0, which indicates that Stokes Q and V are anti-correlated (180° out of phase). This is consistent with the anti-correlation that is visible in the QV scatter plots in Fig. 4 and the upper panels of Fig. A.1. Furthermore, the oscillations of Stokes U and V are 90° out of phase, as indicated by the antisymmetry of R UV (ϕ,τ) $ R^\prime_{{UV}}(\phi,\tau) $. Both the anticorrelation of Stokes Q and V and the antisymmetry of Stokes U and V are consistent with an axis of revolution of the torus that lies primarily in the QV plane, such that there is an apparently elliptical path projected into the UV plane.

Similarly, at later pulse longitudes within this region, e.g., 0° ≲ ϕ ≲ 1°, R UV (ϕ,τ) $ R^\prime_{{UV}}(\phi,\tau) $ is negative and symmetric about τ = 0, which indicates that Stokes U and V are anti-correlated. This is consistent with the anti-correlation that is visible in the UV scatter plots in Fig. 4 and the lower panels of Fig. A.1. Furthermore, the oscillations of Stokes Q and V are 90° out of phase, as indicated by the antisymmetry of R VQ (ϕ,τ) $ R^\prime_{{VQ}}(\phi,\tau) $. Both the anticorrelation of Stokes U and V and the antisymmetry of Stokes Q and V are consistent with an axis of revolution of the torus that lies primarily in the UV plane, such that there is an apparently elliptical path projected into the QV plane. That is, as a function of pulse longitude, the axis of revolution of the torus precesses about the Stokes V axis, rotating from primarily in the QV plane to primarily in the UV plane.

3.4. Average drift band

As in Taylor et al. (1971) and Edwards et al. (2003), we estimated the polarisation state of the mean drift bands by averaging the single-pulse Stokes parameters and quantities derived from them synchronously with the modulation period, P3 ≃ 4.2P1, which was determined via short-time Fourier transform analysis. The results of P3 folding the first 500 pulses in the data set are shown in Fig. 8, which exhibits several noteworthy features.

thumbnail Fig. 8.

P3-folded polarisation state of the first 500 pulses in the PSR B1919+21 dataset. From top to bottom are the total intensity; the position angle, θ = tan−1(U′/Q′)/2; the fractional linear polarisation, L ̂ = Q 2 + U 2 / I $ \hat{L}=\sqrt{Q^2+U^2}/I $; the ellipticity angle, ϵ = sin−1(V′/L′)/2; and the fractional circular polarisation, C ̂ = | V | / I $ \hat{C}=|V|/I $; where P3-folded quantities are shown with primed symbols (e.g., Q′ and U′) and the symbols for single-pulse quantities are not primed (e.g., Q and U). That is, L ̂ $ \hat{L} $ and C ̂ $ \hat{C} $ are derived from the single-pulse Stokes parameters before P3 folding and θ and ϵ are derived from the P3-folded Stokes parameters. Furthermore, L ̂ $ \hat{L} $ and C ̂ $ \hat{C} $ are integrated only when the single-pulse total intensity exceeds 4 times the standard deviation of the of noise. In each panel, the modulation period (P3 ≃ 4.2 P1) is plotted twice to better illustrate the continuity of the pattern.

In all panels, pulse longitude ϕ ≃ −5.4° marks several significant changes in the sub-pulse modulation pattern. First, there is a 180° phase shift in the modulation of the total intensity associated with the transition from region 1 to region 2 of the mean drift band. The overlap of drift regions 1 and 2 leads to an apparent doubling of the modulation frequency, which is also detected in Fig. A.21 of Weltevrede et al. (2007), where it is described as “P3 = 2P0 flickering”. Overlapping and offset drift regions also explain the reduced the longitude-resolved modulation index observed over this region of pulse longitude (Weltevrede et al. 2006).

The transition from drift region 1 to drift region 2 is accompanied by a longitudinal transition between orthogonally polarised modes, which is also clearly visible in the position angle histogram plotted in panel b of Fig. 1. Before this transition (over −6.3° ≲ ϕ ≲ −5.4°), the position angle exhibits little modulation and has a mean value close to 0°. Immediately after this transition (over −5.4° ≲ ϕ ≲ −4.7°), the position angle varies between ∼ − 60° and ∼75° as a function of modulation phase. The minimum degrees of polarisation (both fractional linear and fractional circular) in Fig. 8 occur just before ϕ ≃ −5.4°. Shortly after this transition (over −5.0° ≲ ϕ ≲ −3.7°) the ellipticity angle exhibits minimal modulation and remains greater than 15°; correspondingly, the fractional circular polarisation also reaches its maximum value in this pulse longitude range.

In the longitudes shortly before the OPM switch at ϕ ≃ −5.4° (between approximately −7.6° and −6.3°) the P3-folded position angle and ellipticity angle indicate variation between orthogonally polarised states as a function of modulation phase. However, as first noted for PSR B0809+74 (Taylor et al. 1971), the position angle does not jump discontinuously from one mode to the other, as expected for incoherent OPM superposition; rather, there is a continuous decrease followed by a continuous increase in position angle as a function of modulation phase. With reference to Fig. A.2, the apparent continuous oscillation of the polarisation state could be explained by periodic switching along an axis that does not pass through the origin in the QU plane. A likely origin for such off-center switching is the incoherent transition between slightly non-orthogonal polarisation modes (e.g., McKinnon 2003). Similarly, continuous oscillation of position angle as a function of modulation phase is observed in the longitude range following the OPM transition, between approximately −5.4° and −3.6°. Here, there is also a clear systematic swing of position angle through each sub-pulse as a function of longitude and, as first noted for PSR B0809+74 (Taylor et al. 1971), this swing follows the sub-pulse structure as it drifts in longitude as a function of modulation phase. This swing gives rise to the continuous oscillation of position angle at a fixed longitude as a function of modulation phase.

Over the pulse longitudes that exhibit the toroidal and partially toroidal distributions (approximately −3.6° ≲ ϕ ≲ 1.3°) the position angle no longer oscillates back and forth; rather, it cycles monotonically through the full spectrum of visible colours, clearly indicating the elliptical motion of the polarisation vector in the QU plane. Here, the degree of linear polarisation is also a maximum. Owing to the high degree of linear polarisation and the high signal-to-noise ratio of the individual pulses, the continuous cycling of the position angle over this range of pulse longitudes can also be clearly seen in the animations provided online.

4. Discussion

In light of the complexities outlined in the introduction, it would be premature to interpret our results in terms of magnetospheric emission physics without a more comprehensive overview of single-pulse polarimetry across the radio spectrum. Instead, we aim to add valuable pieces to the puzzle and focus our discussion on two key features of the polarisation of the drifting sub-pulses: (1) the dramatic change in the modulation-resolved polarisation pattern seen in Fig. 8 that coincides with the longitudinal transition between orthogonally polarised modes near ϕ ≃ −5.4°; and (2) the toroidal distributions and cyclic drifting of polarisation state around an elliptical path as a function of modulation phase, observed over the region −3.6° ≲ ϕ ≲ 1.3°.

First, the change in P3-modulated polarisation near a pulse longitude of ϕ ∼ −5.4° can be explained as a transition between primarily incoherent (ϕ ≲ −5.4°) and partially coherent (ϕ ≳ −5.4°) mode superposition. At ϕ ≲ −5.4°, the distribution of polarisation states in the Poincaré sphere has, to first order, a prolate spheroidal shape as expected for the incoherent superposition of OPMs (e.g., Fig. A.2, van Straten & Tiburzi 2017). Over the region roughly defined by −5.0° ≲ ϕ ≲ −3.7°, the polarisation state varies as a function of modulation phase between a low degree of polarisation (both linear and circular), as expected for an incoherent superposition of OPMs, and a high degree of circular polarisation, as expected for a coherent superposition of linearly polarised natural modes in a highly relativistic plasma (e.g., Kennett & Melrose 1998). Further evidence of coherent superposition is displayed on either side of this region. Over −5.4° ≲ ϕ ≲ −4.7°, the position angle varies between approximately −60° and 75°; and over −3.7° ≲ ϕ ≲ −3.6°, the position angle varies between approximately 30° and 75°. These approximately ±45° offsets can also be interpreted as transitions between incoherent OPM superposition, dominated by one mode, and a coherent superposition of linearly polarised modes (Dyks 2019). The sign of the position angle offset is opposite on either side of the region where the degree of circular polarisation reaches its maximum, which can be interpreted as monotonic variation of the phase of the coherent superposition of linearly polarised natural modes as a function of drift phase (and hence pulse longitude), as also considered for PSR B1451−68 (Dyks et al. 2021).

Over −0.8° ≲ ϕ ≲ 0.8°, the single-pulse polarisation state exhibits remarkable longitude-resolved toroidal distributions that are clear in both the integrated scatter plot shown in Fig. 2 and the normalised scatter plots shown in Figs. 4 and A.1. Furthermore, the polarisation state cycles around the elliptical path of the torus in a continuous periodic manner, rather than a stochastic one, modulated in phase with the drifting sub-pulse pattern observed in the total intensity. Cycling of the polarisation state around an elliptical path is demonstrated by the cross-correlations of Stokes Q and U in Fig. 7 and in the P3-folded position angle in Fig. 8. The longitude-resolved cross-correlations between the Stokes parameters also indicate that the axis of elliptical revolution precesses about the Stokes V axis as a function of pulse longitude.

To date, continuous cyclic variation of the Stokes polarisation vector around an elliptical locus in Poincaré space has been observed in only three other sources, PSR B0320+39, PSR B0809+74, and PSR B0818−13 (observed at ∼360 MHz, 1380 MHz, and 330 MHz, respectively; E04). For PSR B0320+39, E04 reports significant periodic fluctuations of the polarisation vector about an elliptical locus over a range of pulse longitude (∼18° to ∼19° in Fig. 4 of E04) where periodic modulations of the total intensity are significantly diminished. In our analysis of PSR B1919+21, barely detectable modulation of total intensity is observed in the “saddle region” (−2.3° ≲ ϕ ≲ −0.3°; cf. panels c and d of Fig. 1 just prior to the range of pulse longitudes where the toroidal distributions are observed (−0.8° ≲ ϕ ≲ 0.8°). As noted by Backer (1973), the decreased intensity modulation in this saddle region could be due to overlapping sub-pulses from the drift band regions on either side, especially if they differ by 180° in modulation phase. For PSR B0818−13, the elliptical locus reported by E04 appears to occur in a plane that is perpendicular to the Stokes QU plane and above it, as inferred from the positive ellipticity angles at pulse longitudes between ∼38° and ∼40° in Fig. 2 of E04 and the meridional distribution of polarisation states in Fig. 5 of E04. This is quite different to the toroidal distribution of polarisation states observed in our analysis of PSR B1919+21, for which the angle between the axis of revolution and the Stokes V axis varies between ∼15° and ∼45° (cf. Fig. A.1).

4.1. Origin of the toroidal distribution of polarisation state

We consider two plausible explanations for the observed toroidal distributions of polarisation state. These include (1) generalised Faraday rotation; and (2) four-mode mixing.

Generalised Faraday rotation. As in previous studies (e.g., Cheng & Ruderman 1979; Kennett & Melrose 1998; Edwards & Stappers 2004; Jones 2016; Dyks 2017), consider a coherent superposition of orthogonally polarised natural modes, which have anti-parallel Stokes polarisation vectors that define an axis in Poincaré space. Phase-coherent superposition of the natural modes produces a new polarisation state that lies in a plane that is normal to this axis. The azimuthal direction of the new polarisation state in this normal plane depends on the phase relation between the orthogonal modes. (In the special case of Faraday rotation, the natural modes have left-handed and right-handed circular polarisation, and birefringence-induced differential phase between these modes rotates the polarisation vector about the Stokes V axis.) If the phase difference between the modes varies smoothly between 0 and 2π during one drift cycle, then the polarisation state sweeps through a circle in the normal plane, creating a toroidal distribution.

Four-mode mixing. Consider the incoherent superposition of polarised radiation from two emission regions, ℳ1 and ℳ2. The emission from each region oscillates periodically between two orthogonally polarised states: M1 and −M1 from region ℳ1, and M2 and −M2 from region ℳ2, such that the radiation consists of an incoherent superposition of four modes in total (in contrast to the coherent superposition of four modes considered by Dyks et al. 2021). The primary axes defined by these pairs of OPMs in the Poincaré sphere are offset by some angle close to 90° (e.g., the emission from region ℳ1 could oscillate between positive and negative Stokes Q and the emission from region ℳ2 could oscillate between positive and negative Stokes U). On its own, amplitude modulation of a single pair of incoherently superposed OPMs would create a prolate distribution of polarisation states with its primary axis aligned with the axis defined by the OPMs in Poincaré space.

If the two pairs of OPM oscillations are 90° out of phase with respect to each other, then their incoherent superposition will generate a toroidal distribution of polarisation vectors in the Poincaré sphere (as depicted in Fig. A.3).

4.2. Origin of drifting sub-pulses

To explain the drifting of the observed sub-beams, we consider two physical models. These include (1) the refractive steering model; and (2) the rotating carousel model.

Refractive steering model. In this model, sub-beams are steered across the primary beam owing to refraction in the pulsar magnetosphere (e.g., Barnard & Arons 1986; Lyubarskii & Petrova 1998; Petrova & Lyubarskii 2000). Each sub-beam consists of a coherent superposition of natural modes and, as the beam is steered, the relative phase difference between the modes slowly varies, causing the resulting polarisation state to move in an arc around the primary axis defined by the OPM pair. Because the magnetospheric conditions that determine the refracted ray path of each sub-beam (e.g., Fussell & Luo 2004; Petrova 2000) may also determine the polarisation state (e.g., through birefringent effects Cheng & Ruderman 1979; Lyubarskii & Petrova 1999), it seems plausible that periodic modulation of these magnetospheric conditions would result in temporal correlation between the angle at which the refracted sub-beam leaves the pulsar magnetosphere and the polarisation state of the sub-beam. However, more detailed modeling is required to explain how the relative phase between the natural modes cycles through approximately 2π per drift cycle. Furthermore, at pulse longitudes −0.8° ≲ ϕ ≲ 0.8°, the polarisation state rotates in a plane that is offset by only a small angle from the QU plane. This seems to suggest that the natural modes of the magnetospheric plasma that is responsible for the birefringence are highly circularly polarised, which is inconsistent with the properties of a highly relativistic plasma for which the natural modes are expected to be linearly polarised (e.g., Kennett & Melrose 1998). This inconsistency could be resolved by introducing two different birefractive regions in the pulsar magnetosphere, both with linearly polarised natural modes, such that the first region introduces the torus and the second region rotates the primary axis of the torus out of the QU plane.

Rotating carousel model. In the rotating carousel model, charged particle streams circulate around the magnetic poles owing to E×B drift, causing sub-beams to drift across the primary beam (e.g., Ruderman & Sutherland 1975; Rankin 1993; Rankin et al. 2006). Rankin & Ramachandran (2003) extend this model such that each sub-beam is split (e.g., owing to birefringence) into two spatially separated beams; as our line of sight passes from one sub-beam to the next, the polarisation state switches between OPMs. For elliptical cycling of polarisation state to arise requires a fortuitous incoherent superposition of two distinct rotating carousels, one carousel comprised of the OPM pair from region ℳ1 and the other comprised of the OPMs from region ℳ2 (as described in the four mode mixing model). The sub-beams in the ℳ1 carousel must also be offset in magnetic longitude from the sub-beams in the ℳ2 carousel, such that the polarisation state periodically cycles through the four modes and completes a full 360° rotation with each drift cycle. It is important to note that, if the sub-beams are as well resolved as shown in Fig. 6 of Rankin & Ramachandran (2003), then the total intensity would cycle through four maxima in the time required for the polarisation state to complete one circle. This would not be consistent with our observations; e.g., Fig. 6 shows that the total intensity is only weakly modulated compared to the much deeper modulations of Stokes Q, U, and V. This can be explained if the sub-beams observed in the toroidal region are unresolved (i.e. the width of each sub-beam is comparable to the spacing between them) such that any remaining modulation in intensity arises from the (possibly small) differences in the peak intensities of the beams along the line of sight that cuts through them (as depicted in Fig. A.3).

In principle, four-mode mixing by a pair of unresolved or partially resolved carousels is consistent with either the barely resolved conal quadruple (Olszanski et al. 2019) or the five-component multiple (Rankin et al. 1989) classifications that have been applied to describe the mean profile of PSR B1919+21. In this interpretation, the first and third drift-bands represent the leading and trailing edges of the outer cone, the second drift-band represents the leading edge of the inner cone, the trailing edge of the inner cone is unresolved, and the inner and outer cones represent the ℳ1 and ℳ2 carousels of orthogonally polarised sub-beams. As noted in the introduction, quantitatively testing this interpretation of the beam geometry is challenging.

First, it is not possible to model the entire P.A. curve, plotted in panel b of Fig. 1, using the RVM. Therefore, to apply the RVM it is necessary to apply judgement and model only limited regions of pulse longitude (e.g., Everett & Weisberg 2001). We identify two candidate regions: A) the steepest, negative slope of the P.A. profile ranging from −5° to −2.5° in pulse longitude; and B) the shallow, positive slope from −13° to +2°, excluding −7.5° to −2°. Region A exhibits the highest degree of circular polarisation and the P.A. in this region most likely represents a distortion produced by coherent superposition of natural modes rather than a regular RVM swing (Dyks 2019). Moreover, assuming that the leading component centred at −12° can be neglected as anomalous, the steepest slope of P.A. swing in region A leads the centre of the total intensity profile, which conflicts with the expected effects of aberration and retardation (Blaskiewicz et al. 1991). Regardless of the selected region, the best-fit value for the colatitude of the magnetic axis α approaches 0°, which is inconsistent with the observed narrow pulse width. Therefore, we abandon any further interpretation of the beam geometry based on the P.A. profile.

To measure the apparent longitudinal offset between the centres of the putative inner and outer cones of emission, we model the mean total intensity pulse profile as a sum of six components, each described by a von Mises function (e.g., Weltevrede & Johnston 2008). The derived difference of 0.36° in relativistic phase shifts corresponds to an offset in the inner and outer cone emission heights of ∼200 km (Dyks et al. 2004); therefore, relativistic aberration is a plausible explanation for the overlapping trailing components of the inner and outer cones.

5. Conclusions

The brightness of PSR B1919+21, combined with the sensitivity of Arecibo, enables a variety of approaches to studying the polarisation of its drifting sub-pulses. The methods adopted in this study facilitated the discovery of toroidal and partially toroidal distributions of the Stokes parameters. Further investigation showed that, over the range of pulse longitudes that exhibit toroidal distributions, the polarisation state rotates continuously and synchronously with the sub-pulse drift period.

This analysis is one step forward in the development of a framework that can be used to interpret the fourth-order statistics of the Stokes parameters in weaker sources. Eigenvalues can be used to identify rare oblate distributions of the Stokes polarisation vectors, and the cross-correlations between the Stokes parameters can reveal the continuous cycling of polarisation state around an elliptical path, synchronous with the intensity modulation of drifting sub-pulses. In the Arecibo observations of PSR B1919+21, these interpretations are supported by visualisation of the single pulses using methods that are applicable only at high S/N. For example, the scatter plots of the normalised Stokes parameters clearly show the toroidal distribution of polarisation states and the P3-folded position angle confirms the continuous cycling of polarisation state.

Our study demonstrates that novel statistical methods can facilitate the discovery of relatively rare phenomena, even when analysing decades-old archival observations of the most studied pulsars. Although only 6 pulsars are currently known to exhibit drift-synchronous modulation of polarisation state, we speculate that more such sources remain to be discovered. To obtain a large sample of uniformly studied pulsars, we are conducting the Pulsar Radio Emission Statistics Survey (PRESS)4, an observing campaign that will record single-pulse observations of around 200 pulsars using the ultra-wide bandwidth low-frequency (UWL) receiver at the Parkes Observatory (Hobbs et al. 2020). These data will be made public as soon as they are calibrated and various quality assurance checks have been performed. All of the software required to perform the analysis presented in this paper is freely available as part of PSRCHIVE, an open-source C++ development library for the analysis of pulsar astronomical data (Hotan et al. 2004; van Straten et al. 2012).

Supplementary Material

Animations Access here


1

These data are located at http://www.uvm.edu/~pulsar/P1260

2

Manchester et al. (1975) noted the possibility of correlation between polarisation state and the drifting sub-pulse structures.

3

The offset varies from pulse longitude to pulse longitude.

Acknowledgments

The authors dedicate this paper to Professor Jocelyn Bell-Burnell, who discovered the first pulsar, PSR B1919+21, in 1967. The authors thank Joanna Rankin for providing the software required to interpret the public data of Hankins & Rankin (2010); this software is now incorporated into PSRCHIVE. We also thank Joanna Rankin for providing valuable comments that greatly improved this paper. C. Tiburzi is supported by Veni grant (project number 016.Veni.192.086) awarded by the Dutch Research Council (NWO). J. Dyks was supported by grant 2017/25/B/ST9/00385 of National Science Centre, Poland. Figures 5 through 8 of this paper use the perceptually uniform colour maps developed by Kovesi (2015) (https://colorcet.com/). At the time of recording the observations presented in this paper, the Arecibo Observatory was part of the National Astronomy and Ionosphere Center, operated by Cornell University for the National Science Foundation. The Arecibo Observatory is currently a facility of the National Science Foundation operated under cooperative agreement by the University of Central Florida and in alliance with Universidad Ana G. Mendez, and Yang Enterprises, Inc. We regret that, while drafting this manuscript, the Arecibo radio telescope’s instrument platform crashed into its primary reflector after the cables supporting the platform catastrophically snapped.

References

  1. Backer, D. C. 1970, Nature, 228, 42 [NASA ADS] [CrossRef] [Google Scholar]
  2. Backer, D. C. 1973, ApJ, 182, 245 [NASA ADS] [CrossRef] [Google Scholar]
  3. Backer, D. C., Rankin, J. M., & Campbell, D. B. 1975, ApJ, 197, 481 [NASA ADS] [CrossRef] [Google Scholar]
  4. Barnard, J. J., & Arons, J. 1986, ApJ, 302, 138 [NASA ADS] [CrossRef] [Google Scholar]
  5. Basu, R., Mitra, D., Melikidze, G. I., & Skrzypczak, A. 2019, MNRAS, 482, 3757 [NASA ADS] [CrossRef] [Google Scholar]
  6. Bhat, N. D. R., Rao, A. P., & Gupta, Y. 1999, ApJS, 121, 483 [NASA ADS] [CrossRef] [Google Scholar]
  7. Blaskiewicz, M., Cordes, J. M., & Wasserman, I. 1991, ApJ, 370, 643 [Google Scholar]
  8. Bondonneau, L., Grießmeier, J. M., Theureau, G., et al. 2021, A&A, 652, A34 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  9. Cheng, A. F., & Ruderman, M. 1979, ApJ, 229, 348 [NASA ADS] [CrossRef] [Google Scholar]
  10. Cordes, J. M. 1975, ApJ, 195, 193 [NASA ADS] [CrossRef] [Google Scholar]
  11. Cordes, J. M. 1986, ApJ, 311, 183 [NASA ADS] [CrossRef] [Google Scholar]
  12. Deshpande, A. A., & Rankin, J. M. 2001, MNRAS, 322, 438 [NASA ADS] [CrossRef] [Google Scholar]
  13. Drake, F. D., & Craft, H. D. 1968, Nature, 220, 231 [NASA ADS] [CrossRef] [Google Scholar]
  14. Dyks, J. 2017, MNRAS, 472, 4598 [NASA ADS] [CrossRef] [Google Scholar]
  15. Dyks, J. 2019, MNRAS, 488, 2018 [NASA ADS] [CrossRef] [Google Scholar]
  16. Dyks, J., Rudak, B., & Harding, A. K. 2004, ApJ, 607, 939 [NASA ADS] [CrossRef] [Google Scholar]
  17. Dyks, J., Weltevrede, P., & Ilie, C. 2021, MNRAS, 501, 2156 [NASA ADS] [CrossRef] [Google Scholar]
  18. Edwards, R. T. 2004, A&A, 426, 677 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Edwards, R. T., & Stappers, B. W. 2002, A&A, 393, 733 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  20. Edwards, R. T., & Stappers, B. W. 2004, A&A, 421, 681 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  21. Edwards, R. T., Stappers, B. W., & van Leeuwen, A. G. J. 2003, A&A, 402, 321 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Ekers, R. D., & Moffet, A. T. 1969, ApJ, 158, L1 [NASA ADS] [CrossRef] [Google Scholar]
  23. Everett, J. E., & Weisberg, J. M. 2001, ApJ, 553, 341 [NASA ADS] [CrossRef] [Google Scholar]
  24. Fussell, D., & Luo, Q. 2004, MNRAS, 349, 1019 [NASA ADS] [CrossRef] [Google Scholar]
  25. Hankins, T. H. 1973, ApJ, 181, L49 [NASA ADS] [CrossRef] [Google Scholar]
  26. Hankins, T. H., & Rankin, J. M. 2010, AJ, 139, 168 [NASA ADS] [CrossRef] [Google Scholar]
  27. Hankins, T. H., Kern, J. S., Weatherall, J. C., & Eilek, J. A. 2003, Nature, 422, 141 [CrossRef] [Google Scholar]
  28. Hewish, A., Bell, S. J., Pilkington, J. D. H., Scott, P. F., & Collins, R. A. 1968, Nature, 217, 709 [NASA ADS] [CrossRef] [Google Scholar]
  29. Hobbs, G., Manchester, R. N., Dunning, A., et al. 2020, PASA, 37, e012 [NASA ADS] [CrossRef] [Google Scholar]
  30. Hotan, A. W., van Straten, W., & Manchester, R. N. 2004, PASA, 21, 302 [Google Scholar]
  31. Ilie, C. D., Weltevrede, P., Johnston, S., & Chen, T. 2020, MNRAS, 491, 3385 [Google Scholar]
  32. Jones, P. B. 2016, MNRAS, 455, 3814 [NASA ADS] [CrossRef] [Google Scholar]
  33. Kennett, M., & Melrose, D. 1998, PASA, 15, 211 [NASA ADS] [CrossRef] [Google Scholar]
  34. Knight, H. S., Bailes, M., Manchester, R. N., Ord, S. M., & Jacoby, B. A. 2006, ApJ, 640, 941 [NASA ADS] [CrossRef] [Google Scholar]
  35. Kovesi, P. 2015, ArXiv e-prints [arXiv:1509.03700] [Google Scholar]
  36. Lyubarskii, Y. E., & Petrova, S. A. 1998, A&A, 333, 181 [NASA ADS] [Google Scholar]
  37. Lyubarskii, Y. E., & Petrova, S. A. 1999, Ap&SS, 262, 379 [Google Scholar]
  38. Manchester, R. N., Taylor, J. H., & Huguenin, G. R. 1975, ApJ, 196, 83 [NASA ADS] [CrossRef] [Google Scholar]
  39. McKinnon, M. M. 2003, ApJ, 590, 1026 [NASA ADS] [CrossRef] [Google Scholar]
  40. McKinnon, M. M. 2004, ApJ, 606, 1154 [CrossRef] [Google Scholar]
  41. Melrose, D. B., & Yuen, R. 2016, J. Plasma Phys., 82, 635820202 [CrossRef] [Google Scholar]
  42. Melrose, D. B., Rafat, M. Z., & Mastrano, A. 2021, MNRAS, 500, 4530 [Google Scholar]
  43. Mitra, D., & Rankin, J. M. 2002, ApJ, 577, 322 [NASA ADS] [CrossRef] [Google Scholar]
  44. Mitra, D., Arjunwadkar, M., & Rankin, J. M. 2015, ApJ, 806, 236 [NASA ADS] [CrossRef] [Google Scholar]
  45. Olszanski, T. E. E., Mitra, D., & Rankin, J. M. 2019, MNRAS, 489, 1543 [NASA ADS] [CrossRef] [Google Scholar]
  46. Petrova, S. A. 2000, A&A, 360, 592 [NASA ADS] [Google Scholar]
  47. Petrova, S. A., & Lyubarskii, Y. E. 2000, A&A, 355, 1168 [NASA ADS] [Google Scholar]
  48. Prószyński, M., & Wolszczan, A. 1986, ApJ, 307, 540 [CrossRef] [Google Scholar]
  49. Radhakrishnan, V., & Cooke, D. J. 1969, Astrophys. Lett., 3, 225 [NASA ADS] [Google Scholar]
  50. Ramachandran, R., Rankin, J. M., Stappers, B. W., Kouwenhoven, M. L. A., & van Leeuwen, A. G. J. 2002, A&A, 381, 993 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  51. Rankin, J. M. 1983, ApJ, 274, 359 [Google Scholar]
  52. Rankin, J. M. 1993, ApJ, 405, 285 [Google Scholar]
  53. Rankin, J. M., & Ramachandran, R. 2003, ApJ, 590, 411 [NASA ADS] [CrossRef] [Google Scholar]
  54. Rankin, J. M., Stinebring, D. R., & Weisberg, J. M. 1989, ApJ, 346, 869 [NASA ADS] [CrossRef] [Google Scholar]
  55. Rankin, J. M., Ramachandran, R., van Leeuwen, J., & Suleymanova, S. A. 2006, A&A, 455, 215 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  56. Rankin, J. M., Olszanski, T. E. E., & Wright, G. A. E. 2020, ApJ, 890, 151 [NASA ADS] [CrossRef] [Google Scholar]
  57. Romani, R., & Johnston, S. 2001, ApJ, 557, L93 [NASA ADS] [CrossRef] [Google Scholar]
  58. Ruderman, M. A., & Sutherland, P. G. 1975, ApJ, 196, 51 [Google Scholar]
  59. Stinebring, D. R., Cordes, J. M., Rankin, J. M., Weisberg, J. M., & Boriakoff, V. 1984, ApJS, 55, 247 [NASA ADS] [CrossRef] [Google Scholar]
  60. Taylor, J. H., Huguenin, G. R., Hirsch, R. M., & Manchester, R. N. 1971, Astrophys. Lett., 9, 205 [NASA ADS] [Google Scholar]
  61. Taylor, J. H., Manchester, R. N., & Huguenin, G. R. 1975, ApJ, 195, 513 [NASA ADS] [CrossRef] [Google Scholar]
  62. van Straten, W. 2009, ApJ, 694, 1413 [NASA ADS] [CrossRef] [Google Scholar]
  63. van Straten, W., & Tiburzi, C. 2017, ApJ, 835, 293 [NASA ADS] [CrossRef] [Google Scholar]
  64. van Straten, W., Demorest, P., & Osłowski, S. 2012, Astron. Res. Technol., 9, 237 [Google Scholar]
  65. Weisberg, J. M., Cordes, J. M., Lundgren, S. C., et al. 1999, ApJS, 121, 171 [NASA ADS] [CrossRef] [Google Scholar]
  66. Weltevrede, P., & Johnston, S. 2008, MNRAS, 391, 1210 [NASA ADS] [CrossRef] [Google Scholar]
  67. Weltevrede, P., Edwards, R. T., & Stappers, B. W. 2006, A&A, 445, 243 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  68. Weltevrede, P., Stappers, B. W., & Edwards, R. T. 2007, A&A, 469, 607 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  69. Wolszczan, A. 1978, A&A, 63, 425 [NASA ADS] [Google Scholar]

Appendix A: Additional Figures

thumbnail Fig. A.1.

Distributions of the normalised Stokes parameters in the Poincaré sphere at pulse longitudes of approximately −0.7° (upper left), −0.3° (upper right), 0.3° (lower left), and 0.7° (lower right). In each quadrant, the two upper panels and the lower-left panel show the projected, 2D distribution of the normalised Stokes parameters in the major planes of Poincaré space, and the 3D distribution is displayed in the lower-right panel.

thumbnail Fig. A.2.

Prolate spheroidal distribution of the normalised Stokes parameters in the Poincaré sphere at pulse longitude ∼−6.8°. The two upper panels and the lower-left panel show the projected, 2D distribution of the normalised Stokes parameters in the major planes of Poincaré space, and the 3D distribution is displayed in the lower-right panel.

thumbnail Fig. A.3.

Four-mode mixing model. Top row: signal from emission region ℳ1; middle row: signal from emission region ℳ2; and bottom row: incoherent superposition of ℳ1 and ℳ2. Upper two rows: from right to left, the right-most panels depict the OPM pairs as linearly polarised states in the x − y plane orthogonal to the line of sight of the antenna; the middle panels depict the OPM pairs in the Poincaré sphere; and the left-most panels depict the modulated time series produced by the OPM pairs over two drift cycles, where the dashed lines represent the positive OPM state, the dash-dot lines represent the negative OPM state, the upper thick lines represent the total intensity, and the lower thick lines represent the superposed polarisation state. Bottom row: the right panel depicts the four modes in the equatorial plane of the Poincaré sphere and the left panel depicts the modulated time series produced by the incoherent superposition the OPM pairs over two drift cycles (colours and line types are the same as for the upper two rows, and the resultant total intensity is shown in red.

All Figures

thumbnail Fig. 1.

Longitude-resolved mean and variability of the polarised emission from PSR B1919+21 at 1414 MHz. Panel a: both the ellipticity angle of the average Stokes parameters (points with error bars) and a histogram of the ellipticity angles of the individual pulses (greyscale intensity). Similarly, panel b: both the position angle of the average Stokes parameters and a histogram of the position angles of the individual pulses. Panel c: integrated profile in total intensity (black, upper solid), linear (red, dashed) and circular polarisation (blue, dot-dashed). These profiles were obtained by averaging together all the available data. In panel d, the black solid line shows the standard deviation of the total intensity, and the red dashed, green dot-dashed, and blue dotted lines display the square roots of the eigenvalues (from largest to smallest, respectively) of the 3 × 3 partition of the covariance matrix that describes the fluctuations of the polarisation vector S = (Q, U, V).

In the text
thumbnail Fig. 2.

Integrated number of single-pulse polarisation states projected onto a two-dimensional plane that is normal to the minor axis of the oblate spheroidal distribution observed at pulse longitude 0°. The greyscale intensity represents the count of instances in each pixel after integrating 1100 single pulses; contours are spaced at intervals of 5 counts. The angle between the Stokes V axis and the normal to the Q′−U′ plane is 7.2°.

In the text
thumbnail Fig. 3.

Total intensity histogram at pulse longitude 0°, where the Stokes parameters are toroidally distributed. This example demonstrates the validity of Stokes parameters normalisation. Stokes I has a similar behaviour throughout all pulse longitudes. The vertical dashed line shows the maximum value of the median off-pulse noise.

In the text
thumbnail Fig. 4.

Distribution of the normalised Stokes parameters in the Poincaré sphere at ϕ 0.°. The two upper panels and the lower-left panel show the projected, 2D distribution of the Stokes parameters in the three planes normal to the basis vectors, and the 3D distribution is displayed in the lower-right panel.

In the text
thumbnail Fig. 5.

Longitude-offset cross-correlations of each Stokes parameter as a function of lag τ and longitude offset Δϕ with respect to the reference pulse longitude ϕ0 ≃ −6.8°. In each panel, a horizontal black line marks zero lag and a vertical black line marks zero longitude offset.

In the text
thumbnail Fig. 6.

As for Fig. 5, but for reference longitude ϕ0 = 0°, in the middle of the region where the toroidal distribution of the polarisation vector in the Poincaré sphere is observed.

In the text
thumbnail Fig. 7.

Longitude-resolved Stokes Q, U, V cross-correlations as a function of lag τ and pulse longitude ϕ. From top to bottom are plotted the correlations between Stokes Q and U, Stokes U and V, and Stokes V and Q. In each panel, a horizontal black line marks zero lag.

In the text
thumbnail Fig. 8.

P3-folded polarisation state of the first 500 pulses in the PSR B1919+21 dataset. From top to bottom are the total intensity; the position angle, θ = tan−1(U′/Q′)/2; the fractional linear polarisation, L ̂ = Q 2 + U 2 / I $ \hat{L}=\sqrt{Q^2+U^2}/I $; the ellipticity angle, ϵ = sin−1(V′/L′)/2; and the fractional circular polarisation, C ̂ = | V | / I $ \hat{C}=|V|/I $; where P3-folded quantities are shown with primed symbols (e.g., Q′ and U′) and the symbols for single-pulse quantities are not primed (e.g., Q and U). That is, L ̂ $ \hat{L} $ and C ̂ $ \hat{C} $ are derived from the single-pulse Stokes parameters before P3 folding and θ and ϵ are derived from the P3-folded Stokes parameters. Furthermore, L ̂ $ \hat{L} $ and C ̂ $ \hat{C} $ are integrated only when the single-pulse total intensity exceeds 4 times the standard deviation of the of noise. In each panel, the modulation period (P3 ≃ 4.2 P1) is plotted twice to better illustrate the continuity of the pattern.

In the text
thumbnail Fig. A.1.

Distributions of the normalised Stokes parameters in the Poincaré sphere at pulse longitudes of approximately −0.7° (upper left), −0.3° (upper right), 0.3° (lower left), and 0.7° (lower right). In each quadrant, the two upper panels and the lower-left panel show the projected, 2D distribution of the normalised Stokes parameters in the major planes of Poincaré space, and the 3D distribution is displayed in the lower-right panel.

In the text
thumbnail Fig. A.2.

Prolate spheroidal distribution of the normalised Stokes parameters in the Poincaré sphere at pulse longitude ∼−6.8°. The two upper panels and the lower-left panel show the projected, 2D distribution of the normalised Stokes parameters in the major planes of Poincaré space, and the 3D distribution is displayed in the lower-right panel.

In the text
thumbnail Fig. A.3.

Four-mode mixing model. Top row: signal from emission region ℳ1; middle row: signal from emission region ℳ2; and bottom row: incoherent superposition of ℳ1 and ℳ2. Upper two rows: from right to left, the right-most panels depict the OPM pairs as linearly polarised states in the x − y plane orthogonal to the line of sight of the antenna; the middle panels depict the OPM pairs in the Poincaré sphere; and the left-most panels depict the modulated time series produced by the OPM pairs over two drift cycles, where the dashed lines represent the positive OPM state, the dash-dot lines represent the negative OPM state, the upper thick lines represent the total intensity, and the lower thick lines represent the superposed polarisation state. Bottom row: the right panel depicts the four modes in the equatorial plane of the Poincaré sphere and the left panel depicts the modulated time series produced by the incoherent superposition the OPM pairs over two drift cycles (colours and line types are the same as for the upper two rows, and the resultant total intensity is shown in red.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.