Abstract
The average pulse shape of a pulsar is typically stable over decadal timescales, enabling estimation of pulse times of arrival to better than a small fraction of the pulse width using matched filtering techniques. However, in North American Nanohertz Observatory for Gravitational Waves (NANOGrav) observations of PSR J1713+0747, three discrete timing events that depart from the prevailing timing model have been seen in the last 20 yr. All three correspond to morphological changes in pulse shape. Using principal component analysis, we analyze the pulse profiles of nine NANOGrav pulsars, including seven with profiles from the 15 yr dataset and two with additional profiles from the forthcoming 20 yr dataset. We recover the three known pulse shape change events in PSR J1713+0747 and another previously known event in PSR J1643 - 1224. We implement a ranking metric for candidate events and address four highly ranked candidates in this nine-pulsar sample. We also recover known slow pulse shape variations in PSR J1643 - 1224, PSR J1903+0327, and PSR B1937+21 and report an unexpected recurrence after 10 yr of one such variation in PSR B1937+21.