Folded NICER XTI profiles from Figure 2 but grouped in 256 (for PSRs J0437−4715 and J0030+0451) or 128 phase bins (for PSRs J1231−1411 and J2124−3358). The solid red lines show the best fit with a model of the profile constructed from the empirical Fourier coefficients, given by Equation (1), of the set of photon phases. The dashed lines show the sinusoids corresponding to the first four harmonic components of the fit (blue, orange, purple, and cyan for the first, second, third, and fourth harmonic, respectively). The bottom panel for each pulsar shows the residuals from the fit.