the study of the temperature anisotropies in the cosmic microwave background ( cmb ) can provide crucial information on the origin of the universe and help us to discriminate between different cosmological models .
recent observations in many cosmological fields tend to support the standard inflationary model , which predicts the anisotropies to represent an isotropic gaussian random field . the wilkinson microwave anisotropy probe ( wmap ) 1-year data ( bennett et al .
2003 ) , measured the cmb fluctuations with high accuracy . in a first approach ( komatsu et al .
2003 ) , these were found to be consistent with the gaussian predictions .
further analyses revealed asymmetries or non - gaussian features which have been confirmed in the 3-year wmap data ( hinshaw et al . 2007 ) .
some of these anomalies are : low multipole alignments ( de oliveira - costa et al .
2006 , copi et al .
2007 , land & magueijo 2007 , chiang et al . 2007 ) , north - south asymmetries ( eriksen et al . 2007 ) , structure alignment ( vielva et al . 2007 ) , non - gaussian features found with steerable ( wiaux et al . 2008 ) and directional wavelets ( mcewen et al .
2006 ) , and a very cold spot in the southern hemisphere ( cruz et al .
this cold spot ( cs ) is found at galactic coordinates ( @xmath1 ) , and has an angular radius in the sky of @xmath2 .
the spherical mexican hat wavelet ( smhw ) was used to detect this structure .
this spherical wavelet is an optimal tool to enhance non - gaussian features on the sphere .
the first smhw analysis ( vielva et al . 2004 ) was performed on 15 different wavelet scales and two estimators , skewness and kurtosis , were used .
the kurtosis deviated significantly from the range of values expected from gaussian simulations at scales around @xmath3 .
the spot was found to be the main cause of the deviation from gaussianity , using the area estimator ( cruz et al .
cruz et al . ( 2006 ) showed the incompatibility of the frequency independent signal of the cs with galactic foregrounds , and cayn , jin & treaster ( 2005 ) introduced two new estimators , namely @xmath4 and higher criticism , redetecting the cs .
cruz et al .
( 2007a ) confirmed the detection in the 3-year wmap data , using all the mentioned estimators . in the latter paper ,
the probability of finding an at least as high deviation in gaussian simulations was calculated without any _ a posteriori _ assumption .
a conservative estimate of this probability was 1.85% , and hence this spot is a very special feature whose origin deserves further investigation .
rth et al.(2007 ) detect the cs in the 3-year wmap data using scalar indices instead of wavelets .
several models were proposed in order to explain this anomaly : voids ( inoue & silk 2006 , 2007 , rudnick , brown & williams 2007 ) , second order gravitational effects ( tomita 2005 , tomita & inoue 2007 ) , anisotropic bianchi vii@xmath5 model ( jaffe et al . 2005 ) , finite cosmology model , ( adler , bjorken & overduin 2006 ) , asymptotically flat lemaitre - tolman - bondi model ( garcia - bellido & haugbolle 2008 ) and brane - world model ( cembranos et al . 2008 ) are some of them .
also a large - scale , non - gaussian angular modulation ( naselsky et al . 2007 ) has been suggested although this hypothesis is not based on any physical model .
the bianchi vii@xmath5 model was proven by jaffe et al .
( 2006 ) to be ruled out at high significance , whereas there is still no further evidence for the validity of the other explanations .
a full bayesian evidence analysis of that bianchi model can be found in bridges et al .
( 2007a ) .
bridges et al .
( 2007b ) show that the cold spot was likely to be driving any bianchi vii@xmath5 template detection .
recently , cruz et al . ( 2007b ) showed that the cs could be caused by a cosmic texture .
the amplitude and scale of the spot were consistent with that interpretation and the kurtosis of the data was compatible with the gaussian cmb plus texture model at all scales .
the predicted number of spots at a scale of @xmath3 or larger , is @xmath6 which is consistent with the observed spot .    in this paper
we study two other possibilities , namely the rees - sciama ( rees & sciama 1968 ) effect caused by voids as proposed by inoue & silk ( 2006 , 2007 ) , rudnick , brown & williams ( 2007 ) and the sunyaev - zeldovich effect produced by clusters ( hereafter sz effect , sunyaev & zeldovich 1970 ) already considered in cruz et al .
both of them could in principle produce temperature decrements as the cs . in addition
a dip ( rudnick , brown & williams 2007 ) in the nrao vla sky survey data ( nvss , condon et al .
1998 ) and the eridanus group of galaxies ( brough et al . 2006 ) match approximately the size and position in the sky of the cs
. however , the values required by the parameters of either model are not very realistic .
a bayesian framework is very convenient to decide which model best describes the data .
the prior theoretical and observational knowledge on the parameters is included in the bayesian analysis and can be decisive in some cases .
our analysis is performed on the recently released 5-year wmap data ( hinshaw et al .
the large - scale structure of the cmb anisotropies remains almost unchanged respect to the 3-year release ( komatsu et al .
2008 ) .
we will first review the data and the template fitting method in section 2 , then we recall in section 3 the results obtained for the texture model in order to compare with the sz effect in section 4 and the void model in section 5 .
finally we present our conclusions in section 6 .
we use the combined , foreground - cleaned q - v - w map ( hereafter wcm ) of the 5-year wmap data release ( hinshaw et al . 2008 ) .
the wcm is in the healpix pixelization scheme ( gorski et al .
2005 ) with resolution parameter nside = 64 .
since the scale of the spot is around @xmath3 this resolution is sufficient and reduces the number of pixels used in the template fitting .
the template fitting is performed in a circular area of @xmath7 radius centered at galactic coordinates ( @xmath1 ) .
although the angular radius of the cs is about @xmath3 we have to consider at least a @xmath7 radius patch to take into account the whole neighborhood of the spot for the fit .
the smhw convolves all the pixels in this region and they contribute in an important way to the detected structure .    in our template fitting we apply bayesian statistics to the data , @xmath8 , in order to find the best hypothesis , @xmath9 , describing the data through a set of parameters , @xmath10 .
bayes theorem states : @xmath11 where @xmath12 is the posterior distribution , @xmath13 the likelihood , @xmath14 the prior and @xmath15 the bayesian evidence .
the null hypothesis , @xmath16 , describes the data as an isotropic gaussian random field ( cmb ) plus the noise of the wmap data . the alternative hypothesis
, @xmath17 describes the data as cmb plus noise and an additional template , @xmath18 , given by a physical model as for instance a texture spot , the sz effect produced by a cluster or the rees - sciama effect generated by a local void . @xmath19
( with @xmath20 ) depends on the parameter set @xmath21 .    for @xmath19 ,
the probability density associated with the observed data is @xmath22 the evidence is hence the average of the likelihood with respect to the prior and has been extensively used for model selection in cosmology ( e.g. jaffe 1996 , hobson , bridle & lahav 2002 , hobson & mclachlan 2003 , mukherjee , parkinson & liddle 2005 , parkinson , mukherjee & liddle 2006 ) . a model having a large evidence shows a high likelihood for its allowed parameter space given the data .
we redefine our notation , with the likelihood @xmath23 , and the prior @xmath24 so the previous integral for the evidence , @xmath25 , reads : @xmath26    the likelihood function is @xmath27 , where @xmath28 is defined as : @xmath29 and @xmath30 is the generalised noise matrix including cmb and noise .
the template @xmath18 depends on the parameter set @xmath10 .
the calculation of the noise contribution to @xmath30 is straightforward since the noise in the wcm is uncorrelated and well known . in order to obtain the cmb part
, we calculate the correlation function for the wcm taking into account the pixel and beam effects . comparing the resulting evidence ratios obtained from the analytical expression for the wcm correlation function with those obtained
using simulations , the errors are negligible .    the posterior probability ratio : @xmath31 can be used to decide beween both hypotheses .
the alternative hypothesis is the more probable one when @xmath32 . if @xmath33 the null hypothesis is the preferred one to describe the data .
the a priori probability ratio for the two models , @xmath34 is set to the fraction of the sky , @xmath35 , covered by either texture , clusters or voids . in this manner
we compensate the _ a posteriori _ selection of the pixel where the template is centred , namely the center of the cs , ( @xmath1 ) .
unified theories of high energy physics predict the production of topological defects ( kibble 1976 , vilenkin & shellard 2000 ) during a symmetry - breaking phase transition in the early universe . depending on the order of the broken symmetry ,
different types of defects arise , such as domain walls , strings , monopoles or texture .
each of these extremely energetic events leave a characteristic pattern in the cmb .
for instance cosmic texture ( turok 1989 ) produce hot and cold spots in the cmb sky .
cruz et al . ( 2007b ) showed recently that the cs could be due to a collapsing texture .
a bayesian template fit was performed on the 3-year wmap data , using an approximated temperature profile based on the analytic result given in turok & spergel ( 1990 ) .
the a priori probability ratio for the two models , @xmath34 , was estimated in this case by the fraction of sky covered by texture , @xmath36 .
the evidence ratio was @xmath37 , hence the resulting posterior probability ratio was @xmath38 favoring the texture hypothesis versus the null hypothesis .
the scale and amplitude of the spot ( @xmath39 and @xmath40 ) were consistent with having being caused by a texture at redshift @xmath41 and energy scale @xmath42 .
the number of texture spots of scale @xmath39 or above is around @xmath43 consistent with the single observed spot .
the amplitude of the texture spot was shown to be overestimated by a selection bias .
texture spots with amplitude @xmath44 were simulated and added to cmb gaussian simulations . performing a template fit on these single extreme events ,
the inferred amplitude was in average @xmath45 , i.e. almost twice the true amplitude .
this means that the amplitude required by the cs is not in conflict with the upper limits on the amplitude calculated from the power spectrum ( urrestilla et al . 2007 ) .
moreover , bevis et al . ( 2008 ) and urrestilla et al . ( 2007 ) find a @xmath46 preference for a cosmological model including texture or strings , when considering only cmb data .
however this preference is reduced when including hubble key project ( freedman et al .
2001 ) and big bang nuclesosynthesis ( kirkman et al . 2003 ) constraints in the bayesian power spectrum anlysis .    in order to estimate the significance of the @xmath38 result
, the same fitting procedure was performed on 10,000 gaussian cmb simulations .
the most prominent spot was selected for each simulation and a texture profile was fitted to it . around @xmath47 of the simulations showed @xmath48 and the typical values where @xmath49 .
this result was found to be almost independent on the priors or the exact shape of the approximated temperature profile .    using the 5-year data release we obtain @xmath50 and @xmath51 .
the lower noise level allows a slightly better detection of the texture template in the new data release .
we have not repeated the frequentist estimation of the significance because the result was very similar to the 3-year one .
the probability of having an at least as extreme value ( i.e. @xmath52-value ) of @xmath53 assuming that the null hypothesis is true , should be again around @xmath54 .
the presence of texture would explain the excess of kurtosis observed at scales around @xmath3 .
the @xmath55 texture predicted by the model would be confused with gaussian cmb spots since the latter have their maximum power at these scales .
a number of follow - up tests can be performed to check the texture hypothesis .
for instance an all - sky template fitting looking for several texture can confirm or refute the alternative hypothesis .
also a polarization analysis can help to decide , since gaussian cmb spots have a correlated polarization pattern , whereas late texture do not show associated polarization .
in addition the b - mode texture power spectrum is expected to be larger than the typical inflationary one , having also a different shape ( urrestilla et al .
the hot gas trapped in clusters modifies the cmb spectrum through the sunyaev - zeldovich effect .
the high - energetic electrons present in the gas collide with the cmb photons , changing their original frequency .
two different components can be identified .
the former is produced by random thermal motion of the electrons ( thermal effect ) , and the second one , by the peculiar motion of the cluster with respect to the cmb ( kinetic effect ) .
the thermal sz effect is typically one order of magnitude higher than the kinetic one , and therefore we focus on the first one .
the effect on the cmb photons is frequency dependent and at frequencies below 217 ghz it produces a temperature decrement .
hence this effect is a possible candidate to explain an anomalously cold spot at the wmap frequencies which range from 23 ghz to 94 ghz .
the temperature decrement is given by : @xmath56 where @xmath57 is the frequency in dimensionless units , with @xmath58 and @xmath59 planck and boltzmann constants respectively , and @xmath60 .
the cluster comptonization parameter is defined as @xmath61 where @xmath62 and @xmath63 are the intracluster electron density and temperature , @xmath64 is the thomson cross section , and @xmath65 is the electron mass .
the integral is performed along the line of sight through the cluster . in hot clusters ( @xmath66 kev ) , the relativistic electrons lead to a slightly modified spectral behaviour of the thermal sz effect , but this is negligible in our analysis . also neglected is the kinetic sz effect of the cluster .
cruz et al . (
2006 ) showed that the flat frequency dependence of the cs is incompatible with being caused by a sz signal alone .
however a combination of cmb plus sz effect may explain the spot and could have a sufficiently flat frequency spectrum .
cruz et al .
( 2005 ) already reported that a group of local galaxies were present in the direction of the cs .
they argued that the amount of gas needed to explain the whole signal was too large for a small group of galaxies .
however the group of galaxies in the direction of the spot has been shown recently to be larger than previously thought .
brough et al .
( 2006 ) , define it as the eridanus _ supergroup _ with mass @xmath67 .
the angular size and position in the sky of eridanus are similar to those of the cs , as can be seen in figure  ( [ fig : patch ] ) .
the top panel shows the wcm data and the bottom one , the nearby galaxies catalog ( tully , 1987 )
. both are @xmath68 patches centered at ( @xmath1 ) and were convolved with the smhw at scale @xmath3 .
patch of the wmap combined map ( wcm , top ) and the nearby galaxies catalog ( bottom ) , centered at ( @xmath1 ) and convolved with the smhw at scale @xmath3 .
the units are wavelet coefficients of temperature in mk and number of counts for the wcm and nearby galaxies catalog respectively.,width=317,height=566 ]    however there are other , more prominent nearby groups which do not have an associated cmb spot , hence this could be a statistical coincidence , but it is worth testing through a template fit whether there is a hidden sz - like template in the wcm data .    the temperature profile as a function of the angle in the sky , @xmath69 , can be calculated assuming a @xmath70-model ( with @xmath71 ) for the gas - density profile ( see e.g. atrio - barandela & mcket , 1999 ) : @xmath72 where @xmath73 is the central electron density and @xmath74 is the core radius of the cluster . expressing the virial radius of the cluster @xmath75 in terms of the core radius as @xmath76 , eq .  [
[ yc ] ] can be rewritten as : @xmath77 where @xmath78r_c t n_c$ ] , and @xmath79 .
\label{phi}\ ] ] the parameter @xmath80 is the angle subtended by the core radius of the cluster , whereas @xmath52 is not a parameter but a function of mass @xmath81 and redshift @xmath82 .
typical values for clusters yield @xmath83 .
summarising , our temperature profile with parameters @xmath80 and @xmath84 reads : @xmath85 .
\label{dttfinal}\ ] ] eq .  [
[ dttfinal ] ] gives the effect produced by a single cluster through a particular line of sight . as the wcm is a combination of maps at different frequencies , we choose the central frequency , @xmath86 ghz .
we consider two alternative hypotheses in this context . in the first one
the spot would be caused by a rich cluster and in the second one by a nearby group of the nearby galaxies catalog .
we simulate the entire population of clusters with flux above 10 mjy ( at 353 ghz ) , using a monte carlo algorithm where masses and redshifts were assigned to each cluster according to the press - schechter prescription ( press & schechter 1974 ) .
we assume @xmath87 , @xmath88 and @xmath89 ( where @xmath90 and @xmath91 are the matter and vacuum energy density parameters , and @xmath92 the dispersion of the density field at scales of 8 mpc ) .
the cut in flux removes nearby small groups of galaxies and distant small clusters .
the redshifts of the simulated clusters range from 0.005 to 3.4 . for each cluster with
known mass and redshift the core radius is determined from the scaling relation : @xmath93 where @xmath94 is the mass of the cluster in @xmath95 @xmath96 .
it is important to note that the cut in flux does not have any negative effect in the results of this paper since the fluxes we are trying to explain are much larger than 10 mjy .
since the redshift and core radius are known for each cluster , we calculate the angle subtended by the core , @xmath80 .
we find the maximum angle of the 300,000 clusters to be @xmath97 .
therefore the cluster hypothesis has to be discarded for generating a @xmath3 spot .
returning to the particular case of the eridanus supergroup found in the nearby galaxies catalog , we can consider larger values of @xmath80 , due to its proximity and lower values of @xmath84 since it is not a cluster but a big group or an ensemble of groups .
in addition we assume lower values of @xmath52 because the gas could be rather dispersed since the gravitational attraction is lower than in big clusters . now
the alternative hypothesis is that a big group from the nearby galaxies catalog is causing the cs .
although the center of eridanus does not match exactly the center of the spot , we still center the template at ( @xmath1 ) , since the gas does not need to fit exactly the positions of the galaxies .
this assumption is the most favorable for the alternative hypothesis . if the cs was due to the sz effect we would expect high x - ray emission in the considered region .
we can use the rosat all - sky survey ( snowden et al . 1997 ) to estimate the prior on @xmath84 .
the rosat survey excludes some small regions of the sky , and one of them crosses the considered patch .
however the estimated average counts in the rosat r6 band are around @xmath98 cts s@xmath99 arcmin@xmath100 .
the @xmath70-models which are compatible with this background and at the same time show an extended sz signal of several degrees , range between models with a core radius of @xmath101 mpc and bolometric luminosities of @xmath102 erg / s and models with @xmath103 mpc and bolometric luminosities of @xmath104 erg / s .
these models would predict at most @xmath105 .
the nearest groups in the nearby galaxies catalog are at @xmath106 mpc distance and as we have shown , the upper limit for the core radius is @xmath103 mpc . in this case
we get @xmath107 rad , which we set as the upper limit for @xmath80 .
we set the lower limit on @xmath80 to @xmath108 , since the spot has @xmath3 radius and the pixel size is around @xmath108 for the chosen resolution .
hence we consider flat priors for parameter @xmath84 in the interval @xmath109 $ ] and for @xmath80 in @xmath110 $ ] radians .
setting @xmath111 we obtain a relative evidence of @xmath112 .
the fraction of the sky covered by nearby galaxies in the catalog , is @xmath113 therefore @xmath114 and the sz hypothesis is discarded again .
the best fit parameters are the upper limits of the priors , @xmath115 and @xmath107 rad .
the amplitude of the temperature decrement does not exceed @xmath116 and is hence negligible compared to the deep decrement .
therefore the @xmath112 result is not suggesting that the sz hypothesis is viable for explaining the spot .
as the template has a negligible effect on the data , @xmath117 , the likelihood is almost constant in the whole parameter range and as the priors integrate to unity we have @xmath118 .
this is a well - known feature of bayesian model selection . in any problem , if @xmath17 has additional parameters compared to @xmath16 , and the likelihood is insensitive to these parameters , then @xmath119 .
the results do not significantly depend on the value of @xmath52 , the chosen priors or the resolution .
the gravitational effect due to a supercluster can not strongly affect large angular fluctuations because these massive structures are very concentrated in a particular region .
on the contrary , voids are better candidates to generate large angular spots .
the gravitational effect produced by cosmic voids can produce temperature decrements in the cmb and has been proposed as a possible source of the cs ( inoue & silk 2006 , 2007 ) . in analogy to the cluster case with the eridanus supergroup , there is a dip in the nvss radio source catalogue that approximately matches the position of the spot .
rudnick , brown & williams ( 2007 ) , claim that the dip could be associated with a void at redshift @xmath120 which could be the origin of the cs , but affirm that the significance of this coincidence is unclear .
as the fraction of the universe filled by voids is around @xmath121 in volume at redshift @xmath122 ( colberg et al . 2005 ) , the fraction of the sky covered in projection , @xmath35 , is close to one and hence such coincidences are not rare .
the dip and the cs appear at different scales as can be seen in figure [ fig : wmap_nvss ] where we present a @xmath68 patch of the wmap and nvss data convolved with the smhw at scales @xmath3 and @xmath123 . at scale @xmath3 the cs reaches its maximum power in the wcm , whereas there is no significant dip visible in the nvss data . at scale @xmath123 where the cs is not an outstanding feature , a dip becomes visible but there are other similar dips with no associated cmb spots , therefore the association between dip and spot is unclear and could be a coincidence .
moreover , there are several other dips with similar depth outside the selected region hence the dip does not seem to be an outstanding feature of the nvss data .
mcewen et al . ( 2007 ) found several regions where both data sets were highly correlated and one of them was located close to the cs .
this could be due to the very cold temperature of the cs which gives a high correlation with any underdense region in the nvss data at scale @xmath3 and not to the relation between dip and spot since they have different scales ( see figure [ fig : wmap_nvss ] ) .
patch of the wcm ( left column ) and the nvss ( right column ) , centered at ( @xmath1 ) and convolved with the smhw at scale @xmath123 ( top row ) and scale @xmath3 ( bottom row ) .
the units are wavelet coefficients of temperature in mk and radio source number counts for the wcm and nvss respectively.,width=317 ]    rudnick , brown & williams ( 2007 ) also performed an order of magnitude calculation to estimate the radius of the void , @xmath75 , needed to generate such a cmb spot through the rees - sciama effect obtaining @xmath124 mpc in concordance with previous calculations ( inoue & silk 2007 ) .
the probability of observing such a void in standard theory is very small .
for instance , a void with @xmath125 mpc will occur with a probability of @xmath126 in an optimistic case ( rudnick , brown & williams 2007 ) .
typical observed values are @xmath127 mpc for a matter density contrast of @xmath128 ( hoyle & vogeley 2004 , patiri et al .
similar results are found in n - body simulations ( colberg et al . 2005 ) .
however it is worth performing a bayesian template fit for the void hypothesis .
martnez  gonzlez et al .
( 1990 ) and martnez  gonzlez & sanz ( 1990 ) , performed a calculation of the temperature anisotropy generated by a compensated void in an einstein - de sitter universe , @xmath129 where @xmath130 is the angle formed between the direction of observation and the direction defined by the centre of the structure , and @xmath131 is the distance to the centre of the structure .
they adopt the thin shell approximation , where the void is an empty spherical region surrounded by a thin shell containing the matter which would have been inside the void . outside the void the universe is represented by an einstein - de sitter model and the shell propagates approximately as @xmath132 and @xmath133 .
the units are @xmath134 , where @xmath135 is the horizon distance .
[ [ eq : dt_t_void ] ] applies to photons crossing the void .
the photons that do not cross the void are unaffected since there is no net mass excess .
introducing the above expressions for @xmath131 and @xmath141 in eq .  [ [ eq : dt_t_void ] ] we have our temperature template with parameters @xmath82 , @xmath75 in a flat @xmath136-cdm universe .
we choose the prior on @xmath75 in concordance with the void radius distribution found in the two - degree field galaxy redshift survey data ( fig . 3 in hoyle & vogeley 2004 ) .
we extrapolate these values for large @xmath75 through an exponential decay as shown in figure  ( [ fig : pr_rv1 ] ) .
this extrapolation is generous with big voids since the derived decay is faster than the exponential used .        for the redshift
, we choose a flat prior between 0 and 1 .
this assumption is again rather generous with large radii , since according to n - body simulations ( colberg et al .
2005 ) the void sizes should be smaller at high , unobserved redshifts .
we set the a priori probability ratio for the two models to unity since voids cover almost all the sky .
the obtained evidence ratio is @xmath118 and therefore the posterior probability ratio is @xmath142 .
this result is not suggesting that the void model is viable as the best fit template corresponds to @xmath143 mpc and @xmath144 which produces a temperature anisotropy of @xmath145 .
as in the cluster case the template is completely negligible , the likelihood almost constant , and the priors integrate to unity , therefore @xmath118 .
the contribution of voids with @xmath146mpc to the evidence integral is negligible and the results remain unchanged varying the upper limit of parameter @xmath75 for values above @xmath147mpc .
the big voids ( @xmath148mpc and @xmath149 ) described in rudnick , brown & williams ( 2007 ) which do produce an important temperature decrement are strongly penalized by the prior , despite showing higher likelihood values .
we redefine @xmath17 in order to test the big void hypotesis .
now only voids with @xmath148mpc and @xmath150 are considered .
we renormalise the priors in the new parameter ranges and recalculate the prior probability ratio since the big voids ( in case they exist ) are expected to cover only a small part of the sky . as shown above a void with @xmath125 mpc is expected to occur with a probability of @xmath126 ( rudnick , brown & williams 2007 ) . hence as a valid upper limit
we can assume that there is at most one big void with @xmath125mpc and @xmath151 .
the fraction of the sky covered by this void is @xmath152 .
therefore we assume @xmath153 . recalculating the evidence for the new parameter range
, we obtain @xmath154 which yields @xmath155 .
this result discards the possibility of explaining the cs through such a void .
hence the void hypothesis has to be rejected .
inoue & silk ( 2006 , 2007 ) proposed slightly different voids with very high values for the radius of the void , @xmath156 mpc and quasi linear density contrast values @xmath157 , while rudnick , brown & williams ( 2007 ) used @xmath158 .
inoue & silk ( 2007 ) recognize that the required values for @xmath75 are inconsistent with observations at low @xmath82 and n - body simulations .
such a void is a @xmath159 result in the standard theory with @xmath160 .
they consider that a percolation process of the typically observed small voids could form the large quasi linear voids .
however these voids have not been observed yet and as already mentioned , according to n - body simulations ( colberg et al .
2005 ) the void sizes should be smaller at high , unobserved redshifts .
in this paper we consider different possibilities to explain a non - gaussian cold spot found in the wmap data .
this spot is located at galactic coordinates ( @xmath1 ) , covers around @xmath161 in the sky , has been shown to be incompatible with gaussian simulations and is unlikely to be produced by foregrounds ( vielva et al .
2004 , cruz et al .
2005 , cruz et al .
2006 , cruz et al .
we compare three alternative hypotheses to explain the cs , namely texture , clusters and voids versus the null hypothesis .
the latter describes the data as an isotropic gaussian random field ( cmb ) plus the noise of the wmap data .
a bayesian analysis takes into account the prior observational and theoretical knowledge on the parameters and can be used to find the hypothesis best describing the data .
the posterior probability ratio , @xmath53 , is the quantity we use to compare the different hypotheses .
the alternative hypothesis is preferred if @xmath32 and disfavored otherwise .
cruz et al . (
2007b ) performed such an analysis for the texture model , finding that this model is favored versus the hypothesis of gaussian fluctuations with @xmath162 .
the predicted size , amplitude , temperature profile , and number of spots generated by texture is consistent with the observed spot . repeating this analysis for the 5-year wmap data we find @xmath163 due to the lower noise level .
the frequentist @xmath52-value associated to this result is around @xmath54 .
the finding of a local supergroup ( brough et al .
2006 ) and a dip in the nvss survey ( rudnick , brown & williams 2007 ) , which match the position of the cs could be a hint for the cluster or void hypotheses , but our detailed bayesian analysis reveals that these are statistical coincidences .
the fraction of the sky covered by clusters and voids is included in our analysis .
the large angular size and temperature decrement of the cs are hard to explain with conventional clusters or voids .
we would need either a local , rich cluster or a void with radius @xmath164 mpc located at redshift @xmath149 . in the cluster case
only a local group with compton parameter @xmath165 is observed in the required direction .
the angular scale is marginally consistent with the spot but the temperature decrement generated is only a few @xmath166 and hence can not account for the cs .
our bayesian analysis gives @xmath167 , discarding the cluster hypothesis .    for conventional voids with @xmath127
mpc we obtain @xmath142 because the likelihood is insensitive to the extra parameters .
a detailed analysis reveals that these voids produce negligible temperature decrements below @xmath116 .
renormalising the priors for voids with radius @xmath164 mpc , we obtain a posterior probability ratio of @xmath168 .
therefore the void hypothesis is not viable either . explaining the cold spot by a void with @xmath148 mpc ,
would introduce a @xmath159 anomaly in the standard model to solve the @xmath169 anomaly of the spot .
the solution would be worse than the problem .
hence we can conclude that the texture hypothesis is the most plausible explanation for the cs .
many scientists are uncomfortable with cosmic defects , because these failed to explain structure formation . however , cosmic defects are a generic consequence of symmetry breaking , which is an essential ingredient in unified theories of particle physics .
hence it is very plausible that they exist .
they might be formed at any scale below the grand unified scale , and would be consistent with current cosmological observations .
they are even compatible with inflation .
the only constraint is that they should make a subdominant contribution to structure formation ( see for instance urrestilla et al .
2007 ) .
the virtue of the texture hypothesis is that it is fully testable .
a texture would not lead to a correlated polarization pattern around the cold spot , but it would lead to a very specific gravitational lensing pattern .
in addition there should be more cmb texture spots at scales around @xmath170 .
although they are more difficult to detect because the gaussian cmb reaches its maximum power at these scales , we may be able to detect them through an all - sky bayesian analysis , which is already in progress .
mc , emg and pv acknowledge financial support from the spanish mcyt project aya2007 - 68058-c03 - 02 .
pv acknowledges financial support from the ramn y cajal project .
we acknowledge the use of the legacy archive for microwave background data analysis ( lambda ) .
support for lambda is provided by the nasa office of space science .
this work has used the software package healpix ( hierarchical , equal area and iso - latitude pixelization of the sphere , http://www.eso.org/science/healpix ) , developed by k.m .
grski , e. f. hivon , b. d. wandelt , j. banday , f. k. hansen and m. barthelmann ; the visualisation program univiewer , developed by s.m .
mingaliev , m. ashdown and v. stolyarov ; and the camb and cmbfast software , developed by a. lewis and a. challinor and by u. seljak and m. zaldarriaga respectively .
adler r. j. , bjorken j. d. , overduin j. m. , 2006 , gr - qc/0602102 .
atrio - barandela f. , mcket j.  p. , 1999
, apj 515 , 465 .
garcia - bellido j. , haugbolle t. , 2008 , jcap , 04 , 003 .
bennett c.l . ,
et al . , 2003 , apjs , 148 , 1 .
bevis n. , hindmarsh m. , kunz m. , urrestilla j. , 2008 , phys .
lett . , 100 , 021301 .
bridges m. , mcewen j.d .
, lasenby a.n . ,
hobson m.p . , 2007 , mnras , 377 , 1473 bridges m. , mcewen j.d . , cruz m. , hobson m.p . , lasenby a.n .
, vielva p. , martnez
gonzlez e. , 2008 , submitted to mnras ( arxiv:0712.1789 ) .
brough , s. , forbes d.  a. , kilborn v.  a. , couch w. , colless m. , 2006 , mnras , 369 , 1365 .
cayn l. , jin l. & treaster a. , 2005 , mnras , 362 , 826 .
cembranos j. a. r. , de la cruz - dombriz a. , dobado a. , maroto , a. l. , 2008 , ( arxiv:0803.0694 ) condon j.j .
, cotton w.d . ,
greisen e.w . ,
yin q.f . , perley r.a . , taylor g.b
. , broderick j.j . , 1998 , aj , 115 , 1693 .
chiang l .- y . , naselsky p.  d. , coles p. , 2007
, apj , 664 , 8 .
cruz m. , martnez
gonzlez e. , vielva p. , cayn l. , 2005 , mnras , 356 , 29 .
cruz m. , tucci m. , martnez
gonzlez e. , vielva p. , 2006 , mnras , 369 , 57 .
cruz m. , cayn l. , martnez
gonzlez e. , vielva p. , jin j. , 2007a , apj , 655 , 11 .
cruz m. , turok n. , vielva p. , martnez
gonzlez e. , hobson m. , 2007b , science , 318 , 1612 .
colberg j.m . ,
sheth r.k .
, diaferio a. , liang g. , yoshida n. , 2005 mnras , 360 , 216 .
copi c. j. , huterer d. , schwarz d. j. , starkman g. d. , 2007 , phys . rev .
d 75 , 023507 .
de oliveira - costa a. , tegmark m. , 2006 , phys . rev .
d 74 , 023005 .
eriksen h. k. , banday a. j. , grski k. m. , hansen f. k. , lilje p. b. , 2007 , apj , 660 , 81 .
freedman et al . , 2001 , ,
apj , 553 , 47 .
hinshaw et al . , 2007 , apjs , 170 , 288 .
hinshaw et al . , 2008 ,
submitted to apj , ( arxiv:0803.0732 ) hobson m , bridle s. , lahav o. , 2002 , 335 , 377 .
hobson m. , mclachlan , 2003 , mnras , 338 , 765 .
m. s. , 2004 , apj , 607,751 .
inoue k. t. , silk j. , 2006 , apj , 648 , 23 .
inoue k. t. , silk j. , 2007 , apj , 664 , 650 .
jaffe a. , 1996 , apj , 471 , 24 .
jaffe t. r. , banday a. j. , eriksen h. k. , grski k. m. , hansen f. k. , 2005 , apj , 629 , 1 .
jaffe t. r. , hervik s. , banday a. j. , grski k. m. , 2006 , apj 644 , 701 .
kibble , j. phys . , 1976 ,
a9 , 1387 .
kirkman d. , tytler d. , suzuki n. , omeara j.m .
, lubin d. , 2003 , apjs , 149 , 1 .
komatsu e. et al . , 2003 ,
apjs , 148 , 119 .
komatsu e. et al . , 2008 ,
submitted to apj , ( arxiv:0803.0547 ) .
land k. , magueijo j. , 2007 , mnras.378 , 153 .
gonzlez e. & sanz j. l. , 1990 , mnras , 247 , 473 . martnez
gonzlez e. , sanz j. l. & silk j. , 1990 , apj , 355 , 5 .
mcewen j.d .
, hobson m.p .
, lasenby a.n . ,
mortlock d.j . , 2006 ,
mnras , 371 , l50 mcewen j.d .
, vielva p. , hobson m.p .
, martnez  gonzlez e. , lasenby a.n . , 2007 , mnras , 376 , 1211 .
mukherjee p. , parkinson d. , liddle a. , 2005 , apj , 635 , 51 .
naselsky p.  d. , christensen p.  r. , coles p. , verkhodanov o. , novikov d. , kim j. , 2007 , apj submitted , arxiv:0712.1118 .
parkinson d. , mukherjee p. , liddle a. , 2006 , phys rev .
d 73 , 123523 .
patiri s.g .
, betancort - rijo j. , prada f. , klypin a. , gottlber s. 2006 , mnras , 369 , 335 .
press w.h , schechter p. , 1974 ,
apj , 187 , 425 .
rth c. , schuecker p. , banday a.  j. , 2007 , mnras , 380 , 466 .
rees m.j . ,
. , 1968 , nature , 217 , 511 .
rudnick l. , brown s. , williams l. , 2007 , apj accepted , astro - ph/0704.0908 .
tomita k. , 2005 , phys .
d , 72 , 10 .
tomita k. , inoue k.t .
, 2007 , arxiv:0712.1291 .
tully r. b. , 1988 , nearby galaxies catalog .
cambridge university press , 1988 , 221 p. turok n. & spergel d. n. , 1990 , phys .
letters , 64 , 2736 .
urrestilla j. , bevis n. , hindmarsh m. , kunz m. , liddle a. , 2007 , arxiv:0711.1842 .
vielva p. , martnez - gonzlez e. , barreiro r.  b. , sanz j.l .
, cayn l. , 2004 , apj , 609 , 22 .
vielva p. , wiaux y. , martnez - gonzlez e. , vandergheynst p. , 2007 ,
mnras , 381 , 932 .
vilenkin a. , shellard e.p.s . , 2000 , cosmic defects and other topological defects ( cambridge univ . press ) .
wiaux y. , vielva p. , barreiro r.  b. , martnez - gonzlez e. , vandergheynst p. , 2008 , mnras , 385 , 939 .


the influence of environmental degrees of freedom has been identified as the key concept in explaining the classical behavior of macroscopic systems in a quantum framework @xcite . according to this notion a preferred set of localized system states  called the pointer basis @xcite
is induced in the course of the interaction of the system with its surrounding .
most characteristically , any initial superposition of these pointer states gets rapidly mixed , while the only states who retain their purity for a long time are the pointer states themselves .
while the basic ideas behind the decoherence process seems to be settled , it still remains an open problem to understand the emergence , the dynamics , and the main properties of the pointer states for microscopic realistic environments .
several strategies have been proposed so far for determining the pointer basis given the environmental coupling .
in @xcite the suggestion was made to sort all pure states in the hilbert space according to their linear entropy production rate .
the pointer states are then identified with the states having minimal loss of purity .
similar results are obtained by the approach of @xcite which is based on a time evolution equation whose solitonic solutions are identified with the pointer states .
so far , this concept has been applied to the damped harmonic oscillator @xcite and to a free quantum particle coupled linearly to a bath of harmonic oscillators @xcite . there , the solitonic solutions of the corresponding nonlinear equation are coherent states and gaussian wave packets , respectively .
moreover , the decoherence to gaussian pointer states was proved to be generic for linear coupling models @xcite .    in this paper , we extend the analysis from linear models to a non - perturbative treatment of the interaction .
we focus on the model of collisional decoherence , which provides a realistic description of the decoherence process generated by an ideal gas environment .
notably , experiments with interfering fullerene molecules displayed a reduction of interference visibility in agreement with this model @xcite .
we derive the corresponding pointer states which are shown to form an overcomplete , exponentially localized set of basis states , and we prove decoherence to these states using the orthogonal unraveling of the master equation @xcite . this stochastic process on the one hand provides the statistical weights of the pointer basis , and on the other hand represents an efficient way of solving master equations which exhibit pointer states .
moreover , it explains the emergence of classical hamiltonian dynamics .
the main result was already announced in @xcite . here , we provide a more detailed explanation of the proofs and derivations . as an extension of @xcite , we prove the decoherence dynamics for a more general situation where the localization rate of collisional decoherence is in a non - saturated regime , and we utilize the relative entropy in order to illustrate the emergence of the statistical weights .
while the above results are derived within the framework of decoherence theory , they can also be applied to dynamic reduction models which propose a modification of the schrdinger equation by means of nonlinear and stochastic terms .
in fact , the observational consequences of the ghirardi - rimini - weber ( grw ) spontaneous localization model @xcite are equivalent to the ones of collisional decoherence , since they are described by the same master equation @xcite .
the present work therefore applies to the grw model . in particular
, it provides the corresponding pointer basis .
the structure of the article is as follows . in sect .  2 ,
we briefly review the notion of pointer states , and we summarize the method for determining the pointer states discussed in @xcite . in order to motivate the approach
, we consider a two level system and a dephasing process .
this method is then applied to collisional decoherence in sect .  3 , which provides a set of solitonic states to be regarded as ` candidate ' pointer states .
we show in sect .
4 that these solitons form an overcomplete basis of exponentially localized states , and we give an expression for their spatial extension .
moreover , we demonstrate that the ` candidate ' pointer states move on classical phase space trajectories if they are sufficiently localized . in sect .  5 , we briefly review the method of quantum trajectories , focusing in particular on the orthogonal unraveling of the master equation .
this stochastic process is then applied to collisional decoherence which allows us to show that the ` candidate ' states are indeed pointer states in the sense of the definition given in sect .  2 .
furthermore , we use the orthogonal unraveling to show that the statistical weights of the pointer states are given by the overlap with the initial state .
we present our conclusion in section 6 .
to motivate the definition of pointer states , let us consider the quantum dynamics of the damped harmonic oscillator .
its evolution can be described by a master equation in lindblad form defined by the standard hamiltonian @xmath0 , and a single lindblad operator @xmath1 , with associated rate @xmath2 @xcite .
we take as initial state a superposition of two quasi - orthogonal coherent states @xmath3 which satisfy @xmath4 , with @xmath5 .
it is then easy to show that for times larger than the decoherence time @xmath6 the solution of the master equation is well approximated by @xcite @xmath7 with @xmath8 .
thus , any coherent state remains pure during the damped time evolution , while any superposition of distinct coherent states decays into a mixture ( with a decay rate @xmath9 ) whose statistical weights are determined by the initial overlaps @xmath10 and @xmath11 . due to this property ,
the coherent states are to be identified with the pointer states of the damped harmonic oscillator .
the above observation serves as the starting point for the following definition of the pointer states for an open quantum system evolving according to a lindblad master equation @xmath12 .
one says that the system exhibits a pointer basis if its dynamics involves a separation of time scales , characterized by a fast decoherence time @xmath13 , such that for any time much larger than @xmath13 , the evolved state is well approximated by a mixture of uniquely defined pure states @xmath14 which are independent of the initial state @xmath15 , @xmath16 with @xmath17 .
following the above example , we further demand that for initial states @xmath15 , which are superpositions of mutually orthogonal pointer states @xmath18 , the probability distribution @xmath19 is given by the initial projections @xmath20 the pointer states @xmath21 initially form an ( overcomplete ) basis , and they may evolve in time , though slowly compared to @xmath13 .    the name _ pointer state _ was coined in @xcite due to its relevance for the physical description of a measurement apparatus . a measurement device which probes an observable @xmath22
is constructed such that macroscopically distinct positions of the pointer or indicator are obtained for the different eigenstates of @xmath22 . for a quantum system initially prepared in an eigenstate of @xmath22
the read - out will display the corresponding eigenvalue with certainty provided these pointer states remain pure during the time evolution . on the other hand , if the quantum system is prepared in a superposition of eigenstates of @xmath22 , we expect the pointer not to end up in a superposition of different read - out states , but rather to be at a definite position , though probabilistically , with probabilities given by the overlap ( [ eq : bornrule ] ) .
we emphasize that the importance of pointer states goes beyond the physics of measurement devices and the quantum - to - classical transition since they are also a practical tool for the solution of master equations . knowing the pointer states @xmath21 , their time evolution @xmath23 , and their probabilities @xmath24
, one can immediately specify the solution of the master equation for any initial state and times greater than the decoherence time .
since the decoherence time is generically much shorter than the system and dissipation time scales of the pointer state motion , this allows one to capture a large part of the system evolution without solving the master equation .      a practical way to obtain the pointer states @xmath21 , for a given environmental coupling ,
was discussed in @xcite .
we will illustrate this method by means of a two level system , suspect to a dephasing environment . the corresponding master equation in interaction picture ,
@xmath25 is characterized by the lindblad jump operator @xmath26 @xmath27 . in bloch representation , @xmath28 with bloch vector @xmath29 and pauli matrices @xmath30 , the solution reads @xmath31 thus , the bloch sphere @xmath32 , which represents the set of all pure states , is projected onto the z - axis in the course of the dephasing process .
this implies that the decohered state @xmath33 is a mixture of the eigenstates of @xmath34 ( denoted by @xmath35 and @xmath36 respectively ) , @xmath37 \mathsf{p_{\uparrow } } + { { \rm tr } }    \left [ \mathsf{p_{\downarrow } } _ { } \rho_0 \right ] \mathsf{p_{\downarrow } \ , }    .
\label{eq : decoherencedephasing}\end{aligned}\ ] ] the comparison with ( [ eq : decoherence ] ) shows that the north and the south pole of the bloch sphere ( corresponding to @xmath38 and @xmath39 ) form the pointer basis of the pure dephasing process .
since the solution of the master equation will not be at hand in general , one requires a method which yields the pointer states without the knowledge of @xmath40 . to motivate this
, we notice that the north pole of the bloch sphere is the asymptotic end point of the trajectory illustrated by the thick line in fig .  1 .
this trajectory is generated by an equation of motion with the following properties .
( _ i _ ) it is nonlinear because it must distinguish pointer states from their superpositions .
( _ ii _ ) it preserves the purity of pure initial states , i.e. an initial state which lies on the bloch sphere remains on the surface .
( _ iii _ ) the generated trajectory follows the exact solution of the master equation as close as possible . in order to find such an equation of motion ,
it is suggestive to minimize the distance of the initial increments , i.e. @xmath41 here , @xmath42 denotes the generator of the master equation ( [ eq : dephasing ] ) in bloch representation and @xmath43 is subject to the condition that the generated trajectory remains on the surface of the bloch sphere .
the solution of the above optimization problem reads , in spherical coordinates , @xmath44 since the sine is positive for @xmath45 , the solutions of these equations tend asymptotically towards a pointer state of the system , see fig .  1 .
the equator of the bloch sphere forms a set of unstable fixed points of ( [ eq : kugelkoordinaten ] ) .
let us now generalize the above argument to general markovian master equations . replacing the euclidean norm by the hilbert - schmidt norm @xmath46 in the space of operators , the generalization of ( [ eq : minimum ] ) to higher dimensional systems reads @xmath47 where the minimization is with respect to all evolution equations @xmath48 which propagate @xmath49 within the set of pure states , such that @xmath50 .
it can be shown that the general structure of an equation that evolves state vectors @xmath51 preserving their normalization has the structure @xmath52 , with @xmath53-dependent , hermitian and anti - hermitian mappings @xmath54 and @xmath55 .
this implies that the equation of motion for the projector @xmath56 must be of the form @xmath57 \right ] \,,$]where @xmath58 $ ] .
using this form , the optimization problem ( [ eq : minimum2 ] ) reduces to @xmath59 \right ] \|^{}_{\rm{hs } } \,_{}$ ] .
as shown in @xcite the solution is determined by the generator of the master equation , @xmath60 .
hence , the generalization of ( [ eq : kugelkoordinaten ] ) reads as @xcite @xmath61 \ , .
\label{eq : ng } \right .
\right.\end{aligned}\ ] ] motivated by the example in sect .
2.2 , one conjectures that the asymptotic solutions of ( [ eq : ng ] ) provide the pointer states in more complex systems as well .
it will be important below that eq .
( [ eq : ng ] ) is known also in another context : it corresponds to the deterministic part of the orthogonal unraveling of the master equation . as we will demonstrate in section 5
, one can use this specific unraveling to prove for a specific model that the asymptotic solutions of ( [ eq : ng ] ) indeed provide the pointer states .
in order to assess the nonlinear equation ( [ eq : ng ] ) in the context of a nontrivial environmental coupling we now apply it to the model of collisional decoherence @xcite . the latter describes the motion of a quantum test particle in an ideal gas environment and it accounts for the quantum effects of the scattering dynamics in a non - perturbative fashion .
the corresponding master equation has lindblad form , @xmath62 +    \int { \mathrm{d}}q { } \left ( \mathsf{l}_q \rho \mathsf{l}^{^{\dag}}_q -    \frac{1}{2 } \mathsf{l}^{^{\dag}}_q \mathsf{l}_q \rho - \frac{1}{2 } \rho
\mathsf{l}^{^{\dag}}_q \mathsf{l}_q \right ) \ , ,   \label{eq : lindblad}\end{aligned}\ ] ] where the jump operators are proportional to momentum kick operators , @xmath63 ( with position operator @xmath64 ) . the continuous label @xmath65 has the meaning of a momentum transfer experienced by the test particle with @xmath66 the corresponding distribution , @xmath67 ; @xmath2 is the collision rate of the gas environment .
the 1d equation of motion thus reads @xmath68 +    \gamma \int_{- \infty}^{\infty } { \mathrm{d}}q { } \ , g \left ( q \right )    e^{iq \mathsf{x } / \hbar } \rho e^{- iq \mathsf{x } / \hbar } - \gamma \rho \ ,    .
\label{eq : cd}\end{aligned}\ ] ] it leads to a localization in position space , i.e. to a loss of spatial coherence , as can be seen by switching to the interaction picture , @xmath69 , and the position representation , @xmath70 the decay rate of the spatial coherences is thus characterized by a _ localization rate _
@xmath71 which is related to the momentum transfer distribution @xmath72 by @xmath73 since the fourier transform of the distribution @xmath74 tends to zero for large distances @xmath75 , the localization rate saturates for large @xmath75 at the maximum value given by the collision rate , @xmath76 , which can be interpreted as the limit where one collision is sufficient to reveal the particles ` which path ' information .
this behavior is in sharp contrast to linear models , where the localization rate grows quadratically , and thus approaches infinity in the limit of a large separation @xmath75 .      in order to apply the nonlinear equation ( [ eq : ng ] ) to collisional decoherence ( [ eq : cd ] ) of a free particle , that is @xmath77 , we rewrite the projector equation ( [ eq : ng ] ) in vector representation and choose for @xmath42 the lindblad form ( [ eq : lindblad ] ) , which gives @xcite @xmath78 the expectation value @xmath79 is disregarded in the following , since it contributes only an additional phase .
now , we choose the jump operator of collisional decoherence , @xmath80 , and switch to position representation , which yields @xmath81 \left ( x \right ) \ , ,     \label{eq : ngcd}\\ \fl     \lambda \left [ \left| \psi_t \right|^2 \right ] \left ( x \right ) & = & \gamma    \left ( | \psi_t |^2 \ast \hat{g } \left ( x \right ) - \int^{\infty}_{- \infty }    { \mathrm{d}}y| \psi_t |^2 \left ( y \right )   \left ( | \psi_t |^2 \ast \hat{g }    \right ) \left ( y \right ) \right ) \ , .
\label{eq : lambda}\end{aligned}\ ] ] here , @xmath82 denotes the convolution and @xmath83 is the fourier transform of @xmath72 , i.e. @xmath84 .
the two summands in ( [ eq : ngcd ] ) have counteractive effects on the temporal evolution of the wave function : the coherent term leads to its dispersion , whereas the second , incoherent summand tends to localize the solution . in order to explain this localization
, we note that the second summand in ( [ eq : lambda ] ) is independent of @xmath85 .
this implies that the centered parts of the wave function , where the convolution @xmath86 exceeds the constant term in ( [ eq : lambda ] ) , get amplified , i.e. @xmath87 , whereas the tails of the wave function get damped , i.e. @xmath88 . as a consequence of these competing effects , solutions of ( [ eq : ngcd ] ) evolve towards solitonic states where both effects are in equilibrium , such that the state moves with fixed shape and constant velocity . as discussed above
, these solitons are candidates for the pointer states of collisional decoherence .    assuming the momentum transfer distribution @xmath74 to be a centered gaussian with variance @xmath89 , we can rewrite ( [ eq : ngcd ] ) in dimensionless form .
@xmath90 here we use the dimensionless variables @xmath91 and @xmath92 to define the dimensionless wave function @xmath93 .
notably , eq .  ( [ eq : dimensionless ] ) depends only on the single dimensionless parameter @xmath94 .
figure  2 shows a numerical solution of ( [ eq : dimensionless ] ) , where we choose a superposition of two counter propagating localized states @xmath95 as the initial state , @xmath96 .
as expected from the above discussion , the ( modulus of the ) solution converges to a soliton .
moreover , we find that the soliton inherits its initial position and momentum expectation value from that localized component @xmath97 of the initial state which has the greatest weight @xmath98 , @xmath99 .
similar observations are found for various other initial states .
we proceed to characterize the solitonic solutions of ( [ eq : ngcd ] ) . in sect .  4.1 ,
the consequences of the conservation of probability on the phase of the solitons are analyzed , allowing us to predict the asymptotic shape of the solitons in sect .
4.2 . in sect .
4.3 we estimate the spatial extension of the solitons , followed by the proof that they form a basis of the hilbert space in sect
.  4.4 . finally , in sect .  4.5 , we discuss the dynamics of the solitonic solutions in the presence of an external potential .      as demonstrated in the previous section , the nonlinear equation ( [ eq : ngcd ] ) exhibits solitonic solutions @xmath100 in the sense that the modulus of @xmath100 moves with constants shape and velocity , i.e. @xmath101 with @xmath102 and @xmath103 real . in this section
, we analyze the general structure of the phase @xmath104 , which will be relevant subsequently .
the time derivative of a solution @xmath105 of ( [ eq : ngcd ] ) , yields the continuity equation for @xmath106 , @xmath107 \left ( x \right )   \ , .
\label{eq : continuity}\end{aligned}\ ] ] plugging the solitonic solution ( [ eq : shape ] ) into ( [ eq : continuity ] ) , gives @xmath108 \left ( x - vt \right ) - v \partial_x \log f^2    \left ( x - vt \right ) & = & - \frac{\hbar}{m } ( \partial_x^2 g \left ( x , t    \right ) \nonumber\\    &   & \left .
+ \partial_x^ { } g \left ( x , t \right ) \partial_x \log f^2    \left ( x - vt \right )
\right ) \ , .
\label{eq : intermediate1}\end{aligned}\ ] ] here we have used that @xmath109 \left ( x \right ) = \lambda \left [ f^2 \right ] \left ( x - vt \right)$ ] , which follows from @xmath110 .
the time dependence of the left hand side of ( [ eq : intermediate1 ] ) corresponds to a spatial shift .
thus , also the right hand side of ( [ eq : intermediate1 ] ) must exhibit such a simple time dependence , which implies that @xmath111 where @xmath112 denotes the right hand side of ( [ eq : intermediate1 ] ) .
it follows that @xmath113 since this equation must hold for all @xmath114 and @xmath115 , we may assume that the equality holds already for the summands , such that @xmath116 \ , .
\label{eq : ableitung1}\end{aligned}\ ] ] therefore , the temporal and spatial dependence of the phase has the general structure @xmath117 with unknown functions @xmath118 and @xmath119 .      to explore the tails of the solitonic states
@xmath100 let us consider the form of ( [ eq : ngcd ] ) for asymptotically large positions .
it reads @xmath120 with @xmath121 a @xmath53-dependent , positive constant . inserting the solitonic form ( [ eq : shape ] ) into ( [ eq : asymptotic ] ) yields @xmath122 + vf    \left ( x - vt \right ) \nonumber   \\    &   & - \frac{\hbar}{m }
\left [ \partial_x f \left ( x - vt \right ) \partial_x    g \left ( x , t \right ) + f \left ( x - vt \right ) \partial_x^2 g \left ( x , t    \right ) \right]\nonumber   \\ & & - \gamma a_{\psi } f \left ( x - vt \right ) \ , .
\label{eq : long}\end{aligned}\ ] ] using ( [ eq : phase ] ) , we find that both @xmath123 and @xmath124 are only a function of @xmath125 , and accordingly , that also the left hand side of ( [ eq : long ] ) must be a function of @xmath126 . if follows that @xmath127 is at most linear in @xmath128 ( that is @xmath129 , with unknown constants @xmath130 and @xmath131 ) .
considering the real and imaginary part of ( [ eq : long ] ) separately , one obtains two coupled ( second order ) differential equations @xmath132 ^ 2 \right ) \
, ,     \label{eq : set2}\end{aligned}\ ] ] where @xmath133 and @xmath134 .
this set of equations has two unique solutions @xmath135 where the constant @xmath136 depends on the boundary condition for ( [ eq : set1 ] ) ( which can be determined only by solving the full nonlinear equation ( [ eq : ngcd ] ) ) .
the solution with the positive exponent in ( [ eq : f ] ) is irrelevant , since it is not normalizable .
figure 3 confirms that the tails of the numerically obtained solitonic solutions of ( [ eq : dimensionless ] ) are in agreement with the functional form ( [ eq : f ] ) ; they are straight lines in the semi - logarithmic plot .
this shows that , unlike in linear models @xcite where the pointer states are gaussian , the pointer states of collisional decoherence are exponentially localized .
an important characteristic of the pointer states is their spatial extension .
as explained in @xcite , the latter can be related to the experimentally accessible one - particle coherence length of a thermal gas .
we will determine the pointer width in this section , and apply the result later , when studying the dynamics of pointer states in an external potential .
as a first step , consider the standard deviation @xmath137 of the numerically obtained dimensionless solitonic solution @xmath138 of ( [ eq : dimensionless ] ) as a function the dimensionless parameter @xmath139 . as shown by the solid line in fig .  4 , the size @xmath137 increases linearly with @xmath140 .
this observation can be reproduced by a simplified model which has the practical advantage that it can be applied to more complex systems , such as 3d gases with a microscopically realistic localization rate @xmath141 @xcite .
the idea of the model goes as follows : the ideal gas environment consists of particles which collide with the system at a rate @xmath2 . at each collision , the ambient particles gain position information , such that the wave function gets spatially localized to a length scale @xmath142 determined by the localization rate @xmath143 , see ( [ eq : decaycoherences ] ) .
after the scattering event , the particle disperses freely , until it gets localized again by a subsequent collision .
the pointer width @xmath144 is then obtained by averaging the time - dependent width of the wave function over the waiting time distribution of a poisson process .
more specifically , we assume that the length scale @xmath145 is characterized by the free parameter @xmath146 using ( [ eq : localization ] ) and taking the momentum transfer distribution @xmath147 to be a centered gaussian with variance @xmath148 , we obtain @xmath149 with @xmath150 . the free dispersion after the collision yields the time dependent size @xmath151 finally , the average over the waiting time distribution @xmath152 gives @xmath153 where we use a linearization of @xmath154 in the second line . the dimensionless version of ( [ eq : size ] ) reads @xmath155 the dashed line in figure 4 shows that the form of this equation agrees with the numerical solution of ( [ eq : dimensionless ] ) ; the fit yields a value of @xmath156 for the parameter characterizing the localization length scale .
our next aim is to show that the solitonic solutions of ( [ eq : ngcd ] ) , which are interpreted below as the pointer states of collisional decoherence , form an overcomplete basis .
for this purpose , we first present a general method to construct a whole manifold of solutions of ( [ eq : ng ] ) given a specific one .
it relies on the symmetry properties of the corresponding master equation . since collisional decoherence exhibits galilean ( i.e. translation and boost ) invariance , it is then easy to show that the pointer states of this model form an overcomplete basis .
suppose there is a family of unitary operators @xmath157 , satisfying @xmath158 ,   \label{eq : condition2}\end{aligned}\ ] ] where we denote by @xmath159 the incoherent part of the master equation , @xmath160 / \left ( i \hbar \right ) + \mathcal{d } \left ( \rho \right)$ ] .
then , given a solution @xmath161 of the nonlinear equation ( [ eq : ng ] ) , also @xmath162 constitutes a solution of ( [ eq : ng ] ) .
this can be verified easily : @xmath163 \right ]   & = & \frac{1}{i    \hbar } \left ( \mathsf{h } \mathsf{u } \mathsf{p } \mathsf{u}^{\dag } -    \mathsf{u } \mathsf{p } \mathsf{u}^{\dag } \mathsf{h } \right ) + \mathsf{u }    \left [ \mathsf{p } , \left [ \mathsf{p } , \mathcal{d } \left ( \mathsf{p } \right )     \right ] \right ] \mathsf{u}^{\dag }   \nonumber\\    & = & { \frac{1}{i \hbar } ( \mathsf{h } \mathsf{u } \mathsf{p }    \mathsf{u}^{\dag } - \mathsf{u } \mathsf{p } \mathsf{u}^{\dag } \mathsf{h } +    \mathsf{u } \mathsf{h } \mathsf{p } \mathsf{u}^{\dag } - \mathsf{u } \mathsf{h }    \mathsf{p } \mathsf{u}^{\dag } + \mathsf{u } \mathsf{p } \mathsf{h }    \mathsf{u}^{\dag } } \nonumber\\    &   & - \mathsf{u } \mathsf{p } \mathsf{h } \mathsf{u}^{\dag } ) + \mathsf{u }    \left [ \mathsf{p } , \left [ \mathsf{p } , \mathcal{d } \left ( \mathsf{p } \right )     \right ] \right ] \mathsf{u}^{\dag }   \nonumber\\    & = & \frac{1}{i \hbar } \left [ \mathsf{h } , \mathsf{u } \right ] \mathsf{p } _ { }    \mathsf{u}^{\dag } + \mathsf{u } \left [ \mathsf{p } _ { } , \left [ \mathsf{p } _ { } ,    \frac{1}{i \hbar } \left [ \mathsf{h } , \mathsf{p } \right ] + \mathcal{d } \left (    \mathsf{p } _ { } \right ) \right ] \right ] \mathsf{u}^{\dag } \nonumber\\    &   & - \frac{1}{i \hbar } \mathsf{u } \mathsf{p } _ { } \left [ \mathsf{u}^{\dag } ,    \mathsf{h } \right ]   \nonumber\\    & = & \partial_t \left ( { \mbox{\rm $ \mathsf{u } \mathsf{p } \mathsf{u}^{\dag}$ } }    \right ) \
, ,   \label{eq : upu}\end{aligned}\ ] ] where we dropped the time argument for brevity . here
, the first equality makes use of ( [ eq : condition1 ] ) and the unitarity of @xmath164 , and the third equality is due to the relation @xmath165 = \left [ \mathsf{p } , \left [ \mathsf{p } , \left [ \mathsf{h } , \mathsf{p } \right ] \right ] \right]$ ] . in the last line ,
we use ( [ eq : ng ] ) and ( [ eq : condition2 ] ) .
let us now apply this to the galilean invariance of master equation ( [ eq : cd ] ) .
we will see that the phase space translations @xmath166 satisfy the symmetry conditions ( [ eq : condition1 ] ) and ( [ eq : condition2 ] ) provided the time dependence of @xmath167 and @xmath168 has the particular form @xmath169 the latter enact a phase space translation in accordance with the free shearing motion .
let us first verify condition ( [ eq : condition1 ] ) using @xmath170 .
@xmath171 in order to verify ( [ eq : condition2 ] ) , use the campbell - hausdorff formula to rewrite the translation operator as @xmath172 the time derivative thus yields @xmath173 \ , ,   \label{eq : timederivativeoft}\end{aligned}\ ] ] where the shearing transformation ( [ eq : time1 ] ) and ( [ eq : time2 ] ) is required in the second line .
this confirms ( [ eq : condition2 ] ) for @xmath174 .
we conclude that the nonlinear equation ( [ eq : ngcd ] ) exhibits a family of solitonic solutions @xmath175 , parameterized by the phase space coordinate @xmath176 . in order to verify that this family forms an overcomplete basis , let us consider a specific class of phase space representations . according to @xcite ,
any hilbert - schmidt operator @xmath22 can be represented as @xmath177 provided @xmath178 is a trace - class operator , i.e. @xmath179 . here
, @xmath180 denotes a phase space integral and @xmath181 is a function of the phase space coordinate @xmath182 .
choosing for @xmath22 the identity @xmath183 , and for @xmath178 the solitonic solution @xmath184 of ( [ eq : ngcd ] ) with vanishing position and momentum expectations , we obtain a resolution of the identity in terms of the solitons @xmath185 , @xmath186 this demonstrates that the pointer states of collisional decoherence form an overcomplete basis .
* * * * * *      so far , we have characterized the solitonic solutions of the nonlinear equation ( [ eq : ngcd ] ) which applies in the absence of an external force .
if an additional potential is present the corresponding nonlinear equation contains an additional term @xmath187 on the right hand side of ( [ eq : ngcd ] ) .
the numerical treatment shows that the solutions still converge to localized wave packets , which , however , change their shape and velocity in the course of the evolution .
we find that the center of these wave packets moves on the corresponding classical phase space trajectory provided the collision rate is sufficiently large .
we first summarize our numerical findings and then proceed with an analytic explanation .
figure 5 shows the position and momentum expectation values of the numerical solution of the nonlinear equation , in case of an anharmonic external potential of the form @xmath188 ( starting from an gaussian initial state ) .
the panel on the left hand side of fig .
5 was obtained in the limit of a vanishing collision rate @xmath189 ( i.e. @xmath190 ) , which turns ( [ eq : dimensionless ] ) into the schrdinger equation .
the solution therefore disperses , and the solid line shows a typical evolution of the phase space expectation values . the dashed line , on the other hand , gives the classical trajectory of the phase space point where the initial state is localized .
the result for a large collision rate @xmath2 ( or small @xmath191 ) is shown on the right hand side of fig .
, the initial state turns rapidly into a soliton whose expectation values move on the corresponding classical trajectory .
this illustrates that the temporal evolution turns from quantum to classical dynamics with increasing collision rate @xmath2 ( i.e. decreasing @xmath191 ) .
we made similar observations with various other potentials .     in order to explain the numerical observation ,
first consider a particle in a linear potential @xmath192 .
the corresponding nonlinear equation reads as @xmath193    \ , .
\label{eq : linear}\end{aligned}\ ] ] as discussed in sect .
4.1 , the field - free version of this equation @xmath194 exhibits uniformly moving solitonic solutions of the form @xmath195 \right ) \ , .
\label{eq : movingsoliton}\end{aligned}\ ] ] this implies that ( [ eq : linear ] ) has solitonic solutions of the form @xmath196 with @xmath197 and @xmath198 . in order to verify this statement
, we evaluate the left hand side of ( [ eq : linear ] ) with the ansatz given by ( [ eq : accelerated ] ) .
this yields @xmath199 with @xmath200 , @xmath201 and @xmath202 .
the expression can be further simplified by noting that the free soliton ( [ eq : movingsoliton ] ) is a solution of the field - free version of ( [ eq : linear ] ) , implying that @xmath203 \left ( x \right ) & = & if \partial_t    \chi \left ( t \right ) - v \left ( \partial_x f + if \partial_x \phi \right )    \nonumber\\    &   & + \frac{\hbar}{m } \left ( \partial_x f \partial_x \phi - \frac{i}{2 }    \partial_x^2 f + \frac{i}{2 } f \left ( \partial_x \phi \right)^2 +    \frac{1}{2 } f \partial_x^2 \phi \right ) \ , ,   \label{eq : freefieldcase}\end{aligned}\ ] ] with @xmath204 \left ( x \right)$ ] defined in ( [ eq : lambda ] ) . using ( [ eq : derivatives ] ) , ( [ eq : freefieldcase ] ) and the above definition of @xmath205 , one finds that @xmath206 \left ( x \right ) \,,\end{aligned}\ ] ] which confirms that @xmath106 evolves according to ( [ eq : linear ] ) .
we conclude that in a linear potential the pointer states have the same shape as in the field - free case and they are uniformly accelerated like a classical particle . for general potentials ,
this implies that the pointer states follow the corresponding classical motion , provided the spatial width of the solitons is sufficiently small , such that the linearization of the potential is justified over their spatial extension . since the size of the pointer states decreases with the collision rate ( see section 4.3 ) , the pointer states must exhibit classical dynamics in the limit of large collision rates .
thus far , we have calculated ` candidate ' pointer states as the solitonic solutions of ( [ eq : ngcd ] ) , and we have studied their properties and dynamics .
next , it will be shown that these ` candidates ' are genuine pointer states @xmath21 in the sense of ( 1 ) .
moreover , we find that the statistical weights @xmath207 of the pointer states are given by the overlap of the initial state @xmath15 with the initial pointer state @xmath208 , i.e. @xmath209 \ , .
\label{eq : initaloverlap}\end{aligned}\ ] ] in order to verify the above conjectures , let us now make use of the formalism of quantum trajectories @xcite to solve the master equation ( [ eq : cd ] ) .
more precisely , a specific quantum trajectory method , the _ orthogonal unraveling _ , is distinguished by the physics of pointer states because the deterministic part of the associated stochastic differential equation coincides with the nonlinear equation ( [ eq : ng ] ) .
we start with a general description of the quantum trajectory approach and the orthogonal unraveling in section 5.1 .
the latter will be applied to collisional decoherence in section 5.2 .
this will allow us to evaluate the statistical weights of the pointer states in section 5.3 .        in the quantum trajectory approach ,
the wave function corresponding to a pure initial state @xmath210 is propagated stochastically to generate pure state trajectories @xmath211 , whose ensemble average recovers the solution of the master equation ( [ eq : lindblad ] ) , i.e. @xmath212 such a stochastic process is called an unraveling of the master equation .
a common unraveling is provided by the quantum monte carlo method @xcite which is based on a piecewise deterministic process . here , one realization of a trajectory consists of smooth deterministic pieces generated by an effective ( non - hermitian ) hamiltonian , in our case @xmath213 which are interrupted by random jumps .
the jumps occur with the rate @xmath214 an expectation value with respect to @xmath215 , and they are effected by the operators @xmath216 the quantum monte carlo method is not the only stochastic process satisfying ( [ eq : average ] ) .
in fact , there exists an infinite set of these stochastic processes , because the convex decomposition of the density matrix on the left hand side of ( [ eq : average ] ) is not unique .
for instance , other unravelings can be obtained from the quantum monte carlo method , by noting that the generator @xmath42 does not uniquely fix the lindblad operators @xmath217 and the hamiltonian @xmath218 @xcite .
this is due to the fact that the master equation ( [ eq : lindblad ] ) is invariant under certain transformations of the lindblad operators , such as the addition of a complex multiple @xmath219 of the identity @xmath220 , @xmath221 in the latter case , also the hamiltonian must be transformed as @xmath222 in order to assure the invariance of the master equation .      to obtain a different , but
again piecewise deterministic unraveling , we now make the choice @xmath223 .
it then follows that the deterministic pieces of a sample path are given by the solution of the nonlinear equation ( [ eq : nonlinear ] ) ( or equivalently the corresponding projector equation ( [ eq : ng ] ) ) .
the jumps occur with the rate @xmath224 and are caused by the nonlinear operators @xmath225 as a distinctive feature , the states @xmath226 into which the system may jump are orthogonal to the original state @xmath227 , thus justifying its naming .
( the states @xmath228 are not necessarily mutually orthogonal , though . ) to our knowledge , this unraveling was first noted by rigo and gisin @xcite , although , it has not been studied numerically so far .
a related unraveling , which is also referred to as the ` orthogonal unraveling ' , was introduced by disi @xcite . here , the deterministic pieces of the evolution are as well generated by ( [ eq : nonlinear ] ) .
however , the states @xmath228 into which the system may jump are obtained differently , as the eigenvectors of the hermitian operator @xmath229 where @xmath230 . as a consequence ,
these states are also _ mutually _ orthogonal ( in finite dimensional systems ) . since the orthogonal unraveling of @xcite requires the diagonalization of the operator ( [ eq : rateoperator ] )
, it is much more involved than the one defined by ( [ eq : nonlinear ] ) , ( [ eq : rate2 ] ) and ( [ eq : jump ] ) , which is why we will use the latter in the following .      as mentioned above , all unravelings are equivalent in the same sense as the different convex decompositions of the density matrix @xmath231 .
note , however , that a preferred set of pure states  the pointer basis  may be singled out through the environmental coupling . in that case , those unravelings are distinguished which generate for any initial state an ensemble of these state independent projectors .
an unraveling will do this job if ( a ) its deterministic part exhibits stable fixed points or solitons @xmath232 , which ( b ) are characterized by a vanishing jump rate , @xmath233 . in that case , the sample paths of the process will end up in one of the states @xmath21 , by all means .
hence , the ensemble mean is of the form ( [ eq : decoherence ] ) , such that the fixed points or solitons @xmath21 can be identified with the pointer states of the system .
for the case of collisional decoherence the orthogonal unraveling fulfills the aforesaid conditions .
( a ) its deterministic part is given by ( [ eq : ngcd ] ) which exhibits the solitonic solutions @xmath234 shown in fig .  2 .
we shall see explicitly in section 5.2.1 that these states are attractive fixed points .
( b ) we will show in section 5.2.2 that the jump rate ( [ eq : rate2 ] ) vanishes for the solitons @xmath234 , which finally demonstrates that our ` candidate ' pointer states @xmath235 are genuine pointer states @xmath21 in the sense of ( 1 ) .
moreover , this shows that the orthogonal unraveling is a very efficient numerical scheme for the long time solution of the master equation ( [ eq : cd ] ) , since the state is no longer affected by the stochastic part , once it has turned into the soliton , and the trajectory is therefore more easy to integrate .
we note that the orthogonal unraveling is not the only stochastic process which generates the ensemble of pointer states .
in fact , there is a _ diffusive _ unraveling @xcite , which also involves ( [ eq : ng ] ) as its deterministic part .
it was applied in @xcite to investigate pointer states in a linear model .
we now apply the described orthogonal unraveling to collisional decoherence ( [ eq : cd ] ) , first evaluating the deterministic part ( [ eq : nonlinear ] ) of the stochastic process in sect .
5.2.1 and then the stochastic one ( [ eq : rate2 ] ) and ( [ eq : jump ] ) in sect .  5.2.2 .      applying ( [ eq : nonlinear ] ) to the case of collisional decoherence yields the soliton equation ( [ eq : ngcd ] ) discussed in sect .  3.2 .
we will now further simplify this equation , by considering initial states @xmath236 which are superpositions of non - overlapping wave functions @xmath237 , @xmath238 the latter are assumed to be localized in the sense that @xmath239 where @xmath240 and @xmath241 denote the variances of the distributions @xmath242 and @xmath243 , respectively . under this assumption , which will be justified at the end of this section
, one can extract a system of evolution equations for the time evolution of the coefficients in ( [ eq : initialstate ] ) , @xmath244 here , the matrix @xmath245 is obtained from the localization rate ( [ eq : localization ] ) , where the @xmath246 denote the mean positions of the constituent wave functions @xmath247 .
the latter evolve according to @xmath248 \left ( x , t \right ) \nonumber\\    &   & + \phi_{^ { } i } ( x , t ) \sum^n_{j = 1 , j \neq i } { } \left| c_j    \left ( t \right ) \right|^2   \tilde{\gamma}_{ij } \left ( x , t \right ) \ , ,     \label{eq : basis}\end{aligned}\ ] ] where @xmath249 is defined in ( [ eq : lambda ] ) and @xmath250 is a rate of the order of @xmath2 , @xmath251 let us now verify that @xmath252 , with @xmath253 and @xmath254 solutions of ( [ eq : coefficients ] ) and ( [ eq : basis ] ) , evolves according to ( [ eq : ngcd ] ) .
first , we note that the assumption ( [ eq : a2a ] ) leads to the approximation @xmath255 for all @xmath65 contributing appreciably to integrals weighted with the momentum transfer distribution @xmath74 .
this , in turn , implies @xmath256 hence , we find that @xmath257 , and one thus obtains @xmath258 \left ( x , t \right ) & = & -    \left ( \left| \phi_i \right|^2 \ast f \right ) \left ( x , t \right ) \ , .
\label{eq : lambdaa2}\end{aligned}\ ] ] now , consider the time derivative of @xmath259 , which gives @xmath260 where we dropped the arguments for brevity .
in ( [ eq : beweis1 ] ) , we used ( [ eq : coefficients])-([eq : basisrate ] ) , ( [ eq : lambdaa2 ] ) , the fact that @xmath261 , and the normalization condition @xmath262 . now , we replace @xmath263 in ( [ eq : beweis1 ] ) by the right hand side of ( [ eq : fij ] ) , and use ( [ eq : a1 ] ) , which yields @xmath264 finally , by using ( [ eq : localization ] ) , we obtain @xmath265 which confirms that @xmath259 evolves according to ( [ eq : ngcd ] ) .
let us now discuss the solution of eq .
( [ eq : coefficients ] ) for the time evolution of the coefficients .
we first consider situations where the wave packets @xmath266 are sufficiently far apart such that the localization rate is saturated , i.e. @xmath267 under this assumption , ( [ eq : coefficients ] ) reduces to the equation @xmath268 which was already studied in @xcite in the context of a discrete model for quantum measurement .
it is shown there that all stable fixed points of ( [ eq : coefficientssaturated ] ) have the form @xmath269 , and that the particular fixed point @xmath270 , with @xmath271 is approached monotonically , i.e. the component with the largest initial weight wins .
this behavior is visualized in fig .
6 which was obtained by solving ( [ eq : coefficientssaturated ] ) numerically for the case @xmath272 . here ,
the @xmath85- and the @xmath273-axis indicate the weights @xmath274 and @xmath275 , respectively .
the plot on the left hand side shows various trajectories , illustrating in particular the fixed points .
the plot on the right displays the regions of attraction of the stable fixed points @xmath276 , in agreement with the criterion ( [ eq : argmax ] ) .
for instance , area @xmath277 highlights the region of attraction of the fixed point @xmath278 .
figure 7 , on the other hand , depicts a scenario where the wave packets @xmath266 are close together such that the localization rate is unsaturated , i.e. @xmath279 . here
, we choose @xmath272 wave packets with non - equidistant position expectations , @xmath280 .
similarly to the saturated case , we observe that all stable fixed points of ( [ eq : coefficients ] ) have the form @xmath281 .
however , the regions of attraction are deformed such that criterion ( [ eq : argmax ] ) is no longer valid , and the fixed points are not necessarily approached monotonically .   to see that @xmath281 are stable fixed points of ( [ eq : coefficients ] ) , assume that the coefficients are close by , i.e. @xmath282 with @xmath283 .
it follows from ( [ eq : coefficients ] ) that @xmath284 and hence , @xmath285 .
the knowledge of the fixed points of the coefficients allows one to discuss the asymptotic evolution of the initial state shown in eq .
( [ eq : initialstate ] ) .
since the coefficients @xmath286 with @xmath287 tend to zero asymptotically , it follows that @xmath288 for a specific @xmath289 ( which is given by ( [ eq : argmax ] ) in the saturated case ) .
the asymptotic behavior of the wave packets @xmath290 can , in turn , be predicted from eq .
( [ eq : basis ] ) .
since the @xmath291 vanish for large times , the coupling term given by the last summand in * * ( [ eq : basis ] ) vanishes as well , implying that the time evolution ( [ eq : basis ] ) of @xmath290 is asymptotically equal to the soliton equation ( [ eq : ngcd ] ) .
therefore , in the absence of stochastic jumps , the state @xmath292 evolves into that solitonic solution @xmath293 of ( [ eq : ngcd ] ) which is associated to the initial wave packet @xmath294 .
it should be mentioned that eqs .
( [ eq : coefficients ] ) and * * ( [ eq : basis ] ) for the coefficients @xmath98 and the constituent wave packets @xmath97 are not completely decoupled , since ( [ eq : coefficients ] ) depends on the matrix @xmath295 which contains the position expectations @xmath296 of the wave packets @xmath97 .
however , the position expectation follows the classical trajectory for sufficiently large @xmath191 s , such that ( [ eq : coefficients ] ) can be solved without knowing the solution of * * ( [ eq : basis ] ) .
let us now discuss the validity of the assumption of small position variance ( [ eq : a2a ] ) , and the ensuing approximation ( [ eq : a2b ] ) .
it can be justified by our observation in section  4.3 , that the dimensionless pointer width @xmath297 is a function of the parameter @xmath298 only , @xmath299 thus , for all @xmath300 we find that the position variance @xmath301 of a pointer state is one order of magnitude smaller than the reciprocal width of the momentum transfer distribution @xmath74 , @xmath302 the above relation for the width of the pointer state is sufficient to justify the approximation @xmath303 , as we checked numerically , by using the solitonic solution of ( 6 ) . the relative error is less than @xmath304 for @xmath305 $ ] and @xmath306 .      upon inserting the lindblad operator @xmath307 into ( [ eq : jump ] ) we find that the jump operator takes the form @xmath308 with normalization @xmath309 .
we again consider states of the form ( [ eq : initialstate ] ) which are superpositions of non - overlapping ( [ eq : a1 ] ) and localized ( [ eq : a2a ] ) wave packets @xmath310 . under this assumption
, one can evaluate the expectation value in ( [ eq : jumpoperator ] ) @xmath311 such that the state @xmath312 into which the system may jump takes the form @xmath313    later we will choose the initial wave packets @xmath266 to be solitons @xmath314 .
let us therefore assume that the @xmath97 s form a basis , such that @xmath315 can be represented as @xmath316 .
then the transformed coefficients @xmath317 can be evaluated by the overlap @xmath318 similarly , one can evaluate the rate ( [ eq : rate2 ] ) associated to the jump operator of collisional decoherence , @xmath319 the above approximation ( [ eq : expectation ] ) further simplifies this expression , @xmath320 one observes that this rate vanishes for the stable fixed points @xmath321 , indicating that the quantum trajectories of the orthogonal unraveling evolve into the pointer states , the solitonic solutions of ( [ eq : ngcd ] ) .
moreover , we note that the original stochastic process ( [ eq : nonlinear ] ) , ( [ eq : rate2 ] ) , ( [ eq : jump ] )  which is defined in the infinite dimensional hilbert space of the system  has been reduced to a stochastic process in @xmath322 , demonstrating the efficiency of the orthogonal unraveling . however , due to the finite pointer width ( [ eq : sizedimensionless ] ) the _ exact _ expression for the jump rate ( [ eq : exactrate ] ) does not vanish identically , although it is very small compared to @xmath2 . for instance , the numerically obtained soliton displays a strongly suppressed total jump rate @xmath323 of @xmath324 for @xmath325 , while the superposition state decays with the rate @xmath326 .
this implies that the solitons are not perfect pure state solutions of the master equation ( [ eq : cd ] ) , though the loss of purity is small .
the previous section showed that the orthogonal unraveling of an initial superposition state subject to collisional decoherence can be reduced to a stochastic process with respect to the corresponding coefficients .
in particular , this applies to the case where the initial state is a superposition of pointer states , @xmath327 thus we can now verify , by using the discrete process defined by eqs .
( [ eq : coefficients ] ) , ( [ eq : process2 ] ) and ( [ eq : process3 ] ) , that after decoherence the statistical weights of the pointer states are given by the overlap of the initial state with the initial pointer states .
more specifically , this demonstrates that the initial state @xmath328 evolves into the mixture @xmath329 where the statistical weights are given by the overlap @xmath330 we first present an analytic proof of the above for @xmath331 . the general case , @xmath332 ,
is then treated numerically in the following section .
we consider the expectation value for the coefficients after a jump , that is @xmath333 upon inserting ( [ eq : process2 ] ) one obtains @xmath334 , with the normalization constant @xmath335 , see ( [ eq : jumpoperator ] ) . using @xmath336 , we find @xmath337 which implies @xmath338 this shows that after an average jump the moduli of the coefficients are simply interchanged .
this property ( which does not hold for @xmath339 ) makes the stochastic process analytically tractable , not least because the dynamics is independent of the phases of the coefficients .
since the deterministic part ( [ eq : coefficients ] ) of the evolution is monotonic for @xmath340 , a trajectory starting from @xmath341 will end up in the state @xmath342 if and only if an odd number of jumps occurs in the process .
this is demonstrated in fig .
crucially , the jump rate @xmath343 , @xmath344 is unaffected by the jump ( [ eq : averagejump ] ) at all times , since it is invariant under interchanging the coefficients .
hence , the time dependence of the jump rate is identical for all trajectories , which , in turn , means that the number of jumps follows an inhomogeneous poisson process .
therefore , the probability for an odd number of jumps , which is equal to the statistical weight @xmath345 of the pointer state @xmath346 , is given by @xmath347 with @xmath348 the integrated jump rate .
the latter can easily be evaluated by noting that ( [ eq : coefficients ] ) can be written for @xmath340 as @xmath349 upon inserting this result into ( [ eq : totalrate ] ) we obtain the integrated jump rate : @xmath350 noting ( [ eq : probodd ] ) we thus find the probability for an odd number of jumps @xmath351 this finally confirms that the statistical weights of the pointer states are indeed given by the expected overlap ( [ eq : distribution ] ) .
the stochastic process is much more involved if the initial superposition consists of more than two pointer states .
our numerical implementation of the stochastic process defined by ( [ eq : coefficients ] ) , ( [ eq : process2 ] ) and ( [ eq : process3 ] ) is based on a metropolis - hastings algorithm to draw the momentum transfer @xmath65 in accordance with the rate ( [ eq : process3 ] ) , with @xmath147 a gaussian . each of the generated trajectories ends asymptotically in one of the fixed points corresponding to a pointer state , and we thus obtain a numerical estimate of the statistical weights by means of the relative frequencies @xmath352 , of the asymptotic states . to confirm the expected probability distribution @xmath353 , we evaluate the relative entropy @xmath354 between these two distributions .
figure 9 shows the result for a random initial state with @xmath355 as a function of the number of trajectories @xmath356 , indicating convergence to zero .
i**n addition , we found for @xmath357 random initial states , with random @xmath358 , based on @xmath359 trajectories that the relative entropy was always less than @xmath360 .
this holds both for cases where the initial wave packets @xmath361 are far apart such that the localization rate is saturated , @xmath362 , and for situations where the wave packets are located close together such that @xmath363 .
we conclude that the asymptotic trajectories are indeed distributed according to the expected overlap ( [ eq : distribution ] ) .   as an alternative confirmation of the statistical weights
, we performed a @xmath364-test .
similar to the treatment above , @xmath357 random initial states @xmath365 , with random @xmath366 , were drawn by the simplex picking method .
for each random state , @xmath367 trajectories were generated , each of which ends asymptotically in one of the pointer states . using the observed relative frequencies @xmath368 , of the pointer states , we evaluate @xmath369 for each random state . in order to verify that the pointer states are distributed according to @xmath370
, the set @xmath371 must be shown to be sampled from a @xmath372-distribution with @xmath373 degrees of freedom . comparing the set @xmath371 with the @xmath374-quantiles * *
quantile is the value such that @xmath375 of the samples lie below @xmath376 . ]
( denoted by @xmath377 ) of the corresponding @xmath364-distribution , a typical run shows ten cases where @xmath378 , one case where @xmath379 , but not a single case where @xmath380 , as one expects if the @xmath371 are @xmath364-distributed .
like above , this confirms statistically that the asymptotic trajectories are distributed according to the expected overlap ( [ eq : distribution ] ) .
in this article , we related the nonlinear pure state equation discussed in @xcite to a specific orthogonal unraveling of the collisional decoherence master equation .
this puts into evidence that the dynamics of a particle in an ideal gas environment can be represented by an ensemble of pure state trajectories which evolve into spatially localized pointer states . for sufficiently strong collisions with the background gas ,
these solitonic wave packets move according to the classical equations of motion , thus explaining the emergence of classical dynamics within the quantum framework .
once the pointer state is reached by an individual quantum trajectory , the latter is no longer affected by the stochastic part of the unraveling , such that the integration of the trajectory is reduced to the solution of the classical equations of motion .
this suggests that the orthogonal unraveling is an efficient algorithm for the long time solution of master equations which exhibit a pointer basis . on the other hand ,
also the short time solution turns out to be efficient , since the orthogonal unraveling can be reduced ( under appropriate assumptions ) from an infinite dimensional unraveling to a stochastic process in @xmath322 .
future studies might consider the emergence and dynamics of pointer states in dissipative quantum systems .
we note that the present work relies on the model of pure collisional decoherence which does not describe long time effects such as dissipation or thermalization .
it would certainly be worth to determine the pointer states of a more involved model such as the quantum linear boltzmann equation presented in @xcite . for large mass ratios between the test particle and the gas
, one expects that the pointer states then evolve according to a langevin equation , thus explaining the emergence of classical brownian motion within the quantum framework .
k.  hornberger , _ introduction to decoherence theory _ , in : _ entanglement and decoherence .
foundations and modern trends _ , ed . by a.
buchleitner , c.  viviescas , and m.  tiersch , lect .
notes phys .
* 768 * , 221276 , springer , berlin , 2009 .