Using the thermal field dynamics theory to convert the thermal state into a "pure" state in doubled Fock space, we find that the average value of e f a † a under squeezed thermal state (STS) is just the generating function of Legendre polynomials. Based on this remarkable result, the normalization and photon-number distributions of m-photon added (or subtracted) STSs are conviently obtained as the Legendre polynomials. This new concise method can be expanded to the entangled case.OCIS codes: 270.0270, 270.5290. doi: 10.3788/COL201210.082701.The nonclassicality of optical fields is helpful in understanding the fundamentals of quantum optics and has many applications in quantum information processing [1] . The subtraction or addition of photons from/to traditional quantum states or Gaussian states has been proposed to generate and manipulate various nonclassical optical field [2−11] . For example, photon addition and subtraction have been experimentally demonstrated to probe quantum commutation rules [9] . Recently, photon-added (-subtracted) Gaussian states have received more attention from both experimentalists and theoreticians [12−21] , because these states exhibit numberous nonclassical properties and may provide access to a complete engineering of quantum states and fundamental quantum phenomena.Theoretically, the normalization factors of such quantum states are essential for studying their nonclassical properties. Very recently, Fan et al.[22] presented a new concise approach for normalizing m-photon-added (-subtracted) squeezed vacuum state (pure state) by constructing a generating function. However, most systems are not isolated, are immersed in a thermal reservoir, and thus, we often have no enough information to specify completely the state of a system. In such situations, the system only can be described by mixed states, such as thermal states. In addition, the squeezed thermal states (STSs) can be considered as the generalized Gaussian states.In this letter, we shall extend this case to the mixed state, i.e., by using the thermal field dynamics (TFD) theory to convert the thermal state into a "pure" state in doubled Fock space. We present a new concise method for normalizing photon-added (-subtracted) STSs (PASTSs, PSSTSs) and deriving their photon-number distributions (PNDs), which have been a major topic of studies on quantum optics and statistics. The normalization factors and PNDs were found to be related to the Legendre polynomials in compact form.We begin by briefly reviewing the properties of a thermal state. For a single mode with frequency ω in a thermal equilibrium state corresponding to absolute temperature T , the density operator iswhere n c = {exp[ ω/(kT )] − 1} −1 is the average photon number of the thermal state ρ th and k is the Boltzmann's constant. |n = a †n / √ n! |0 and the normally ordering form of vacuum projector |0 0| = : exp(−a † a) : (the symbol : : denotes normal ordering). One can express Eq. (1) aswhere in the last step, the operator identity exp λa † a = : exp e λ − 1 ...