# Unified quantum density matrix description of coherence and polarization

###### Abstract

The properties of coherence and polarization of light has been the subject of intense investigations and form the basis of many technological applications. These concepts which historically have been treated independently can now be formulated under a single classical theory. Here, we derive a quantum counterpart for this theory, with basis on a density matrix formulation, which describes jointly the coherence and polarization properties of an ensemble of photons. The method is used to show how the degree of polarization of a specific class of mixed states changes on propagation in free space, and how an interacting environment can suppress the coherence and polarization degrees of a general state. This last application can be particularly useful in the analysis of decoherence effects in optical quantum information implementations.

## I Introduction

Coherence and polarization are undoubtedly two of the most important properties of light. In general terms, the coherence of an optical field can be understood as the ability to produce interference, as remarkably demonstrated by Young in his famous double-slit experiment, and theoretically developed by the works of Fresnel in the context of waves hecht (). Another important development in the coherence theory was the one made by Glauber and Sudarshan, which established the connections about the coherence properties of light with the concept of photon statistics in a quantum mechanical scenario glauber (); glauber2 (); sud (). Conversely, the modern study of the polarization properties was introduced by Stokes, who proposed a set of parameters to completely describe the polarization state of a random electromagnetic wave; the so-called Stokes parameters mc (); mc2 (), that can also be extended to the quantum realm agrawal (). Together, these two concepts form the basis of numerous applications of light in microscopy moerner (), criptography bennet (); ekert (), metrology katori (), astronomy hanbury (); abbott (), as well as in future quantum information technologies nielsen (); kok ().

Although the importance of the theories of coherence and polarization, their theoretical descriptions have historically been developed independently born (); brosseau (); mandel (). However, since the last decade the study of these two apparently distinct properties could be established into a single formulation through the unified theory of coherence and polarization introduced by Wolf wolf (). In this seminal work, it was shown that both coherence and polarization of a random electromagnetic beam could be understood as manifestations of the correlations between fluctuations of the optical field. In this respect, coherence manifests itself from correlations between fluctuations of the electric field of a light beam at two or more points in space, whereas polarization arises from the correlations of the optical field components at a single point in space wolf2 ().

Since the publication of the unified theory, many other advances have been made towards a complete understanding of this problem. For example, the introduction of the generalized Stokes parameters koro (), the description of the polarization change of partially coherent electromagnetic beam upon propagation in free space koro2 (); salem (), and in the turbulent atmosphere koro3 (); roy (), just to mention a few. Nevertheless, almost all these works have been limited to the scope of the classical electromagnetic theory lind (); lepp (). In fact, there have been some recent works extending the classical unification theory to the realm of quantum mechanics by direct quantization of the electromagnetic field lahiri (); lahiri2 (). So far, this extension did not provide a significant clarification of the problem, when compared to the classical counterpart, maybe because the state of the field is characterized in the Fock space, which sometimes makes the physical intuition less precise and, depending on the environment in which the system is inserted, it is difficult to write an appropriate Hamiltonian to account for the time evolution of the system agrawal2 ().

In this work, we derive a unified quantum mechanical description of coherence and polarization from first principles, that is to say, without direct reference to the classical theory. As we shall see, the central element in this formalism is the density matrix of the system written directly in terms of the position and polarization Hilbert spaces. This last point is the responsible for making the method relatively simple when describing the behavior of a general ensemble of photons on propagation in free space, as well as under the action of an interacting environment. Indeed, we provide some applications of the model to demonstrate how a partially coherent ensemble of photons change the degree of polarization when propagating in free space, and how decoherence and depolarization take place when photons are subjected to an environment whose constituents can be refractive and birefringent. Since all these examples are presented by means of simple quantum-mechanical arguments, the present description can be particularly valuable in the study of environmental disturbance in optical quantum information processes, in which the properties of coherence and polarization play a fundamental role.

## Ii Theory

To start with, we derive an expression for the degree of spatial coherence of light in a context similar to the one used to derive the classical theory wolf (). In doing so, let us consider a Young’s double-slit experiment which consists in an ensemble of photons propagating close to the z-axis which are mostly blocked by a mask with two small openings on it. After this stage, the positions of the photons that passed through the slits are permanently registered by a distant detection screen, as shown in Fig. 1. Let and denote the quantum states of the photons which passed through the slits and , respectively, and and the states of the photons linearly polarized along the horizontal (x-axis) and vertical (y-axis) directions, respectively. In this scenario, we can write the general quantum state of the photons in the form

(1) |

with , in order to specify simultaneously both the slit in which the photon passes through and the state of polarization. Also, we can write the density matrix for this system as , which provides a matrix in the following format:

(2) |

where the asterisk denotes complex conjugation.

Now, if we are interested in computing the probability density to find a photon at a point on the detection screen, keeping in mind that horizontally polarized photons do not interfere with vertically polarized ones, we have that

(3) |

where and represent the states of photons localized at with horizontal and vertical polarizations, respectively. Assuming that the size of the slits is much smaller than the wavelength of the photons, we can consider that after passing through a given slit the wavefunction of the photons are spherical waves. Therefore, the probability amplitudes of finding a photon at with horizontal (vertical) polarization which passed through the slit () are given, respectively, by

(4) |

and

(5) |

with and being the imaginary unity and the wavenumber, respectively. The parameters and are the distances from the slits and to the point , respectively. By substitution of Eqs. (4) and (5) into Eq. (3), and cancelling out terms with inner products between horizontal and vertical polarization states, we obtain that

(6) | |||||

where we used the fact that is Hermitian, , and denotes the real part.

Let us visualize Eq. (6) under a different perspective. Observe that if the slit is closed, the amplitudes and are null in Eq. (1), therefore, Eq. (6) reduces to

(7) |

which represents the probability density of finding a photon that emerged exclusively from at . Similarly, the probability density of finding a photon that emerged from at is given by

(8) |

In this context, Eq. (6) can be rewritten as

(9) |

where the parameter is given by

(10) |

The first two terms in Eq. (9) correspond to the sum of the individual probability densities of the photons which passed through each slit, and the last term is responsible for the interference pattern on the detection screen. Note that the parameter that dictates the prominence of the interference pattern in this system is , which we define as the degree of coherence. Therefore, this parameter can be measured by detecting the patterns due to the photons which emerge from each slit separately, and the pattern formed when both slits are open, by means of Eq. (9).

By substitution of the amplitudes , , and in Eq. (10), and using the Cauchy-Schwarz inequality, it is easy to show that . Since the interference term is maximum when , we say that the ensemble of photons is completely coherent with respect to the slits and . On the other hand, for , the interference pattern is completely destroyed and we say that the ensemble of photons is completely incoherent with respect to the slits. In the intermediate cases, , we say that the photons are partially coherent. Despite the fact that we used a purely quantum mechanical method to derive the expression for the degree of coherence, Eq. (10), it has an interesting mathematical similarity with the one found in the classical theory (See Eq. (8) in Ref. wolf ()).

Now we turn to the analysis of the polarization in this system, and to the derivation of an expression for the degree of polarization. Initially, let us concentrate only on the photons emerging from the slit . At this point, we can define the quantum version of the Stokes parameters as given by (see a similar analysis in Ref. alte ())

(11) |

(12) |

(13) |

(14) |

where denotes the trace. Observe that, contrary to the second quantization formalism, in which the Stokes parameters are operators agrawal (); luis (); luis2 (), here they appear simply as numbers playing a role similar to that of the classical electromagnetic theory. Physically, the Stokes parameters above represent the ensemble average of the identity operator and the three Pauli operators in the basis brosseau1 (), but only for the photons which emerge from slit . In this form, they can be measured with an appropriate combination of linear polarizers and a phase shifter (See Refs. schaef () and alte () for a classical and quantum approach to this problem, respectively).

Accordingly, the degree of polarization at can be defined, and measured, in agreement with the following relation brosseau2 (); bjork ():

(15) |

which is a real quantity and its range is wolf2 (). After substitution of Eqs. (11), (12), (13) and (14) into Eq. (15), and some algebra, it can be written as

(16) |

which is our final expression for the degree of polarization. With this definition, we have that: (i) if , the ensemble is said to be partially polarized; (ii) if , the ensemble is unpolarized, and (iii) if , the ensemble is totally polarized. In a similar fashion, one can show that the degree of polarization of the photons which pass through the slit is given by

(17) |

Again, we call attention to the mathematical similarity between the expression of the degree of polarization derived here, with basis only on a quantum mechanical background, and the one found by means of classical methods wolf ().

So far, we have discussed the properties of coherence and polarization based on the pure state of Eq. (1). However, for a general mixed state

(18) |

with as the fractional populations of each pure state contained in the ensemble of photons, we can write the density matrix for the system as

(19) |

with each element given by . In this form, all the derivations developed above are equally valid for a mixed density matrix, including the expressions for the degree of coherence and polarization of Eqs. (10), (16) and (17).

Now we present some applications of the equations derived above in order to elucidate their physical meaning. Consider, for example, a pure ensemble of horizontally polarized photons with equivalent probabilities of passing through both slits. This is described by the state , such that the degrees of coherence as well as polarization are found to reach their maximum value, namely, , and . These are expected results since we have no information about which slit the photons passed, and all of them have a well defined polarization. In fact, for a general pure state described by Eq. (1), the two aspects that determine the degree of coherence are: (i) the relation between the probabilities of a photon to emerge from and , and (ii) the similarity between horizontally and vertically polarized photons with respect to the phase relation of the passage through the two slits. In the first case, the more distributed the probabilities of the photon to pass through each slit, the larger the degree of coherence. In the second case, the closer the relative phases between the passage through and for photons with polarization horizontal and vertical, the larger the degree of coherence. In regards to the degrees of polarization, it is easy to see that they reach their maximum for a pure state, . This is also expected since, independent of the amplitudes , , and , we always have a well defined polarization in the H-V basis for the state in Eq. (1).

Note that the present discussion is also valid for states whose coherent properties are entangled with the polarization properties. For instance, the state provides , and . The reason why the degree of coherence is null is because the polarization gives information about the path of the photons, eliminating the interference pattern. The degree of polarization is maximum because the polarization is completely defined at each opening, namely, horizontal at and vertical at .

The situation is much richer for mixed states. Here, we want to analyze two particular examples. First, consider the mixed state

which represents a completely random ensemble of photons with equal probability of being horizontally and vertically polarized, and passing through the slits and . Then, for obvious reasons, the degrees of coherence and polarization are null: , and . The next example is particularly interesting. It represents a mixed state whose coherence and polarization degrees of freedom are separable:

The degrees of coherence and polarization for this state satisfy , and . That is, the ensemble of photons is completely unpolarized at each slit; however, it is completely coherent with respect to them. Wolf called attention to this unusual behavior when the classical theory was developed wolf ().

## Iii Polarization change on propagation in free space

In this section, we use the theory that we developed to show how the degree of polarization of a mixed ensemble of photons is changed upon propagation along the direction. Let us consider the ensemble initially prepared (at ) in the following mixed state:

with and , which is composed of two groups of pure states that we assume to be, in the beginning, equally probable to be detected at the points and . Namely, note (). The first group is horizontally polarized, with equal probability of detecting a photon at and . The second is vertically polarized, but also with equal probability of detecting a photon at the two referred points. Moreover, let us assume that the spatial wavefunctions of each group perpendicular to the direction of propagation, , are Gaussians with different widths, which spread on propagation due to diffraction (See Fig. 2):

(23) |

and

(24) |

where and are the respective widths of the wavefunctions. Note that, for a given point along the direction the relative probabilities to find a photon described by (with horizontal polarization) or (with vertical polarization) change on propagation. This, as we shall see, is the physical reason for which the degree of polarization varies with .

The -dependent density matrix of this system in the basis is given by:

(25) |

Thus, the degrees of polarization at the points and are given by:

(26) |

At this point we call attention to the fact that the fractional populations and to find each subensemble at and , at , introduced in Eq. (III) are numbers which are determined at the moment of the creation of the ensemble. Then, it should be clear that the dependence on shown in Eq. (25) reflects that the probabilities to obtain a photon from each group changes upon propagation because of the independent evolution of and .

In order to estimate the probabilities to obtain each group of photons, we will assume that the subensembles comprise two Gaussian beams whose intensity profiles (probability density) are given by mandel (); siegman ():

(27) |

and

(28) |

where and are the minimum widths of and , respectively, which we assume to be at . We also have that , with . The parameters are the Rayleigh lengths of the beams described by and .

For the sake of simplicity, let us assume that the distance between the points and and the -axis are much smaller than the minimum width of the beams, such that we can write

(29) |

which after normalization we find that the matrix elements in Eq. (25) are given by

(30) |

Observe that , which are the fractional populations that we defined at the ensemble creation. Finally, if we substitute Eq. (30) into Eq. (26), we can find how the degree of polarization changes as the photons propagate along the -axis. Fig. 3 shows the degree of polarization of Eq. (26), with as given in Eq. (30), as a function of the propagation direction , for . The behavior is similar to the one found by using the classical theory applied to the Gaussian Schell-model james (). In this regard, it is important to emphasize that despite the similarity between the results obtained from the classical formalism and the quantum treatment, only the latter was able to clarify the physical reason of the polarization change effect in free space. Namely, the different evolutions of two independent subensembles contained in a photonic mixed state. This result is being presented here for the first time.

At this stage, we can point out that the unified theory of polarization and coherence was fundamental in the present description in the sense that, to completely describe the dynamics of the degree of polarization of the photons upon propagation, information about the influence of the dynamics of the phase relation between the two subensembles, which is related the (spatial) coherence properties of the whole ensemble, has to be taken into account. In this regard, the density matrix of Eq. (25) is the mathematical entity capable of providing such complete knowledge, in a similar fashion to that expected from the cross-spectral density matrix in the unified classical theory wolf ().

## Iv Decoherence and depolarization due to environmental interactions

In this section we analyze the effects of the correlations naturally created between the ensemble of photons and the environment constituents that may take place upon propagation. This effect, so-called environment induced decoherence, has been recognized as the responsible for the emergence of the classical behavior of light and matter from the underlying quantum substrate zurek (); schloss (); schloss2 (); bert (). In this sense, as we shall see, the decay of the coherence and polarization degrees due to the interaction with the environment appears to be a natural and irreversible process. In mathematical terms, the characteristic trait of decoherence is the decay of the off-diagonal elements of the reduced density matrix of the system obtained from the partial trace of the system-environment density matrix with respect to the environmental states zurek2 (); zurek3 ().

In order to account for this effect on the photonic states described here, we will describe the decoherence process in the formalism of quantum channels kraus (); nielsen (). One of the advantages of the method, also called operator-sum formalism, is that the influence of the environment on the photons can be described without specific reference to the interaction Hamiltonian. Thus, we initiate this study by defining the system-environment density matrix under the assumption that they are initially uncorrelated. In this form,

(31) |

where the density matrix of the system contains information of both coherence and polarization, as in Eq. (19), and the density matrix of the environment can be written in the diagonal decomposition as , with as the fractional populations of the environmental states in the basis . If we consider that the photons and the environment form an isolated system, they necessarily evolve under an unitary operation . Then, the evolution of the reduced density matrix of the photons is given by

(32) |

where denotes the partial trace over the states of the environment. This equation can be written as

(33) |

where are the so-called Kraus operators, which are given by

(34) |

It is easy to show that the Kraus operators obey the relation , where is the identity matrix in the Hilbert space of the system .

Now, for us to proceed, it is necessary to describe the characteristics of the environment to be experienced by the photons. As a first example, let us assume that the photons propagate through a region in space composed of many small transparent particles, randomly distributed in space, with an index of refraction different from that of vacuum. This type of environment, which could simulate the lower atmosphere weich (), as well as impurities in optical fibers wai (), causes random, uncorrelated phase shifts in the photon quantum state at both points and . Each possible single interaction of this type can be described by a unitary transformation that modifies uniquely the state of the environment in the following form:

(35) |

(36) |

(37) |

(38) |

The parameter in the above equations represent the probability for an interaction between a photon and an environment constituent to occur during a given time interval, . The states , and represent the initial state of the environment, and the states of the environment after the interaction with one photon at and at , respectively. Note that, in the present case, the change of the states of the environment does not depend on the polarization of the photon, but only on the localization of the interaction. Physically, we can attribute the change of the environmental state upon interaction to the momentum imparted to the atoms of the environment due to either the scattering of the photon or absorption and reemission of it. Also, note that this kind of interaction is unable to cause a transition in the basis , i.e., a photon at cannot be sent to due to the environment, and vice-versa.

Before evaluating the Kraus operators of Eq. (34), we call attention to the fact that in the present example . Then, from Eqs. (33) and (34) we have that the time-dependent density matrix can be simplified to

(39) |

with

(40) |

Accordingly, by using Eqs. (35) to (38) we find that the three possible Kraus operators in the basis are given by

(41) | |||||

(42) |

(43) |

Then, by using these results in Eq. (39) we have that the evolution of the density matrix after a time is given by

(44) |

If this operation is applied times in succession, the terms in the matrix above become . Also, if we assume that the interaction probability in the time interval is of the form , with as the probability of an interaction between a photon and an environment constituent per unit time, then, after a time , we have that . Thus, for we obtain preskill (). Therefore, the time evolution of the density matrix can be written as

(45) |

Now, from Eq. (10), we obtain that the degree of coherence of the ensemble of photons under the interaction with the environment decays in the following form:

(46) |

On the other hand, it is easy to see that the degrees of polarization, and , given by Eqs. (16) and (17) remain constant. This was expected, since the environment only disturbs the relative phase of the state of the photons with respect to the points and .

As a last example, we analyze the decoherence effect on the photons due to an environment whose constituents, besides causing random phase shifts to the photonic states as in the previous case, are also birefringent, i.e., the phase shifts now depend on the polarization state of the photon. In this case, after interaction, the photons modify the environment state in a form that depends both on the localization, and , and the polarization state, and . Under these assumptions, each possible single interaction can be described by a unitary transformation that changes the initial state of the environment as

(47) |

(48) |

(49) |

(50) |

Again, is the interaction probability. In this form, from Eq. (40) we can find the five possible Kraus operators related to this case

(51) | |||||

(52) |

(53) |

(54) |

(55) |

Thus, by substitution of Eqs. (52) to (55) into the time dependent expression for the density matrix,

(56) |

we have that

(57) |

which, similar to the previous example, if we assume that the interaction probability is linear with time, , after many interactions the terms becomes approximately , with as the time elapsed by the interactions. In this case, the evolution of the density matrix is given by

(58) |

Therefore, given the temporal evolution of the density matrix, with the formalism introduced here we can evaluate the evolution of the degrees of coherence and polarization with time. In this case, one can easily see from equation Eq. (10) that the degree of coherence has an exponential decay, , similar to the previous case, Eq. (46). Nevertheless, from Eqs. (16) and (17), we can verify that the degrees of polarization also decay with time according to

(59) |

and

(60) |

In this case, in which the light-environment interaction is ruled by Eqs. (47) to (50), we have a polarization-dependent decoherence, i.e., contrary to case of the last section, now the polarization properties cause influence on the coherence of the photons. This is where the importance of the unified theory comes into play. It would also be interesting to analyze the case in which an ensemble of photons in a mixed state like that of Eq. (III) propagates in the medium described by Eqs. (47) to (50), instead of in free space. In such scenario, one has interplay between the coherence and polarization properties of light. Such cross influence could be depicted under the perspective of the present unified framework.

We want to call attention to the fact that we have provided a simplified model for the interactions by assuming ideal (orthonormal) environmental states to illustrate the validity of the present model. In fact, as indicated above, the simplicity of the method lies in accounting for the evolution of the reduced density matrix of the system without specifying the interaction Hamiltonian with the environment. This approach delivers a compact and practical description for the dynamics of open quantum systems schloss2 (). However, it is not difficult to imagine a situation in which the states obtained after the interactions in Eqs. (35) to (38), as well as in Eq. (47) to (50), are not orthonormal as we assumed. If we relax this condition, it can be shown that the decay rate in the elements of Eqs. (45) and (58) are not necessarily the same, which would provide more interesting time evolutions for the degrees of coherence and polarization.

The knowledge about the dynamical properties of a light beam which propagates through a disturbing environment is very important and finds application in many fields, such as optical communications, remote sensing and radar systems wang (). However, the usual classical description of this problem may be cumbersome, especially when the environment is turbulent salem2 (). We believe that the quantum density matrix approach developed in this section opens a new avenue for investigations of the coherence and polarization properties of light under the action of many types of environments, once we know the quantum state transformations which rule the interactions of the photons with the environment constituents. Furthermore, as well known from decoherence theory, depending on the type of environment interacting with the system, the master equation formalism can also be applied bru (); rivas (). Here, we also have this option since our formalism stand on the density matrix of the system as the fundamental element.

Finally, we want to point out that, contrary to the classical framework, our quantum description of the problem can depict the coherence and polarization properties of a subensemble of photons which composes, for example, one party of an entangled multipartite system, e.g., one of the constituents emitted from a EPR or a GHZ source pan (). In this case, our four-dimensional density matrix used to obtain information about coherence and polarization would be the reduced density matrix of the subensemble of photons, which provides all the measurement statistics. As a matter of fact, the quantum density matrix method accounts for entanglement in the multipartite system both in the position of the photons and the polarization degree of freedom.

## V Conclusion

In conclusion, we have proposed a new unified theory of coherence and polarization based solely on first-principles quantum mechanical arguments. The theory relies on a density matrix written in terms of position and polarization states of an ensemble of photons, from which we derived expressions for the degrees of coherence and polarization of the system. To confirm the validity and efficiency of the model, it was applied to show how the degree of polarization of a mixed ensemble of photons varies on propagation in free space; a problem that, to our knowledge, has been studied only with basis on the classical electromagnetic theory. Furthermore, we successfully used our method to describe the behavior of the coherence and polarization properties of a generic ensemble of photons subjected to interactions with an external environment. In this case, we showed two examples of interacting environments: one causing random phase shifts at two different points perpendicular to the propagation direction, and another causing polarization-dependent phase shifts. In this context, we used the operator-sum representation to unveil the temporal evolution of the system. However, it is important to emphasize that, depending on the type of interaction in which the photonic system is submitted, the master equation formalism can also be used to describe the dynamics. In future works, we intend to use the present study to investigate the action of other types of environment by using master equations.

###### Acknowledgements.

The author is grateful to A. S. L. Gomes and Cid B. de Araújo for their hospitality at Universidade Federal de Pernambuco, and to the financial support from the Brazilian funding agency CNPq, Grant Number 309292/2016-6.## References

- (1) E. Hecht, Optics, Addison-Wesley, 4th Ed. (San Francisco, 2002).
- (2) R. J. Glauber, Phys. Rev. 130 2529 (1963).
- (3) R. J. Glauber, Phys. Rev. 131 2766 (1963).
- (4) E. C. G. Sudarshan, Phys. Rev. Lett. 10 277 (1963).
- (5) W. H. McMaster, Am. J. Phys. 22 351 (1954).
- (6) W. H. McMaster, Rev. Mod. Phys. 8 33 (1961).
- (7) G. Agarwal, Quantum Optics, Cambridge University Press (Cambridge, 2013).
- (8) W. E. Moerner, Rev. Mod. Phys. 87 1183 (2015).
- (9) C. H. Bennett and G. Brassard, Proceedings of IEEE International Conference on Computers, Systems and Signal Processing, 175, 8 (1984).
- (10) A. Ekert, Phys. Rev. Lett. 67, 661 (1991).
- (11) H. Katori, Nature Photon. 5, 2013 (2011).
- (12) R. H. Brown and R. Q. Twiss, Nature 177, 27 (1956).
- (13) B. P. Abbott et al., Phys. Rev. Lett. 116, 061102 (2016).
- (14) M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information, 1st Ed. (Cambridge University Press, 2000).
- (15) P. Kok and B. W. Lovett, Introduction to Optical Quantum Information Processing, Cambridge University Press, (Cambridge, 2010).
- (16) M. Born and E. Wolf, Principles of Optics, Cambridge University Press, (Cambridge, 1999).
- (17) C. Brosseau, Prog. Opt. 54, 149 (2010).
- (18) L. Mandel and E. Wolf, Optical Coherence and Quantum Optics, Cambridge University Press, (Cambridge, 1995).
- (19) E. Wolf, Phys. Lett. A 312, 263 (2003).
- (20) E. Wolf, Introduction to the Theory of Coherence and Polarization of Light, Cambridge University Press (Cambridge, 2007).
- (21) O. Korotkova and E. Wolf, Opt. Lett. 30 198 (2005).
- (22) O. Korotkova and E. Wolf, Opt. Commun. 246 35 (2005).
- (23) M. Salem and E. Wolf, Opt. Lett. 30 198 (2005).
- (24) O. Korotkova, M. Salem, and E. Wolf, Opt. Commun. 233 255 (2004).
- (25) H. Roychowdhurya, S. A. Ponomarenkob and E. Wolf, J. Mod. Opt. 52 1611 (2005).
- (26) J. Lindberg, T. Setälä, M. Kaivola, A. T. Friberg, Opt. Commun. 283, 4452 (2010).
- (27) L.-P. Leppänen, K. Saastamoinen, A. T. Friberg, and T. Setälä, New. J. Phys. 16, 113059 (2014).
- (28) M. Lahiri and E. Wolf, Phys. Rev. A 82, 043805-1 (2010).
- (29) M. Lahiri and E. Wolf, Phys. Rev. A 82, 043837-1 (2010).
- (30) G. S. Agarwal and R. R. Puri, Phys. Rev. A 40, 5179 (1989).
- (31) J. B. Altepeter, E. R. Jeffrey, P. G. Kwiat, Adv. At. Mol. Opt. Phys. 52, 105 (2005).
- (32) A. Luis, Prog. Opt., 61 283 (2016).
- (33) A. Luis, Phys. Rev. A, 66 013806 (2002).
- (34) C. Brosseau and A. Dogariu, Prog. Opt. 49, 315 (2006).
- (35) B. Schaefer, E. Collett, R. Smyth, D. Barrett, and B. Fraher, Am. J. Phys. 75, 163 (2007).
- (36) C. Brosseau, Fundamentals of Polarized Light: A Statistical Optics Approach, Wiley (New York, 1998).
- (37) G. Björk, J. Söderholm, L. L. Sánchez-Soto, A. B. Klimov, I. Ghiu, P. Mariand, T. A. Marian, Opt. Commun. 283 4440 (2010).
- (38) These fractional populations concern the probabilities to find photons from each subensemble at the points and (at z=0), not the probabilities to find them in the whole volume which contains the ensemble.
- (39) A. E. Siegman, Lasers, University Science Books (Mill Valley, 1986).
- (40) D. F. V. James, J. Opt. Soc. Am A 11, 1641 (1994)
- (41) W. H. Zurek, Rev. Mod. Phys. 75 715–775 (2003).
- (42) M. Schlosshauer, Rev. Mod. Phys. 76 1267–1305 (2005).
- (43) M. Schlosshauer, Decoherence and the quantum to classical transition, 1st Ed. (Springer, 2007).
- (44) B. de Lima Bernardo, Braz. J. Phys. 44 202 (2014).
- (45) W. H. Zurek, Phys. Today 67 44 (2014).
- (46) W. H. Zurek, Nat. Phys. 5 181 (2009).
- (47) K. Kraus, Effects, and Operations: Fundamental Notions in Quantum Theory, Lect. Notes Phys., Vol. 190 (Springer-Verlag, Berlin, 1983).
- (48) H. Weichel, Laser Beam Propagation in the Atmosphere SPIE Optical Engeneering Press (Bellingham, 1990).
- (49) P. K. A. Wai and C. R. Menyuk, Opt. Lett 19 1517 (1994).
- (50) J. P. Preskill, Quantum Computation Lecture Notes for Physics 219/Computer Science 219 (http://www.theory.caltech.edu/people/preskill/ph229/).
- (51) F. Wang, X. Liu, and Y. Cai, Prog. Electromagn. Res., 150, 123, (2015).
- (52) M. Salem, O. Korotkova, A. Dogariu, and E. Wolf, Waves Random Media, 14 513 (2004).
- (53) H.-P. Breuer, F. Petruccione, The Theory of Open Quantum Systems, (Oxford University Press, Oxford, 2002).
- (54) A. Rivas, A. D. K. Plato, S. F. Huelga, and M. B. Plenio, New J. Phys. 12 113032 (2010).
- (55) J.-W. Pan, Z.-B. Chen, C.-Y. Lu, H. Weinfurter, A. Zeilinger, and M. Zukowski, Rev. Mod. Phys. 84 777 (2012).