Search

(abridged) We run a Faraday structure determination data challenge to benchmark the currently available algorithms including Faraday synthesis (previously called RM synthesis in the literature), wavelet, compressive sampling and $QU$-fitting. The frequency set is similar to POSSUM/GALFACTS with a 300 MHz bandwidth from 1.1 to 1.4 GHz. We define three figures of merit motivated by the underlying science: a) an average RM weighted by polarized intensity, RMwtd, b) the separation $\Delta\phi$ of two Faraday components and c) the reduced chi-squared. Based on the current test data of signal to noise ratio of about 32, we find that: (1) When only one Faraday thin component is present, most methods perform as expected, with occasional failures where two components are incorrectly found; (2) For two Faraday thin components, QU-fitting routines perform the best, with errors close to the theoretical ones for RMwtd, but with significantly higher errors for $\Delta\phi$. All other methods including standard Faraday synthesis frequently identify only one component when $\Delta\phi$ is below or near the width of the Faraday point spread function; (3) No methods, as currently implemented, work well for Faraday thick components due to the narrow bandwidth; (4) There exist combinations of two Faraday components which produce a large range of acceptable fits and hence large uncertainties in the derived single RMs; in these cases, different RMs lead to the same Q, U behavior, so no method can recover a unique input model., Comment: 15 pages, 8 figures, 5 tables, accepted by AJ. Some figures have been degraded, for a full resolution version, see http://physics.usyd.edu.au/~xhsun/datachallengeI.pdf

We present RESOLVE, a new algorithm for radio aperture synthesis imaging of extended and diffuse emission in total intensity. The algorithm is derived using Bayesian statistical inference techniques, estimating the surface brightness in the sky assuming a priori log-normal statistics. RESOLVE not only estimates the measured sky brightness in total intensity, but also its spatial correlation structure, which is used to guide the algorithm to an optimal reconstruction of extended and diffuse sources. For a radio interferometer, it succeeds in deconvolving the effects of the instrumental point spread function during this process. Additionally, RESOLVE provides a map with an uncertainty estimate of the reconstructed surface brightness. Furthermore, with RESOLVE we introduce a new, optimal visibility weighting scheme that can be viewed as an extension to robust weighting. In tests using simulated observations, the algorithm shows improved performance against two standard imaging approaches for extended sources, Multiscale-CLEAN and the Maximum Entropy Method., Comment: 22 pages, 8 figures, submitted to Astronomy & Astrophysics

We develop a method to infer log-normal random fields from measurement data affected by Gaussian noise. The log-normal model is well suited to describe strictly positive signals with fluctuations whose amplitude varies over several orders of magnitude. We use the formalism of minimum Gibbs free energy to derive an algorithm that uses the signal's correlation structure to regularize the reconstruction. The correlation structure, described by the signal's power spectrum, is thereby reconstructed from the same data set. We show that the minimization of the Gibbs free energy, corresponding to a Gaussian approximation to the posterior marginalized over the power spectrum, is equivalent to the empirical Bayes ansatz, in which the power spectrum is fixed to its maximum a posteriori value. We further introduce a prior for the power spectrum that enforces spectral smoothness. The appropriateness of this prior in different scenarios is discussed and its effects on the reconstruction's results are demonstrated. We validate the performance of our reconstruction algorithm in a series of one- and two-dimensional test cases with varying degrees of non-linearity and different noise levels., Comment: 18 pages, 12 figures

We present a sequence of five deep observations of SS433 made over the summer of 2007 using the VLA in the A configuration at 5 and 8 GHz. In this paper we study the brightness profiles of the jets and their time evolution. We also examine the spectral index distribution in the source. We find (as previously reported from the analysis of a single earlier image) that the profiles of the east and west jets are remarkably similar if projection and Doppler beaming are taken into account. The sequence of five images allows us to disentangle the evolution of brightness of individual pieces of jet from the variations of jet power originating at the core. We find that the brightness of each piece of the jet fades as an exponential function of age (or distance from the core), exp(-tau/tau'), where tau is the age at emission and tau' = 55.9 +- 1.7 days. This evolutionary model describes both the east and west jets equally well. There is also significant variation (by a factor of at least five) in jet power with birth epoch, with the east and west jets varying in synchrony. The lack of deceleration between the scale of the optical Balmer line emission (10^15 cm) and that of the radio emission (10^17 cm) requires that the jet material is much denser than its surroundings. We find that the density ratio must exceed 300:1., Comment: 26 pages, 13 Figures, Accepted for publication in the Astrophysical Journal