The MgSiO3 post-perovskite phase is the most abundant silicate phase in a super-Earth's mantle, although it only exists within the Earth's lowermost mantle. In this study, we established the thermal equation of state (EoS) of the MgSiO3 post-perovskite phase, which were determined by using both laser-heated diamond anvil cell and density-functional theoretical techniques, within a multi-megabar pressure range, corresponding to the conditions of a super-Earth's mantle. The Keane and AP2 EoS models were adopted for the first time to extract meaningful physical properties. The experimentally determined Gruneisen parameter, which is one of the thermal EoS parameters, and its volume dependence were found to be consistent with their theoretically obtained values. This reduced the previously reported discrepancy observed between experiment and theory. Both the experimental and theoretical EoS were also found to be in very good agreement for volumes at pressures and temperatures of up to 300 GPa and 5000 K, respectively. Our newly developed EoS should be applicable to a super-Earth's mantle, as well as the Earth's core-mantle boundary region.