the universe started from a very hot , dense phase about 15 billion years ago in the big bang . the very early universe was opaque due to the constant interchange of energy between matter and radiation . the universe then cooled to a temperature of  4000k through its expansion about 300,000 years after the big bang . at this stage the electrons combine with the protons to form atoms , and the cross section for compton scattering with photons is dramatically reduced . the radiation from this point in time has been traveling towards us for 15 billion years and has now cooled to a blackbody temperature of 2.7 degrees kelvin . , the planck spectrum has its peak at microwave frequencies (  1  1000 ghz ) , and its study forms a branch of astronomy called cosmic microwave background astronomy ( hereafter cmb astronomy ) . in 1965 , arno penzias and robert wilson   were the first to detect this radiation . there must have been some bumps in the early universe that could grow to create the structures we see today . as the universe grows older , the observable universe gets bigger , due to the finite speed of light . the particle horizon of an observer is the distance to the farthest object that could have affected that observer . any objects further than this point are not , and never have been , in causal contact with the observer . at the last scattering surface the particle horizon corresponds to   2 as seen from earth today . no physical processes will act on scales larger than this . at the epoch of recombination , fluctuations on scales larger than 2 must have been produced by matter perturbations already present at this time . it is thought these fluctuations were probably laid down only 10 seconds after the big bang , and have their origin in quantum fluctuations of a scalar field . on smaller scales , corresponding to regions in causal contact at recombination , we should be able to see in the cmb the effects of physical processes occurring at recombination , such as acoustic oscillation of the coupled photon - baryon fluid , which thus also gives us a direct link with the physics of galaxy formation . in 1992 , the nasa cosmic microwave background explorer ( cobe ) satellite was the first experiment to detect the bumps . these initial measurements were in the form of a statistical detection rather than individual physical features ( some features in the first year maps were real cmb anisotropies . but it was not possible to distinguish these from the noise except in a statistical way . ) . today , experiments all around the world are finding these bumps , both causally connected and non - causally connected , that eventually grew into galaxies and clusters of galaxies . the main principle behind all of these experiments is that , instead of measuring the actual brightness , they measure the difference in brightness between different regions of the sky . the experiments at tenerife produced what was probably the first detection of the real , individual cmb fluctuations on scales comparable to the beam size of the experiment . there are many different theories of how the universe began its life and how it evolved into the structures seen today . each of these theories make slightly different predictions of how the universe looked at the very early stages which up until now have been impossible to prove or disprove . knowing the structure of the cmb , within a few years it should be possible for astronomers to tell us where the universe came from , how it developed and what will happen to it in the future . useful overall reviews of the physics of the production of cmb fluctuations , which will complement the informal presentation given in the next section , are contained in white et al . , 1994 , scott et al . , 1995   and hu et al . , 1995 this review is intended to be an extension and update of lasenby and jones , 1997 , and parts of that review are reproduced here so that the arguments are easily followed . we start from einstein s equation for the evolution of the universe 
 ( 1)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\dot r^2}(t ) = \frac{8}{3}\pi g{\rho _ { \rm{o}}}{r^ { - 1}}(t ) + \frac{\lambda } { 3}{r^2 } - k{c^2}$$\end{document } 
 where r describes the size of the universe , g is the gravitational constant , o is the present density of the universe , k is a measure of the curvature of space and  is the cosmological constant which can be thought of as the zero energy of a vacuum . if the cosmological term dominates ( as the scalar field is expected to at very high temperatures ) , then the other two terms become negligible and it is possible to solve and find 
 ( 2)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$r(t ) \propto \exp \left [ { { { ( \lambda /3)}^{\frac{1}{2}}}t } \right].$$\end{document } 
 therefore , the inflationary theory describes an exponential expansion of space in the very early universe . amplification of initial quantum irregularities then results in a spectrum of long wavelength perturbations on scales initially bigger than the horizon size . central to the theory of inflation is the potential v( ) , which describes the self - interaction of the scalar inflaton field . due to the unknown nature of this potential , and the unknown parameters involved in the theory , inflationary theory is bad ( at the moment ) at predicting the overall amplitude of the matter fluctuations at recombination . however , there is a reasonable agreement that the form of the fluctuation spectrum coming out of inflation should be given by 
 ( 3)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$|{\delta _ k}{|^2 } \propto { k^n},$$\end{document } 
 where k is the comoving wavenumber and n is the  tilt of the primordial spectrum . the latter is predicted to lie close to 1 ( the case n = 1 being the harrison - zeldovich , or an overdensity in the early universe does not collapse under the effect of self - gravity until it enters its own particle horizon when every point within it is in causal contact with every other point . the perturbation will continue to collapse until it reaches the jean s length , at which time radiation pressure will oppose gravity and set up acoustic oscillations . since overdensities of the same size will pass the horizon size at the same time , they will be oscillating in phase . these acoustic oscillations occur in both the matter field and the photon field and so will induce the level of the doppler peaks in the power spectrum depend on the number of acoustic oscillations that have taken place since entering the horizon . for overdensities that have undergone half an oscillation , there will be a large doppler peak ( corresponding to an angular size of  1 ) . other peaks occur at harmonics of this . as the amplitude and position of the primary and secondary peaks are intrinsically determined by the number of electron scatterers and by the geometry of the universe , they can be used as a test of the density parameter of baryons and dark matter , as well as other cosmological constants . prior to the last scattering surface , the photons and matter interact on scales smaller than the horizon size . through diffusion the electrons are coupled to the protons through coulomb interactions , and so the matter will move from high density regions to low density regions . this diffusion has the effect of damping out the fluctuations and is more marked as the size of the fluctuation decreases . it occurs when collisionless particles ( e.g. neutrinos ) move from high density to low density regions . if these particles have a small mass , then free streaming causes a damping of the fluctuations . the exact scale this occurs on depends on the mass and velocity of the particles involved . slow moving particles will have little effect on the spectrum of fluctuations as silk damping already wipes out the fluctuations on these scales , but fast moving , heavy particles ( e.g. a neutrino with 30 ev mass ) , can wipe out fluctuations on larger scales corresponding to 20 mpc today . putting this all together , we see that on large angular scales ( > 2 ) we expect the cmb power spectrum to reflect the initially near scale - invariant spectrum coming out of inflation ; on intermediate angular scales we expect to see a series of peaks , and on smaller angular scales ( < 10 arcmin ) we expect to see a sharp decline in amplitude . these expectations are borne out in the actual calculated form of the cmb power spectrum in what is currently the standard model for cosmology , namely inflation together with cold dark matter ( cdm ) . the spectrum for this , assuming  = 1 and standard values for other parameters , is shown in figure 2 . parameters assumed are  = 1 , n = 1 , h0 = 50 km s 
 mpc 
 and a baryon fraction of b = 0.04 . parameters assumed are  = 1 , n = 1 , h0 = 50 km s 
 mpc 
 and a baryon fraction of b = 0.04 . the quantities plotted are ( + 1)c , versus  where c is defined via 
 ( 4)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${c_\ell } = \left\langle { { { \left| { { a_{\ell m } } } \right|}^2 } } \right\rangle , \,\:\:\:\:\frac{{\delta t}}{t}(\theta , \,\phi ) = \sum\limits_{\ell , m } { { a_{lm } } } { y_{\ell m}}(\theta , \,\phi ) , $ $ \end{document } 
 and the ym are standard spherical harmonics . the reason for plotting (+1)c is that it approximately equals the power per unit logarithmic interval in . increasing  corresponds to decreasing angular scale  , with a rough relationship between the two of   2/ radians . in terms of the diameter of corresponding proto - objects imprinted in the cmb , a rich cluster of galaxies corresponds to a scale of about 8 arcmin , while the angular scale corresponding to the largest scale of clustering we know about in the universe today corresponds to 1/2 to 1 degree . the first large peak in the power spectrum , at  s near 200 , and therefore angular scales near 1 , is known as the , the inflationary cmb power spectrum plotted in figure 2 is that predicted by assuming the standard values of the cosmological parameters for a cdm model of the universe . in order for an experimental measurement of the angular power spectrum to be able to place constraints on these parameters , we must consider how the shape of the predicted power spectrum varies in response to changes in these parameters . in general , the detailed changes due to varying several parameters at once can be quite complicated . however , if we restrict our attention to the parameters  , h0 and b , the fractional baryon density , then the situation becomes simpler . perhaps most straightforward is the information contained in the position of the first doppler peak , and of the smaller secondary peaks , since this is determined almost exclusively by the value of the total  and varies as peak  . ( this behaviour is determined as mentioned above by the linear size of the causal horizon at recombination , and the usual formula for angular diameter distance . ) this means that if we were able to determine the position ( in a left / right sense ) of this peak , and we were confident in the underlying model assumptions , then we could read off the value of the total density of the universe . ( in the case where the cosmological constant was non - zero , we would effectively be reading off the combination matter + . ) this would be a determination of  free of all the usual problems encountered in local determinations using velocity fields etc . the height of the doppler peak is controlled by a combination of h0 and the density of the universe in baryons , b . we have a constraint on the combination bh02 from nucleosynthesis , and thus using this constraint and the peak height we can determine hg within a band compatible with both nucleosynthesis and the cmb . alternatively , if we have the power spectrum available to good accuracy covering the secondary peaks as well , then it is possible to read off the values of tot , b and h0 independently , without having to bring in the nucleosynthesis information . the overall point here is that the power spectrum of the cmb contains a wealth of physical information , and that once we have it to good accuracy and have become confident that an underlying model ( such as inflation and cdm ) is correct , then we can use the spectrum to obtain the values of parameters in the model , potentially to high accuracy . this will be discussed further below , both in the context of the current cmb data , and in the context of what we can expect in the future . the formation of topological defects is a generic property of symmetry breaking in unified field theories , and formation of such defects is expected to occur during phase transitions in the early universe . such phase transitions lead to spatial gradients in the value of the field which can not be eliminated , and the resulting field defects contain energy concentrations that can gravitationally perturb the surrounding matter and thus induce structure formation . depending on the properties of the field involved , the topological defects are strings , monopoles , walls and textures . if the seeds for galaxy formation are indeed provided by topological defects , we can expect the resulting power spectrum of cmb fluctuations to differ from that predicted by inflation . the full calculation of the power spectrum predicted by defects is much harder than that for inflation , and consequently the progress in this field has been relatively slow . qualitatively , however , the main difference from the inflationary scenario is that all perturbations have to be generated causally , that is , they must be within the horizon volume at a given epoch . thus anisotropies on scales above about 2 , as already measured by e.g the cobe satellite , have to be generated after recombination , and correspond to late - time effects . ( the horizon size grows after inflation , roughly like ( 1 + z ) . ) this entails , for example in the case of strings , the calculation of the properties of the string network through recombination to late times , which is computationally intensive . figure 3 shows a comparison of predictions for global strings , monopoles and textures made by pen , seljak & turok   with experimental results from present cmb experiments . the defect model predictions were each normalised to cobe at  = 10 . in practice , one is allowed to slide each of the curves up and down so as to best match a range of criteria , rather than ( as here ) matching a single large - scale c value . however , even with this freedom it is easily seen that the models do not fit the data very well . there appears to be evidence for a much larger acoustic peak in the data than predicted by defect theories . the same models can be used to predict the matter power spectrum ; this is shown in figure 4 . again , the models have been normalised to cobe , but a change in the normalisation does not appear to be sufficient to produce a good agreement between the predictions and experimental results . therefore , at present , the models considered by pen , seljak & turok appear to be ruled out by current experimental data , and the case for a topological defect origin of cmb fluctuations is less strong . the defects generally have more power on small scales than large scales relative to the adiabatic scdm model . the defects generally have more power on small scales than large scales relative to the adiabatic scdm model . however , recent work which uses an analytic expression for the evolution of local strings from the radiation to matter dominated era as well as including a non - zero cosmological constant , shows that it is possible to produce results that contain a broad peak in the angular power spectrum of the cmb ( figure 5 ) . however , this peak occurs at   400  600 which is in disagreement with the cat and ovro detections ( table 2 ) at these scales . the numerical method for calculating the matter and cmb power spectra used in this work differs from 
 that of pen , seljak and turok . are the results from simulations with 
  
 of 0.0 ( short - long dashed ) , 0.3 ( dotted - dashed ) , 0.5 ( short dashed ) , 0.7 ( long dashed ) , 0.9 ( solid line ) and standard cdm ( dot - dashed ) . table 2some recent cmb anisotropy measurements.experimenttypebeamwidth/beamthrow ( )v(ghz)tenerifeground5/810 , 15 , 33mit / firsballoon3.8180 + 3 higheracme / hemtground1.5/2.130/40maxballoon0.5/1.0110 , 180 , 270 , 360msamballoon0.5/0.6180 + 3 higherwhite dishground0.18/0.4790pythonground0.75/2.7590saskatoonground1.5/2.4526 - 36argoballoon0.9/1.8150 + 3 highercatground0.2513 - 16ovroground0.1 - 0.414.5 + 32 
  comparison between the observations and predictions for local strings . shown are the results from simulations with 
  
 of 0.0 ( short - long dashed ) , 0.3 ( dotted - dashed ) , 0.5 ( short dashed ) , 0.7 ( long dashed ) , 0.9 ( solid line ) and standard cdm ( dot - dashed ) . having given a description of the theoretical context , we now consider the experimental considerations that go into making measurements of the cmb . these include both the observing modes that are used to gain the sensitivity required and the various contaminants that emit at the frequencies of interest to cmb astronomy . there are a number of foregrounds that need to be well understood , or at the very least minimised , when making cmb measurements . one of the main foregrounds that is seen in cmb data originates from extra - galactic sources . the best analysis of the contribution by unresolved point sources to cmb experiments has been produced by franceschini it et al .   and de zotti et al . . they used numerous surveys , including vla and iras data , to put limits on the contribution to single beam cmb experiments by a random distribution of point sources . this analysis assumes that there are no unknown sources that only emit radiation in a frequency range between  30 ghz and  200 ghz . this range of frequency has not been properly surveyed and therefore there is still a cause of concern in the cmb community . the problem with using past surveys to predict the effects of point sources is that a large fraction of these sources are highly variable . therefore , by far the best way to subtract the sources from the data is to make simultaneous observations with a higher resolution telescope at the same observing frequency ( the ryle telescope at cambridge is used by cat for this purpose ) . at the higher frequency range of the microwave background experiments , this is the hardest galactic foreground to estimate , as it depends on the properties of the individual dust grains and their environment . at lower frequencies galactic synchrotron and free - free emission each of the galactic foregrounds are extremely difficult to eliminate , and so far the best method has been to observe in regions where their emissions are expected to be small and at frequencies where they are less dominant . with the multi - frequency observations it is possible to subtract the effect of these foregrounds if their spectral dependencies are known ( see section 6.1 ) . the final foreground that is seen with experiments looking at the microwave background is closer to earth than those already discussed . fluctuations in the atmosphere are hard to distinguish from actual extra - terrestrial fluctuations , when limited frequency coverage is available . space missions are the best way to do this , but their main problem is cost . high altitude sites ( either at the top of a mountain or in a balloon ) can reduce the atmospheric contribution as can moving the experiment to a region with a stable atmosphere . a cheaper alternative to physically moving as the atmospheric effects occur on a short time - scale , compared with the life - time of the experiment ( typically of order a few months for each data set taken with ground based cmb experiments ) , and the extra - terrestrial fluctuations are essentially constant , by integrating over a long time the contribution from the extra - terrestrial fluctuations are increased with respect to the atmospheric effects . stacking together n data points ( taken from n separate observations ) will reduce the variable atmospheric signal with respect to the constant galactic or cosmological one by a factor of \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\sqrt n $ $ \end{document } ( providing that they are independent with respect to the atmospheric signal and any long term atmospheric effects that affect the gain have been removed ) . the third way , which can also be combined with both the first and second way , is to design the experiment to be as insensitive as possible to atmospheric variations . the first obvious design consideration is to make the telescope sensitive to frequencies at which the atmospheric contribution is a minimum . by avoiding various bands in the spectrum , where much emission is expected ( for example water lines ) , the atmospheric effect is too large to allow useful observations from a ground based telescope . taken with the increasing foreground contamination from the galaxy at low frequencies ( it is expected that the galaxy dominates over the cmb signal at frequencies below 10 ghz . ) , this reduces the observable frequencies for ground based cmb experiments to between 10 and 100 ghz . this narrow observable range results in the need for balloon or satellite experiments so that a larger frequency coverage can be made to check the consistency of the results and to check the contamination from the various foregrounds that are expected . figure 6 shows the expected level of the various foregrounds for a typical cmb experiment . figure 7 shows a typical region of the sky over a range of frequencies covered by cmb experiments ( no atmospheric emission has been added ) . 
 figure 6the expected level of the various anisotropic foregrounds for a typical cmb experiment ( resolution  1 ) observing in the best low galactic flux regions . figure 7a typical region of sky over a range of frequencies covered by cmb experiments . synchrotron dominates at 1 ghz , free - free , synchrotron and cmb at 10 ghz , cmb at 100 ghz and dust at 1000 ghz . flux is in t equivalent temperature in k . no atmospheric emission has been added . the cmb model shown here is for a standard cold dark matter model with h = 0.5 and b = 0.05 . the expected level of the various anisotropic foregrounds for a typical cmb experiment ( resolution  1 ) observing in the best low galactic flux regions . synchrotron dominates at 1 ghz , free - free , synchrotron and cmb at 10 ghz , cmb at 100 ghz and dust at 1000 ghz . the cmb model shown here is for a standard cold dark matter model with h = 0.5 and b = 0.05 . the largest atmospheric variations occur mainly on relatively long time scales , compared to the integration time of telescopes ( typically of order a few minutes ) , as the variations are produced by pockets of air moving over the telescope . if an experiment could be insensitive to these long term variations , then it should effectively see through the atmosphere . an interferometer ( e.g. cat ) extracts a small range of fourier coefficients from the sky , reducing any incoherent signal ( short time scale variations ) or any signal that is coherent on large angular scales ( long time scale variations ) , and so should see through the atmosphere very well . similarly , an experiment that switches between two positions on the sky relatively quickly will also reduce the long term atmospheric variations . a variation on the beamswitching technique is to scan a single beam backwards and forwards across the sky ( e.g. saskatoon , python ) . in addition to the primary cmb anisotropies discussed so far , we also expect secondary anisotropies due to the interaction of cmb photons with clusters along the line of sight . the best known ( and observed ) is the sunyaev zeldovich ( sz ) effect , which is due to the upscattering of cmb photons by electrons in the hot gas ( te  10 - 10 k ) at the centre of clusters . other secondary anisotropies include , for example , the rees - sciama effect , which is a result of a cluster collapsing during its encounter with a cmb photon . in this review , we shall consider only the sz effect , since this the only secondary anisotropy to have been observed to date . the sz effect is in fact made up of two separate effects : one is due to the bulk velocity of the cluster , and the other due to the thermal velocities of the electrons in the cluster gas . the kinematic effect measures the cluster peculiar velocity , whereas the thermal effect can be used , in conjunction with images and spectra of its x - ray emission , to study the cluster gas . in particular , for a dynamically relaxed cluster , we can use the thermal sz effect and x - ray data to estimate the physical size of the cluster , and hence its distance . if the cluster has a peculiar velocity vp along the observer s line of sight , then the temperature of a cmb photon is doppler - shifted by an amount 
 ( 5)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${(\delta t)_{{\rm{kin } } } } = { t_0}\frac{{{v_p}}}{c}\int { { n_e } } { \sigma _ { \rm{t}}}{\rm{d}}l,$$\end{document } 
 where t0 is the temperature of the cmb , ne is the electron number density , t is the thomson scattering cross - section and the integral is taken along the line of sight . in terms of the rayleigh - jeans brightness temperature this becomes 
 ( 6)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${(\delta { t_{{\rm{rj}}}})_{{\rm{kin } } } } = { t_0}\frac{{{v_p}}}{c}\int { { n_e } } { \sigma _ { \rm{t}}}{\rm{d}}l\left\ { { \frac{{{x^2}{e^x}}}{{{{({e^x } - 1)}^2 } } } } \right\},$$\end{document } 
 or as an intensity 
 ( 7)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${(\delta { i_\nu } ) _ { { \rm{kin } } } } = \frac{{2{k^3}t_0 ^ 3}}{{{h^2}{c^2}}}\frac{{{v_p}}}{c}\int { { n_e } } { \sigma _ { \rm{t}}}{\rm{d}}l\left\ { { \frac{{{x^4}{e^x}}}{{{{({e^x } - 1)}^2 } } } } \right\},$$\end{document } 
 where we have set x = hv / kt0 . the effect due to thermal motions of the electrons is second - order in the electron velocity , and does not preserve the blackbody shape of the cmb spectrum . for the thermal sz effect the change in the rayleigh - jeans brightness temperature is given by 
 ( 8)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${(\delta { t_{{\rm{rj}}}})_{{\rm{thermal } } } } = { t_0}\int { { n_e } } { \sigma _ { \rm{t}}}\frac{{k{t_e}}}{{{m_e}{c^2}}}{\rm{d}}l\left\ { { \frac{{{x^2}{e^x}}}{{{{({e^x } - 1)}^2}}}[x\coth ( \frac{x}{2 } ) - 4 ] } \right\},$$\end{document } 
 where te is the electron temperature and me the electron mass . in terms of intensity this becomes 
 ( 9)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${(\delta { i_\nu } ) _ { { \rm{thermal } } } } = 2\frac{{{k^3}t_0 ^ 3}}{{{h^2}{c^2}}}\int { { n_e } } { \sigma _ { \rm{t}}}\frac{{k{t_e}}}{{{m_e}{c^2}}}{\rm{d}}l\left\ { { \frac{{{x^4}{e^x}}}{{{{({e^x } - 1)}^2}}}[x\coth ( \frac{x}{2 } ) - 4 ] } \right\}.$$\end{document } 
 it is usual to describe the magnitude of the thermal effect in terms of the y parameter , which is given by 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$y = \int { { n_e } } { \sigma _ { \rm{t}}}\frac{{k{t_e}}}{{{m_e}{c^2}}}{\rm{d}}l,$$\end{document } the frequency dependencies of the kinetic and thermal effects ( i.e. those functions in curly brackets in equations 69 ) , are shown in figure 8 . figure 8the frequency dependences of the thermal and kinetic sz effects expressed as a brightness temperature change ( top ) and intensity change ( bottom ) . the frequency dependences of the thermal and kinetic sz effects expressed as a brightness temperature change ( top ) and intensity change ( bottom ) . note that at v = 210 ghz , the maximum change in intensity due to the kinematic effect coincides with the null of the thermal effect . the magnitude of the thermal effect for a hot , dense cluster is ( trj)thermal  1 mk , and for reasonable cluster velocities the kinematic effect is an order of magnitude smaller . observations of the sz effect have been made in cambridge using the ryle telescope , which is an 8-dish interferometer operating at 15 ghz , in caltech using the owen s valley 5.5 m telescope , the owen s valley 40 m telescope   and the owen s valley millimeter array ( ovmmi )   , at nasa using the msam balloon experiment , at the caltech submillimeter observatory using a purpose built instrument called the sunyaev - zeldovich infrared experiment ( suzie )   and by various other groups . the magnitude of the observed sz effect in these clusters can be combined with x - ray data from rosat and asca to place limits on h0 . another important feature of the sz effect is that the decrement does not change with redshift . therefore , it should be possible to detect clusters out to very high redshift . to test this figure 9 shows the ryle telescope map of the sky towards the quasar pair pc1643 + 4631 ( jones et al . these two quasars are at a redshift of  3.8 and are  3 apart on the sky . they are a strong candidate for a gravitational lensed object due to the similarity in their spectra . as there is no x - ray detection with rosat the cluster responsible for the lensing must be at a redshift greater than 2.5 and from modelling of the gravitational lensing , or from fitting for a density profile in the sz effect , it must have a total mass of about 2  10m. 
 figure 9ryle telescope map of the sky towards quasar pair pc1643 + 4631 ( the positions of which are shown as crosses ) . ryle telescope map of the sky towards quasar pair pc1643 + 4631 ( the positions of which are shown as crosses ) . it is possible to predict the mass density of the universe . using the press - schecter formalism the number of clusters in terms of sz flux counts from this figure it is seen that , if the decrements are really due to the sz effect and there was no bias in selecting the fields ( e.g. the field was chosen because of the magnification of the quasar pair images ) , then the  = 1 model of the universe is ruled out . have been made , it is impossible to say how accurate these findings are . figure 10sz source counts with observational constraints , as a function of sz flux density expressed at 400 ghz . the two hatched boxes show the 95% onesided confidence limits from the vla and the rt ; due to the uncertain redshift of the clusters , there is a range of possible total sz flux density , which has for a minimum the value observed in each beam and a maximum chosen here to correspond to z > 1 . from the suzie blank fields , one can deduce the 95% upper limit shown as the triangle pointing downwards ( church et al . ) . we also plot the predictions of our fiducial open model (  = 0.2 ) for all clusters ( dashed line ) and for those clusters with z > 4 . sz source counts with observational constraints , as a function of sz flux density expressed at 400 ghz . the two hatched boxes show the 95% onesided confidence limits from the vla and the rt ; due to the uncertain redshift of the clusters , there is a range of possible total sz flux density , which has for a minimum the value observed in each beam and a maximum chosen here to correspond to z > 1 . from the suzie blank fields , one can deduce the 95% upper limit shown as the triangle pointing downwards ( church et al . ) . we also plot the predictions of our fiducial open model (  = 0.2 ) for all clusters ( dashed line ) and for those clusters with z > 4 . the above treatment of the sz effect is purely non - relativistic . in clusters where kbte > 10kev , the relativistic effect can be included either by extending the equations to include relativistic terms ( , , ) or by including multiple scattering descriptions of the comptonization process ( , , , , ) . order in temperature the correction for the thermal sz effect is given by 
 ( 10)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\begin{array}{l } 
 { \left ( { \delta { t_{{\rm{rj } } } } } \right)_{{\rm{thermal } } } } = { t_0}y\left [ { \frac{{{x^2}{e^x}}}{{{{({e^x } - 1)}^2 } } } } \right.[x\coth \left ( { \frac{x}{2 } } \right ) - 4 + \frac{{k{t_e}}}{{{m_e}}}\left (   -   \right.10 + \frac{{47}}{2}x\coth \left ( { \frac{x}{2 } } \right)\\ 
 - \frac{{42}}{5}{x^2}{\coth ^2}\left ( { \frac{x}{2 } } \right ) + \frac{7}{{10}}{x^3}{\coth ^3}\left ( { \frac{x}{2 } } \right ) + \frac{{7{x^2}}}{{5{{\sinh } ^2}\left ( { \frac{x}{2 } } \right)}}\left . { \left . { \left . { \left ( { x\coth \left ( { \frac{x}{2 } } \right ) - 3 } \right ) } \right ) } \right ] } \right ] 
 \end{array}$$\end{document } 
 and in the rayleigh - jeans limit ( small x ) we find 
 ( 11)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${(\delta { t_{{\rm{rj}}}})_{{\rm{thermal } } } } =   - 2{t_0}y\left [ { 1 - \frac{{17}}{{10}}\frac{{k{t_e}}}{{{m_e } } } + o\left ( { { { \left ( { \frac{{k{t_e}}}{{{m_e } } } } \right)}^2 } } \right ) } \right]$$\end{document } 
 and so it is seen that the inclusion of the relativistic treatment tends to lead to a small decrease in the sz effect . the hubble constant is inferred from combined sunyaev - zeldovich and x - ray data by a relation of the form 
 ( 12)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${h_{\rm{o } } } \propto \;\delta { t^ { - 2}}.$$\end{document } 
 the reduction in t in the rayleigh - jeans region , for given cluster parameters , leads to a decrease in the constant of proportionality in equation 12 , and hence a small reduction in the determined values of ho . for example , it is found that if the cluster temperature is 8kev the reduction in ho due to relativistic effects is 5% . since the discovery of a statistical anisotropy four years ago by the cobe satellite , there has been a great increase in the number of ground and balloon - based measurements of cmb anisotropy . many groups around the world have been engaged in projects which aim to detect individual features in the cmb , and to establish their statistical properties . table 1 lists all of the major experiments in the field and includes links to their home pages . the focus of these has tended to be towards smaller angular scales than measured by cobe , and this holds with it the exciting prospect of being able to detect structure in the primordial power spectrum . table 1 
 the main cmb experiments with links to their www home pages where available 
 space basedcobeplanckmapground basedtenerifesouth polesaskatooncatatcapythonovrosuziejodrell bankiac bartolwhite dishvipercobramatdasivsacbipolarbrown / wisc polarizationballoon bornefirsargomaxbeastmsambamaceqmapboomerangmaximatop hat 
 
 the main cmb experiments with links to their www home pages where available 
  space based cobe planck map ground based tenerife south pole saskatoon cat atca python ovro suzie jodrell bank iac bartol white dish viper cobra mat dasi vsa cbi polar brown / wisc polarization balloon borne firs argo max beast msam bam ace qmap boomerang maxima top hat a summary of the properties and results of several experiments , covering the period through to 1996 , is given in lasenby & hobson . since that time , i.e. in the last year , there have been significant developments with regard to a number of experiments , and we concentrate on these new results here . table 2 shows in a convenient form the main parameters of several recent experiments . the cosmic background explorer satellite ( cobe ) primordial anisotropy measurements are made using the dmr experiment , which consists of six differential microwave radiometers , two at each of 31.5 ghz , 53.0 ghz and 90.0 ghz . the first - year cobe observations provided convincing statistical evidence for the existence of cmb fluctuations . it was not however , possible to see individual cmb features , on the scale of the beam size , in the dmr maps , because even combining all of the maps together , the noise level per beam area was  45k and the signal to noise remained less than one . the results of the analysis of all four years of dmr data have now become available . , the results can be used to constrain the normalisation of a power law primordial spectrum . for a given slope n , normalisation is usually expressed via the implied amplitude of the quadrupole component of the power spectrum , c2 , as 
 ( 13)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${q_{rms - ps } } = t\sqrt { \frac{{5{c_2}}}{{4\pi } } } , $ $ \end{document } 
 where t is the mean cmb temperature . ( note that this qrms - ps value need not be the same as the actual quadrupole component . the fit is to a whole power spectrum as parameterised by a given n ) . for an assumed value of n = 1 , ( the harrison - zeldovich value ) , bennett et al . the joint best fit values of n and qrms - ps are n = 1.2  0.3 and qrms - ps = 15.3k . this restriction on the value of n is of course of great interest in the context of inflationary predictions that n = 1 . it is also of interest that inflation predicts gaussian fluctuations , and while this is much harder to test for than finding the amplitude and slope of the spectrum , the data are also consistent with this prediction statistical tests prefer gaussian over other toy statistical models by a factor of  5. with the accumulation of four years of data , the individual anisotopy features within the maps on the scale of the beam size are now becoming statistically significant . figure 11 shows the all - sky maps at each frequency taken from bennett et al . . some of the features in these maps away from the galactic plane are expected to be real cmb fluctuations , since the signal to noise in these regions is now about 2 sigma per 10 degree sky patch . indeed , features which repeat well between the different frequencies are now clearly visible . [ project collaborators : r. davies ( jodrell bank - principal investigator ) ; r. rebolo , c. gutierrez , r. watson , r. hoyland ( iac , tenerife ) ; a. lasenby , a. jones , g. rocha , s. hancock ( mrao , cambridge ) . ] a description of these experiments was given in lasenby & hancock . the tenerife experiments constitute a suite of three instruments , working at 10 , 15 and 33 ghz , designed and built at jodrell bank ( davies et al . 1992 , 1996 ) and operated by the iac in - situ on tenerife island . data are taken by drift scanning in right ascension at fixed declination and by sampling at 2.5 intervals in declination with a beamwidth of 5.1 ( fwhm ) it is possible to build up a fully sampled two - dimensional map of the sky at each frequency . initially , observations were concentrated on the strip of sky at declination + 40.0 where the low frequency surveys of haslam et al . ( 408 mhz )   and reich and reich ( 1420 mhz )   indicate a minimum in galactic foreground emission . analysis of the data at declination + 40.0 has been made and reported in hancock et al . 35 data scans at 10 and 15 ghz have recently been reported in gutierrez et al .   and are consistent with the dec . 40 data . a level of qrms - ps = 20  6 kwas found in a high galactic latitude region at dec . the comparison between this prediction and the tenerife 15 ghz data is shown in figure 12 , and it is seen that there is very good agreement with the main features in the cobe scan . this is evidence for features that are constant in amplitude over a frequency range of 15 ghz to 90 ghz and is a very strong candidate for a real cmb anisotropy . 35 data at 15 ghz ( solid line ) and the cobe dmr predictions of bunn et al . 35 data at 15 ghz ( solid line ) and the cobe dmr predictions of bunn et al . the tenerife programme is continuing , with the objective of mapping some 4000 square degrees of sky at 10 , 15 and 33 ghz . in conjunction with the cobe four - year data set , these experiments will continue to offer a useful source of large - scale cmb anisotropy measurements and hence to directly probe cosmological theories . in terms of power spectrum results , a revised estimate of the fluctuation amplitude for  values near  20 has been constructed from the dec 40 data , making an allowance for an atmospheric contribution ( hancock et al . + 30 to 45 ) has been made and will be presented in a future paper . for an assumed value of n = 1 a full two dimensional likelihood analysis gives qrms - ps = 22  4 k . a cooled hemt receiver with six channels is used to span the frequency range 26 to 46 ghz . the beam sizes used range from  1.5 at the lowest frequency to  0.5 at the highest . the analysis of observations of a 24 hour ra strip at declination + 85.1 is described in wollack et al . ( 1997 ) , exciting results have been presented which show not just a detection , but for the first time for a switched - beam instrument , evidence for the form of the power spectrum itself on the angular scales probed . the final map after 3 years of data from this experiment is shown in figure 13 and the power spectrum from this map will be used below , in comparison with theoretical predictions . we note here that there is currently an overall scaling uncertainty in the saskatoon results of 14% , due to calibration uncertainties . recent analysis of the saskatoon data ( knox , in prep . ) appear to show that the previous calibration is an underestimate of the true level of the saskatoon data . this would make the doppler peak even higher and lower the value of ho found below . figure 13saskatoon 3 year map showing region analysed as compared to the cobe full sky coverage . saskatoon 3 year map showing region analysed as compared to the cobe full sky coverage . the owen s valley radio observatory ( ovro ) is a ground - based , 40 m single dish , telescope . a hemt reciever operating at two frequencies ( 14.5 ghz and 32 ghz ) has recently been added to improve the sensitivity of the telescope . the resolution is 0.12. the original experiment put very strong constraints on galaxy formation scenarios   . 88 give t = 56 - 65 + 8.6k at   589 ( leitch , private communication ) . it consists of a single bolometer on a 75 cm telescope , operating at a frequency of 90 ghz with a 0.75 fwhm beam . three seasons of observations have already been analysed ( python i , python ii   and python iii ) . the python results can be separated into two distinct bins in  space ( one centred at   87 and the other at   170 ) . the results for these two bins are t = 60 - 13 + 15k and 66 - 16 + 17k respectively . python iv and v data have already been taken and the analysis should be published very shortly . duffett - smith , m. hobson , m. jones , a. lasenby , c. osullivan , g. pooley , r. saunders and p. scott . ] the cosmic anisotropy telescope ( cat ) is a three element , ground - based interferometer telescope , of novel design . horn - reflector antennas mounted on a rotating turntable , track the sky , providing maps at four ( non - simultaneous ) frequencies of 13.5 , 14.5 , 15.5 and 16.5 ghz . the interferometric technique ensures high sensitivity to cmb fluctuations on scales of 0.5 , ( baselines  1 m ) whilst providing an excellent level of rejection to atmospheric fluctuations . despite being located at a relatively poor observing site in cambridge , the data is receiver noise limited for about 60% of the time , proving the effectiveness of the interferometer strategy . the first observations were concentrated on a blank field ( called the cat1 field ) , centred on ra 08 20 , dec . + 68 59 , selected from the green bank 5 ghz surveys under the constraints of minimal discrete source contamination and low galactic foreground . recently observations of a new blank field ( called the cat2 field ) , centred on ra 17 00 , dec . accurate information on the point source contribution to the cat2 field maps , which contain sources at much lower levels , has been obtained by surveying the fields with the ryle telescope at cambridge , and the multi - frequency nature of the cat data can be used to separate the remaining cmb and galactic components . some preliminary results from cat2 have been presented in baker ( 1997 )   and the 16.5 ghz map is shown in figure 14 . clear structure is visible in the central region of this map , and is thought to be actual structure , on scales of about 1/4 , in the surface of last scattering . when interpreting this map , however , it should remembered that for an interferometer with just three horns , the  synthesised beam of the telescope has large sidelobes , and it is these sidelobes that cause the regular features seen in the map . in the full analysis of the data , these sidelobes must be carefully taken into account . figure 1416.5 ghz cat image of 60 x 60 area centred on the cat2 geld , after discrete sources have been subtracted excess power can be seen in the central 2  2 primary beam ( because the sensitivity drops sharply outside this area , the outer regions are a good indicator of the noise level on the map ) . 16.5 ghz cat image of 60 x 60 area centred on the cat2 geld , after discrete sources have been subtracted excess power can be seen in the central 2  2 primary beam ( because the sensitivity drops sharply outside this area , the outer regions are a good indicator of the noise level on the map ) . the flux density range scale spans 40 mjy per beam . for an interferometer , visibility space correlates directly with the space of spherical harmonic coefficients  discussed earlier , and the data may be used to place constraints directly on the cmb power spectrum in two independent bins in . these constraints , along with those from the other experiments , are shown in figure 15 . the solid line is the prediction ( normalised to cobe ) for standard cdm with o = 1.0 , b = 0.1 and h0 = 45 km s 
 mpc . the solid line is the prediction ( normalised to cobe ) for standard cdm with o = 1.0 , b = 0.1 and h0 = 45 km s 
 mpc . the next few years should bring great improvements in the quality of cmb data available , both as regards to accuracy , and to range of coverage of angular scales . for balloon experiments , these improvements will come through the use of long - duration balloon flights , launched e.g. in antarctica , and circling the pole , and from the use of arrays of detectors covering several frequency bands . specific projects already in the pipeline using these techniques are tophat , beast , ace and boomerang . on the ground , interferometers will play an increasingly important rle , with three new arrays now under development , the very small array ( vsa ) , the degree angular scale interferometer ( dasi ) , the cosmic background interferometer ( cbi ) and the millimeter anisotropy experiment imaging array ( maxima ) . we will concentrate here briefly on a new interferometer array to be built by cambridge and jodrell bank in the u.k . , and to be sited in tenerife ( the design of the three new interferometer arrays are similar ) , and then discuss the two new satellite projects recently selected . although the cat has already provided maps of cmb anisotropy on scales  0.25 , these are relatively poor as images due to the limited number of baseline lengths and pixels available . in fact , the cat is a prototype for a considerably more advanced instrument , the very small array ( vsa ) . the objectives of the vsa are to obtain detailed maps of the cmb with a sensitivity approaching 5k and covering a range of angular scales from 10 to 2. an artists impression of the vsa is shown in figure 16 . using two interchangeable t - shaped configurations of 1015 horn elements , simulations have shown that it is possible to obtain maps of suitable sensitivity over the desired range of angular scales . the planned instrument would be sited at the teide observatory , tenerife , at an altitude of 2400 m and would make observations between 28 and 38 ghz , to enable the galactic component to be estimated and removed . the good accuracy available over a scale range that is well - matched to the positions of the first and secondary doppler peaks in the power spectrum , should enable measurements of  and h0 to be made to an accuracy of better than 10% after 12 months of observations . we believe that this will be refined somewhat by better array configuration design and the use of proper models and secondary peak information ( all work in progress ) . in addition , simulations have shown that the proposed observing strategy will be quite sensitive to the non - gaussian features expected on these angular scales if ( e.g. ) textures or monopoles are the seed perturbations for galaxy formation ( maisinger , hobson , lasenby & turok ) . the instrument is currently under construction at cambridge and jodrell bank in the u.k . , and is hoped it will be operational by the year 2000 . two new satellite experiments to study the cmb have recently been selected as future missions . these are map , or microwave anisotropy probe , which has been selected by nasa as a midex mission , for launch in august 2000 , and the planck surveyor , which has been selected by esa as an m3 mission , and will be launched hopefully soon after 2005 . an artist s impression of the map satellite , which has five frequency channels from 30 ghz to 100 ghz , is shown in figure 17 . an artist s impression of the planck surveyor satellite , which combines both hemt and bolometer technology in 10 frequency channels covering the range 30 ghz to 850 ghz , is shown in figure 18 . a crucial feature of a satellite experiment is the potential all - sky coverage that it affords , and the ability to map features on large angular scales ( > 10 ) . 
 figure 17artist s impression of the map satellite . figure 18 
 artist s impression of the planck surveyor satellite ( formally cobras / samba ) 

  artist s impression of the map satellite . with such a large increase in the amount of data available from cmb experiments it is becoming increasingly important to improve the analysis techniques used . with multifrequency observations of the same patch of sky to high precision ( the satellite will cover the full sky ) it should be possible to extract information on the foregrounds and leave a one technique that has been shown to be robust in the analysis of a number of different experiments is maximum entropy ( ,   and ) . if we start with bayes theorem which states , given a hypothesis h and some data d the posterior probability pr(h|d ) is the product of the likelihood pr(d|h ) and the prior probability pr(h ) , normalised by the evidence pr(d ) , 
 ( 14)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\rm{pr}}(h|d ) = \frac{{{\rm{pr}}(h){\rm{pr}}(d|h)}}{{{\rm{pr}}(d)}}.$$\end{document } 
 if the instrumental noise on each frequency channel is gaussian - distributed , then the probability distribution of the noise is a multivariate gaussian . assuming the expectation value of the noise to be zero at each observing frequency , the likelihood is therefore given by 
 ( 15)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\rm{pr}}(d|h ) \propto \exp [ - { \chi ^2}(h)]$$\end{document } 
 where the  misfit statistic has been introduced . let us consider a discretised image hj consisting of l cells , so that j = 1 , , l ; we may consider the hj as the components of an image vector h. if we base the derivation of the prior on purely information theoretic considerations ( subset independence , coordinate invariance and system independence ) we are naturally led to the maximum entropy method ( mem ) . it may be shown   the prior probability takes the form 
 ( 16)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\rm{pr}}(h ) \propto \exp [ \alpha s(h,\:m)],$$\end{document } 
 where the dimensional constant a depends on the scaling of the problem and may be considered as a regularising parameter , and m is a model vector to which h defaults in the absence of any data . in standard applications of the maximum entropy method , nevertheless , the mem approach can be extended to images that take both positive and negative values by considering them to be the difference of two pads , so that 
 ( 17)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$h = u - v$$\end{document } 
 where u and v are the positive and negative parts of h respectively . in this case , the cross entropy is given by 
 ( 18)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$s(h,\:{m_u},\:{m_v } ) = \sum\limits_{j = 1}^l { \left\ { { { \psi _ j } - { m_{uj } } - { m_{vj } } - { h_j}\ln \left [ { \frac{{{\psi _ j } + { h_j}}}{{2{m_{uj } } } } } \right ] } \right\ } } , $ $ \end{document } 
 where \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\psi _ j } = { [ h_j^2 + 4{m_{uj}}{m_{vj}}]^{1/2}}$$\end{document } and mu and mv are separate models for each pad . the most probable image h is then just the result from finding the maximum probability or , equivalently , the minimum of  - s . it can be shown   that the lagrange multiplier  is completely defined in a bayesian way and any prior correlation information can also be incorporated into the analysis . also , the assignment of errors is straightforward in the fourier domain where all the pixels in the discretised image will be independent . simulated data taken by the planck surveyor satellite and used mem to reconstruct the underlying cmb and foregrounds . they used six input maps ( the cmb , thermal and kinetic sz , dust emission , free - free emission and synchrotron emission ) to make up the data and then added gaussian noise to each frequency . after using mem with the bayesian value for a and giving the algorithm the average power spectra of each channel , it was found that features in all six maps were recovered . without any prior power spectrum information it was found that only the kinetic sz was not recovered and all others were recovered to some degree ( the cmb and dust were almost indistinguishable from the input maps with residual errors of 6k and 2k per pixel respectively ) . figure 20 shows the results from mem as compared to the input maps for the case with assumed average power spectrum . it is easily seen that mem reconstructs both the gaussian cmb and the non - gaussian thermal sz effect very well . figure 20the left hand side shows the input maps used in the planck simulations for a cdm simulation of the cmb and the thermal sz effect . the right hand side shows the reconstructions obtained by mem . it is easily seen that mem does a very good job at reconstructed these two components . for comparison , the grey scales on the input maps are the same as on the reconstructed maps . the left hand side shows the input maps used in the planck simulations for a cdm simulation of the cmb and the thermal sz effect . it is easily seen that mem does a very good job at reconstructed these two components . for comparison , the grey scales on the input maps are the same as on the reconstructed maps . in the absence of non - gaussian sources it is possible to simplify the maximum entropy method using the quadratic approximation to the entropy ( hobson et al . ) . this has the same affect as chosing a gaussian prior so that equation 16 becomes 
 ( 19)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\rm{pr}}(h ) \propto \exp [ - \:\frac{1}{2}{h^\dagger } { c^ { - 1}}h]$$\end{document } 
 where c is the covariance matrix of the image vector h given by 
 ( 20)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$c = { h^\dagger } h.$$\end{document } 
 the solution to the analysis with this bayesian prior is called wiener filtering ( see , for example ,   and ) and has been applied to many data sets in the past when non - gaussianity could be ignored . the data d can be written as 
 ( 21)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$d = b \star h +   \epsilon $ $ \end{document } 
 where the convolution of the image vector h is with b , the beam response of the instrument and frequency dependance of h.  is the noise vector . in this case , the best reconstructed image vector ,  , is given by 
 ( 22)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\hat h = wd$$\end{document } 
 where the wiener filter , w , is given by 
 ( 23)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$w = c{b^\dagger } { ( bc{b^\dagger } + n)^ { - 1}}$$\end{document } 
 and n is the noise covariance matrix given by n = . it is important to note that not only the cmb signal but all the foregrounds are implicitly assumed to be non - gaussian in this method (   and ) . it will have become apparent in the preceding sections that the cmb data are approaching the point where meaningful comparison between theory and prediction , as regards the shape and normalisation of the power spectrum , can be made . this is particularly the case with the new availability of the recent cat and saskatoon results , where the combination of scales they provide is exactly right to begin tracing out the shape of the first doppler peak . ( if this exists , and if tot = 1 . ) before embarking on this exercise , some proper cautions ought to be given . first , the current cmb data is not only noisy , with in some cases uncertain calibration , but will still have present within it residual contamination , either from the galaxy , or from discrete radio sources , or both . experimenters make their best efforts to remove these effects , or to choose observing strategies that minimise them , but the process of getting really clean cmb results , free of these effects to some guaranteed level of accuracy , is still only in its infancy . secondly , in any comparison of theory and data where parameters are to be estimated , the results for the parameters are only as good as the underlying theoretical models and assumptions that went into them . if cdm turns out not to be a viable theory for example , then the bounds on  derived below will have to be recomputed for whatever theory replaces it . many of the ingredients which go into the form of the power spectrum are not totally theory - specific ( this includes the physics of recombination , which involves only well - understood atomic physics ) , so that one can hope that at least some of the results found will not change too radically . bearing these caveats in mind , it is certainly of interest to begin this process of quantitative comparison of cmb data with theoretical curves . figure 21 shows a set of recent data points , many of them discussed above , put on a common scale ( which may effectively be treated as \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\sqrt { \ell ( \ell   + 1){c_\ell } } ) $ $ \end{document } , and compared with an analytical representation of the first doppler peak in a cdm model . the work required to convert the data to this common framework is substantial , and is discussed in hancock et al . the analytical version of the power spectrum is parameterised by its location in height and left / right position , and enables one to construct a likelihood surface for the parameters  and apeak , where apeak is the height of the peak , and is related to a combination of b and h0 , as discussed above . the dotted and dashed extreme curves in figure 21 indicate the best fit curves corresponding to varying the saskatoon calibration by 14% . the central fit yields a 68% confidence interval of 
 ( 24)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$0.2 < \omega   < 1.5,$$\end{document } 
 with a maximum likelihood point of  = 0.7 after marginalisation over the value of apeak . incorporating nucleosynthesis information as well , as sketched above ( specifically the copi et al . bounds of 0.009  bh  0.02 are assumed ) , a 68% confidence interval for h0 of 
 ( 25)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\rm{30km}}\;{{\rm{s}}^ { - 1}}\;{\rm{mp}}{{\rm{c}}^ { - 1 } } < { h_0 } < 50\;{\rm{km}}\;{{\rm{s}}^ { - 1}}\;{\rm{mp}}{{\rm{c}}^ { - 1}}$$\end{document } 
 is obtained . generally , in the range of parameters of current interest , increasing h0 lowers the height of the peak . thus taking the saskatoon calibration to be lower than nominal , for example by the 14% figure quoted as the one - sigma error , enables us to raise the allowed range for h0 . by this means , an upper limit closer to 70 km s mpc is obtained . ) the best angular resolution offered by map is 12 arcmin , in its highest frequency channel at 90 ghz , and the median resolution of its channels is more like 30 arcmin . this means that it may have difficulty in pining down the full shape of the first and certainly secondary doppler peaks in the power spectrum . on the other hand , the angular resolution of the planck surveyor extends down to 5 arcmin , with a median ( across the six channels most useful for cmb work ) of about 10 arcmin . this means that it will be able to determine the power spectrum to good accuracy , all the way into the secondary peaks , and that consequently very good accuracy in determining cosmological parameters will be possible . figure 19 , taken from the planck surveyor phase a study document , shows the accuracy to which  , h0 and b can be recovered , given coverage of 1/3 of the sky with sensitivity 2  10 in t / t per pixel . we can see that the good angular resolution of the planck surveyor should mean a joint determination of  and h0 to  1% accuracy is possible in principle . these figures do not , however , take into account any reduction in sensitivity as a result of the need to separate galactic foregrounds from the cmb . nevertheless , simulations using a maximum entropy separation algorithm ( hobson , jones , lasenby & bouchet , in press ) suggest that for the planck surveyor the reduction in the final sensitivity to the cmb is very small indeed , and that the accuracy of the cosmological parameters estimates indicated in figure 19 may be attainable . figure 19expected capability of a satellite experiment as a function of resolution . the percentage error in recovering cosmological parameters from the cmb power spectrum is shown versus the resolution available . figure 22the contours show 50 , 5 , 2 and 0.1 percentile likelihood contours for pairs of parameters determined from fits to the cmb power spectrum . the figures to the left show results for an experiment with resolution fwhm = 1. those to the right for a higher resolution experiment with fwhm = 10 plotted on the same scale ( central column ) and with expanded scales ( rightmost column ) . this figure is taken from bersanelli et al . the percentage error in recovering cosmological parameters from the cmb power spectrum is shown versus the resolution available . the contours show 50 , 5 , 2 and 0.1 percentile likelihood contours for pairs of parameters determined from fits to the cmb power spectrum . the figures to the left show results for an experiment with resolution fwhm = 1. those to the right for a higher resolution experiment with fwhm = 10 plotted on the same scale ( central column ) and with expanded scales ( rightmost column ) . it is possible to formulate two models with similar power spectra , but different underlying physics . for example , standard cdm and a model with a non zero cosmological component and a gravity wave component can have almost identical power spectra ( to within the accuracy of the map satellite ) . to break the degeneracy more accuracy is required ( like the planck surveyor ) or information about the polarisation of the cmb photons can be used . this extra information on polarisation is very good at discriminating between theories but requires very sensitive polarimeters .
we present a brief review of current theory and observations of the cosmic microwave background ( cmb ) . 
 new predictions for cosmological defect theories and an overview of the inflationary theory are discussed . 
 recent results from various observations of the anisotropies of the microwave background are described and a summary of the proposed experiments is presented . 
 a new analysis technique based on bayesian statistics that can be used to reconstruct the underlying sky fluctuations is summarised . 
 current cmb data is used to set some preliminary constraints on the values of fundamental cosmological parameters  and ho using the maximum likelihood technique . 
 in addition , secondary anisotropies due to the sunyaev - zeldovich effect are described .