Simulation of Astronomical Images from Optical Survey Telescopes using a Comprehensive Photon Monte Carlo Approach

Abstract

We present a comprehensive methodology for the simulation of astronomical images from optical survey telescopes. We use a photon Monte Carlo approach to construct images by sampling photons from models of astronomical source populations, and then simulating those photons through the system as they interact with the atmosphere, telescope, and camera. We demonstrate that all physical effects for optical light that determine the shapes, locations, and brightnesses of individual stars and galaxies can be accurately represented in this formalism. By using large scale grid computing, modern processors, and an efficient implementation that can produce 400,000 photons/second, we demonstrate that even very large optical surveys can be now be simulated. We demonstrate that we are able to: 1) construct kilometer scale phase screens necessary for wide-field telescopes, 2) reproduce atmospheric point-spread-function moments using a fast novel hybrid geometric/Fourier technique for non-diffraction limited telescopes, 3) accurately reproduce the expected spot diagrams for complex aspheric optical designs, and 4) recover system effective area predicted from analytic photometry integrals. This new code, the photon simulator (PhoSim), is publicly available. We have implemented the Large Synoptic Survey Telescope (LSST) design, and it can be extended to other telescopes. We expect that because of the comprehensive physics implemented in PhoSim, it will be used by the community to plan future observations, interpret detailed existing observations, and quantify systematics related to various astronomical measurements. Future development and validation by comparisons with real data will continue to improve the fidelity and usability of the code.

arXiv commentAccepted for publication in ApJS; 34 pages, 15 figures

References in corpus (20)

• Observational Evidence from Supernovae for an Accelerating Universe and a Cosmological Constant
• Measurements of Omega and Lambda from 42 High-Redshift Supernovae
• The cosmological simulation code GADGET-2
• The Dust Content and Opacity of Actively Star-Forming Galaxies
• Detection of the Baryon Acoustic Peak in the Large-Scale Correlation Function of SDSS Luminous Red Galaxies
• The USNO-B Catalog
• Athena: A New Code for Astrophysical MHD
• The Fermi Large Area Telescope On Orbit: Event Classification, Instrument Response Functions, and Calibration
• SDSS Standard Star Catalog for Stripe 82: the Dawn of Industrial 1% Optical Photometry
• The brighter-fatter effect and pixel correlations in CCD sensors
• Precision simulation of ground-based lensing data using observations from space
• Synthetic Modeling of Astronomical Closed Loop Adaptive Optics
• UBVRI twilight sky brightness at ESO-Paranal
• Testing turbulence model at metric scales with mid-infrared VISIR images at the VLT
• Atmospheric image blur with finite outer scale or partial adaptive correction
• Smoothed Particle Inference: A Kilo-Parametric Method for X-ray Galaxy Cluster Modeling
• Astronomical Image Simulation for Telescope and Survey Development
• Image Ellipticity from Atmospheric Aberrations
• Application of an XMM-Newton EPIC Monte Carlo Technique to Analysis and Interpretation of Data for the Abell 1689, RX J0658-55 and Centaurus clusters of galaxies
• Multivariate Monte Carlo Methods for the Reflection Grating Spectrometers on XMM-Newton

Cited by in corpus (41)

• LSST: from Science Drivers to Reference Design and Anticipated Data Products
• Models and Simulations for the Photometric LSST Astronomical Time Series Classification Challenge (PLAsTiCC)
• Lensing substructure quantification in RXJ1131-1231: A 2 keV lower bound on dark matter thermal relic mass
• Deblending and Classifying Astronomical Sources with Mask R-CNN Deep Learning
• The Prospects of Observing Tidal Disruption Events with the LSST
• The LSST Data Management System
• Approximate Bayesian Computation for Forward Modeling in Cosmology
• The LSST DESC DC2 Simulated Sky Survey
• Curvature Wavefront Sensing for the Large Synoptic Survey Telescope
• Detection, Instance Segmentation, and Classification for Astronomical Surveys with Deep Learning (DeepDISC): Detectron2 Implementation and Demonstration with Hyper Suprime-Cam Data
• A computationally efficient autoregressive method for generating phase screens with frozen flow and turbulence in optical simulations
• Euclid preparation XXVI. The Euclid Morphology Challenge. Towards structural parameters for billions of galaxies
• Generative deep fields: arbitrarily sized, random synthetic astronomical images through deep learning
• Euclid preparation. XXV. The Euclid Morphology Challenge -- Towards model-fitting photometry for billions of galaxies
• Wavelength Dependent PSFs and their impact on Weak Lensing Measurements
• HIDE & SEEK: End-to-End Packages to Simulate and Process Radio Survey Data
• Point spread function modelling for astronomical telescopes: a review focused on weak gravitational lensing studies
• The LSST DESC Data Challenge 1: Generation and Analysis of Synthetic Images for Next Generation Surveys
• Effects of overlapping sources on cosmic shear estimation: Statistical sensitivity and pixel-noise bias
• A Study of the Point Spread Function in SDSS Images
• Structured Variational Inference for Simulating Populations of Radio Galaxies
• The Brighter-Fatter and other Sensor Effects in CCD Simulations for Precision Astronomy
• The Point Spread Function Reconstruction by Using Moffatlets - I
• Enabling catalog simulations of transient and variable sources based on LSST cadence strategies
• A framework for modeling the detailed optical response of thick, multiple segment, large format sensors for precision astronomy applications
• Testing the Impact of Satellite Anisotropy on Large and Small Scale Intrinsic Alignments using Hydrodynamical Simulations
• Sensor Distortion Effects in Photon Monte Carlo Simulations
• Deformation of Optics for Photon Monte Carlo Simulations
• PhoSim-NIRCam: Photon-by-photon image simulations of the James Webb Space Telescope's Near-Infrared Camera
• Next-generation telescopes with curved focal surface for ultra-low surface brightness surveys
• Weak Lensing analysis of Abell 2390 using short exposures
• The analysis of effective galaxies number count for Chinese Space Station Optical Survey(CSS-OS) by image simulation
• Masked areas in shear peak statistics: a forward modeling approach
• Prompt Detection of Fast Optical Bursts with the Vera C. Rubin Observatory
• Euclid preparation. Spatially resolved stellar populations of local galaxies with Euclid: a proof of concept using synthetic images with the TNG50 simulation
• A hybrid Fourier--Real Gaussian Mixture method for fast galaxy--PSF convolution
• Searching for Sub-Second Stellar Variability with Wide-Field Star Trails and Deep Learning
• Using AI for Wavefront Estimation with the Rubin Observatory Active Optics System
• Neural Posterior Estimation for Cataloging Astronomical Images with Spatially Varying Backgrounds and Point Spread Functions
• Simulating Photometric Images of Moving Targets with Photon-mapping
• The WEAVE acquisition and guiding software: pattern recognition-based acquisition and multi-fibre guiding