Thus, the basic model for smoking dynamics is given by Eqs. (1)-(4), (see flowchart diagram in Fig. 1a). It is an extension of the model in [4,10], by incorporating a class (and dynamics) of smokers who temporarily quit smoking (Qt(t)). Note that this model does not include smoking-related mortality.