A rate equation based VCSEL thermal model is established, which can simulate the features in the non dc operating regimes, namely, small signal and transient modulating condition. The model is apt to implement in conventional SPICE like circuit simulators, including PSPICE, and used to simulate key features of VCSEL. The results compare favorably to experimental data from two devices reported in the literature.