arxiv:astro-ph/ v1 12 May 2006 Accepted. Received ; in original form 15 May 2006

Transcription

1 Mon. Not. R. Astron. Soc. 000, 1 8 (2005) Printed 15 May 2006 (MN LATEX style file v2.2) Markov chain Monte Carlo analysis of ianchi VII h models M. ridges, 1 J.D. McEwen, 1 A.N. Lasenby, 1 M.P. Hobson 1 1 Astrophysics Group, Cavendish Laboratory, Madingley Road, Cambridge C3 0HE, UK arxiv:astro-ph/ v1 12 May 2006 Accepted. Received ; in original form 15 May 2006 ASTRACT We have extended the analysis of Jaffe et al. (2005a, 2006) to a complete Markov chain Monte Carlo (MCMC) parameter space study of the ianchi type VII h models including a dark energy density, using Wilkinson Microwave Anisotropy Probe (WMAP) cosmic microwave background (CM) data from the 1-year and 3-year releases. Since we perform the analysis in a ayesian framework our entire inference is contained in the multidimensional posterior distribution from which we can extract marginalised parameter constraints and the comparative ayesian evidence. Treating the ianchi CM anisotropy as a template, the parameter estimates derived for the total energy density, tightness and vorticity from 1-year data are found to be: Ω tot = , h = , ω = with Euler angles ( , , ), preferred by a factor of roughly unity in log-evidence over a concordance cosmology alone. For 3-year data this becomes Ω tot = , h = , ω = centred on ( , , ), preferred by a log-evidence difference of 1.9. This confirms that a ianchi type template is supported by the data. The low total energy density of this template, however, implies a geometry that is incompatible with cosmologies inferred from recent CM observations. Jaffe et al. (2005b) found that by extending the ianchi model to include a term in Ω Λ creates a degeneracy in the Ω m Ω Λ plane. We explore this region fully by MCMC and find that the degenerate likelihood contours do not intersect areas of parameter space that 1 or 3 year WMAP data would prefer at any significance above 2σ. Thus we can confirm the conclusion that a physical ianchi model is not responsible for this signature, which we have treated in our analysis as merely a template. Key words: cosmological parameters cosmology:observations cosmology:theory cosmic microwave background 1 INTRODUCTION It has recently been proposed that some of the current anomalous results seen in observations of the CM sky by WMAP, namely non-gaussian structure centered on the so called cold spot in the southern hemisphere (Vielva et al. 2005; Cruz et al. 2005), low quadrupole (Efstathiou 2004) and multipole alignments (de Oliveira-Costa et al. 2003), could be removed by allowing for large scale vorticity and shear. Components such as these arise in ianchi type VII h models (arrow et al. 1985) due to large scale rotation of the universe distorting the CM since the time of last scattering. This assumption has the unfortunate side-effect of violating universal isotropy and hence the cosmological principle, so any such result should be examined very carefully. Such a finding was indeed made by Jaffe et al. (2005a) [hereafter referred to as Jaffe] where a significant correlation was found between a class of ianchi VII h models and WMAP observations. However the resultant best-fit ianchi component was derived di- rectly from arrow et al. (1985) for a universe with a large negative curvature with total energy density Ω tot = 0.5 and including no dark energy component. Such a cosmology cannot be reconciled with cosmologies inferred from either current CM or other astronomical observations. In Jaffe et al. (2005b) the authors extended their study to include a cosmological constant and searched for a morphologically identical template to their first analysis in the increased parameter space. We perform a new study to explore fully the parameter space of both ianchi and cosmological parameters in order to select the best fitting template. We carry out this investigation using Markov chain Monte Carlo (MCMC) sampling in a fully ayesian manner. Hence we also find the comparative ayesian evidence allowing an assessment of whether the inclusion of such a template is supported by the data. m.bridges@mrao.cam.ac.uk

2 2 M. ridges et al. 2 MCMC ANALYSIS AND MODEL SELECTION FRAMEWORK Several recent studies have been performed using a ayesian approach to cosmological parameter estimation and model selection (Jaffe 1996, Drell et al. 2000, John & Narlikar 2002, Slosar et al. 2003, Saini et al. 2004, Marshall et al. 2003, Niarchou et al. 2004, asset et al. 2004, Mukherjee et al. 2005, Trotta 2005, eltran et al. (2005), ridges et al. 2005, Parkinson et al. 2006). Here we are interested in determining the best fitting ianchi VII h sky from WMAP data using a selection of all sky maps available; the 1-year WMAP Internal Linear Combination (WILC) map (ennett et al. 2003), the Lagrange Internal Linear Combination (LILC) map (Eriksen et al. 2004), the Tegmark, de Oliveira-Costa & Hamilton (TOH) map (Tegmark et al. 2003), the Maximum Entropy Method (MEM) map (Stolyarov et al. 2005) and the new 3-year WILC map (Hinshaw et al. 2006). We have incorporated ianchi sky simulations into an adapted version of the COSMOMC package (Lewis & ridle 2002) which we show via simulated sky maps is capable of extracting both the underlying cosmology and the ianchi signal, while also determining when the addition of such a component is favoured by the data. The ayesian framework allows the examination of a set of data D under the assumption of a model M defined by a set of parameters Θ with ayes Theorem: P(Θ D, M) = P(D Θ, M)P(Θ M), (1) P(D M) where P(Θ D, M) is the posterior distribution, P(D Θ, M) the likelihood, P(Θ M) the prior and P(D M) the ayesian evidence. We explore the parameter space by MCMC sampling to maximise the likelihood utilising a reasonable convergence criteria from the Gelman & Rubin R statistic (Gelman & Rubin 1992) so that once the Markov chain has reached its stationary distribution, samples from it reliably reflect the actual posterior distribution of the model. The final result will be parameter estimates marginalised from the multi-dimensional posterior and a value of the evidence for the model being studied. 2.1 Likelihood A CM experiment records the temperature anisotropy on the sky T (ˆp). We wish to find a cosmological model that can optimally describe these fluctuations by maximising the probability of observing the data T (ˆp) given a model described by a set of parameters Θ C. Then for a Gaussian model, from ayes Theorem we can express this probabilty, the likelihood, as: P(T (ˆp) Θ C) 1 C e 1 2 (T T C 1 T ) where C is the covariance matrix containing the predictions of the model Θ C and C its determinant. Expanding the temperature sky in spherical harmonmics T (ˆp) = lm a lmy lm, in a universe that is globally isotropic means C lm,l m = a lma l m = C l δ ll δ mm where the variance of the harmonic coefficients a lm is C l = 1 2l+1 m a lm 2. Thus Eqn. 2 becomes 1 : P({â lm } Θ C) lm (2) e ( â lm 2 /2C l) Cl. (3) 1 Where we have used the notation that Ĉl and a lm ˆ refer to observed quantities while C l and a lm refer to their theoretical counterparts. If one assumes universal isotropy the likelihood is then independent of m which can be summed over to give (for a full sky survey) the form used in most conventional CM analyses: ) (Ĉl ln P(Ĉl Θ C) = (2l + 1) + ln C l, (4) C l where Ĉ l is the observed quantity. However ianchi VII h models are anisotropic so one cannot simply compare observed Ĉl s with theoretical C l s, instead the full set of a lm s must be used. In this analysis we make the assumption that there is both a background cosmological contribution and ianchi VII h contribution to the CM. We take the cosmological component to be that of a standard ΛCDM universe, producing a set of a lm s, while the ianchi model component (defined by parameters Θ ) is given by a lm. The ianchi component can then be subtracted from the observed â lm s so that Eqn 3 becomes: P({â lm } Θ,Θ C) l 1 Cl e ( (â l0 a ) l0 )2 C l l m=1 2 C l e ( ) â lm a lm 2 2C l. (5) For numerical convenience we take the natural logarithm to yield the final likelihood function used in this analysis, hereafter referred to as L: ln L l (l + 1/2) ln(c l ) (a l0 a l0) 2 2C l 1 C l l m=1 ( a lm a lm 2), (6) which would reduce to Eqn. 4 if the ianchi signal were zero. 2.2 Evidence The ayesian evidence is the average likelihood over the entire prior parameter space of the model: L(Θ C,Θ )P(Θ C)P(Θ )d N Θ Cd M Θ, (7) where N and M are the number of cosmological and ianchi parameters respectively. Those models having large areas of prior parameter space with high likelihoods will produce high evidence values and vice versa. This effectively penalises models with excessively large parameterisations, thus naturally incorporating Ockam s razor. We have followed the analysis performed in our previous paper on primordial power spectrum determination (ridges et al. 2005) and so we only give a brief overview of the method here. In calculating the evidence our MCMC sampling procedure is fundamentally different from that used to produce the posterior distribution for parameter estimation purposes, where we were interested soley in tracing the areas of high likelihood. Now we require sufficient mobility over areas of low likelihood in order to accurately determine the average over the full prior. To do this we perform a reverse annealing schedule by drawing samples from: L(Θ C,Θ ) λ P(Θ C)P(Θ ), where λ is the inverse temperature and is lowered from 1 to 0. With this method log evidence values varying to within only one unit were successfully obtained from independent chains, while an error was calculated from the variance

3 MCMC analysis of ianchi VII h models 3 between chains. The total number of sampler calls made during the evidence burn in was typically 2000, deliberately large to ensure the chain was at a stationary point. During the reverse annealing schedule typically calls were made. 3 IANCHI VII H SIMULATIONS Jaffe implemented solutions of the geodesic equations developed by arrow et al. (1985) to compute ianchi induced temperature maps in type VII h cases. Jaffe et al. (2005b) extended these to include cosmologies with dark energy allowing an almost flat geometry, in keeping with the so called concordance cosmology. One of us (ANL, in preparation) independently performed the same extension and found results agreeing well with Jaffe et al. (2005b). The models are parameterised by six independent quantities; the total energy density Ω tot = Ω m + Ω Λ, the current vorticity ω, h which determines the tightness of the spiral and the Euler angles α and β determining the position of the centre of the spiral on the sky with γ allowing rotation of the pattern about this point 2. In addition one must specify the direction of rotation or handedness of the spiral which we have assumed in this analysis to be that found in Jaffe et al. (2005a). The ratio of the current shear to the Hubble parameter (σ/h) 0 (differential expansion) is determined by the combination of these parameters: ( σ H ) 0 = 6hΩ tot 2(1 + h) 1/2 (1 + 9h) 1 /2(1 Ω tot) ( ω H ) Jaffe found a best fit bianchi template T by a χ 2 analysis of the observed CM sky D for a selection of WMAP data variants by minimising: χ 2 = (D T ) T M 1 (D T ) (8) under the assumption of an embedded cosmology M given by the best fit 1-year WMAP angular power spectrum (Spergel et al. 2003). There is no contribution from noise in this expression as the dominant uncertainty on large scales is from cosmic variance. Their result is shown in Fig. 1 3 computed in our formalism with parameters Ω tot = 0.5, ω = , h = and position (42, 28, 310 ) in Euler angles. The primary feature is a central core positioned on the so called cold spot, assumed to be the source of a large area of non-gaussianity. The corrected map (see Fig. 2), produced after removal of this ianchi contamination, provides even by cursory visual inspection a more statistically homogeneous sky. In Section 4 we will perform our analysis on simulated data in order to test our parameter recovery and model selection algorithm before going on to examine real sky maps. The data was naturally created under the assumption that the same model universe (defined by Ω m, Ω Λ and Ω tot) was responsible for both the background ΛCDM cosmology and the ianchi VII h structure. Thus in sampling the parameter space we can describe the two elements as being coupled together as they always share the same set of energy densities and geometry. Of course our analysis method does not require this coupling, we also have the facility to vary both components independently, with seperate cosmological (Ω m, Ω Λ and 2 We adopt the active zyz Euler convention corresponding to the rotation of a physical body in a fixed coordinate system about the z, y and z axes by γ, β and α respectively. 3 See jdm57 for an animation of this template on the sphere. 0. Figure 1. est fit bianchi template as found by Jaffe et al. (2006) computed from parameters Ω tot = 0.5, ω = , h = and position (42, 28, 310 ) in Euler angles. Figure 2. Raw LILC map (left) with the ianchi corrected map (right) Ω tot) and ianchi energy densities (Ω m, Ω Λ and Ω tot) 4, allowing a powerful and efficient way to examine the CM sky for contaminant templates over the ΛCDM cosmology. We will adopt this method for the majority of the analysis in Section 5 where it will be shown that there is no ianchi VII h model that can be reconciled with a viable background cosmology. This of course does not render nugatory the examination of a signature on the CM sky that appears similar to a ianchi sky if one exists its significance should be quantified. When performed on real sky data the former method effectively extends the Jaffe study to fit simultaneously for a coupled ianchi and background cosmological component, while the latter method provides two main advantages over Jaffe s original template analysis: i) confirmation of the robustness of the detection within a varying cosmology thus accounting for possible degenerate parameters and ii) estimation of the ratio of ayesian evidences with and without a ianchi component providing an effective model selection to favour or disfavour the inclusion of a ianchi 4 SIMULATED DATA We produced a set of simulated high resolution CM sky maps with the HEALPIX 5 software routines. These contained a cosmic variance limited concordance cosmological component and a ianchi VII h component with fixed shear, matter and dark energy density and geometry (coupled to the cosmology) but varying the vorticity ω to study the sensitivity of the method to the total amplitude of the ianchi signal (of which ω is a good tracer). Eight separate 4 When we refer to decoupled ianchi energy densities we will always use a superscript. 5

4 4 M. ridges et al. Figure 3. ianchi component (enhanced by factor of 10 for illustration): ω = , h = 0.005, α = 42, β = 28, γ = 310, Ω m = 0.3, Ω Λ = 0.69 (left) and background cosmology: Ω m = 0.3, Ω Λ = 0.69, τ = 0.16, H 0 = 70 (right). Markov chains (yielding eight estimates of the evidence) were produced for two parameterisations: ianchi + cosmology (model A); and cosmology only (model ). The ianchi simulations were created with parameters: h = 0.005, Ω m = 0.3, Ω Λ = 0.69, α = 42, β = 28 and γ = 310 (see Fig. 3 left) combined with the cosmology (see Fig. 3 right) with eight equally spaced vorticities incremented between 3 and (see Fig. 4). Ultimately we aim to ascertain via these simulations the lowest amplitude of ianchi signal we can succesfully extract via our method. We cannot however directly relate the amplitudes of the Jaffe template with our simulated components using the vorticity since they have inconsistent values of shear, matter and dark energy density 6. It is however illustrative to compare the amplitude of their respective angular power spectra C l. Fig. 5. places the Jaffe component roughly between the ω = and ω = maps in overall amplitude. Marginalised parameter constraints (see Fig. 6) illustrate we are able to successfully extract the input cosmology and ianchi parameters for amplitudes at and above ω = , which suggests our method should be able to determine parameters of any ianchi component similar to the Jaffe template, should one exist. As for the question of whether such an inclusion is actually warranted by the data, we consider the evidence. A useful guide has been given by Jeffreys (1961) to rank relative evidence differences in order of significance: lne < 1 inconclusive, 1 < lne < 2.5 significant, 2.5 < lne < 5 strong and lne > 5 decisive. While there is marginal evidence in these simulations for a component at ω = (see Table 1) a definitive model selection can only be made at ω = , also roughly at the point where the ianchi template becomes visibly apparent in the sky maps (see Fig. 4). Since a template with ω = has an amplitude considerably larger than the template found by Jaffe we must be sceptical about making a decisive conclusion for the Jaffe template with real data. However the form of this template is sufficiently different to our simulated data for the study to be warranted. Figure 4. Simulated data of ianchi component in an embedded WMAP type cosmology (vorticity ω increases from top to bottom in column 1 then column 2) l(l+1)c l /2π l 10 ω = 4x10-10 ω = 5x10-10 ω = 6x10-10 Jaffe Template Figure 5. Comparison of the relative amplitudes of the angular power spectra of three simulated ianchi skies of varying vorticity and the Jaffe template. 5 APPLICATION TO REAL DATA We use the three full sky ILC maps available, downgraded to a HEALPIX resolution N side of 32 and smoothed with a large 330 beam as Jaffe did. This only allows the analysis to probe structure up to multipoles of 64, not limiting any study of ianchi structure, which contains few features beyond l = 20 (see Fig. 5), but precluding a full examination of the background cosmology. Thus we also analyse a map cleaned by the maximum-entropy 6 Our simulations were not carried out with a Jaffe template due to the inconsistency between the cosmology it required and the concordance model. method (MEM) (Stolyarov et al. 2005), which is capable of effectively separating and removing physical component contaminants. MEM produces an all-sky set of a lm s that are accurate and relatively noise free to l 300, thus allowing an improved examination of the cosmological parameter space. On the basis of Jaffe s previous work we will focus on three seperate classes of model: (i) A full ianchi VII h component coupled with a ΛCDM cosmology. (ii) A full ianchi VII h component varied independently of the cosmology.

5 MCMC analysis of ianchi VII h models 5 Table 2. Summary of full and reduced ianchi VII h and cosmology parameters and priors. Full ianchi Reduced ianchi Cosmology τ h ω 10 [x10 ] Ω Λ = [ 0.8,0.8] Ω Λ = [0,0] Ω b( H )2 = [0.005,0.1] Ω m = [0.1,0.99] Ω m = [0.1,0.99] Ω dm( H )2 = [0.01, 0.99] h = [0.01, 0.2] h = [0.01,0.2] H 0 = [40,100] ω = [0,20] ω = [0,20] τ = [0.01, 0.8] α = [0.1,1.8]rads α = [0.1,1.8]rads A s = [15, 55] β = [0.1,1.8]rads β = [0.1,1.8]rads n s = [0.5, 1.5] γ = [3.5,6.3]rads γ = [3.5,6.3]rads α β [Radians] γ Ω Λ Ω m H 0 Figure 6. Extracted cosmological and ianchi parameters from a simulated map with ω = (described in text) with input values (dotted lines). Table 1. Differences of log evidence with (model A) and without (model ) a ianchi component with vorticity incremented between ω [ ] A ±0.8-5 ± ±0.7-3 ± ±0.8-1 ± ± ± ± ± ± ± ± ± ± ± 2 (iii) A reduced ianchi VII h component with no dark energy density. The full ianchi VII h components are parameterised by Ω m,ω Λ, h, ω,α, β, γ, where the energy densities are identical to those used for the background cosmology (as discussed in Sec. 3). The independent component is parameterised by Ω m,ω Λ, h, ω,α, β, γ while the reduced component has Ω Λ = 0. The last case was motivated by the positive detection of such a component by Jaffe despite its alarmingly open cosmology Ω tot = 0.5 and their subsequent discussion of the viability of ianchi models in a concordance cosmology. The MCMC method provides the ideal medium in which to study the multidimensional parameter space of the former, where any degeneracies, such as those found by Jaffe, can be mapped efficiently. In addition we must make a distinction between full and reduced cosmological components depending on whether we use the high resolution MEM map or the low resolution ILC type maps. In the reduced case we describe the large scale CM fluctuations using the scalar amplitude A s, the spectral index n s, and the optical depth to reionisation τ, these being the only cosmological parameters that can be constrained at large scales. The full cosmological parameterisation requires in addition: the physical baryonic and dark matter energy densities Ω b (H 0/100) 2, Ω dm (H 0/100) 2 and the Hubble parameter H 0 7. We have summarised the parameters and priors in Table 2. In section 5.1 we search for any correlated full ianchi structure within a viable background cosmology, a full ianchi component decoupled from the cosmology and finally we attempt to reproduce the template found by Jaffe. In section 5.2 we use the ayesian evidence to select between competing ianchi components and to determine whether its inclusion is preferred by the data. 5.1 Parameter Constraints The findings of Jaffe et al. (2005a), Jaffe et al. (2005b) and Jaffe et al. (2006)) all point to a ianchi VII h component with a total energy density Ω tot = 0.5, this value being required to focus the centre of the spiral, thus alleviating, at least partially, the observable non-gaussianity (McEwen et al. 2005; McEwen et al. 2006) if centered in the southern hemisphere. However both 1-year and 3-year WMAP data would suggest a value of Ω tot much closer to unity, perhaps larger, motivating the search for a ianchi model within a more viable cosmology. The extension of the ianchi parameter space to include a dark energy density Ω Λ should thus be explored thoroughly to alleviate this tension. Jaffe et al. (2005b) found that a degeneracy is introduced with Ω Λ in the Ω m Ω Λ plane, similar to the geometric degeneracy which exists in the CM. Along this region, ianchi skies appear identical only being distinguished by their relative amplitudes. Exploration of this feature by MCMC sampling (Fig. 7.) shows that despite the increased degree of freedom in the ianchi model, at the 1σ level no overlap with either 1 or 3-year WMAP likelihood contours is found. Attempts to find harmonious areas of coupled Ω m Ω Λ space proved futile. ased on these results we can effectively rule out ianchi VII h models as the physical origin of this effect. Of course the corrections that a ianchi component makes to the WMAP data are no less interesting, whether we can model them or not. This is especially true given that at least the areas of non-gaussianity have been mostly preserved into the 3-year release (McEwen et al. 2006). In the remainder of this paper we examine the ianchi components decoupled from the cosmological parameters which weakens their significance, so that we assume that they 7 Where Ω m = Ω dm + Ω b.

6 6 M. ridges et al h ω [ ] Ω Λ Ω m α [rad] β [rad] Figure 7. Comparison of the Ω m Ω Λ ianchi degeneracy (shaded with 1 and 2σ contours) with the familiar CM geometric degeneracy from WMAP first (blue) and third year + polarisation (green) data (with 1 and 2σ contours) γ [rad] Ω m Figure 9. Extracted reduced ianchi VII h parameters from 1 (black) and 3 (red) year WMAP data with the Jaffe best fit values (dotted lines) h ω [ ] Table 3. Cosmological and ianchi parameterisations for each model tested α [rad] γ [rad] β [rad] Ω m Model Cosmology ianchi A τ, A s, n s - τ, A s, n s Ω m, Ω Λ, h, ω, α, β, γ C τ, A s, n s Ω m, h, ω, α, β, γ D τ, A s, n s, Ω c( H )2, Ω dm ( H )2, H 0 - E τ, A s, n s, Ω c( H )2, Ω dm ( H )2, H 0 Ω m, Ω Λ, h, ω, α, β, γ F τ, A s, n s, Ω c( H )2, Ω dm ( H )2, H 0 Ω m, h, ω, α, β, γ ates slightly broader constraints perhaps indicating that some of the 1-year ianchi template signal may have been noise driven Ω tot Ω Λ Figure 8. Extracted full ianchi VII h parameters from 1 (black) and 3 (red) year data with the Jaffe best fit values (dotted lines). provide a template, approximating some as yet unexplained structure on the sky. Exploring the full ianchi parameter space independently of the cosmology did allow reasonable constraints on the general morphology of the template (defined by h, ω and the Euler angles) though it left both Ω m and Ω Λ poorly constrained (see Fig. 8) illustrating the difficulty introduced by the degeneracy. We finally consider the case, as Jaffe did, of a reduced template with Ω Λ = 0. Marginalised estimates (see Fig. 9.) confirm very similar template fits among all four 1-year maps confirming the results of the Jaffe analysis. The 3-year data curiously gener- 5.2 Model Selection The conclusion that the ianchi VII h models may only be treated as templates does not preclude a consideration of the necessity of such a template to explain the data. Indeed, removal of areas of non-gaussianity and correction of large-scale power would explain serious anomalies in the current data, but at the cost of removal of universal isotropy. The ayesian evidence provides an ideal way to assess quantitatively whether such action is necessary. The models chosen are listed in Table 3 always using a ianchi component decoupled from the cosmology. We now explore the implications on model selection of the inclusion of the dark energy density in addition to the basic model template recovered above. Using the four downgraded ILC maps based on both 1 and 3-year WMAP data we compare models A, and C, while the higher resolution MEM map facilitates a comparison of models D, E and F. All models with a ianchi component (, C, E and F) are

7 MCMC analysis of ianchi VII h models 7 Table 4. Differences of log evidence using first and third year all sky maps with (models, C, E and F) and without (model A and D) a ianchi template component. Data \ Model A C 3-year WILC 0.0 ± ± ± year LILC 0.0 ± ± ± year TOH 0.0 ± ± ± year WILC 0.0 ± ± ± 0.7 D E F MEM 0.0 ± ± ± 1.1 supported by the both 1 and 3-year data over a concordance cosmology alone (see Table 4). However the significance of each result, as judged by Jeffreys scale varies somewhat over the data and model combinations. In general models and E show lower significance log-evidence differences when compared with C and F. This suggests that we have over complicated the model by the addition of Ω Λ, quantitatively supporting the poor constraints obtained and degeneracies introduced by this extension. It is interesting to note the MEM data provides higher significance results even within the weaker constraints of the full cosmological parameter space, substantiating the robustness of the template fit. The mean log-evidence for 3-year data is slightly higher for both models and C, albeit well within 1σ uncertainty. This is unlikely to be due to any improvement in the quality of data between years 1 and 3 as on the scales of interest (l < 10) the data is almost cosmic variance limited, however it could be due to improvements in the data pipeline leading to the ILC map. The 3 year WILC method (Hinshaw et al. 2006) differs from the 1 year method in a number of important ways. Firstly they now use the LILC approach (Eriksen et al. 2004) minimising variance by Lagrange multiplers, so comparison of the 1 year LILC and 3 year WILC are the most sensible. Secondly there is far better removal of galactic plane contamination which would have an effect on very large scales. On this basis it is conceivable that we are making a slightly more significant detection of a ianchi VII h template in 3-year data 8. 6 CONCLUSIONS We have investigated the possibility of signatures of universal shear and vorticity in 1 and 3-year WMAP data. We implemented ianchi VII h simulations in a fully ayesian MCMC analysis to extract the best fitting parameters of the model and to determine whether its inclusion is warranted by the data. We have made a significant detection for a ianchi template with parameters; Ω tot = , h = , ω = with Euler angles ( , , ) from 1-year data. Third year data prefers the marginally different parameters Ω tot = , h = , ω = centred on ( , , ), though detected at a slightly higher significance of 1.9 units in log-evidence. We conclude this discrepancy could be due to improvements in the Internal Linear Combination method used by the WMAP team for 3-year data. Since the open universal geometry 8 Our best fitting ianchi sky maps for 1 and 3 year data are shown in Appendix A and can be downloaded in FITS format from implied by these templates effectively rules out the ianchi VII h models as the origin of this structure we can only treat these models as template fittings, morphologically similar to some as yet unknown physical cause. Nevertheless our analysis shows that without their inclusion we cannot describe the observed CM structure satisfactorily, using only the concordance cosmology. In other words, at present the data suggests such an addition. ACKNOWLEDGEMENTS We acknowledge A. Slosar for permission to use his ayesian evidence code. This work was carried out largely on the COSMOS UK National Cosmology Supercomputer at DAMTP, Cambridge and we would like to thank S. Rankin and V. Treviso for their computational assistance. M was supported by a enefactors Scholarship at St. John s College, Cambridge and an Isaac Newton Studentship. JDM was supported by a Commonwealth (Cambridge) Scholarship. REFERENCES arrow J.D., Juszkiewicz R., Sonoda, D.H., 1985, MNRAS, 213, 917 asset.a., Corasaniti P.S., Kunz, M., Astrophys. J. Lett., 2004, 617, L1 eltran M., Garcia-ellido J., Lesgourgues J., Liddle A., Slosar A., 2005, Phys. Rev. D, 71, ennett C.L. et al. 2003, ApJS, 148, 97 ridges M., Lasenby A.N., Hobson, M.P., in press MNRAS (astro-ph/ ) Cruz M., Martinez-Gonzalez E., Vielva P., Cayon L., MNRAS, 356, 29 Drell P.S., Loredo T.J., Wasserman, I., Astrophys. J., 2000, 530, 593 Efstathiou G., MNRAS, 356, 29 Eriksen H.K., anday A.J., Gorski K.M., Lilje P.., 2004, ApJ, 612, 633 Gelman A., Rubin D., Statistical Science, 1992, 7, 457 Hinshaw, G. et al., 2006, submitted to ApJ Jaffe A., 1996, Astrophys. J., 471, 24 Jaffe T.R., anday A.J., Eriksen H.K., Gorski K.M., Hansen F.K., 2005a, submitted to ApJL (astro-ph/ ) Jaffe T.R., Hervik, S., anday A.J., Gorski K.M., 2005b, submitted to ApJ (astro-ph/ ) Jaffe T.R., anday A.J., Eriksen H.K., Gorski K.M., Hansen F.K., 2006, submitted to ApJL (astro-ph/ ) Jeffreys H., 1961, Theory of Probability, 3rd ed., Oxford University Press John M.V., Narlikar J.V., 2002, Phys. Rev. D, 65, Lewis A. and ridle S., 2002, Phys. Rev. D, 66, Marshall P.J., Hobson M.P., Slosar A., 2003, MNRAS, 346, 489 McEwen, J.D, Hobson, M.P., Lasenby, A.N., Mortlock, D.J., in press MNRAS (astro-ph/ ) McEwen, J.D, Hobson, M.P., Lasenby, A.N., Mortlock, D.J., submitted to MNRAS (astro-ph/ ) Mukherjee, P. Parkinson D. Liddle, A., 2005, ApJ, 638, L51 Niarchou A., Jaffe A., Pogosian L., 2004, Phys. Rev. D, 69, de Oliveira-Costa A., Tegmark M., Zaldarriaga M., Hamilton A., 2004, Phys. Rev. D, 69,

8 8 M. ridges et al. Figure A1. est fit bianchi template from WMAP 1 year data with parameters Ω tot = 0.47, ω = , h = 0.13 and position (42, 32, 301 ) in Euler angles. Figure A2. est fit bianchi template from WMAP 3 year data with parameters Ω tot = 0.45, ω = , h = 0.12 and position (44, 33, 298 ) in Euler angles. Parkinson D., Mukherjee P., Liddle A.R., astro-ph/ Saini T.D., Weller J., ridle S.L., 2004, MNRAS, 348, 603 Slosar A. et al., 2003, MNRAS, 341, L29 Spergel D.N. et al., 2003, Astrophys. J. Suppl., 148, 175 Stolyarov, V., Hobson, M.P., Lasenby, A.N., arreiro, R.., 2005, MNRAS, 357, 145 Tegmark M., de Oliveira-Costa A., Hamilton A., 2003, Phys. Rev. D, 68, Trotta R., 2005, astro-ph/ Vielva P., Martinez-Gonzalez E., arreiro R.,., Sanz J., L., Cayon L., 2004, ApJ, 609, 22 APPENDIX A: IANCHI VII H TEMPLATES We illustrate here our best fitting ianchi VII h templates found from 1 and 3 year WMAP data. They are both in agreement within 1σ, with that found by Jaffe Fig. 1 with some minor differences in the location and tightness of the spiral centre.