Cytometry by time-of-flight (CyTOF) is a proteomic analysis technique capable of providing information on millions of cells. Along with technological development, researchers have devised many data analysis tools to offer deeper insights into the data at the single-cell level. However, objective evaluations of these methods remain absent. In this paper, we develop Cytomulate: an accurate, efficient, and well-documented simulation model and algorithm for reproducible and accurate simulation of CyTOF data, which could serve as a foundation for future development and evaluation of new methods. With its two simulation modes: creation mode and emulation mode, Cytomulate is capable of generating simulated datasets that capture various characteristics of CyTOF data under different circumstances. We also provide the users with methods to simulate batch effects, temporal effects, and cell differentiation paths. Empirical experiments suggest that Cytomulate is superior in learning the overall data distribution than other single-cell RNA-seq-oriented methods such as Splatter.