Image Ellipticity from Atmospheric Aberrations

We investigate the ellipticity of the point-spread function (PSF) produced by imaging an unresolved source with a telescope, subject to the effects of atmospheric turbulence. It is important to quantify these effects in order to understand the errors in shape measurements of astronomical objects, such as those used to study weak gravitational lensing of field galaxies. The PSF modeling involves either a Fourier transform of the phase information in the pupil plane or a ray-tracing approach, which has the advantage of requiring fewer computations than the Fourier transform. Using a standard method, involving the Gaussian weighted second moments of intensity, we then calculate the ellipticity of the PSF patterns. We find significant ellipticity for the instantaneous patterns (up to more than 10%). Longer exposures, which we approximate by combining multiple (N) images from uncorrelated atmospheric realizations, yield progressively lower ellipticity (as 1 / sqrt(N)). We also verify that the measured ellipticity does not depend on the sampling interval in the pupil plane using the Fourier method. However, we find that the results using the ray-tracing technique do depend on the pupil sampling interval, representing a gradual breakdown of the geometric approximation at high spatial frequencies. Therefore, ray tracing is generally not an accurate method of modeling PSF ellipticity induced by atmospheric turbulence unless some additional procedure is implemented to correctly account for the effects of high spatial frequency aberrations. The Fourier method, however, can be used directly to accurately model PSF ellipticity, which can give insights into errors in the statistics of field galaxy shapes used in studies of weak gravitational lensing.Comment: 9 pages, 5 color figures (some reduced in size). Accepted for publication in the Astrophysical Journa

The Astrophysical Multipurpose Software Environment

We present the open source Astrophysical Multi-purpose Software Environment (AMUSE, www.amusecode.org), a component library for performing astrophysical simulations involving different physical domains and scales. It couples existing codes within a Python framework based on a communication layer using MPI. The interfaces are standardized for each domain and their implementation based on MPI guarantees that the whole framework is well-suited for distributed computation. It includes facilities for unit handling and data storage. Currently it includes codes for gravitational dynamics, stellar evolution, hydrodynamics and radiative transfer. Within each domain the interfaces to the codes are as similar as possible. We describe the design and implementation of AMUSE, as well as the main components and community codes currently supported and we discuss the code interactions facilitated by the framework. Additionally, we demonstrate how AMUSE can be used to resolve complex astrophysical problems by presenting example applications.Comment: 23 pages, 25 figures, accepted for A&

Investigating the interstellar dust through the Fe K-edge

The chemical and physical properties of interstellar dust in the densest regions of the Galaxy are still not well understood. X-rays provide a powerful probe since they can penetrate gas and dust over a wide range of column densities (up to $10^{24}\ \rm{cm}^{-2}$). The interaction (scattering and absorption) with the medium imprints spectral signatures that reflect the individual atoms which constitute the gas, molecule, or solid. In this work we investigate the ability of high resolution X-ray spectroscopy to probe the properties of cosmic grains containing iron. Although iron is heavily depleted into interstellar dust, the nature of the Fe-bearing grains is still largely uncertain. In our analysis we use iron K-edge synchrotron data of minerals likely present in the ISM dust taken at the European Synchrotron Radiation Facility. We explore the prospects of determining the chemical composition and the size of astrophysical dust in the Galactic centre and in molecular clouds with future X-ray missions. The energy resolution and the effective area of the present X-ray telescopes are not sufficient to detect and study the Fe K-edge, even for bright X-ray sources. From the analysis of the extinction cross sections of our dust models implemented in the spectral fitting program SPEX, the Fe K-edge is promising for investigating both the chemistry and the size distribution of the interstellar dust. We find that the chemical composition regulates the X-ray absorption fine structures in the post edge region, whereas the scattering feature in the pre-edge is sensitive to the mean grain size. Finally, we note that the Fe K-edge is insensitive to other dust properties, such as the porosity and the geometry of the dust.Comment: 11 pages, 10 figures. Accepted for publication in Astronomy and Astrophysic

Long Term Variability of SDSS Quasars

We use a sample of 3791 quasars from the Sloan Digital Sky Survey (SDSS) Early Data Release (EDR), and compare their photometry to historic plate material for the same set of quasars in order to study their variability properties. The time base-line we attain this way ranges from a few months to up to 50 years. In contrast to monitoring programs, where relatively few quasars are photometrically measured over shorter time periods, we utilize existing databases to extend this base-line as much as possible, at the cost of sampling per quasar. Our method, however, can easily be extended to much larger samples. We construct variability Structure Functions and compare these to the literature and model functions. From our modeling we conclude that 1) quasars are more variable toward shorter wavelengths, 2) their variability is consistent with an exponentially decaying light-curve with a typical time-scale of ~2 years, 3) these outbursts occur on typical time-scales of ~200 years. With the upcoming first data release of the SDSS, a much larger quasar sample can be used to put these conclusions on a more secure footing.Comment: 16 pages, accepted for publication in AJ, Sept issu

Molecular Line Profile Fitting with Analytic Radiative Transfer Models

We present a study of analytic models of starless cores whose line profiles have ``infall asymmetry,'' or blue-skewed shapes indicative of contracting motions. We compare the ability of two types of analytical radiative transfer models to reproduce the line profiles and infall speeds of centrally condensed starless cores whose infall speeds are spatially constant and range between 0 and 0.2 km s-1. The model line profiles of HCO+ (J=1-0) and HCO+ (J=3-2) are produced by a self-consistent Monte Carlo radiative transfer code. The analytic models assume that the excitation temperature in the front of the cloud is either constant (``two-layer'' model) or increases inward as a linear function of optical depth (``hill'' model). Each analytic model is matched to the line profile by rapid least-squares fitting. The blue-asymmetric line profiles with two peaks, or with a blue shifted peak and a red shifted shoulder, can be well fit by the ``HILL5'' model (a five parameter version of the hill model), with an RMS error of 0.02 km s-1. A peak signal to noise ratio of at least 30 in the molecular line observations is required for performing these analytic radiative transfer fits to the line profiles.Comment: 48 pages, 20 figures, accepted for publication in Ap

Lower bounds on photometric redshift errors from Type Ia supernovae templates

Cosmology with Type Ia supernovae heretofore has required extensive spectroscopic follow-up to establish a redshift. Though tolerable at the present discovery rate, the next generation of ground-based all-sky survey instruments will render this approach unsustainable. Photometry-based redshift determination is a viable alternative, but introduces non-negligible errors that ultimately degrade the ability to discriminate between competing cosmologies. We present a strictly template-based photometric redshift estimator and compute redshift reconstruction errors in the presence of photometry and statistical errors. With reasonable assumptions for a cadence and supernovae distribution, these redshift errors are combined with systematic errors and propagated using the Fisher matrix formalism to derive lower bounds on the joint errors in $\Omega_w$ and $\Omega_w'$ relevant to the next generation of ground-based all-sky survey.Comment: 23 pages, 6 figure

The Cape Triage Score - a triage system for South Africa

The Cape Triage Score (CTS) has been derived by the Cape Triage Group (CTG) for use in emergency units throughout South Africa. It can also be used in the pre-hospital setting, although it is not designed for mass casualty situations. The CTS comprises a physiologically based scoring system and a list of discriminators, designed to triage patients into one of five priority groups for medical attention. Three versions have been developed, for adults, children and infants. As part of the ongoing assessment process the CTG would value feedback from the readers of this Journal