Searches for gravitational waves from known pulsars at two harmonics in the second and third LIGO-Virgo observing runs

We present a targeted search for continuous gravitational waves (GWs) from 236 pulsars using data from the third observing run of LIGO and Virgo (O3) combined with data from the second observing run (O2). Searches were for emission from the l = m = 2 mass quadrupole mode with a frequency at only twice the pulsar rotation frequency (single harmonic) and the l = 2, m = 1, 2 modes with a frequency of both once and twice the rotation frequency (dual harmonic). No evidence of GWs was found, so we present 95% credible upper limits on the strain amplitudes h (0) for the single-harmonic search along with limits on the pulsars' mass quadrupole moments Q (22) and ellipticities epsilon. Of the pulsars studied, 23 have strain amplitudes that are lower than the limits calculated from their electromagnetically measured spin-down rates. These pulsars include the millisecond pulsars J0437-4715 and J0711-6830, which have spin-down ratios of 0.87 and 0.57, respectively. For nine pulsars, their spin-down limits have been surpassed for the first time. For the Crab and Vela pulsars, our limits are factors of similar to 100 and similar to 20 more constraining than their spin-down limits, respectively. For the dual-harmonic searches, new limits are placed on the strain amplitudes C (21) and C (22). For 23 pulsars, we also present limits on the emission amplitude assuming dipole radiation as predicted by Brans-Dicke theory.