the HB simulator rarely converges with BJT models. It easily converges with diode and JFET. I tried with several settings and number of harmonic but the convergence with BJT is very hard to obtain.
This is a known problems. It appears due to the transit times, namely the parameters “Tf” and “Tr” of the BJT. With zero values here, it should work.
It would be great, if someone can help to solve this problem, i.e. check if it works in simulators like XYCE (https://xyce.sandia.gov/), fREEDA (https://www.freeda.org/) or similar. Or maybe someone finds a publication about this topic, knows a friend who has a solution etc.