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 sim
ulating 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.
We present a large X-ray selected serendipitous cluster survey based on a novel joint analysis of archival Chandra and XMM-Newton data. The survey provides enough depth to reach clusters of flux of $approx 10^{-14} {ergs} {cm}^{-2} {s}^{-1}$ near $z$
$approx$ 1 and simultaneously a large enough sample to find evidence for the strong evolution of clusters expected from structure formation theory. We detected a total of 723 clusters of which 462 are newly discovered clusters with greater than 6$sigma$ significance. In addition, we also detect and measure 261 previously-known clusters and groups that can be used to calibrate the survey. The survey exploits a technique which combines the exquisite Chandra imaging quality with the high throughput of the XMM-Newton telescopes using overlapping survey regions. A large fraction of the contamination from AGN point sources is mitigated by using this technique. This results in a higher sensitivity for finding clusters of galaxies with relatively few photons and a large part of our survey has a flux sensitivity between $10^{-14}$ and $10^{-15} {ergs} {cm}^{-2} {s}^{-1}$. The survey covers 41.2 square degrees of overlapping Chandra and XMM-Newton fields and 122.2 square degrees of non-overlapping Chandra data. We measure the log N-log S distribution and fit it with a redshift-dependent model characterized by a luminosity distribution proportional to $e^{-frac{z}{z_0}}$. We find that $z_0$ to be in the range 0.7 to 1.3, indicative of rapid cluster evolution, as expected for cosmic structure formation using parameters appropriate to the concordance cosmological model.
