Source author record

Paul Alexander

Paul Alexander 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

29works
8topics
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

29 published item(s)

preprint2022arXiv

Automated Detection of Antenna Malfunctions in Large-N Interferometers: A Case Study with the Hydrogen Epoch of Reionization Array

We present a framework for identifying and flagging malfunctioning antennas in large radio interferometers. We outline two distinct categories of metrics designed to detect outliers along known failure modes of large arrays: cross-correlation metrics, based on all antenna pairs, and auto-correlation metrics, based solely on individual antennas. We define and motivate the statistical framework for all metrics used, and present tailored visualizations that aid us in clearly identifying new and existing systematics. We implement these techniques using data from 105 antennas in the Hydrogen Epoch of Reionization Array (HERA) as a case study. Finally, we provide a detailed algorithm for implementing these metrics as flagging tools on real data sets.

preprint2021arXiv

Demonstrations of Cooperative Perception: Safety and Robustness in Connected and Automated Vehicle Operations

Cooperative perception, or collective perception (CP) is an emerging and promising technology for intelligent transportation systems (ITS). It enables an ITS station (ITS-S) to share its local perception information with others by means of vehicle-to-X (V2X) communication, thereby achieving improved efficiency and safety in road transportation. In this paper, we present our recent progress on the development of a connected and automated vehicle (CAV) and intelligent roadside unit (IRSU). We present three different experiments to demonstrate the use of CP service within intelligent infrastructure to improve awareness of vulnerable road users (VRU) and thus safety for CAVs in various traffic scenarios. We demonstrate in the experiments that a connected vehicle (CV) can "see" a pedestrian around the corners. More importantly, we demonstrate how CAVs can autonomously and safely interact with walking and running pedestrians, relying only on the CP information from the IRSU through vehicle-to-infrastructure (V2I) communication. This is one of the first demonstrations of urban vehicle automation using only CP information. We also address in the paper the handling of collective perception messages (CPMs) received from the IRSU, and passing them through a pipeline of CP information coordinate transformation with uncertainty, multiple road user tracking, and eventually path planning/decision making within the CAV. The experimental results were obtained with manually driven CV, fully autonomous CAV, and an IRSU retrofitted with vision and laser sensors and a road user tracking system.

preprint2021arXiv

First Results from HERA Phase I: Upper Limits on the Epoch of Reionization 21 cm Power Spectrum

We report upper-limits on the Epoch of Reionization (EoR) 21 cm power spectrum at redshifts 7.9 and 10.4 with 18 nights of data ($\sim36$ hours of integration) from Phase I of the Hydrogen Epoch of Reionization Array (HERA). The Phase I data show evidence for systematics that can be largely suppressed with systematic models down to a dynamic range of $\sim10^9$ with respect to the peak foreground power. This yields a 95% confidence upper limit on the 21 cm power spectrum of $Δ^2_{21} \le (30.76)^2\ {\rm mK}^2$ at $k=0.192\ h\ {\rm Mpc}^{-1}$ at $z=7.9$, and also $Δ^2_{21} \le (95.74)^2\ {\rm mK}^2$ at $k=0.256\ h\ {\rm Mpc}^{-1}$ at $z=10.4$. At $z=7.9$, these limits are the most sensitive to-date by over an order of magnitude. While we find evidence for residual systematics at low line-of-sight Fourier $k_\parallel$ modes, at high $k_\parallel$ modes we find our data to be largely consistent with thermal noise, an indicator that the system could benefit from deeper integrations. The observed systematics could be due to radio frequency interference, cable sub-reflections, or residual instrumental cross-coupling, and warrant further study. This analysis emphasizes algorithms that have minimal inherent signal loss, although we do perform a careful accounting in a companion paper of the small forms of loss or bias associated with the pipeline. Overall, these results are a promising first step in the development of a tuned, instrument-specific analysis pipeline for HERA, particularly as Phase II construction is completed en route to reaching the full sensitivity of the experiment.

preprint2021arXiv

Validation of the HERA Phase I Epoch of Reionization 21 cm Power Spectrum Software Pipeline

We describe the validation of the HERA Phase I software pipeline by a series of modular tests, building up to an end-to-end simulation. The philosophy of this approach is to validate the software and algorithms used in the Phase I upper limit analysis on wholly synthetic data satisfying the assumptions of that analysis, not addressing whether the actual data meet these assumptions. We discuss the organization of this validation approach, the specific modular tests performed, and the construction of the end-to-end simulations. We explicitly discuss the limitations in scope of the current simulation effort. With mock visibility data generated from a known analytic power spectrum and a wide range of realistic instrumental effects and foregrounds, we demonstrate that the current pipeline produces power spectrum estimates that are consistent with known analytic inputs to within thermal noise levels (at the 2 sigma level) for k > 0.2 h/Mpc for both bands and fields considered. Our input spectrum is intentionally amplified to enable a strong `detection' at k ~0.2 h/Mpc -- at the level of ~25 sigma -- with foregrounds dominating on larger scales, and thermal noise dominating at smaller scales. Our pipeline is able to detect this amplified input signal after suppressing foregrounds with a dynamic range (foreground to noise ratio) of > 10^7. Our validation test suite uncovered several sources of scale-independent signal loss throughout the pipeline, whose amplitude is well-characterized and accounted for in the final estimates. We conclude with a discussion of the steps required for the next round of data analysis.

preprint2020arXiv

Absolute Calibration Strategies for the Hydrogen Epoch of Reionization Array and Their Impact on the 21 cm Power Spectrum

We discuss absolute calibration strategies for Phase I of the Hydrogen Epoch of Reionization Array (HERA), which aims to measure the cosmological 21 cm signal from the Epoch of Reionization (EoR). HERA is a drift-scan array with a 10 degree wide field of view, meaning bright, well-characterized point source transits are scarce. This, combined with HERA's redundant sampling of the uv plane and the modest angular resolution of the Phase I instrument, make traditional sky-based and self-calibration techniques difficult to implement with high dynamic range. Nonetheless, in this work we demonstrate calibration for HERA using point source catalogues and electromagnetic simulations of its primary beam. We show that unmodeled diffuse flux and instrumental contaminants can corrupt the gain solutions, and present a gain smoothing approach for mitigating their impact on the 21 cm power spectrum. We also demonstrate a hybrid sky and redundant calibration scheme and compare it to pure sky-based calibration, showing only a marginal improvement to the gain solutions at intermediate delay scales. Our work suggests that the HERA Phase I system can be well-calibrated for a foreground-avoidance power spectrum estimator by applying direction-independent gains with a small set of degrees of freedom across the frequency and time axes.

preprint2020arXiv

Detection of Cosmic Structures using the Bispectrum Phase. II. First Results from Application to Cosmic Reionization Using the Hydrogen Epoch of Reionization Array

Characterizing the epoch of reionization (EoR) at $z\gtrsim 6$ via the redshifted 21 cm line of neutral Hydrogen (HI) is critical to modern astrophysics and cosmology, and thus a key science goal of many current and planned low-frequency radio telescopes. The primary challenge to detecting this signal is the overwhelmingly bright foreground emission at these frequencies, placing stringent requirements on the knowledge of the instruments and inaccuracies in analyses. Results from these experiments have largely been limited not by thermal sensitivity but by systematics, particularly caused by the inability to calibrate the instrument to high accuracy. The interferometric bispectrum phase is immune to antenna-based calibration and errors therein, and presents an independent alternative to detect the EoR HI fluctuations while largely avoiding calibration systematics. Here, we provide a demonstration of this technique on a subset of data from the Hydrogen Epoch of Reionization Array (HERA) to place approximate constraints on the brightness temperature of the intergalactic medium (IGM). From this limited data, at $z=7.7$ we infer "$1σ$" upper limits on the IGM brightness temperature to be $\le 316$ "pseudo" mK at $κ_\parallel=0.33$ "pseudo" $h$ Mpc$^{-1}$ (data-limited) and $\le 1000$ "pseudo" mK at $κ_\parallel=0.875$ "pseudo" $h$ Mpc$^{-1}$ (noise-limited). The "pseudo" units denote only an approximate and not an exact correspondence to the actual distance scales and brightness temperatures. By propagating models in parallel to the data analysis, we confirm that the dynamic range required to separate the cosmic HI signal from the foregrounds is similar to that in standard approaches, and the power spectrum of the bispectrum phase is still data-limited (at $\gtrsim 10^6$ dynamic range) indicating scope for further improvement in sensitivity as the array build-out continues.

preprint2020arXiv

Foreground modelling via Gaussian process regression: an application to HERA data

The key challenge in the observation of the redshifted 21-cm signal from cosmic reionization is its separation from the much brighter foreground emission. Such separation relies on the different spectral properties of the two components, although, in real life, the foreground intrinsic spectrum is often corrupted by the instrumental response, inducing systematic effects that can further jeopardize the measurement of the 21-cm signal. In this paper, we use Gaussian Process Regression to model both foreground emission and instrumental systematics in $\sim 2$ hours of data from the Hydrogen Epoch of Reionization Array. We find that a simple co-variance model with three components matches the data well, giving a residual power spectrum with white noise properties. These consist of an "intrinsic" and instrumentally corrupted component with a coherence-scale of 20 MHz and 2.4 MHz respectively (dominating the line of sight power spectrum over scales $k_{\parallel} \le 0.2$ h cMpc$^{-1}$) and a baseline dependent periodic signal with a period of $\sim 1$ MHz (dominating over $k_{\parallel} \sim 0.4 - 0.8$h cMpc$^{-1}$) which should be distinguishable from the 21-cm EoR signal whose typical coherence-scales is $\sim 0.8$ MHz.

preprint2020arXiv

Measuring HERA's primary beam in-situ: methodology and first results

The central challenge in 21~cm cosmology is isolating the cosmological signal from bright foregrounds. Many separation techniques rely on the accurate knowledge of the sky and the instrumental response, including the antenna primary beam. For drift-scan telescopes such as the Hydrogen Epoch of Reionization Array \citep[HERA, ][]{DeBoer2017} that do not move, primary beam characterization is particularly challenging because standard beam-calibration routines do not apply \citep{Cornwell2005} and current techniques require accurate source catalogs at the telescope resolution. We present an extension of the method from \citet{Pober2012} where they use beam symmetries to create a network of overlapping source tracks that break the degeneracy between source flux density and beam response and allow their simultaneous estimation. We fit the beam response of our instrument using early HERA observations and find that our results agree well with electromagnetic simulations down to a -20~dB level in power relative to peak gain for sources with high signal-to-noise ratio. In addition, we construct a source catalog with 90 sources down to a flux density of 1.4~Jy at 151~MHz.

preprint2019arXiv

Mitigating Internal Instrument Coupling II: A Method Demonstration with the Hydrogen Epoch of Reionization Array

We present a study of internal reflection and cross coupling systematics in Phase I of the Hydrogen Epoch of Reionization Array (HERA). In a companion paper, we outlined the mathematical formalism for such systematics and presented algorithms for modeling and removing them from the data. In this work, we apply these techniques to data from HERA's first observing season as a method demonstration. The data show evidence for systematics that, without removal, would hinder a detection of the 21 cm power spectrum for the targeted EoR line-of-sight modes in the range 0.2 < k_parallel < 0.5\ h^-1 Mpc. After systematic removal, we find we can recover these modes in the power spectrum down to the integrated noise-floor of a nightly observation, achieving a dynamic range in the EoR window of 10^-6 in power (mK^2 units) with respect to the bright galactic foreground signal. In the absence of other systematics and assuming the systematic suppression demonstrated here continues to lower noise levels, our results suggest that fully-integrated HERA Phase I may have the capacity to set competitive upper limits on the 21 cm power spectrum. For future observing seasons, HERA will have upgraded analog and digital hardware to better control these systematics in the field.

preprint2016arXiv

Contamination of the Epoch of Reionization power spectrum in the presence of foregrounds

We construct foreground simulations comprising spatially correlated extragalactic and diffuse Galactic emission components and calculate the `intrinsic' (instrument-free) two-dimensional spatial power spectrum and the cylindrically and spherically averaged three-dimensional k-space power spectra of the Epoch of Reionization (EoR) and our foreground simulations using a Bayesian power spectral estimation framework. This leads us to identify a model dependent region of optimal signal estimation for our foreground and EoR models, within which the spatial power in the EoR signal relative to foregrounds is maximised. We identify a target field dependent region, in k-space, of intrinsic foreground power spectral contamination at low k_perp and k_parallel and a transition to a relatively foreground-free intrinsic EoR window in the complement to this region. The contaminated region of k-space demonstrates that simultaneous estimation of the EoR and foregrounds is important for obtaining statistically robust estimates of the EoR power spectrum; biased results will be obtained from methodologies that ignore their covariance. Using simulated observations with frequency dependent uv-coverage and primary beam, with the former derived for HERA in 37-antenna and 331-antenna configuration, we recover instrumental power spectra consistent with their intrinsic counterparts. We discuss the implications of these results for optimal strategies for unbiased estimation of the EoR power spectrum.

preprint2015arXiv

SKA synergy with Microwave Background studies

The extremely high sensitivity and resolution of the Square Kilometre Array (SKA) will be useful for addressing a wide set of themes relevant for cosmology, in synergy with current and future cosmic microwave background (CMB) projects. Many of these themes also have a link with future optical-IR and X-ray observations. We discuss the scientific perspectives for these goals, the instrumental requirements and the observational and data analysis approaches, and identify several topics that are important for cosmology and astrophysics at different cosmic epochs.

preprint2014arXiv

A Bayesian blind survey for cold molecular gas in the Universe

A new Bayesian method for performing an image domain search for line-emitting galaxies is presented. The method uses both spatial and spectral information to robustly determine the source properties, employing either simple Gaussian, or other physically motivated models whilst using the evidence to determine the probability that the source is real. In this paper, we describe the method, and its application to both a simulated data set, and a blind survey for cold molecular gas using observations of the Hubble Deep Field North taken with the Plateau de Bure Interferometer. We make a total of 6 robust detections in the survey, 5 of which have counterparts in other observing bands. We identify the most secure detections found in a previous investigation, while finding one new probable line source with an optical ID not seen in the previous analysis. This study acts as a pilot application of Bayesian statistics to future searches to be carried out both for low-$J$ CO transitions of high redshift galaxies using the JVLA, and at millimeter wavelengths with ALMA, enabling the inference of robust scientific conclusions about the history of the molecular gas properties of star-forming galaxies in the Universe through cosmic time.

preprint2014arXiv

Bayesian estimation of non-Gaussianity in pulsar timing analysis

We introduce a method for performing a robust Bayesian analysis of non-Gaussianity present in pulsar timing data, simultaneously with the pulsar timing model, and additional stochastic parameters such as those describing red spin noise and dispersion measure variations. The parameters used to define the presence of non-Gaussianity are zero for Gaussian processes, giving a simple method of defining the strength of non-Gaussian behaviour. We use simulations to show that assuming Gaussian statistics when the noise in the data is drawn from a non-Gaussian distribution can significantly increase the uncertainties associated with the pulsar timing model parameters. We then apply the method to the publicly available 15 year Parkes Pulsar Timing Array data release 1 dataset for the binary pulsar J0437$-$4715. In this analysis we present a significant detection of non-Gaussianity in the uncorrelated non-thermal noise, but we find that it does not yet impact the timing model or stochastic parameter estimates significantly compared to analysis performed assuming Gaussian statistics. The methods presented are, however, shown to be of immediate practical use for current European Pulsar Timing Array (EPTA) and International Pulsar Timing Array (IPTA) datasets.

preprint2014arXiv

Radio Source Evolution on Galactic Scales

There is mounting evidence that mechanical radio source feedback is important in galaxy evolution and in order to quantify this feedback, detailed models of radio source evolution are required. We present an extension to current analytic models that encompasses young radio sources with physical sizes on sub-kiloparsec scales. This work builds on an existing young source dynamical model to include radiative losses in a flat environment, and as such, is the best physically-motivated Compact Symmetric Object model to date. Results predict that young radio sources experience significant radiative loss on length scales and spectral scales consistent with observed Compact Steep-Spectrum sources. We include full expressions for the transition to self-similar expansion and present this complete model of radio source evolution from first cocoon formation to end of source lifetime around 10^8 years within the context of a simplified King profile external atmosphere.

preprint2013arXiv

Hyper-efficient model-independent Bayesian method for the analysis of pulsar timing data

A new model independent method is presented for the analysis of pulsar timing data and the estimation of the spectral properties of an isotropic gravitational wave background (GWB). We show that by rephrasing the likelihood we are able to eliminate the most costly aspects of computation normally associated with this type of data analysis. When applied to the International Pulsar Timing Array Mock Data Challenge data sets this results in speedups of approximately 2 to 3 orders of magnitude compared to established methods. We present three applications of the new likelihood. In the low signal to noise regime we sample directly from the power spectrum coefficients of the GWB signal realization. In the high signal to noise regime, where the data can support a large number of coefficients, we sample from the joint probability density of the power spectrum coefficients for the individual pulsars and the GWB signal realization. Critically in both these cases we need make no assumptions about the form of the power spectrum of the GWB, or the individual pulsars. Finally we present a method for characterizing the spatial correlation between pulsars on the sky, making no assumptions about the form of that correlation, and therefore providing the only truly general Bayesian method of confirming a GWB detection from pulsar timing data.

preprint2013arXiv

TempoNest: A Bayesian approach to pulsar timing analysis

A new Bayesian software package for the analysis of pulsar timing data is presented in the form of TempoNest which allows for the robust determination of the non-linear pulsar timing solution simultaneously with a range of additional stochastic parameters. This includes both red spin noise and dispersion measure variations using either power law descriptions of the noise, or through a model-independent method that parameterises the power at individual frequencies in the signal. We use TempoNest to show that at noise levels representative of current datasets in the European Pulsar Timing Array (EPTA) and International Pulsar Timing Array (IPTA) the linear timing model can underestimate the uncertainties of the timing solution by up to an order of magnitude. We also show how to perform Bayesian model selection between different sets of timing model and stochastic parameters, for example, by demonstrating that in the pulsar B1937+21 both the dispersion measure variations and spin noise in the data are optimally modelled by simple power laws. Finally we show that not including the stochastic parameters simultaneously with the timing model can lead to unpredictable variation in the estimated uncertainties, compromising the robustness of the scientific results extracted from such analysis.

preprint2012arXiv

A new connection between the opening angle and the large-scale morphology of extragalactic radio sources

In the case of an initially conical jet, we study the relation between jet collimation by the external pressure and large-scale morphology. We first consider the important length-scales in the problem, and then carry out axisymmetric hydrodynamic simulations that include, for certain parameters, all these length-scales. We find three important scales related to the collimation region: (i) where the sideways ram-pressure equals the external pressure, (ii) where the jet density equals the ambient density, and (iii) where the forward ram-pressure falls below the ambient pressure. These scales are set by the external Mach-number and opening angle of the jet. We demonstrate that the relative magnitudes of these scales determine the collimation, Mach-number, density and morphology of the large scale jet. Based on analysis of the shock structure, we reproduce successfully the morphology of Fanaroff-Riley (FR) class I and II radio sources. Within the framework of the model, an FR I radio source must have a large intrinsic opening angle. Entrainment of ambient gas might also be important. We also show that all FR I sources with radio lobes or similar features must have had an earlier FR II phase.

preprint2012arXiv

Future Science Prospects for AMI

The Arcminute Microkelvin Imager (AMI) is a telescope specifically designed for high sensitivity measurements of low-surface-brightness features at cm-wavelength and has unique, important capabilities. It consists of two interferometer arrays operating over 13.5-18 GHz that image structures on scales of 0.5-10 arcmin with very low systematics. The Small Array (AMI-SA; ten 3.7-m antennas) couples very well to Sunyaev-Zel'dovich features from galaxy clusters and to many Galactic features. The Large Array (AMI-LA; eight 13-m antennas) has a collecting area ten times that of the AMI-SA and longer baselines, crucially allowing the removal of the effects of confusing radio point sources from regions of low surface-brightness, extended emission. Moreover AMI provides fast, deep object surveying and allows monitoring of large numbers of objects. In this White Paper we review the new science - both Galactic and extragalactic - already achieved with AMI and outline the prospects for much more.

preprint2011arXiv

Interaction of Fanaroff-Riley class II jets with a magnetised intra-cluster medium

We present 3-D MHD and synthetic numerical simulations to follow the evolution of randomly magnetized intra-cluster medium plasma under the effects of powerful, light, hypersonic and bipolar jets. We prescribe the cluster magnetic field (CMF) as a Gaussian random field with power law energy spectrum tuned to the expectation for Kolmogorov turbulence. We investigate the power of jets and the viewing angle used for the synthetic Rotation Measure (RM) observations. We find the model radio sources introduce and amplify fluctuations on the RM statistical properties; the average RM and the RM standard deviation are increased by the action of the jets. This may lead to overestimations of the CMFs' strength up to 70%. The effect correlates with the jet power. Jets distort and amplify CMFs especially near the edges of the lobes and the jets' heads. Thus the RM structure functions are flattened at scales comparable to the source size. Jet-produced RM enhancements depend on the orientation of the jet axis to the line of sight and are thus more apparent in quasars than in radio galaxies.

preprint2011arXiv

Interaction of Fanaroff-Riley class II radio jets with a randomly magnetised intra-cluster medium

A combination of three-dimensional (3D) magnetohydrodynamics (MHD) and synthetic numerical simulations are presented to follow the evolution of a randomly magnetised plasma that models the intra-cluster medium (ICM), under the isolated effects of powerful, light, hypersonic and bipolar Fanaroff-Riley class II (FR II) jets. We prescribe the cluster magnetic field (CMF) as a Gaussian random field with a Kolmogorov-like energy spectrum. Both the power of the jets and the viewing angle that is used for the synthetic Rotation Measure (RM) observations are investigated. We find the model radio sources introduce and amplify fluctuations on the RM statistical properties which we analyse as a function of time as well as the viewing angle. The average RM and the RM standard deviation are increased by the action of the jets. Energetics, RM statistics and magnetic power spectral analysis consistently show that the effects also correlate with the jets' power, and that the lightest, fastest jets produce the strongest changes in their environment. We see jets distort and amplify the CMFs especially near the edges of the lobes and the jets' heads. This process leads to a flattening of the RM structure functions at scales comparable to the source size. The edge features we find are similar to ones observed in Hydra A. The results show that jet-produced RM enhancements are more apparent in quasars than in radio galaxies. Globally, jets tend to enhance the RM standard deviation which may lead to overestimations of the CMFs' strength by about 70%. This study means to serve as a pathfinder for the SKA, EVLA and LOFAR to follow the evolution of cosmic magnetic fields.

preprint2011arXiv

Three-dimensional magnetohydrodynamic simulations of the evolution of magnetic fields in Fanaroff-Riley class II radio sources

Radio observations of Fanaroff-Riley class II sources often show correlations between the synchrotron emission and the linear-polarimetric distributions. Magnetic position vectors seem to align with the projected emission of both the radio jets and the sources' edges. Using statistics we study such relation as well as its unknown time evolution via synthetic polarisation maps of model FR II sources formed in 3D-MHD numerical simulations of bipolar, hypersonic and weakly magnetised jets. The magnetic field is initially random with a Kolmogorov power spectrum, everywhere. We investigate the structure and evolution of magnetic fields in the sources as a function of the power of jets and the observational viewing angle. Our synthetic polarisation maps agree with observations, showing B-field vectors which are predominantly aligned with the jet axis, and show that magnetic fields inside sources are shaped by the jets' backflow. Polarimetry is found to correlate with time, the viewing angle and the jet-to-ambient density contrast. The magnetic structure inside thin elongated sources is more uniform than inside more spherical ones. We see jets increase the magnetic energy in cocoons in proportion to the jet velocity and the cocoon width. Filaments in the synthetic emission maps suggest turbulence develops in evolved sources.

preprint2010arXiv

3D-MHD simulations of the evolution of magnetic fields in FR II radio sources

3D-MHD numerical simulations of bipolar, hypersonic, weakly magnetized jets and synthetic synchrotron observations are presented to study the structure and evolution of magnetic fields in FR II radio sources. The magnetic field setup in the jet is initially random. The power of the jets as well as the observational viewing angle are investigated. We find that synthetic polarization maps agree with observations and show that magnetic fields inside the sources are shaped by the jets' backflow. Polarimetry statistics correlates with time, the viewing angle and the jet-to-ambient density contrast. The magnetic structure inside thin elongated sources is more uniform than for ones with fatter cocoons. Jets increase the magnetic energy in cocoons, in proportion to the jet velocity. Both, filaments in synthetic emission maps and 3D magnetic power spectra suggest that turbulence develops in evolved sources.

preprint2010arXiv

Magnetic fields in galaxies: I. Radio disks in local late-type galaxies

We develop an analytical model to follow the cosmological evolution of magnetic fields in disk galaxies. Our assumption is that fields are amplified from a small seed field via magnetohydrodynamical (MHD) turbulence. We further assume that this process is fast compared to other relevant timescales, and occurs principally in the cold disk gas. We follow the turbulent energy density using the Shabala & Alexander (2009) galaxy formation and evolution model. Three processes are important to the turbulent energy budget: infall of cool gas onto the disk and supernova feedback increase the turbulence; while star formation removes gas and hence turbulent energy from the cold gas. Finally, we assume that field energy is continuously transferred from the incoherent random field into an ordered field by differential galactic rotation. Model predictions are compared with observations of local late type galaxies by Fitt & Alexander (1993) and Shabala et al. (2008). The model reproduces observed magnetic field strengths and luminosities in low and intermediate-mass galaxies. These quantities are overpredicted in the most massive hosts, suggesting that inclusion of gas ejection by powerful AGNs is necessary in order to quench gas cooling and reconcile the predicted and observed magnetic field strengths.

preprint2010arXiv

Multipath Parameter Estimation from OFDM Signals in Mobile Channels

We study multipath parameter estimation from orthogonal frequency division multiplex signals transmitted over doubly dispersive mobile radio channels. We are interested in cases where the transmission is long enough to suffer time selectivity, but short enough such that the time variation can be accurately modeled as depending only on per-tap linear phase variations due to Doppler effects. We therefore concentrate on the estimation of the complex gain, delay and Doppler offset of each tap of the multipath channel impulse response. We show that the frequency domain channel coefficients for an entire packet can be expressed as the superimposition of two-dimensional complex sinusoids. The maximum likelihood estimate requires solution of a multidimensional non-linear least squares problem, which is computationally infeasible in practice. We therefore propose a low complexity suboptimal solution based on iterative successive and parallel cancellation. First, initial delay/Doppler estimates are obtained via successive cancellation. These estimates are then refined using an iterative parallel cancellation procedure. We demonstrate via Monte Carlo simulations that the root mean squared error statistics of our estimator are very close to the Cramer-Rao lower bound of a single two-dimensional sinusoid in Gaussian noise.

preprint2009arXiv

Estimation of Magnetic Field Strength in the Turbulent Warm Ionized Medium

We studied Faraday rotation measure (RM) in turbulent media with the rms Mach number of unity, using isothermal, magnetohydrodynamic turbulence simulations. Four cases with different values of initial plasma beta were considered. Our main findings are as follows. (1) There is no strong correlation between the fluctuations of magnetic field strength and gas density. So the magnetic field strength estimated with RM/DM (DM is the dispersion measure) correctly represents the true mean strength of the magnetic field along the line of sight. (2) The frequency distribution of RMs is well fitted to the Gaussian. In addition, there is a good relation between the width of the distribution of RM/$\bar{\rm RM}$ ($\bar{\rm RM}$ is the average value of RMs) and the strength of the regular field along the line of sight; the width is narrower, if the field strength is stronger. We discussed the implications of our findings in the warm ionized medium where the Mach number of turbulent motions is around unity.

preprint2009arXiv

Measurements of the Cosmological Evolution of Magnetic Fields with the Square Kilometre Array

We investigate the potential of the Square Kilometre Array (SKA) for measuring the magnetic fields in clusters of galaxies via Faraday rotation of background polarised sources. [...] We find that about 10 per cent of the sky is covered by a significant extragalactic Faraday screen. Most of it has rotation measures between 10 and 100 rad/m/m. We argue that the cluster centres should have up to about 5000 rad/m/m. We show that the proposed mid frequency aperture array of the SKA as well as the lowest band of the SKA dish array are well suited to make measurements for most of these rotation measure values, typically requiring a signal-to-noise of ten. We calculate the spacing of sources forming a grid for the purpose of measuring foreground rotation measures: it reaches a spacing of 36 arcsec for a 100 hour SKA observation per field. We also calculate the statistics for background RM measurements in clusters of galaxies. We find that a first phase of the SKA would allow us to take stacking experiments out to high redshifts (>1), and provide improved magnetic field structure measurements for individual nearby clusters. The full SKA aperture array would be able to make very detailed magnetic field structure measurements of clusters with more than 100 background sources per cluster up to a redshift of 0.5 and more than 1000 background sources per cluster for nearby clusters, and could for reasonable assumptions about future measurements of electron densities in high redshift clusters constrain the power law index for the magnetic field evolution to better than dm=0.4, if the magnetic field in clusters should follow B ~ (1+z)^m.

preprint2009arXiv

Radio source feedback in galaxy evolution

We present a galaxy evolution model which incorporates a physically motivated implementation of AGN feedback. Intermittent jets inflate cocoons of radio plasma which then expand supersonically, shock heating the ambient gas. The model reproduces observed star formation histories to the highest redshifts for which reliable data exists, as well as the observed galaxy colour bimodality. Intermittent radio source feedback also naturally provides a way of keeping the black hole and spheroid growth in step. We find possible evidence for a top-heavy Initial Mass Function (IMF) for $z>2$, consistent with observations of element abundances, and sub-mm and Lyman break galaxy counts.

preprint2009arXiv

The relationship between star formation rate and radio synchrotron luminosity at 0 < z < 2

We probe the relationship between star formation rate (SFR) and radio synchrotron luminosity in galaxies at 0 < z < 2 within the northern Spitzer Wide-area Infrared Extragalactic survey (SWIRE) fields, in order to investigate some of the assumptions that go into calculating the star formation history of the Universe from deep radio observations. We present new 610-MHz Giant Metrewave Radio Telescope (GMRT) observations of the European Large-Area ISO Survey (ELAIS)-N2 field, and using this data, along with previous GMRT surveys carried out in the ELAIS-N1 and Lockman Hole regions, we construct a sample of galaxies which have redshift and SFR information available from the SWIRE survey. We test whether the local relationship between SFR and radio luminosity is applicable to z = 2 galaxies, and look for evolution in this relationship with both redshift and SFR in order to examine whether the physical processes which lead to synchrotron radiation have remained the same since the peak of star formation in the Universe. We find that the local calibration between radio luminosity and star formation can be successfully applied to radio-selected high-redshift, high-SFR galaxies, although we identify a small number of sources where this may not be the case; these sources show evidence for inaccurate estimations of their SFR, but there may also be some contribution from physical effects such as the recent onset of starburst activity, or suppression of the radio luminosity within these galaxies.

preprint1997arXiv

Evolutionary tracks of FRII sources through the P-D diagram

This is the second in a series of papers presenting an analytical model for the evolution of FRII radio sources. In this paper we evaluate the expected radio emission from a radio source incorporating energy loss processes for the relativistic electrons. By combining these results with our earlier dynamical model we calculate evolutionary tracks through the Power-Linear size diagram. These tracks are in good agreement with the observed distribution of sources in this diagram. The effects of different forms for the evolution of the magnetic field in the cocoon, the redshift of the source, the environment of the source and the defining parameters of the jet are investigated. The evolutionary tracks are found to be insensitive to the assumed form of the magnetic field evolution. Some evidence against protons as a major constituent of the jet material is also found.