Continuum study of various susceptibilities within thermal QED$_3$
  In this paper, the relations of four different susceptibilities (i.e., the chiral susceptibility, the fermion number susceptibility, the thermal susceptibility and the staggered spin susceptibility) are investigated both in and beyond the chiral limit. To this end, we numerically solve the finite temperature version of the truncated Dyson-Schwinger equations for fermion and boson propagator. It is found that, in the chiral limit, the four susceptibilities give the same critical temperature and signal a typical second order phase transition. But the situation changes beyond the chiral limit: the critical temperatures from the chiral and the thermal susceptibilities are different, which shows that to define a critical region instead of an exclusive point for crossover might be a more suitable choice; meanwhile, both the fermion number and the staggered spin susceptibilities have no singular behaviors any more, this may mean that they are no longer available to describe the crossover properties of the system. 
Pei-lin Yin, Yuan-mei Shi, Zhu-fang Cui, Hong-tao Feng, and Hong-shi   Zong
hep-ph