given some radiation field , a prediction can be made about the outcome of an experiment where a suitable detector is placed at some position @xmath0 and probes the radiation field for some time @xmath1 .
even when the radiation field is known as precisely as possible , there is still a finite uncertainty in the outcome .
this statement can also be formulated in the inverse way : it is not possible to build an apparatus that provides a `` better '' ( = more defined ) illumination at the point @xmath0 than this limit .
we will demonstrate in this paper that the uncertainty of the photocount , quantified by its variance , depends on both the number of modes , including the energy distribution among the modes , and on the spectral properties and quantum correlations of the radiation field .
another contribution to the uncertainty will appear if the radiation field is known at the source but the mode structure is then perturbed in some uncontrolled way before it reaches the location of the detector .
the latter effect is frequently referred to as static speckle  @xcite named after the pattern it creates on a camera .
the opposing term `` dynamic speckle '' unfortunately is used to denote two completely different concepts in the literature .
frequently it refers to a static speckle pattern that is changing over time because the source of the perturbation is moving  @xcite but in this paper we refer to `` dynamic laser speckle ''  @xcite  an effect that is due to the quantum dynamics of the radiation field .
quantum effects are most prominent on small length scales .
modern optical lithography operates precisely in this regime , printing structures smaller than the wavelength of the light with a precision of a few nanometres  @xcite . to this end
, photo resists have been developed that are more sensitive to changes of the light intensity than most technical sensors , thereby acting as efficient ( albeit unintentional ) detectors for uncertainties of the photocount . later in this paper
we will show that the effects of multi - mode multi - frequency radiation are most pronounced neither at the limit of large or short times but in the intermediate regime .
recent studies have confirmed that this is precisely the regime used in modern optical lithography  @xcite .
classically a radiation field is described by its electric field @xmath2 as function of position @xmath3 and time @xmath4 . in a quantum treatment , the electric field @xmath5
is replaced by a suitable operator . in both treatments ,
the electric field is not a perfectly determined quantity , and the maximum knowledge possible is contained in the mutual - coherence function  @xcite @xmath6 we restrict ourselves to stationary fields such that @xmath7 does not depend on @xmath8 and @xmath9 but only on the time - difference @xmath10 .
furthermore , for ease of writing we will only treat one component @xmath11 of the electric field @xmath5 but the extension to cover the other components is straight - forward .
the fourier transform @xmath12 of @xmath13 with respect to the time difference is a nonnegative hermitian operator , and thus possesses an eigenrepresentation  @xcite @xmath14 the eigenvectors @xmath15 are called the modes of the electric field and form a complete orthonormal set .
we assume that the frequency spectrum is small enough such that @xmath15 is independent of @xmath16 .
whenever the frequency spectrum is equal to the natural linewidth of the light source , i.e. , the spectrum is due to the finite lifetime of some excited light - emitting medium such as in a laser , this assumption is always fulfilled . for other light sources ,
this assumption can sometimes be problematic for cavity - like systems but for open system or in a waveguide geometry , this is less of an issue .
if the condition that @xmath15 is independent of @xmath16 , should be violated , one can split the mode into several discrete modes ( one for each frequency interval in which the shape of the mode can assumed to be constant ) and the remainder of this paper be applied nonetheless . the electric field
can then be written as @xmath17 since @xmath18 is independent of @xmath4 , the @xmath18 are the modes instead of being some arbitrary base of the electric field .
comparison of eqs .
( [ eqgamma2 ] ) and  ( [ eqgamma3 ] ) shows that the magnitude of @xmath19 can be computed from @xmath20 but not its phase .
a semiclassical treatment of the photocount statistics would assume that @xmath20 is fluctuating in time whereas within quantum optics , which we will apply in this paper , fluctuations are inherent to the operator description .
this paper is organised as follows . in sec .
[ seczeitunschaerfe ] we compute the variance of the photocount on a quantum - optical level when the radiation field is completely known .
we will find that there is a shot noise term , a quantum correlation term and a term describing the beating of different frequencies .
frequently , by design or unintentionally , the mode structure emitted by a known light source is completely changed and unknown when it reaches the point of the detector . as will be demonstrated in sec
.  [ secstaticspeckle ] , random - matrix theory allows a compact and exact treatment of this problem .
the resulting `` static speckle '' becomes the smaller the more modes of the radiation field are excited , and the more uniform the energy distribution among the modes is .
depending on the `` size '' of the radiation field , there is thus a minimum amount of static speckle , which is computed in sec .
[ secetendue ] . in sec .
[ secstatickorrelation ] we extend this question to computing the most likely amount of static speckle .
since we will demonstrate this quantity to be self - averaging , the computed average is more than just an average in that it is characteristic for almost all individual speckle values .
while in `` general '' electrodynamics the electric field @xmath2 can be probed directly , this is not possible in the realm of optics as @xmath2 changes too quickly in time and space to allow a direct measurement .
rather , the radiation field is probed by means of photodetection , i.e. , by absorbing photons inside some device and counting the number of photons absorbed .
this principle applies to technical machines as well as to the human eye or to photo resists .    on small length and time scales , quantum effects become important . for this reason , and
because all effects can then be treated in a more compact way , we will use the quantum theory of photodetection in the following .
this theory  @xcite is usually formulated for single - mode detection .
since multi - mode fields  @xcite are at the core of this paper , we will present an extension of the theory to multi - mode photodetection here .
we will treat the multi - frequency aspect explicitly and not express the different frequencies by different modes as is frequently done in textbooks .
we label the modes of the electromagnetic field by @xmath21 .
the annihilation operator associated with this mode is @xmath19 . using this notation ,
the quantum operator for the electric field at some position @xmath0 is given by  @xcite @xmath22 later it will prove helpful to switch from the time representation @xmath23 to the spectral representation @xmath24 , @xmath25 since in time representation , operators taken at different times do not necessarily commute whereas they do in the spectral representation , @xmath26=\delta_{nm}\delta(\omega-\omega ' )      \;.
\label{eqvertauschungsrelation}\ ] ] when we label the number of photons in the @xmath27-the mode with @xmath28 , which basically amounts to the total intensity of that mode , this gives the expectation value @xmath29 where @xmath30 is the ( normalised ) spectrum of the radiation in the @xmath27-th mode , and the prefactor @xmath31 has been introduced for later convenience .
photodetection within some time interval @xmath1 is then described by the quantity  @xcite @xmath32 where @xmath33 marks the detection efficiency and includes information on the size of the detector .
this amounts to a perturbative description of the interaction between the detector and the radiation field .
the perturbative approach neglects the effect that every detected photon decreases the number of photons remaining in the radiation field  @xcite but this is mainly an issue for microcavities where only a small number of photons are excited at one point in time and the dynamics are then studied .
for the purpose of this paper , this is no relevant restriction .
the factorial moment @xmath34 of the photodetection count distribution is given by @xmath35 where the colons denote normal - ordering of the operator within , and the brackets @xmath36 denote the average which has to be taken over both the quantum fluctuations of the operators @xmath19 and the frequency distribution @xmath30 .
mean and variance of the photo count then follow from @xmath37 ^ 2\;.      \label{eqdefvar}\ ] ]    for the first moment , this yields the result @xmath38 this is the same result as expected by classical theory since @xmath39 corresponds to the intensity of the @xmath27-th mode at position @xmath0 .
the second factorial moment is    @xmath40 }          \;.\end{gathered}\ ] ]    taking the average over the quantum mechanical expectation operators gives nonzero contributions for three cases of the indices @xmath27 : @xmath41 , @xmath42 , and @xmath43 , thus @xmath44 ^ 2          [ a_{n}(\omega)]^2          \rangle   \\          + \frac{\eta^2}{4\pi^2}\int\!\!\!\int_0^t\mathrm{d}t_{1,2 }          \int\mathrm{d}\omega_{1,3 } \sum_{n_1\ne n_3 }           { \lvert\phi_{n_1}(\vec{r})\rvert}^2{\lvert\phi_{n_3}(\vec{r})\rvert}^2           \langle a^{\dagger}_{n_1}(\omega_1 )          a_{n_1}(\omega_1)a^{\dagger}_{n_3}(\omega_3)a_{n_3}(\omega_3 )          \rangle            \\      + \frac{\eta^2}{4\pi^2}\int\!\!\!\int \mathrm{d}\omega_{1,2 }          \int\!\!\!\int_0^t\mathrm{d}t_{1,2 }          \sum_{n_1\ne n_2 }           { \lvert\phi_{n_1}(\vec{r})\rvert}^2{\lvert\phi_{n_2}(\vec{r})\rvert}^2                  \langle a^{\dagger}_{n_1}(\omega_1 )          a_{n_1}(\omega_1)a^{\dagger}_{n_2}(\omega_2)a_{n_2}(\omega_2 )          \rangle           \mathrm{e}^{\mathrm{i}\bigl[(\omega_1-\omega_2)t_1          + ( \omega_2-\omega_1)t_2\bigr ] }       \;.      \label{eqlang1}\end{gathered}\ ] ]    this equation can be simplified by inserting the expectation values from eq .
( [ eqerwart1 ] ) .
the double integral over @xmath8 and @xmath9 in the last line of eq .
( [ eqlang1 ] ) can be reduced to a single integral via @xmath45\;.      \label{eqttransform}\ ] ] the remaining two integrations over @xmath46 and @xmath47 in that line each yield the fourier transform @xmath48 of @xmath49 , @xmath50 this thus gives @xmath51 ^ 2          [ a_{n}(\omega)]^2          \rangle   \\
+ \eta^2 t^2   \sum_{n_1\ne n_3 }          { \lvert\phi_{n_1}(\vec{r})\rvert}^2{\lvert\phi_{n_3}(\vec{r})\rvert}^2          i_{n_1 } i_{n_3 } \\
+ \frac{\eta^2}{2\pi }          \int_{-t}^t\mathrm{d}t          \sum_{n_1\ne n_2 }           { \lvert\phi_{n_1}(\vec{r})\rvert}^2{\lvert\phi_{n_2}(\vec{r})\rvert}^2              \\          \times              i_{n_1 } i_{n_2 } \hat{g}_{n_1}(t ) \hat{g}^*_{n_2}(t )           [ t-{\lvertt\rvert } ]      \;.      \label{eqlang4}\end{gathered}\ ] ] the summation over @xmath52 is rather inconvenient .
thus , we include the missing terms @xmath53 in that summation , yielding the value @xmath54 ^ 2 $ ] , and correct for this by subtracting in the first summation the terms just added .
there now appears the quantity @xmath55 ^ 2          [ a_{n}(\omega)]^2          \rangle            - \langle a^{\dagger}_{n}(\omega ) a_{n}(\omega ) \rangle^2           \ ; ,           \label{eqdeff}\ ] ] that quantifies the deviation of the radiation in the @xmath27-th mode from a coherent state .
for coherent radiation , @xmath56 since then @xmath57 ^ 2 [ a_{n}(\omega)]^2 \rangle=\langle a^{\dagger}_{n}(\omega ) a_{n}(\omega ) \rangle^2 $ ] whereas for any gaussian light source , in particular any object emitting thermal radiation , @xmath57 ^ 2 [ a_{n}(\omega)]^2 \rangle=2 \langle a^{\dagger}_{n}(\omega ) a_{n}(\omega ) \rangle^2 $ ] , and @xmath58 . for classical radiation , @xmath59 whereas for certain nonclassical radiation @xmath60 is possible .
collecting results , using eq .
( [ eqdefvar ] ) , gives the variance of the photocount @xmath61      \;.
\label{eqlang3}\end{gathered}\ ] ]    equation  ( [ eqlang3 ] ) constitutes a core result of this paper .
it gives the general expression for the noise that is intrinsic to any measurement of a multi - mode multi - frequency radiation field over some time @xmath1 .
the first term is the well - known shot noise term that appears naturally in the quantum treatment : electric field and photocount are quantum - optically related via the factorial moment , see eq .
( [ eqfactorialmoment1 ] ) , whereas in classical optics the relation is via the `` plain '' moment , with the difference between factorial and plain moments being precisely the shot - noise term .    the second term , quantified by @xmath62 from eq .
( [ eqdeff ] ) , marks the excess noise when the radiation field is not in a coherent state .
as mentioned above , @xmath56 for coherent radiation , and @xmath63 for thermal radiation .
even when @xmath62 is small , any nonzero @xmath62 will eventually become the dominating effect for the photonoise as its contribution increases quadratically with time , in contrast to all other contributions .
the cross - over point , where this term becomes larger than the shot - noise term , is easily estimated from eq .
( [ eqlang3 ] ) as @xmath64 . for any thermal radiation
this regime thus is already entered once the measurement time is long enough for one photon per mode to be counted on average .
@xmath28 can be computed from the bose - einstein factor , evaluated at the temperature @xmath65 of the source , @xmath66 this value depends on the ratio of frequency and temperature , and thus has the same value @xmath67 for both a light bulb emitting visible light and an euv plasma source operating at up to @xmath68 to emit radiation with @xmath69 .    the third and final contribution in eq .
( [ eqlang3 ] ) quantifies the beating of different frequency components when the measurement time is finite .
it is sometimes referred to as dynamic laser speckle since it is most relevant for measurement times @xmath1 not much larger than the coherence time @xmath70 of the radiation , @xmath71 hence the name `` temporal degree of coherence '' frequently given to @xmath48 . however , the factor @xmath72 $ ] is also relevant , and it would be wrong to replace the integral simply by @xmath73 , as we will now demonstrate .
the extra factor @xmath72 $ ] , which seems to be largely ignored in phenomenological literature , gives less weight to the frequently surprisingly large tails of @xmath48 .    two frequently encountered spectral distributions are gaussian and lorentzian , @xmath74 the curves above are already normalised to @xmath75 , i.e. , their integral ( not their square integral ) is @xmath75 .
the widths of the spectral distributions have been chosen such that the coherence time , computed from eq .
( [ eqdeftau ] ) , is the same value @xmath70 , allowing for a direct comparison of these two curves .
their fourier transforms are @xmath76   \;,\\
\hat{g}_{\mathrm{lorentz}}(t ) & = \frac{1}{\sqrt{2\pi } } \exp\bigl [          -\frac{{\lvertt\rvert}}{2\pi\tau } \bigr ] \;.
\label{eqkurvenformen}\end{aligned}\ ] ] the results for entering these expressions into the integral in eq .
( [ eqlang3 ] ) are shown in fig .
[ figkurvenformen ] .
there is a significant dependence on the spectral shape , with the longer tails of the lorentzian leading to a slower approach to saturation .
please note that the graph corresponds to the variance of the photocount .
when the square root of the variance is scaled by the mean intensity  this quantity is called the contrast  the saturating curve shown in the figure turns into a curve decreasing to zero as @xmath1 becomes larger .
@xmath77 is a plain exponential and thus conforms best to the simple model of how to model finite temporal coherence .
we thus have used this curve for contrasting eq .
( [ eqlang3 ] ) by an approximation where the factor @xmath78 is replaced by @xmath1 , thereby ignoring the additional dependence on @xmath79 . for all values of @xmath1 ,
the integral increases by about @xmath80 , making the factor @xmath78 important for correct modelling of the photon counting statistics .    ) that quantifies the strength of the dynamic speckle .
the two curves at the bottom follow by inserting the expressions for a gaussian and a lorentzian spectrum from eq .
( [ eqkurvenformen ] ) . when the term @xmath81 in eq .
( [ eqlang3 ] ) is replaced by @xmath1 , and an exponentially decreasing @xmath48 is assumed , the third curve follows . ]
a final word about the cross - over from `` dynamics '' to `` statics '' seems to be in order .
since the temporal effect begins to saturate at time @xmath70 it is frequently assumed that the nontemporal regime is entered once @xmath82 .
while this statement is correct for single - mode radiation fields , it is incorrect here since the shot - noise term @xmath83 , setting the reference value , scales linearly with the number @xmath84 of modes whereas the prefactor of the temporal term scales as @xmath85 such that the condition @xmath82 has to be replaced by @xmath86 .
incidentally , since it is the regime @xmath87 that current high - power excimer lasers are operating in  @xcite , this distinction is of actual importance .
similar to classical mechanics , theoretical optics knows the concept of phase space  @xcite .
the optical phase space is spanned by position and spatial frequency , informally called @xmath88-vector , and the radiation field is completely described by the ( pseudo)-density function @xmath89 known as the wigner function . any radiation field then occupies a certain volume in phase space . unfortunately , there is no good mathematical metric to actually measure or at least define such a volume .    in experimental and technological areas of optics , in contrast , the volume of phase space occupied is well - defined and is referred to as etendue .
the energy density @xmath90 of the light , given as a function of position and angle , is assumed to be larger than zero only inside a finite and well - defined region .
this allows for an easy and direct definition of the etendue of the radiation field .
etendue is an important quantity since , at least within geometrical optics , it can only be changed by vignetting the radiation field .
the etendue supplied by some radiation field incident on some optical apparatus should thus be no larger than the etendue that this apparatus can accept .    convolving the wigner function @xmath89 with a small kernel , its size given by an uncertainty - relation  @xcite , yields @xmath90 when the relation @xmath91 between angle and @xmath88-vector is utilised .
thus on first sight it might seem that the relation between microscopic and macroscopic definition would be straight - forward .
the essential difference is that @xmath89 has no finite support whereas @xmath90 is assumed to have .
it can be shown that due to the uncertainty principle , every mode @xmath18 needs to have infinite tails in at least either real space ( position ) or @xmath88-space ( angle )  @xcite , and this is reflected in @xmath89 .
of course , also @xmath90 has the same infinite tails but in the macroscopic world , the function values inside those tails drop down so quickly that they are `` assumed '' to be zero once they are so low that they have become irrelevant for practical purposes .
we thus need a definition for phase space volume that is based on a microscopic radiation field but reduces to the macroscopic etendue for `` large '' radiation fields . for the purpose of this paper ,
the relevant way of quantifying phase space volume is by the number of modes that can be fit into it .
we will tackle the problem in the opposite way : we pick a certain number of `` sensible '' base functions and compute all possible radiation fields that can be formed from them .
we will then identify for the intensity profiles @xmath92 as function of position and @xmath93 as a function of direction the properties that allow to deduce the number of base functions used .    to have shorter mathematical expressions , we will express directions in terms of the @xmath88-vector instead of the angle @xmath94 .
furthermore we make two assumptions .
first , we restrict ourselves to one spatial dimension but the extension is straight - forward and becomes trivial when the radiation field factorises in its spatial dimensions .
second , we assume that the intensity distribution @xmath95 factorises into a spatial term @xmath96 and an angular term @xmath97 . since the etendue accepted by practically any real - life apparatus factorises , these assumptions do not limit the applicability of this approach .
we start with @xmath98 orthonormal base functions @xmath99 ,
the orthogonality is essential since this prevents the base functions from being too similar and thus collapsing into a very small volume of phase space .
furthermore , all modes of the radiation field have to be mutually orthogonal anyhow as they are eigenfunctions of a hermitian operator , cf .  eq .  ( [ eqgamma2 ] )
we assume them to be normalised for the ease of the calculation .
for every @xmath99 , there also exists its fourier transform @xmath101    to progress , we need to choose a certain set of base functions .
the aim is to pack the base functions as tightly as possible in phase space , and without loss of generality we try to pack them around the origin .
if all @xmath18 , @xmath100 , are assumed to fulfil the conditions    [ eqhardy1 ] @xmath102    then it can be shown  @xcite that @xmath103 is the most condensed situation for which there are solutions .
furthermore , all solutions can then be written as @xmath104 where @xmath105 , @xmath106 are polynomials of order not higher than @xmath84 .
only the product @xmath107 is fixed by this but not @xmath108 or @xmath109 on their own .
this freedom is equivalent to changing @xmath110 with the simultaneous change @xmath111 .
this also reflects that , while phase space volume is well - defined and preserved as volume , its projection onto a particular axis is allowed to change .
the set of solutions @xmath99 of eq .
( [ eqhardy2 ] ) is @xmath98 dimensional , and a convenient orthonormal set consists of the hermite functions @xmath112 , defined as @xmath113 h_n(x / r)\ ; ,      \label{eqdefhermite1}\ ] ] where @xmath114 are the hermite polynomials , @xmath115 the prefactor in eq .
( [ eqdefhermite1 ] ) ensures normalisation of @xmath112 .
the hermite functions are eigenfunctions of the fourier transform , @xmath116 such that the derivation presented in the following for real space also applies to @xmath88-space .
all normalised functions @xmath117 inside the solution space of eq .
( [ eqhardy1 ] ) can be written as @xmath118 the maximum intensity @xmath119 possible for such a state can be computed using the method of lagrange multipliers . finding an extremum of @xmath117 is equivalent to finding one of @xmath119 , and we will pick either of these quantities depending on which one is more convenient .
the condition @xmath120 of the function @xmath121 yields the intermediary result @xmath122 . by multiplying this with @xmath123 respectively @xmath124 and summing
, this gives the two conditions    [ eqlagrange3 ] @xmath125    solving for @xmath126 by eliminating @xmath117 and inserting @xmath126 into eq .
( [ eqlagrange3b ] ) , remembering that @xmath127 , one arrives at the result for the maximum intensity , namely @xmath128 the method of lagrange multipliers only gives necessary but not sufficient conditions for extrema , meaning that there can not be any additional extrema different from eq .
( [ eqlagrange4 ] ) but this equation might not describe an extremum in the first place .
equation  ( [ eqlagrange4 ] ) actually describes two solutions , @xmath129 and @xmath130 .
it is obvious from physics reasons that at least two extrema , one being the most positive and one the most negative , exist .
hence , eq .  ( [ eqlagrange4 ] ) uniquely describes these two extrema .     and @xmath131 , computed from eq .  ( [ eqlagrange4 ] ) , as a function of @xmath84 where @xmath84 is the number of basis functions minus 1 .
depicted are the curves for @xmath132 .
if the corresponding curves for a normalised radiation field are below a given line , the radiation field can not have more modes than indicated by the number @xmath84 on the curve . ]    a graphical display of the solutions @xmath133 and @xmath134 can be found in fig .
[ figlagrange1 ] as a function of @xmath84 .
the right edge of @xmath133 becomes progressively steeper as @xmath84 increases .
this means that in the macroscopic limit @xmath135 there exist a sharp edge that can be associated with the macroscopic etendue . to define the position of this edge
, one could pick the @xmath136-coordinate where @xmath137 or where it is equal some other particular value but this would lead to an arbitrary result .
it can be shown that the number of basis functions is then proportional to the square of the width at which the `` cut '' is done  @xcite such that the result depends more on this arbitrary decision than on the radiation field .
also the method of `` almost bandwidth - limited functions '' ( which is usually formulated for time - dependent electric signals , hence the term `` bandwidth '' instead of `` angular range '' ) suffers from a similar problem  @xcite .
luckily , there exists a strong link between the microscopic and macroscopic worlds because the hermite functions @xmath112 are the solutions of the eigenvalue equation @xmath138 that can be found in any quantum mechanics textbook as it describes the harmonic oscillator .
the eigenfunction @xmath139 has eigenvalue ( = energy ) of @xmath140 , and a classical oscillator with the same energy is restricted to the interval @xmath141 .
the extremal curve @xmath117 computed above thus is the extremal superposition of harmonic - oscillator solutions with energy level @xmath27 up to @xmath84 .
hence , the classical edge beyond which the oscillator can not be found is given by @xmath142 .
summarising , a microscopic radiation field has at most @xmath98 modes when for all of its modes @xmath143 the conditions    @xmath144    are fulfilled for some value @xmath145 identical for all modes .
since @xmath146 and @xmath147 are the intensities as a function of position and angle , respectively , these are quantities that can in principle be measured .
a macroscopic radiation field with a half - diameter @xmath148 in real space and a half - diameter @xmath149 in angle space is thus spanned by modes @xmath139 with @xmath100 determined by @xmath150 the scale factor @xmath145 drops out when expressing results in terms of the macroscopic etendue . for given etendue @xmath151 ,
the number @xmath84 of modes ( minus 1 ) is thus @xmath152 if the etendue is specified in angle ( radians ) instead of a @xmath88-vector , this becomes @xmath153
the variance of the photocount computed in sec .  [ seczeitunschaerfe ] describes the uncertainty of the photocount if the radiation field at the detector is completely known .
frequently , the mode structure is known at the light source but not at the detector .
the standard undergraduate example is a hene - laser pointed at the wall of the lecture room where the observer can see a speckle pattern . when the eyes are then moved , the direction of motion of the observed speckle pattern depends on whether the viewer is near - sighted or far - sighted  @xcite .
the hene - laser emits a very regular mode structure but the modes are then changed by scattering at the rough wall .
the influence of the eye movement demonstrates that scattering at the wall does not result in local intensity changes at the place of scattering but rather in a deformation of the mode structure that translates into intensity changes only by propagation to and focusing in the eye .
the change of the mode structure on its way from the radiation source to the detector can be unintentional , just as in the example above , or it can be intentional due to some apparatus designed to mix the modes of the radiation field . unless optical systems are manufactured to the highest standards , they will have surface roughness in excess of the wavelengths .
this results in uncontrolled  thus in a certain sense random
 mode deformations .
in addition , any apparatus designed to provide a uniform and stable illumination at its exit can only do this by means of `` mixing '' the incoming light .
illumination systems used in optical lithography are the supreme example of such an apparatus as the local intensity at the output varies only by about @xmath154 even when the form of the illumination of its entrance is completely changed . while the light paths inside such an apparatus might be very controlled , for an outside observer there ( intentionally ) is no recognisable connection between the light at the entrance and the exit . in other words ,
an apparently random change of the mode structure occurs in this case as well .    in sec .
[ seczeitunschaerfe ] we have used the modes @xmath21 of the radiation field to expand the operator for the electric field , cf .
( [ eqfexpansion ] ) .
any orthonormal basis could have been used instead but using the modes offered the advantage that the associated annihilation operators @xmath24 become uncorrelated , cf .
( [ eqerwart1 ] ) .
we will now switch to a different basis set , at the expense of replacing eq .
( [ eqerwart1 ] ) by something more complicated .    ,
one switches from the modes @xmath18 that are optimal to describe the light source , to new base functions @xmath155 that are optimal to describe the photodetector . along with the transformation @xmath156
, the annihilation operators also need to be transformed , @xmath157 .
the additional action of scattering or mixing is then described by the matrix @xmath158 , transforming @xmath159 into @xmath109 . ]    considering only the radiation source and the detector , hence ignoring the scattering between them for the moment , we pick new basis functions @xmath160 and associated annihilation operators @xmath161 such that at the location @xmath0 of the detector all @xmath155 except for @xmath162 vanish . such a basis set can always be found and offers the advantage that photodetection at the point @xmath0 has been reduced to the photocount of the operator @xmath163 .
since both bases are orthonormal , the transformation between @xmath164 and @xmath165 can be described by a unitary matrix @xmath166 , @xmath167 with @xmath168 .
computing @xmath166 explicitly would be a very formidable task but it will turn out that knowledge of @xmath166 is not needed .
we will fix the number @xmath169 of new basis functions @xmath155 shortly but need to allow for the case that @xmath169 is larger than the number @xmath84 of excited modes of the radiation field .
this is easily achieved by simply adding @xmath170 vacuum states to the input for eq .
( [ secsmatrix ] ) .
we now assume that the radiation field is perturbed between the light source and the place of the detector by means of some `` virtual device ''  either intentionally by a properly designed apparatus or , usually unintentionally , by quasi - random scattering .
this can be described in two opposite but physically equivalent approaches .
either , one assumes that the mode structure is perturbed , thus changing @xmath171 while keeping @xmath172 unchanged , or one assumes that the mode structure @xmath171 is unchanged while energy is transferred between modes , thus modifying @xmath172 .
we use the second approach , also known as the method of input - output relations  @xcite , since the annihilation operators @xmath173 describing the radiation leaving the device can then be expressed in terms of the annihilation operators @xmath172 entering device . in the absence of nonlinear media , this relation is linear , and can be described by a matrix @xmath158 @xmath174 if there is no absorption , @xmath158 is unitary .
apart from energy conservation , already the commutation relation  ( [ eqvertauschungsrelation ] ) demands the unitarity of @xmath175 .
in the presence of absorption , eq .
( [ secsmatrix2 ] ) would need to be supplemented by noise sources  @xcite to ensure these commutation relations .    in the absence of further information
, the unitary matrix @xmath166 describing the distortion of the mode structure is uniformly distributed in the space of unitary matrices  @xcite .
this concept of `` maximum uncertainty '' is known from other areas of physics as well , most notably from thermodynamics where the basic lemma is that every microstate compatible with the macroscopic boundary conditions is equally likely .
the photocount distribution at some point @xmath0 has thus been reduced to the photocount distribution for the mode @xmath176 , related to the original incident mode structure of the radiation source by @xmath177 since for a given position of detector and light source @xmath0 , @xmath166 is a fixed unitary matrix , the product @xmath178 is also distributed uniformly  @xcite , and so no knowledge of @xmath166 is needed , as promised above .    the only information about what happens between light source and detector enters in the form of the number @xmath169 of basis functions .
this number is not arbitrary but is determined by the physics of the mode deformation process described by the matrix @xmath158
. please remember that @xmath169 is the number of base functions that are mixed by the transformation @xmath158 .
if the mode deformation is due to some intentional mixing in an apparatus accepting and emitting some etendue @xmath11 , the number @xmath169 follows from the results of sec .
[ secetendue ] and has a finite , well - defined value .
if the mode deformation is due to scattering or a diffuser , @xmath169 is determined by the number of basis functions that could in principle end up in the detector due to scattering or diffusing . depending on maximum scattering angle
, this can be a very large number and is always larger than @xmath84 .
computing the first factorial moment analogue to eq .
( [ eqafrequenz2 ] ) yields @xmath179 where in the final equality sign all averages over the fluctuations of electromagnetic field have been taken , just as described in sec .
[ seczeitunschaerfe ] , and an average over the ensemble of random unitary matrices remains .
the average is easily computed from @xmath168 , hence the average in eq .
( [ eqafrequenz3 ] ) is equal to @xmath180 where @xmath169 is the order of the matrix , yielding @xmath181 this is the expected result , namely that the mean intensity at the detector is equal to the properly scaled total intensity of the original light field .
computing the first moment demonstrated that to translate the equations from sec .
[ seczeitunschaerfe ] , all terms of the form @xmath182 have to be replaced by @xmath183 , and all terms @xmath184 by @xmath185 .
the correct starting point is eq .
( [ eqlang4 ] ) , which transforms into @xmath186 ^ 2          [ a_{n}(\omega)]^2          \rangle   \\          + \eta^2 t^2   \sum_{n_1\ne n_3 }          \langle { \lverts_{1n_1}\rvert}^2 { \lverts_{1n_3}\rvert}^2 \rangle          i_{n_1 } i_{n_3 } \\
+ \frac{\eta^2}{2\pi }          \int_{-t}^t\mathrm{d}t          \sum_{n_1\ne n_2 }           \langle { \lverts_{1n_1}\rvert}^2 { \lverts_{1n_2}\rvert}^2 \rangle           \\          \times              i_{n_1 } i_{n_2 } \hat{g}_{n_1}(t ) \hat{g}^*_{n_2}(t )           [ t-{\lvertt\rvert } ]      \;.      \label{eqlang6}\end{gathered}\ ] ] it only remains to compute the necessary averages over the unitary group in a similar spirit to above where we found that @xmath187 .
the average of the square of this expression is frequently needed and can be found in many papers , e.g. in ref .  , with the result @xmath188 ^ 2 \rangle = \frac{2}{m(m+1 ) } \;.      \label{eqsquadratmittel}\ ] ] the average @xmath189 for @xmath190 is not identical to @xmath191 as the unitary condition @xmath192 introduces correlations among the elements of @xmath175 .
starting from the exact relation @xmath193 allows us to write @xmath194 which , together with eq .
( [ eqsquadratmittel ] ) , gives us the desired average , @xmath195 this allows us to transform eq .
( [ eqlang6 ] ) into @xmath196 ^ 2          [ a_{n}(\omega)]^2          \rangle   \\          + \frac{\eta^2 t^2}{m(m+1 ) }   \sum_{n_1\ne n_3 }          i_{n_1 } i_{n_3 } \\      +          \frac{\eta^2}{2\pi m(m+1 ) }          \int_{-t}^t\mathrm{d}t          \sum_{n_1\ne n_2 }               i_{n_1 } i_{n_2 } \hat{g}_{n_1}(t ) \hat{g}^*_{n_2}(t )           [ t-{\lvertt\rvert } ]      \;.      \label{eqtemp1exakt}\end{gathered}\ ] ] when computing the variance in sec
.  [ seczeitunschaerfe ] , the term with @xmath197 dropped out when @xmath54 ^ 2 $ ] was subtracted from @xmath198 .
this is no longer the case here , and one arrives at @xmath199      \;.
\label{eqtemp2exakt}\end{gathered}\ ] ]    before discussing this result , we will analyse an approximation that is frequently done : the exact averages from eq .  ( [ eqmittelwertes ] ) are replaced by their gaussian approximation , equivalent to a phasor representation  @xcite , @xmath200 the denominator in eq .
( [ eqmittelwertes ] ) was @xmath201 so that the mixing of the modes would conserve energy .
equation  ( [ eqmittelwertesnaeherung ] ) neglects this constraint . while this error might seem to be negligible for large @xmath169
, we will demonstrate that its effect is surprisingly large .
first , we compute @xmath202 ^ 2          [ a_{n}(\omega)]^2          \rangle
\\          + \frac{\eta^2 t^2}{m^2 }   \sum_{n_1\ne n_3 }          i_{n_1 } i_{n_3 } \\
+          \frac{\eta^2}{2\pi m^2 }          \int_{-t}^t\mathrm{d}t          \sum_{n_1\ne n_2 }               i_{n_1 } i_{n_2 } \hat{g}_{n_1}(t ) \hat{g}^*_{n_2}(t )           [ t-{\lvertt\rvert } ]      \;.\end{gathered}\ ] ] in contrast to eq .
( [ eqtemp1exakt ] ) , the term with @xmath203 now cancels against @xmath54 ^ 2 $ ] .
the variance then becomes @xmath204      \;.\end{gathered}\ ] ] this is basically eq .
( [ eqlang3 ] ) with an additional contribution @xmath205 that is called `` static '' as its scaling with @xmath206 implies that its strength , when scaled by the mean squared intensity , is independent of measurement time @xmath1 .
we arrived at the ensemble of random @xmath175 for calculating the variance @xmath207 by assuming a fixed position of the detector while the mode deformations between radiation source and detector is different for a every member of the ensemble .
however , from eq .
( [ eqsstatischgesamt ] ) it follows that this variance can equally well be computed or measured by keeping the mode formations constant while changing the position of the detector .
the static speckle contrast @xmath208 , @xmath209^{1/2}}{\sum_n i_n } \;,\ ] ] can thus equally well be determined by doing a ( long ) measurement at different detector positions @xmath0 . in practise ,
one uses a camera instead of a point - like photodetector , and @xmath208 is calculated from the observed contrast of the image taken by the camera .
the image resembles the speckle pattern found on the fur or skin of many animals , hence the name `` speckle '' .
the static speckle contrast can , in this approximation , conveniently be expressed in terms of the effective number of modes @xmath210 , @xmath211 ^ 2}{\sum_n i_n^2 } \;,\ ] ] that is equal to the actual number of modes if all intensities are equal , and smaller ( implying more static speckle ) otherwise .
for the exact result  ( [ eqtemp2exakt ] ) , static speckle is not so easy to define , as the difference between eq .
( [ eqtemp2exakt ] ) and the results from sec .
[ seczeitunschaerfe ] is not a simple term proportional to @xmath206 .
however , taking the limit @xmath212 , hence living up to the name `` static '' , the difference between these two terms gives @xmath213^{1/2 } } {          \overline{n } } \nonumber\\      & = \left [ \frac{m}{m+1 } \frac{\sum_n i_n^2}{\bigl[\sum_n i_n\bigr]^2 }       + \frac{m-1}{m+1 } \frac{\sum_n f_n i_n^2}{\bigl[\sum_n i_n\bigr]^2 }      -\frac{1}{m+1}\right]^{1/2 } \nonumber\\      & = \left [ \frac{m}{m+1 } \bigl(\frac{1 } { n_{\mathrm{eff } } }          -\frac{1}{m } \bigr )      + \frac{m-1}{m+1 } \frac{\sum_n f_n i_n^2}{\bigl[\sum_n i_n\bigr]^2 }      \right]^{1/2 } \;.      \label{eqcexakt}\end{aligned}\ ] ]    there thus is a very significant difference between the exact result and the approximation . in the latter ,
the contribution  ( [ eqvarnstatisch ] ) is minimal but still nonzero if all modes carry the same intensity @xmath28 whereas for the exact result eq .
( [ eqtemp2exakt ] ) , assuming coherent radiation , this contribution vanishes if all @xmath28 are identical since then @xmath214 .
the approximation becomes valid in the limit that there are many modes but only a small number of them are actually excited : equation  ( [ eqvarnstatisch ] ) depends on the average squared intensity @xmath215 of the modes of the light source , where for a fixed light source the average is taken over its modes .
this kind of average is indicated by the notation @xmath216 .
in contrast , the exact result eq .
( [ eqtemp2exakt ] ) depends on @xmath217 $ ] but the term @xmath218 scales with the square of the fraction of excited modes and may thus be neglected in the limit stated above .
this important conclusion can also be formulated in terms of the etendue introduced in sec .
[ secetendue ] .
the approximation is valid in the limit that the light source fills only a small fraction of the etendue that the `` device '' used for mixing can accept .
it also means that mixing will not introduce static speckle as long as the entire etendue of the mixer is uniformly filled by the light source ( assuming coherent radiation ) .
this might seem to contradict phenomenological theories that claim that there should aways be some pattern in the detector plane because light from the same point of the light source can reach the detector along different paths , thereby resulting in uncontrolled interference that inevitably creates intensity fluctuations .
however , the unitary condition @xmath168 ensures that , whenever a mode of the light source is brighter at the detector , another mode has to be dimmer .
this property is fulfilled for any individual @xmath175 since @xmath219 , and no averaging over @xmath175 is necessary to arrive at this conclusion .
technologically , mixing apparatuses have a well - defined etendue , hence a well - defined value of @xmath169 .
when the etendue that they can accept , is filled to a large extend by the incident radiation , the effect of the unitary condition is strong , and the exact formula  ( [ eqtemp2exakt ] ) has to be used . for random scattering , on the other hand ,
in particular by a diffuser , frequently @xmath220 , and therefore the approximate formula  ( [ eqvarnstatisch ] ) might suffice .
in sec .  [ secstaticspeckle ] we demonstrated that the static speckle depends on the mean squared intensity in each mode of the radiation source , hence on @xmath221 , cf .
( [ eqtemp2exakt ] ) and  ( [ eqvarnstatisch ] ) .
the effect is minimal if all modes carry equal weight , and the maximum number of modes possible follows from the etendue via sec .
[ secetendue ] . the minimum amount of static speckle possible
can then be computed  it is zero if the etendue is not increased by the scattering and mixing processes , and finite otherwise .
this raises the following question : is it possible not only to compute the minimum amount of static speckle for given etendue of the radiation field but also the `` typical '' amount of speckle ? we will demonstrate that static speckle contrast is a self - averaging quantity , i.e. , in the limit of many modes @xmath169
, the relative fluctuation while taking the average over the unitary group becomes small .
this implies that the ensemble average is characteristic for almost all realisation of the ensemble .    when @xmath21 are the modes of the electric field , the electric field @xmath222 is uniquely described by coefficients @xmath19 , @xmath223 and the static speckle contrast is determined by the effective number @xmath210 of modes in a more or less complicated way depending on whether the exact result or the approximation is used , @xmath224 ^ 2 } \;.
\label{eqneff1}\ ] ]    in sec .
[ secetendue ] the linear span corresponding to some given etendue @xmath11 was computed . in the absence of dynamic information
, one can not select the modes out of all possible sets of basis vectors , and any set of basis vectors could thus represent the modes .
we assume that the decomposition  ( [ eqfeld1 ] ) represents the actual modes @xmath99 of the radiation field . since , as just explained , we can not know @xmath99 , the equations for the static speckle contrast are instead evaluated using some other orthonormal basis @xmath160 .
the two bases are related by a unitary transformation , @xmath225 , and consequently @xmath226 , but @xmath227 is unknown for the reasons just explained .
if no additional information is available , the best available estimate is the average of @xmath228 ^ 2 }          \nonumber \\      & =   \frac{\sum_{n , i_1,\ldots , i_4 }          u^\dagger_{n , i_1 } u_{n , i_2 } u^\dagger_{n , i_3 } u_{n , i_4 }          a^*_{i_1}a_{i_2}a^*_{i_3}a_{i_4 }          } {          \bigl[\sum_n { \lverta_n\rvert}^2\bigr]^2 } \ ; ,      \label{eqneffmittel}\end{aligned}\ ] ] after integrating over the unitary group , i.e. , over all possible unitary matrices .
we study the average of @xmath229 instead of @xmath230 since it both is easier to compute and is the relevant quantity for the speckle contrast .
the variance @xmath231 ^ 2\rangle -\langle 1/n'_{\mathrm{eff}}\rangle^2 $ ] then specifies how well the average is representative for all possible members of the ensemble .
averaging over the unitary group , the nonvanishing terms are of the form @xmath232 where @xmath233 and @xmath234 are permutations of the numbers @xmath235 , and @xmath236 denotes the cycle structure of a permutation .
the coefficients @xmath237 depend only on the cycle structure of @xmath238 and have been tabulated  @xcite .    labelling the numerator in eq .
( [ eqneffmittel ] ) as @xmath239 , @xmath240 and counting all the permutations of the indices in it yields @xmath241 for the square , @xmath242 one finds @xmath243    the terms @xmath244 appearing in the averages above cancel against the denominator of eq .
( [ eqneffmittel ] ) .
this is as expected since , as explained , these coefficients depend on the base that was used to determine them , i.e. , they depend on the choice of one particular unitary matrix .
averaging over the unitary group , this choice must be irrelevant , thus the coefficients drop out , and the results depends only on the number @xmath169 of modes .
inserting eq .
( [ equmittel2 ] ) into eq .
( [ eqneffmittel ] ) and using the tabulated coefficients from ref .   gives the result @xmath245 which has an easy interpretation .
the term @xmath246 instead of the naive term @xmath169 in the denominator is due to the correlations imposed on @xmath227 as it is unitary , cf . the factor @xmath246 in eq .
( [ eqcexakt ] ) .
the @xmath247 in the numerator implies that every mode is occupied with an effective intensity of @xmath248 relative to the maximum intensity of any mode  a very plausible result since the intensity is a random quantity between @xmath249 and that maximum value .    a similar calculation using eq .
( [ equmittel4 ] ) gives for the square @xmath250 so that the variance becomes @xmath251    the variance decreases much faster with @xmath169 than @xmath229 or @xmath252 as the lowest - orders in @xmath180 cancel when computing the variance .
this implies that @xmath229 is a self - averaging quantity , and its value computed from the number @xmath169 and hence from the etendue of the radiation field gives a good estimation of the actual value of @xmath229 .
one should be aware , however , that this is a statistical statement : in `` almost all cases '' this statement is correct . by use of particular light sources it still is possible to have a very different value for @xmath229 .
this is also reflected in the calculation presented here .
equation  ( [ eqneff1 ] ) gives the actual value of @xmath253 which , of course , depends on the values of the @xmath19 .
the average value from eq .
( [ eqeffmittelwert ] ) no longer depends on the @xmath19 and can thus deviate .
the modes @xmath21 of the electromagnetic field are mutually orthogonal and stationary , cf .
( [ eqgamma2 ] ) , and all the time - dependence of the electromagnetic field is included in the frequency - dependence of the annihilation operators .
however , frequently it is rather inconvenient to describe time - dependence in this way .
an extreme example are pulsed lasers where it is more appropriate to describe the mode structure of each pulse separately .
this comes at a price , though , as the modes of one pulse are then not necessarily orthogonal to the modes of another pulse .
this is primarily an issue for the static speckle computed in sec .
[ secstaticspeckle ] as the mutual orthogonality of the modes was essential for the derivation presented there .
the static speckle contrast in eq .
( [ eqcexakt ] ) has an easy interpretation : except for a minor correction due to the unitary of the scattering matrix , every mode creates a mutually uncorrelated speckle pattern in the detector plane . the effect of a superposition of different modes is easily understood from this , and for the static part it is irrelevant if all the modes are excited simultaneously or rather one after the other . in this section
, we will show that the speckle patterns of two nonorthogonal modes ( i.e. , necessarily at two different times ) are correlated .
we will , without loss of generality , refer to these two times as first and second pulse , respectively .
the first pulse consists of modes @xmath18 and annihilation operators @xmath19 . for a fixed arrangement of mode structure of the light source , scattering between source and detector , and detector position ,
thus only averaging over the fluctuations of the radiation field , the mean photocount is , completely analogue to eq .
( [ eqafrequenz3 ] ) , @xmath254 we now consider a second set @xmath155 of modes of the light source but keep the scattering between source and detector and the detector position fixed .
the transformation from the @xmath18 to the @xmath155 is linear and can thus be described by a matrix @xmath166 , @xmath255 where the notation @xmath256 denotes the overlap integral .
the annihilation operators transform as @xmath257 and the mean photocount thus becomes @xmath258 the essential quantity is the cross - correlation between @xmath83 and @xmath259 after averaging over all possible scattering configurations , hence over all @xmath175 , and only the matrix @xmath166 may remain in the result .
a positive correlation means that the light patterns caused by the two radiation fields are similar .
when the photocounts of two such fields are subsequently integrated on the same detector , the static speckle contrast will be higher than naively expected by the inverse square - root of number - of - modes law .
the cross - correlation between @xmath83 and @xmath259 is given by @xmath260 and the average has to be taken over the matrix @xmath175 keeping the matrix @xmath166 fixed .
writing out the first term explicitly , @xmath261 the average yields nonzero terms only for @xmath262 , and the two cases @xmath263 and @xmath264 have to be distinguished .
the necessary averages have already been given in eq .
( [ eqmittelwertes ] ) , and the result is @xmath265 this can be simplified by noting that @xmath266 . the cross - correlator then becomes @xmath267 i_n i'_m \;.
\label{eqcrosscorrelator}\ ] ] this equation gives the expected result in the two extreme situations .
if all intensities are equal , the cross - correlation becomes zero since , as explained earlier in the context of eq .
( [ eqcexakt ] ) , the static speckle itself vanishes , and thus also the cross - correlation does .
if only a single mode @xmath27 is excited , the cross - correlation becomes proportional to @xmath268 .
the average value of @xmath269 when assuming a random @xmath166 is equal to @xmath180 , and the cross - correlator is again zero .
if @xmath18 and @xmath270 are identical , the cross - correlator becomes maximal .
the cross - correlator thus correctly quantifies if @xmath18 and @xmath155 are stronger correlated than in a random configuration , and the overlap integral @xmath271 is directly related to the static speckle pattern .    when a photodetector sums over @xmath84 pulses , the variance @xmath272 of the integrated photocount follows from the variances @xmath273 of the @xmath88-th pulse , given by eq .
( [ eqtemp2exakt ] ) , and from the cross - correlator @xmath274 , given by eq .
( [ eqcrosscorrelator ] ) , between the @xmath88-th and the @xmath275-th pulse , @xmath276    apart from the rather obvious relation  ( [ eqvarmehrere ] ) , the results from this section have another application . frequently , radiation fields are expanded not in their modes but rather in a set of `` convenient '' functions , such as gaussian beams or plane waves .
drawback of this approach is that this set then usually has to be overcomplete , hence it has more elements than there had been modes , and not all its elements are mutually orthogonal .
for example , two gaussian beams have an overlap that decreases exponentially with separation but is always larger than zero .
two plane waves restricted to a finite interval have an overlap that decreases exponentially with the angle between the two waves but is always finite .
depending on the way that this expansion is done , some coherence information might be lost , namely when doing the expansion such that only the time - averaged intensities at all field - points are matched , and each @xmath18 is then assigned an intensity @xmath28 .
this might sound like a bad approximation but effectively amounts to using the assumptions from sec .
[ secstatickorrelation ] . inserting the @xmath28 into eq .
( [ eqtemp2exakt ] ) would yield a low static speckle contrast since the number of @xmath28 s is large .
however , the correct equation to use is eq .
( [ eqvarmehrere ] ) where every term is evaluated for a single function @xmath18 , and the nonorthogonality of the @xmath18 implies that the cross - correlation terms are large .
this is as expected since introducing additional functions @xmath18 to describe the same system should not change its computed properties .
the purpose of this paper has been to give a consistent and concise description of multi - frequency and multi - mode effects on photodetection on a quantum level trying to avoid the use of uncontrolled approximations .
its focus lies on exactness , not on ease of presentation or application to some particular device .
it thus complements previous work , for example by joseph w. goodman or the group around christer rydberg  @xcite .    in sec .
[ seczeitunschaerfe ] we derived how the uncertainty of the photocount depends on the measurement time and the temporal coherence function . among others , it was shown that a description via the coherence time of the radiation is insufficient to arrive at an exact result .
similarly , in sec .  [ secstaticspeckle ] not only previously - known approximations of the static speckle contrast were retrieved , namely an inverse square - root law with the effective number of modes , but in addition also effects of finite etendue accepted by the `` mixing device '' or of noncoherent , e.g. thermal , radiation followed from the chosen mathematical formalism without additional effort .    some of the results of this paper can be explained , albeit only on a qualitative level , by the `` brick wall '' model : the radiation field is like a brick wall ( horizontally : time , vertically : position ) , and only bricks in the same column or row can cause interference effects . counting the fraction of such bricks
gives an intuitive explanation of the scaling of dynamic and static speckle contrast upon changing system parameters .
if this is sufficient , the methods and results presented in this paper are a bit of an overkill .
if , on the other hand , exact agreement between predication and performance of some device is essential , this paper offers the advantage of giving results using only well - defined input parameters , thus not relying on effective parameters that need to be tuned until the desired result is retrieved .
main function of any theory is to predict whether some device or experimental setup will work as expected and/or required . for state - of - the - art technological applications , the performance might improve by about @xmath277 from one generation of a device or setup to the next  but rarely more .
to predict if something actually is an improvement , the error due to approximations thus has to be much smaller this number . as an example , fig .
[ figkurvenformen ] gives an indication of the error made when using the coherence time instead of the exact shape of the frequency spectrum .
if improvement of dynamic speckle is a topic , one should thus use the exact formulae and refrain from approximations .    depending on the desired application , either the present paper or one of the previous works by other authors
will thus be more appropriate .
furthermore , we restrict us to universal problems .
for example , we study in sec .
[ secstaticspeckle ] only the variance of the photocount , hence the contrast of the dark and bright `` spots '' in the plane of the detector , but not their spatial correlations , i.e. , the size of these spots .
reason is that , while the contrast is universal , the size of the spots is not and depends on the setup .
the book of goodman  @xcite is to a large extend devoted to studying this question ( on the classical level ) for a large number of setups of technological or scientific importance .    in sec .
[ seczeitunschaerfe ] the uncertainty of the photocount for known radiation field has been computed , yielding in addition to the well - known shot noise term two terms depending on the frequency spectrum and on the photon correlations of the radiation source , respectively .
the importance of the frequency spectrum is obvious , in particular in view of the two examples presented in this paper . for the importance of the effects of noncoherent
radiation , the situation is less clear . for traditional thermal light sources
, the necessary information can be found in this paper ; for traditional lasers , the emitted radiation is coherent , and the question becomes trivial . for modern light sources , such as excimer lasers  @xcite that are increasingly used for high - power applications ,
insufficient data is available .
technological advances are fast , and no good characterisation of the photon correlations for modern devices seems to have been published . since this contribution to the noise inevitably becomes the dominant term once the measurement time is sufficiently long , this lack of knowledge is an actual issue .
the static speckle , computed in sec .
[ secstaticspeckle ] , depends on the number of excited modes of the radiation field ( and the energy distribution among them ) but , in addition , also on the number of modes that could in principle be excited .
this is closely related to the concept of etendue , cf .
[ secetendue ] .
even though etendue is a very basic quantity , its microscopic definition does not seem to have been addressed before .
the relation presented between number of basis function and macroscopic etendue has application beyond this paper as any modelling of a radiation field  for the purpose of simulation or for an analytic study  needs to start with a given number of modes , and knowing this number a priori makes this process much more efficient .    while the applicability of the relation between etendue and maximum number of modes is obvious , this is less the case for sec .
[ secstatickorrelation ] discussing the relation between etendue and average effective number of modes .
we demonstrated that this average is identical to the actual number for almost all radiation fields .
however , that is `` only '' a statistical statement : it assumes that all possible radiation fields with given etendue are equally likely , and it is still allowed that a few cases ( albeit of measure zero ) deviate strongly .
different light sources have different properties , and for any technological application , one will choose the most appropriate light source . selecting a single - mode laser that operates in a high - order hermite - gauss mode
would give the most extreme disagreement between the actual number of modes ( = one ) and the prediction from sec .
[ secstatickorrelation ] .
the results from this section can thus only be applied when the choice for a particular light source does not adversely impair the freedom of the mode structure of the generated light field .