The first part of this work is devoted to the study of higher derivatives of pressure functions of Hölder potentials on shift spaces with finitely many symbols. By describing the derivatives of pressure functions via the central limit theorem for the associated random processes, we discover some rigid relationships between derivatives of various orders. The rigidity imposes obstructions on fitting candidate convex analytic functions by pressure functions of Hölder potentials globally, which answers a question of Kucherenko and Quas. In the second part of the work, we consider fitting candidate analytic germs by pressure functions of locally constant potentials. We prove that all 1-level candidate germs can be realised by pressures of some locally constant potentials, as long as the number of symbols in the symbolic set is large enough. There are also some results on fitting 2-level germs by pressures of locally constant potentials obtained in the work.
