Single-cell sequencing technologies have evolved and expanded rapidly in recent years, accompanied by the successful development of numerous computational tools. Single-cell Assay for Transposase-Accessible Chromatin sequencing (scATAC-seq) is a powerful method for genome-wide analysis of open (accessible) chromatin regions within individual cells derived from heterogeneous populations. Robust evaluation is critical for the development of scATAC-seq data analysis workflows, which are currently far from stable or mature. Notwithstanding important existing software tools, workflows remain in flux with many open questions regarding the best ways to optimize all elements of data analysis. Simulations are crucial for generating reproducible datasets with known characteristics that enable successful development, testing, and benchmarking of analysis methods. However, current scATAC-seq simulation frameworks do not fully account for the biological and technical characteristics of scATAC-seq data and are incapable of simulating population-scale scATAC-seq data with realistic genetic effects on chromatin accessibility and overall population structure. We present the scATACPop framework, a model and software tool for flexible and reproducible simulation of population-scale scATAC-seq data. We are developing a new data generative model for scATAC-seq data that enables simulation of datasets that closely resemble real scATAC-seq datasets in library size, sparsity, and chromatin accessibility signals. scATACPop can also mimic complex batch, cell group, and conditional effects between individuals from different cohorts. The inclusion of genetic effects in our simulation framework will establish scATACpop as a key tool for the simulation of population-scale scATAC-seq data, and, specifically, the development of methods for single-cell multi-omics quantitative trait locus mapping.