<p>In this work, we develop a new implicit Monte Carlo method to treat the time evolution of the coupled radiation, electron and ion radiative transfer equations. This method deals with all the coupling of electron and ion energies in the Monte Carlo simulation. It is shown that the method has the property of total energy conservation, and is asymptotic preserving in the limit that the electron-ion coupling coefficient tends to zero or infinity. Moreover, it involves nonlinear iterations for the electron and ion temperatures, which allows for an implicit estimate of the electron and ion temperatures dependent quantities. We demonstrate that when the electron temperature satisfies the maximum principle, then the radiation and ion temperatures will satisfy the maximum principle automatically. Some numerical simulation results are given to demonstrate the effectiveness of the proposed method.</p>