Issue |
A&A
Volume 507, Number 2, November IV 2009
|
|
---|---|---|
Page(s) | 1087 - 1105 | |
Section | Astronomical instrumentation | |
DOI | https://doi.org/10.1051/0004-6361/200912539 | |
Published online | 15 September 2009 |
A&A 507, 1087-1105 (2009)
Simulated square kilometre array maps from Galactic 3D-emission models
X. H. Sun1,2 - W. Reich2
1 - National Astronomical Observatories, CAS,
Jia 20, Datun Road, Chaoyang District, Beijing 100012, PR China
2 - Max-Planck-Institut für Radioastronomie, Auf dem Hügel 69,
53121 Bonn, Germany
Received 20 May 2009 / Accepted 13 August 2009
Abstract
Context. Planning of the Square Kilometre Array (SKA)
requires simulations of the expected sky emission at arcsec angular
resolution to evaluate its scientific potential, to constrain its
technical realization in the best possible way, and to guide the
observing strategy.
Aims. We simulate high-resolution total intensity, polarization,
and rotation measure (RM) maps of selected fields based on our recent
global 3D-model of Galactic emission.
Methods. Simulations of diffuse Galactic emission were conducted using the HAMMURABI
code modified for arcsec angular resolution patches towards various
Galactic directions. The random magnetic field components are set to
follow a Kolmogorov-like power-law spectrum. We analysed the simulated
maps in terms of their probability density functions (PDFs) and
structure functions.
Results. We present maps for various Galactic longitudes and
latitudes at 1.4 GHz, which is the frequency where deep SKA
surveys are proposed. The maps are about
in size and have an angular resolution of about
.
Total intensity emission is smoother in the plane than at high
latitudes because of the different contributions from the regular and
random magnetic field. The high-latitude fields show more extended
polarized emission and RM structures than those in the plane, where
patchy emission structures dominate on very small scales. The RM PDFs
in the plane are close to Gaussians, but clearly deviate from that at
high latitudes. The RM structure functions show smaller amplitudes and
steeper slopes towards high latitudes. These results emerge from much
more turbulent cells being passed through by the line-of-sights in the
plane. Although the simulated random magnetic field components
distribute in 3D, the magnetic field spectrum extracted from the
structure functions of RMs conforms to 2D in the plane and approaches
3D at high latitudes. This is partly related to the outer scale of the
turbulent magnetic field, but mainly to the different lengths of the
lines-of-sight.
Conclusions. The significant scatter of the simulated RM
distributions agrees with the large scatter of observed RMs of pulsars
and extragalactic sources in the Galactic plane and also at high
Galactic latitudes. A very dense grid of RMs from extragalactic sources
is required to trace and separate Galactic RM foreground fluctuations.
Even at high latitudes, total intensity and polarized emission is
highly structured, which will contaminate sensitive high-resolution
extragalactic observations with the SKA and have to be separated in an
appropriate way.
Key words: magnetic fields - polarization - radiation mechanisms: non-thermal - ISM: magnetic field - ISM: structure
1 Introduction
Intensive radio emission from the Galaxy originates from ionized thermal gas, cosmic rays, and magnetic fields. All-sky surveys yield an overview of the Galactic structures; however, to reveal distributions of the Galactic emission components, modelling is required because of our position inside the Galaxy. Galactic emission is a severe foreground contamination for observations of fluctuations of the cosmic microwave background radiation and of all kinds of distant, faint extragalactic objects, which will be investigated by the Square Kilometre Array (SKA) with unprecedented sensitivity. Therefore it is crucial to properly account for Galactic emission in all directions of the sky. Unfortunately, most of the properties of these constituents of the Galaxy from large scales to small scales are insufficiently constrained up to now.
The large-scale Galactic magnetic field is of particular interest. It has been
intensively studied by using rotation measures (RMs) of pulsars
(e.g. Han et al. 2006; Noutsos et al. 2008) and extragalactic sources (EGSs) in the Galactic
plane (e.g. Brown et al. 2007,2003). Pulsar dispersion measures (DMs) are used to
find their distances based on a thermal gas distribution model. Recent pulsar
results indicate several magnetic field reversals within the Galaxy
(e.g. Han et al. 2006; Noutsos et al. 2008), which is difficult to understand in terms of a
simple axisymmetric or bisymmetric magnetic field configuration. A denser
RM grid is needed to derive a clear large-scale pattern, because complications
arise from local magnetic field disturbances or excessive thermal gas
contributions to the observed pulsar DMs and RMs (Mitra et al. 2003). Such effects
must be separated to derive a reliable large-scale magnetic field configuration.
Moreover, to decompose the observed RM information into a magnetic field
contribution, detailed knowledge of the thermal gas distribution is required.
The currently widely accepted distribution of Galactic diffuse thermal gas, the
NE2001 model, as well as its filling factor modelled from pulsar and
H
observations was presented by Cordes & Lazio (2002) and Berkhuijsen et al. (2006),
respectively. Galactic magnetic field models based on pulsar RM observations
have to agree with the observed RMs of extragalactic sources shining through
the entire Galaxy. A firm picture of the large-scale magnetic
field needs to account for all the effects above.
Recently, Sun et al. (2008) have derived a Galactic 3D-emission model that
reproduces all-sky surveys over a wide frequency range and also the all-sky RM
distribution of EGSs. The simulated maps have an angular resolution of about
15
,
which is suitable for comparing them with all-sky maps having
a lower angular resolution. Based on these global results for the Galaxy,
simulated arcsec angular resolution maps are the next step in getting an
idea of what might be observed with the SKA and other future high-resolution
radio telescopes. It is important to understand the origin of the observed
structures at high angular resolution to estimate the contamination of
extragalactic observations by Galactic small-scale structures as the foreground.
One of the key science projects for the SKA is ``Cosmic Magnetism'',
where RMs of about
polarized EGSs are proposed to be observed in a region of 104 deg2 (Beck & Gaensler 2004). The survey is proposed for
1.4 GHz with a resolution of about
.
Simulations of the expected
number density of faint sources and their polarization properties from SKA
observations were made by e.g. Wilman et al. (2008). A feasibility study of the
decomposition of the source RMs from the Galactic foreground RMs and its
removal need high-resolution Galactic RM maps. This is the basic motivation for
extending our simulations to selected Galactic patches with arcsec angular
resolution.
Fluctuations in the components of the interstellar medium (ISM) dominate the emission structures at high angular resolution. The thermal electron density distribution has been found to follow a Kolmogorov-like power law from about 108 m to about 1018 m (Armstrong et al. 1995). The fluctuations of the Galactic magnetic field are mainly studied based on the analysis of the structure function of RMs, where an additional effort is needed to decouple the magnetic field and the electron density fluctuations (Minter & Spangler 1996). The outer scale of turbulence was reported to be in the range of about 100 pc (Haverkorn et al. 2008; Armstrong et al. 1995) down to a few pc (Haverkorn et al. 2008; Minter & Spangler 1996).
We summarise our global Galactic 3D-modelling in Sect. 2. The high angular resolution simulation method and the treatment of small-scale magnetic fields are described in Sect. 3. Results and analysis of the simulated RM maps are presented in Sect. 4, including a discussion of the consequences of the simulated RM fluctuations. In Sect. 5 we discuss the results of the simulated total intensity and polarization maps. Some general conclusions are given in Sect. 6.
2 A new Galactic 3D-emission model
Sun et al. (2008) present a new Galactic 3D-emission model, which agrees with a wide range of observations:
- galactic synchrotron radiation, which is best represented by the 408 MHz all-sky survey (Haslam et al. 1982) in total intensity and by the WMAP 22.8 GHz all-sky survey in polarized intensity (Page et al. 2007);
- optically thin thermal emission by the WMAP 22.8 GHz free-free emission template (Hinshaw et al. 2007), and optically thick thermal emission causing absorption, which is reflected in a spectral flattening at low frequencies along the Galactic plane;
- the all-sky 1420 MHz polarization survey combined from northern (Wolleben et al. 2006) and southern sky data (Testori et al. 2008) constraining depolarization properties of the interstellar medium;
- RM data from EGSs along the Galactic plane taken from the Canadian (CGPS) and the Southern Galactic Plane Survey (SGPS) (Brown et al. 2007,2003), and EGSs out of the plane from Han et al. (1999), including preliminary RMs from the new Effelsberg L-band survey (Han et al., in prep.) and RMs from the CGPS high latitude extension (Brown et al., in prep.).
The quantitative model parameters of the three constituents of the Galaxy are
those of the NE2001 thermal electron density model, an ASS magnetic field
configuration with a local strength of 2 G, asymmetric toroidal halo
fields with respect to the plane with a maximum strength of 10
G,
isotropic and homogeneous random magnetic fields with a strength of 3
G,
and a cosmic-ray electron distribution truncated at 1 kpc above and below the
plane. The details are given by Sun et al. (2008).
Although this Galactic 3D-model agrees with the observed all-sky maps,
some of its parameters seem unrealistic. To fit the observed high-latitude RMs,
a strong regular halo magnetic field is required, because the scale height of
the thermal electrons of the NE2001 model is just about 1 kpc. The strong halo
field in turn requires a cosmic-ray electron distribution limit of about 1 kpc
to avoid excessive high-latitude synchrotron emission. These parameters are
unrealistic for both the halo field strength (Moss & Sokoloff 2008) and the distribution
of cosmic-ray electrons. Sun et al. (2008) note that, by increasing the scale
height of the thermal electron density to about 2 kpc, the halo magnetic field
strength reduces to about 2 G and the truncation of the cosmic-ray
distribution turns obsolete. Recently, Gaensler et al. (2008) have again critically analysed the
thermal electron scale height and quote a revised value of 1.8 kpc, very close
to the proposed 2 kpc by Sun et al. (2008). This makes the Galactic properties
similar to those observed for nearby galaxies. It should be noted that the
simulated maps are not affected by these modifications of the parameters.
3 High-resolution simulations
3.1 The general simulation concept
The simulations presented in this paper make use of the HAMMURABI code
(Waelkens et al. 2009; Waelkens 2005), where the sky is pixelized according to the HEALPix
scheme (Górski et al. 2005). Each pixel contains information from a cone-shaped
volume centred at the observer. Each surface element is quadrilateral and
covers the same area. The number of pixels
of a sphere
are calculated by
,
where
is
an integer power of 2. Therefore d
,
the linear size or the angular
resolution, can be calculated as
For each pixel of the map, the Stokes parameters (I, U, and Q) and the RM are integrated in the radial direction from the observer up to a maximal distance



3.2 Zoom-in for specified patches
According to Eq. (1) the possible resolutions depend on
.
To achieve arcsec resolution and to optimize the amount of
computation time, we use
of 131 072, which corresponds to a
resolution of about
.
The pixel number for an all-sky map then
amounts to about
,
which is far beyond any computer memory
capacity today. Therefore high-resolution simulations are only possible for a
small patch of sky. Thus we adapted the HAMMURABI code to be able to
simulate parts of the sky.
In the HAMMURABI code, the NEST ordering (Górski et al. 2005)
facilitates the zoom-in of specified patches as explained below. For
the entire sky is split into 12 patches as shown in
Fig. 1. They are numbered from 0 to 11. For
,
each
pixel is further split into four child pixels, also shown in Fig. 1
for patch ``5''. This patch is centred at a longitude of
in the plane with a diagonal size of about
.
Subsequently, when
is increased by 2, each pixel is further split into four pixels.
In this scheme one can determine a target patch and make proper splits
depending on the resolution. To reach a smaller patch size and higher angular
resolution a much higher value is required for
.
![]() |
Figure 1: The HEALPix pixel scheme (Górski et al. 2005) for Galactic coordinates used in this paper. |
Open with DEXTER |
3.3 Realisation of the random magnetic field component
It is generally accepted that the power spectrum of the random magnetic field
component should follow a power law as
.
Here
is the wave vector and l the spatial scale. For a Kolmogorov-like
turbulence spectrum, the spectral index is
,
which is used in all the present simulations. Independent fluctuations of the electron
density are not included in the simulations. However, fluctuations of the
electron density are also expected to follow a Kolmogorov-like spectrum
(Armstrong et al. 1995).
We introduced the coupling of electron density and random magnetic
fields in our all-sky simulations to model the observed depolarization (Sun et al. 2008), which results in Kolmogorov-like fluctuations for both components.
Random magnetic field components realised in the simulations presented by
Sun et al. (2008) are generated in situ for every volume unit by Gaussian
random numbers. To obtain the random
magnetic field components following a given power-law spectrum, an inverse
Fourier transformation must be performed. All available fast Fourier
transformation (FFT) algorithms are based on a Cartesian coordinate system. We
therefore setup a cubic box with a size of L3 and divided it into N3
small cubes with a common volume of
.
Within the box we constructed a
Cartesian coordinate system and perform the FFT to derive the random field
strength and direction in each small cube. Here L corresponds to the outer
scale of the turbulence and
corresponds to the inner scale.
An important constraint for
is that
.
This
ensures that several neighbouring line-of-sights do not pass the same small
cubes. Otherwise spurious stripes are seen in the maps. For r of about 1 kpc
and
of
,
the size of a small cube should not be larger than about 0.01 pc. If our Galaxy is treated as one box, meaning L = 40 kpc, the total number of cubes is about
.
Again there is no way to process this huge amount of cubes by current computing facilities.
In our simulation, we devised a scheme to generate the random magnetic field
components in an adequate way. We first conducted an FFT to create a magnetic
field distribution in a box. We then rotate this box to acquire new boxes,
which are glued together to fill the entire Galaxy. The box rotation angles are
integer multiples of ,
which avoids cutting boxes. The turbulence status of the Galaxy is realised through this scheme. This process is
explained in Appendix A.
Actually the size of the box cannot be too small, or else there are too many
boxes to handle. On the other hand, L cannot be too large to get a
reasonable number of cubes of about 0.01 pc in size. As a compromise, we take
a size of 10 pc for the box and of 0.0125 pc for the cubes. Although
0.0125 pc is much larger than the inner scale of the electron density
fluctuations, which are of the order of 108 cm (e.g. Armstrong et al. 1995), it adapts to the angular resolution of our study. In total
there are about
cubes, requiring a memory of about 12 giga-byte for the random
magnetic field components when stored as double precision numbers. It
is thus obvious that the number of cubes is a compromise between high
resolution, which requires small cube sizes, hence a very large number
of cubes, and the available computer memory, which is thus able to
handle a limited number of cubes.
![]() |
Figure 2: 2D random magnetic field distribution in a slice parallel to one surface of a 10 pc box. The upper left panel shows magnetic field components following a power law index of -11/3, while the upper right panel shows random field components realised by random Gaussian numbers corresponding to a spectral index of zero. Their corresponding power spectra are shown in the lower panels. |
Open with DEXTER |
In Fig. 2, we show the 2D-distribution of random magnetic field
components for a slice parallel to one surface of a 10 pc box. For
comparison, we show the random magnetic fields generated in situ in each
grid of the same plane as random Gaussian number, which is similar to the one used by Sun et al. (2008). The power spectra of the magnetic fields are also shown in Fig. 2.
It is apparent that the distribution of random magnetic fields
following a power-law spectrum shows coherent structures of various
extents, whereas the Gaussian random magnetic field distribution is
featureless. The spectral index of
for Kolomogorov magnetic field distributions and
for Gaussian field distributions agree with the expectations and prove
that the random magnetic field components are properly distributed in
the box.
4 Modelling results on rotation measures
Except for the random magnetic fields, the simulations presented below use the
same Galactic 3D-emission models as derived by Sun et al. (2008), where the
``ASS+RING'' model for the disk magnetic field is selected here. The random
magnetic fields are introduced as described above and follow a Kolmogorov-like
spectrum. The strength of the random fields is scaled to 3 G, which is the
same as in Sun et al. (2008). All simulations are made for 1.4 GHz
observations, as this is the proposed SKA RM survey frequency. As mentioned
above, the simulated patch size is limited by available computer memory, so
that just relatively small patches of about
in size are simulated at a
resolution of
.
This is compatible with the expected SKA angular
resolution.
Simulations are presented for only a handful of patches, because each run takes
a very long computation time. We selected the first patch around
,
which is in the direction of the
DRAO Planck Deep Fields (Taylor et al. 2007), where soon high-resolution
continuum and polarization data will be available. This patch is named MN
standing for Medium-latitude-North. More patches with roughly the same
longitude, but with different latitudes, are also selected. These are
labelled HN (High-latitude-North), PNAC1 (Plane-North-Anti-Centre 1), PSAC
(Plane-South-Anti-Centre), MS (Medium-latitude-South) and HS
(High-latitude-South). To compare with PNAC1, three additional patches in the
plane, PNAC2 (Plane-North-Anti-Centre 2), PNC1, and PNC2
(Plane-North-Centre 1 and 2) are included. The centre coordinates of these
patches are listed in Table 1. We cannot set the patch
coordinates to exactly the same longitudes or latitudes because of HEALPix
constraints (Górski et al. 2005). In Table 1 the maximal integral
lengths or the length of lines-of-sight (
)
and the integral steps
(
)
are also listed. The integral length is determined by taking
the cosmic-ray scale height of about 1 kpc and the disk cutoff radius
of 17.5 kpc into account for the electron density in the NE2001 model (Cordes & Lazio 2002).
Integrating further out results in negligible contributions. For the
patches with large integral length, the step size is larger than for
those with a small length to limit the computation time. The selected
patches allow study of variations in emission and RM structures
related to turbulent magnetic-field components versus Galactic
longitudes and latitudes. Variations in the structure functions for RMs
from EGSs as a function of Galactic direction have already been noted
by Sun & Han (2004) for angular scales of several degrees.
Table 1: Simulated patches.
All simulations were made by using the computer cluster of the Max-Planck-Institut für Radioastronomie operated by the VLBI group.
4.1 Overview of the results
As described by Sun et al. (2008) the HAMMURABI simulations result in
Stokes I, U, and Q maps at the specified frequency, as well as the RM
maps. Subsequently, the polarization angle (PA) and polarized intensity
(PI) maps were calculated according to
and
.
In this paper, we do not discuss the separation process of total and polarized intensities of EGSs from the foreground, which requires adding polarized EGSs and noise to the simulated maps. The simulated RMs represent the contribution of the whole Galaxy to the RM from an EGS in case the source size does not exceed the pixel size of the simulated map. Thus, the RM maps are directly related to observations. Below we focus on analysing RM maps, which are directly connected to the planned SKA RM survey. All simulated RM maps are shown in Fig. 3. The other maps are shown in Figs. B.1-B.4, and B.5 in Appendix B.
The RM maps clearly show structures ranging from small to large scales. The variation of the structures versus latitudes is also evident. Below, some statistical methods are applied to analyse these structures.
![]() |
Figure 3:
RM maps for all simulated patches ordered from high to low latitudes ( top to bottom)
in gnonomic projection. Their designation and Galactic coordinates are
quoted on top of each individual panel (also listed in Table 1).
Different simulations for the same patch and for patches at the same
latitude, but different longitudes, are shown in one row (see text for
details). Each field has about
|
Open with DEXTER |
4.2 Statistical description of the simulations
The rotation measure, RM, is defined as
where K is a constant,




4.2.1 The probability distribution function
We use the probability distribution function (PDF) to describe the structure
properties of a map as discussed by Mizeva et al. (2007). For the RM distribution,
quantifies the relative number of
pixels falling between
and
.
A Gaussian
shape of the PDF indicates that no extended features are present in the map.
To measure the deviation from a Gaussian distribution, two more parameters,
skewness and kurtosis, are used (e.g. Kowal et al. 2007),
![]() |
(3) |
where


The RM PDFs for all simulated patched are shown in Figs. 4, 6, 7, and 8 (upper panels). Their average values, together with variance, skewness, and kurtosis are listed in Table 2.
4.2.2 The structure function
For turbulence studies, it is customary to investigate the power spectrum of
the magnetic field, which can be retrieved from the power spectrum,
autocorrelation function, or structure function of RMs
(e.g. Enßlin & Vogt 2003; Simonetti et al. 1984; Minter & Spangler 1996). Observationally, uneven sampled RMs of
EGSs are always obtained. In this case the structure function is the most suitable (Sun & Han 2004). Since RM maps are available as simulation results, the
autocorrelation function is also investigated here as a crosscheck for the
structure function results. The structure function (
)
is defined as
![]() |
(4) |
where


It would take a huge amount of computing time to obtain the structure
functions from all pixels of our simulated maps. We simplified this by randomly selecting
pixels to conduct the calculation for structure functions.
These pixels are randomly distributed across the maps. To accurately assess
the structure functions for small angular separations, the number
should
be as large as possible. As a compromise we selected
pixels
from the maps. The structure functions are shown in Figs. 4,
6, 7, and 8 (lower panels).
Table 2: Statistical parameters for the simulated RM patches.
Table 3: Fitting results for structure functions and autocorrelation functions of RMs.
The theoretical deduction of how structure functions are related to the
power spectrum of magnetic fields has been given by Minter & Spangler (1996), where the
structure function of RMs is found to follow a power law as
with A a constant and the
spectral index m related to the power-law index
of P(k) as
.
This relation holds for the inertial range, i.e.,
for the scale
,
where
and
are
inner and outer scales, respectively. The angular inertial scale can thus be
defined as
,
where
and
.
In principle the inertial range
can be estimated by checking where the linear fit of the structure functions
for a logarithmic scale fails. For the structure functions from our simulated
results, we are able to infer the outer scale, but the inner scale is smaller
than the pixel size. Below we name the outer angular scale
as the transition angle. The fit results for structure
functions are listed in Table 3.
Enßlin & Vogt (2003) have noted that the very inner part of the autocorrelation
function of RMs (
)
can be written as
,
where B is a
constant, and
is the spectral index of the magnetic energy
spectrum. This implies that the power-law index of the structure function
m and the power-law index of the autocorrelation function s should be
equal. We obtain the autocorrelation functions by a Fourier transform of the
RM power spectra and show them in Fig. 5. The RM fit results for each
patch are listed in Table 3. As can be seen, the spectral
indices of the autocorrelation functions do not deviate very much from that of
the structure functions. We explain the differences mainly by computational
effects when obtaining the autocorrelation functions. Unwanted components might
be introduced during the pixelization of the maps when Fourier transformed.
Therefore we are confident that the results based on the two methods are
basically consistent with each other.
![]() |
Figure 4: PDFs and structure functions of the patches HN and MN at high and medium latitudes. See Tables 1 and 2 for details of the patch parameters. |
Open with DEXTER |
The spectral index of 5/3 is expected from the power-law fit of the RM structure functions for a Kolmogorov-like turbulent magnetic field. This is, however, not seen in Table 3. Possible explanations are discussed in Sect. 4.7.1.
4.3 Gaussian versus Kolmogorov-like random magnetic fields
In the simulations presented by Sun et al. (2008), the random magnetic field components were realised in situ for each volume unit by Gaussian random numbers. The random fields created that way are called ``Gaussian fields'' below. In that scheme, the simulated RMs are entirely uncorrelated between the map pixels, which implies that no small-scale structures are expected. When random magnetic fields with a Kolmogorov-like spectrum are generated, as described in Sect. 3.3, many turbulent boxes are placed along the line-of-sight. For arcsec angular resolutions, several neighbouring line-of-sights pass the same boxes, which results in extended coherent features in the simulated RM maps.
For the northern patches HN and MN, we also performed simulations with Gaussian random magnetic field component for comparison. Both RM maps are shown in the first two rows in Fig. 3 and prove the above expectations. As seen from Fig. 4, the PDFs of the RM simulations with Gaussian random magnetic fields conform to Gaussian distributions nicely, while the PDFs from the simulations with Kolmogorov-like random magnetic fields do not. Quantitatively, this can be verified by the skewness and kurtosis in Table 2, where the values of these two parameters are much higher for Kolmogorov-like magnetic field spectra than Gaussian ones. The structure functions shown in the bottom panel of Fig. 4, as well as the autocorrelation functions in Fig. 5, reinforce the argument. The structure functions for Gaussian random magnetic fields are flat, while the structure functions for a Kolmogorov-like magnetic field spectrum have slopes, which is discussed below.
![]() |
Figure 5: Autocorrelation functions for RMs. |
Open with DEXTER |
4.4 RM patches from different realizations
The RM structures in the present simulations are no predictions for sky emission structures towards the quoted directions, but are just one possible realization of the random magnetic fields following a Kolmogorov-like power-law spectrum in context with other model parameters as described before.
We generated the random magnetic field components in the following way. Their amplitudes are determined according to the slope of the power-law spectrum with a spectral index of -11/3, and the phase angles are realised by random numbers. The rotation angle of the numerous 10 pc boxes placed along the line-of-sight are also determined from random numbers. These random numbers are actually pseudo random and are initiated by an integer called ``seed''. When the random ``seed'' number is changed, a different realization of the same random magnetic field distribution is obtained.
To check structure variations and structure functions for different
realizations, we ran two more simulations for patch MS with different ``seed''
numbers. The RM maps are displayed in the fifth row of Fig. 3 and
show different distributions of RM features as expected. Also the PDFs for the
RM maps (Fig. 6) vary considerably, as do the parameters
listed in Table 2. The structure functions of RMs, however,
are almost identical as shown in the bottom panel of Fig. 6,
except for regions with large angular separations exceeding about
,
where the structure function becomes less accurate due to the decreasing number
of selected pixels.
The mean RM for the three realizations of patch MS is between -20 rad m-2 and -28 rad m-2 (Table 2), and the majority of RMs have negative values. However, because of the variance of up to 19 rad m-2 (Table 2), all three simulations show a number of small extended features in the maps with positive sign. Naively this could be interpreted as reversals of the large-scale magnetic field in case the available number of observed RMs from EGSs is small. However, in contrast to large-scale reversals in the Galactic plane as discussed for instance by Han et al. (2006), the reversals seen here are caused entirely by the distribution of random magnetic fields. For this medium-latitude area, we roughly estimate the probability for picking up a ``reversal'' to be about 7%. This should be reflected in the distribution of observed RMs of EGSs and will prove the relevance of the simulations and the parameters used.
![]() |
Figure 6: PDFs and structure functions of the patch MS for the same amplitudes but different realizations of the phases of the random magnetic field components created in a 10 pc box (see text for details). |
Open with DEXTER |
4.5 Galactic-plane RM patches at different longitudes
In this section we investigate RM variations in the plane versus Galactic
longitudes. The four simulated patches have latitudes of about
.
Two
of them (PNC1 and PNC2) are in the inner Galaxy direction and the other two
(PNAC1 and PNAC2) are in the anti-centre direction (see Table 1).
The visible structures in the maps (the third row in Fig. 3) are less extended than those at higher latitudes. There is more fine structure for the inner Galaxy regions than for those in the anti-centre direction. The PDFs of all four simulated patches closely resemble Gaussian shapes (Fig. 7). The inner Galaxy patches PNC1 and PNC2 have lower skewness values than the anti-centre patches (Table 2), indicating that their PDFs are closer to Gaussian shapes. The skewness and kurtosis for all patches are very small, which indicates that it is impossible to infer whether the turbulence spectrum is Gaussian or Kolmogorov-like based on their PDF. The structure functions as discussed before, however, clearly reveal the properties of the underlying turbulence spectrum. The average RM values are consistent with the RM results obtained by Sun et al. (2008).
As displayed in Fig. 7, all Galactic plane structure functions are
well-fitted by power laws. Compared to the medium- and high-latitude
fields, a flattening of the structure functions occurs at smaller angular
distances. The fit results for the power laws are listed in
Table 3. From the parameters listed in Table 1 we obtain
a transition angle of about 1
for PNC1 and PNC2 and about
for PNAC1 and PNAC2. Figure 7 shows a smooth transition from a
linear slope to a flat curve in the logarithmic scale instead of a sharp change in the
structure function. The turnover of the structure functions ranges from
about
to
for patches PNC1 and PNC2, from
to
for patch PNAC1, and from
to
for patch PNAC2,
which are roughly consistent with the aforementioned expectations. The spectral
indices are consistent with theoretical predictions by Minter & Spangler (1996)
for 2D turbulence, which results in a spectral index of about 0.7
for the structure functions. However, our input turbulent magnetic
fields are for 3D. The amplitudes of the structure functions increase
towards the inner Galaxy regions, which is expected because the
line-of-sights pass more turbulent cells or magnetic field boxes than
in the anti-centre or high-latitude direction. The amplitude of the
structure function for PNC1 is larger than that for PNC2, which is
mainly caused by a higher electron density in this direction according
to the NE2001 model.
![]() |
Figure 7: PDFs and structure functions of the five simulated Galactic plane patches (see Table 1 for field parameters). |
Open with DEXTER |
4.6 RM patches at different latitudes
We finally compared the different results towards different latitudes. The
patches HN/HS, MN/MS,and PNAC1/PSAC have longitudes between
and
.
The latitudes are
,
and
.
The ``seed1'' case was used for patch MS. The six simulated RM maps are all shown
in the first column of Fig. 3. A clear trend toward more extended
structures is clearly seen for high-latitude regions. The PDFs for the
high-latitude regions (HN/HS and MN/MS) all deviate significantly from Gaussian (Fig. 8 and Table 2).
The structure functions derived for the four high-latitude patches are shown in
Fig. 8, as well as for the five plane patches in Fig. 7. All
of them have similar properties. The transition angle is calculated to be
about
for HN/HS and MN/MS regions, which is consistent with the
results displayed in Fig. 8. The amplitudes of the structure
functions are smaller for the high-latitude patches compared to those in the
plane. This can be understood because the length along the line-of-sight is
short compared to that in the Galactic plane. The number of turbulent cells
passed by the line-of-sights is smaller as is the amplitude of the
structure functions. We also note that the slopes of the structure functions
increase for the high-latitude patches. These spectral changes are
discussed below.
![]() |
Figure 8: PDFs and structure functions of four patches (see first column of Fig. 3) at roughly the same Galactic longitude located at medium and high latitudes (see Table 1). |
Open with DEXTER |
4.7 Discussion of RM simulations
4.7.1 Theoretical structure functions revisited
From the structure functions of simulated RM maps, we found evidence of latitude-dependent variations in their spectral slope, which gets steeper towards high latitudes. The regular magnetic field and the electron density have spatial irregularities, which have a larger scale than the sizes of the patches. Therefore the large-scale magnetic field cannot cause variations in the spectral slopes, but it does influence the amplitude of the structure functions. The coupling between the electron density and random magnetic fields introduced by Sun et al. (2008) augments the scattering of RMs, hence the amplitude of the structure function rather than changes its spectral index. For the present simulations, the random magnetic field is comparable to the regular magnetic field strength, and it dominates the regular field towards high-latitude regions. We conclude that all RM structures across the simulated patches are caused by the random magnetic field properties.
Table 4: Fit parameters for structure functions from different simulation settings.
Table 5: Statistical parameters for the simulated I, PI, and PA patches. The patch designations are listed in Table 1.
The variations versus latitude must be ascribed to geometric effects, which
reflects in the integral length (
)
as listed in Table 1.
The size of the magnetic field box (L)
is fixed in the simulations presented; therefore, the integral length
determines the number of turbulent cells or boxes passed by the
line-of-sight. When there are a large enough number of cells, the
fluctuations along the line-of-sight are averaged out. This results in
a spectral slope of
for the structure function.
This spectral index was observed and interpreted as a 2D turbulence in the
interstellar medium by Minter & Spangler (1996). The 2D turbulence is confined to thin sheets or filaments. However, the present simulations are based on a 3D distribution
of the turbulent magnetic field components, which do not have the shape of sheets
or filaments. When there are not enough cells for the integral
length to be comparable to the outer scale, the structure function tends to have
the index of
as derived by (Minter & Spangler 1996) for 3D
Kolmogorov-like turbulence. Therefore we conclude that, in the inertial range
(
), the derived
structure functions behave like
with 2/3<m<5/3
depending on the number of cells. When the angular separation is
greater than the transition angle, the structure function spectrum
flattens and m gradually turns to zero.
To prove these arguments, we made several simulations for a patch of about
in the plane with
corresponding to
a resolution of about
.
Although the resolution is much lower than
the high-resolution simulations for the patches in Table 1, the
dependence of the structure functions should be the same. These simulations can
be realised in much less computing time than the arcsec resolution
simulations. We fit the structure functions by using
and list the results
including errors of C and m in Table 4. We first fix the outer
scale (L) and the integral length and vary the size of the small cubes
(
). The results in Table 4
show that the size of the small cubes has little influence on the
structure function parameter. We then change the outer scale, but keep
the size of the small cubes and the integral length constant. As
expected, the spectral index of the structure functions tends to be
about 0.7 for very small outer scales, which means there are many
turbulent cells. Finally, we increase the integral length from
10 pc to 10 kpc by using the same values of cube size and
outer scale. The spectral index of the structure functions varies
between 1.3 and 1.7 until the integral length exceeds about
1600 pc, becoming smaller for longer lengths. These test
simulations support our conclusion that the spectral index of the
structure function primarily depends on the number of turbulent cells,
which is mainly determined by the length of the line-of-sight. Thus,
when assessing the turbulent magnetic field properties from the
structure functions of RMs, the line-of-sight dependence must be taken
into account. The integral length is always large in the Galactic
plane, which assures that enough turbulent cells are passed through by
the line-of-sight; therefore, the relation
holds in general. Towards higher Galactic latitudes, however, the length of the
line-of-sight decreases, so the number of turbulent cells also decreases.
The spectral index of the random fields approaches
.
![]() |
Figure 9: Structure functions for the total intensity I ( top panels) and the polarized intensity PI maps ( bottom panels). |
Open with DEXTER |
4.7.2 Diagnosis of the galactic foreground
It has been proposed to conduct a very dense RM survey of EGS with the SKA to
study magnetic field structures in the Milky Way, in galaxies, in clusters of
galaxies and beyond (Beck & Gaensler 2004). The key issue to achieve relevant information
for extragalactic objects is to properly extract and separate the foreground
RM contributions from the observations. We tested the minimum number
of pixel needed to calculate the correct structure function for the simulated
maps and found that about 10 000 pixel are required. For our map size, this
corresponds to a mean angular separation of about
in case the EGS
are not larger than the pixel size of
and their intrinsic RMs
average out.
In the plane the PDFs are close to Gaussian. Therefore the average of RMs is a good estimate of the foreground contribution caused by the large-scale magnetic field and the thermal electron distribution as modelled by Sun et al. (2008). However, large variances (Table 2) and amplitudes of the structure functions (Table 3) imply that, after having subtracted the average RM in a certain direction, the RM scattering of the foreground emission is still preserved. Thus patches in the plane are certainly not well suited to studying extragalactic magnetic fields.
At high latitudes the PDFs deviate significantly from a Gaussian distribution
and also more extended structures show up (Fig. 3). To take the
foreground RMs appropriately into account a dense sample of EGS is required.
In case an extended extragalactic object is seen towards an excessive RM patch
of the same size, additional information is required to solve the RM foreground
problem. On the other hand, the variation in the RMs are much less than in
the plane, thus high-latitude regions seem better suited to studying extragalactic
magnetic fields. However, according to Kronberg et al. (2008), cosmological magnetic
fields of a few 10-7 G at a redshift of 2 have an expected RM excess
of up to 20 rad m-2. This is within the same order of magnitude as the RM
variance of the Galaxy at high latitudes obtained from the present simulations,
which means a detection of cosmological magnetic fields is quite challenging.
Reducing the coupling factor between electron density and random
fields introduced by Sun et al. (2008) to explain the observed depolarization of
Galactic synchrotron emission at low angular resolution will in turn also
reduce the foreground RM variance. Therefore future high-resolution
observations of the small-scale structures of the magnetized ISM are needed to
investigate the coupling factor in detail. Such observations are planned for
LOFAR, which will operate at metre-wavelengths and will thus be
able to detect very small RM-variations from small-scale Galactic emission
features.
5 Simulated total intensity and polarization maps
We briefly describe the properties of the simulated total intensity and linear polarization maps emerging from the same simulation runs where the already discussed RM distributions were calculated. These simulated Stokes I, U, Q maps are all for the frequency of 1.4 GHz, whereas the RM distribution is frequency-independent.
The simulated total intensity I, U, Q, polarized intensity PI, and polarization angle PA maps are shown in Figs. B.1-B.4, and B.5 for the same patches and simulation variations as displayed in Fig. 3 for RMs. The statistical parameters for the I, PI, and PA maps are calculated and listed in Table 5. The structure functions for I and PI are shown in Fig. 9.
Similar to RM maps, the patches in the plane are dominated by structures on very small scales, while the patches at high latitudes exhibit more extended features. These structures are all related to the properties of the turbulent magnetic fields. The variation in the sizes of the features primarily stems from the much larger integral length along the line-of-sights in the plane than at high latitudes.
The large-scale cosmic-ray and magnetic field also have imprints on the simulations. For example, the clear gradient of the total intensity, particularly in the plane maps towards the inner Galaxy, traces the large-scale magnetic field and the cosmic-ray distribution. The asymmetric halo field adds up to the disk field and enhances the field strength below the Galactic plane, which produces higher intensities there. This can be seen from the averages of I and PI intensities listed in Table 5 for the PSAC and the PNAC1 patch.
The average polarization percentage is about 10% for the HN and MN patches and about 16% for the HS and MS patches. Towards these directions, both RM and its variance are small so that Faraday depolarization is low as well (Sokoloff et al. 1998). The low degree of polarization thus comes from the random magnetic fields, which dominate the field components. The percentage becomes about 7% for patch PNAC1, 5% for patch PNAC2, and 2% for the inner Galaxy patches, which is caused by strong Faraday depolarization, as both the RM and the RM scattering are large. There is no indication of a correlation among the variances of RM, polarized intensity PI, and polarization angle PA.
The structure functions for total intensity I (Fig. 9) have a linear slope on the logarithmic scale, but these are generally more shallow than those for the corresponding RMs. The total intensity I and polarized intensity PI maps at high latitudes have quite similar structure functions, while they are fairly distinct in the plane. The structure functions of PI are nearly flat. This indicates that strong depolarization in the plane means the PI structure functions no longer contain the fluctuation information of the turbulent magnetic field.
6 Conclusions
Galactic 3D-models of the distribution of thermal electrons, cosmic-ray
electrons, and magnetic fields were used to calculate total intensity,
polarization, and RM maps for selected patches of the Galaxy via the
HAMMURABI code. The angular resolution of the maps is about
,
close to that of the planned SKA. The simulated patches were
distributed at different longitudes and latitudes, which enabled us to study
spatial variations of the results. The simulated RM maps are of particular
interest, as they are directly related to the planned deep RM EGS survey with
the SKA. Complete SKA simulations in addition need simulated sky patches
with distributed polarized sources. Similar simulations were made by
Wilman et al. (2008).
In the present simulations the random magnetic field component is assumed to follow a Kolmogorov-like power spectrum. Compared to purely random Gaussian magnetic fields, structures on different scales show up in the RM, as well as in the total intensity I and the polarized intensity PI maps. The simulations predict more extended structures for high-latitude regions.
We study the PDFs and structure functions for the RM maps. The PDFs of the patches in the plane are close to Gaussian, whereas at high-latitudes clear deviations from Gaussian shapes are noted. The structure functions have a power-law slope for angular separations smaller than the transition angle and flatter for larger angular separations. The transition angle is determined by the outer scale of the turbulence and by the length of the line-of-sight. The amplitudes of the structure functions become smaller and the slopes of the structure functions steepen towards higher-latitude regions. All this is entirely caused by geometric effects. In the plane, the integral length exceeds that at high latitudes, therefore more turbulent cells are passed by the lines-of-sight. This significantly smears out fluctuations along the lines-of-sight, which makes the 3D turbulence observationally a 2D turbulence. Thus the structure functions are flatter in the plane than at high latitudes. We notice that the PDFs of RMs have a wide distribution in all directions, which potentially leads to depolarization of Galactic radio emission, when observed with larger beams. It also strongly influences the interpretation of discrete features sampled by RMs of EGS. It is obvious that only a very dense grid of extragalactic RMs will help to separate Galactic foreground influence from extragalactic contributions. This is planned for the SKA.
The simulations at 1.4 GHz presented in this paper can be easily run at any other frequency. While the RM maps are frequency independent, total intensity and polarization simulations are of particular interest for low-frequency simulations, where LOFAR is expected to do sensitive high-resolution observations in the near future. At low frequencies, thermal absorption adds to depolarization effects and must be properly taken into account. LOFAR will make a major attempt to trace the very weak signal from the Epoch of reionization, where the selection of a high-latitude Galactic area with very low foreground signal is important for the success of these ambitious low-frequency observations as well as the subtraction of all kinds of extragalactic foregrounds as discussed by Jelic et al. (2008).
AcknowledgementsX. H. Sun acknowledges support by the European Community Framework Programme 6, Square Kilometre Design Study (SKADS), contract no. 011938. We are very grateful to Andrè Waelkens and Torsten Enßlin for providing the HAMMURABI code and for many detailed discussions on its application, improvements and results. We thank Michael Kramer and Patricia Reich for critical reading of the manuscript and helpful comments. We like to thank Walter Alef for support in using the MPIfR PC-cluster operated by the VLBI-group, where all HAMMURABI based simulations were performed.
Appendix A: Details for the realization of random magnetic field components
The general concept is to generate random magnetic fields in a box by a
fast Fourier transformation (FFT), randomly rotate the box around its
centre to obtain a different view on the box, and finally fill the
volume of the Galaxy with these randomly rotated boxes. The size of a
box is L3 and the size of the Galaxy is
.
We first define several Cartesian coordinate systems:
- (X, Y, Z): with XY-plane parallel to the Galactic plane and the Galactic centre located at (DX/2, DY/2, DZ/2).
- (x, y, z): based on the box with the origin at its vertex. When the Galaxy is filled with boxes the origin of the first box is the same as for the system (X, Y, Z) and the xy-plane is parallel to the XY-plane.
-
or
: the same as (x, y, z), but with the origin at the centre of the box.
We first consider the random magnetic field distribution in a box. The random
magnetic field is assumed to be isotropic and homogeneous, and its three
components do not correlate with each other. The three components thus follow
the same power law spectrum
.
Here Cb2 is the
amplitude,
is the wave vector with l being the spatial scale and
is the spectral index. For the simulations presented in this paper a Kolmogorov-like turbulence is assumed, which means
.
We first generate the 1D component of the random magnetic field. The other two
components were obtained by simply repeating the process for the first one.
This 1D component denoted as
b(x, y, z) is the inverse Fourier
transformation of the power spectrum as the following,
where






where



![]() |
(A.3) |
with the same for km and kn. Generally klmn is limited to the inertial range





-
(A.4)
-
is a real number, if the following is satisfied,
(A.5)
![[*]](/icons/foot_motif.png)
![]() |
Figure A.1:
2D-scheme for the realization of random magnetic fields. The sizes of the box and the cube are L and |
Open with DEXTER |
To calculate the random magnetic field in any position (X, Y, Z) in the Galaxy b(X, Y, Z), the following procedure is executed (Fig. A.1),
- 1.
- transform (X, Y, Z) to
as,
(A.6)
where,
and
are the integer part.
- 2.
- rotate
by Euler angles
around the centre of the box to get
as,
(A.7)
and(A.8)
(A.9)
(A.10)
Here the Euler angles are obtained as,(A.11)
(A.12)
and(A.13)
Here we use three random numbers,
, and
to determine the Euler angles. For simplicity the Euler angles are set to be integer multiple of
to avoid complex projections.
- 3.
- transform
to (x, y, z) then to
,
(A.14)
and(A.15)
Finally.
Appendix B: Simulated total intensity and polarization maps
We display the total intensity I maps in Fig. B.1, Stokes U maps in Fig. B.2, Stokes Q maps in Fig. B.3, polarized intensity PI maps in Fig. B.4 and polarization angle PA maps in Fig. B.5.
![]() |
Figure B.1: Total intensity I maps for all the simulated patches in gnomonic projection. The centre coordinates are given on top of each panel together with its label (see Table 1 for details). Intensities are in mK brightness temperature at 1.4 GHz. The minimum and maximum of each map is shown together with the wedge below each panel. |
Open with DEXTER |
![]() |
Figure B.2: As Fig. B.1, but for the simulated Stokes U maps. |
Open with DEXTER |
![]() |
Figure B.3: As Fig. B.1, but for the simulated Stokes Q maps. |
Open with DEXTER |
![]() |
Figure B.4: As Fig. B.1, but for PI maps calculated from the simulated U and Q maps. |
Open with DEXTER |
![]() |
Figure B.5: As Fig. B.1, but for the polarization angles PA calculated from the simulated U and Q maps. |
Open with DEXTER |
References
- Armstrong, J. W., Rickett, B. J., & Spangler, S. R. 1995, ApJ, 443, 209 [CrossRef] [NASA ADS]
- Beck, R., & Gaensler, B. M. 2004, New Astron. Rev., 48, 1289 [CrossRef] [NASA ADS]
- Berkhuijsen, E. M., Mitra, D., & Müller, P. 2006, AN, 327, 82 [NASA ADS]
- Brown, J. C., Taylor, A. R., & Jackel, B. J. 2003, ApJS, 145, 213 [CrossRef] [NASA ADS]
- Brown, J. C., Haverkorn, M., Gaensler, B. M., et al. 2007, ApJ, 663, 258 [CrossRef] [NASA ADS]
- Cordes, J. M., & Lazio, T. J. W. 2002, preprint [arXiv:astro-ph/0207156]
- Enßlin, T. A., & Vogt, C. 2003, A&A, 401, 835 [EDP Sciences] [CrossRef] [NASA ADS]
- Gaensler, B. M., Madsen, G. J., Chatterjee, S., & Mao, S. A. 2008, PASA, 25, 184 [NASA ADS]
- Górski, K. M., Hivon, E., Banday, A. J., et al. 2005, ApJ, 622, 759 [CrossRef] [NASA ADS]
- Han, J. L., Manchester, R. N., & Qiao, G. J. 1999, MNRAS, 306, 371 [CrossRef] [NASA ADS]
- Han, J. L., Manchester, R. N., Lyne, A. G., Qiao, G. J., & van Straten, W. 2006, ApJ, 642, 868 [CrossRef] [NASA ADS]
- Haslam, C. G. T., Salter, C. J., Stoffel, H., & Wilson, W. E. 1982, A&AS, 47, 1 [NASA ADS]
- Haverkorn, M., Brown, J. C., Gaensler, B. M., & McClure-Griffiths, N. M. 2008, ApJ, 680, 362 [CrossRef] [NASA ADS]
- Hinshaw, G., Nolta, M. R., Bennett, C. L., et al. 2007, ApJS, 170, 288 [CrossRef] [NASA ADS]
- Jelic, V., Zaroubi, S., Labropoulos, P., et al. 2008, MNRAS, 389, 1319 [CrossRef] [NASA ADS]
- Kowal, G., Lazarian, A., & Beresnyak, A. 2007, ApJ, 658, 423 [CrossRef] [NASA ADS]
- Kronberg, P. P., & Perry, J. J. 1982, ApJ, 263, 518 [CrossRef] [NASA ADS]
- Kronberg, P. P., Bernet, M. L., Miniati, F., et al. 2008, ApJ, 676, 70 [CrossRef] [NASA ADS]
- Minter, A. H., & Spangler, S. R. 1996, ApJ, 458, 194 [CrossRef] [NASA ADS]
- Mitra, D., Wielebinski, R., Kramer, M., & Jessner, A. 2003, A&A, 398, 993 [EDP Sciences] [CrossRef] [NASA ADS]
- Mizeva, I., Reich, W., Frick, P., Beck, R., & Sokoloff, D. 2007, AN, 328, 80 [NASA ADS]
- Moss, D., & Sokoloff, D. 2008, A&A, 487, 197 [EDP Sciences] [CrossRef] [NASA ADS]
- Noutsos, A., Johnston, S., Kramer, M., & Karastergiou, A. 2008, MNRAS, 386, 1881 [CrossRef] [NASA ADS]
- Page, L., Hinshaw, G., Komatsu, E., et al. 2007, ApJS, 170, 335 [CrossRef] [NASA ADS]
- Simonetti, J. H., Cordes, J. M., & Spangler, S. R. 1984, ApJ, 284, 126 [CrossRef] [NASA ADS]
- Sokoloff, D. D., Bykov, A. A., Shukurov, A., et al. 1998, MNRAS, 299, 189 [CrossRef] [NASA ADS]
- Sun, X. H., & Han, J. L. 2004, in The Magnetized Interstellar Medium, ed. B. Uyaniker, W. Reich, & R. Wielebinski (Copernicus GmbH), 25
- Sun, X. H., Reich, W., Waelkens, A., & Enßlin, T. A. 2008, A&A, 477, 573 [EDP Sciences] [CrossRef] [NASA ADS]
- Taylor, A. R., Stil, J. M., Grant, J. K., et al. 2007, ApJ, 666, 201 [CrossRef] [NASA ADS]
- Testori, J. C., Reich, P., & Reich, W. 2008, A&A, 484, 733 [EDP Sciences] [CrossRef] [NASA ADS]
- Waelkens, A. 2005, Diploma Thesis, Ludwig-Maximilian-Universität München
- Waelkens, A., Jaffe, T., Reinecke, M., Kitaura, F. S., & Enßlin, T. A. 2009, A&A, 495, 697 [EDP Sciences] [CrossRef] [NASA ADS]
- Wilman, R. J., Miller, L., Jarvis, M. J., et al. 2008, MNRAS, 388, 1335 [NASA ADS]
- Wolleben, M., Landecker, T. L., Reich, W., & Wielebinski, R. 2006, A&A, 448, 441 [CrossRef] [NASA ADS]
Footnotes
- ... numbers
- Throughout all the simulations, we use the random number generator from the GNU Scientific Library. The algorithm is MT19937, which has a Mersenne prime period of 219 937-1.
- ... library
- www.fftw.org
All Tables
Table 1: Simulated patches.
Table 2: Statistical parameters for the simulated RM patches.
Table 3: Fitting results for structure functions and autocorrelation functions of RMs.
Table 4: Fit parameters for structure functions from different simulation settings.
Table 5: Statistical parameters for the simulated I, PI, and PA patches. The patch designations are listed in Table 1.
All Figures
![]() |
Figure 1: The HEALPix pixel scheme (Górski et al. 2005) for Galactic coordinates used in this paper. |
Open with DEXTER | |
In the text |
![]() |
Figure 2: 2D random magnetic field distribution in a slice parallel to one surface of a 10 pc box. The upper left panel shows magnetic field components following a power law index of -11/3, while the upper right panel shows random field components realised by random Gaussian numbers corresponding to a spectral index of zero. Their corresponding power spectra are shown in the lower panels. |
Open with DEXTER | |
In the text |
![]() |
Figure 3:
RM maps for all simulated patches ordered from high to low latitudes ( top to bottom)
in gnonomic projection. Their designation and Galactic coordinates are
quoted on top of each individual panel (also listed in Table 1).
Different simulations for the same patch and for patches at the same
latitude, but different longitudes, are shown in one row (see text for
details). Each field has about
|
Open with DEXTER | |
In the text |
![]() |
Figure 4: PDFs and structure functions of the patches HN and MN at high and medium latitudes. See Tables 1 and 2 for details of the patch parameters. |
Open with DEXTER | |
In the text |
![]() |
Figure 5: Autocorrelation functions for RMs. |
Open with DEXTER | |
In the text |
![]() |
Figure 6: PDFs and structure functions of the patch MS for the same amplitudes but different realizations of the phases of the random magnetic field components created in a 10 pc box (see text for details). |
Open with DEXTER | |
In the text |
![]() |
Figure 7: PDFs and structure functions of the five simulated Galactic plane patches (see Table 1 for field parameters). |
Open with DEXTER | |
In the text |
![]() |
Figure 8: PDFs and structure functions of four patches (see first column of Fig. 3) at roughly the same Galactic longitude located at medium and high latitudes (see Table 1). |
Open with DEXTER | |
In the text |
![]() |
Figure 9: Structure functions for the total intensity I ( top panels) and the polarized intensity PI maps ( bottom panels). |
Open with DEXTER | |
In the text |
![]() |
Figure A.1:
2D-scheme for the realization of random magnetic fields. The sizes of the box and the cube are L and |
Open with DEXTER | |
In the text |
![]() |
Figure B.1: Total intensity I maps for all the simulated patches in gnomonic projection. The centre coordinates are given on top of each panel together with its label (see Table 1 for details). Intensities are in mK brightness temperature at 1.4 GHz. The minimum and maximum of each map is shown together with the wedge below each panel. |
Open with DEXTER | |
In the text |
![]() |
Figure B.2: As Fig. B.1, but for the simulated Stokes U maps. |
Open with DEXTER | |
In the text |
![]() |
Figure B.3: As Fig. B.1, but for the simulated Stokes Q maps. |
Open with DEXTER | |
In the text |
![]() |
Figure B.4: As Fig. B.1, but for PI maps calculated from the simulated U and Q maps. |
Open with DEXTER | |
In the text |
![]() |
Figure B.5: As Fig. B.1, but for the polarization angles PA calculated from the simulated U and Q maps. |
Open with DEXTER | |
In the text |
Copyright ESO 2009
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.