A&A 467, 411-419 (2007)
DOI: 10.1051/0004-6361:20065086
M. Maturi^{1,2,3,4} - T. Enßlin^{4} - C. Hernández-Monteagudo^{4,5} - J. A. Rubiño-Martín^{4,6}
1 -
Dipartimento di Astronomia, Università di Padova, Vicolo
dell'Osservatorio 2, 35120 Padova, Italy
2 -
ITA, Universität Heidelberg,
Albert-Überle-Str. 2, 69120 Heidelberg, Germany
3 -
Dipartimento di Astronomia, Università di Bologna, via
Ranzani 1, 40127 Bologna, Italy
4 -
Max-Planck-Institut für Astrophysik, PO Box 1523,
85740 Garching, Germany
5 -
Dept. of Physics & Astronomy, University of Pennsylvania,
209 South 33rd Str. Philadelphia, PA 19104-6396, USA
6 -
Instituto de Astrofísica de Canarias, C/Vía Láctea s/n,
38200 Tenerife, Spain
Received 24 February 2006 / Accepted 8 January 2007
Abstract
We discuss the prospects for detecting the Rees-Sciama
(RS) effect associated with merging clusters of galaxies in millimetric
observations. Using the extended Press-Schechter formalism, we characterize
the distribution of the parameters describing the mergers, and use these
results to generate a mock catalogue of cluster mergers. We use realistic
N-body simulations to confirm the reliability of the mock catalogue, and
build a new filter that attempts to optimize the extraction of the RS
signal. This is an extension of the spatial-filtering method described in
Haehnelt & Tegmark (1996, MNRAS, 279, 545), that has a general definition and hence can be
applied in other contexts, such as gravitational lensing of the CMB or
lensing of background galaxies. When applied to our mock cluster merger
catalogue, we show that with the announced sensitivities of future
experiments like the Atacama Cosmology Telescope (ACT), the South Pole
Telescope (SPT), or the Atacama Large Millimeter Array (ALMA), a detection
of the signal will be possible if we consider roughly 1000 cluster
mergers.
Key words: cosmology: theory - galaxies: clusters: general - cosmic microwave background
Observations of the fluctuations of the cosmic microwave background (CMB) can give us information about the formation of structures in the Universe, because the evolution of the gravitational potentials leaves its imprint on the CMB photons as they travel along them. This physical effect is usually split into two terms. One is the integrated Sachs-Wolfe effect (ISW; Hu & Sugiyama 1994; Sachs & Wolfe 1967), which is produced by the linear evolution of the potentials. The other is the Rees-Sciama effect (RS; Seljak 1996; Rees & Sciama 1968; Martinez-Gonzalez et al. 1990), in which the density contrast producing the gravitational potential is in its non-linear regime.
In this work we discuss the RS effect associated with the non-linear regime of galaxy cluster mergers. In a previous work (Rubiño-Martín et al. 2004, hereafter Paper I), this problem was examined using a simplified model of the physics of the merger. The amplitude, morphology, and symmetries of the effect were discussed and characterized in terms of the physical parameters of the merger. However, the RS effect is not the only signal being generated by clusters. The presence of hot gas in the intra-cluster medium (ICM) and the peculiar velocity of the cluster with respect to the Hubble flow introduce temperature fluctuations that are proportional to the integral of pressure (thermal Sunyaev-Zel'dovich effect, [tSZ]) and the momentum density (kinematic Sunyaev-Zel'dovich effect, [kSZ]) along the line of sight (Sunyaev & Zel'dovich 1980; Sunyaev & Zeldovich 1972). In addition, galaxy clusters may host radio-galaxies or infrared galaxies whose flux might be relevant to the small-amplitude effects we are studying here. Furthermore, the intrinsic CMB field and its deflection caused by the matter distribution in the cluster produce different anisotropies that must be accounted for. In Paper I, it was proposed to coherently co-add the RS signals from a sample of cluster mergers in order to statistically detect the effect against these contaminating anisotropy signals.
Here, we have extend this work in three ways. First, we used hydrodynamical simulations in order to compute realistic maps of the RS effect and to show that the morphology of the signal has the characteristic quadrupole structure described in Paper I. We also showed that N-body simulations contain an additional signal that originates during the collapse of the environmental dark-matter in the potential well of the system, although the morphology of this signal is completely different from the one associated with the one induced by merging. Second, we performed a detailed study of the expected distribution of the parameters describing a merger of galaxy clusters (masses of the components, distances and velocities) and presented semi-analytical expressions for the probability of finding a cluster merger with a given parameter set. These equations can be used to produce a realistic "catalogue'' of cluster mergers, which we used in order to forecast the detectability of the effect. Finally, we investigated the prospects for the RS signal extraction from cluster mergers by future CMB experiments, like the South Pole Telescope^{} or the Atacama Cosmology Telescope^{}. To this end, we extended the filtering method described in Haehnelt & Tegmark (1996). Our new filter has a general mathematical definition, so it could be applied to extracting the signal in many other fields, e.g. gravitational lensing of the CMB (Seljak 1996; Maturi et al. 2005a) or lensing of background galaxies (e.g. Maturi et al. 2005b). We used this extended filter to separate the RS signal from the following three components: the primordial CMB, the gravitational lensing imprint, and the kinetic SZ effect. In doing that we analysed a simple and optimistic scenario, showing that with upcoming instruments the RS signal may be detected if we consider roughly 1000 observed cluster mergers. As we will see below, the main limiting "noise'' component in such a measurement will be the lensed CMB.
We underline that a number of simplifying and optimistic assumptions were made in this work. It could well be that observational imperfections drive our suggested sigma detection level in the ideal case into only a weakly or even non-significant signal in a more realistic case. In this sense, this work presents the difficulties found in estimating the RS measurements and a possible way to overcome them.
In Paper I, the Rees-Sciama effect associated with mergers of galaxy clusters was discussed using an analytical model for the merger. In addition, an analytical prescription to prepare Rees-Sciama maps from hydrodynamical simulations was introduced. Here, we implement this prescription and present an RS map resulting from a numerical simulation. The morphology of the effect can then be directly compared to what is predicted in Paper I.
In order to prepare a map of the RS effect from a hydrodynamic simulation, we
used the analytic recipe proposed in Paper I to compute the CMB
temperature change along the line-of-sight (LOS) direction ,
Due to its gravitational nature, the RS effect can also be related to the
deflection angle, yielding a very useful expression for analytical
computations:
The dot product has a dipolar structure elongated along and, according to Eq. (2), the RS anisotropy is proportional to the halo mass and to its transverse velocity. Due to these dependencies, forming galaxy clusters with two major merging haloes seem to be the best candidates for RS detection; first, because of the large masses and velocities involved, on the order of and respectively; second, because the contributions of the merging components add constructively due to the converging cluster motions.
A simple description of this scenario can be built by considering two dark matter haloes, having a NFW profile (Navarro et al. 1995), in free fall. It provides a simple derivation for the lensing and the RS effects, from which it is straightforward to obtain the templates necessary for the data reduction (Sect. 4). The RS effect can be computed according to Eq. (2), where the deflection angle can be computed analytically (Bartelmann 1996). As will be shown, the fact that the deflection field is axially symmetric on large scales justifies the assumption of the haloes being spherical.
As an example, we present the RS effect of a major merger event in Fig. 1. In the left panel, the RS signal was computed with our simple analytical recipe, whereas in the right panel we used the numerically simulated cluster g24+200 described in Rasia et al. (2004). This cluster consists of two merging haloes, of masses and , with an infall velocity of , provided by G. Tormen. The field of view is 2 degrees on a side. This figure demonstrates that our simple analytical model can describe the RS signal generated in major cluster mergers with reasonable accuracy. The main differences between our description and the output obtained from the numerical simulations are related to the presence of minor merging substructures within the main object and to the roughly spherical collapse of environmental dark matter (DM) onto the potential. This collapse-induced signal has the shape of a negative peak with circular symmetry centered on the potential minimum. This contribution adds coherently to the signal of the two main halo mergers and is strictly related to the merging dynamic we want to estimate. Since we are ignoring this term, the resulting signal amplitude is smaller, making this choice conservative. Also this extra RS component is related to the dynamical status of DM and could be considered part of the signal we are searching for. The inclusion of this contribution in our procedure can be easily implemented by simply changing the filter's RS template (see Sect. 4.2). In this case we would measure the full DM momentum instead of only the DM momentum of the main mergers. Since we ignore this signal, our estimate of the detectability of the RS signal is conservative. The effect of minor mergers can be neglected, since their amplitude is small and their contribution vanishes after applying our filter (see Sect. 4.2 below).
Figure 1: A comparison between our analytical description of the RS signal generated in a major merger ( left panel) and a numerical simulation ( right panel). The analytical model represents the same merger of the numerical simulation. | |
Open with DEXTER |
We investigated the probability distribution of the physical parameters
describing a merger. We used the approach of Lacey & Cole (1993) based on the
extended Press-Schechter formalism to model the distribution of cluster
mergers in the Universe. This description basically relies on the probability
of a region of mass M_{1} having over-density
provided that,
at the same point, the over-density measured on a larger scale (corresponding
to a mass M_{2}) is .
In other words, our starting point is the
conditional probability of having
on scales of M_{1}, given
that on scales corresponding to M_{2} the density contrast is
We follow the model given in Sarazin (2002) when assigning distances and
velocities to the merging clusters. The typical initial distance between a
pair of clusters, following Eq. (10) in Sarazin (2002), may be approximated as
The millimetric observations of upcoming experiments will provide extensive catalogues of galaxy clusters detected through their thermal SZ effect (tSZ), together with estimates of their masses and their geometrical properties, such as the merger components' centers, and hence the axes joining them. If the main merger components cannot be resolved, the tSZ detection would show a strong elongation anyway, which can be uniquely related to the un-relaxed clusters we are looking for. On the other hand, the dipole kSZ pattern expected for merging clusters could be used as a test in order to distinguish real cluster mergers from spurious ones.
The tSZ effect can be measured separately from the other effects due to its unique spectral dependence. The data of CMB temperature fluctuations will contain signatures of the RS effect we aim to detect, together with the following components that will be regarded here as contaminants: the primordial CMB field, the instrumental noise, the lensing of the CMB anisotropies induced by the matter present along the line of sight, and the kinematic SZ (kSZ) effect generated by the cluster's peculiar velocity with respect to the Hubble flow. As long as the contaminants show different spatial characteristics compared to the signal we are trying to unveil (the RS effect), it is possible to define a filter that optimally reduces the impact of those sources of noise. Regarding this, instrumental noise will be assumed to be scale-free white noise, i.e., to have the same power on all scales. The intrinsic CMB fluctuations usually have (if clusters are not too close) scales larger than those of the RS signal. On the other hand, the kSZ and the tSZ residuals are confined to the region of the cluster where the gas is located. Compared to these scales, the RS is usually broader.
However, if a contaminant shows spatial power on scales that are too close to those of the signal, only a more direct approach can be performed. This turns out to be the case for the anisotropies introduced by the lensing of the intrinsic CMB fluctuations. This component will require a separate processing (Sect. 4.1): from a template of the lensing-induced deflection field, we shall "deflect back'' the data, in an attempt to subtract the effect of gravitational lensing on the CMB. This then allows us to minimize the impact of this contaminant.
Once the cluster lensing signal is reduced by this "de-lensing procedure'', an optimized filter will be applied to suppress the other noise terms and to estimate the RS effect amplitude. This amplitude depends on the cluster infall velocity, which is essentially the physical parameter that can be measured with the RS effect (see Eq. (2)). This filter incorporates two templates in its definition: one for the RS effect described in Sect. (2), parameterized by the mass which we assumed to be estimated from observations of the tSZ effect, and one for the kSZ effect. The second one could be directly derived from the millimetric observations by observing at and filtering away the CMB features through a high-pass filter.
All templates are defined with the same hypotheses used when building the catalogues (spherical NFW haloes) and the cluster orientations are supposed to be perfectly known. This simplifies our analyses to a great extent, but leads to signal-to-noise ratios that do not account for errors in the morphological description of every system.
As we shall see, the resulting signal-to-noise ratio of a single merger is too small to provide a detectable signal. Thus it is necessary to average multiple measurements of many mergers. The detection of the cluster through their tSZ effect is already described well in other papers (e.g. Schaefer et al. 2004), so it will not be discussed here. We now describe in detail the steps of the observational strategy summarized above.
The CMB anisotropy produced by the gravitational lensing effect of galaxy clusters has a typical amplitude ten times larger than the RS anisotropy induced by the same clusters. These two effects arise from the same gravitational potential, so that their typical spectral scales are comparable and correlated in position and, although their patterns are not identical^{}, the similarity of their power spectra complicates the use of any Wiener filter.
It is therefore necessary to reduce the lensing effect before the application
of any optimized filter. To achieve this, we propose to remap the observed
data by applying a distortion that compensates the cluster lensing effect as
accurately as possible. To do this, it is necessary to invert the lens
equation
by introducing the new coordinate
.
In this way it is
possible to derive the de-lensed map
from the observed data Tthanks to the deflection field template
,
according to
In order to build a deflection field template, we define the merger using two NFW haloes according to Eq. (2). We assume that their masses are derived by the tSZ measurements, and for this reason we consider some uncertainty in the mass determination. The assumption of spherical haloes is justified since the deflection field is roughly symmetric on large scales. Furthermore, the final estimate of the signal will be produced after stacking all available merger events.
Of course, this procedure will - apart from the CMB - also distort the RS signal, the kSZ signal, and the distribution of the instrumental noise. Because of this, the RS effect and the kSZ templates used in the optimal filter definition have to be corrected through the same process. This procedure leaves some randomly distributed residuals in the field and some non-Gaussian features in the instrumental noise, which will be ignored in the following analysis.
Once the cluster lensing contribution is minimized, the remaining data contaminants in the maps are the CMB primary anisotropy, the kSZ effect, and the instrumental noise.
The Gaussian contribution of CMB and instrumental noise are fully described by their corresponding power spectra and can be filtered out by optimized filters. In particular, the CMB power spectrum shows an exponential cutoff around the arcminute scale due to Silk damping of primordial fluctuations. This angular scale is close to the typical one for distant galaxy clusters, and hence the CMB intrinsic fluctuations should make small contributions at the spatial frequencies we are interested in. For CMB experiments with very high sensitivity on small scales, such as ACT, SPT, or the Atacama Large Millimeter Array^{}, we assume that the noise is spatially white, i.e., that it introduces the same amount of power on all scales.
We present here an optimally matched filter (derived in detail in
Appendix A), which maximizes the signal-to-noise ratio by
processing the signal in the Fourier domain. A similar filter construction was
proposed by Haehnelt & Tegmark (1996) in order to extract the kSZ signal from single
clusters. However, that method considered the CMB and instrumental noise the
only sources of confusion when trying to recover the kSZ signal. Our filter
is defined in a slightly more general context. It provides an estimate of the
amplitude of a signal with a known spatial template and is to be measured
despite the presence of homogeneous noise with a known spatial power spectrum
and an arbitrary number of other contaminants following known spatial
patterns. Hence, the total observed data can be modeled by
(10) |
The Gaussian noise contribution
will contain the CMB temperature
fluctuations, the instrumental noise, and all the noise sources that can be
modeled as Gaussian random fields. In this case, the filter derived in
Appendix A takes the form
(12) |
(13) |
Figure 2: Two simple examples where the signal and a contaminant of known spatial pattern both have Gaussian profiles with the same width. In the upper panel, both signal and noise are positive (solid line). In the bottom panel the signal is positive and the contaminant is negative. | |
Open with DEXTER |
It is clear from Eq. (11) that the filter can be split into two components. The first, (T^{*}/P_{n}), maximizes the sensitivity at the spatial frequencies where the signal T is strong and the noise power spectrum P_{n} is small, as opposed to the second ( ), which underweights the regions where the Z_{i} components are important and correlated with the filter itself. Note that, since the filter is present in the definition of the coefficients, the final form of the filter must be computed iteratively.
We give two simple examples in Fig. 2 of how our filter behaves. Here the signal and the contaminant have the same Gaussian profile, and add in a constructive way (top panels) or in a destructive way (bottom panels). The estimator defined in Eq. (15) is the integral of the solid line times the dashed one. In the first (upper panel) case, the filter is negative where the contaminant is positive, so that their product is negative and the contaminant contribution is subtracted. In the second case (bottom panel), both filter and contaminant are negative so that their product is positive and the negative contribution of the noise is compensated. This actually shows that the same filter provides an unbiased estimator also if only the variance of the contaminant amplitude is known. That is, the sign of the contaminant is irrelevant to the filter, and this reflects the statistical nature of our filter.
As a more complex example, we show one of the realizations taken from our mock catalogue in Fig. 3. The merger was modeled by two NFW haloes at z=0.3, with masses and and with an infall velocity of at 45 degrees with respect to the line of sight (see Sect. 2). Since the two haloes are colliding, the RS effect (top left) of both components adds constructively, amplifying the resulting signal in the region between them. The kSZ effect (top right) shows the typical dipolar pattern reflecting opposite velocities along the line of sight. The anisotropy introduced by the lensing effect has a complex shape, since it depends on the configuration of both haloes and the CMB gradient (bottom left). The bottom right panel shows our filter. It is high-pass in order to suppress the CMB primary anisotropy noise contribution.
In Fig. 4, we plot the temperature amplitudes along the line containing the haloes displayed in Fig. 3. The top left panel shows the de-lensed RS effect superimposed on the filter. Note the sudden change in amplitude and sign of the filter close to the cluster cores. The top right panel displays the kSZ effect and the filter, whereas the bottom left one shows how the de-lensing procedure acts on the kSZ effect, shrinking its profile and shifting the peaks closer to the center. The bottom right panel compares our filter (new filter) with the filter originally proposed by Haennelt & Tegmark (1996, H&T 1996). The difference between these two filters is negligible when the contaminant given by kSZ effect (the Z template) is orthogonal to the signal, otherwise it becomes relevant.
The contribution from small merging substructures should be negligible, since our filter is a weighted integral over the whole field (Eq. (15)), and these minor mergers are odd functions of typical scales that are much smaller than the one of the filter.
Figure 3: Example taken from our mock catalogue, showing the spatial distribution of the RS signal ( top left), the kSZ effect ( top right), the anisotropy introduced by the lensing effect ( bottom left), and our filter ( bottom right). The temperature amplitudes of these panels along the line containing the haloes are displayed in Fig. 4. | |
Open with DEXTER |
The kSZ masks are applied as weight templates and are not directly subtracted from the input data. This makes the accuracy to which the templates are known uncritical. For this reason it is possible to use the observations at 217 GHz (the frequency at which the tSZ vanishes to a good approximation) as a good template for the kSZ. There might be some tSZ residuals, together with the lensing and RS signals, but their contribution is negligible for the filter construction. Again, this is because the filter does not subtract this template directly from the maps, which of course would erase all the RS signal we are searching for, but it only optimally weights down the regions contaminated by the kSZ effect (see "Case B'' in Sect. 6).
To test the described procedure, we performed simulations of the secondary anisotropies induced on the CMB by binary cluster mergers in a cosmological context. For the two merging haloes, we adopted a NFW profile for the dark matter component and a spherical isothermal sphere (SIS) for the gas distribution. We ignored tidal effects and the hydrodynamic processes involved in the interaction between the intra-cluster gas (ICM) of the two haloes. For our purpose, the model adopted gives a sufficiently realistic representation of RS (see Sect. 2), kSZ, lensing effect, and instrumental noise.
We adopted a standard CDM cosmology with a density contribution from dark matter, baryons and cosmological constant of , , and , respectively. The Hubble constant was set to with h=0.7(Bennett et al. 2003).
We simulated the CMB primary anisotropies as Gaussian random fields with a resolution of 1024^{2} pixels and a field of view of 120'. The CMB power spectrum used was computed using CMBEASY (Doran 2003), assuming a re-ionization fraction of 0.1 at redshift z=6.2 and a present helium abundance of 0.24. The field size ensures an adequate sampling of the CMB multipoles, and permits the retention of only the central quarter of the whole field, avoiding the boundary effects of the Fourier transforms.
For an isothermal cluster with negligible internal motions, the kSZ and the
tSZ effects are proportional to one another:
Figure 4: Comparison between the filter profile and the different components under consideration (see Fig. 3). The top left panel shows the de-lensed RS effect superimposed on the filter. The top right panel displays the kSZ effect and the filter, whereas the bottom left one shows how the de-lensing procedure acts on the kSZ effect. The bottom right panel compares our filter (new filter) with the filter originally proposed by Haennelt & Tegmark (1996, H&T 1996). | |
Open with DEXTER |
The lensing effect on the CMB is described by the lens equation:
The RS and kSZ effects were computed after introducing this deflection angle into Eq. (2). The distribution of cluster mergers in the universe is described according to the equations presented in Sect. 3. The number of clusters N_{M,z} was computed with that formalism on a finite number of mass-redshift cells. Then, for each cell, we randomly distributed N_{M,z} clusters with redshifts and masses enclosed in the cell range.
We included in our simulation the noise from a single, dish telescope with a
resolution of 1 arcmin and a sensitivity of
per beam. This
noise was modeled as a Gaussian field with the following power spectrum
(18) |
Table 1: Main characteristics of the upcoming millimetric observatories ACT, SPT, and our simulated instrument.
Finally, one of the inputs required to build our templates is the cluster mass, and therefore an error in this quantity will yield errors in the produced templates and errors in the deduced RS amplitudes. Therefore, the way we propagate the error in masses is by simply assigning wrong masses to the input of our procedure and letting this input be processed by our algorithm. These errors in cluster masses are introduced under a normal distribution in our Monte Carlo simulations with an uncertainty of .
In order to probe the ability of our extended filter to separate the RS signal, we prepared three different mock catalogues of cluster mergers following the recipe described in Sect. 4. These catalogues were generated using three minimum mass thresholds of , and and contain 6807, 266, and 28 cluster mergers, respectively.
The model (spherical NFW haloes) used to build the mock catalogues is the same assumed in defining the filter. This simplified approach leads to optimistic signal-to-noise ratios that should be considered as upper limits. In fact it could well be that a practical application of this procedure to real measurements leads to a weak or even non-significant signal.
We carried out two different analyses of the same set of mock catalogues, "Case A'' and "Case B'', where the filter was defined with the same model for the RS and lensing templates, but with two different templates for the kSZ contamination.
Even with a moderately sized cluster merger catalogue (266 clusters), the significance of our method approaches the 2- level in the detection of the RS signal. Since we are observing different mergers, and the errors in different mergers are uncorrelated, our Poissonian error-bar scale is the inverse square root of the number of mergers.
Table 2: RS amplitude estimates for three mass-selected cluster subsamples.
The results for Case B are also presented in Table 2. As one can see, RS signal detection is still possible. Comparing these results with those of the previous table, it is clear that the filter efficiency is not significantly affected. This is mainly due to the statistical nature of our filter which relies not on the accuracy of our templates for each merger, but on the accuracy of our templates when describing the average spatial properties of the kSZ and RS signals. Indeed, one should also be able to use the same average RS and kSZ templates for all mergers. The precision of our method is hence limited by the number of available mergers according to Poisson statistics.
We recall that these results have been obtained by ignoring the component of the RS signal of galaxy clusters given by the collapse of the environmental dark-matter in the potential well of the system. The inclusion of this component in the simulations would increase the signal level.
We have presented a method for extracting the Rees-Sciama effect of the major merger events of galaxy clusters when observing at centimeter/millimeter wavelengths. This would provide important kinematical information about the processes involved in structure formation. The RS effect is a secondary anisotropy of the CMB produced by the time variation in the gravitational potential along the line of sight. In the our case, the main contribution comes from the infall motion of the merging components, so that the total signal measures the convergence of the projected perpendicular momentum of the system (see Paper I).
These observations will be contaminated by the primordial CMB fluctuations, the kinematic SZ signal and the lensing effect from the clusters, together with the instrumental noise. We assumed that the thermal component of the SZ effect can be subtracted by multi-frequency observations, given that we know its frequency dependence.
The RS effect amplitude in one single merger event is much smaller than all the other components and thus only a statistical detection of the signal may be achieved. To maximize the signal-to-noise ratio and decrease the number of mergers that we need to co-add, we can apply a filter which enhances the RS signal above the other components. Therefore, we proposed in Sect. (4.2) an extended version of a matched filter presented by Haehnelt & Tegmark (1996). A complete derivation of our filter is given in Appendix A. It requires as input the power spectra of the CMB and the instrumental noise, plus two templates, one for the RS signal and another for the contaminating kSZ effect. These templates can be simple analytic models parameterized with the mass estimates provided by the tSZ measurements and with the expected average velocities. We would like to note that the filter construction is general so that it can be applied to other problems, like the extraction of gravitational lensing on CMB maps (e.g. Haehnelt & Tegmark 1996; Maturi et al. 2005a) or on background galaxies (Maturi et al. 2005b).
We showed that the developed filter can control the contamination from all signals except CMB lensing. This is because the RS and the lensing effect originate from the same gravitational deflection field, and thus their spatial frequencies are similar. Therefore, we included an additional step in our pipeline before the filter application, i.e. a "de-lensing'' procedure. We propose remapping the input data according to the inverse template deflection field, again parameterized by the tSZ mass estimate.
This method was applied to mock catalogues of cluster mergers where the merging components were assumed to be NFW spherical haloes. This was also assumed for the templates used in defining the filter. With this simplification, the signal-to-noise ratios achieved have to be considered as upper limits, because it could well be that ignored contaminants drive our suggested sigma detection level into a only weakly or even non-significant signal.
Our results show that the RS signal of merging clusters could be measured with this method after observing of the order of 1000 cluster mergers. These numbers will be achieved by combining the expected yields of the upcoming hig-resolution millimetric surveys (e.g. ACT, SPT, ALMA). Assuming a fraction of 30% for the number of mergers over the number of clusters, a RS detection will be feasible if future SZ surveys are able to perform sensitive observations of the order of 10^{4} clusters.
Acknowledgements
This work was supported in part by an EARA fellowship spent in the Max-Planck-Institut für Astrophysik and the COFIN 2001 fellowship provided by Bologna University. C.H.M. and J.A.R.M acknowledge the financial support of the European Community through the Human Potential Program under contract HPRN-CT-2002-00124 (CMBNET). C.H.M. is currently supported by NASA grants ADP03-0000-0092 and ADP04-0000-0093 and NSF grant PIRE 0507768. M.M. thanks Giuseppe Tormen for providing the simulation used in Fig. 1.
We aim to build an optimized filter
to extract from a noisy data
set
the best estimate of amplitude A of any signal with a known
spatial shape
:
We define the integral of a general function
as
This estimate is required for a bias
.
It has
a variance
given by
Since we are interested in deriving a filter that minimizes the estimate
variance, maintaining the unbias condition b=0, we introduce the Lagrangian
multiplier
and search for the filter function
that minimizes
the action
,
obtaining
The normalization factor
is obtained by substituting
Eq. (A.6) into Eq. (A.4), yielding
As is clear from Eq. (A.6), the filter has to be computed through an iterative procedure, where the starting point can be the filter with the variance matrix equal to the null matrix, which is essentially the filter originally proposed by Haehnelt & Tegmark (1996). This iterative procedure is required since the filter evaluation requires the correlation between the filter itself and the noise sources.
The simple intuitive interpretation of this filter is that the first term T/P_{n} maximizes the sensitivity on the spatial frequencies where the signal T is strong and the noise power spectrum P_{n} is small. At the same time, the second term introduces a weighted mask that takes into account the correlation between the filter and all the Z_{i}noise components, as well as the correlation between the different Z_{i} noise components. It is important to note that the filter not applies a mask by only subtracting the Z_{i} templates from the input data, leading to spurious residuals, but it down-weights the influence of the regions where the Z_{i} components are dominant and correlated with the filter itself. This means it is not critical to have a very accurate template for the noise sources, in order to apply the filter properly. It is sufficient to have a simple model that represents their average shapes.
The derived filter satisfies a generalized problem, where the noise sources
v_{i}Z_{i} can also be correlated, but the noise sources will be uncorrelated in
the majority of cases, either due to their different positions or to their
amplitudes
being mutuallyuncorrelated. In this case, the matrix
is diagonal and the filter function simplifies to
(A.8) |
(A.9) |
(A.10) |
(A.11) |