ever since the foundation of modern cosmology based on general relativity , there have been attempts to determine the cosmological parameters of our universe .
gravitational lensing has been recognized as one of several such tools .
for example , the results of radio surveys for gravitational lenses have recently been used to place significant constraints on cosmological models ( falco et al .
it has also been realised that observations of supernovae at high redshift can be used for this purpose ( goobar & perlmutter @xcite ) , in particular for determining the value of the cosmological constant .
in fact , several collaborations with this in mind are in progress , and the first sets of data show an intriguing hint of a non - vanishing cosmological constant in a universe consistent with having a flat geometry ( see , e.g. , riess et al .
@xcite ; perlmutter et al .
although the two groups which have published results mutually agree on the best - fit parameters , it is important to note that the effects of geometry are small ( on the order of half a magnitude ) , and the need to go to even higher redshift to get larger effects is obvious .
when observing such distant sources , at redshift greater than unity , it is necessary to estimate the effects of lensing due to inhomogeneities in the matter distribution .
this will be of growing importance when supernovae at very high redshifts become accessible , e.g. , with ngst , the next generation space telescope ( miralda - escud & rees @xcite ; dahln & fransson @xcite ) or with the dedicated supernova search satellite , snap ( perlmutter et al .
, private communication ) .
the literature on gravitational lensing is quite rich .
much of the history and developments up to the early 1990s can be found in the excellent textbook by schneider , ehlers & falco ( @xcite ) . a new method for examining lensing effects
has recently been proposed by holz & wald ( @xcite ; hw ) , one which has the virtue of lending itself easily to numerical calculations . in that work
it is also shown that given a large and deep enough sample of standard candle supernovae , one may in principle distinguish between different matter distributions in the galactic halos responsible for the lensing . in particular the case of point masses gives a quite different distribution of magnifications ( or de - magnifications ) compared to the case of singular isothermal spheres ( see also metcalf & silk @xcite ) .
the method of hw can be summarized as follows : first , a friedmann - lematre ( fl ) background geometry is selected .
inhomogeneities are accounted for by specifying a matter distribution in cells which have an average energy density equal to that of the underlying fl model .
following the history of observed light from a distant source , a light ray is traced backwards to the desired redshift by being sent through a series of cells , each time with a randomly selected impact parameter with respect to the matter distribution in the cell . between each cell ,
the fl background is used to update the scale factor and the expansion rate . by using monte carlo techniques to trace a large number of such light rays , statistics for the apparent luminosity of an ensemble of sources at a given redshift
is obtained .
the purpose of this paper is to generalize the method of hw in a number of ways and to investigate some issues related to the different lensing signatures in the luminosity distribution of supernovae caused by different mass distributions of the intervening lens population .
the outline is as follows . in sect .
[ sec : cosmodel ] it is demonstrated how the hw method can be generalized from matter in dust form to perfect fluids with pressure .
[ sec : area ] is concerned with the interpretation of results as probability distributions .
[ sec : devi ] presents the geodesic deviation and in sect .
[ sec : hwsummary ] , we summarize and discuss some of the conceptual issues in the analysis of hw . in sect .
[ sec : nfw ] , we include the possibility of using the more realistic halo matter distribution proposed by navarro et al .
( @xcite ; nfw ) based on their n - body simulation results , and in sect .
[ sec : massdist ] we discuss the mass distribution of lensing galaxies . by using realistic density profiles and mass distributions of dark matter halos ,
we can obtain high accuracy results on gravitational lensing without the need of using the extensive full data sets of n - body simulations , as in some other methods ( wambsganss et al .
@xcite ; premadi et al .
@xcite ; jain et al .
one strength of our ray - tracing monte carlo method is further that it can be continuously refined as more observational information , e.g. , on galaxy distributions is obtained .
also , effects such as absorption by dust and other possibly @xmath0-dependent effects can be straightforwardly added to the algorithm . in sect .
[ sec : results ] , some results from the generalizations of the method of hw are presented as luminosity distributions . in sect .
[ sec : limit ] , we investigate some limitations of the numerical implementation of the method and in sect .
[ sec : mi ] , we discuss multiple images .
[ sec : cc ] investigates different consistency checks , sect .
[ sec : formula ] contains convenient analytical fits to our numerical results and the paper is concluded with a discussion in sect .
[ sec : disc ] .
the starting point of hw is a newtonianly perturbed fl universe .
the line element of robertson - walker form is then @xmath1 where @xmath2 is the scale factor , and @xmath3 is a constant - curvature three - geometry , for which @xmath4 indicates the sign of the curvature .
( we use geometrized units such that @xmath5 ) .
the perturbation @xmath6 is assumed to satisfy a number of properties that will be employed when determining dominant terms .
firstly , @xmath7 , reflecting that @xmath6 indeed is a perturbation .
secondly , inhomogeneities are assumed to be significant , so that @xmath8 is small compared to spatial derivatives .
thirdly , second - order spatial derivatives will be assumed to dominate over products of first - order ones ( e.g. , @xmath9 ) . in this paper , we have generalized the treatment of hw so that the matter content of the model can be taken to be a number of perfect fluids : @xmath10 where @xmath11 is the stress - energy tensor , @xmath12 is the energy density , @xmath13 the pressure of fluid @xmath14 , @xmath15 is the ( normalized ) fluid four - velocity , and @xmath16 is the metric corresponding to the line element ( [ eq : line - element ] ) .
cosmologically reasonable equations of state are linear barotropic ones : @xmath17 . in this paper ,
the @xmath18 s are constants with @xmath19 .
however , it is straightforward to allow for the equations of state to evolve with redshift , i.e. , @xmath20 .
hw considered dust models with a cosmological constant , which corresponds to a two - fluid model , the first fluid having @xmath21 , and the second fluid homogeneously distributed with @xmath22 and @xmath23 .
the field equations , corresponding to eqs .
( hw7,8 ) , become + _ the friedmann equation : _
@xmath24 _ the raychaudhuri equation : _
@xmath25 where the definitions of the hubble scalar @xmath26 and deceleration parameter .
] @xmath27 have been made explicit , and where @xmath28 denotes the spatial projection of contracted covariant derivatives .
the pressure does not enter the friedmann equation , i.e. , the hubble scalar is entirely given by the energy densities and the curvature .
the deceleration , on the other hand , is affected by the presence of pressure . comparing the above field equations with the field equations for the corresponding fl background ( with averaged densities and pressures @xmath29 , @xmath30 ) it can be shown [ eq .  ( hw11 ) ] that @xmath6 satisfies a poisson equation @xmath31 , where @xmath32 , while @xmath33 vanishes .
thus , the pressures result in homogeneous contributions , and the specification of an equation of state only makes sense in terms of the averaged energy density : @xmath34 .    from the field equations for the background model
follows + _ the energy conservation equation : _
@xmath35 if the fluids are assumed to couple only to gravity , they will be separately conserved , and each component in the sum of eq .
( [ eq : ec ] ) vanishes .
note that it is only for dust ( @xmath36 ) that the mass within a volume co - moving with the hubble flow is constant ( @xmath37 ) . for other types of matter
, the mass varies with time because of the work done by the pressure .
also note that @xmath38 results in a constant energy density .
thus , a homogeneously distributed fluid with @xmath38 is equivalent to a cosmological constant , or in terms of quantum field models , to vacuum energy .
it is convenient to introduce density parameters @xmath39 which depend on time and therefore on redshift . assuming the above discussed equations of state the time dependence is governed by the following set of conservation equations ( see , e.g. , goliath & ellis @xcite ) @xmath40\omega_ih .\ ] ]    hw proceed by performing a number of coordinate transformations to co - moving coordinates on isotropic form , eq .
the inclusion of pressure results in a modification of this line element : @xmath41    dt^2 \nonumber \\    & + & \left[1 - 2\phi\right](dx^2+dy^2+dz^2 ) , \end{aligned}\ ] ] where @xmath42 , and the effective potential @xmath43 satisfies the poisson equation @xmath44    the lookback time of an event at redshift @xmath45 is the time difference ( with respect to the coordinate @xmath46 ) between the event and the present . by using @xmath47 , we get a relation between the lookback time and the redshift for a fl universe : @xmath48 to obtain @xmath49 , the friedmann equation , eq .
( [ eq : fe ] ) , is expressed in terms of redshift and the present densities @xmath50 .
to do so , first note that the energy conservation equation , eq .
( [ eq : ec ] ) , together with the equation of state gives @xmath51 the hubble scalar can then be expressed @xmath52 where we introduced the curvature density parameter @xmath53 the luminosity distance to the event at redshift @xmath45 is given by a similar expression ( see , e.g. , bergstrm & goobar @xcite ) @xmath54 note that , when @xmath55 , the present scale factor can be expressed [ using ( [ eq : ook ] ) ] as @xmath56 .
consider a narrow ( strictly speaking ; infinitesimal ) ray bundle connecting a source and an observer as depicted in fig .
[ fig : focusing ] .
the ray bundle can be focused either at the source or the observer .
the relevant quantity when determining luminosity distances is the ratio of the area of a beam and the solid angle it subtends .
of course , the smaller the distance , the brighter the source .    in the method of hw ,
a beam focused at the observer with infinitesimal @xmath57 is integrated backwards using the geodesic deviation equation , and the resulting area at the source sphere is considered . in the lower panel of fig .
[ fig : focusing ] , the beam has undergone additional focusing .
a finite aperture at the observer will detect light from a larger @xmath57 of an isotropically emitting source , and a finite source will subtend a larger solid angle on the sky . since surface brightness is conserved , the source will be magnified .
the magnification , @xmath58 , is defined as the ratio of the _ observed _ flux and the flux of the corresponding image in absence of any additional focusing .
@xmath59 where @xmath60 denotes the solid angle subtended by the source in absence of lenses .
hence the magnification will be proportional to the inverse of the area of the beam . in absence of any lenses
, the ray bundle will have an area corresponding to the so called _ empty - beam _ area ( see sect .
[ sec : devi ] ) , and we will therefore define the magnification as @xmath61 for further reference , see schneider , ehlers & falco ( @xcite ) .    as noted in hw ,
some care is needed when interpreting results as probability distributions .
this is due to the fact that individual ray bundles do not correspond to random source positions .
this is illustrated in the left panel of fig .
[ fig : focusing3 ] where we consider the situation of a ray bundle originating from an observer at @xmath62 at the center of a source sphere at some fixed redshift .    from the right panel of fig .
[ fig : focusing3 ] we see that focused beams sample a smaller fraction of the source sphere .
since the fraction of the source sphere being sampled by a beam is proportional to the area of the beam , i.e. , inversely proportional to the magnification , magnified sources will be overrepresented if not compensated for . in hw
, this is done by convoluting the obtained probability distribution , @xmath63 , with a factor proportional to @xmath64 , thus obtaining a new distribution @xmath65 .
the method we use is very similar in that we designate a probability , proportional to the area of the ray bundle , to keep every event .
the advantage of this rejection method is that each selected event represents a statistically `` real '' event .
further discussions of the statistical weight given to random line - of - sights can be found in ehlers & schneider ( @xcite ) .
the geodesic deviation equation for a light ray is ( see , e.g. , wald @xcite ) @xmath66 where @xmath67 is the deviation vector and @xmath68 is a null vector tangent to the light ray .
the affine parameter @xmath69 is a parametrization along the light ray given by @xmath70 , where @xmath71 is the frequency of the photons in the beam .
it is instructive to calculate the deviation in an exact ( homogeneous ) fl universe , analogous to eq .
( hw21 ) : @xmath72 thus , the pressure will cause additional focusing of the light beam . note that a cosmological constant , which corresponds to @xmath38 , does not affect the deviation .    to derive the deviation experienced when propagating through a cell in the inhomogeneous model , hw start by writing the deviation equation in terms of a @xmath73 matrix @xmath74 , such that @xmath75 where @xmath76 ( hw28 ) .
the properties of the geodesic deviation are contained in @xmath74 .
for example , the area ( with orientation ) of the beam is given by @xmath77 .
in the case of the homogeneous background friedmann - lematre cosmology , we denote the corresponding reference area @xmath78 , which will of course depend on the redshift , but also on the cosmological parameters of the fl background model . also , the area of a light beam that propagates without encountering any matter will be denoted @xmath79 .
sometimes @xmath78 and @xmath79 are referred to as the _ filled - beam _ and _ empty - beam _ areas , respectively .    the deviation equation , eq .
( [ eq : devia ] ) , can now be written as a set of difference equations ( hw35,36 ) @xmath80 with @xmath81 where the integration is along a straight line through the cell .
next , @xmath82 is calculated from the line element eq .
( [ eq : lineco ] ) @xmath83        \!\eta_{\alpha\beta}\!\right\ } , \ ] ] where @xmath84 is the minkowski metric associated with the coordinates @xmath85 .
thus , with non - vanishing pressure , there will be an additional contribution in @xmath86 that takes the form @xmath87 with @xmath88 where @xmath89 is an impact parameter , and @xmath90 and @xmath91 are the radius and volume of the cell , respectively . assuming linear barotropic equations of state , @xmath92 , eq .
( [ eq : jp ] ) becomes @xmath93 where @xmath94 is the mass of matter component @xmath14 in the cell .
let us now examine the above result for the filled - beam case .
( hw39 ) gives the quantity @xmath86 for matter uniformly distributed in a ball of radius @xmath95 at the center of a cell of radius @xmath90 .
setting @xmath96 in this equation , and adding our derived pressure contribution @xmath97 to this yields @xmath98 from which it is seen that a cosmological constant ( @xmath99 ) does not contribute , in consistency with eq .
( [ eq : fldev ] ) .    as an interesting example , consider a flat cosmological model in which ordinary matter ( @xmath100 ) constitutes only 10 % of the total matter content , i.e. , @xmath101 , while the rest is made up of some homogeneously distributed `` x - matter '' with @xmath102 .
the ordinary matter is assumed to be in the form of point masses ( we will return later to more realistic mass distributions ) .
for this cosmology , garnavich et al .
( @xcite ) have put an upper limit on the parameter of the equation of state for an x - matter component using type ia supernova data .
their limit is @xmath103 at 95 % confidence level .
similarly , the supernova cosmology project ( perlmutter et al .
@xcite ) has deduced a limit of @xmath104 ( 95 % cl ) for any value of @xmath105 .    in fig .
[ fig : xmatter ] , we compare the cases of @xmath106 , corresponding to a cosmological constant , @xmath107 , which may either arise from `` quintessence '' ( i.e. , a slowly evolving scalar field ) ( wang et al .
@xcite ) or from a domain wall network ( battye et al .
@xcite ) . for reference
, we also show the result for @xmath108 , corresponding to dust ( i.e. , equivalent to an einstein  de sitter model with @xmath109 , of which 90 % is homogeneously distributed ) .
the choice of @xmath110 is made to allow for a direct comparison of the curve corresponding to @xmath106 in fig .
[ fig : xmatter ] with fig .  7 in hw ,
with which it agrees well .
note that the normalization differs between the three cases since @xmath78 depends on the cosmology . nevertheless , from the slope of the curves
, one can conclude that lensing effects become more prominent for large negative @xmath111 .
this is _ not _ what one would expect from studying the geodesic deviation equation , eq .
( [ eq : fldev ] ) , since negative pressure causes less focusing of a beam .
however , for negative @xmath111 , distances will get larger for a fixed redshift , which means a larger number of lenses between the source and the observer ( see , e.g. , turner @xcite ) .
the method of holz & wald ( @xcite ) which we adopt has a number of properties , which we summarize and comment upon here :    * * local nature of lensing * + if @xmath6 is small within distances much smaller than the hubble radius @xmath112 , then the curvature ( and thereby the lensing effects ) is determined locally within @xmath113 .
likewise , if we assume that there is no strong correlation on scales greater than a clustering scale @xmath114 , it follows that the curvature is determined locally within @xmath114 . * * relevant scales * + the question now is what relevant clustering scale to use .
hw find it convenient to employ the typical separation between galaxies , and argue that galaxy clustering should have a negligible effect on magnification , as shown in some of their simulations .
they also argue that the smallest relevant scale for matter clumping is of the order @xmath115 , since objects smaller than that would not give a coherent lensing effect over the angular size of most of the interesting sources . *
* generality * + for the case of point masses , hw demonstrate that their results are fairly insensitive to the size of individual point masses
. furthermore , their results are rather insensitive to clustering .
this is tested by cutting out a tube in a uniform matter distribution and distributing `` stars '' within this tube . *
* largest effects * + based upon their results , hw conjecture that randomly distributed point masses should cause the largest lensing effects of all possible mass distributions for a given average mass density . note that the point - mass case corresponds to the limit in which all matter clumps lie well within their einstein radii [ see eq .
( [ eq : re ] ) ] .
* comment on relevant scales : * + the einstein radius is a characteristic length scale in the lens plane given by @xmath116 where @xmath117 is the schwarzschild radius of the lens and @xmath118 , @xmath119 and @xmath120 are angular diameter distances between lens and source , observer and lens , and observer and source , respectively .
the corresponding length scale in the source plane , @xmath121 , is @xmath122    for lensing to be important , the source size , @xmath123 , has to be smaller than the length scale @xmath121 introduced by a point mass lens , @xmath124 .
the source size of a type  i supernova is @xmath125 cm at maximum luminosity . putting the source at @xmath126 and the lens at @xmath127 in a flat universe with @xmath128 and the dimensionless hubble parameter @xmath129
, we obtain @xmath130 i.e. , only compact objects with @xmath131 should be of importance through microlensing of supernovae .
the work of hw has as a first treatment been concerned with the @xmath86 expressions for idealised cases like point sources ; singular , truncated isothermal spheres ( sis ) ; uniform spheres ; and uniform cylinders .
however , another often - used matter distribution is the one based on the results of detailed n - body simulations of structure formation given by navarro et al .
( @xcite ) .
they parametrize the density profile by @xmath132^{(\beta-\gamma)\delta } } , \ ] ] where @xmath133 is the critical density , @xmath134 is a dimensionless density parameter and @xmath135 is a characteristic radius .
nfw found that halos found in n - body simulations ranging in mass from dwarf galaxies to rich galaxy clusters can be fit by eq .
( [ eq : nfw ] ) with @xmath136 , @xmath137 , @xmath138 . the potential for this density profile
is given by @xmath139 where @xmath140 .
the matrix @xmath86 can then be obtained analytically , see [ app : jnfw ] .
the mass inside radius @xmath141 of a nfw halo is given by @xmath142,\ ] ] since we want to have the average mass density in each cell equal to that of the background cosmology , we get the relation @xmath143 setting @xmath144 , we obtain @xmath145},\ ] ] where @xmath146 .
that is , for a given mass @xmath147 , @xmath148 is a function of @xmath149 . from the numerical simulations of nfw we also get a relation between @xmath148 and @xmath149
this relation is computed numerically by a slight modification of a fortran routine kindly supplied by julio navarro .
of course , one wants to find an @xmath149 which gives both the required average density in each cell and compatibility with the numerical results of nfw .
generally , @xmath149 will be a function of mass @xmath147 , the hubble parameter @xmath150 , @xmath105 , @xmath151 and @xmath0 .
however , we use the result from del popolo ( @xcite ) and bullock et al .
( @xcite ) that @xmath149 is approximately constant with redshift .
a comparison of luminosity distributions between point masses , isothermal spheres ( sis ) and the nfw matter distribution is presented in sect .
[ sec : results ] .
as discussed above , hw have shown that when studying luminosity distributions and other quantities where the null geodesics only differ by infinitesimal amounts , their method is insensitive to the individual masses and clustering of point mass lenses . for point mass lenses ,
the mass distribution and number density will therefore be of small importance .
however , realistic modeling of galaxies does call for realistic mass distributions and number densities , i.e. , one has to allow for the possibility of @xmath90 to reflect the actual distances between galaxies , and for galaxy masses to vary in agreement with observations .
one of the advantages of the method of hw is that any mass distribution and number density , including possible redshift dependencies , can easily be implemented and used as long as the average energy density in each cell agrees with the underlying fl model .    to exemplify this
, we derive a galaxy mass distribution , @xmath152 , by combining the schechter luminosity function ( see , e.g. , peebles @xcite , eq .  5.129 ) @xmath153 with the mass - to - luminosity ratio ( see , e.g. , peebles @xcite , eq .  3.39 ) @xmath154 normalizing to a `` characteristic '' galaxy with @xmath155 and @xmath156 , this can be written @xmath157 using eq .
( [ eq : schechter ] ) , we find that @xmath158 in fig .
[ fig : mdistr ] , mass distributions for various values of @xmath159 are displayed .
assuming that the entire mass of the universe resides in galaxy halos to @xmath160 at the truncation radius .
] we can write @xmath161 using the schechter luminosity function and the mass - to - luminosity ratio we get @xmath162 thus , by supplying values for @xmath163 , reasonably well - determined by observations , and @xmath164 and @xmath165 , on which the dependence of @xmath166 is weak , together with parameters @xmath159 and @xmath167 we can obtain an @xmath166 consistent with @xmath105 .    note that assuming that the schechter luminosity function and the mass - to - luminosity ratio is valid over the entire luminosity range ( i.e. , putting @xmath168 and @xmath169 ) , we get @xmath170 where @xmath171 .    a compilation of parameter values can be found in [ app : parametervalues ] . for @xmath172 , @xmath173 and @xmath174 , eq .
( [ eq : mref2 ] ) can be written @xmath175      the sis model is built on the assumption that stars ( and other mass components ) behave like particles in an ideal gas , confined by their spherically symmetric gravitational potential .
the density profile is given by @xmath176 where @xmath177 is the line - of - sight velocity dispersion of the mass particles .
[ note that this corresponds to eq .
( [ eq : nfw ] ) with @xmath178 , @xmath179 and @xmath180 . ] the mass of a sis halo truncated at radius @xmath95 is then given by @xmath181 we want this to be equal to the mass given by the schechter distribution @xmath147 , @xmath182 where , in addition to @xmath183 , we have introduced a characteristic velocity dispersion @xmath184 . combining the faber - jackson relation @xmath185 where @xmath186 is defined in eq .
( [ eq : schechter ] ) , with the mass - to - luminosity ratio , eq .
( [ eq : masstolum ] ) , we can substitute for @xmath177 in eq .
( [ eq : d ] ) , @xmath187 for values of @xmath184 and @xmath69 , see [ app : parametervalues ]
. using eq .
( [ eq : mrefestimate ] ) , we can write the truncation radius for a halo with mass @xmath188 as @xmath189    a comparison between point masses and isothermal spheres for sources at @xmath190 in a @xmath109 , @xmath191 universe is presented in fig .
[ fig : iso - nfw ] , where the full line corresponds to matter in the form of point masses .
this plot agrees well with fig .  5 of hw .
the long - dashed line corresponds to matter clumped in isothermal spheres with cut - off radius @xmath192 kpc .
we choose this fix @xmath95 in order to be able to compare directly with fig .
both lines agree very well with the results in hw .
finally , the short - dashed line corresponds to isothermal spheres for which the cut - off radius @xmath95 has been calculated according to eq .
( [ eq : dcalc ] ) ( see sect . [
sec : siscutoff ] ) .
for the presented model , this results in a typical cut - off of @xmath193 mpc , which explains why this case comes closer to the filled - beam value .
in figs .  [ fig : model-1][fig : model-9 ] , we compare the luminosity distributions obtained for nfw halos with the point - mass and sis matter distributions . comparing fig .
[ fig : model-1 ] with fig .
18 and 20 of hw , we see that the point - mass cases agree very well .
the sis luminosity distributions differs in that our results is shifted towards the filled beam value .
this is due to the fact that the truncation radii used [ see eq .
( [ eq : destimate ] ) ] are an order of magnitude larger than the ( constant ) value of @xmath192 kpc used in hw .    lensing effects are governed by the lensing potential of the galaxy model used , and the number of lenses between the source and the observer . as noted in the discussion of fig .
[ fig : xmatter ] , a larger number of lenses will tend to make lensing effects for a given galaxy model more prominent . for a fixed source redshift and a fixed number density of lenses ,
the number of lenses will only be a function of the cosmology .
we get the largest number of lenses in the @xmath194 cosmology , followed by the @xmath195 case , and the smallest number of lenses for @xmath196 .    for point masses
, the lensing potential will be unaffected by the cosmology . as discussed in sect .
[ sec : hwsummary ] , point masses cause the largest lensing effects for a given average mass density .    for a sis halo , the truncation radius , @xmath95 , is proportional to @xmath105 [ see eq .
( [ eq : mref ] ) and eq .
( [ eq : dcalc ] ) ] . since the potential for a truncated sis halo is equal to the point - mass potential for @xmath197 ( see fig .
[ fig : sispotentials ] ) , the sis luminosity distribution will be shifted towards the point - mass distribution when lowering @xmath105 .
in fact , the agreement between our results and those in hw is better for lower @xmath105 , as can be seen by comparing the curve corresponding to the sis luminosity distribution in fig .
[ fig : model-2 ] with fig .
taking also the number of lenses into consideration , we expect the largest lensing effects for the @xmath194 cosmology , followed by the @xmath195 case , and the smallest effects in the @xmath196 cosmology .
for nfw halos , the value of the characteristic radius @xmath135 will depend on the cosmology used , the largest value obtained with the @xmath195 cosmology , followed by @xmath194 and the @xmath196 case .
now , it can be shown that for a fixed halo mass , @xmath147 , the nfw potential will approach the point - mass potential when lowering @xmath135 , see fig .
[ fig : nfwpotentials ] , @xmath198 . ] .
therefore , we expect the closest resemblance to the point - mass case for @xmath196 , followed by @xmath194 and @xmath195 .
the method of hw , or rather the numerical implementation of it , has some limitations due to the approximations used .
since the area of the light beam when exiting a cell is computed by linear interpolation from the value before entering the cell [ see eq .
( [ eq : ada ] ) ] , we need very small cells or very many cells to get robust results .
small cells means shorter interpolations while a large number of cells tends to average out the error in every interpolation .    in order to quantify this effect
we have made some tests with empty and filled beams using cell sizes obtained from the schechter matter distribution . in fig .
[ fig : error ] , we have plotted the mean fractional error in area vs. redshift . to get a mean fractional error less than one percent , one would need to have @xmath199 . note that this is the average error of individual measurements ; the average value will be very close to the real value even at very small redshifts . for the cosmology used ( @xmath200 ) , a redshift of @xmath201 corresponds to @xmath202 cells
another drawback with the hw method is that there is no way to keep track of multiple images from the same source .
this poses no problem when the time delay between a primary image and secondary images is big enough to allow observational separation of the images .
however , if image separation is not observationally feasible , we have to compensate for multiple images . this will be the case when considering microlensing , see sect .
( [ subsec.ola ] ) .    a typical case of multiple imaging is shown in fig .
[ fig : focusing4 ] .
note that the secondary image has undergone a caustic crossing and that its parity has been reversed ( this corresponds to @xmath203 ) .
also , the secondary image is demagnified compared to the primary image ( remember that it is the relation between @xmath57 and the area that determines the magnification ) .
since the quantity @xmath57 is more or less constant in our method , areas of secondary images often becomes very big , i.e. , very dim .    for a point mass lens ,
the time delay between two images with flux ratio @xmath141 is given by @xmath204 where @xmath205 for a lens at redshift @xmath127 and a flux ratio of two , we get @xmath206 as an illustration , to get a time delay of the order of weeks , we need @xmath207 .
the image separation in the same case ( assuming a source at @xmath126 ) is given by @xmath208 i.e. , a mass of @xmath207 will give an image separation of @xmath209 .
in the realistic case that a galaxy is causing the lensing , the surface brightness of the lens may make the detection of secondary images impossible ( porciani & madau @xcite ) .      now , let us discuss what kind of information we can obtain from secondary images . in fig .
[ fig : points ] we see what a source would look like in a universe filled with point - mass lenses .
since the deflection angle from a point mass can be very large , there will , apart from the primary image , be one secondary image next to every point mass lying in the vicinity of the line - of - sight to the source ( not to mention all images with two or more caustic crossings ) .. ] most of these images will be _ very _ dim , i.e. not observable .    in our implementation of the hw method , we can happen to follow any one of these .
however , the vast majority of the flux will almost always be in the primary image , the only exception being the case where we have a point mass very close to the line - of - sight ( see fig .  [
fig : points2 ] ) . in this case , the primary and the brightest secondary image will have almost equal flux . if we want to consider the total flux from a source we need to correct for these images .
of course , the correction would have to be larger in the case of fig .
[ fig : points2 ] .
multiple - image corrections are discussed in sects .
[ subsec.ola ] and [ subsec.olsis ] below .
note that if we happen to follow any of the weak secondary images in fig .
[ fig : points2 ] , we do not have any information on the flux of the primary image , i.e. , we would not be able to discriminate between the case of fig .
[ fig : points ] and fig .
[ fig : points2 ] . from this
we conclude that weak secondary images does not give any valuable information concerning the total flux of a source .      filling the universe with sis - type galaxies
, the situation would be quite different .
no longer can we have very large deflection angles , and multiple images will therefore be less common .
, where @xmath177 is the velocity dispersion of the lens and @xmath210 is the impact parameter , i.e. , the magnitude of the deflection does not depend on @xmath210 . ] since we are typically dealing with galaxy masses , images will be separable and we need only to consider secondary images bright enough to be observable . since for ( almost ) every source
, there will be at most one secondary image , it is a more straightforward procedure to associate to every secondary image a primary image .      here , we will discuss reasonable corrections when images are non - separable .
since we have a very large number of images with two or more caustic crossings contributing to the flux of the total image , we can not correct for each one of them .
instead , we will make a first correction for the brightest secondary image using the approximation that lensing effects from a single lens ( cell ) is dominant .
this is the so called one - lens approximation .
possible tests of the validity of this assumption is discussed in sect .
[ sec : cc ] .
its plausibility is also reinforced by figs .  8 and 9 in hw .
we can now use standard analytical expressions for one - lens systems for this correction : @xmath211 in fact , this formula is valid also for secondary images but one then has to include the sign in @xmath58 , i.e. , @xmath212 .
however , eq .  ( [ eq : mucorr ] ) is only valid for the brightest secondary image , and one should therefore not apply it to very weak , secondary images . also note , that to make the correction self - consistent , the magnification should be defined as the ratio of the expected area if the most prominent cell is empty and the real area .
it is easy to see that the most important effect of this correction is to make the high magnification tail more pronounced .
it is as straightforward to compute the total magnification using the one - lens approximation with a sis - lens . for a primary image
we get @xmath213 for a secondary image we get @xmath214 since images are separable in the majority of cases , this is perhaps most interesting to use as a consistency check for the total luminosity .
in order to validate our implementation of the method , we have performed a number of tests .    *
* filled and empty beam distances * + first  in order to get an independent check  we have compared results from our implementation of the hw method using empty cells and cells with a homogeneous dust component with analytical empty - beam and filled - beam results ( using the angsiz routine , see kayser et al .
this has been done for a variety of cosmologies and the agreement is good even at very high @xmath0 ( less than 1 % discrepancy up to @xmath215 ) . * * comparison with hw * + second , we have reproduced most of the plots presented by hw .
we have also tested that primary images satisfy @xmath216 , see fig .
[ fig : areasum ] . * * average luminosities * + imagine a source situated at the center of the sphere in the lower panel of fig .
[ fig : focusing3 ] and observers sitting on the sphere at some constant redshift .
since the number of photons emitted by the source , as well as the area of the observers sphere is constant , i.e. independent of any details of the distribution of matter as long as the average matter density is the same , observers will on average receive the same flux , regardless of the matter distribution .
this average will be equal to the fl value , @xmath217 when computing @xmath58 for a source , one has to consider the total flux if there are multiple images .
if we consider primary images only , without any corrections , we would expect to obtain a @xmath58 slightly lower than the fl value ( since we lose some flux in secondary images ) .
+ note that this does not mean that the average area should equal the fl area . instead , @xmath218 i.e. , the fl area is the harmonic mean of the event areas , not the arithmetic ( remember that @xmath219 ) .
since we lose some flux in secondary images , we have ( considering primary images only ) @xmath220 see fig .
[ fig : lumtest ] for a plot of the average magnification using point mass lenses and nfw lenses .
the result for sis halos is not plotted since the result is indistinguishable from the nfw result
. note that since the deviation from the filled - beam value is due to the absence of the flux contained in secondary images , fig .
[ fig : lumtest ] indicates that the flux in secondary images is very low when using nfw or sis halos , even at very high redshifts .
+ the relationship between this test and the area test of hw is discussed in [ app : meanarea ] .
since the numerical work involved in computing the results presented in this paper is quite extensive , it may be convenient to use simple analytical fitting formulas for the probability distributions for the different cosmologies , which we now give .
the probability distribution of the deviations from filled - beam magnitudes induced by sis and nfw lenses , @xmath221mag(filled beam)@xmath222mag(lensed ) , can be parametrized by an exponential function folded with a gaussian distribution : @xmath223    thus , for any combination of cosmological parameters , the corrections due to gravitational lensing can be described by three parameters : the mean of the gaussian distribution , @xmath58 , the standard deviation @xmath224 , and the magnification tail , @xmath225 .
the expression in eq .
( [ eq : gausexp1 ] ) can be rewritten so that ( after normalization ) a simpler integral is invoked :    @xmath226    where @xmath227    in fig .
[ fig : sisdis ] the fits to the magnification distributions are shown for sis lenses for source locations @xmath126 , @xmath228 and @xmath229 and for three different cosmologies . in fig .
[ fig : nfwdis ] the same source locations and cosmologies are used but the lenses are of nfw type .    the parametrization of eqs .
( [ eq : gausexp1 ] ) and ( [ eq : gausexp2 ] ) may be used to study the shape of the magnification distribution as a function of redshift and cosmological parameters , as shown in figs .  [
fig : sispar ] and [ fig : nfwpar ] .
the estimates of the cosmological parameters may be further refined and the understanding of the mass distribution in galaxy halos could be deduced by fitting the parameters @xmath225 , @xmath224 and @xmath58 to the residual magnitude distributions in future high - redshift supernova searches , such as the snapsat project or at ngst .
the accuracy of this method will be discussed in a forthcoming paper ( bergstrm et al .
, in preparation ) .
in this paper , the method of holz & wald ( @xcite ) has been generalized to include the effect of non - vanishing pressure .
a motivation for this exercise is to use this method as part of a large package for a detailed simulation of most aspects high - redshift supernova determination of the geometry of the universe ( bergstrm et al .
, in preparation ) .
the possibility to allow for various matter contributions is then desirable .
another improvement has been to allow for a more realistic mass profile of the galaxy halos , the navarro - frenk - white distribution , where we have given analytical formulas for the geodesic deviation integrals . with this improvement , and by adding matter in our model universes according to a schechter function with parameters fixed by observation , we think we have achieved a sufficient degree of sophistication to make realistic predictions of the lensing effects of present and future deep supernova searches .
our results show , in agreement with previous work ( holz & wald @xcite ; metcalf & silk @xcite ) that the ( unlikely ) case of having most of the dark matter distributed in the form of point masses can be quite easily distinguished from more realistic halo models , if a sufficient number of supernovae can be detected at large redshifts .
it will be more difficult to distinguish between a singular isothermal sphere model and the n - body results of navarro , frenk and white , unless @xmath105 is large and one can measure supernovae out to @xmath230 ( see fig .  [ fig : model-7 ] ) .
however , this also means that our method of determining the properties of the lens population by matching the measured schechter function with the required @xmath105 gives quite robust results for the predicted magnifications of type ia supernovae over a broad range of reasonable halo models .
this will be helpful when making realistic predictions for future deep supernova searches .
the authors would like to thank daniel holz and bob wald for helpful comments , and julio navarro for providing his numerical code relating the parameters of the nfw model .
in this appendix , we give analytic expressions for the integrals ( [ eq : jab ] ) needed for the evolution of the optical matrix for a bundle of light which passes near a mass distribution of the nfw form .
we assume a spherically symmetric matter distribution , the impact parameter is @xmath89 , the characteristic radius of the nfw model is @xmath135 and the cell radius is @xmath90 .
when @xmath231 , it follows that @xmath232    \nonumber \\ & &     -\frac{\sqrt{r_c^2-b^2}(2r_c^2+b^2)}{b^2r_c^3 }    \ln\left(1+\frac{r_c}{r_s}\right )     \nonumber \\ & & \left .
+ \frac{2r_s(r_s^2 - 2b^2)}{b^2|r_s^2-b^2|^{3/2}}f(b , r_s , r_c)\right\ } , \end{aligned}\ ] ] @xmath233    \nonumber \\ & &    -\frac{\sqrt{r_c^2-b^2}(2r_c^2+b^2)}{b^2r_c^3 }    \ln\left(1+\frac{r_c}{r_s}\right )    \nonumber \\ & & \left .
+ \frac{2r_s\,{\rm sgn}(r_s - b)}{b^2|r_s^2-b^2|^{1/2}}f(b , r_s , r_c)\right\ } , \end{aligned}\ ] ] @xmath234 and the limit @xmath235 is : @xmath236    \right.\nonumber\\    & & -\frac{\sqrt{r_c^2-b^2}(2r_c^2+b^2)}{b^2r_c^3 }    \ln\left(1+\frac{r_c}{b}\right)-   \nonumber\\    & & \left.\left(\frac{r_c - b}{r_c+b}\right)^{3/2 }    \frac{r_c^2 + 3(r_c+b)^2}{3r_c^2b^2 } \right\ } ,    \nonumber\\    ( j^{\rm nfw})^y\!_y & = & -2 c \left\ {    \frac{\sqrt{r_c^2-b^2}}{r_c^2(r_c+b ) }    \right .
\nonumber \\ & &    + \frac{2}{b^2 }    \ln\left[\frac{r_c}{b}+\sqrt{\left(\frac{r_c}{b}\right)^2 - 1}\right ]    \nonumber \\ & &    -\frac{\sqrt{r_c^2-b^2}(2r_c^2+b^2)}{b^2r_c^3 }    \ln\left(1+\frac{r_c}{b}\right )    \nonumber \\ & & \left .
-2\frac{\sqrt{r_c^2-b^2}}{b^2(r_c+b ) } \right\ } .\end{aligned}\ ] ] the constant @xmath237 is given by @xmath238 where @xmath239 .
here , we have compiled various parameter values , relevant for the mass distribution discussed in sect .  [ sec : massdist ] . in
what follows , absolute magnitudes are denoted by @xmath240 , while masses are denoted by @xmath147 .
* lin et al . * : :    see eq .
( 32 ) of ( lin et al .
@xcite ) .
+    * @xmath241    * @xmath242    * @xmath243    +    lin et al . gives its range of validity as    @xmath244 , which is about    @xmath245 magnitudes around their characteristic magnitude    @xmath246 .
this translates into a luminosity range    @xmath247 , and a mass range    @xmath248 , which , with    @xmath249 , gives @xmath250 $ ] .
* kirshner et al . * : :    ( kirshner et al .
@xcite ; ki )    +    * @xmath251 ( ki19 ) +    @xmath252 ( ki22 )    * @xmath253 ( ki18 ) +    @xmath254 ( ki21 )    * @xmath255 ( see binney &    tremaine @xcite ) .
* @xmath256 ( ki23 ) * peebles * : :    ( peebles @xcite ; pe )    +    * @xmath174 ( pe3.39 )    * @xmath257 ( pe5.130 )    * @xmath258 in the    @xmath259 magnitude system where    @xmath260 .
( pe5.139,5.140 )    * @xmath261    ( pe5.141 )    * @xmath262    ( pe5.142 )    * @xmath263 ( pe3.35 )    * @xmath264 ( pe3.35 )    in the numerical work in this paper , we have used    * @xmath174 * @xmath173 * @xmath265 * @xmath266 * @xmath267 * @xmath263 * @xmath264
holz and wald ( @xcite ) use the fact that the area of the boundary of the past of an observer should be very nearly equal to the area of the past light cone in the underlying fl model . from this
they deduce that adding up all areas  counting beams with an even number of caustics as positive and beams with an odd number as negative  one should very nearly obtain the fl result , i.e. , @xmath268 , where the zero - subscript denote unweighted quantities . here
we investigate whether this claim is consistent with the fact that the average magnification should be equal to the fl magnification .
imagine that , without weighting , we get an area distribution @xmath63 , where we assume that lensing effects are weak enough to only give primary images . according to hw , we then have @xmath269 in our simulations , the area distribution will be a convolution of @xmath63 and a probability distribution @xmath270 , @xmath271 since we want the harmonic mean of areas to be equal to the fl value , we get @xmath272 that is , the two consistency checks at least show consistency in the case where we only have primary images .
the inclusion of non - primary images is highly non - trivial since the concept of negative areas is difficult to interpret as luminosities .
one possibility is to combine this approach with the one - lens approximation method for multiple image correcting .
however , this is difficult to do analytically since one would then need the form of @xmath63 .


massive multiple - input multiple - output ( mimo ) ( also known as large - scale mimo ) is considered as one of the most promising technologies for next generation wireless communication systems @xcite , since it can provide very high spectral and energy efficiencies . besides , for massive mimo systems , the effects of small - scale fading and interference can be significantly reduced with linear signal processing , such as maximal - ratio transmission ( mrt ) and zero - forcing ( zf ) @xcite .
the performance of massive mimo single / multi - cell downlink systems has been widely studied in the literature . in @xcite , the asymptotic signal - to - interference - plus - noise - ratio ( sinr ) for linear precoding in single - cell _ frequency - division - duplex _ ( fdd ) systems was derived based on which the optimization of the user number was investigated . in @xcite , the power scaling ,
i.e. , how fast the transmission power can decrease with the antenna number while maintaining a certain sinr level , for multi - cell _ time - division - duplex _
( tdd ) systems was given based on asymptotic sinr analysis , where the pilot contamination was shown to be the only performance limit . in @xcite , capacity lower bounds and energy efficiency of mrt and zf in single - cell tdd systems were studied and compared with each other . in the parallel field of uplink systems , similar results were derived in @xcite , @xcite .
these works mainly focus on the derivations of asymptotically deterministic expression for the sinr , and the performance analysis only consider some special cases , e.g. , the constant or linearly increasing user number with respect to the base station ( bs ) antenna number , the power scaling for perfect or imperfect channel state information ( csi ) with equal training and transmission power and constant user number .    recently , along this route , the performance improvement for using coherent detection based on the dedicated downlink channel estimation was shown in @xcite , especially for large but finite antenna numbers .
meanwhile , the impact of other practical factors such as channel aging @xcite , hardware limits @xcite , and rician fading @xcite on the performance of massive mimo systems was studied . however , a comprehensive scaling law analysis for multi - cell multi - user massive mimo is still missing in existing works .
another important problem in the performance analysis of massive mimo is when the derived asymptotic results are applicable in practical systems .
while many works derived deterministic sinr approximations via calculating the expected value of the sinr or its compenents , the analysis on whether the instantaneous sinr becomes asymptotically deterministic and converges to the derived expressions is sometimes missing . in @xcite , this problem was studied through the so called massive mimo efficiency for the ideal situation that the noise , channel estimation error , and interference are negligible compared with pilot contamination , where the dependence of the massive mimo efficiency on the ratio of spatial degrees of freedom / antenna number to the user number and the product of antenna number and transmit power was shown by simulations . in @xcite , via simulations , the convergence speed of the strong law of large numbers with respect to the antenna number and the convergence of the sinr of mrt and zf to their expected values were both analyzed for fixed ratio of antenna number to user number .
it was shown that the sinr of mrt has slower convergence than that of zf , while the explicit relationship between the convergence speed and important parameters , i.e. , antenna number , user number , transmit power , and channel training energy ( csi quality ) , is missing .    in this work , by drawing lessons from the analyzing method and the mean square convergence definition in @xcite , we provide a comprehensive performance scaling law analysis for multi - cell multi - user tdd massive mimo downlink systems with channel spatial correlation , csi error , and pilot contamination and investigate the applicability of the derived asymptotic scaling law in practical systems with large but finite antenna numbers .
extensive simulations are also conducted to validate the analytical results .
our main new contributions are summarized as follows .
* we consider a multi - cell multi - user massive mimo network with correlated channel model , csi error , and pilot contamination . in dealing with pilot contamination ,
recent advancements in pilot contamination elimination ( pce ) @xcite is also considered , while in existing works it is treated as a constant bottleneck of massive mimo .
an approximate model for the pce effect is proposed .
* a lower bound on the system sum - rate is derived based on which a general performance scaling law is obtained .
the result reveals the effect and joint interaction of extensive network parameters , i.e , the channel training energy , the transmission power , and the user number . in the scaling - law analysis ,
a general model is used where the parameters have continuous scaling exponents with respect to the number of bs antennas , while in most existing work , only a few discrete values for the power scaling exponent , e.g. , @xmath0 ; @xmath1 ; @xmath2 , are allowed . * to understand the applicability of the derived scaling law in practical systems with large but finite antenna number
, the effect of antenna number on its accuracy is analyzed and evaluated by simulations .
* we conduct quantitative analysis on the asymptotically deterministic property of random instantaneous sinr using the mean square convergence .
sufficient conditions for the instantaneous sinr to be asymptotically deterministic are given , from which the tradeoff among different parameter scales is discovered .    in @xcite ,
the performance scaling law for multi - user massive mimo relay networks was extensively studied .
while our work adopts the mean square convergence idea and the sinr analysis method from @xcite , it is fundamentally different in both network model and channel model .
we consider the multi - cell multi - user massive mimo downlink networks with a general spatially correlated channel model .
furthermore , as crucial issues for multi - cell systems , the pilot contamination and pce are taken into account in our work . as a result , the corresponding derivation procedures and analytical results are significantly different . in addition , the analysis on the effectiveness of the derived scaling law for practically large but finite antenna numbers is completely new .
the remaining of the paper is organized as follows . in section
ii , the channel model and system model including both the multi - cell channel estimation with pilot contamination and the downlink transmission with mrt are introduced .
section iii shows the performance scaling law result , discussions on typical network scenarios , and the analysis on the applicability to practical systems with large but finite antenna numbers . in section iv ,
sufficient conditions for the sinr to be asymptotically deterministic are derived .
section v shows simulations and conclusions are provided in section vi .    in this paper , bold upper case letters and bold lower case letters
are used to denote matrices and vectors , respectively . for a matrix @xmath3
, its conjugate transpose , transpose , and trace are denoted by @xmath4 , @xmath5 and @xmath6 , respectively . @xmath7 and @xmath8 are the mean and variance operators .
@xmath9 means equivalence in probabilistic distribution .
the notation @xmath10 means that when @xmath11 , @xmath12 and @xmath13 have the same scaling with respect to @xmath14 , in other words , for large @xmath14 , there exists positive constants @xmath15 and @xmath16 , such that @xmath17 .
consider a multi - cell multi - user massive mimo network consisting of @xmath18 cells with one @xmath14-antennas bs and @xmath19 scheduled single - antenna users in each cell .
the bss and users are assumed to operate a tdd protocol with universal frequency reuse .
the uplink spatially correlated channel from user @xmath20 in cell @xmath21 to the antennas of bs @xmath22 in cell @xmath22 can be written as @xmath23 where @xmath24 is the channel correlation matrix and @xmath25 is the independent fast - fading channel vector .
the following assumptions on @xmath24 are made by taking into consideration the pilot contamination , pce , and tractability analysis .
the set @xmath26 denotes the index set of all cells having pilot contamination to cell @xmath22 while @xmath27 denotes the index set of remaining cells whose pilot contamination is eliminated .
we assume that for all @xmath28 and @xmath29 , @xmath30 where the channel direction matrix @xmath31 is an @xmath32 unitary matrix and @xmath33 with @xmath34 $ ] indicates the channel correlation level .
for example , elements of the channel vector become independent and identically distributed ( i.i.d . )
when @xmath35 .
also , @xmath36 is the pass - loss fading factor where @xmath37 and @xmath38 .
the parameter @xmath39 denotes the large - scale fading between different cells and the intra - cell one is normalized to @xmath1 .    for @xmath40 , the pilot contamination from cell @xmath21 to cell
@xmath22 is eliminated by using the spatial information in channel correlation matrix .
thus we have @xmath41 this channel model considers both the channel correlation and the inter - cell pass - loss fading .
it is similar to the one used in ( * ? ? ?
* section iii ) where motivations and justifications for this model can be referred to .
new feature of the modeling in and is that the possible effect of pce , i.e. , the bss schedule users with relatively less pilot contamination for channel training and transmission , can be embodied . in @xcite , several efficient methods by exploiting the statistical channel information have been proposed to alleviate the pilot contamination due to the pilot reuse among @xmath42 cells @xcite , @xcite . generally speaking , the overlap among users channel directions can be in any degree leading to different pilot contamination level @xcite . for explicit analytical results , in
, we assume that channels from both users in cell @xmath22 and users in other cells with pilot contamination to cell @xmath22 have the same @xmath3 matrix , while in other cells , pilot contamination is eliminated by having the orthogonal channel directions as shown in .
this is a discretization approximation to show the average effect of pce and the extension to the case of more general channel direction distribution is left for future research .      in the uplink training phase ,
users in cell @xmath21 transmit mutually orthogonal pilot sequences which are reused by users in other cells .
denote the length of the pilot sequence as @xmath43 ( @xmath44 for reliable performance @xcite ) .
assume that all @xmath19 users in each cell use the same training power @xmath45 .
define @xmath46 which is the channel training energy . by correlating the received training signals with the pilot sequence of user @xmath47 ,
bs @xmath21 has the observation @xmath48 as @xmath49 where @xmath50 is the local noise . with the minimum - mean - squared - error ( mmse ) estimation @xcite , the estimate @xmath51 is given as @xmath52 we assume that all the @xmath42 cells have the same pilot contamination level , i.e. , @xmath53 $ ] for all @xmath21 . by defining @xmath54 and utilizing ,
, we have @xmath55 , for @xmath56 where @xmath57 from , it can be shown that @xmath58 this will be utilized to simplify the derivations in later parts .
define @xmath59 , which is the estimation error . due to the feature of mmse estimation and gaussian distribution , @xmath60 and @xmath61
are mutual independent and @xmath62 where @xmath63 note that @xmath64 denotes the csi quality , which is dependent on both channel training energy and pilot contamination . when @xmath65 , we have @xmath66 and the csi quality is only affected by pilot contamination .
if further perfect pilot contamination elimination is achieved , we have @xmath67 , @xmath68 and perfect csi can be obtained .
the mrt precoding is considered due to its computational simplicity , robustness , and high asymptotic performance @xcite . by assuming channel reciprocity , the received signal @xmath69 of user @xmath20 in cell @xmath21 is given as @xmath70 where @xmath71 denotes the average downlink transmission power or the signal - to - noise - ratio ( snr ) , @xmath72\in\mathbb{c}^{m \times k}$ ] is the mrt precoding matrix , @xmath73^t\in\mathbb{c}^{k}\sim\mathcal{cn}({\bf0},{\bf i}_k)$ ] contains the data symbols for the @xmath19 users in cell @xmath22 , and @xmath74 is the receiver noise .
further , @xmath69 can be rewritten as @xmath75 in this formulation , the instantaneous effective gains of the estimated channel vector , i.e. , @xmath76 , are assumed to be known at the user side for coherent detection . in several recent works , downlink training for the estimation of the instantaneous effective gains is considered @xcite . in this work , we use @xmath77 as an approximate estimate of @xmath78 to indicate the effect of instantaneous effective channel gain and at the same time to delete the insignificant details of the downlink training to better focus on the main purpose of this paper .
this assumption has been used in recent works on tdd massive mimo relay networks @xcite .
since the bs contaminated by one user will create interference to it in the downlink transmission , we define @xmath79 and @xmath80 which are the index sets of cells or bss which are contaminated and not contaminated by users in cell @xmath21 in the uplink training period , respectively . for simplicity , we assume @xmath81 , for all @xmath21 which is reasonable in the average sense .
define @xmath82which are respectively the normalized power of the desired signal , intra - cell interference , inter - cell interference , and effective noise created by unknown csi error .
their different normalization factors are used to guarantee that the mean of each term is bounded as shown in the next section . in the derivations of and we utilize @xmath83 for @xmath84 based on , and . with the results in - , the sinr of user @xmath85 in cell @xmath21 can be written as : @xmath86 the achievable ergodic rate is thus @xmath87
this section studies the general performance scaling law of the massive mimo network .
first , the means and scvs of all sinr components are calculated .
consequently , a sum - rate lower bound based on jensens inequality is derived via exploiting the asymptotically deterministic property of the desired signal power .
then , the general performance scaling law is obtained and typical network scenarios for the sinr to be non - decreasing are discussed . finally , conditions for the scaling law to be applicable for large but finite antenna numbers and extended discussions on the possible effect of pce are both given .
different from directly using the almost sure convergence based lemmas @xcite , @xcite , @xcite to calculate the deterministic equivalence of , we consider the mean square convergence @xcite and use the scv ( defined as the square of the ratio of the standard deviation over the mean of a random variable ) in the analysis @xcite .
specifically , for a random variable sequence @xmath88 with a bounded mean , it converges in mean square to a deterministic value @xmath89 , denoted as , @xmath90 , if @xmath91 .
further , @xmath88 is said to be asymptotically deterministic in the sense of mean square convergence if @xmath92 @xcite . to utilize the mean square convergence ,
we first calculate the means and scvs of all sinr components . the following lemma is obtained .
[ lemma1 ] @xmath93    see appendix [ appendix a ] .
[ remark 1 ] by noticing that @xmath94 $ ] , and @xmath95 $ ] , it can be seen that the random variables @xmath96 , @xmath97 and @xmath98 , @xmath99 , all have bounded means . moreover , apparently the scv of @xmath100 approaches 0 as @xmath11 with a fast linear convergence rate .
thus for large @xmath14 , @xmath96 can be approximated with its mean .
however , for @xmath97 , @xmath98 and @xmath101 , their scvs depend on the scaling of @xmath19 , @xmath102 and @xmath103 which are not necessarily @xmath104 .
thus , at this point , we do not use their mean values for further derivations , but keep their random nature for more careful analysis .
therefore , the sinr expression in becomes @xmath105 based on which the following result on the achievable rate can be obtained .
[ lemma 2 ] the achievable rate of user @xmath20 in cell @xmath21 in the massive mimo network has the following lower bound : @xmath106 where @xmath107    as @xmath108 is a convex function of @xmath89 @xcite , according to jensen s inequality , we have @xmath109 by applying the sinr approximation in , it can be shown that @xmath110 by defining the effective asymptotic sinr as @xmath111 and directly using the results in - , the lower bound in can be obtained .    in existing works
@xcite , closed - form approximations of the ergodic rate were derived by using ( * ? ? ?
* lemma 1 ) , i.e. , @xmath112 .
however , in lemma [ lemma 2 ] we utilize the definite asymptotically deterministic property of the desired signal power term only in the sinr expression to derive the lower bound of the ergodic rate . due to ,
understanding the scaling law of the achievable sum - rate is transformed to understanding the scaling law of @xmath113 where the subscript of user index and cell index in is omitted due to homogeneous network assumptions .      in this subsection ,
the scaling law of the asymptotically effective sinr is analysed to show how the system performance is affected by network parameters and pilot contamination level . for all system parameters ,
including the number of users @xmath19 , the transmission power @xmath71 and the channel training energy @xmath102 , a general scaling model with respect to the bs antenna number @xmath14 is used .
note that other system parameters , i.e. , the channel spatial correlation metric @xmath114 , the normalized inter - cell large scale fading @xmath115 are constant with respect to @xmath14 .
assume that @xmath116 where the exponents @xmath117 , @xmath118 , and @xmath119 represent the scales of @xmath19 , @xmath120 , and @xmath121 with respect to @xmath14 ; and @xmath122 represents the scale of the effective sinr .
for practical ranges of the network parameters , it is assumed that @xmath123 .
the reasons are given as follows .
firstly , in typical applications of massive mimo , @xmath19 either increases with @xmath14 or keeps constant .
thus @xmath124 . on the other hand
, @xmath19 can not exceed @xmath14 since the maximum multiplexing gain is @xmath14 .
thus , @xmath125 . secondly , due to the high energy efficiency requirement of massive mimo@xcite , the transmission power @xmath71 should not increase with @xmath14 .
but it can decrease as @xmath14 increases with the condition that its decreasing rate is lower than the increasing rate of @xmath14 .
this is because that the maximum achievable array gain for compensating the loss of receiving energy is @xmath14 .
thus @xmath126 .
for the same reasons , @xmath127 is considered .
another important parameter for the performance analysis is @xmath103 which quantifies the pilot contamination level .
existing works @xcite , @xcite show that pilot contamination does not vanish as @xmath14 grows for i.i.d .
channel @xcite or spatially correlated channel with uniform channel directions @xcite . for more practical channel models , by exploiting the large spatial degrees of freedom provided by the massive antenna array at the bs , several efficient methods have been proposed to reduce pilot contamination @xcite .
specifically , the difference among users channel directions @xcite or locations @xcite and smart pilot allocation based on them @xcite were utilized to avoid pilot contamination , and corresponding simulations showed that as @xmath14 increases , lower pilot contamination level can be achieved . in our work , the effect of pce is reflected by the @xmath128 parameter in our modeling , which is assumed to be a constant .
the special case of @xmath129 represents perfect pce , i.e. , highly correlated channels .
] , while @xmath130 represents that pilot contamination is not fully annihilated and on average there are @xmath128 cells with pilot contamination to the cell of interest . with the above scaling model , the general performance scaling law
is given as follows .
[ theorem 1 ] for the massive mimo network with mmse channel estimation and mrt , if perfect pce is achieved , the performance scaling law is @xmath131 for imperfect pce , the performance scaling law is @xmath132    the scaling exponent of @xmath133 in is determined by the maximal scaling exponent of the terms in its denominator . directly from , it can be seen that given aforementioned parameter range , either @xmath134 or @xmath135 or both have the highest order with respect to @xmath14 .
further , since @xmath64 and @xmath102 have the same scaling from , ( [ eq30 ] ) and ( [ eq31 ] ) can be obtained from ( [ eq28 ] ) .    since having decreasing performance with respect to @xmath14 contradicts the motivations of massive mimo , we present the condition for _ non - decreasing _ sinr ( i.e. , @xmath136 ) in the following corollary .    [ corollary 1 ] for both perfect and imperfect pce , the necessary and sufficient condition for the massive mimo network with mmse channel estimation and mrt to have non - decreasing sinr is @xmath137\ ] ]    the result is a straightforward extension from theorem [ theorem 1 ] .
the scaling law in ( [ eq30 ] ) and ( [ eq31 ] ) illustrates quantitatively the relationship between the performance scaling and scalings of all important parameters .
the condition in ( [ eq32 ] ) , on the other hand , provides useful guidelines for the design of the massive mimo network .
equation ( [ eq30 ] ) shows that with perfect pce , the scaling of the sinr @xmath122 is determined by the scaling exponent of the channel training energy @xmath119 and per - user transmission power @xmath138 .
moreover , ( [ eq30 ] ) also shows that the scaling of the sinr is a decreasing function of both .
thus higher per - user transmission power and training energy result in improved performance , and one can compensate for the other in performance .
( [ eq31 ] ) shows that with imperfect pce , @xmath122 is upper bounded by zero due to the pilot contamination bottleneck and it can be seen from that @xmath139 .
this means that the asymptotic performance will be constant with respect to @xmath14 , while the concrete value is dependent on the efficiency of pce , which is more general than the result in @xcite .      in this subsection
, the scaling law in ( [ eq30 ] ) and ( [ eq31 ] ) and the condition for non - decreasing sinr in ( [ eq32 ] ) are elaborated for typical network settings .    * [ d1 ] first , we consider the case of @xmath140 , i.e. , the training energy @xmath102 is constant . for perfect pce , @xmath141 . for imperfect pce , @xmath142 .
to have non - decreasing sinr , @xmath143 is needed .
this shows a tradeoff between the scaling of the user number @xmath19 and the transmission power @xmath71 .
the most power - saving design is to make the per - user transmission power decrease linearly with @xmath14 , i.e. , @xmath144 . with a larger @xmath14 ,
the network can serve more users or have less power consumption , while maintaining certain performance .
however , improvements in the two aspects both have a limit : 1 ) @xmath145 , @xmath146 and 2 ) @xmath147 , @xmath148 .
case 1 ) means that when @xmath19 increases linearly with @xmath14 , to achieve non - decreasing sinr , @xmath71 must remain constant , and thus the goal of reducing @xmath71 can not be achieved .
case 2 ) means that when @xmath71 is scaled inversely proportional to @xmath14 , the goal of serving more users can not be achieved .
the latter case is the major power scaling scenario considered in the literature @xcite .
obviously , our results cover this case , and shows more insights for general scalings of @xmath19 and @xmath71 .
* [ d2 ] then we consider the case of @xmath149 , i.e. , @xmath150 . for perfect pce , @xmath151 . for imperfect pce , @xmath152 .
to have non - decreasing sinr , @xmath153 is needed , i.e. , the transmission power and the user number should both remain constant .
this shows that the training energy is key to the performance , i.e. , with low training energy , the promising features of massive mimo can not be achieved . *
[ d3 ] for the general case where @xmath154 , non - decreasing sinr requires @xmath155 .
that is , the scale of the per - user transmission power should be no less than @xmath156 .
this shows the trade - off between the training phase and the data transmission phase . *
[ d5 ] in previous discussions , the scaling of the training energy @xmath119 is treated as a free parameter , next , the special but commonly used training case where @xmath157 and @xmath158 is considered , i.e. , each user uses the same power in training as that for data transmission . since in this case , @xmath159 , we have @xmath160 for perfect pce and @xmath161 for imperfect pce .
non - decreasing sinr requires @xmath162 .
if further @xmath19 is constant , @xmath163 is needed , i.e. , the transmission power should scale as @xmath164 or higher .
this is the same as the results in @xcite , @xcite . *
[ d6 ] another typical setting from the perspective of supporting a large number of users is to have @xmath19 increase linearly with @xmath14 , i.e. , @xmath165 .
for this case , @xmath166 for perfect pce and @xmath167 for imperfect pce .
non - decreasing sinr requires @xmath168 .
thus , to support such a number of users , the training energy must be constant or higher and the transmission power can not decrease with @xmath14 simultaneously .
the performance scaling law is derived for @xmath169 . to be practically useful
, its effectiveness for large but finite antenna numbers needs to be understood . in our scaling law derivation
, the term with the highest scaling exponent in the denominator of is kept while other lower order terms are neglected .
thus for the interested range of antenna number , the result is said to be applicable if the term with the highest order is dominant , in other words , significantly larger ( denoted by the symbol `` @xmath170 '' ) than the rest is said to be dominated by another positive number @xmath171 , i.e. , @xmath172 , when @xmath171 is at least one - order high than @xmath173 , or in other words , @xmath174 .
this common practice is used in our work . ] . in the following ,
we derive conditions for the scaling law to be applicable in practical systems with large but finite antenna numbers .
[ corollary 2 ] for networks with perfect pce , the scaling law in ( [ eq30 ] ) is applicable for large but finite @xmath14 when @xmath175 .    in this case , @xmath129 and the term with the highest order with respect to @xmath14 is @xmath176 . for it to be dominant in the denominate of ( [ eq25 ] )
, it can be seen straightforwardly that @xmath175 is sufficient .
[ corollary 3 ] for networks with imperfect pce , the scaling law in ( [ eq31 ] ) is applicable for large but finite @xmath14 when the following holds : @xmath177 where @xmath178 .
the derivation is almost the same as the one for corollary [ corollary 2 ] .
thus it is omitted .
[ remark 5 ] corollary [ corollary 2 ] shows that for large but finite @xmath14 , the scaling law in ( [ eq30 ] ) for perfect pce is more accurate for less correlated channel ( larger @xmath114 ) and/or lower transmission power .
corollary [ corollary 3 ] shows that for large but finite @xmath14 , the scaling law in ( [ eq31 ] ) for imperfect pce is also more accurate for less correlated channel , while the effect of the transmission power on the scaling law accuracy is unclear .    the sufficient conditions .
thus , another case for the scaling law in ( [ eq30 ] ) to be applicable is @xmath146 and @xmath179 . ] in corollaries [ corollary 2 ] and [ corollary 3 ] for the scaling laws to be applicable provide guidance for using the scaling law and in practical system design with large but finite antenna number . for systems that satisfy the conditions in corollaries
[ corollary 2 ] and [ corollary 3 ] , the derived scaling law can accurately reflect the network performance ; while when the conditions are not guaranteed , the results becomes less accurate but can still provide rough suggestions on the performance .      in the scaling law analysis in section [ scalinglawres ] , @xmath128 ,
the average number of cells with pilot contamination is assumed to be a constant . as also mentioned in that section , a decreasing @xmath103 with respect to @xmath14
is possible by pce . to take into consideration
the possible variation of @xmath128 with respect to @xmath14 , the following approximate scaling model for @xmath128 can be used : @xmath180 where @xmath181 .
the value of the exponent @xmath182 is dependent on the concrete algorithm .
for example , @xmath183 corresponds to the imperfect pce case discussed previously , where no pce is available or the pilot contamination is only reduced to a constant value . since the highest achievable scaling of pilot contamination can not surpass the increasing rate of antenna number , we have @xmath184 .    with this more general modeling , the corresponding performance scaling law by following similar derivation procedure can be obtained as @xmath185    [ remark 6_170211 ] the expressions in and are approximations since @xmath103 is modeled as an integer in this paper .
its variation with respect to @xmath14 can be roughly fitted with the function @xmath186 for a constant @xmath171 .
when @xmath103 decreases with @xmath14 ( meaning high pce efficiency and a positive @xmath182 ) , an increasing effective sinr rather than the bounded one in can be obtained if @xmath187 , which will be evaluated by simulations in section [ simulation ] .
one important concept in massive mimo analysis is asymptotically deterministic . for example , when the desired signal power and interference power which are random in finite - dimension converge to deterministic values as @xmath14 increases to infinity @xcite , their expected values can be used to replace the random counterparts to simplify the analysis .
then a natural question is when the massive mimo system has asymptotically deterministic sinr for the corresponding performance analysis to be reliable . in existing literature ,
the deterministic equivalence is mainly based on the almost sure convergence@xcite .
while many results have been reported based on the aforementioned analysis , the corresponding convergence condition tends to be vague and explicit analysis of the effect of parameter scalings on the convergence is missing .    in this section , following the work in @xcite , we adopt the mean square convergence @xcite for the definition of asymptotically deterministic property , and use the scv scaling with respect to @xmath14 to derive sufficient conditions for asymptotically deterministic sinr .
finally , typical network scenarios are discussed .
[ definition 2 ] the random variable sequence @xmath88 is said to be asymptotically deterministic if its scv decreases linearly with @xmath14 or faster .
strictly speaking , there is no constraint on the scaling of the scv for a random sequence to converge to a constant value in the mean square sense other than that the scv converges to 0 as @xmath188 .
notice that for any positive number @xmath115 , @xmath189 when @xmath188 .
however , the convergence can be very slow for small @xmath115 , in which case , the derived asymptotic results may not be applicable for practical range of large but finite @xmath14 .
thus , for practical system applications , the sinr is only considered to be asymptotically deterministic when the scv decreases linearly with @xmath14 or faster . with this definition ,
we actually derive sufficient conditions for the sinr to be asymptotically deterministic .
[ pro1]when @xmath190 , for perfect pce , a sufficient condition for the sinr in to be asymptotically deterministic is @xmath191 for imperfect pce , the sinr is always asymptotically deterministic
.    please see appendix [ appendix b ] .
note that the above sufficient conditions are based on the scaling law in and .
they can be utilized to justify whether the cases discussed in d1-d5 ( including the power scaling analysis in existing works @xcite ) satisfy the condition for asymptotically deterministic sinr . for perfect pce , implies @xmath192 , meaning that to make the sinr asymptotically deterministic , the highest possible sinr scaling is @xmath193 . without the constraint for sinr to be asymptotically deterministic ,
the maximum possible value for @xmath122 is @xmath1 . for imperfect pce ,
proposition [ pro1 ] implies that no new limit is casted to the sinr scaling for the sinr to be asymptotically deterministic .
next , typical scenarios with asymptotically deterministic sinr are investigated in which @xmath119 is allowed to take values from @xmath194 only .
the tradeoff between parameters will be revealed .    *
[ d7 ] for perfect pce , to achieve both @xmath195 ( the sinr increasing linearly with @xmath193 ) and asymptotically deterministic sinr , the sufficient condition reduces to @xmath196 , @xmath197 .
it means that the user number should be constant , the product of training energy and the transmission power must scale as @xmath164 where the tradeoff between the energy efficiency in training and transmission period can be found .
two typical cases with constant user number may happen : a ) @xmath198 ; and b ) @xmath199 .
for case a ) , when the training energy is constant , the transmission power should scales as @xmath164 . for case
b ) , when the training energy scales as @xmath164 , the transmission power should be constant . *
[ d8 ] for imperfect pce , to achieve both the highest scaling @xmath200 and asymptotically deterministic sinr , the sufficient condition is the same as the condition for non - decreasing sinr in , i.e. , @xmath201 .
it means that the product of the training energy and the per - user transmission power should scale as @xmath202 or higher .
typical cases are referred to the former discussions in section [ dis_sec_1 ] . *
[ d9 ] for perfect pce , to achieve constant and asymptotically deterministic sinr , the sufficient condition reduces to @xmath203 .
it means that the product of training energy and the per - user transmission power must scale as @xmath202 .
three typical cases may happen : a ) @xmath204 ; b ) @xmath205 and c ) @xmath206 . for case
a ) , when the training energy is constant , the per - user transmission power should scales as @xmath202 . for case
b ) , when the training energy scales as @xmath164 , the per - user transmission power should also scales as @xmath164 . for case
c ) , when the training energy scales as @xmath202 , the per - user transmission power should be constant .    in existing works ,
only constant sinr case ( @xmath207 ) has been considered @xcite .
discussion d6-d8 show that for perfect pce , the asymptotically deterministic sinr can scale as @xmath193 . for imperfect pce , at most constant asymptotically deterministic sinr can be achieved without any extra conditions besides the one for non - decreasing sinr in .
in this section , main analytical results in this paper will be verified by simulations .
we assume that the total number of cells is @xmath208 , the spatial correlation metric is @xmath209 and the inter - cell normalized large scale fading is @xmath210 .
the matrix @xmath3 in and ( [ eq5 ] ) is assumed to consist of the first @xmath211 columns of the @xmath14 by @xmath14 discrete fourier transfer ( dft ) matrix . due to the main focus of this paper ,
we only consider systems with non - decreasing sinr scenario as given in corollary [ corollary 1 ] .
first , we consider the case of perfect pce , i.e. , @xmath212 for all simulated antenna numbers .
eleven typical parameter settings as specified in table i are considered .
[ table 1 ]    .
parameter settings for perfect pce [ cols="^,^,^,^,^,^,^",options="header " , ]         in fig .
4 , the simulated effective sinr with respect to @xmath14 is shown for the network settings given in table ii . for the constant pilot contamination @xmath183 case , the limiting effect of pilot contamination on the performance scaling is shown in cases 4 and 5 .
note that these two cases do not match the scaling law as well as that of perfect pce since the pilot contamination term is not dominant enough for the antenna number interval in both cases .
however , for the simulated finite antenna numbers , the scaling law for case 4 is more reliable compared with that for case 5 .
this is because that the former case has a better compliance with the condition in corollary [ corollary 3 ] .
further , as @xmath14 increases , the accuracy of the scaling law for both cases improves . on the other hand ,
cases 1 - 3 show that for decreasing pilot contamination with respective to @xmath14 , an increasing sinr can be obtained ( case 3 ) , which is in accordance with remark [ remark 6_170211 ] .        in fig .
5 , it can been shown that the lower bound given in lemma [ lemma 2 ] also match the simulated achievable ergodic sum - rate well for all parameter settings in table ii .
for tdd mutli - cell multi - user massive mimo downlink systems with csi error , channel spatial correlation and mrt precoding , this paper provided a comprehensive general scaling law based performance analysis , which quantitatively shows the effect of important network parameters and the pilot contamination elimination and the tradeoff among them .
the application of the scaling law for large but finite antenna number was also studied .
moveover , a sufficient condition for the sinr to be asymptotically deterministic in the sense of mean square convergence was derived .
our results cover existing ones on such analysis as special cases and show the effect of pilot contamination more explicitly rather than only considering it as the performance limit in existing works . in the future work , more general spatial correlation distribution and more practical modeling on the effect of pilot contamination elimination
should be studied .
due to the relationship in , the key step for calculating the means and scvs of the sinr components is the calculation of @xmath213 , @xmath214 , and @xmath215 , @xmath216 .
the following results are derived for all @xmath217 and @xmath218 .    first , with some tedious calculations which are omitted here for space limit , we have @xmath219 @xmath220 @xmath221 and @xmath222 where - are derived from properties of gamma distribution and in , the central limit theorem is used to approximate the distribution of @xmath223 as @xmath224 for all @xmath225 based on . based on ( [ eq38])-([eq41 ] ) and the relationship in ( [ eq11 ] ) , the next step for calculating the means and
scvs of the sinr components is to follow their definitions and find the highest order terms with respect to @xmath14 .
details are omitted here due to space limit .
the sinr expression in ( [ eq17 ] ) can be reformulated as @xmath226 the sinr is asymptotically deterministic when its scv approaches zero as @xmath227 .
however , the complex structure of the sinr expression makes it challenging to obtain its scv directly . alternatively , since @xmath96 is asymptotically deterministic as noted in remark [ remark 1 ] , for the sinr to be asymptotically deterministic , the sufficient and necessary condition is that the denominator in ( [ eq42 ] ) is asymptotically deterministic .
one sufficient condition is the scv of the denominator in ( [ eq42 ] ) scales no larger than @xmath202 .
when @xmath228 is asymptotically deterministic , @xmath228 asymptotically converges to @xmath229 as defined in ( [ eq27 ] ) .
thus from the definition in ( [ eq28 ] ) , i.e. , @xmath230 , and @xmath231 ( refer to lemma [ lemma1 ] ) , we have @xmath232 therefore , the sufficient condition can be rewritten as @xmath233 for some constant @xmath234 .    by using the fact that a sufficient condition for ( [ eq45 ] ) is that the variance of each term in ( [ eq45 ] ) scales no larger than @xmath202 ( details are omitted for space limit ) and combining the scaling law in ( [ eq30 ] ) and ( [ eq31 ] ) ,
the sufficient condition in proposition [ pro1 ] can be obtained via some tedious calculations based on ( [ eq20])-([eq22 ] ) .
99 t. l. marzetta ,  noncooperative cellular wireless with unlimited numbers of base station antennas , " in _ ieee trans .
wireless commun .
_ , vol . 9 , no . 11 , pp . 3590 - 3600 , nov
e. bjrnson , e. g. larsson and t. l. marzetta ,  massive mimo : ten myths and one critical question , " in _ ieee commun .
114 - 123 , feb . 2016 .
f. rusek et al . ,  scaling up mimo : opportunities and challenges with very large arrays , " in _ ieee signal process . mag .
_ , vol . 30 , no .
40 - 60 , jan . 2013 .
s. wagner , r. couillet , m. debbah et al . ,
 large system analysis of linear precoding in correlated miso broadcast channels under limited feedback , " in _ ieee trans .
inf . theory _ , vol .
58 , no . 7 , pp . 4509 - 4537 , july 2012 .
j. hoydis , s. ten brink and m. debbah ,  massive mimo in the ul / dl of cellular networks : how many antennas do we need ? , " in _
ieee j. sel .
areas commun .
160 - 171 , feb . 2013 .
h. yang and t. l. marzetta ,  performance of conjugate and zero - forcing beamforming in large - scale antenna systems , " in _ ieee j. sel .
areas commun .
172 - 179 , feb . 2013 .
h. q. ngo , e. g. larsson and t. l. marzetta ,  energy and spectral efficiency of very large multiuser mimo systems , " in _ ieee trans .
1436 - 1449 , apr . 2013 .
q. zhang , s. jin , k. k. wong et al . ,  power scaling of uplink massive mimo systems with arbitrary - rank channel means , " in _ ieee j. sel .
topics signal process . _ , vol . 8 , no .
966 - 981 , oct . 2014 .
h. q. ngo , e. g. larsson and t. l. marzetta ,  massive mu - mimo downlink tdd systems with linear precoding and downlink pilots , " _
2013 51st annual allerton conference on commun . , control , and comput .
( allerton ) , monticello , il , _ 2013 , pp .
a. khansefid and h. minn ,  achievable downlink rates of mrc and zf precoders in massive mimo with uplink and downlink pilot contamination , " in _ ieee trans .
4849 - 4864 , dec . 2015 .
j. zuo , j. zhang , c. yuen et al .
 multicell multiuser massive mimo transmission with downlink training and pilot contamination precoding , " in _ ieee trans . veh .
65 , no . 8 , pp . 6301 - 6314 , aug
a. papazafeiropoulos ,  impact of general channel aging conditions on the downlink performance of massive mimo , " in _ ieee trans
1428 - 1442 , feb .
e. bjrnson , m. matthaiou and m. debbah .
 massive mimo with non - ideal arbitrary arrays : hardware scaling laws and circuit - aware design , " in _ ieee trans .
wireless commun .
4353 - 4368 , aug .
h. falconet , l. sanguinetti , a. kammoun et al .
 asymptotic analysis of downlink miso systems over rician fading channels , " _ 2016 ieee international conference on acoust . , speech and signal process .
( icassp ) , shanghai , _ 2016 , pp .
3926 - 3930 . p.
j. smith , c. neil , m. shafi et al .
,  on the convergence of massive mimo systems , " _ 2014 ieee international conference on commun .
( icc ) , sydney _ , nsw , 2014 , pp . 5191 - 5196
. q. wang and y. jing .
 performance analysis and scaling law of mrc / mrt relaying with csi error in massive mimo systems . "
arxiv preprint arxiv:1606.07480 ( 2016 ) .
h. yin , d. gesbert , m. filippou et al .
,  a coordinated approach to channel estimation in large - scale multiple - antenna systems , " in
_ ieee j. sel .
areas commun .
264 - 273 , feb . 2013
p. zhao , z. wang , c. qian et al .
 location - aware pilot assignment for massive mimo systems in heterogeneous networks , " in _ ieee trans .
65 , no . 8 , pp . 6815 - 6821 , aug
b. tomasi and m. guillaud ,  pilot length optimization for spatially correlated multi - user mimo channel estimation , " _
2015 49th asilomar conference on signals , syst . and comput . ,
pacific grove , ca , _ 2015 , pp .
1237 - 1241 .
o. elijah , c. y. leow , t. a. rahman et al .
 a comprehensive survey of pilot contamination in massive mimo - 5 g system , " in _ commun .
surveys tuts .
905 - 923 , second quarter 2016 .
j. jose , a. ashikhmin , t. l. marzetta et al . ,  pilot contamination and precoding in multi - cell tdd systems , " in _ ieee trans .
wireless commun .
10 , no . 8 , pp .
2640 - 2651 , aug . 2011 .
 fundamentals of statistical signal processing , volume i : estimation theory . " ( 1993 ) .
y. wang , s. li , c. li et al .
 ergodic rate analysis for massive mimo relay systems with multi - pair users under imperfect csi , " in _ 2015 ieee global conf . on sig . and
info . process .
( globalsip ) , orlando , fl , _ 2015 , pp .
t. v. t. le and y. h. kim ,  power and spectral efficiency of multi - pair massive antenna relaying systems with zero - forcing relay beamforming , " in _ ieee commun . letters _ , vol .
2 , pp . 243 - 246 , feb .
z. zhang , z. chen , m. shen et al .
spectral and energy efficiency of multipair two - way full - duplex relay systems with massive mimo , " in _ ieee j. sel .
areas in commun .
848 - 863 , apr . 2016 .
b. everitt , the cambridge dictionary of statistics , cambridge university press , 1998 .
s. boyd and l. vandenberghe , convex optimization , cambridge univ .
press , 2004 .
h. kobayashi , b. l. mark , and w. turin , probability , random processes , and statistical analysis , cambridge univ .
press , 2012 .