skysurvey: a pure python package to simulate the transient sky
Pith reviewed 2026-06-29 20:36 UTC · model grok-4.3
The pith
skysurvey structures simulations of astronomical transients around Target, Survey and DataSet objects plus modeldag to enable complex population modeling with minimal code.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The skysurvey package relies on three core objects—a Target that models how an astrophysical target exists in nature, a Survey that specifies how the sky has been observed, and a DataSet that combines these two to generate data as they would have been acquired—along with the modeldag package for parameter modeling; this structure allows users to implement complex populations such as Type Ia Supernovae with varying beta term and to replicate the ZTF SNe Ia DR2 dataset.
What carries the argument
Three core objects Target, Survey and DataSet together with the modeldag structure for parameter modeling.
If this is right
- Complex populations such as Type Ia Supernovae with varying color-brightness beta term can be implemented straightforwardly.
- The package can replicate the rate and redshift distribution of the ZTF SNe Ia DR2 dataset.
- Simulations can account for selection effects and correct biases in analyses.
- The package supports use in simulation-based inference.
- The package is ready for general usage in scientific publications.
Where Pith is reading between the lines
- The same three-object structure could be applied to transients other than supernovae.
- Users might combine skysurvey outputs directly with existing machine-learning detection pipelines.
- The modeldag component could be extended to include additional physical parameters without rewriting core simulation logic.
Load-bearing premise
The three core objects Target, Survey, DataSet together with modeldag suffice to implement arbitrary complex astrophysical population models with minimal additional user code.
What would settle it
A concrete complex population model that requires substantial additional user code beyond the three core objects and modeldag to implement.
Figures
read the original abstract
Accurate simulation of astronomical observations is a critical element for any modern analyses, be it to measure event rates, analyses population properties, validate or train pipelines, account for selection effects, or correct biases. We present a novel pure python package named skysurvey made to enable the user to quickly simulate astrophysical transients as observed by a survey. The package is structured to make the implementation of any complex population modeling fast and easy. The skysurvey package relies on three core objects: a Target, that models how an astrophysical target exists in nature, a Survey, that specifies how the sky has been observed and, a DataSet that combine these two to generate data as they would have been acquired. In addition, we present a side stand-alone package named modeldag that contains the core structure that simplifies the parameter modeling. We present in this paper how skysurvey is structured and we clearly illustrate how the code can straightforwardly be used to simulate complex populations, such as Type Ia Supernovae with varying color-brightness $\beta$ term. We also illustrate how the package can be made to replicate the rate and redshift distribution of the ZTF SNe Ia DR2 dataset. The skysurvey package, already used in recent scientific publications, is now ready for general usage and paves the way for future use of simulations such as simulation based inference.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript describes skysurvey, a pure Python package for simulating astronomical transients as observed by surveys. It introduces three core objects (Target for astrophysical populations, Survey for observational specifications, DataSet for combining them) plus the standalone modeldag package for parameter DAGs. The central claim is that this structure enables fast and easy implementation of arbitrary complex population models, illustrated via a Type Ia supernova example with varying color-brightness eta and a replication of the ZTF SNe Ia DR2 rate/redshift distribution. The package is stated to be already used in publications and suitable for simulation-based inference.
Significance. If the core objects and modeldag indeed allow complex models with minimal additional code, the package would be a useful contribution to transient astronomy for rate measurements, bias corrections, and pipeline validation. The pure-Python implementation, open availability, and existing scientific usage are positive factors that could facilitate reproducible simulations.
major comments (3)
- [Abstract] Abstract and the section on complex population modeling: the assertion that Target/Survey/DataSet + modeldag suffice to implement 'any complex population modeling' with 'minimal additional user code' is load-bearing for the paper's main claim but is supported only by the two provided illustrative examples (varying-eta Type Ia and ZTF DR2 replication); no test cases involving multi-population correlations, time-varying selection functions, or non-standard DAG extensions are shown, leaving the generality unverified.
- [ZTF replication section] ZTF DR2 replication illustration: the claim of replicating the observed rate and redshift distribution is presented without any quantitative validation (e.g., Kolmogorov-Smirnov statistic, χ^{2}, or residual plots comparing simulated vs. DR2 data), which is required to assess whether the simulation is accurate enough for the scientific applications listed in the abstract.
- [Core objects description] Section describing the core objects: no code-complexity metrics, line-count comparisons, or user-effort benchmarks are supplied against direct Python implementations or alternative frameworks, which is needed to substantiate the 'fast and easy' claim for arbitrary models.
minor comments (2)
- The manuscript would benefit from embedding short, self-contained code snippets for the two main examples directly in the text (rather than only pointing to external documentation) to improve immediate reproducibility.
- A brief comparison table listing skysurvey against other transient simulation tools (e.g., SNANA, simsurvey) would help readers place the new package in context.
Simulated Author's Rebuttal
We thank the referee for their constructive comments. We address each major comment below, making revisions where the points identify clear gaps in the presented evidence.
read point-by-point responses
-
Referee: [Abstract] Abstract and the section on complex population modeling: the assertion that Target/Survey/DataSet + modeldag suffice to implement 'any complex population modeling' with 'minimal additional user code' is load-bearing for the paper's main claim but is supported only by the two provided illustrative examples (varying-β Type Ia and ZTF DR2 replication); no test cases involving multi-population correlations, time-varying selection functions, or non-standard DAG extensions are shown, leaving the generality unverified.
Authors: We agree that the generality claim would be stronger with explicit discussion of additional use cases. The two examples already cover intra-population parameter variation and replication of an observed distribution, both of which rely on the modular Target/Survey/DataSet structure and modeldag. In the revision we will add a concise paragraph in the discussion section outlining how the same objects extend to multi-population correlations and time-varying selection functions, without introducing new full test cases. revision: partial
-
Referee: [ZTF replication section] ZTF DR2 replication illustration: the claim of replicating the observed rate and redshift distribution is presented without any quantitative validation (e.g., Kolmogorov-Smirnov statistic, χ², or residual plots comparing simulated vs. DR2 data), which is required to assess whether the simulation is accurate enough for the scientific applications listed in the abstract.
Authors: We accept that quantitative metrics are needed. The revised manuscript will include a Kolmogorov-Smirnov test and a residual plot comparing the simulated and DR2 redshift distributions in the ZTF replication section. revision: yes
-
Referee: [Core objects description] Section describing the core objects: no code-complexity metrics, line-count comparisons, or user-effort benchmarks are supplied against direct Python implementations or alternative frameworks, which is needed to substantiate the 'fast and easy' claim for arbitrary models.
Authors: We will add a short comparison of lines of user code required for the Type Ia example versus an equivalent direct implementation, placed in the core-objects section of the revision. revision: yes
Circularity Check
No circularity: software description with no derivation chain
full rationale
This is a software package paper describing skysurvey and modeldag. It defines three core objects (Target, Survey, DataSet) and illustrates their use via code examples for Type Ia supernovae and ZTF replication. No equations, predictions, fitted parameters, or first-principles derivations are present. The central claim that the objects enable complex modeling is a statement about code structure and API design, supported directly by the presented implementation and examples rather than any self-referential reduction or self-citation load-bearing argument. No patterns from the enumerated circularity kinds apply.
Axiom & Free-Parameter Ledger
Reference graph
Works this paper leans on
-
[1]
Amenouche, M., Rosnet, P., Smith, M., et al. 2025, A&A, 694, A3. doi:10.1051/0004-6361/202452134 Andreoni, I., Coughlin, M. W., Perley, D. A., et al. 2022, Nature, 612, 7940,
-
[2]
doi:10.1038/s41586-022-05465-8 Barbary, K., Barclay, T., Biswas, R., et al. 2016, Astrophysics Source Code Li- brary. ascl:1611.017 Barjou-Delayre, C., Rosnet, P., Raboux, C., et al. 2026, arXiv:2601.11139. doi:10.48550/arXiv.2601.11139 Bellm, E. C., Kulkarni, S. R., Graham, M. J., et al. 2019, PASP, 131, 995, 018002. doi:10.1088/1538-3873/aaecbe Bernstei...
-
[3]
doi:10.1088/0004-637X/753/2/152 Betoule, M., Kessler, R., Guy, J., et al. 2014, A&A, 568, A22. doi:10.1051/0004- 6361/201423413 Boyd, B. M., Mandel, K. S., Grayling, M., et al. 2026, arXiv:2603.11165. doi:10.48550/arXiv.2603.11165 Brout, D., & Scolnic, D. 2021, The Astrophysical Journal, 909,
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1088/0004-637x/753/2/152 2014
-
[4]
Brout, D., Scolnic, D., Popovic, B., et al
doi:10.3847/1538-4357/abd69b. Brout, D., Scolnic, D., Popovic, B., et al. 2022, ApJ, 938, 2,
-
[5]
The Pantheon+ Analysis: Cos- mological Constraints
doi:10.3847/1538-4357/ac8e04 Bulla, M. 2019, MNRAS, 489, 4,
-
[6]
doi:10.1093/mnras/stz2495 Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 345,
-
[7]
LightCurveLynx: Forward Modeling of Time-Domain Surveys with Application to ZTF SN Ia DR2
doi:10.1086/167900 Carreres, B., Bautista, J. E., Feinstein, F., et al. 2023, A&A, 674, A197. doi:10.1051/0004-6361/202346173 Chotard, N., Gangler, E., Aldering, G., et al. 2011, A&A, 529, L4. doi:10.1051/0004-6361/201116723 Cranmer, K., Brehmer, J., & Louppe, G. 2020, Proceedings of the National Academy of Science, 117, 48, 30055. doi:10.1073/pnas.191278...
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1086/167900 2023
-
[8]
doi:10.1088/1475-7516/2019/10/005 Freeburn, J., Andreoni, I., de Soto, K. M., et al. 2025, ApJ, 994, 1, L24. doi:10.3847/2041-8213/ae1cba Frohmaier, C., Sullivan, M., Nugent, P. E., et al. 2019, MNRAS, 486, 2,
-
[9]
Illuminating the Local Universe: Large-Scale Structure from ZTF Type Ia Supernovae
doi:10.1093/mnras/stz807 Gilles Lordet, A., Goobar, A., Jasche, J., et al. 2026, arXiv:2604.12714. Gillies, S., van der Wel, C., Van den Bossche, J., et al. 2025, Zenodo. doi:10.5281/zenodo.5597138. Ginolin, M., Rigault, M., Smith, M., et al. 2025, A&A, 695, A140. doi:10.1051/0004-6361/202450378 Ginolin, M., Rigault, M., Copin, Y ., et al. 2025, A&A, 694,...
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1093/mnras/stz807 2026
-
[10]
doi:10.1051/0004-6361:20020930 Green, G. M. 2018, The Journal of Open Source Software, 3, 26,
-
[11]
doi:10.21105/joss.00695 Guy, J., Sullivan, M., Conley, A., et al. 2010, A&A, 523, A7. doi:10.1051/0004- 6361/201014468 Hinton, S. & Brout, D. 2020, The Journal of Open Source Software, 5, 47,
-
[12]
doi:10.21105/joss.02122 Hložek, R., Malz, A. I., Ponder, K. A., et al. 2023, ApJS, 267, 2,
-
[13]
doi:10.3847/1538-4365/accd6a Ho, M., Bartlett, D. J., Chartier, N., et al. 2024, The Open Journal of Astro- physics, 7,
-
[14]
doi:10.33232/001c.120559 Hounsell, R., Scolnic, D., Foley, R. J., et al. 2018, ApJ, 867, 1,
-
[15]
doi:10.3847/1538-4357/aac08b Ivezi´c, Ž., Kahn, S. M., Tyson, J. A., et al. 2019, ApJ, 873, 2,
-
[16]
doi:10.3847/1538-4357/ab042c Ishiyama, T., Fukushige, T., & Makino, J. 2009, PASJ, 61,
-
[17]
doi:10.1093/pasj/61.6.1319 Jasche, J. & Lavaux, G. 2019, A&A, 625, A64. doi:10.1051/0004- 6361/201833710 Karambelkar, V . R., Kasliwal, M. M., Blagorodnova, N., et al. 2023, ApJ, 948, 2,
-
[18]
doi:10.3847/1538-4357/acc2b9 Karchev, K., Trotta, R., & Weniger, C. 2023, MNRAS, 520, 1,
-
[19]
doi:10.1093/mnras/stac3785 Kelsey, L., Sullivan, M., Smith, M., et al. 2021, MNRAS, 501, 4,
-
[20]
doi:10.1093/mnras/staa3924 Kessler, R., Bernstein, J. P., Cinabro, D., et al. 2009, PASP, 121, 883,
-
[21]
doi:10.1086/605984 Kessler, R., Bassett, B., Belov, P., et al. 2010, PASP, 122, 898,
-
[22]
doi:10.1086/657607 Kessler, R. & Scolnic, D. 2017, ApJ, 836, 1,
-
[23]
Spinning Primordial Black Holes and Scalar In- duced Gravitational Waves from Single Field Inflation
doi:10.3847/1538- 4357/836/1/56 Kessler, R., Narayan, G., Avelino, A., et al. 2019, PASP, 131, 1003, 094501. doi:10.1088/1538-3873/ab26f1 Leaman, J., Li, W., Chornock, R., et al. 2011, MNRAS, 412, 3,
-
[24]
doi:10.1111/j.1365-2966.2011.18158.x Li, W., Chornock, R., Leaman, J., et al. 2011, MNRAS, 412, 3,
-
[25]
doi:10.1111/j.1365-2966.2011.18162.x Masci, F. J., Laher, R. R., Rusholme, B., et al. 2019, PASP, 131, 995, 018003. doi:10.1088/1538-3873/aae8ac Nissanke, S., Kasliwal, M., & Georgieva, A. 2013, ApJ, 767, 2,
-
[26]
doi:10.1088/0004-637X/767/2/124 Pain, R., Fabbro, S., Sullivan, M., et al. 2002, ApJ, 577, 1,
-
[27]
A., Fremling, C., Sollerman, J., et al
doi:10.1086/342129 Perley, D. A., Fremling, C., Sollerman, J., et al. 2020, ApJ, 904,
-
[28]
doi:10.3847/1538-4357/abbd98 Planck Collaboration, Aghanim, N., Ashdown, M., et al. 2016, A&A, 596, A109. doi:10.1051/0004-6361/201629022 Planck Collaboration, Aghanim, N., Akrami, Y ., et al. 2020, A&A, 641, A6. doi:10.1051/0004-6361/201833910 Popovic, B., Brout, D., Kessler, R., et al. 2021, ApJ, 913, 1,
-
[29]
doi:10.3847/1538-4357/abf14f Popovic, B., Brout, D., Kessler, R., et al. 2023, ApJ, 945, 1,
-
[30]
doi:10.3847/1538-4357/aca273 Popovic, B., Shah, P., Kenworthy, W. D., et al. 2025, arXiv:2511.07517. doi:10.48550/arXiv.2511.07517 Riess, A. G., Press, W. H., & Kirshner, R. P. 1996, ApJ, 473,
work page internal anchor Pith review Pith/arXiv arXiv doi:10.3847/1538-4357/aca273 2025
-
[31]
G., Casertano, S., Yuan, W., et al
doi:10.1086/178129 Riess, A. G., Casertano, S., Yuan, W., et al. 2019, ApJ, 876, 1,
-
[32]
doi:10.3847/1538-4357/ab1422 Rigault, M., Copin, Y ., Aldering, G., et al. 2013, A&A, 560, A66. doi:10.1051/0004-6361/201322104 Rigault, M., Brinnel, V ., Aldering, G., et al. 2020, A&A, 644, A176. doi:10.1051/0004-6361/201730404 Rigault, M., Smith, M., Goobar, A., et al. 2025, A&A, 694, A1. doi:10.1051/0004-6361/202450388 Rigault, M., Smith, M., Regnault...
work page internal anchor Pith review doi:10.3847/1538-4357/ab1422 2013
-
[33]
doi:10.3847/1538-4357/ade1d6 Sagués Carracedo, A., Bulla, M., Feindt, U., et al. 2021, MNRAS, 504, 1,
-
[34]
doi:10.1093/mnras/stab872 Sagués Carracedo, A., Goobar, A., Mörtsell, E., et al. 2024, arXiv:2406.00052. doi:10.48550/arXiv.2406.00052 Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 2,
-
[35]
doi:10.1086/305772 Scolnic, D. & Kessler, R. 2016, ApJ, 822, 2, L35. doi:10.3847/2041- 8205/822/2/L35 Scolnic, D. M., Jones, D. O., Rest, A., et al. 2018, ApJ, 859, 2,
work page internal anchor Pith review doi:10.1086/305772 2016
-
[36]
doi:10.3847/1538-4357/aab9bb Shvartzvald, Y ., Waxman, E., Gal-Yam, A., et al. 2024, ApJ, 964, 1,
-
[37]
2024, ApJ, 964, 74, doi: 10.3847/1538-4357/ad2704
doi:10.3847/1538-4357/ad2704 Smith, M., Sullivan, M., Wiseman, P., et al. 2020, MNRAS, 494, 3,
-
[38]
doi:10.1093/mnras/staa946 Sullivan, M., Le Borgne, D., Pritchet, C. J., et al. 2006, ApJ, 648, 2,
-
[39]
doi:10.1086/506137 Sullivan, M., Conley, A., Howell, D. A., et al. 2010, MNRAS, 406, 2,
-
[40]
doi:10.1111/j.1365-2966.2010.16731.x Tripp, R. 1998, Astronomy and Astrophysics, 331,
-
[41]
E., et al
Vincenzi, M., Sullivan, M., Firth, R. E., et al. 2019, MNRAS, 489, 4,
2019
-
[42]
doi:10.1093/mnras/stz2448 Vincenzi, M., Brout, D., Armstrong, P., et al. 2024, ApJ, 975, 1,
-
[43]
doi:10.3847/1538-4357/ad5e6c Wiseman, P., Vincenzi, M., Sullivan, M., et al. 2022, MNRAS, 515, 3,
-
[44]
doi:10.1093/mnras/stac1984 1 UniversitÃl’ Lyon 1, CNRS, IP2I Lyon, UMR 5822, Villeurbanne, France 2 Institute of Astronomy and Kavli Institute for Cosmology, Mading- ley Road, Cambridge CB3 0HA, UK 3 School of Physics and Astronomy, University of Southampton, Southampton, SO17 1BJ, UK 4 The Oskar Klein Centre, Department of Physics, AlbaNova, SE-106 91 St...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.