Crack propagation due to thermal expansion is one of the multi-physics problems which be a concern to many researchers. Therefore, we study thermal fracturing using the phase field model (or, PFM). Here, the new models to study thermal fracturing in linear thermoelasticity solids are proposed through PFM. Herein, the damage evolution equation is directly derived through the microforce balance approach or the Gurtin concept, while the heat evolution equation is derived through the first principle of thermodynamics. Furthermore, the thermodynamic consistency of the model is shown by Clausius-Duhem inequality. In particular, the Gurtin concept and the first principle of thermodynamics follow the entropy, internal energy, Helmholtz free energy, and energy dissipation functions which are based on the Biot of thermoelasticity model and the Ambrosio-Tortorelli regularization. Since our models are based on the microforce balance approach, we also derive the PFM for crack propagation which was proposed by Kimura and Takaishi via this approach. In the present study, we validate our proposed PFMs through several numerical experiments. Herein, we solve the numerical experiments using the adaptive finite element method. From these, good agreements are achieved between our models and the previous model.