Do not waste your time to simulate THD with these models. The nonlinearity is not modeled in these simple models. If it is modeled as good as the chips, anyone understand it can reverse engineer it with ease. The model needs every transistor and parasitic components to fully model the nonlinearity, more than full schematic, it is likely kept as secret. These models can only do small signal AC response, nothing else.
I do not think we talked about THD measurements, all Spice models are a crude approximation, for both discrete components, even more for analog IC, it's just a way to have a general idea of how it tends to change the behavior in amplitude, gain, frequency, phase,
or estimate the behavior of the power supply according to the current delivered and its internal impedance etc..
I agree, the only reliable way is real measure, especially if you want to be sure that the circuit is stable, but not everyone has a suitable equipment, at this point, the only way to have a generic idea is the simulation.
If not, what is the usefulness of the manufacturers to provide spice models?