Discharges exhibit a bit of negative resistance (higher current leads to higher cathode spot temperature, so lower cathode drop),
but also the dynamics of the ionization may cause the system to become unstable and oscillate: In an equilibrium state the free charge avalanche multiplication is suppossed to generate equal amount of free electron/ion pairs, as they recombine.
But both effects have some time delay, so cause a phase lag (when expressed as how they respond to a small superimposed sinewave like changes). So we have a system with two phase lags in the loop (one is voltage so multiplication rate -> amount of charges generated means an integration, so 90deg phase lag, other is gradual recombination, forming a kind of exponential decay, so could be described as a low pass filter with some time constant, so again up to 90deg phase lag), which could lead to phase lag approaching 180deg.
That means in the loop (voltage -> charge multiplication + charge decay -> amount of charges -> conductivity -> current -> voltage) which is has normally 180deg phase shift (it is a negative feedback, at least per static behavior), the extra 180deg phase shift means the originally negative feedback becomes positive
(see Bode stability criteria) when dynamic behavior taken into account. The remaining few degrees lag to the full 360deg may come from either the electrical response (e.g. a capacitance adding extra lag current -> voltage) or even within the tube when taking the finite speed of sound waves within the plasma (how fast the electrons and ions distribute, so causing extra delay), you may loose all of the remaining phase margin so get an oscillator.
And for the power factor:
The PF=~0.9 comes just from the fact at mains frequency AC supply, the current is nearly a sinewave, but the voltage is nearly a square wave, so they have quite a mismatch. And if you calculate just an ideal sinewave vs an ideal square wave, you get PF = (2/pi) / sqrt(1/2) = 0.900316...