Next Article in Journal
Correction: Godzieba, D.A.; Radice, D. High-Order Multipole and Binary Love Number Universal Relations. Universe 2021, 7, 368
Next Article in Special Issue
Neutron Stars in the Symmetron Model
Previous Article in Journal
The Effect of Charge, Isospin, and Strangeness in the QCD Phase Diagram Critical End Point
Previous Article in Special Issue
Quasi-Periodic Oscillatory Motion of Particles Orbiting a Distorted, Deformed Compact Object
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Magnetospheric Dichotomy for Pulsars with Extreme Inclinations

1
Gravitational Wave and Cosmology Laboratory, Department of Astronomy, Beijing Normal University, Beijing 100875, China
2
Advanced Institute of Natural Sciences, Beijing Normal University, Zhuhai 519087, China
Universe 2021, 7(12), 455; https://doi.org/10.3390/universe7120455
Submission received: 11 October 2021 / Revised: 17 November 2021 / Accepted: 18 November 2021 / Published: 23 November 2021
(This article belongs to the Special Issue Compact Astrophysical Objects)

Abstract

:
In this work, we expand on a comment by Lyne et al. (2017), that intermittent pulsars tend to congregate near a stripe in the logarithmic period versus period-derivative diagram. Such a stripe represents a small range of polar cap electric potential. Taking into account also the fact (already apparent in their Figure 7, but not explicitly stated there) that high-fraction nulling pulsars also tend to reside within this and an additional stripe, we make the observation that the two stripes further match the “death lines” for double- and single-pole interpulses, associated with nearly orthogonal and aligned rotators, respectively. These extreme inclinations are known to suffer from pair production deficiencies, so we propose to explain intermittency and high-fraction nulling by reinvigorating some older quiescent (no pulsar wind or radio emission) “electrosphere” solutions. Specifically, as the polar potential drops below the two threshold bands (i.e., the two stripes), corresponding to the aligned and orthogonal rotators, their respective magnetospheres transition from being of the active pair-production-sustained-type into becoming the electrospheres, in which charges are only lifted from the star. The borderline cases sitting in the gap outside of the stable regime of either case manifest as high-fraction nullers. Hall evolution of the magnetic field inside orthogonally rotating neutron stars can furthermore drive secular regime changes, resulting in intermittent pulsars.

1. Introduction

Some pulsars show systematic long term variations in their emission behaviors, in the form of intermittency [1,2,3,4], whereby the emission ceases for a significant period of time before turning back on again. Although rather preliminary in terms of statistics due to the small number of such systems so far observed, it has been noted by [4] that these pulsars appear to congregate along a straight stripe in the period versus period-derivative or log 10 P log 10 P ˙ diagram (see Figure 1). We also point out that the a grouping of high-fraction nullers (we distinguish them from intermittent pulsars by their much shorter nulling periods—comparable in order of magnitudes to pulse intervals; see also Section 4 for their difference with low-fraction nullers) around that stripe is also quite discernible, with furthermore, the remaining high-fraction nullers seen (with significant scatter) to congregate near another similar stripe. It is noted by [4] that such lines/stripes are likely significant because they correspond to constant potential drops across polar caps (abbreviated as P below), in the rudimentary vacuum dipole model of the pulsar magnetosphere (summarized in Appendix A.1), which also serves as a baseline parameter for constructing more sophisticated models (summarized in Appendix A.2 and Appendix A.3).
Zooming further out to encompass all pulsars, it is particularly noteworthy that not all pulsars near the proposed lines share such variability, so the variable population must correspond to special configurations possessing particular but not abnormal (otherwise the pulsars likely will not behave normally in their “on” states) parameters. Recalling that those lines of constant P are meaningful only when comparing pulsars of similar inclinations (angle between the rotational and magnetic axes, see Appendix A.1 below), it is then natural to suspect that these pulsars are nearly (not precisely, thus some scattering about the lines are to be expected) orthogonal or aligned rotators, respectively for the two lines, since such special inclinations are known to be problematic for secondary pair plasma production1. The corresponding magnetospheres could then plausibly begin to fall into comatose states when the electromotive force across the polar caps drop below critical values (signified by the two stripes), for which radio emission and pulsar wind start to become suppressed.
In this paper, we investigate this possible association of some variable pulsar populations with extreme inclinations, beginning with an assessment of the observational evidences in Section 2, followed by theoretical speculations on the underlying physics in Section 3, before concluding in Section 4.

2. Statistics

2.1. Clustering

Figure 1a2,b are visually suggestive as to the clustering of intermittent pulsars and nullers in the log 10 P log 10 P ˙ diagram, but eye-balling is not adequate, since, e.g., the underlying overall pulsar population density (plotted in Figure 23) may already be abnormally high in those P neighborhoods, such that similar clustering patterns would arise if we just randomly drawing some samples out of the general population. We need a more quantitative gauge of how likely this feature arises purely by chance, especially since the data points are scarce and scattering is not negligible (recall that, to complicate matter, some scattering is unavoidable because the relevant pulsars inclination angles will not all be precisely 0 or π / 2 ). For example, we could compute the silhouette values for the members of these populations against the generic pulsar population displayed in Figure 2. Let a i denote the average distance in P (computed according to Equation (A6) below) of a sample point i amongst, e.g., nullers, to all other members of the same collection, and b i its average distance to all pulsars, then the silhouette
s i = b i a i max ( a i , b i )
satisfies 1 s i 1 , and a value close to 1 indicates vary tight clustering, i.e., a i b i , while a negative value means the clustering is phantom, in that the sample point is more closely associated with the general population than with other points in its assigned collection. The results are listed in Table 1, and indeed all show large positive values. Note that we have also listed the silhouette values for stripe O nullers against the A grouping (i.e., b i are average distances to nullers on A rather than the general population), and vice versa. The fact all these values are positive indicates that we have not misassigned affiliation of points (reassign any nuller will result in negative values).
This simple (thus robust) silhouette would already constitute a quantitative measure, but we must also correct for the influence of our small sample sizes. To achieve “normalization” (and also to provide a more intuitive interpretation of what the silhouette values imply), we can further compare the mean silhouette value s ¯ for a observational sample against its counterparts computed for the same-sized mock up samples drawn randomly via Monte Carlo simulations from the general pulsar population (i.e., 8 general pulsars are randomly picked to mock up a sample whose silhouette value can be compared with that of nullers in stripe O ; similarly 22 are plucked for comparison with nullers in stripe A and 5 for the intermittents). The simulations as such provide a distribution of s ¯ appropriate for this particular sample size, and where the physical ensemble lands within this distribution (specifically, what ratio of Monte Carlo outcomes end up with greater or equal s ¯ than the observational sample) tells us the likelihood of the stripe-like-grouping arising purely by chance—i.e., the probability that the intermittent pulsars and high-fraction nullers have no special preference for P , as compared to the general population.
The simulations are trivial to carry out. One merely needs to draw random numbers from a uniform distribution over the index set of all pulsars, and then the pulsars whose numbers get called join the simulated ensemble. This is analytically equivalent (but numerically more efficient and accurate) to the more spelled-out procedure of extracting the joint distribution of pulsars over the log 10 P log 10 P ˙ plane, collapsing it into a marginal distribution over P by integrating out the directions of constant P on the log 10 P log 10 P ˙ plane, and then drawing random samples by inverting the resulting cumulative distribution function over a uniform distribution over the interval [ 0 , 1 ] . The results of our simulation are displayed in Table 2, which show that very small percentages of samples have greater or equal s ¯ than the observed populations, thus demonstrating that the observational clusterings are unlikely chance occurrences. Converting to the familiar Gaussian sigma terminology, the odds of the three clusters arising by chance are ≈2.7 σ for intermittent pulsars, ≈4.3 σ for nullers near O and >5.7 σ for nullers near A (up to the maximum 10 8 samples that our computing resources permit, there has been no instance of equally tight or tighter clustering arising by chance, which sets an upper limit of 10 8 , and thus an inequality). Note the progression of increasing number of sigmas is due to the differences in population sizes, matching the intuition that it should be more unlikely for larger populations to cluster by chance. This is what we meant by “normalization by sample size” earlier.
It should also be noted that we have split the nullers into two subpopulations that are examined separately, an action that has thus far been based on our theoretical understanding that there exists of two classes of problematic magnetospheres—the aligned and orthogonal rotators. However, if one prefers statistics extracted from unadulterated purely observational data, without introducing theoretical priors, then they perhaps should consider all nullers as a single collection, and count all Monte Carlo samples that can be subdivided into a small number of very tightly clustered subgroups as being equally good or superior. We expect this alternative approach to yield slightly less extreme number of sigmas, depending on how many sub-clusters we allow and how we average their silhouettes. This fact that there are many tunable knobs is a significant deficiency though. In particular, we do not see how theoretical priors can be prevented from sneaking back in, through whichever way we penalize larger numbers of sub-clusters (at the very least we need to decide on how many is too many, inevitably evoking a theoretical expectation of how many there should be). Failing this prevention would then defeat the purpose of introducing this more opaque and technically more complicated approach in the first place, so we do not adopt it here. Furthermore, as we now discuss in the next section, there in fact exists observational evidence for the two stripes being associated with aligned and orthogonal rotators, respectively, thus our division of the high-fraction nuller population in two is not solely theoretical.

2.2. Inclination Angles

For a small number of pulsars, relatively precise measurements of the inclination angle α , by fittings to polarization position angles or PPA (as per [31]), is available. Near O , ref. [32] showed that the intermittent J1933+2421 is indeed a nearly orthogonal rotator. In addition, measurements for J1841-0500 had been done by [2], giving an impact parameter of β 3 . Although a measurement on α is not explicitly given in that work, we can infer from the very steepest slope of the PPA curve (the steepest part equals sin α / sin β ) that sin α should be quite close to saturating out at sin π / 2 . For nullers near line A on the other hand, explicit position angle studies by [33] (≲2 for J1944+1755) and [34] (for J1107-5907) have also confirmed that the inclination angle α is indeed very small. The high-fraction ( 70 % ) nuller J0826-3417 likewise has a very small α since the emission covers the entire pulse longitude (we are always inside the emission cone) [11].
Another supporting evidence for associating O and A stripes with inclinations is that since they are essentially the “death lines” for orthogonal and aligned rotators within our proposal, pulsars (not necessarily nulling or intermittent) that exhibit robust characteristics confirming their extremal inclinations should exist only to the upper-left regions of these stripes. One particular feature of orthogonal and aligned rotators is that they often (but not always, depending on the detection threshold among others, but for our purpose we only need this effect to be sufficient and not necessary) exhibit interpulses, of the double and single-pole variety, respectively [35]. With orthogonal rotators, it is because the emission cones from both poles could sweep across Earth, so interpulses sit roughly half-way in-between “regular” pulses. With nearly aligned rotators, the “regular pulse”-interpulse pair would be due to the Earth traversing the two walls of a hollow conic emission beam [36], usually with a bridge connecting the two (see, e.g., [11] for J0826-3417). This is essentially an extreme version of the double-poled pulse profile, whereby the two peaks are extra-widely separated due to the small α , to the extent that they are (mis)identified as separate pulses. Alternatively, the single-pole interpulses can be due to the line of sight being always within the beam (possible for nearly aligned rotators), and threading through two nested cones or a core-conal duo [37,38]. With either interpretation, single-pole interpulses flag nearly-aligned rotators. In Figure 34, we plot the log 10 P log 10 P ˙ locations of the pulsars exhibiting double and single-pole interpulses (as compiled by [35]), and they indeed appear to reside mostly to the upper left of the relevant stripes.
Once again, we could assess this observation more quantitatively, by defining a quantity that we term the “penumbra”, that measures the amount of spillover beyond the relevant death lines:
p D = 1 n i max P D P i , 0 ,
where i indexes the single or double-pole population, n is their population size, D = O or A , and P D is associated with the center of the relevant stripe (average of the P values for the intermittent and/or nulling pulsars clustered around the stripe, see caption to Figure 1). For any prescribed death line, the penumbra vanishes when all data points lie to the upper-left of it. It does not however, by itself, tells us whether the death line is appropriate, in the sense that one could obviously place it very far to the bottom-right of the log 10 P log 10 P ˙ diagram, so that the interpulse (and indeed all) pulsars lie to its top-left trivially. Such a situation conjures little statistical significance.
To handle this issue and to provide better intuition, we can once again evoke the Monte Carlo simulation, and record what percentage of simulated sample populations, of the same size as the single- or double-poled pulsars, possess smaller or equal penumbras than the actual observed population. Or in other words, to what likelihood a population of this size, residing mostly to the top-left of the death lines, would arise purely by chance. This gauge also penalizes over-generous death lines, since pushing them too far towards the bottom-right corner would lead to greater portions of simulated samples possessing equal or smaller penumbras, reducing (as should be) the statistical significance readings. This in fact causes trouble with the A stripe, which is unfortunately very close to the overall pulsar death line (see Figure 2), so the probability of the single-pole population residing mostly to its top-left purely by chance is quite high (see the “Single-Poled” rows in Table 3), and we cannot assert this stripe’s role as a death line for the aligned rotators based on interpulse population statistics. Beyond the available α measurements for the few nullers discussed at the beginning of this section, we will have to base our confidence in this hypothesis mostly on its agreement with theoretical predictions. Note though, this feature of A sitting far to the bottom right does not translate into aligned rotators being easy to turn on, and is instead simply the result of small α causing large Δ Φ vac values to map into small P (i.e., for similar stripe placements, much larger potential drops are demanded for aligned rotators than for inclined rotators, to which the overall death line refer), see Equation (A7) below. Fortunately, the O stripe sits much further up the P spectrum, so meaningful statistics (≈ 4.4 σ or ≈5.6 σ depending on ambiguous population grouping, see Table 3) can be obtained for it, validating O as a likely death line for orthogonal rotators.
So far, our conclusions have been based on observational statistics, but theoretical arguments can be brought in to further enhance our confidence, that O is associated with orthogonal rotators. Specifically, with intermittent pulsars, it is quite natural to expect that the “on” and “off” states correspond to two modes of the pulsar magnetosphere, one with active pulsar winds and another with only vacuum-like dipole radiation. The simplest baseline estimate for the ratio of spin-down rates between them is [53]
3 2 1 + sin 2 α sin 2 α ,
which while unfortunately always overshoots the observationally measured values that sit in the range of 1.5–2.5 [54,55]. Nevertheless, this expression attains its minimal value (closest to the observed values) of 3 at α = π / 2 (this ratio diverges for aligned or anti-aligned rotators in contrast), meaning it would be much easier to tweak the model (specifically in our case, the “off” state are electrospheres rather than actual vacuum, while the “on” state deficiencies of orthogonal rotators that we discuss in Section 3.3 have not been included in Equation (3)) to match observations if intermittent pulsars are orthogonal rotators.
Lastly, for completeness, we note that due to the scarcity of precise inclination angle estimates based on position angle sweeps, indicative estimates have frequently been produced in literature using pulse width measurements, with [56]
W ( α , ρ , β ) = 4 arcsin sin ρ + β 2 sin ρ β 2 sin α sin ( α + β ) ,
where the emission beam width ρ (recall that β is the impact parameter) is commensurate with the emission level (e.g., 50 % or 10 % ) relevant for the pulse width W . This method suffers from great uncertainty though. The sensitive dependence on β in particular, causes considerable variation in the estimate of α from W . Most crudely but conveniently, one may take β = 0 , and obtain the commonly evoked approximation W 2 ρ / cos α [35]. A more complex empirical rule have also been proposed [57] based on the existence of a Lower Boundary Line, where
α arcsin 2 . 45 P 0.5 W core ,
with W core being pulse width at 50 % strength of the core component in degrees. Even Equation (5) is too broad-stroked though (only useful when averaging over a large population of pulsars), since, e.g., for J1944+1755 the position angle analysis gives α 2 [33], while Equation (5) gives 14 (assuming the whole width is core). To avoid misleading conclusions then, one must consider a range of α values achieved by allowing β to vary. To this end, it is more convenient to replace β with ζ = α + β , i.e., the angle between the spin axis and the observer, yielding [58]
W ( α , ρ , ζ ) = 2 arccos cos ρ cos α cos ζ sin α sin ζ ,
which is relevant only when | α ρ | ζ α + ρ , ensuring that the pulsar is seen on Earth as pulsed emissions [59]. We plot in Figure 4 the W values as functions of α and ζ for a couple of example choices ρ = 10 and 20 . The contours of constant W shows that α can span vast ranges, so α derived from such considerations tend not to be very constraining. In particular, the contour plot of Figure 4b shows that α = 0 is always an admissible limiting solution in the ζ ρ limit, while maximum α value occurs close to the diagonal α = ζ line. Therefore, for any width W 0 corresponding to a particular pulsar, the value α diag as given by W ( α diag , ρ , α diag ) = W 0 provides an indication of the maximum α achievable for that pulsar.
Note also that in addition to β or equivalently ζ , ρ is also variable. One can in principal take the empirical ρ 10 % 4 . 9 P 1 / 2 to 6 . 3 P 1 / 2 [60,61,62] for the general pulsar population, but since we are examining geometrically special aligned and orthogonal rotators here, their beam widths may not exactly follow this rule. Furthermore, for many of the high-fraction nullers depicted in Figure 1, it turns out that there is not a solution for α diag , since the corresponding contour does not intersect the diagonal line in Figure 4b for the given ρ (i.e., the contours corresponding to those crowding near the boundary of the rectangular region; equivalently the ones at lowest W values in Figure 4a), meaning that there is no constraints on α at all. In summary, constraints on the inclination angle from pulse width is not particularly informative in our context, and we shall not pursue this avenue further.

3. Proposed Explanation

3.1. Magnetospheric Dichotomy

Finer details within the log 10 P log 10 P ˙ diagram can clue us in on interesting features of pulsar magnetospheres, so it is worthwhile coming up with a theoretical model that explains the stripes. In other words, it is interesting to see what surprising structural properties of the magnetosphere they may be hinting at. We propose that charge-abundant active magnetospheres and charge-starved quiescent electrospheres (summarized in Appendix A.2 and Appendix A.3) describe the live and dead pulsars that are separated by the O and A stripes. These stripes thus serve as the death lines for (nearly) orthogonal and aligned rotators, respectively.
The discriminating factor determining which of these magnetospheric states is realized is whether pair cascade occurs to provide an abundance of charged particles. Typically, cascading behavior tends to emerge only when a threshold is crossed in the parameter space. It is quite plausible that aligned and orthogonal rotators, that have been known to suffer from different deficiencies with regard to pair production, will cross below said threshold once the electromotive force across the polar caps drops to their respective critical levels, as signified by the stripes A or O . Indeed, numerical work by, e.g., [63] (for 1-D) and [64,65] (for aligned rotators) had indicated that the pulsar magnetospheres would revert to the charge-separated electrospheres when the pair production threshold exceeds the full polar cap potential.
In other words, for pulsars with extreme inclinations and residing to the lower right of the relevant stripe, their e ± cascade shuts down, leaving only those charges lifted directly out of the stars available to short out the parallel (to the magnetic field) electric field E and achieve force-freeness5 in the charge regions. This in turn reduces plasma multiplicity, driving the magnetospheres into the quiescent charge-separated electrospheres that those earlier studies had seen6. This subsequently shuts down pulsar winds, since there are no current flows, as well as radio emission, which is assumed to emerge from near the separatrix current sheets (missing in electrospheres) in some models (see, e.g., [67]), and otherwise modeled in neutral and not charge-separated (as in the electrosphere case) plasma in the rest (see, e.g., [68]).
In essence then, we propose that both sides of the aforementioned historical debate on the magnetospheric structure are in fact correct, with both types of magnetospheres realized in nature, reigning their respective realms (we will slight abuse terminology in the interest of brevity, and use “pulsars” to refer to both cases, even though the electrospheres give off no radio pulses). This paper is of course not the first to touch on such a bimodal possibility, e.g., [69,70] motivated the study of electrospheres by noting that they would be useful if for some reason pair cascade ceases. The novelty of our work is mainly in the proposed role played by extreme inclinations.
There are also the more marginal cases where the nearly extreme inclination pulsars reside close to the relevant death lines, and hop back-and-forth between active and quiescent states. The magnetospheres for these pulsars evolve dynamically and thus should differ from either of the stationary cases described in the appendix. Absent full 3-D first-principle (including pair production microphysics) particle-in-cell numerical simulations, an accurate picture for such dynamical states is difficult to obtain (although already quite sophisticated simulations as those done by [71,72] may have already captured many essential qualitative aspects; note that with the first reference, although their modeling is aimed at intermittent pulsars, the description could plausibly also be applied to the nulling cases), but we know such states can reside within the gap between the stability thresholds, P O / A a and P O / A e for the active and electrospheric magnetospheres, respectively, where both those stationary configurations become unstable. Specifically, if
P O / A a > P O / A e ,
and if the actual potential P O / A for a nearly aligned or orthogonal rotator satisfies
P O / A a > P O / A > P O / A e ,
then it cannot be steadily “off”, nor can it remain consistently active. It would be stuck in a cycle of “somewhat on” and “somewhat off”, manifesting as a high-fraction nuller.

3.2. Short Circuiting the Electrospheres

A physically significant feature of the electrospheres is that E is shorted out only within the charge clouds that are confined to regions close to the star (mostly inside the trap surfaces). All of the regions outside of the clouds are vacuum “gaps” whereby E is not shielded, and are thus possible sites of pair production, as has been highlighted by, e.g., [70]. Such additional sources of charges are not included in the computations of the electrospheres, and could thus in principle invalidate them. However, the electrospheres should be stable against small amounts of charge injections [73], given that the direction of the electric field in the gaps is such that it transports the particles to the clouds of the same charges, thus maintain charge separation.
However, if these “gap” regions host strong enough E to drive runaway pair cascade, then we will witness copious flows of particles along certain magnetic field lines, esp. along the would-be separatrix current sheets [72], turning on global current circulation within the magnetosphere [74], with current density reaching such intensity that the magnetic field is no longer dipolar, and a transition into the active high-multiplicity magnetosphere occurs. To narrow down the scope for a look into finer details, we concentrate on the single photon conversion process, since the thermal X-ray photons required for the two-photon interactions (serving as the abundant targets for the less numerous curvature γ -ray photons to hit) would unlikely be sufficiently populous, because the stellar temperature of around 10 6 K [75] would be too low [70] (the two-photon process can however be important for more powerful pulsars when current sheets are present, see next section for more discussions).
For aligned and orthogonal rotators, the site for pair cascading by the single photon absorption process would likely be substantially divergent. With orthogonal rotators, if Δ Φ vac is large enough, pair cascade could occur far away beyond the trap surfaces, since the rotation induced electric field as given by (A12) and (A13) contains slowly declining 1 / r 2 pieces when α = π / 2 . In contrast, the leading contribution to the electric field is 1 / r 4 (quadrupolar) for aligned rotators. Indeed, computation by [70] shows that this process is likely operative only very close to the star, within around tens of stellar radii. One possible site for aligned rotator pair cascade is then in the crevice between the dome and the disk (close to the orange thread between the red and blue clouds of Figure A3a). Such a case is indeed seen in the axisymmetric simulations by [72], where an intermediate Δ Φ vac corresponds to a cyclic behavior of the magnetosphere, switching between activity and comatose, bearing obvious resemblance to the nulling behavior. The electronics analogy of such a magnetosphere is perhaps an LC circuit, with a capacitor in the form of a charge cloud residing at the Y point, as opposed to an open circuit as electrospheres or properly closed circuit as the active case.
Although not yet numerically investigated, we can expect that if we further relax axisymmetry, another cascade site would become available. The reason being that the diocotron/slipping-stream instability (see Appendix A.3) would become activated in this case, causing the equatorial disk to develop into high charge density vortices in the azimuthal direction, which subsequently merge and either collapse to hug the stellar surface closely, turning the entire large radii regions vacuum, or, when particle injection is introduced, form complex radially-extended patterns containing prominent vacuum gaps (see Figures 9 and 11 of [76]). These gaps cannot be closed up by the pairs produced within (they are not closed even by artificially injected single-signed charges), since the walls bounding a gap are all made up of the same charges. Specifically if, say, electrons produced by pair production, re-enforce and thus draw one wall in, the positrons of the pairs would inevitably annihilate and push back the border on the other side. These long-lasting gaps, residing within the disk where the magnetic field is strong and single photon conversion is efficient, would then keep injecting charges into the magnetosphere, turning the pulsar on (see, e.g., [77,78] for active magnetospheres sustained by particle injected from close to the star). Furthermore, the diocotron instability grows on timescales comparable to pulsar rotation periods [76,79,80], implying that the transition from “off” to “on” states can be accomplished within timescales as short as a few pulse intervals, which would indeed be necessary if it is to explain the nulling behavior. Furthermore, the cyclic behavior resulting from pair cascading only in the “crevice” is also of periods comparable to rotation period [72], thus either site (“crevice” or “disk opening”), or indeed both, can be relevant for nulling.
For a more quantitative assessment, we note that the approximate criteria for single photon cascading is given by [81] (see their Equation (1.4) and discussion below Equation (1.5), also their Figure 12) as
χ ϵ γ 2 m e c 2 B B cr 10 1 ,
where ϵ γ is the characteristic energy associated with the physical process under consideration (photon energy in our context), B is the magnetic field strength and
B cr m e 2 c 3 e = 4.414 × 10 13 G .
The electron rest energy on the other hand is
m e c 2 = 0.51 MeV = 8.2 × 10 7 erg .
When applied to magnetospheres, we would want nearly full conversion of photons as they traverse the gap. Let s be the full conversion shielding distance, then Equation (3.6) of [81] yield the requirements (these expressions are derived under the assumption B B cr , which is reasonably valid for typical pulsar magnetospheres of O ( 10 12 ) G)
B s 4.6 , ϵ γ m e c 2 6 × 10 7 s ,
where B is in units of MG, and s is in cm. For an active gap with characteristic dimension g , we would need s g , we have thus the conditions
C 1 e : g 4.6 B , C 2 e : E m e c 2 6 × 10 7 ,
where we have used ϵ γ E g (we are considering the minimal requirement that there exists sufficiently energetic photons, these photons may not be the most populous among those produced by synchrotron radiation), with E being the electric field strength in the gap with a unit of statV/cm. The condition C 1 e is quite simple to satisfy since neutron star vicinity has typical length scales of kilometers, so g 10 5 cm, thus B∼100 G would already be sufficient, meaning that C 1 e is satisfied for regular pulsars all the way out to the light cylinder. The second condition C 2 e is E∼50 statV/cm, which is a little more stringent. From vacuum approximations for which we have analytical estimates (A11)–(A13), we obtain that
E align B 0 Ω 2 c R 5 r 4 , E ortho B 0 Ω 2 c R 3 r 2 .
Substituting in the typical parameters B 0 10 12 G, Ω ∼1, and R∼6 km, we have that the maximum distances at which we could still have sufficient E strength are at around r = 127 km and r = 2700 km for the aligned and orthogonal rotators, respectively.
These numbers are in agreement with our naive expectation that the electrosphere-destroying pair cascades can occur far outside of the trap surfaces for the orthogonal rotators, but can only occur quite close to the star for the aligned rotator. Although we will not be able to predict the exact location of the transition lines due to the limitations imposed by the various approximations we have to adopt (chiefly that the field strength expressions are not those of actual electrospheres), the fact that these numbers, computed with realistic pulsar parameters, are greater than neutron star radius and smaller than the light cylinder radius, suggests that the electrospheres should not be uniformly stable nor unstable for the entire regular pulsar population (if both numbers are smaller than stellar radius, we will not have any place with a sufficiently strong E to break down the electrospheres; and if both numbers are larger than typical light cylinder radius, then there will always be gaps powerful enough to destroy electrospheres, so this quiescent configuration can never exist in reality)—it is indeed plausible for there to be two subpopulations possessing very different magnetospheric structures.

3.3. Replenishing an Active Magnetosphere

The quiescent charge-separated electrospheres has a muted discharge with all the plasma supplied by the star, which is sufficient since there is no pulsar wind and thus no voracious demand for replenishment. The high multiplicity active magnetosphere on the other hand, must be underpinned by an active pair cascade inside gaps where the electric field parallel to the magnetic field is not shorted out. These gaps are formed afresh when the magnetosphere turns active and are not necessarily associated with the gaps in the electrospheres. Once such gaps open up, the same single-photon conversion consideration of the previous section applies, but there is now the prerequisite that such gaps would not be quickly shorted out by the surrounding plasma (as opposed to the electrosphere case whereby stable gaps are known to always be present). We turn now to this additional condition, on the stable existence of the necessary gaps that continuously supply plasma to maintain a charge neutral force-free magnetosphere—if the gaps cannot persist, the unbalanced loss of charged particles through pulsar winds will cause the magnetosphere to become lower multiplicity, evolving towards the charge separated electrospheres.
Previous studies (see, e.g., [63,64,65,82]) have found that the main criteria for efficient pair production in the polar cap is that C a : the local current density j along the magnetic field lines should exceed the Goldreich-Julian or GJ value
j GJ = c ρ GJ = 2 ϵ 0 c Ω · B 1 Ω × r c 2 ,
so that the particle number density as inferred from j still exceeds the GJ density value7, even if the charges are moving with a maxed-out speed vc (thus the charge density requirement is at a minimum). Such over-charging prevents the proper shielding of the E field, instead resulting in an(other) E 0 gap region that keeps accelerating charged particles while they move away from the star, and the photons emitted via curvature radiation or energized via inverse-Compton scattering can trigger pair cascade via single photon conversion [9]. The potential difference that develops in such a gap is approximated by [83]
Δ Φ 4 π R ζ 1 c ρ GJ c h 2 ,
topping out at Δ Φ vac , where h < r pc is the height of the gap ( r pc being the radius of the polar cap), and R is the Ramp function.
It is worth noting that many studies (e.g., [84]) would assume force-freeness a priori for mathematical tractability (i.e., within the formalism adopted, already assuming that the charge density is at the GJ value and E is shielded), then condition C a is equivalently stated for them as there being regions where the 4-D current J a becomes spacelike J a J a > 0 (because the temporal component J 0 = c ρ , equaling c ρ GJ by the force-free assumption, is too small as compared to the spatial component j , exactly the same situation that C a describes), flagging a consistency problem (charges need to move superluminally to produce such a current, which is impossible) indicating that the force-free assumption breaks down, and we would instead witness a gap there. Worth noting also is the fact that we have thus far described only a single-charge-species simplification serving as a time-averaged toy model. The full picture is more complex and dynamic, with the charged particle number density possibly many times that of the GJ value (but since both charges are present, average charge density can be much lower), and with pair production likely occurring in bursts (counter-streaming opposite charges can momentarily keep j high while retaining charge density at ρ GJ in localized regions, but cannot maintain such a balance everywhere in a stationary manner, see [9,63,65]; so spacelike currents still flag the need for pair production, but the temporal and spatial features could be much richer).
There is furthermore an alternative site for pair production relevant for young and fast spinning pulsars [64,85], in the current sheets where ζ j / j GJ < 0 [9], so that the charges have the wrong sign for shielding the electric field. The accelerated charges can emit synchrotron photons (e.g., γ -ray signals from pulsars can be produced this way [86]), and some of which can then produce pair cascade via γ γ interactions. The demands on the magnetic field strength is relaxed as compared to the single photon conversion, since the magnetic field does not participate directly in this process, which can then occur in the current sheets further from the star. However, hints are that the polar caps are likely the more important pair production sites when the pair-injection rate is low [78] (meaning lower overall pair flux is required if injected at stellar surface only), which is the relevant regime for our present considerations. Indeed, computations based on single photon magnetic conversion predicts a pair creation death valley at Δ Φ vac 10 12 V [87], which matches the observed pulsar radio emission turn-off region [9,88].
The determination of j is nontrivial however, since if we see the magnetosphere as a circuit, then some resistivity comes from the crossing of the twisted magnetic field lines by the particles on the outside [74], where the interaction with interstellar medium is complicated. Nevertheless, it appears reasonable to expect that:
  • The aligned (and anti-aligned) rotators are not effective accelerators, because j GJ is quite large when Ω and B are aligned, making it more difficult for j to exceed it. This leads to the likelihood that the previously derived P A P A e reduces to a necessary but not sufficient condition, namely that while it being satisfied is sufficient for pair cascade if there is a gap, the gap will not in fact be available at this P A yet. A more restrictive P P A a fulfilling Equation (7) is needed instead.
    Evidences for this behavior exist in literature, e.g., analytical models predict ζ < 1 [89], and numerical simulations such as those done by [64,82] confirmed that for nearly aligned and anti-aligned rotators, charge-separated electrosphere solutions are obtained (for the non-juvenile pulsars). Such a complete shut-down for all low inclination pulsars can however be remedied by General Relativistic frame-dragging effects [65], which reduce the effective angular velocity of the star and subsequently j GJ . This fact would be consistent with the A stripe being not far from the overall pulsar death line.
  • With orthogonal rotators, the story is different. Because Ω and B are nearly orthogonal, j GJ nearly vanishes and condition C a is easily satisfied. However, this also translates into the fact that we are working with much smaller numbers and the resulting accelerating electric field may not be sufficiently large to drive pair cascades. Indeed, the indicative vacuum electric field as given by Equation (A11) with α = π / 2 and θ = 0 gives a vanishing time-averaged E r in the polar regions. Alternatively, substituting Equation (15) into Equation (16), we see that since Ω · B 0 for an orthogonal rotator, the electric field inside the gap would be E Δ Φ / g 0 . In other words, ζ > 1 only ensures that the electric field will not be shielded, it does not guarantee that the resulting unshielded field is strong enough to accelerate charges to energies sufficient for pair production, as demanded by Equation (13)8.
    As a result, the existence of the polar gaps is not restrictive and the opening up of the band (where high fractional nullers are found) as required by Equation (7) is not due to the extra requirement C a . Instead, the fact that P O a exceeds P O e is simply due to the polar gaps of active magnetospheres residing close to the magnetic axis, where the rotation-induced electric field is suppressed, while the gaps of the electrosphere occupy other latitudes where the electric field is much larger.

3.4. Intermittency

Beyond the magnetospheric changes that we have considered so far, there is a further complication relating to evolutions within the neutron star. Namely, that the magnetic field configurations within the star can contain higher order multipoles, which will affect the field structure immediately abutting the star (but not further away, since the higher order multipoles drop off much more quickly with radial distance), and significantly boost pair cascade. It does so by two means: ( e 1 ) the polar cap shifts location, which turns the effective inclination less extreme; ( e 2 ) the curvature of the magnetic field lines are enhanced so curvature radiation becomes more efficient.
If this happens, an otherwise borderline “off” extreme-inclination pulsar may be turned “on”. For moderate-inclination pulsars, similar internal magnetic field dynamics would also be present, but with less dramatic impacts—the effective inclination is already non-extremal, so altering it a little will not impart significant qualitative changes to pair cascading. The two extreme-inclination cases also differ significantly, as the following simple argument already shows: referring to Figure A1 for a crude dipolar illustration, we see that with an orthogonal rotator (second row in the figure), the Lorentz force ( e / c ) v e × B due to the star’s rotation, as experienced by those highly mobile electrons within the neutron star near the dipolar cap region, is pointed along the z axis (with the direction being consistent for all electrons near a polar cap), meaning that the electrons and thus currents sustaining the magnetic field will be pushed around as compared to the lowest order dipole approximation (thus offsetting from the positively charged immobile ions, producing an electric field, to achieve force balance as per usual Hall effect), so the centers of the polar caps will likely shift from being in the equatorial plane to somewhere with a nonvanishing z value; in contrast, this force is axisymmetric against the z axis for an aligned rotator, so the polar caps’ centers will remain on that axis, and the polar caps will thus only change in size, but will not shift in location. In other words, effect ( e 1 ) is suffered more acutely by orthogonal rotators.
Furthermore, if there are multiple equilibrium states that are similar in energy, but some with greater higher order multipoles giving rise to abundant pair production and thus active magnetospheres, and some less, leading to subdued pair production and thus quiescent magnetospheres, then we would observe the internal field and current configuration switching between these states on Hall evolution timescales9, driving the magnetosphere to hop between active and comatose states. In other words, we have intermittency for orthogonal rotators (more pedestrian radio emission mode-changing may be due to the same physics, but for moderately oblique and aligned rotators), which is a different thing from nulling, since the pulsar jumps between truly (quasi-stationary) “on” and “off” states, rather than being stuck in the middle band between P O a and P O e (although the “on” state for some cases could obviously also land in the middle zone).
Turning to more details, we have the Hall evolution equation
B t = c 4 π e n e × × B × B ,
which is simply the Faraday induction law supplemented by the Hall-drifting condition on the particles (i.e., ignoring the non-Hall terms in the generalized Ohms law for plasmas, such as pressure gradient and the small resistivity [93]; see [94] for the physical conditions that the negligence of those terms corresponds to)
E = v e × B c ,
as well as the neglecting of the displacement current (thanks to the high conductivity in the star) in the Ampère’s law so v e j × B , where v e is the electron velocity and j is the current density. Equation (17) is nonlinear so we will not be able to solve it to provide quantitative descriptions. Nevertheless, it has the following qualitative/semiquantitative properties that argue for its relevance to the intermittency phenomenon:
  • Because the Hall drift conserves energy, there is no reason (on the short timescales we are interested in, when Ohmic dissipation is negligible) for the internal magnetic structure to be stuck in some equilibrium state. It would instead be able to wonder around along equi-energy surfaces in the space of such structures (i.e., the shifted polar caps will likely move about).
    Specifically, Equation (17) is hyperbolic, with component form resembling the well-known Burgers’ equation (after introducing simplifying symmetries) [92]. As such, it generically exhibit wave-like periodic behavior—manifesting as the low frequency helicons/whistlers [95,96] in the linearized limit. Such oscillatory behavior is indeed seen in the numerical simulations of, e.g., [97] (although the star is simplified into a uniform one rather than having a crust-core division for this study, and the setup is axisymmetric).
    Note that while the underlying internal magnetic field structure evolves continuously, the external pair cascade should turn on and off in a more clear-cut binary manner, because the relevant microphysics is characterized by sharper threshold conditions, such as C 2 e .
  • Note that when deriving Equation (17), we have neglected inertial forces such as the Coriolis force that would be seen in the co-rotating frame, since they are negligible in strength (contributing to the force balance at equivalent to around only 200 G of magnetic field strength). This amounts to the zero electron mass limit (the strong magnetic field means that the cyclotron radius is small and so the microscopic cyclotron motion of the electrons will not complicate the dynamics), which removes the Trivelpiece–Gould mode and disables the associated potential enhancement on plasma production at the conductor surface. As a consequence, the motion of the electrons is determined completely by the electromagnetic field, and so just as with the force-free case, we have that the particle motions contributing to a nonlinear modification to an otherwise vacuum electromagnetism.
    Unlike most discussions in helicon literature (see introductory texts such as [96]), our physical setup is not one where we can approximate the magnetic field as a small perturbation on top of a constant background. In other words, we cannot linearize the problem, and the periodicity in our nonlinear setting cannot be expected to be very regular (indeed, the rigorous mathematical definition of “waves” in the nonlinear context is subtle). This nonlinearity would generically cause the evolution of the internal magnetic field to appear rather stochastic.
    Note for theoretical studies such as [97] though, one has to impose artificial symmetries for tractability, thus will see more organized periodic behavior. In other words, one should expect observed intermittency to exhibit less clean periodicity than seen naively in theoretical predictions. Observationally, the intermittent pulsars are indeed seen to be quasi-periodic, see, e.g., [98].
  • The Hall drifting timescale is given by [99] as (robust and not sensitive to the state of matter in the neutron star; also dimensionally consistent with Equation (18))
    t Hall 5 × 10 8 L km 2 B 12 ρ ρ nucl yr ,
    where the characteristic magnetic field strength B 12 = B / 10 12 G O ( 1 ) for regular pulsars, while L km = L / 10 5 cm is the length scale at which the magnetic field density changes, which is again ∼ O ( 10 1 ) to O ( 1 ) , being limited from above by the thickness of the neutron star crust, and below by the size of the polar cap (we are mostly interested in the lowest few higher order multipoles that would change the overall structure of the polar caps if ( e 1 ) is the most important effect; if ( e 2 ) is more important, higher order multipoles may also become relevant). The constant ρ nucl = 2.8 × 10 14 g / cm 3 is the nuclear density, and ρ is the density of the neutron star in the upper crustal region of interest, which is around 8.0 × 10 6 g / cm 3 (specifically, the density of 56 Fe in the ground state of cold dense matter10 [101]). Substituting all these numbers into Equation (19), we obtain t Hall on the order of a few months or years, in good agreement with observed intermittency timescales.
    This t Hall being quite long is a result of the helicon waves being much more slowly moving than vacuum electromagnetic waves. We essentially have a helicon standing wave in the crust, the phase of which the higher order multipoles of the magnetic field depends on. Our range of t Hall obtained above is also consistent with the oscillation period seen in the numerical simulations of [97].

4. Discussion and Conclusions

We highlighted an apparent feature seen in the distribution of the high-fraction nulling and intermittent pulsars on the log 10 P log 10 P ˙ diagram, namely their congregating near two stripes of constant polar electric potential. We analyzed the statistical significance of this feature and proposed an explanation.
Our discussions are exploratory. In particular, we caution that due to limited amounts of availability data (in, e.g., the number of intermittent pulsars, the polarization measurements, etc.), which are further contaminated by the complication of the inclination angle not being perfectly zero or π / 2 , we do not claim that the existence of the claimed stripes is firmly established. We aim instead to motivate further observational efforts towards their confirmation or rejection, ideally involving observational campaigns on the relevant pulsars, using highly sensitive radio telescopes that would be able to see weak interpulses. Such observations have the potential to directly verify the extreme inclination hypothesis.
The theory part of the note are also only semiquantitative, due to the fact that we are discussing finer details in the pulsar behavior, when consensus on even the more rudimentary pulsar emission mechanism is yet elusive. We have thus to evoke various subtitles that must combine to provide us with the observations, while each still remaining an active debate within the pulsar theory community. Under such circumstances, firm conclusions are not prudent. Nevertheless, simple order of magnitude estimates appear to offer up reasonable matches with observations, so we have adventured to discuss the explanation, with the aim to offer an example of how the aforementioned observational campaigns, if carried out, could benefit the theoretical modeling efforts.
Lastly, we wish to clarify that we concentrated only on the high-fraction nullers because the low-fraction ones may just be due to our sight line missing the carousel of sub-beams, or the extinguishing of some sub-beams [12,102,103,104,105,106] (the formation of such beams may be quite dynamic and thus capricious; see, e.g., [67] for an example), and are thus less mysterious. In particular, observational hints exist that point to differences in the underlying mechanisms driving low- and high-fraction nulling. For example, the debate as to whether nulling occurs simultaneously across frequencies could possibly be resolved if one notes that studies yielding an affirmative result appear to have examined high-fraction nullers (e.g., 75 % for B0826-34 [107]), while those answering in the negative looked at lower fraction ones instead (e.g., 1.4 % for B0809+74 [108], 15 % for B1133+16 [5]). Indeed, our proposed existence of two separate bands, for only the high-fraction nullers11, if true, represents a somewhat complicated structure and could possibly help explain the apparently conflicting conclusions reached by different studies (see, e.g., [10,12,13,14]) regarding correlations among parameters of nulling pulsars. That is, there is indeed an underlying structure, so correlations do exist, yet different relationships are obeyed by separate subpopulations; thus, varying conclusions can be drawn when studying different samples.

Funding

This work is supported by the National Natural Science Foundation of China grants 12073005, 12021003, 11503003 and 11633001, the Interdiscipline Research Funds of Beijing Normal University, and the Strategic Priority Research Program of the Chinese Academy of Sciences Grant No. XDB23000000.

Acknowledgments

We would like to thank Duncan Lorimer for bringing to our attention the existence of the O stripe that intermittent pulsars congregate around.

Conflicts of Interest

The author declares no conflict of interest.

Appendix A. Models of Pulsar Magnetospheres

Appendix A.1. Vacuum Magnetosphere

In a pulsar magnetosphere, the vacuum dipole induced potential drop Δ Φ vac across the polar caps (the base region of the field lines that reach beyond the light cylinder) is always a useful baseline parameter, providing an estimate on the typical energy of charged particles lifted out of the star that might be available to seed cascades, regardless of whether the polar cap is the actual site of cascading (such as proposed by [87,112,113]). For example, even with young pulsars where cascading happens predominantly in the current sheets, a transition of the magnetosphere into a charged-starved electrosphere state is observed by [64] to occur when the threshold photon energy for pair production becomes comparable to e Δ Φ vac .
Following the treatment of [114,115] for an aligned rotator, we have that inside of the neutron star seen as a perfect conductor, the charged particles must not experience any net force, so
E ind + 1 c Ω × r × B = 0 ,
where E ind is an electric field produced by charge redistribution in the conductor (it is not an induction field generated by the variation of the magnetic field; e.g., for an aligned dipole, we have B / t = 0 but this E ind does not vanish). Substituting in the expression for a dipolar magnetic field, one can solve for the vacuum electric scalar potential outside of the star due to the induced charge redistribution, yielding the quadrupolar
Φ vac = B 0 Ω R 5 6 c r 3 3 cos 2 θ 1 ,
where R is the stellar radius, B 0 is the magnetic field strength at the equator, and θ is measured against the rotation axis. One can then further evaluate the potential drop across the polar cap, which would be responsible for lifting charges off of the star to serve as seeds for pair production cascades ( e ± discharge). To this end, one begins by demarcating the boundary of the polar cap, by noting that the boundary field lines will touch the light cylinder tangentially at the equatorial plane (given the reflection symmetry for an aligned rotator), or θ = π / 2 . Since sin 2 θ / r is a constant along a dipolar magnetic field line, we have that the angle θ P at which the boundary line strikes the stellar surface is given (for aligned rotators) by
sin 2 θ P = R Ω c sin 2 π 2 .
Substituting into Equation (A2), we then have that
Δ Φ vac Φ θ = θ P vac Φ θ = 0 vac = B 0 Ω R 5 6 c r 3 3 sin 2 θ p = B 0 Ω 2 R 3 2 c 2 .
One then further replace B 0 by its estimate from dipole radiation (only roughly valid if dipole radiation contributes a roughly constant proportion of energy loss for any given fixed inclination angle α )
B 0 = 3 c 3 8 π 2 I R 6 sin 2 α P P ˙ ,
where I ( 2 / 5 ) M R 2 is the neutron star moment of inertia (M is stellar mass), resulting in lines of
3 2 log 10 P + 1 2 log 10 P ˙ P ( α )
marking out constant Δ Φ vac values in the log 10 P log 10 P ˙ diagram when α is held constant at a small value (because the treatment leading up to Equation (A4) is valid only for aligned rotators; we have also assumed that the relative stellar mass and radius variation between pulsars is less than order unity, as variations in sin α are). For this case, we further have explicitly that
P = log 10 sin α 5 c Δ Φ vac 3 π M R .
This treatment can also be generalized to oblique rotators, e.g., [75,116,117] give (the vector basis { e r , e θ , e ϕ } is orthonormal)
B r = B 0 R 3 r 3 cos α cos θ + sin α sin θ cos ( φ Ω t ) ,
B θ = B 0 R 3 2 r 3 cos α sin θ sin α cos θ cos ( φ Ω t ) ,
B φ = B 0 R 3 2 r 3 sin α sin ( φ Ω t ) ,
E r = B 0 Ω R 5 2 c r 4 cos α 3 cos 2 θ 1 + 3 2 sin α sin 2 θ cos ( φ Ω t ) + Q 4 π ϵ 0 r 2 ,
E θ = B 0 Ω R 3 2 c r 2 sin α R 2 r 2 cos 2 θ 1 cos ( φ Ω t ) R 2 r 2 cos α sin 2 θ ,
E φ = B 0 Ω R 3 2 c r 2 1 R 2 r 2 sin α cos θ sin ( φ Ω t ) ,
where Q is the overall charge on the star, which is zero for either dipolar or uniform magnetizations inside of the star. The expressions above are visualized in Figure A1, and they ignore the displacement current t E , so back-action on the magnetic field is neglected. This will be a common theme encountered throughout neutron star literature and this paper. It is justified for the near-star region of regular pulsars for which the rotation period is on the order of seconds, so we have R Ω / c 10 5 1 (more intuitively manifesting as the light cylinder being very far away from the star), thus | E | | B | . The displacement current expression in the full Ampère-Maxwell equation
× B = 1 c E t + 4 π c j
further brings in another such factor ( t introduces the | v | dependence), so the back-action is only very mildly perturbative ( × gives a factor of O ( 10 6 ) cm 1 suppression to the left hand side for even the global neutron star scale features, thus the ratio of the perturbations to the original magnetic field is around O ( 10 4 ) ). One should note though that the situation is different with millisecond pulsars.
Figure A1. The magnetic (yellow, light in black and white) and electric (blue, dark) field lines projected onto the x z (panels (a,d,g)), x y (panels (b,e,h)), and y z (panels (c,f,i)) planes, for the aligned (panels (ac), the rotation axis is in the z direction, and the magnetic axis is in the x z plane for the moment plotted), orthogonal (panels (df)) and moderately oblique ( α = π / 4 ; panels (gi)) rotators. The semi-transparent disk represents the stellar surface, and note that the plotting software produces spoke patterns when the field on the plane is vanishing.
Figure A1. The magnetic (yellow, light in black and white) and electric (blue, dark) field lines projected onto the x z (panels (a,d,g)), x y (panels (b,e,h)), and y z (panels (c,f,i)) planes, for the aligned (panels (ac), the rotation axis is in the z direction, and the magnetic axis is in the x z plane for the moment plotted), orthogonal (panels (df)) and moderately oblique ( α = π / 4 ; panels (gi)) rotators. The semi-transparent disk represents the stellar surface, and note that the plotting software produces spoke patterns when the field on the plane is vanishing.
Universe 07 00455 g0a1aUniverse 07 00455 g0a1b
Integrating E r along r direction to infinity then gives Φ vac equaling ( r / 3 ) E r . The second line in Equation (A11) averages to zero over a rotation period thus does not contribute to a time-average Δ Φ vac . For an oblique rotator, the potential difference
Δ Φ vac = Φ vac θ = θ P 1 Φ vac θ = θ P 2 = B 0 Ω R 5 6 c r 3 sin 2 θ P 1 sin 2 θ P 2
between the two ends of the polar cap that differ the greatest in θ (since Φ vac is θ dependent) is what’s interesting to us (in contrast, for an aligned rotator, the rim of the polar cap share the same θ P , which was why we used the θ = 0 reference point to evaluate the maximum potential drop). To evaluate θ P 1 / 2 , we note that the magnetic field outside of the star is essentially a superposition of aligned and orthogonal rotators, including and in particular B r and B θ as given by Equations (A8) and (A9). Assume that a boundary field line strikes the light cylinder tangentially at a location with angle θ LC , and thus
tan θ LC = B θ B r ,
where the right hand side is independent of r, we have then that θ LC must be independent of how far the light cylinder is to the star, or in other words independent of Ω . We thus have that
sin 2 θ P 1 / 2 α = R Ω c F 1 / 2 .
where F 1 / 2 depends on α but not Ω . For α = π / 2 , the left hand side of Equation (A17) becomes cos 2 θ P 1 / 2 which can be substituted directly into Equation (A15) resulting in once again Δ Φ vac Ω 2 B 0 , and subsequently Equation (A6) remains valid.
It is important to remember that the right hand side of Equation (A6) is α dependent, so the lines are only meaningful when comparing pulsars with similar inclinations. Therefore, the observed huddling of intermittent and high-fraction nulling pulsars near the O and A stripes (each has a slope consistent with Equation (A6)) is indicative that they should be grouped into two particular inclinations. Furthermore, at a given P and P ˙ , B 0 for larger inclination angles will be smaller, so Δ Φ vac at that point in the log 10 P log 10 P ˙ diagram will be smaller for larger inclinations. Alternatively and more intuitively, ref. [114] computed the induced charge density inside of the star to be 2 ϵ 0 Ω · B 0 = 2 ϵ 0 Ω B 0 cos α , so the associated electric field is smaller for larger inclinations. In turn, if there is a roughly common shut-down threshold Δ Φ vac value for both nearly-aligned and orthogonal pulsars, one would expect the orthogonal death line to occur further to the top-left of the diagram than the aligned death line, which is indeed consistent with the placement of the O and A stripes.
Finally, we note that the vacuum treatment of this section provides us with only the maximum possible P that would be available. In reality, since we are interested in charge-filled magnetospheres, it is not the entire Δ Φ vac , but its charge-modified version that manages or fails to drive the cascades. We thus turn to plasma filled electrospheres.

Appendix A.2. With Plasma: The Active Magnetospheres

Although [114] (GJ) already painted a semiquantitative picture of a plasma-filled magnetosphere, the precise details of the charge distribution and generation mechanism has been the subject of extensive debate. Earlier analytical and numerical work [73,84,118,119,120] showed that charges lifted from the star itself will form force-free charge-separated “electrosphere” structures (crudely stylized in Figure A3a,b below) that are rather different from the GJ prescription (see Appendix A.3 for more details). These solutions are well constructed but predicts inactive magnetospheres, as they are characterized by charge trapping (thus no current/pulsar wind flowing out to the nebula [118]). Interests in them thus waned and the search went on to look for active magnetospheres.
Later on, numerous studies (e.g., [77,82,121,122,123]) successfully showed that the GJ-like configuration is instead achieved when, in addition to charges from the star, abundant e ± pair production [113,124] by cascading is assumed, turning the magnetosphere into a high multiplicity (neutral plasma) one. In these solutions, currents flow out along the open field lines, and return through the current sheets (styled in Figure A2, see also e.g., [125] and references therein), forming an active circuit that maintains charge neutrality, thereby answering a fundamental critique leveled against the original GJ model. This is the canonical pulsar magnetosphere that has nowadays become widely accepted. We do not summarize much further details about its properties, since the review literature on it is extensive and easy to find. We do emphasize however, that the magnetic field structure for the active magnetosphere is dipolar only in the near zone (potentially containing higher multipoles immediately abutting the star), and becomes split monopolar at far away due to the strong impact of the high current density within the current sheets. On the other hand, the magnetic field for the electrospheres that we describe next do not differ much from the vacuum case of Equations (A8)–(A10) (the electrical field become shielded within charge clouds, and are no longer described by Equations (A11)–(A13), but since the back-reaction on the magnetic field is small, this change in the electric field does not alter the magnetic field).
Figure A2. The structural backbone of an active magnetosphere for an aligned rotator (sliced in half). The semi-transparent yellow surface is the light cylinder, beyond which plasma cannot co-rotate with the star without going superluminal. The white meshed surfaces are the current sheets, with the two pieces within the light cylinder enclose the closed field lines region in which the plasma particles co-rotate. The field lines residing outside will penetrate the light cylinder and sweep back in the azimuthal direction, along which the particles stream (thus not co-rotating with the star). These out-flowing pulsar wind particles could cross the magnetic field lines in the outer regions where the magnetic field strength is weak, and return through the current sheets, completing a full circuit that maintains charge neutrality. The neutron star is depicted as an orange sphere at the center, whose size is, for illustrative purposes, massively enlarged from what should be for common pulsars.
Figure A2. The structural backbone of an active magnetosphere for an aligned rotator (sliced in half). The semi-transparent yellow surface is the light cylinder, beyond which plasma cannot co-rotate with the star without going superluminal. The white meshed surfaces are the current sheets, with the two pieces within the light cylinder enclose the closed field lines region in which the plasma particles co-rotate. The field lines residing outside will penetrate the light cylinder and sweep back in the azimuthal direction, along which the particles stream (thus not co-rotating with the star). These out-flowing pulsar wind particles could cross the magnetic field lines in the outer regions where the magnetic field strength is weak, and return through the current sheets, completing a full circuit that maintains charge neutrality. The neutron star is depicted as an orange sphere at the center, whose size is, for illustrative purposes, massively enlarged from what should be for common pulsars.
Universe 07 00455 g0a2

Appendix A.3. With Plasma: The Charge-Separated Electrospheres

We summarize now the aforementioned charge-separated electrospheres. These are of a dome plus torus structure for aligned rotators [73,118,119,120], and quadrants of alternating charges for orthogonal rotators [122,123] (see Figure A3). More quantitatively, it was noted by [123] that the sharp boundaries of the charge distribution tracks closely the vacuum solution’s force-free trapping surfaces, i.e., places where E = 0 is already satisfied in the vacuum solution (A8)–(A13) (and continuity dictates that E is of opposite signs on either side of it). Such surfaces are natural discontinuous boundaries between charge clouds and vacuum gaps. On the force-free side, the charges ensure that E = 0 (and if any lapse in shielding occurs, the electric field direction is such that it drives charges towards the trap surface, which is why the cloud fills at least up to this surface), but on the vacuum side, E is non-vanishing and points in the expedient direction to collect and transport stray particles of the same (sign of) charge into the force-free side, while repelling particles of the undesirable opposite charge. This ensures that the charge clouds, and in turn the electrosphere solutions, are stable [120,126]12. Such a self-maintenance feature also prevents clouds from dissipating if they spill over beyond the trap surfaces, as would be the case when there is too much charge and the electric repulsion is severe. In other words, the trap surfaces mark only the minimal sizes of the clouds, but do not set the upper bound.
Figure A3. Schematic drawings of the electrospheres. (a) Electrosphere for aligned rotators, resembling the rotating Terrela configuration familiar to plasma physicists. (b) Electrosphere for orthogonal rotators. The meshed semi-transparent surfaces are the sharp boundaries of the charged regions. The plasmas with charges indicated by the color (red for negative, blue for positive) fill the space between these surfaces and the stellar surface (the orange surface faintly visible behind the meshed surface). There are no currents flowing in any region of the electrospheres. (c,d) The trap surfaces in vacuum magnetospheres as given by Equation (A18). The orange sphere signifies the neutron star, the red arrow is the rotation axis and the black arrow the magnetic axis.
Figure A3. Schematic drawings of the electrospheres. (a) Electrosphere for aligned rotators, resembling the rotating Terrela configuration familiar to plasma physicists. (b) Electrosphere for orthogonal rotators. The meshed semi-transparent surfaces are the sharp boundaries of the charged regions. The plasmas with charges indicated by the color (red for negative, blue for positive) fill the space between these surfaces and the stellar surface (the orange surface faintly visible behind the meshed surface). There are no currents flowing in any region of the electrospheres. (c,d) The trap surfaces in vacuum magnetospheres as given by Equation (A18). The orange sphere signifies the neutron star, the red arrow is the rotation axis and the black arrow the magnetic axis.
Universe 07 00455 g0a3
Explicitly, these trapping surfaces are (derived from vacuum—the actual boundaries of the charge clouds are deformed from these due to the presence of charges, but this higher order correction does not admit a simple analytical description)
θ = π 2 or r = 3 R | cos θ | , for α = 0 , θ = π 2 or r = R 2 cos 2 θ + 2 cos 2 ϕ sin 2 θ , for α = π 2 ,
and plotted in Figure A3c,d.
For aligned rotators, we can further write out the expression
E · B = 4 R 2 cos θ ( 3 cos 2 θ + 1 ) ,
which changes sign at
θ c = ( 1 / 2 ) cos 1 1 / 3 55 ,
meaning that charges of opposite signs are lifted out of the star in regions θ < θ c and θ > θ c , so that θ c is the angle at which the trap surface intersects the star. When combined with the trapping surfaces of Equation (A18), and the fact that charges can only move along magnetic field lines, one then arrives at the following (leading order) qualitative expectation: charges of one sign (negative if the rotation and magnetic axes are aligned, positive if anti-aligned) are lifted from the polar ( θ < θ c ) regions of the star and moved to the dome-like structures in Figure A3c along the polar magnetic field lines; charges of the opposite sign are lifted from the equatorial ( θ > θ c ) regions of the star and channeled to the equatorial plane along those dipolar magnetic field lines striking the star at θ > θ c . Because sin 2 θ / r is conserved along the dipolar field lines, there is a maximum radial extend (where the magnetic field line hitting the star at θ c crosses the equatorial plane)
r max = R sin 2 π 2 sin 2 θ c = 3 2 R ,
that the latter charges can reach along the equatorial plane. In other words, these charges do not permeate the entire plane in Figure A3c, forming instead only a torus quite close to the star as shown in Figure A3a. More specifically, one can envision that the torus is bounded by the aforementioned field line with shape r = ( 3 R / 2 ) sin 2 θ . It is observed that the charge density in the torus and dome are similar to the GJ values [73]. Such a configuration is depicted in Figure A3a, but is only a leading order approximation. In particular, while the domes co-rotate with the star, and the part of the equatorial thick disk/torus abutting the star is in co-rotation, the further away regions are in differential (even faster) rotation [70] (the differential rotation is needed to evade the restrictions laid down by [127], which only applies to axisymmetry, so no such complication needs to arise for the orthogonal rotator). The differentially rotating disk is shown to suffer from diocotron instability [77,79,80], thus exhibits non-axisymmetric charge modulations and radial expansion [77], and the actual equatorial charge distribution likely is more complicated than just laminar flow around a simple torus, with macroparticles/vortices likely showing up [76] (also note that the disk is not bound by a trap surface, so there is no aforementioned self-maintenance mechanism at work to regulate the shape of the disk; nevertheless, the disk remains confined, at least in the absence of pair injection, and cannot escape to infinity [128]).
For the orthogonal rotators on the other hand, ref. [123] noted that the region within the peanut-shaped force-free trapping surfaces agree roughly with the analytical prescription given by [122] (derived assuming a particular form for the current, among other simplifications), with the charge density being (note unlike below, the original Equation (31) in the reference did not explicitly write out the ϕ dependent factors, but its Equation (25) correspond to l = 1 and solves Equation (23) only when m = ± 1 ; also it was in SI units, which we translate into cgs)
ρ peanut = 3 Ω B 0 4 π c R 3 r 3 sin 2 θ 1 Ω 2 r 2 6 c 2 1 9 sin 2 θ cos ϕ ,
but confined to regions much closer to the stellar surface, rather than reaching outwards to regions close to the light cylinder as this equation would predict (see Figure A3b). The corresponding nonvanishing current density component is
j peanut ϕ = 3 Ω 2 B 0 4 π c R 3 r 2 sin θ sin 2 θ 1 Ω 2 r 2 6 c 2 1 9 sin 2 θ cos ϕ .
Similar to the electric field in Equations (A11)–(A13), this current is also of the order R Ω / c O ( 10 5 ) times B 0 near the stellar surface, but without the benefit of another R Ω supplementing it in Equation (A14), so its back-action on the magnetic field is even smaller than the displacement current, thus also negligible.

Notes

1
e.g., they have both been proposed as terminal quiescent states that pulsars age toward [6,7,8]; note though that whether pulsars actually evolve towards either is more subtle, depending on details of the stellar shape and higher magnetic multipoles [7,9], but this is irrelevant to the present discussion—we are only concerned about the present values of the inclination angles, regardless of whether the pulsars evolved there over time or were born into them.
2
The data for the intermittent pulsars are taken from [1,2,3,4]. The nulling data are taken from [10] and references therein, namely [11,12,13,14,15,16,17,18,19,20,21,22], and we also supplement the collection with more recent discoveries (including updates on null fractions to previously known ones) from [23,24,25,26]. We have also made use of the ATNF catalog [27,28].
3
One might notice in Figure 2 that the O stripe appears to coincide with the death line for millisecond pulsars to the bottom left corner, a feature that may have relevance to the large fraction of such pulsars being identifiable as orthogonal rotators [29,30]. We will however not pursuit this line of investigation further here, and focus on the regular population containing the high-fraction nullers and intermittent pulsars.
4
The peculiar fact that the double-pole pulsars are located much further to the top-left has already been noticed by the cited references, although a connection with nullers were not made. Unfortunately, these real-world pulsars are not perfectly aligned or orthogonal, so the death lines will not be sharp (some overshooting is to be expected; for the same reason, data point scattering in Figure 1 is unavoidable), and their inclination angles are furthermore not in fact precisely known, so a quantitative assessment of the statistical significance of this feature is difficult.
5
Arising from negligible particle inertia as compared to electromagnetic energy density, so Lorentz forces must vanish or else infinite accelerations result. This is only relevant where there are charged particles, and not required in vacuum gaps. Note also that some literature took E · B = 0 to mean force-free, even though it is only a necessary and not sufficient condition.
6
These studies suitably assumed that the star is the sole source of charges. Heuristically, ref. [66] explained that the charge-separated magnetospheres would naturally arise from turning up the magnetic field and thus the induced electric potential to strengths beyond the work function of the neutron star, so that the charges are gradually lifted out of the star but confined to nearby regions (because, e.g., they cannot cross magnetic field lines easily, due to the strong magnetic field inducing synchrotron radiation that push the particles into fundamental Landau levels).
7
A charge deficiency on the other hand would not lead to a lapse in shielding, as it only causes a charge-density wave to set up but the mean flow is still the GJ flow [74].
8
This subtler extremal case has not been scrutinized by numerical studies of generic oblique rotators, which are more interested in identifying all possible potential sites of active pair formation of a generic oblique rotator. In particular, numerical simulations tend to pick an unrealistically (at least for non-millisecond pulsars) large ratio ξ = R / R LC between stellar and light cylinder ( R LC ) radii, in order to avoid having to simultaneously handle two drastically different scales, meaning small time steps are needed to resolve the short scale but the simulation has to run for a long time (thus an enormous number of evolution steps) to see any changes in the slow scale ( 1 / ξ is roughly the number of time steps needed). This in effect massively boost the rotation angular speed of the star, which is 1 / R LC , thus boosts Δ Φ vac and masks the potential problem with orthogonal rotators. For example, ref. [82] adopted a value of ξ = 0.3 (in contrast, for a realistic regular pulsar with a period on the order of a second, ξ 10 5 ), and as such did not see any issue with pair cascade with orthogonal rotators.
9
The Hall drifting is the advection of magnetic field by free electrons, which is the most important dynamics [90] in the crust consisting of a neutral fluid with ions pinned down. Ohmic dissipation can also occur but conductivity is quite high even in the crust so it only operates on much longer timescales for pulsars of regular O ( 10 12 ) G magnetic field [91]—essentially all previous studies on the crustal magnetic field concentrate on pulsar magnetic field decay over thousands to millions of years, but even for them, Ohmic decay is subdominant (but can be non-negligible). It is nevertheless interesting to note that hypothetically, should Ohmic decay becomes important on our timescale, it would tend to dissipate higher order multipoles [90], and the “on” state would become less viable. Finally, ambipolar diffusion is expected to occur deeper down where more particle species are mobile, but, as per common industry practice, we adopt the simplest Meissner condition that the superconducting core expels the magnetic field so only the crust is important for us. This is justified because the strength of ambipolar diffusion declines as magnetic field strength cubed, so is expected to be important only for young magnetars [92].
10
Note simulation works studying magnetic field decays do not normally include this solid upper crust (see, e.g., [100]; simulations are also restricted to the axisymmetric case), because it introduces much smaller timescales than the one they are interested in, and thus requires too many time steps, but it is precisely such short timescales that are relevant for our intermittency modeling effort.
11
The Rotating Radio Transients could also be very high-fraction nullers, but they could just as well be due to the giant pulses of young pulsars that rise above the detection threshold [109]. Since we cannot really distinguish between the various possibilities (unfortunately, they cannot all be nullers, or else the population synthesis does not work out [110,111]), we exclude this class of objects from our considerations in this note and include only confirmed nulling pulsars.
12
It is to be noted that some literature claims that the GJ magnetosphere is unstable and would collapse into an electrosphere. The GJ magnetosphere studied there is not the modern version with current sheets though. The plasma examined was also not of high multiplicity/nearly-neutral. Nevertheless, they demonstrate the stability of the electrospheres under non-drastic (no turning-on of new sources of particles via pair cascade) perturbations. In any case, this mechanism only works for charge-separated, not neutral plasma [73], and is therefore irrelevant for the active magnetospheres

References

  1. Kramer, M.; Lyne, A.G.; O’Brien, J.T.; Jordan, C.A.; Lorimer, D.R. A Periodically Active Pulsar Giving Insight into Magnetospheric Physics. Science 2006, 312, 549–551. [Google Scholar] [CrossRef] [Green Version]
  2. Camilo, F.; Ransom, S.M.; Chatterjee, S.; Johnston, S.; Demorest, P. PSR J1841-0500: A Radio Pulsar That Mostly is Not There. Astrophys. J. 2012, 746, 63. [Google Scholar] [CrossRef] [Green Version]
  3. Lorimer, D.R.; Lyne, A.G.; McLaughlin, M.A.; Kramer, M.; Pavlov, G.G.; Chang, C. Radio and X-ray Observations of the Intermittent Pulsar J1832+0029. Astrophys. J. 2012, 758, 141. [Google Scholar] [CrossRef]
  4. Lyne, A.G.; Stappers, B.W.; Freire, P.C.C.; Hessels, J.W.T.; Kaspi, V.M.; Allen, B.; Bogdanov, S.; Brazier, A.; Camilo, F.; Cardoso, F.; et al. Two Long-Term Intermittent Pulsars Discovered in the PALFA Survey. Astrophys. J. 2017, 834, 72. [Google Scholar] [CrossRef] [Green Version]
  5. Bhat, N.D.R.; Gupta, Y.; Kramer, M.; Karastergiou, A.; Lyne, A.G.; Johnston, S. Simultaneous single-pulse observations of radio pulsars. V. On the broadband nature of the pulse nulling phenomenon in PSR B1133+16. Astron. Astrophys. 2007, 462, 257–268. [Google Scholar] [CrossRef]
  6. Davis, L.; Goldstein, M. Magnetic-Dipole Alignment in Pulsars. Astrophys. J. Lett. 1970, 159, L81–L86. [Google Scholar] [CrossRef]
  7. Goldreich, P. Neutron Star Crusts and Alignment of Magnetic Axes in Pulsars. Astrophys. J. Lett. 1970, 160, L11. [Google Scholar] [CrossRef]
  8. Beskin, V.S.; Gurevich, A.V.; Istomin, Y.N. Physics of the Pulsar Magnetosphere; Cambridge University Press: Cambridge, UK, 1993. [Google Scholar]
  9. Philippov, A.; Tchekhovskoy, A.; Li, J.G. Time evolution of pulsar obliquity angle from 3D simulations of magnetospheres. Mon. Not. R. Astron. Soc. 2014, 441, 1879–1887. [Google Scholar] [CrossRef] [Green Version]
  10. Wang, N.; Manchester, R.N.; Johnston, S. Pulsar nulling and mode changing. Mon. Not. R. Astron. Soc. 2007, 377, 1383–1392. [Google Scholar] [CrossRef] [Green Version]
  11. Durdin, J.M.; Large, M.I.; Little, A.G.; Manchester, R.N.; Lyne, A.G.; Taylor, J.H. An unusual pulsar—PSR 0826-34. Mon. Not. R. Astron. Soc. 1979, 186, 39P–41P. [Google Scholar] [CrossRef] [Green Version]
  12. Ritchings, R.T. Pulsar single pulse intensity measurements and pulse nulling. Mon. Not. R. Astron. Soc. 1976, 176, 249–263. [Google Scholar] [CrossRef] [Green Version]
  13. Rankin, J.M. Toward an empirical theory of pulsar emission. III—Mode changing, drifting subpulses, and pulse nulling. Astrophys. J. 1986, 301, 901–922. [Google Scholar] [CrossRef]
  14. Biggs, J.D. An analysis of radio pulsar nulling statistics. Astrophys. J. 1992, 394, 574–580. [Google Scholar] [CrossRef]
  15. Huguenin, G.R.; Taylor, J.H.; Troland, T.H. The Radio Emission from Pulsar MP 0031-07. Astrophys. J. 1970, 162, 727. [Google Scholar] [CrossRef]
  16. Lyne, A.G.; Ashworth, M. The effect of nulls upon subpulse drift in PSRs 0809+74 and 0818-13. Mon. Not. R. Astron. Soc. 1983, 204, 519–536. [Google Scholar] [CrossRef]
  17. Weisberg, J.M.; Armstrong, B.K.; Backus, P.R.; Cordes, J.M.; Boriakoff, V.; Ferguson, D.C. High-sensitivity observations of 28 pulsars. Astron. J. 1986, 92, 621–626. [Google Scholar] [CrossRef]
  18. Lorimer, D.R.; Camilo, F.; Xilouris, K.M. Arecibo Timing Observations of 17 Pulsars along the Galactic Plane. Astron. J. 2002, 123, 1750–1756. [Google Scholar] [CrossRef]
  19. Lewandowski, W.; Wolszczan, A.; Feiler, G.; Konacki, M.; Sołtysiński, T. Arecibo Timing and Single-Pulse Observations of Eighteen Pulsars. Astrophys. J. 2004, 600, 905–913. [Google Scholar] [CrossRef]
  20. Faulkner, A.J.; Stairs, I.H.; Kramer, M.; Lyne, A.G.; Hobbs, G.; Possenti, A.; Lorimer, D.R.; Manchester, R.N.; McLaughlin, M.A.; D’Amico, N.; et al. The Parkes Multibeam Pulsar Survey—V. Finding binary and millisecond pulsars. Mon. Not. R. Astron. Soc. 2004, 355, 147–158. [Google Scholar] [CrossRef] [Green Version]
  21. Redman, S.L.; Wright, G.A.E.; Rankin, J.M. Pulsar PSR B2303+30: A single system of drifting subpulses, moding and nulling. Mon. Not. R. Astron. Soc. 2005, 357, 859–872. [Google Scholar] [CrossRef] [Green Version]
  22. Young, N.J.; Weltevrede, P.; Stappers, B.W.; Lyne, A.G.; Kramer, M. Long-term observations of three nulling pulsars. Mon. Not. R. Astron. Soc. 2015, 449, 1495–1504. [Google Scholar] [CrossRef] [Green Version]
  23. Burke-Spolaor, S.; Johnston, S.; Bailes, M.; Bates, S.D.; Bhat, R.; Burgay, M.; Champion, D.; D’Amico, N.; Keith, M.J.; Kramer, M.; et al. The High Time Resolution Universe Pulsar Survey—V. Single-pulse energetics and modulation properties of 315 pulsars. Mon. Not. R. Astron. Soc. 2012, 423, 1351–1367. [Google Scholar] [CrossRef] [Green Version]
  24. Gajjar, V.; Joshi, B.C.; Kramer, M. A survey of nulling pulsars using the Giant Meterwave Radio Telescope. Mon. Not. R. Astron. Soc. 2012, 424, 1197–1205. [Google Scholar] [CrossRef] [Green Version]
  25. Gajjar, V.; Joshi, B.C.; Wright, G. On the long nulls of PSRs J1738-2330 and J1752+2359. Mon. Not. R. Astron. Soc. 2014, 439, 221–233. [Google Scholar] [CrossRef] [Green Version]
  26. Naidu, A.; Joshi, B.C.; Manoharan, P.; Krishnakumar, M. Detection of long nulls in PSR B1706-16, a pulsar with large timing irregularities. Mon. Not. R. Astron. Soc. 2018, 475, 2375–2382. [Google Scholar] [CrossRef] [Green Version]
  27. Manchester, R.N.; Hobbs, G.B.; Teoh, A.; Hobbs, M. The Australia Telescope National Facility Pulsar Catalogue. Astron. J. 2005, 129, 1993–2006. [Google Scholar] [CrossRef]
  28. ATNF. Australia Telescope National Facility Pulsar Catalogue. 2005. Available online: http://www.atnf.csiro.au/research/pulsar/psrcat/ (accessed on 16 December 2020).
  29. Jayawardhana, R.; Grindlay, J.E. Pulse profiles of millisecond pulsars. In IAU Colloquium 160: Pulsars: Problems and Progress; Johnston, S., Walker, M.A., Bailes, M., Eds.; Astronomical Society of the Pacific Conference Series; Cambridge University Press: Cambridge, UK, 1996; Volume 105, p. 231. [Google Scholar]
  30. Chen, K.; Ruderman, M. Origin and radio pulse properties of millisecond pulsars. Astrophys. J. 1993, 408, 179–185. [Google Scholar] [CrossRef]
  31. Radhakrishnan, V.; Cooke, D.J. Magnetic Poles and the Polarization Structure of Pulsar Radiation. Astrophys. Lett. 1969, 3, 225. [Google Scholar]
  32. Malov, I.F. The pulsar PSR B1931+24 as an orthogonal rotator. Astron. Rep. 2007, 51, 477–480. [Google Scholar] [CrossRef]
  33. Kloumann, I.M.; Rankin, J.M. On the long and short nulls, modes and interpulse emission of radio pulsar B1944+17. Mon. Not. R. Astron. Soc. 2010, 408, 40–52. [Google Scholar] [CrossRef] [Green Version]
  34. Young, N.J.; Weltevrede, P.; Stappers, B.W.; Lyne, A.G.; Kramer, M. On the apparent nulls and extreme variability of PSR J1107-5907. Mon. Not. R. Astron. Soc. 2014, 442, 2519–2533. [Google Scholar] [CrossRef] [Green Version]
  35. Maciesiak, K.; Gil, J.; Ribeiro, V.A.R.M. On the pulse-width statistics in radio pulsars—I. Importance of the interpulse emission. Mon. Not. R. Astron. Soc. 2011, 414, 1314–1328. [Google Scholar] [CrossRef] [Green Version]
  36. Rickett, B.J.; Lyne, A.G. Further Impulsive Emission from CP 0950. Nature 1968, 218, 934–935. [Google Scholar] [CrossRef]
  37. Gil, J. PSR 0950 + 08—A model for interpulse emission. Astron. Astrophys. 1983, 127, 267–270. [Google Scholar]
  38. Gil, J. Interpulse beams and profile components. Astrophys. J. 1985, 299, 154–160. [Google Scholar] [CrossRef]
  39. Kolonko, M.; Gil, J.; Maciesiak, K. On the pulse-width statistics in radio pulsars. Astron. Astrophys. 2004, 428, 943–951. [Google Scholar] [CrossRef] [Green Version]
  40. Taylor, J.H.; Manchester, R.N.; Lyne, A.G. Catalog of 558 pulsars. Astrophys. J. 1993, 88, 529–568. [Google Scholar] [CrossRef]
  41. Weltevrede, P.; Johnston, S. Profile and polarization characteristics of energetic pulsars. Mon. Not. R. Astron. Soc. 2008, 391, 1210–1226. [Google Scholar] [CrossRef] [Green Version]
  42. D’Amico, N.; Stappers, B.W.; Bailes, M.; Martin, C.E.; Bell, J.F.; Lyne, A.G.; Manchester, R.N. The Parkes Southern Pulsar Survey—III. Timing of long-period pulsars. Mon. Not. R. Astron. Soc. 1998, 297, 28–40. [Google Scholar] [CrossRef] [Green Version]
  43. Morris, D.J.; Hobbs, G.; Lyne, A.G.; Stairs, I.H.; Camilo, F.; Manchester, R.N.; Possenti, A.; Bell, J.F.; Kaspi, V.M.; Amico, N.D.; et al. The Parkes Multibeam Pulsar Survey—II. Discovery and timing of 120 pulsars. Mon. Not. R. Astron. Soc. 2002, 335, 275–290. [Google Scholar] [CrossRef] [Green Version]
  44. Kramer, M.; Bell, J.F.; Manchester, R.N.; Lyne, A.G.; Camilo, F.; Stairs, I.H.; D’Amico, N.; Kaspi, V.M.; Hobbs, G.; Morris, D.J.; et al. The Parkes Multibeam Pulsar Survey—III. Young pulsars and the discovery and timing of 200 pulsars. Mon. Not. R. Astron. Soc. 2003, 342, 1299–1324. [Google Scholar] [CrossRef] [Green Version]
  45. Hobbs, G.; Faulkner, A.; Stairs, I.H.; Camilo, F.; Manchester, R.N.; Lyne, A.G.; Kramer, M.; D’Amico, N.; Kaspi, V.M.; Possenti, A.; et al. The Parkes multibeam pulsar survey—IV. Discovery of 180 pulsars and parameters for 281 previously known pulsars. Mon. Not. R. Astron. Soc. 2004, 352, 1439–1472. [Google Scholar] [CrossRef]
  46. Lorimer, D.R.; Faulkner, A.J.; Lyne, A.G.; Manchester, R.N.; Kramer, M.; McLaughlin, M.A.; Hobbs, G.; Possenti, A.; Stairs, I.H.; Camilo, F.; et al. The Parkes Multibeam Pulsar Survey—VI. Discovery and timing of 142 pulsars and a Galactic population analysis. Mon. Not. R. Astron. Soc. 2006, 372, 777–800. [Google Scholar] [CrossRef] [Green Version]
  47. Manchester, R.; Lyne, A.; Camilo, F.; Bell, J.; Kaspi, V.; D’Amico, N.; McKay, N.; Crawford, F.; Stairs, I.; Possenti, A.; et al. The Parkes multi-beam pulsar survey—I. Observing and data analysis systems, discovery and timing of 100 pulsars. Mon. Not. R. Astron. Soc. 2001, 328, 17–35. [Google Scholar] [CrossRef]
  48. Janssen, G.H.; Stappers, B.W.; Braun, R.; Van Straten, W.; Edwards, R.T.; Rubio-Herrera, E.; Van Leeuwen, J.; Weltevrede, P. Discovery and timing of the first 8gr8 Cygnus survey pulsars. Astron. Astrophys. 2009, 498, 223–231. [Google Scholar] [CrossRef]
  49. Keith, M.J.; Eatough, R.P.; Lyne, A.G.; Kramer, M.; Possenti, A.; Camilo, F.; Manchester, R.N. Discovery of 28 pulsars using new techniques for sorting pulsar candidates. Mon. Not. R. Astron. Soc. 2009, 395, 837–846. [Google Scholar] [CrossRef] [Green Version]
  50. Keith, M.J.; Johnston, S.; Weltevrede, P.; Kramer, M. Polarization measurements of five pulsars with interpulses. Mon. Not. R. Astron. Soc. 2010, 402, 745–752. [Google Scholar] [CrossRef] [Green Version]
  51. Camilo, F.; Ray, P.; Ransom, S.; Burgay, M.; Johnson, T.J.; Kerr, M.; Gotthelf, E.V.; Halpern, J.P.; Reynolds, J.; Romani, R.W.; et al. Radio Detection of LAT PSRs J1741-2054 and J2032+4127: No Longer Just Gamma-ray Pulsars. Astrophys. J. 2009, 705, 1–13. [Google Scholar] [CrossRef]
  52. Malov, I.F.; Nikitina, E.B. The magnetospheric structure of radio pulsars with interpulses. Astron. Rep. 2013, 57, 833–843. [Google Scholar] [CrossRef]
  53. Spitkovsky, A. Time-dependent Force-free Pulsar Magnetospheres: Axisymmetric and Oblique Rotators. Astrophys. J. 2006, 648, L51–L54. [Google Scholar] [CrossRef]
  54. Beskin, V.S.; Nokhrina, E.E. On the role of the current loss in radio pulsar evolution. Astrophys. Space Sci. 2007, 308, 569–573. [Google Scholar] [CrossRef] [Green Version]
  55. Gurevich, A.V.; Istomin, Y.N. The energy loss of a rotating magnetized neutron star. Mon. Not. R. Astron. Soc. 2007, 377, 1663–1667. [Google Scholar] [CrossRef] [Green Version]
  56. Gil, J. The frequency dependence of component separation in pulsar integrated profiles. Part 1. Low frequencies. Acta Phys. Pol. 1981, 12, 1081. [Google Scholar]
  57. Rankin, J.M. Toward an empirical theory of pulsar emission. IV—Geometry of the core emission region. Astrophys. J. 1990, 352, 247–257. [Google Scholar] [CrossRef]
  58. Manchester, R.N.; Taylor, J.H. Pulsars; W. H. Freeman: San Francisco, CA, USA, 1997. [Google Scholar]
  59. Young, M.D.T.; Chan, L.S.; Burman, R.R.; Blair, D.G. Pulsar magnetic alignment and the pulsewidth-age relation. Mon. Not. R. Astron. Soc. 2010, 402, 1317–1329. [Google Scholar] [CrossRef] [Green Version]
  60. Rankin, J.M. Toward an empirical theory of pulsar emission. VI—The geometry of the conal emission region. Astrophys. J. 1993, 405, 285–297. [Google Scholar] [CrossRef]
  61. Gil, J.A.; Kijak, J.; Seiradakis, J.H. On the Two-Dimensional Structure of Pulsar Beams. Astron. Astrophys. 1993, 272, 268. [Google Scholar]
  62. Kramer, M.; Wielebinski, R.; Jessner, A.; Gil, J.A.; Seiradakis, J.H. Geometrical analysis of average pulsar profiles using multi-component Gaussian FITS at several frequencies. I. Method and analysis. Astron. Astrophys. 1994, 107, 515–526. [Google Scholar]
  63. Timokhin, A.N.; Arons, J. Current flow and pair creation at low altitude in rotation-powered pulsars’ force-free magnetospheres: Space charge limited flow. Mon. Not. R. Astron. Soc. 2013, 429, 20–54. [Google Scholar] [CrossRef] [Green Version]
  64. Chen, A.; Beloborodov, A.M. Electrodynamics of Axisymmetric Pulsar Magnetosphere with Electron-Positron Discharge: A Numerical Experiment. Astrophys. J. 2014, 795, L22. [Google Scholar] [CrossRef]
  65. Philippov, A.A.; Cerutti, B.; Tchekhovskoy, A.; Spitkovsky, A. Ab Initio Pulsar Magnetosphere: The Role of General Relativity. Astrophys. J. Lett. 2015, 815, L19. [Google Scholar] [CrossRef] [Green Version]
  66. Michel, F.C.; Smith, I.A. Axisymmetric Rotator Models as Pulsars: Abject Failure. In Revista Mexicana de Astronomia y Astrofisica; Cantó, J., Rodríguez, L.F., Eds.; Revista Mexicana de Astronomia y Astrofisica Conference Series; UNAM: Mexico City, Mexico, 2001; Volume 27, pp. 168–175. [Google Scholar]
  67. Zhang, F. Pulsar current sheet Čerenkov radiation. Astron. Astrophys. 2018, 612, A24. [Google Scholar] [CrossRef] [Green Version]
  68. Philippov, A.; Timokhin, A.; Spitkovsky, A. Origin of Pulsar Radio Emission. Phys. Rev. Lett. 2020, 124, 245101. [Google Scholar] [CrossRef]
  69. Michel, F.C. Theory of Pulsar Action as Revealed by PSR B1931+24 (the “intermittent Pulsar”). In AAS/High Energy Astrophysics Division #11; Bulletin of the American Astronomical Society; American Astronomical Society: Washington, DC, USA, 2010; Volume 42, p. 682. [Google Scholar]
  70. Pétri, J.; Heyvaerts, J.; Bonazzola, S. Global static electrospheres of charged pulsars. Astron. Astrophys. 2002, 384, 414–432. [Google Scholar] [CrossRef] [Green Version]
  71. Li, J.; Spitkovsky, A.; Tchekhovskoy, A. On the Spin-down of Intermittent Pulsars. Astrophys. J. 2012, 746, L24. [Google Scholar] [CrossRef] [Green Version]
  72. Chen, A.Y.; Cruz, F.; Spitkovsky, A. Filling the Magnetospheres of Weak Pulsars. arXiv 2019, arXiv:1911.00059. [Google Scholar] [CrossRef] [Green Version]
  73. Smith, I.A.; Michel, F.C.; Thacker, P.D. Numerical simulations of aligned neutron star magnetospheres. Mon. Not. R. Astron. Soc. 2001, 322, 209–217. [Google Scholar] [CrossRef] [Green Version]
  74. Shibata, S. The field-aligned accelerator in the pulsar magnetosphere. Mon. Not. R. Astron. Soc. 1997, 287, 262–270. [Google Scholar] [CrossRef] [Green Version]
  75. Pétri, J. Theory of pulsar magnetosphere and wind. J. Plasma Phys. 2016, 82, 635820502. [Google Scholar] [CrossRef]
  76. Pétri, J. Non-linear evolution of the diocotron instability in a pulsar electrosphere: Two-dimensional particle-in-cell simulations. Astron. Astrophys. 2009, 503, 1–12. [Google Scholar] [CrossRef] [Green Version]
  77. Philippov, A.A.; Spitkovsky, A. Ab Initio Pulsar Magnetosphere: Three-dimensional Particle-in-cell Simulations of Axisymmetric Pulsars. Astrophys. J. Lett. 2014, 785, L33. [Google Scholar] [CrossRef] [Green Version]
  78. Brambilla, G.; Kalapotharakos, C.; Timokhin, A.N.; Harding, A.K.; Kazanas, D. Electron-Positron Pair Flow and Current Composition in the Pulsar Magnetosphere. Astrophys. J. 2018, 858, 81. [Google Scholar] [CrossRef] [Green Version]
  79. Pétri, J.; Heyvaerts, J.; Bonazzola, S. Diocotron instability in pulsar electrospheres. I. Linear analysis. Astron. Astrophys. 2002, 387, 520–530. [Google Scholar] [CrossRef] [Green Version]
  80. Pétri, J. The diocotron instability in a pulsar cylindrical electrosphere. Astron. Astrophys. 2007, 464, 135–142. [Google Scholar] [CrossRef]
  81. Erber, T. High-Energy Electromagnetic Conversion Processes in Intense Magnetic Fields. Rev. Mod. Phys. 1966, 38, 626–659. [Google Scholar] [CrossRef]
  82. Philippov, A.A.; Spitkovsky, A.; Cerutti, B. Ab Initio Pulsar Magnetosphere: Three-dimensional Particle-in-cell Simulations of Oblique Pulsars. Astrophys. J. Lett. 2015, 801, L19. [Google Scholar] [CrossRef] [Green Version]
  83. Beloborodov, A.M. Polar-Cap Accelerator and Radio Emission from Pulsars. Astrophys. J. Lett. 2008, 683, L41. [Google Scholar] [CrossRef] [Green Version]
  84. Leahy, D.A. Relativistic magnetospheres—Corotating and streaming formulations for the charge filled case. Astron. Astrophys. 1985, 149, 57–64. [Google Scholar]
  85. Cheng, K.S.; Ho, C.; Ruderman, M. Energetic radiation from rapidly spinning pulsars. I—Outer magnetosphere gaps. II—VELA and Crab. Astrophys. J. 1986, 300, 500–539. [Google Scholar] [CrossRef]
  86. Cerutti, B.; A Philippov, A.; Parfrey, K.; Spitkovsky, A. Particle acceleration in axisymmetric pulsar current sheets. Mon. Not. R. Astron. Soc. 2015, 448, 606–619. [Google Scholar] [CrossRef] [Green Version]
  87. Ruderman, M.A.; Sutherland, P.G. Theory of pulsars: Polar caps, sparks, and coherent microwave radiation. Astrophys. J. 1975, 196, 51. [Google Scholar] [CrossRef]
  88. Arons, J. Pulsar Emission: Where to Go. In Astrophysics and Space Science Library; Becker, W., Ed.; Springer: Berlin/Heidelberg, Germany, 2009; Volume 357, p. 373. [Google Scholar]
  89. Timokhin, A. On the force-free magnetosphere of an aligned rotator. Mon. Not. R. Astron. Soc. 2006, 368, 1055–1072. [Google Scholar] [CrossRef]
  90. Gourgouliatos, K.N.; Cumming, A. Hall Attractor in Axially Symmetric Magnetic Fields in Neutron Star Crusts. Phys. Rev. Lett. 2014, 112, 171101. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  91. Gourgouliatos, K.N.; Cumming, A.; Lyutikov, M.; Reisenegger, A. Hall Equilibria: Solutions with toroidal and poloidal magnetic fields in Neutron Star Crusts. arXiv 2013, arXiv:1305.0149. [Google Scholar]
  92. Pons, J.A.; Viganò, D. Magnetic, thermal and rotational evolution of isolated neutron stars. arXiv 2019, arXiv:1911.03095. [Google Scholar] [CrossRef] [Green Version]
  93. Spitzer, L. Physics of fuly ionized gases. In Interscience, 2nd ed.; McGraw-Hill: New York, NY, USA, 1962. [Google Scholar]
  94. Boswell, R.W.; Chen, F.F. Helicons-the early years. IEEE Trans. Plasma Sci. 1997, 25, 1229–1244. [Google Scholar] [CrossRef] [Green Version]
  95. Aigrain, P.R. Proceedings of the International Conference on Semiconductor Physics, Prague, 1960; Academic Press: Cambridge, MA, USA, 1961; p. 224. [Google Scholar]
  96. Maxfield, B.W. Helicon Waves in Solids. Am. J. Phys. 1969, 37, 241–269. [Google Scholar] [CrossRef]
  97. Shalybkov, D.A.; Urpin, V.A. The Hall effect and the decay of magnetic fields. Astron. Astrophys. 1997, 321, 685–690. [Google Scholar]
  98. Cordes, J.M. Pulsar State Switching from Markov Transitions and Stochastic Resonance. Astrophys. J. 2013, 775, 47. [Google Scholar] [CrossRef] [Green Version]
  99. Goldreich, P.; Reisenegger, A. Magnetic field decay in isolated neutron stars. Astrophys. J. 1992, 395, 250–258. [Google Scholar] [CrossRef] [Green Version]
  100. Gourgouliatos, K.N.; Cumming, A. Hall effect in neutron star crusts: Evolution, endpoint and dependence on initial conditions. Mon. Not. R. Astron. Soc. 2014, 438, 1618–1629. [Google Scholar] [CrossRef] [Green Version]
  101. Rüster, S.B.; Hempel, M.; Schaffner-Bielich, J. Outer crust of nonaccreting cold neutron stars. Phys. Rev. C 2006, 73, 035804. [Google Scholar] [CrossRef] [Green Version]
  102. Backer, D.C. Pulsar Nulling Phenomena. Nature 1970, 228, 42–43. [Google Scholar] [CrossRef] [PubMed]
  103. Arons, J. Electron positron pairs in radio pulsars. In Positron-Electron Pairs in Astrophysics; Burns, M.L., Harding, A.K., Ramaty, R., Eds.; American Institute of Physics Conference Series; AIP Publishing LLC: New York, NY, USA, 1983; Volume 101, pp. 163–193. [Google Scholar]
  104. Deshpande, A.A.; Rankin, J.M. The topology and polarization of sub-beams associated with the ‘drifting’ sub-pulse emission of pulsar B0943+10—I. Analysis of Arecibo 430- and 111-MHz observations. Mon. Not. R. Astron. Soc. 2001, 322, 438–460. [Google Scholar] [CrossRef] [Green Version]
  105. Rankin, J.M.; Wright, G.A.E. The ‘periodic nulls’ of radio pulsar J1819+1305. Mon. Not. R. Astron. Soc. 2008, 385, 1923–1930. [Google Scholar] [CrossRef] [Green Version]
  106. Redman, S.L.; Rankin, J.M. On the randomness of pulsar nulls. Mon. Not. R. Astron. Soc. 2009, 395, 1529–1532. [Google Scholar] [CrossRef] [Green Version]
  107. Bhattacharyya, B.; Gupta, Y.; Gil, J. Results from multifrequency observations of PSR B0826-34. Mon. Not. R. Astron. Soc. 2008, 383, 1538–1550. [Google Scholar] [CrossRef] [Green Version]
  108. Bartel, N.; Sieber, W. Simultaneous single pulse observations of radio pulsars at two widely spaced frequencies. Astron. Astrophys. 1978, 70, 307–310. [Google Scholar]
  109. McLaughlin, M.A.; Cordes, J.M. Searches for Giant Pulses from Extragalactic Pulsars. Astrophys. J. 2003, 596, 982–996. [Google Scholar] [CrossRef] [Green Version]
  110. McLaughlin, M.A.; Lyne, A.G.; Lorimer, D.R.; Kramer, M.; Faulkner, A.J.; Manchester, R.N.; Cordes, J.M.; Camilo, F.; Possenti, A.; Stairs, I.H.; et al. Transient radio bursts from rotating neutron stars. Nature 2006, 439, 817–820. [Google Scholar] [CrossRef]
  111. Burke-Spolaor, S. Rotating Radio Transients and their place among pulsars. In Neutron Stars and Pulsars: Challenges and Opportunities after 80 Years; van Leeuwen, J., Ed.; IAU Symposium; Cambridge University Press: Cambridge, UK, 2013; Volume 291, pp. 95–100. [Google Scholar]
  112. Sturrock, P.A. Pulsar Radiation Mechanisms. Nature 1970, 227, 465–470. [Google Scholar] [CrossRef] [PubMed]
  113. Sturrock, P.A. A Model of Pulsars. Astrophys. J. 1971, 164, 529. [Google Scholar] [CrossRef]
  114. Goldreich, P.; Julian, W.H. Pulsar electrodynamics. Astrophys. J. 1969, 157, 869. [Google Scholar] [CrossRef]
  115. Lorimer, D.R.; Kramer, M. Handbook of Pulsar Astronomy; Cambridge University Press: Cambridge, UK, 2012. [Google Scholar]
  116. Deutsch, A.J. The electromagnetic field of an idealized star in rigid rotation in vacuo. Ann. Astrophys. 1955, 18, 1. [Google Scholar]
  117. Michel, F.; Li, H. Electrodynamics of neutron stars. Phys. Rep. 1999, 318, 227–297. [Google Scholar] [CrossRef]
  118. Jackson, E.A. A new pulsar atmospheric model. I—Aligned magnetic and rotational axes. Astrophys. J. 1976, 206, 831–841. [Google Scholar] [CrossRef]
  119. Krause-Polstorff, J.; Michel, F.C. Electrosphere of an aligned magnetized neutron star. Mon. Not. R. Astron. Soc. 1985, 213, 43P–49P. [Google Scholar] [CrossRef] [Green Version]
  120. Krause-Polstorff, J.; Michel, F.C. Pulsar space charging. Astron. Astrophys. 1985, 144, 72–80. [Google Scholar]
  121. Contopoulos, I.; Kazanas, D.; Fendt, C. The Axisymmetric Pulsar Magnetosphere. Astrophys. J. 1999, 511, 351–358. [Google Scholar] [CrossRef] [Green Version]
  122. Parish, J.L. Pulsar near fields. Astrophys. J. 1974, 193, 225–229. [Google Scholar] [CrossRef]
  123. McDonald, J.; Shearer, A. Investigations of the Magnetospheric Plasma Distribution in the Vicinity of a Pulsar. I. Basic Formulation. Astrophys. J. 2009, 690, 13–19. [Google Scholar] [CrossRef] [Green Version]
  124. Cheng, A.; Ruderman, M.; Sutherland, P. Current flow in pulsar magnetospheres. Astrophys. J. 1976, 203, 209–212. [Google Scholar] [CrossRef]
  125. Gralla, S.E.; Lupsasca, A.; Philippov, A. Pulsar Magnetospheres: Beyond the Flat Spacetime Dipole. Astrophys. J. 2016, 833, 258. [Google Scholar] [CrossRef] [Green Version]
  126. Aly, J.J. A stability property of a force-free surface bounding a vacuum gap. Astron. Astrophys. 2005, 429, 779–784. [Google Scholar] [CrossRef] [Green Version]
  127. Pilipp, W.G. On the Electrodynamic Equilibrium of a Space Charge Region around a Rotating Neutron Star with an Aligned Magnetic Field. Astrophys. J. 1974, 190, 391–402. [Google Scholar] [CrossRef]
  128. Aly, J.J. Conservation laws and theorems of confinement and stability for a charged equatorial disk in a pulsar magnetosphere. Astron. Astrophys. 2005, 434, 405–414. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Placement of intermittent (red ⊗) and high-fraction nulling (black dots) pulsars on the log 10 P log 10 P ˙ diagram. (a) Those nulling pulsars with >15% null fraction (chosen due to B1133+16 apparently belonging to the low-fraction category given its non-simultaneous nulling across frequency bands [5]; we also caution that some nulling may be fake—simply due to emission getting into a weaker mode below detection threshold, so a definitive null fraction threshold is unlikely meaningful), as well as the intermittent pulsars (essentially the same as Figure 7 of [4], but with additional data points). The stripes correspond to the one standard deviation band of P associated with all the points clustered around it. (b) Same as (a), but with pulsar names displaced near the data points. Some names have to be offset to avoid overlapping. The labels for the intermittent pulsars are placed at the same vertical height as the corresponding data point. If two labels are similar in height, their horizontal displacements reflect those of the data point. The labels for the nulling pulsars are arranged to coincide with the data points in their horizontal location, and when the labels become too crowded, their vertical height becomes the secondary discriminator that determines the corresponding data point.
Figure 1. Placement of intermittent (red ⊗) and high-fraction nulling (black dots) pulsars on the log 10 P log 10 P ˙ diagram. (a) Those nulling pulsars with >15% null fraction (chosen due to B1133+16 apparently belonging to the low-fraction category given its non-simultaneous nulling across frequency bands [5]; we also caution that some nulling may be fake—simply due to emission getting into a weaker mode below detection threshold, so a definitive null fraction threshold is unlikely meaningful), as well as the intermittent pulsars (essentially the same as Figure 7 of [4], but with additional data points). The stripes correspond to the one standard deviation band of P associated with all the points clustered around it. (b) Same as (a), but with pulsar names displaced near the data points. Some names have to be offset to avoid overlapping. The labels for the intermittent pulsars are placed at the same vertical height as the corresponding data point. If two labels are similar in height, their horizontal displacements reflect those of the data point. The labels for the nulling pulsars are arranged to coincide with the data points in their horizontal location, and when the labels become too crowded, their vertical height becomes the secondary discriminator that determines the corresponding data point.
Universe 07 00455 g001
Figure 2. For reference, the location of the stripes in relation to all known pulsars with a negative P ˙ (a few pulsars are measured to be spinning up, possibly due to the gravitational effects associated with their being inside globular clusters) is shown in this panel, data are from the ATNF catalog.
Figure 2. For reference, the location of the stripes in relation to all known pulsars with a negative P ˙ (a few pulsars are measured to be spinning up, possibly due to the gravitational effects associated with their being inside globular clusters) is shown in this panel, data are from the ATNF catalog.
Universe 07 00455 g002
Figure 3. Pulsars (not necessarily nulling or intermittent) exhibiting interpulses of both double-pole (red crosses; almost orthogonal rotators) and single-pole (black circled crosses; almost aligned rotators) as compiled by [35] (with data originally from [39,40,41,42,43,44,45,46,47,48,49,50,51]), as well as [52]. The red circled crosses are pulsars for which [35,52] made opposite single or double-pole assignments.
Figure 3. Pulsars (not necessarily nulling or intermittent) exhibiting interpulses of both double-pole (red crosses; almost orthogonal rotators) and single-pole (black circled crosses; almost aligned rotators) as compiled by [35] (with data originally from [39,40,41,42,43,44,45,46,47,48,49,50,51]), as well as [52]. The red circled crosses are pulsars for which [35,52] made opposite single or double-pole assignments.
Universe 07 00455 g003
Figure 4. (a) The pulse width W as a function of the inclination angle α and viewing impact angle ζ (see above Equation (6) for its definition and below Equation (4) for the definition of ρ ). All values are in degrees. The semi-transparent surfaces correspond to ρ = 20 (yellow) and ρ = 10 (red). The Contours of constant W span wide α ranges, suggesting that α estimates based off of pulse width can not be very precise. (b) The contours of W on the ( α , ζ ) plane, for the case of ρ = 20 .
Figure 4. (a) The pulse width W as a function of the inclination angle α and viewing impact angle ζ (see above Equation (6) for its definition and below Equation (4) for the definition of ρ ). All values are in degrees. The semi-transparent surfaces correspond to ρ = 20 (yellow) and ρ = 10 (red). The Contours of constant W span wide α ranges, suggesting that α estimates based off of pulse width can not be very precise. (b) The contours of W on the ( α , ζ ) plane, for the case of ρ = 20 .
Universe 07 00455 g004aUniverse 07 00455 g004b
Table 1. Silhouette values indicating the statistical significance of the clustering in P of the intermittent pulsars and high-fraction nullers.
Table 1. Silhouette values indicating the statistical significance of the clustering in P of the intermittent pulsars and high-fraction nullers.
PulsarSilhouetteSilhouette
(Nuller Stripe O )vs. General Pop.vs. Stripe A
J1525-54170.840.88
J1709-16400.770.84
J1255-61310.850.89
J1717-40540.630.60
J1502-56530.850.88
J1820-05090.750.83
J1634-51070.850.88
J1752+23590.800.83
Average0.790.83
PulsarSilhouetteSilhouette
(Nuller Stripe A )vs. General Pop.vs. Stripe O
J2037+19420.800.82
J1725-40430.760.74
J1853+05050.540.39
J1944+17550.670.71
J1107-59070.860.86
J1049-58330.860.87
J1702-44280.850.86
J1727-27390.860.86
J1916+10230.760.78
J1920+10400.850.85
J0034-07210.870.87
J0528+22000.800.79
J0754+32310.860.86
J0826-34170.750.78
J1115+50300.860.86
J1649+25330.860.86
J1744-39220.840.85
J1945-00400.870.87
J1946+18050.850.86
J2113+46440.830.82
J2321+60240.850.84
J1738-23300.700.65
Average0.810.80
PulsarSilhouette
(Intermittent)vs. General Pop.
J1832+00290.76
J1841-05000.66
J1910+05170.59
J1929+13570.71
J1933+24210.72
Average0.69
Table 2. The ratio of Monte Carlo simulated samples (with the same number of pulsars in each sample as in the observed population tagged by the first column) with an average silhouette value equal or above that of the actual observed population. Note that sample size refers to the number of Monte Carlo simulations, not the number of pulsars in the mock up samples.
Table 2. The ratio of Monte Carlo simulated samples (with the same number of pulsars in each sample as in the observed population tagged by the first column) with an average silhouette value equal or above that of the actual observed population. Note that sample size refers to the number of Monte Carlo simulations, not the number of pulsars in the mock up samples.
Pulsar PopulationSample SizeRatio ≥ Observed
Intermittent 10 4 0.69 %
10 5 0.65 %
10 6 0.67 %
10 7 0.67 %
Nuller Stripe O 10 4 0 %
10 5 0.0020 %
10 6 0.0018 %
10 7 0.0020 %
10 8 0.0020 %
Nuller Stripe A 10 4 0 %
10 5 0 %
10 6 0 %
10 7 0 %
10 8 0 %
Table 3. The Monte Carlo simulation results for the penumbra values. The qualification “Uncertain to Single” means the circled red crosses of Figure 3 are grouped with the single-poled population. Sample size here refers to the number of Monte Carlo simulations carried out.
Table 3. The Monte Carlo simulation results for the penumbra values. The qualification “Uncertain to Single” means the circled red crosses of Figure 3 are grouped with the single-poled population. Sample size here refers to the number of Monte Carlo simulations carried out.
Pulsar PopulationSample SizeRatio ≤ Observed
Uncertain to Single
Single-poled 10 4 68 %
10 5 68 %
10 6 68 %
10 7 68 %
Double-poled 10 4 0 %
10 5 0.0030 %
10 6 0.0015 %
10 7 0.0012 %
10 8 0.0012 %
Pulsar PopulationSample SizeRatio ≤ Observed
Uncertain to Double
Single-poled 10 4 88 %
10 5 87 %
10 6 87 %
10 7 87 %
Double-poled 10 4 0 %
10 5 0 %
10 6 0.00030 %
10 7 0.00023 %
10 8 0.00021 %
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, F. A Magnetospheric Dichotomy for Pulsars with Extreme Inclinations. Universe 2021, 7, 455. https://doi.org/10.3390/universe7120455

AMA Style

Zhang F. A Magnetospheric Dichotomy for Pulsars with Extreme Inclinations. Universe. 2021; 7(12):455. https://doi.org/10.3390/universe7120455

Chicago/Turabian Style

Zhang, Fan. 2021. "A Magnetospheric Dichotomy for Pulsars with Extreme Inclinations" Universe 7, no. 12: 455. https://doi.org/10.3390/universe7120455

APA Style

Zhang, F. (2021). A Magnetospheric Dichotomy for Pulsars with Extreme Inclinations. Universe, 7(12), 455. https://doi.org/10.3390/universe7120455

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop