The thermodynamic properties of insensitive high explosive 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) based polymer bonded explosives (PBXs) are investigated by using first principle calculation and molecular dynamics simulation. The results include the phonon dispersion relations, interface thermal conductances, and thermal conductivities of TATB based PBXs. Both TATB and PVDF structures are optimized, in which the optimized lattice constants accord with previous results. The phonon dispersion relation of TATB and PVDF are calculated based on lattice dynamics. All interatomic force constants are calculated by the finite displacement method (numerical derivatives from perturbed supercells). The calculated phonon dispersion relation of TATB and heat capacity are in general agreement with experimental and theoretical results. The imaginary frequencies are observed in both TATB and PVDF dispersion relation. The imaginary frequencies are mainly due to the smaller calculated supercell size and temperature effect. The phonon mode of TATB and PVDF are assigned at Γ point. Based on the calculated phonon dispersion, some information including heat capacity, phonon density of states and phonon mode assignment is derived. The TATB possesses 144 phonon modes including 3 acoustic-phonon modes and 141 optical phonon modes. The anylized phonon mode of TATB shows that -NO₂ dominates the phonon DOS in low frequency zone, phenyl rings dominate in middle frequency zone and -NH₂ dominates in high frequency zone. By analyzing the phonon density of states and capacity, both TATB and PVDF imply that low-frequency vibration dominates the thermal conductivity. The thermal conductivity is determined for TATB by using the equlibrium molecular dynamics method and an established TATB force field. The TATB model is built with 2880 atoms. The structure of TATB is optimized by using molecular mechanics, then this system is relaxed by using a Nose-Hoover thermostat and barostat with a damping factor of 50 fs cin time steps of 0.1 fs. The calcultated thermal conductivity at room temperature shows good agreement with experimental result. The interface thermal conductance of TATB-PVDF is calculated by using a diffusive mismatch model. The interface thermal transport still follows Fourier’s law of heat conduction, and ballistic thermal transport mechanism is not involved. By using the above results, the thermal conductivity of mixture TATB-PVDF system is analized with a simple series model. The particle size smaller than 100 nm significantly suppresses the mixture system thermal conductivity.