First-Principles Polar-Cap Currents in Multipolar Pulsar Magnetospheres

TL;DR

The paper confronts the bias introduced by dipole-only hotspot prescriptions in X-ray pulse-profile modeling of MSPs by deriving analytic expressions for the field-aligned current invariant $\Lambda(\alpha,\beta)$ in mixed dipole–quadrupole magnetospheres within stationary force-free electrodynamics. Using a matched-asymptotic approach with Euler potentials, it develops a cap-interior harmonic closure and provides closed-form $\Lambda$ for oblique dipoles, axisymmetric quadrupoles, and general quadrupoles, including a dipole–quadrupole mixing prescription via a far-zone parameter $\eta$. These analytic results yield self-consistent polar-cap current distributions and hotspot heating maps, which, when fed into relativistic pulse-profile modeling with BB and NSX atmosphere beaming, reveal substantial differences in light curves and energy-resolved profiles compared with dipole-only predictions, especially near pulse peaks (up to ~30% in the NICER band under NSX beaming). The framework enables physically grounded inference of multipolar magnetic geometry in MSPs and provides a foundation for incorporating quadrupole effects into forward models and Bayesian analyses, thereby reducing systematic biases in neutron-star parameter estimation and enabling multiwavelength magnetospheric constraints.

Abstract

X-ray pulse-profile modeling of millisecond pulsars offers a direct route to measuring neutron star masses and radii, thereby constraining the dense-matter equation of state. However, standard analyses typically rely on \emph{ad hoc} hotspot parameterizations rather than self-consistent physical models. While connecting surface heating directly to the magnetospheric geometry provides a more natural physical pathway, computing global magnetospheric solutions is too computationally expensive to perform on-the-fly during parameter inference. In this work, we bridge this gap by deriving fully analytic, first-principles expressions for surface return currents in mixed dipole--quadrupole magnetospheres. Working within force-free electrodynamics, we generalize the field-aligned current invariant $Λ$, the crucial scalar that maps the far-zone magnetic structure to the near-zone heating rate, from the standard dipole approximation to arbitrary quadrupolar configurations. We demonstrate that even when the quadrupole component is sub-dominant in the far zone (the mixing regime), using a dipole-based heating prescription fails to capture the significant enhancement or suppression of the return-current density on the polar cap. Our consistent quadrupole-aware framework reveals that these multipolar currents redistribute the surface heating, leading to systematic discrepancies in predicted pulse profiles that are amplified by atmosphere beaming and can reach $\sim 30\%$ near pulse peaks. These results provide a rigorous analytic foundation for mapping global magnetic geometry to surface heating in multipolar magnetospheres, enabling physically consistent inference beyond the idealized dipole approximation.

Figures (6)

• Figure 1: Mapping between the near-zone quadrupole amplitude $Q_\alpha$ and the far-zone dipole--quadrupole mixing parameter $\eta \equiv \left(B_Q/B_D\right)_{\rm RMS}(R_{\rm LC})$ (evaluated using equation \ref{['eq:eta-Qalpha-lockhart']}). Left: $\eta(Q_\alpha)$ at fixed stellar radius $R=12\,{\rm km}$ for three spin frequencies $\nu=200,300,400\,{\rm Hz}$. Right: $\eta(Q_\alpha)$ at fixed $\nu=400\,{\rm Hz}$ for three radii $R=10,12,14\,{\rm km}$. The representative sequence used in Figures \ref{['fig:jsq-maps']}--\ref{['fig:erpp-nsx']} corresponds to $Q_\alpha=0.50,1.00,1.50$ (i.e., $\eta\simeq 0.041,\,0.082,\,0.123$ for $\nu=400\,{\rm Hz}$ and $R=12\,{\rm km}$).
• Figure 2: Surface maps of the 4-current invariant $J^2$ in magnetic spherical coordinates $(\theta,\phi)$ for three quadrupole amplitudes (columns: $Q_\alpha=0.50,1.00,1.50$). Top row: dipole-only flux-conservation functional $\Lambda$. Middle row: quadrupole-aware (mixed dipole--quadrupole ) $\Lambda$. Bottom row: residual $\Delta J^2 \equiv J^2_{\rm mix}-J^2_{\rm dip}$. Colors use a symmetric-logarithmic scale and are normalized by the maximum absolute value in each panel; red indicates $J^2>0$ (spacelike 4-current) and blue indicates $J^2<0$ (timelike 4-current).
• Figure 3: Surface effective-temperature patterns inferred from the return-current prescription $T\propto |{\bf j}|^{1/4}$ (overall normalization arbitrary), shown in magnetic coordinates $(\theta,\phi)$ for the same $Q_\alpha$ sequence as Fig. \ref{['fig:jsq-maps']}. Top row: dipole-only $\Lambda$. Middle row: mixed dipole--quadrupole $\Lambda$. Bottom row: relative difference $\Delta T \equiv (T_{\rm mix}-T_{\rm dip})/T_{\rm dip}$ (symmetric-logarithmic color scale). In the top and middle rows, each panel is normalized by its own maximum temperature, i.e. colors show $T/T_{\max}$.
• Figure 4: Normalized bolometric light curves integrated over $0.1$--$5.2125\,{\rm keV}$ (NICER band) for three quadrupole amplitudes (columns: $Q_\alpha=0.50,1.00,1.50$), computed for both isotropic blackbody emission (BB) and a beamed NSX atmosphere. Top row: light curves for the dipole-only and mixed dipole--quadrupole $\Lambda$ prescriptions (legend), plotted versus rotational phase. Within each panel, the flux is normalized by the maximum among the four curves shown to show the relative magnitude between different light curves. Bottom row: residuals $\Delta F \equiv F_{\rm mix}-F_{\rm dip}$ shown separately for BB and NSX. All light-curve calculations use the fiducial numerical/geometry setup listed in Table 1.
• Figure 5: Energy-resolved pulse profiles for blackbody (BB) emission, shown as phase--energy maps for the same $Q_\alpha$ sequence (columns: $Q_\alpha=0.50,1.00,1.50$). Top row: dipole-only $\Lambda$. Middle row: mixed dipole--quadrupole $\Lambda$. Bottom row: residual (mix$-$dip), with a zero-centered symmetric-logarithmic color scale. Horizontal axes show rotational phase, vertical axes show photon energy (keV). The color scale shows the photon counts on each energy-phase bin.