2018ApJ...867...68W

Utilizing all the 16 yr RXTE observations, we analyze the X-ray spectra of 32 TeV blazars and perform a systematic investigation of X-ray spectral variability for the five brightest sources during their major flares that lasted several days. We obtain photon spectral index (α), flux and synchrotron radiation peak energy (E_p) from empirical spectral fitting, and electron spectral index (p) from theoretical synchrotron radiation modeling. We find that both α and p generally display a harder-when-brighter trend, confirming the results of many previous works. Furthermore, we confirm and strengthen the result that p must vary in order to explain the observed X-ray spectral variability during flares, which would have useful implications for interpreting the associated higher-energy spectral variability. We see apparent electron spectral hysteresis in many but not all p-flux plots that takes a form of "loop" or oblique "8." We obtain a tight p-hardness ratio (HR) relation and a tighter p-α relation using spectra of flaring periods, both of which are also applicable to stacked data of quiescent periods. We demonstrate that these two empirical relations can be used efficiently to estimate p from HR or α that is readily achieved. Finally, we find that, when considering TeV blazars as a whole, α and X-ray luminosity are positively correlated, E_p is negatively correlated with p and α, and E_p is positively correlated with HR; all these correlations are in line with the blazar sequence. However, after correcting for the Doppler boosting effect, α and intrinsic X-ray luminosity follow an anticorrelation.