A semi-empirical simulation of the extragalactic radio continuum sky for next generation radio telescopes

Monthly Notices of the Royal Astronomical Society 388:3 (2008) 1335-1348

Authors:

RJ Wilman, L Miller, MJ Jarvis, T Mauch, F Levrier, FB Abdalla, S Rawlings, HR Klöckner, D Obreschkow, D Olteanu, S Young

Abstract:

We have developed a semi-empirical simulation of the extragalactic radio continuum sky suitable for aiding the design of next generation radio interferometers such as the Square Kilometre Array (SKA). The emphasis is on modelling the large-scale cosmological distribution of radio sources rather than the internal structure of individual galaxies. Here we provide a description of the simulation to accompany the online release of a catalogue of ≃320 million simulated radio sources. The simulation covers a sky area of 20 × 20 deg2 - a plausible upper limit to the instantaneous field of view attainable with future (e.g. SKA) aperture array technologies - out to a cosmological redshift of z = 20, and down to flux density limits of 10 nJy at 151, 610 MHz, 1.4, 4.86 and 18 GHz. Five distinct source types are included: radio-quiet active galactic nuclei (AGN), radio-loud AGN of the Fanaroff-Riley type I (FR I) and FR II structural classes, and star-forming galaxies, the latter split into populations of quiescent and starbursting galaxies. In our semi-empirical approach, the simulated sources are drawn from observed (or extrapolated) luminosity functions and grafted on to an underlying dark matter density field with biases which reflect their measured large-scale clustering. A numerical Press-Schechter style filtering of the density field is used to identify and populate clusters of galaxies. For economy of output, radio source structures are constructed from point source and elliptical subcomponents, and for FR I and FR II sources an orientation-based unification and beaming model is used to partition flux between the core and extended lobes and hotspots. The extensive simulation output gives users the flexibility to post-process the catalogues to achieve more complete agreement with observational data in the years ahead. The ultimate aim is for the 'idealized skies' generated by this simulation and associated post-processing to be fed to telescope simulators to optimize the design of the SKA itself. © 2008 RAS.

Bayesian galaxy shape measurement for weak lensing surveys - II. Application to simulations

Monthly Notices of the Royal Astronomical Society 390:1 (2008) 149-167

Authors:

TD Kitching, L Miller, CE Heymans, L Van Waerbeke, AF Heavens

Abstract:

In this paper, we extend the Bayesian model fitting shape measurement method presented in Miller et al., and use the method to estimate the shear from the Shear TEsting Programme simulations (STEP). The method uses a fast model fitting algorithm that uses realistic galaxy profiles and analytically marginalizes over the position and amplitude of the model by doing the model fitting in Fourier space. This is used to find the full posterior probability in ellipticity. The shear is then estimated in a Bayesian way from this posterior probability surface. The Bayesian estimation allows measurement bias arising from the presence of random noise to be removed. In this paper, we introduce an iterative algorithm that can be used to estimate the intrinsic ellipticity prior and show that this is accurate and stable. We present results using the STEP parametrization that relates the input shear γT to the estimated shear γM by introducing a bias m and an offset c: γM - γT = mγT + c. The average number density of galaxies used in the STEP1 analysis was 9 per square arcminute, for STEP2 the number density was 30 per square arcminute. By using the method to estimate the shear from the STEP1 simulations we find the method to have a shear bias of m = 0.006 ± 0.005 and a variation in shear offset with point spread function type of σc = 0.0002. Using the method to estimate the shear from the STEP2 simulations we find that the shear bias and offset are m = 0.002 ± 0.016 and c = -0.0007 ± 0.0006, respectively. In addition, we find that the bias and offset are stable to changes in the magnitude and size of the galaxies. Such biases should yield any cosmological constraints from future weak lensing surveys robust to systematic effects in shape measurement. Finally, we present an alternative to the STEP parametrization by using a quality factor that relates the intrinsic shear variance in a simulation to the variance in shear that is measured and show that the method presented has an average of Q ≳ 100 which is at least a factor of 10 times better than other shape measurement methods. © 2008 RAS.

Constraining the quasar population with the broad-line width distribution

Monthly Notices of the Royal Astronomical Society 390:4 (2008) 1413-1429

Authors:

S Fine, SM Croom, PF Hopkins, L Hernquist, J Bland-Hawthorn, M Colless, PB Hall, L Miller, AD Myers, R Nichol, KA Pimbblet, NP Ross, DP Schneider, T Shanks, RG Sharp

Abstract:

In this work, we test the assertion that the scatter in the mass of black holes which drive quasars should be luminosity dependent with less scatter in more luminous objects. To this end, we measure the width of the Mg ii λ2799 line in quasar spectra from the Sloan Digital Sky Survey (SDSS), 2df QSO Redshift survey (2QZ) and 2dF SDSS LRG And QSO (2SLAQ) surveys and, by invoking an unnormalized virial mass estimator, relate the scatter in linewidth to the scatter of mass in the underlying black hole population. We find conclusive evidence for a trend such that there is less scatter in linewidth, and hence black hole mass, in more luminous objects. However, the most luminous objects in our sample show such a low degree of scatter in linewidth that, when combined with measures for the intrinsic scatter in the radius-luminosity relation for the broad-line region (BLR) in active galaxies, an inconsistency arises in the virial technique for estimating black hole masses. This analysis implies that, at least for the most luminous quasars, either there is little-to-no intrinsic scatter in the radius-luminosity relation or the Mg ii broad emission-line region is not totally dominated by virial velocities. Finally, we exploit the measured scatter in linewidths to constrain models for the velocity field of the BLR. We show that the lack of scatter in broad-line widths for luminous quasars is inconsistent with a pure planar/disc-like geometry for the BLR. In the case of a BLR with purely polar flows, the opening angle to luminous quasars must be less than ∼55°. We then explore the effects of adding a random or spherically symmetric component to the velocities of gas clouds in the BLR. Assuming an opening angle to quasars of 45°, a planar field can be made consistent with our results if ∼ 40-50 per cent of the velocities are randomly distributed. © 2008 The Authors.

Fourier-resolved energy spectra of the Narrow-Line Seyfert 1 Mkn 766

Monthly Notices of the Royal Astronomical Society 387:1 (2008) 279-288

Authors:

P Arévalo, IM McHardy, A Markowitz, IE Papadakis, TJ Turner, L Miller, J Reeves

Abstract:

We compute Fourier-resolved X-ray spectra of the Seyfert 1 Markarian 766 to study the shape of the variable components contributing to the 0.3-10 keV energy spectrum and their time-scale dependence. The fractional variability spectra peak at 1-3 keV, as in other Seyfert 1 galaxies, consistent with either a constant contribution from a soft excess component below 1 keV and Compton reflection component above 2 keV or variable warm absorption enhancing the variability in the 1-3 keV range. The rms spectra, which show the shape of the variable components only, are well described by a single power law with an absorption feature around 0.7 keV, which gives it an apparent soft excess. This spectral shape can be produced by a power law varying in normalization, affected by an approximately constant (within each orbit) warm absorber, with parameters similar to those found by Turner et al. for the warm-absorber layer covering all spectral components in their scattering scenario [NH ∼ 3 × 1021 cm-2, log(ξ) ∼ 1]. The total soft excess in the average spectrum can therefore be produced by a combination of constant warm absorption on the power-law plus an additional less variable component. On shorter time-scales, the rms spectrum hardens and this evolution is well described by a change in power-law slope, while the absorption parameters remain the same. The frequency dependence of the rms spectra can be interpreted as variability arising from propagating fluctuations through an extended emitting region, whose emitted spectrum is a power law that hardens towards the centre. This scenario reduces the short time-scale variability of lower energy bands making the variable spectrum harder on shorter time-scales and at the same time explains the hard lags found in these data by Markowitz et al. © 2008 RAS.

Multidimensional modelling of X-ray spectra for AGN accretion disc outflows

Monthly Notices of the Royal Astronomical Society 388:2 (2008) 611-624

Authors:

SA Sim, KS Long, L Miller, TJ Turner

Abstract:

We use a multidimensional Monte Carlo code to compute X-ray spectra for a variety of active galactic nucleus (AGN) disc-wind outflow geometries. We focus on the formation of blueshifted absorption features in the Fe K band and show that line features similar to those which have been reported in observations are often produced for lines of sight through disc-wind geometries. We also discuss the formation of other spectral features in highly ionized outflows. In particular, we show that, for sufficiently high wind densities, moderately strong Fe K emission lines can form and that electron scattering in the flow may cause these lines to develop extended red wings. We illustrate the potential relevance of such models to the interpretation of real X-ray data by comparison with observations of a well-known AGN, Mrk 766. Journal compilation © 2008 RAS.