Source author record

J. R. Peterson

J. R. Peterson appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

11works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

11 published item(s)

preprint2020arXiv

Sensor Distortion Effects in Photon Monte Carlo Simulations

We present a detailed method to simulating sensor distortions using a photon and electron Monte Carlo method. We use three dimensional electrostatic simulations to parameterize the perturbed electric field profile for non-ideal sensor details. We follow the conversion of simulated photons, and the subsequent response of the converted electrons to the electric field pattern. These non-ideal sensor details can be implemented efficiently in a Monte Carlo approach. We demonstrate that the non-ideal sensor distortions have a variety of observable consequence including the modification of the astrometric pattern, the distortion of the electron diffusion size and shape, and the distortion of flats. We show analytic validation of the diffusion physics, reproduce two kinds of edge distortion, and show qualitative validation of field-free regions, lithography errors, and fringing. We also demonstrate that there are two related effects of doping variation having different observable consequences. We show that field distortions from accumulated electrons lead to intensity-dependent point-spread-functions and the sub-linear variance in flats. The method is implemented in the Photon Simulator (PhoSim) and the code is publically available.

preprint2015arXiv

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

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.

preprint2013arXiv

Characterization of ICM Temperature Distributions of 62 Galaxy Clusters with XMM-Newton

We measure the intracluster medium temperature distributions for 62 galaxy clusters in the HIFLUGCS, an X-ray flux-limited sample, with available X-ray data from XMM-Newton. We search for correlations between the width of the temperature distributions and other cluster properties, including median cluster temperature, luminosity, size, presence of a cool core, AGN activity, and dynamical state. We use a Markov Chain Monte Carlo analysis which models the ICM as a collection of X-ray emitting smoothed particles of plasma. Each smoothed particle is given its own set of parameters, including temperature, spatial position, redshift, size, and emission measure. This allows us to measure the width of the temperature distribution, median temperature, and total emission measure of each cluster. We find that none of the clusters have a temperature width, σ_kT, consistent with isothermality. Counterintuitively, we also find that the temperature distribution widths of disturbed, non-cool-core, and AGN-free clusters tend to be wider than in other clusters. A linear fit to σ_kT - kT_med finds σ_kT ~ 0.20kT_med + 1.08, with an estimated intrinsic scatter of ~ 0.55 keV, demonstrating a large range in ICM thermal histories.

preprint2012arXiv

Atmospheric PSF Interpolation for Weak Lensing in Short Exposure Imaging Data

A main science goal for the Large Synoptic Survey Telescope (LSST) is to measure the cosmic shear signal from weak lensing to extreme accuracy. One difficulty, however, is that with the short exposure time ($\simeq$15 seconds) proposed, the spatial variation of the Point Spread Function (PSF) shapes may be dominated by the atmosphere, in addition to optics errors. While optics errors mainly cause the PSF to vary on angular scales similar or larger than a single CCD sensor, the atmosphere generates stochastic structures on a wide range of angular scales. It thus becomes a challenge to infer the multi-scale, complex atmospheric PSF patterns by interpolating the sparsely sampled stars in the field. In this paper we present a new method, PSFent, for interpolating the PSF shape parameters, based on reconstructing underlying shape parameter maps with a multi-scale maximum entropy algorithm. We demonstrate, using images from the LSST Photon Simulator, the performance of our approach relative to a 5th-order polynomial fit (representing the current standard) and a simple boxcar smoothing technique. Quantitatively, PSFent predicts more accurate PSF models in all scenarios and the residual PSF errors are spatially less correlated. This improvement in PSF interpolation leads to a factor of 3.5 lower systematic errors in the shear power spectrum on scales smaller than $\sim13'$, compared to polynomial fitting. We estimate that with PSFent and for stellar densities greater than $\simeq1/{\rm arcmin}^{2}$, the spurious shear correlation from PSF interpolation, after combining a complete 10-year dataset from LSST, is lower than the corresponding statistical uncertainties on the cosmic shear power spectrum, even under a conservative scenario.

preprint2012arXiv

Spurious Shear in Weak Lensing with LSST

The complete 10-year survey from the Large Synoptic Survey Telescope (LSST) will image $\sim$ 20,000 square degrees of sky in six filter bands every few nights, bringing the final survey depth to $r\sim27.5$, with over 4 billion well measured galaxies. To take full advantage of this unprecedented statistical power, the systematic errors associated with weak lensing measurements need to be controlled to a level similar to the statistical errors. This work is the first attempt to quantitatively estimate the absolute level and statistical properties of the systematic errors on weak lensing shear measurements due to the most important physical effects in the LSST system via high fidelity ray-tracing simulations. We identify and isolate the different sources of algorithm-independent, \textit{additive} systematic errors on shear measurements for LSST and predict their impact on the final cosmic shear measurements using conventional weak lensing analysis techniques. We find that the main source of the errors comes from an inability to adequately characterise the atmospheric point spread function (PSF) due to its high frequency spatial variation on angular scales smaller than $\sim10'$ in the single short exposures, which propagates into a spurious shear correlation function at the $10^{-4}$--$10^{-3}$ level on these scales. With the large multi-epoch dataset that will be acquired by LSST, the stochastic errors average out, bringing the final spurious shear correlation function to a level very close to the statistical errors. Our results imply that the cosmological constraints from LSST will not be severely limited by these algorithm-independent, additive systematic effects.

preprint2009arXiv

A direct limit on the turbulent velocity of the intracluster medium in the core of Abell 1835 from XMM-Newton

We examine deep XMM-Newton Reflection Grating Spectrometer (RGS) observations of the X-ray luminous galaxy cluster A1835. For the first time in a galaxy cluster we place direct limits on turbulent broadening of the emission lines. This is possible because the coolest X-ray emitting gas in the cluster, which is responsible for the lines, occupies a small region within the core. The most conservative determination of the 90 per cent upper limit on line-of-sight, non-thermal, velocity broadening is 274 km/s, measured from the emission lines originating within 30 kpc radius. The ratio of turbulent to thermal energy density in the core is therefore less than 13 per cent. There are no emission lines in the spectrum showing evidence for gas below ~3.5 keV. We examine the quantity of gas as a function of temperature and place a limit of 140 Msun/yr (90 per cent) for gas cooling radiatively below 3.85 keV.

preprint2009arXiv

Characterizing the Properties of Clusters of Galaxies as a Function of Luminosity and Redshift

We report the application of a new Monte Carlo method, Smoothed Particle Inference (SPI, described in a pair of companion papers), towards analysis and interpretation of X-ray observations of clusters of galaxies with the XMM-Newton satellite. Our sample consists of publicly available, well-exposed observations of clusters at redshifts z > 0.069, totaling 101 objects. We determine the luminosity and temperature structure of the X-ray emitting gas, with the goal to quantify the scatter and the evolution of the L_X - T relation, as well as to investigate the dependence on cluster substructure with redshift. We confirm that L_X \propto T^3 and we find a weak redshift dependence (\propto (1+z)^(β_LT), β_LT=0.50 +- 0.34), in contrast to some Chandra results. The level of dynamical activity is established using the "power ratios" method, and we find signs of evolution in the P_3/P_0 power ratio. A new method, the "temperature two-point correlation function," is proposed. This method is used to determine the "power spectrum" of temperature fluctuations in the X-ray emitting gas as a function of spatial scale. We show how this method can be fruitfully used to identify cooling core clusters as well as those with disturbed structures, presumably due to on-going or recent merger activity.

preprint2009arXiv

Deep high-resolution X-ray spectra from cool-core clusters

We examine deep XMM-Newton Reflection Grating Spectrometer (RGS) spectra from the cores of three X-ray bright cool core galaxy clusters, Abell 262, Abell 3581 and HCG 62. Each of the RGS spectra show Fe XVII emission lines indicating the presence of gas around 0.5 keV. There is no evidence for O VII emission which would imply gas at still cooler temperatures. The range in detected gas temperature in these objects is a factor of 3.7, 5.6 and 2 for Abell 262, Abell 3581 and HCG 62, respectively. The coolest detected gas only has a volume filling fraction of 6 and 3 per cent for Abell 262 and Abell 3581, but is likely to be volume filling in HCG 62. Chandra spatially resolved spectroscopy confirms the low volume filling fractions of the cool gas in Abell 262 and Abell 3581, indicating this cool gas exists as cold blobs. Any volume heating mechanism aiming to prevent cooling would overheat the surroundings of the cool gas by a factor of 4. If the gas is radiatively cooling below 0.5 keV, it is cooling at a rate at least an order of magnitude below that at higher temperatures in Abell 262 and Abell 3581 and two-orders of magnitude lower in HCG 62. The gas may be cooling non-radiatively through mixing in these cool blobs, where the energy released by cooling is emitted in the infrared. We find very good agreement between smooth particle inference modelling of the cluster and conventional spectral fitting. Comparing the temperature distribution from this analysis with that expected in a cooling flow, there appears to be a even larger break below 0.5 keV as compared with previous empirical descriptions of the deviations of cooling flow models.

preprint2009arXiv

Evidence for Non-Linear Growth of Structure from an X-ray Selected Cluster Survey using a Novel Joint Analysis of the Chandra and XMM-Newton Archives

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$σ$ 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.

preprint2007arXiv

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

We propose a new Monte Carlo method to study extended X-ray sources with the European Photon Imaging Camera (EPIC) aboard {\sl XMM-Newton}. The smoothed particle inference (SPI) technique, described in a companion paper, is applied here to the EPIC data for the clusters of galaxies Abell 1689, Centaurus, and RX J0658-55 (the ``bullet cluster''). We aim to show the advantages of this method of simultaneous spectral and spatial modeling over traditional X-ray spectral analysis. In Abell 1689 we confirm our earlier findings about structure in the temperature distribution and produce a high-resolution temperature map. We also find a hint of velocity structure within the gas, consistent with previous findings. In the bullet cluster, RX J0658-55, we produce the highest resolution temperature map ever to be published of this cluster, allowing us to trace what looks like the trail of the motion of the bullet in the cluster. We even detect a south-to-north temperature gradient within the bullet itself. In the Centaurus cluster we detect, by dividing up the luminosity of the cluster in bands of gas temperatures, a striking feature to the north-east of the cluster core. We hypothesize that this feature is caused by a subcluster left over from a substantial merger that slightly displaced the core. We conclude that our method is very powerful in determining the spatial distributions of plasma temperatures and very useful for systematic studies in cluster structure.

preprint2006arXiv

Smoothed Particle Inference: A Kilo-Parametric Method for X-ray Galaxy Cluster Modeling

We propose an ambitious new method that models the intracluster medium in clusters of galaxies as a set of X-ray emitting smoothed particles of plasma. Each smoothed particle is described by a handful of parameters including temperature, location, size, and elemental abundances. Hundreds to thousands of these particles are used to construct a model cluster of galaxies, with the appropriate complexity estimated from the data quality. This model is then compared iteratively with X-ray data in the form of adaptively binned photon lists via a two-sample likelihood statistic and iterated via Markov Chain Monte Carlo. The complex cluster model is propagated through the X-ray instrument response using direct sampling Monte Carlo methods. Using this approach the method can reproduce many of the features observed in the X-ray emission in a less assumption-dependent way that traditional analyses, and it allows for a more detailed characterization of the density, temperature, and metal abundance structure of clusters. Multi-instrument X-ray analyses and simultaneous X-ray, Sunyaev-Zeldovich (SZ), and lensing analyses are a straight-forward extension of this methodology. Significant challenges still exist in understanding the degeneracy in these models and the statistical noise induced by the complexity of the models.