A circular statistical method for extracting rotation measures

We propose a new method for the extraction of rotation measures from spectral polarization data. The method is based on maximum likelihood analysis and takes into account the circular nature of the polarization data. The method is unbiased and statistically more efficient than the standard χ ² procedure. We also find that the method is computationally much more convenient than the standard χ ² procedure if the number of data points is very large. We find that for most sources the method gives results very close to the standard χ ² minimization procedure. We give results for all the cases for which the method gives significantly different results from χ ² minimization. We also make a λ ³ fit to the data in order to extract non-Faraday rotation behaviour for those sources for which a large number of data points are available.     

A Faraday rotation template for the Galactic sky

Using a set of compilations of measurements for extragalactic radio sources, we construct all-sky maps of the Faraday rotation produced by the Galactic magnetic field. In order to generate the maps, we treat the radio source positions as a kind of 'mask' and construct combinations of spherical harmonic modes that are orthogonal on the masked sky. As long as relatively small multipoles are used, the resulting maps are quite stable to changes in the selection criteria for the sources, and show clearly the structure of the local Galactic magnetic field. We also suggest the use of these maps as templates for cosmic microwave background (CMB) foreground analysis, illustrating the idea with a cross-correlation analysis between the Wilkinson Microwave Anisotropy Probe ( WMAP ) data and our maps. We find a significant cross-correlation, indicating the presence of a significant residual contamination.     

New pulsar rotation measures and the Galactic magnetic field

We measured a sample of 150 pulsar rotation measures (RMs) using the 20-cm receiver of the Parkes 64-m radio telescope. 46 of the pulsars in our sample have not had their RM values previously published, whereas 104 pulsar RMs have been revised. We used a novel quadratic fitting algorithm to obtain an accurate RM from the calibrated polarization profiles recorded across 256 MHz of receiver bandwidth. The new data are used in conjunction with previously known dispersion measures and the NE2001 electron-density model to study models of the direction and magnitude of the Galactic magnetic field.     

Figure rotation of dark haloes in cold dark matter simulations

We investigate the figure rotation of dark matter haloes identified in Λ cold dark matter (CDM) simulations. We find that when strict criteria are used to select suitable haloes for study, five of the 222 haloes identified in our   z = 0  simulation output undergo coherent figure rotation over a  5 h ⁻¹ Gyr  period. We discuss the effects of varying the selection criteria and find that pattern speeds for a much larger fraction of the haloes can be measured when the criteria are relaxed. Pattern speeds measured over a  1 h ⁻¹ Gyr  period follow a lognormal distribution, centred at  Ω_p= 0.2 h rad Gyr⁻¹  with a maximum value of 0.94 h rad Gyr⁻¹. Over a  5 h ⁻¹ Gyr  period, the average pattern speed of a halo is about  0.1 h rad Gyr⁻¹  and the largest pattern speed found is  0.24 h rad Gyr⁻¹  . Less than half of the selected haloes showed alignment between their figure rotation axis and minor axis, the exact fraction being somewhat dependent on how one defines a halo. While the pattern speeds observed are lower than those generally thought capable of causing spiral structure, we note that coherent figure rotation is found over very long periods and argue that further simulations would be required before strong conclusions about spiral structure in all galaxies could be drawn. We find no correlation between halo properties such as total mass and the pattern speed.     

A method for nuclear-component deconvolution in active galaxies

A new method for deconvolution of the nuclear component of active galaxies is presented, along with some initial results with several nuclear-type objects. Unlike the majority of procedures currently in use, this method is based on the following steps: (a) derivation of seeing-free bulge and disc parameters, (b) construction of a two-dimensional model of the galaxy, (c) derivation of a point-spread function (PSF) from star profiles and convolution of the bulge+disc model with the obtained PSF and (d) extraction of the seeing-convolved model radial profile and subtraction from the observed one. The residual gives the nuclear component. Only standard image processing packages are used.     

Objective classification of galaxy spectra using the information bottleneck method

A new method for classification of galaxy spectra is presented, based on a recently introduced information theoretical principle, the information bottleneck . For any desired number of classes, galaxies are classified such that the information content about the spectra is maximally preserved. The result is classes of galaxies with similar spectra, where the similarity is determined via a measure of information. We apply our method to ∼6000 galaxy spectra from the ongoing 2dF redshift survey, and a mock-2dF catalogue produced by a cold dark matter (CDM) based semi-analytic model of galaxy formation. We find a good match between the mean spectra of the classes found in the data and in the models. For the mock catalogue, we find that the classes produced by our algorithm form an intuitively sensible sequence in terms of physical properties such as colour, star formation activity, morphology, and internal velocity dispersion. We also show the correlation of the classes with the projections resulting from a principal component analysis.     

Magnetic and optical spiral arms in the galaxy NGC 6946

The spiral pattern in the nearby spiral galaxy NGC 6946 has been studied using the wavelet transformation technique, applied to galaxy images in polarized and total non-thermal radio emission at λλ 3.5 and 6.2 cm, in broadband red light, in the λ 21.1 cm H  i line and in the optical Hα line. Well-defined, continuous spiral arms are visible in polarized radio emission and red light, where we can isolate a multi-armed pattern in the range of galactocentric distances 1.5–12 kpc, consisting of four long arms and one short spiral segment. The 'magnetic arms' (visible in polarized radio emission) are localized almost precisely between the optical arms. Each magnetic arm is similar in length and pitch angle to the preceding optical arm (in the sense of galactic rotation) and can be regarded as its phase-shifted image. Even details like a bifurcation of an optical arm have their phase-shifted counterparts in the magnetic arms. The average relative amplitude of the optical spiral arms (the stellar density excess over the azimuthal average) grows with galactocentric radius up to 0.3–0.7 at r ≃5 kpc, decreases by a factor of two at r =5–6 kpc and remains low at 0.2–0.3 in the outer parts of the galaxy. By contrast, the magnetic arms have a constant average relative amplitude (the excess in the regular magnetic field strength over the azimuthal average) of 0.3–0.6 in a wide radial range r =1.5–12 kpc. We briefly discuss implications of our findings for theories of galactic magnetic fields.     

A new method for determining the sensitivity of X-ray imaging observations and the X-ray number counts

We present a new method for determining the sensitivity of X-ray imaging observations, which correctly accounts for the observational biases that affect the probability of detecting a source of a given X-ray flux, without the need to perform a large number of time-consuming simulations. We use this new technique to estimate the X-ray source counts in different spectral bands (0.5–2, 0.5–10, 2–10 and 5–10 keV) by combining deep pencil-beam and shallow wide-area Chandra observations. The sample has a total of 6295 unique sources over an area of  11.8 deg²  and is the largest used to date to determine the X-ray number counts. We determine, for the first time, the break flux in the 5–10 keV band, in the case of a double power-law source count distribution. We also find an upturn in the 0.5–2 keV counts at fluxes below about  6 × 10⁻¹⁷ erg s⁻¹ cm⁻²  . We show that this can be explained by the emergence of normal star-forming galaxies which dominate the X-ray population at faint fluxes. The fraction of the diffuse X-ray background resolved into point sources at different spectral bands is also estimated. It is argued that a single population of Compton thick active galactic nuclei (AGN) cannot be responsible for the entire unresolved X-ray background in the energy range 2–10 keV.     

An operator perturbation method of polarized line transfer V. Diagnosis of solar weak magnetic fields

We present an application of the PALI (Polarized Approximate Lambda Iteration) method to the resonance scattering in spectral lines formed in the presence of weak magnetic fields. The method is based on an operator perturbation approach, and can efficiently give solutions for oriented vector magnetic fields in the solar atmosphere.     

Quick Detection System for Planck satellite

Q uick D etection S ystem or qds is a software package that has been developed for detecting point sources in the Planck satellite data as soon as the data become available, a few days after transmission to the Earth. Point sources are detected by filtering the data with a filter defined by the Mexican hat wavelet. An alert is generated on those detections that are found to be interesting, such as prominent flaring, according to the criteria specified to the software. The goal is to detect spectral or flux variability in active galactic nuclei so that instant multifrequency follow-up observations with other instruments could be arranged to study the interesting behaviour.     

An MHD gadget for cosmological simulations

Various radio observations have shown that the hot atmospheres of galaxy clusters are magnetized. However, our understanding of the origin of these magnetic fields, their implications on structure formation and their interplay with the dynamics of the cluster atmosphere, especially in the centres of galaxy clusters, is still very limited. In preparation for the upcoming new generation of radio telescopes (like Expanded Very Large Array, Low Wavelength Array, Low Frequency Array and Square Kilometer Array), a huge effort is being made to learn more about cosmological magnetic fields from the observational perspective. Here we present the implementation of magnetohydrodynamics (MHD) in the cosmological smoothed particle hydrodynamics (SPH) code gadget . We discuss the details of the implementation and various schemes to suppress numerical instabilities as well as regularization schemes, in the context of cosmological simulations. The performance of the SPH–MHD code is demonstrated in various one- and two-dimensional test problems, which we performed with a fully, three-dimensional set-up to test the code under realistic circumstances. Comparing solutions obtained using athena , we find excellent agreement with our SPH–MHD implementation. Finally, we apply our SPH–MHD implementation to galaxy cluster formation within a large, cosmological box. Performing a resolution study we demonstrate the robustness of the predicted shape of the magnetic field profiles in galaxy clusters, which is in good agreement with previous studies.     

A SCUBA scan-map of the Hubble Deep Field: measuring the bright end of the submillimetre source counts

Using the 850-μm SCUBA camera on the JCMT and a scanning technique different from other sub-mm surveys, we have obtained a 125 arcmin² map centred on the Hubble Deep Field . The 1 σ sensitivity to point sources is roughly 3 mJy and thus our map probes the brighter end of the sub-mm source counts. We find six sources with a flux greater than about 12 mJy (>4 σ ) and, after a careful accounting of incompleteness and flux bias, estimate the integrated density of bright sources     (68 per cent confidence bounds).