We present first principles calculations of the thermal and thermal transport properties of Bi2Te3 that combine an ab initio molecular dynamics (AIMD) approach to calculate interatomic force constants (IFCs) along with a full iterative solution of the Peierls-Boltzmann transport equation for phonons. The newly developed AIMD approach allows determination of harmonic and anharmonic interatomic forces at each temperature, which is particularly appropriate for highly anharmonic materials such as Bi2Te3. The calculated phonon dispersions, heat capacity, and thermal expansion coefficient are found to be in good agreement with measured data. .The lattice thermal conductivity, κ l , calculated using the AIMD approach nicely matches measured values, showing better agreement than the κ l obtained using temperature independent IFCs. A significant contribution to κ l from optic phonon modes is found. Already at room temperature, the phonon lineshapes show a notable broadening and onset of satellite peaks reflecting the underlying strong anharmonicity.