paper.md

title	tags	authors	affiliations	date	bibliography
pyuvsim: A comprehensive simulation package for radio interferometers in python.	radio astronomy simulation pyuvdata	name orcid affiliation Adam E. Lanman 0000-0003-2116-3573 4, 7 name orcid affiliation Bryna J. Hazelton 0000-0001-7532-645X 1, 2, 7 name orcid affiliation Daniel C. Jacobs 0000-0002-0917-2269 3, 7 name orcid affiliation Matthew J. Kolopanis 0000-0002-2950-2974 3, 7 name orcid affiliation Jonathan C. Pober 0000-0002-3492-0433 4, 7 name orcid affiliation James E. Aguirre 0000-0002-4810-666X 6, 7 name orcid affiliation Nithyanandan Thyagarajan 0000-0003-1602-7868 5, 7	name index University of Washington, eScience Institute 1 name index University of Washington, Physics Department 2 name index Arizona State University, School of Earth and Space Exploration 3 name index Brown University, Physics Department 4 name index National Radio Astronomy Observatory 5 name index University of Pennsylvania, Physics Department 6 name index Radio Astronomy Software Group 7	29 December 2019	paper.bib

Summary

Data simulations are essential to the progress of low-frequency radio telescope arrays such as HERA, the MWA, and LOFAR, and for future instruments like the SKA. Simulated datasets are used to verify analysis pipelines, to provide models for sky-based calibration, and to test the effects of design choices and environmental factors. Existing simulators (CASA, OSKAR, FHD, PRISim, et al.) make simplifying assumptions to reduce the computational demand of evaluating the measurement equation, compromising accuracy for speed. This can lead to unexpected numerical effects, which can be hard to distinguish from real effects seen in data.

pyuvsim is a simulator designed to be accurate and verifiable up front, with strict control over any approximations being made and including all effects found to be important to 21cm cosmological experiments. The default behavior is to add contributions from each point source above the horizon to each baseline in a full Jones-matrix formalism [@hamaker_understanding_1996], with floating precision source positions and interpolated E-field primary beam values. Results are tested against analytic calculations whenever possible, and compared with other simulators and data for consistency, including PRISim [@nithyanandan_2019_2548117] and FHD [@fhd]. Data handling and primary beam modeling are done using the pyuvdata package [@j_hazelton_pyuvdata:_2017], which supports a variety of output data types and ensures accurate phasing methods. Source motions and coherency rotations are calculated using astropy transformations [@astropy:2013], which takes into account higher order corrections to sky motion.

Currently, pyuvsim supports simulating point sources from the GLEAM catalog [@hurley-walker_galactic_2017] and several mock source catalogs. It includes utilities for making parameter files based on the structure of any data file readable by pyuvdata, allowing users to quickly set up simulations to compare with an observation. It is parallelized using the Message Passing Interface (MPI) [@dalcin2005mpi]. Performance improvements are made by identifying and correcting bottlenecks using built-in line profiling tools. Features to be added include support for spectral models for sources, diffuse source models, and ionospheric effects. These are not expected to change the API. Users can rely on pyuvsim for precise results with steady, long-term improvements in resource usage and runtime.

Acknowledgements

This work was supported in part by NASA award 80NSSC18K0389 and NSF award AST-1613973.

References