Bell Correlations without Entanglement : A Local Wave Model Using Gaussian-Poisson Statistics and Single Count-Pair Selection

Based on the violation of Bell inequalities, it has been believed that the derivation of Bell correlations requires a quantum description that depends on entanglement. However, the present paper computes Bell correlations among polarization analyzer output intensities from two spatially separated sets of superposed random wave pairs. To obtain proper Bell correlations, the general statistics must be modified to represent single event pair selection. The correlations between analyzer output components are then in one-to-one correspondence with those computed from the entanglement formalism.


Introduction
The Bell theorem and experimental violation of Bell's inequalities are believed to make the description of Bell correlations impossible without the use of entanglement.If valid, this view would rule out the possibility of local descriptions of these correlations.However, it is shown in [1] and references therein that conclusions based on the Bell theorem as historically reasoned may not be logically drawn.Thus, it becomes appropriate to investigate possible alternatives to quantum descriptions that are widely believed to imply nonlocality.
Below, it is shown that properly normalized covariances resulting from a statistical model of local "twin beams" from a down converter (see Figure 1) exhibit a one-to-one correspondence with Bell correlations.The experimental selection of single photon pairs from the down-converter source, and imposition of quan- tum source boundary conditions, result in modification of a general wave stochastic model to produce results in agreement with the entanglement description.
In Section 2 below, the correlation of polarization analyzer intensity outputs is first derived from the superpositions of orthogonally polarized light wave output from a type II down-converter as used in Bell experiments (see Figure 1).The correlation has the same form as a Bell correlation but is derived from continuous wave variables rather than individual photon-pair events.In Section 3, a semi-classical assumption is introduced: wave intensity is defined to be a statistic representing temporal count density.The two equal source output intensities become Poisson parameters for corresponding count-pair processes.When single count pairs are selected (per detector time constant) from the multiple possible analyzer output pairs, the resulting modified correlations become equal to those given by the quantum entanglement description.However, the process does not imply nonlocality, since it results from pairs of equal input intensities and associated twin counts originating in one local spatial region, i.e., the correlations result from common initial conditions imposed on each of two subsequently separated twin beams.
The formalism uses statistical assumptions common in the description of chaotic optical waves [2].It does not specify a mechanism for the association of photon counts with such waves.Imposition of count-pair statistics as a source boundary condition at rotation angles equal to zero results in Bell correlations at non-zero rotation angles due to random wave interference.
The overall model may be thought of as a generalization of wave particle duality.In two-slit interference, discrete entities assumed to be emanating from a source are manifested as discrete counts at a detector but exhibit the intervening effects of wave interference.In the model below, pairs of counts at initial analyzer angles equal to zero are inserted into, and determine, the statistical parameters of a wave system that determines the count correlations at non-zero analyzer angles.

Random Wave Description Based on Type II Down-Converter Outputs
Based on the geometry of parametric down converter sources for Bell correlation experiments (see Figure 1) [3], two complex amplitudes 1 U  and 2 U  are introduced corresponding to superpositions of orthogonally polarized light amplitudes in regions denoted as Beams 1 and 2 of the figure. 1 U  and 2 U  are superpositions of random waves: In these equations, horizontal and vertical polarization components are indicated by subscripts H and V respectively, with î and ĵ indicating unit vectors along the x and y axes.The regions of H and V component overlap are used as light sources in Bell experiments.Transient interference of the H and V components is assumed and is indicated by the amplitude superpositions in Equations (2.1) and (2.2).Following common usage in optics [4], such complex amplitudes may be stated in the more explicit form exp u u u iθ = if necessary, and thus contain both amplitude u , and phase information u θ characterizing a light wave.
The action of polarization analyzers placed in each of Beams 1 and 2 is now described.These transmit and reflect components of ( ) where 1H θ and 1V θ are the phases of the fields 1H u and 1V u , respectively, and the original source intensities at 1 0 θ = are 1 Similarly, for analyzer 2 one obtains ( ) To proceed, further conditions on the parameters in Equations (2.6)-(2.9)must be specified.First, imposing source characteristics used in Bell experiments [3] requires that 1 2 1 2 ; This is implied by the fact that photons are produced in pairs of opposite polarization in Beams 1 and 2. By applying Equations (2.10) to Equations (2.6)-(2.9) it follows that , and .
One concludes that for lossless analyzers, the sum of powers in the output polarizations equals the sum of powers in the input polarizations, which by Equation (2.10) are equal for the two beams.
A phase relation between the input beams must be imposed to make the analysis consistent with requirements of the nonlinear optics process [5] used to produce spontaneous parametric down-conversion (SPDC): where the ∆ 's are additional phase shifts implemented by a wave plate used in experiments [3].The difference of phases in the two beams is then θ θ θ π When quantities specified in Equations (2.10) and (2.14) are inserted into Equations (2.6)-(2.9),one obtains ) ( ) ) To compute the correlation of 1 S and 2 S , statistics of the amplitudes and phases are introduced that are consistent with experimental findings [6].These are: 1) The complex amplitude variables introduced in Equations (2.1) and (2.2) are circular complex Gaussian random variables [7].They are assumed to be statistically independent but subject to the phase condition introduced in Equation (2.14).
The properties that will be used below follow from assumption 1.
2) From 1, it follows that the wave intensities 1 are exponentially distributed.They are assumed to be statistically independent, but they have the same mean 0 I .
3) The phases are statistically independent of the amplitudes and are uniformly distributed from π to π − so that this condition also holds for θ in Equations (2.16a) and (2.16b).
Conditions 1 -3 are discussed in [2] and [7], and are commonly used to describe properties of light from thermal sources (and laser speckle patterns).Similar randomness is suggested by photographs of spatial distributions of down converted photons given in [8].
From points 2 and 3 above, two results are obtained that are used repeatedly in the following.First, if I is exponentially distributed, (2.17a) Second, a uniform distribution of θ over 2π implies averages (denoted by { } An important result of this section is now obtained by computing the correlation of 1 S and 2 S .First:

.20)
From the definition of the correlation of any two variables, the correlation of 1 S and 2 S , using Equation (2.20) is The square root variances of 1 S and 2 S in the denominator may be computed from Equation (2.18) by setting 1 2 θ θ = and taking the absolute value.Thus, intensity components, Equations (2.15a)-(2.15d),yield a correlation that is the same as the Bell correlation.However, the above calculation has been performed with random continuous variables, while the Bell correlation is based on measurements of discrete counts.

Modification of the Formalism to Represent Selection of Event Pairs
In spite of Equation (2.21), correlations among individual pairs of 1n I , 1 p I , 2n I , and 2 p I are not yet in a form that allows interpretation in terms of discrete quantum counts.However, the use of normalized covariances rather than the intensites of Equations (2.15) does result in such a form.It will be shown below that the covariances result from modification of the statistics to represent selection of single event pairs in a detector time-constant (as implemented in experiments), rather than arbitrary numbers of event pairs as implied by the statistics thus far.The appropriate modifications are obtained by subtracting their common mean value 0 I from each of the in- tensity components of 1 S and 2 S : The right-hand side of Equation (3.1) equals that of Equation (2.21), and it appears that nothing has been changed.However, Equation (3.1) may be written as a sum of terms resulting from the contributions of products  3) Experimental parameters and observational protocols are adjusted so that single count pairs are selected as assumed in the entanglement formalism.
To understand the reason for the quantization analysis to be given in Section 3.

Selecting Single Photon Pairs from Indefinite Numbers of Photon Pairs
First, consider Equation (3.4c).To measure correlations at the source, analyzer angles are set to 1 2 0 θ θ = = .Writing wave intensities as count averages yields The 0's and 1's in the last two terms are subscripted to indicate to which variable they belong.They represent selected count pair products for which each product equals zero.Note that the values in the right-most double sum occur with low probability, and are deleted under the protocols used in experiments.After averaging over the exponentially distributed variables, the average over the double sum is Thus, the constant on the right side of Equation (3.4c) corresponds to events that are excluded if single event pairs are selected from the range of possibilities.Clearly, the situation is the same for ( ) ( )   The probability product for duplicate photons equals the probability for one photon squared, i.e.

Discussion
It has been shown above that the interference of twin pairs of waves yields polarization analyzer intensity outputs whose appropriate cross correlations may be used to compute Bell correlations.The wave statistical parameters used are set equal to the basic observational count parameters found in SPDC.These results depend on interpreting wave intensity variables as counts per unit time and on the experimental restriction of correlations to represent single event pairs.When the general wave optical statistics are adjusted to represent single source-pair selection at analyzer angles equal to zero, the formalism provides the same results as the entanglement formalism at non-zero angles.This is a principle result of this paper.

Conclusions
It has been believed that the violation of the Bell inequalities from Bell's theorem implies that only the formalism of quantum entanglement can describe the observed correlations of Bell experiments.Although it produces the experimentally observed correlations, the entanglement formalism has the troubling implication of non-locality, or what some call "non-reality".The purpose of the present investigation has been to examine an alternative formulation of the problem that does not have these troubling properties.It is based on the superposition of duplicate pairs of orthogonally polarized waves that, due to quantum source boundary conditions, produce "twin" count pairs and hence "twin" wave intensity pairs.Non-locality does not appear to be a consequence of the formalism.Results of measurements on the spatially separated beams are correlated due to the duplicate wave intensities and counts.
No hidden variables have been used in the calculation.The wave optics formalism and assumed statistics for describing correlations among twin beam polarization components depend on empirical results.These lead to the following assumptions: 1) intensities are defined as counts/time; 2) counts are not divided at polarization analyzers; and 3) photons are created in duplicate at the source, with the two members of each source pair having opposite polarization.Condition 3) provides initial conditions for the physical processes that follow in each down-converter beam.The fact that photons occur in duplicate does not in itself imply nonlocality, since initially, in the down-converter crystal where the photons are created, the output waves are superimposed.
The formalism is silent on the question of the connection of photons to waves other than through the usual intensity definition of count rate.Thus, one obtains a statistical description in both the semi-classical wave optics picture presented, and the quantum entanglement picture.
The possibility of an alternative to the entanglement description has been thought to be ruled out by violation of Bell's inequalities.However, the Bell theorem is fundamentally flawed, as shown in [1] and references therein, so it may no longer be used ab initio to deny efforts to construct alternative descriptions of these phenomena.Some of the analyses given in the present paper were described in a preliminary (unpublished) form in [9].

Figure 1 .
Figure 1.Output of a spontaneous parametric down-converter (SPDC) in type II configuration.Light from the overlap regions of two light-cones constitutes Beams 1 and 2 as used in Bell experiments.Photons are produced in pairs on opposite sides of the arrow with one photon having î , and the other ĵ polarization.In the overlap regions of Beams 1 and 2, each polarization occurs with equal probability, with the two members of a pair having opposite polarization.Correlations of output counts from polarization analyzers (not shown) placed down-beam in each of Beams 1 and 2 are used to measure Bell correlations.
Equations (3.2a) and(3.2b)are equal to the corresponding quantum results, it is important to realize that at this point, little discreteness has been introduced into the model.To achieve a quantum interpretation of Equations (3.1)-(3.2), the continuous variables used are defined in terms of equivalent photon count statistics[2]:1) Wave intensities are defined as count densities in time.2) 1H I and 1V I are the means (Poisson parameters) of two independent processes that each produce duplicate count pairs as seen by analyzers at manner, the other averages may be computed.The relevant averages are 3.4a)-(3.4c)suggest Equations (3.2a) and (3.2b) (and quantum mechanical results), they have an additional constant offset, and an amplitude 2 0 2I multiplying the squared sinusoidal terms.It will now be shown that interpretation of intensities in terms of discrete counts, under the condition of single count pair selection, transforms Equations (3.4a) and (3.4b) into Equations (3.2a) and (3.2b).
(3.4c) holds for all equal analyzer angles so that the boundary condition at at other equal angles as well.Thus, Equation (3.4c) represents the result that follows from the selection of single photon pairs if the constant term is subtracted from both sides: this condition, Equation (3.4a) reduces to Equation (3.4c): is that Equations (3.4a)-(3.4c)represent observations under single count pair selection, if 2 0 I is subtracted from both sides of each of them.Similar considerations apply to the interpretation of the multiplicative factor Consider the average of the left-most expression in Equation (3.4a),

I
Poisson probability densities for photon counts on each side must be the same.Given that the probability of two counts on a side is made small compared to one count,2  Similarly, since the two possible polarization pairs correspond to two statistically independent processes with equal parameters, the second component of Equation (3.4a) yields must be subtracted in Equations (3.9b) and (3.9c) for correction of Equations (3.4a)-(3.4c), the required single-pair correction of Equations (3.9b) and (3.9c) implies ( corrected for single event-pair selection as used in experiments, Equation (3.4a) is modified to Equation (3.2a).From the above analysis, the same procedure modifies Equation (3.4b) to Equation (3.2b).
Two identical beams, "twin beams", have been created with one set rotated 90˚ with respect to the other.The output wave intensities of the analyzers are interchanged due to Equations (2.10) and (2.14).