Direct coupling would have no effect on tube loading, which is responsible for the spectra shown. If we were to re-do the article, we’d try a MOSFET cascode current source load, as well as transformer coupling, SRPP, and RC coupling.
We were surprised that cathode degeneration doesn’t work, and creates some nasty high-order terms instead. Separating the data into even and odd-order terms was essential to unscrambling the chaotic results of the spreadsheet ... a legitimate way of looking at the data, since the underlying transfer curves of odd-order (S-shaped) and even-order (C-shaped) distortion are fundamentally different.
Nowadays, we have the computer power to discover the actual shape of the input/output transfer curve, and exaggerate it enough to be visible. The regrettable drawback of FFT spectral information is that phase is usually discarded, so the underlying transfer shape cannot be found (although it can be inferred).
(What I mean by this is the phase of the distortion harmonics is important. For example, a square wave and a triangle wave look exactly the same on an FFT spectral display; the only difference is the phase of the harmonics. The magnitudes are the same. In real circuits, square waves and triangle waves are created by entirely different mechanisms, so this is important data.)