In high-fluence laser-induced incandescence (LII), current LII models significantly overpredict the soot nanoparticle temperature decay rate compared to that inferred from two-color pyrometry at the first 100 ns after the peak laser pulse in atmospheric pressure flames. One possible cause is the back flow of sublimed species, which to date has been neglected in LII modeling. In this study, the transient direct simulation Monte Carlo (DSMC) method has been used, for the first time, to calculate the temperature evolution of soot particles, taking into account recondensation of sublimed species. In this algorithm, the physical time is discretized into a number of time steps called ensemble time steps, and the heat flux is calculated by performing several DSMC runs in each ensemble time step before proceeding to the next ensemble time step until the variance reaches an acceptable value. This heat flux is then used to update the nanoparticle temperature over the ensemble time step. Using the new algorithm, the temperature evolution of the particle can be predicted by the DSMC code, which is an improvement to previous DSMC simulations in which predetermined temperature decay curves must be prescribed. The results show that recondensation of sublimed species on the originating nanoparticle is not significant. Although accounting for condensation of sublimed species originating from neighboring soot particles enhances the role of recondensation of sublimed species in slowing down the soot particle temperature decay, it is still not sufficient to be considered as a plausible cause for the discrepancy between modeled soot temperature and the two-color pyrometry measured one in high-fluence LII.