The Quantum Regime Operation of Beam Splitters and Interference Filters

The presence of quantum Rayleigh scattering, or spontaneous emission, inside a dielectric medium such as a beam splitter or an interferometric filter prevents a single photon from propagating in a straight line. Modelling a beam splitter by means of a unitary transformation is physically meaningless because of the loss of photons. Additional missing elements from the conventional theory are the quantum Rayleigh-stimulated emission, which can form groups of photons of the same frequency, and the unavoidable parametric amplification of single photons in the original parame-tric crystal. An interference filter disturbs, through multiple internal reflections, the original stream of single photons, thereby confirming the existence of groups of photons being spread out to lengthen the coherence time. The approach of modelling individual, single measurements with probability amplitudes of a statistical ensemble leads to counterintuitive explanations of the experimental outcomes and should be replaced with pure states describing instantaneous measurements whose values are afterwards averaged.


Introduction
Recent developments in the integration of photonic devices for quantum information processing [1,2] are characterized by their capability to generate a single photon per radiation mode in order to give rise to two-photon destructive interference for temporally overlapping indistinguishable photons, which is commonly known as the Hong-Ou-Mandel (HOM) dip [3]. The reduction in the counting rate of coincident detections of photons at two spatially separated photodetectors is explained by opposite sign amplitudes for the probabilities of detecting each photon pair after having been reflected or transmitted in opposite directions by a beam splitter (BS). Yet, with only one pair of photons present in the experimental setup at any given time for an individual measurement, the two types of detection cannot take place simultaneously.
However, a distinction needs to be made between the mathematical formalism of quantum mechanics and the physical processes of quantum optics grounded in optical physics.
The conventional approach of modelling the beam splitter was outlined in ref. [4], concluding that, with identical and synchronized single photons, "the probability of detecting one photon at each output port is zero which is clearly non-classical". However, recent experimental results have clearly demonstrated the possibility of obtaining 100% visibility with no coincidence between two outputs with independent photons [5], thereby contradicting the conventional interpretation mentioned above.
As far as single photons are concerned, the quantum mechanics interpretation would allege that "one photon does not interfere with another one; only the two probability amplitudes of the same photon interfere with each other" [3]. Yet, a probability amplitude can be specified only at the end of a long sequence or an ensemble of measurements consisting of individual measurements. Each of those measurements delivers an instantaneous value created, e.g., through parametric amplification in the nonlinear crystal of the spontaneously down-converted photons [10][11][12] as outlined in the Appendix A. This group of photons will follow a straight-line of propagation by the recapture through stimulated Rayleigh emission of a photon absorbed by an electric dipole. The common experimental combination of a beam splitter and an interference filter, both made of a dielectric material, will also involve the physical process of photon coupling between beams carrying photons of the same frequency [8]. Synchronized groups of identical photons may collide inside a dielectric medium such as a BS or an interference filter (IF), with the possibility of exchanging or coupling of photons from one group to the other through quantum Rayleigh-stimulated emission. The emerging weaker group may not reach the destination photodetector because of quantum Rayleigh scattering. Other physical mechanisms of generating groups of identical photons can be identified as follows.
A quantum dot (QD) placed in a high finesse micro-cavity of a few wavelengths [13,14] and excited with a picosecond pulse, can emit a photon spontaneously and be re-excited within the duration of the same pulse. If the photon was reflected towards the QD, stimulated emission may occur due to the small dimensions of the micro-cavity. This will result in two or more photons leaving the emitter simultaneously, as well as a reduced lifetime of the excited state of the QD, manifesting itself as a higher decay rate overshadowing the Purcell effect.
"Resonators with small mode volume and high-quality factors (Q-factors) enhance the light-matter coupling" [14,15]. However, a high Q-factor indicates a build-up of energy inside the cavity, i.e., a large number of photons do not leave the cavity after being emitted and reaching the output facet for the first time. The early photons emitted by a QD embedded in a resonant cavity of a dielectric structure of distributed Bragg reflectors (DBRs) [13] may be reflected towards the excited quantum dot and amplified, thereby giving rise to a time-varying number distribution which emerges from the "single-photon source" in the form of discrete pulses of identical photons.
The time-varying optical distribution of the monochromatic spectrum of such groups of photons is represented by means of the mixed time-frequency analysis [16] and their physical, spatial distribution was derived in [8].
Additionally, an interference filter will randomly, internally reflect one photon, thereby destroying the needed synchronization. Similarly, the possibility of splitting an incoming group of photons by means of multiple internal reflections will result in a longer pulse or coherence time, which may interfere with the following pulse.

Shortcomings of the Particle-like Approach
The beam splitter is a critical component of any experimental setup for quantum regime measurements. Its role would be to split the input operators of photon creation by propagating them through a unitary transformation across the beam splitter to the output ports. The output quantum state is then generated by applying the creation operators to the zero-photon states, resulting in a mixed quantum state description of possible combinations that preserve the initial number of photons [4].
The physical measuring process of detection involves photons, which are physical entities, as opposed to probability amplitudes, which appear as a mathematical element at the completion of an ensemble of measurements. The wavefunction and operators are mathematical tools or modelling elements. The probability P PD of detecting one photon using a photodetector was derived as the expectation value of the number operatorâ †

PDâ PD
and for an arbitrary state |Ψ out , e.g., [6]: where A PD is the probability amplitude of photodetection through absorption. We recall that the creation operator acting to the left in Equation (1) also indicates an absorption process. The single modes at the input ports a and b of the beam splitter are populated by creation operatorsâ † andb † , respectively. Then, the beam splitter is said to transform the operatorsâ † andb † into the creation operatorsĉ † andd † at the output ports of the beam splitter [4] through a unitary operation: where the reflection and transmission coefficients, respectively, are r ij and t ij , with |r ac | 2 + |t ad | 2 = 1 or |r bc | 2 + |t bd | 2 = 1 as a condition of photon conservation. Yet, these coefficients are evaluated for optical classical waves and may not be relevant for the scattering of single photons.
With the input photonic state written as [4]: the output state becomes [4]: An often-used case would have one input photon randomly exiting from one of two outputs, which leads to n = 1 and m = 0 in Equation (4), that is: Separate detections through absorption would be described by: For the case of the two BS outputs being directed to the same photodetector PD (or joint detection), and as the photon and its operator go together, the annihilation operator of the photodetector acts on both incoming states of Equation (5), leading to the probability of detection for indistinguishable pathways: For |r| = |t| = 1/ √ 2, e.g., a symmetric 1 × 1 Mach-Zehnder interferometer (MZI), and for one input of a cubic BS, r ac = −t ad , the probability vanishes P PD = 0 and photons seem to disappear even though the Hanbury Brown and Twiss (HBT) measurement would indicate the presence of a photon in either arm at any time. Similarly, for the other input with t bc = r bd , the probability P PD = (r + t) 2 > 1, contradicting, on physical grounds, the rule of probability amplitudes that was devised to explain experimental results mathematically. Therefore, the following statement needs to be revaluated: "If the different possible photon paths from source to detector are indistinguishable, then we have to add the corresponding probability amplitudes before squaring to obtain the probability" [3].
As the HBT measurement indicates that the two states |1 c |0 d and |0 c |1 d are mutually exclusive, the physically measured output should be derived from a time-dependent state describing each individual measurement: The case of two input photons exiting randomly from the two outputs of a BS is described from Equation (4) with n = 1 and m = 1, by: after setting r ac = r bd = r, t bc = t ad = t, and r = it . In Equation (9c), the time dependence of one single measurement is denoted by p ij (t M ).
With only one pair of photons, allegedly, present in the experimental setup at any instant of a single measurement, only one component can be measured. The time-dependence of the four terms leads to only one coefficient in Equation (9c) being non-zero. The timeor ensemble-averaged values of the first two terms are p rr (t M ) = r 2 and p tt (t M ) = t 2 . The physical process that implements the mathematical expression of Equation (9) may consist of a collision inside the BS of two groups of identical photons. Coupling between the two groups will be induced by the quantum Rayleigh-stimulated emission [8], depending on the relative phase of the beams.
The first two terms of Equation (9) cannot, mathematically, cancel each other out, if only because the counting consists only of positive numbers. Consequently, only an interference of instantaneous fields could lead to a vanishing value [8]. Additionally, other physical processes are ignored. For example, a photon can be reflected from the upper interface of a glass plate BS (see Figure 1), while the other photon will be reflected from the lower interface. In this case both reflection coefficients as well as the transmission coefficients will be equal with the same sign, making the cancelling of the two amplitudes impossible.
where the time coordinate of one type of measurements is . The conservation of the photon number leads to the condition that ( ) = 1 and ( ) = 0 , or ( ) = 0 and ( ) = 1 for each individual measurement. In this case, two statistical sets are obtained with their ensemble or time averages being ( ) ̅̅̅̅̅̅̅ = and ( ) ̅̅̅̅̅̅̅̅ = . The case of two input photons exiting randomly from the two outputs of a BS is described from Equation (4) with n = 1 and m = 1, by: The physical process that implements the mathematical expression of Equation (9) may consist of a collision inside the BS of two groups of identical photons. Coupling between the two groups will be induced by the quantum Rayleigh-stimulated emission [8], depending on the relative phase of the beams.
The first two terms of Equation (9) cannot, mathematically, cancel each other out, if only because the counting consists only of positive numbers. Consequently, only an interference of instantaneous fields could lead to a vanishing value [8]. Additionally, other physical processes are ignored. For example, a photon can be reflected from the upper interface of a glass plate BS (see Figure 1), while the other photon will be reflected from the lower interface. In this case both reflection coefficients as well as the transmission coefficients will be equal with the same sign, making the cancelling of the two amplitudes impossible. The physical contradictions arise because of the error in the original assumptions: "We characterize interference by the dependence of the resulting light intensities on the optical path length or phase shift, but we need to make a distinction between the measurement of a single realization of the optical field and the average over an ensemble of realizations or over a long time. A single realization may exhibit interference, whereas an The physical contradictions arise because of the error in the original assumptions: "We characterize interference by the dependence of the resulting light intensities on the optical path length or phase shift, but we need to make a distinction between the measurement of a single realization of the optical field and the average over an ensemble of realizations or over a long time. A single realization may exhibit interference, whereas an ensemble average may not. We shall refer to the former as transient interference, because a single realization usually exists only for a short time" [3]. In other words, interference of probability amplitudes was adopted as a rather superficial explanation.
The physical deficiency of Equation (4) is the omission of the photon-dipole interactions as a photon propagates through a BS. A source of photons is required to activate the mathematical expressionâ † 0 a = 1 . The relevant Hamiltonian of inter-action between two modesâ 2 andâ † 1 mediated by an optical susceptibility χ (1) is [8]: , for the same optical frequency ω.

The Physical Operation of Beam Splitters
With, allegedly, only a single photon per radiation mode, a beam splitter is meant to divide the field operators by invoking the Fresnel reflection and transmission coefficients of classical optics, while an interference filter is supposed to shape the spectral distribution of the optical operators [6].
Yet, from a physical perspective, the identical configuration, devices, and methods of operation of two separate systems are the reasons that bring about measurement correlations between the systems [8], rather than entanglement between single photons, which are, in fact, likely to be scattered [10].
An indication of the existence of the quantum Rayleigh scattering can be seen from the extensive loss of photons that has been a constant feature of photon coincidence counting. For example, ref. [17] reports on page 3 of the Supplementary Information: "The success probability of the entanglement generation process, i.e., detection of a photon after an excitation pulse, equals 5.98 × 10 −3 and 1.44 × 10 −3 for Alice's device and Bob's device, respectively. "A typical percentage of lost photons is at least 99.9%, as mentioned independently. The experimental results reported in [13] indicate a very low level of detection probability, with an excitation rate of about 0.5 MHz as shown in its Figure  2c, whereas the maximal coincidence rate was about 175/s as shown in its Figure 3 [13]. Similarly, ref. [15] reports that "single photons were generated at a rate of 740 kHz, chosen to match the path length difference between the EOM fiber and the delay fiber. With this configuration, two independent, subsequently-generated photons can arrive at the beam splitter simultaneously". Yet, the maximum detected number of coincidences per second was only about 75, as displayed in its figures.
Even recent experiments using optically nonlinear crystals for the parametric downconversion of photons report detection probabilities lower than 1%. Post-selection of data for a particular purpose is a typical procedure as "The raw data are sifted" [17]. All these bring to the fore the unavoidable amplification of spontaneously emitted photons as outlined in the Appendix A, which leads to the formation of pulses composed of a number of identical frequency photons. Such a pulse is described, physically, by means of dynamic and coherent number states [8] and, mathematically, by the mixed time-frequency representation [16].
Physical processes taking place inside a glass plate ( Figure 1) or cubic ( Figure 2) beam splitter (BS) are identified as follows:

1.
A stream or sequence of single photons is launched into each input of a BS. In the case of a glass plate, some photons may be reflected at the external interfaces as explained in Section 3.

2.
Inside the BS, the single photons are scattered by quantum Rayleigh spontaneous emission involving photon-dipole interactions.

3.
As the number of photons scattered inside a cubic BS increase, more and more dipoles are excited, so that after an initial build-up, a threshold is reached when a single photon can propagate in a straight line without being absorbed, and may even be amplified though quantum Rayleigh-stimulated emission, creating a group of photons.

4.
At the internal dielectric interface of the cubic BS, a group of photons is split arbitrarily between reflection and transmission. For example, out of 10 photons, three are reflected and seven are transmitted, or vice versa. 5.
The larger groups will survive the propagation to the output facet by recapturing any absorbed photons through quantum Rayleigh-stimulated emission, while the smaller groups will be scattered.  6.
In this way, only one of the two output photodetectors will be triggered, resulting in no coincidence count. 7.
Photon reflections at various dielectric interfaces of a beam splitter will destroy the time synchronization needed for the two-photon HOM dip to appear.
cess probability of the entanglement generation process, i.e., detection of a photon after an excitation pulse, equals 5.98 × 10 −3 and 1.44 × 10 −3 for Alice's device and Bob's device, respectively. "A typical percentage of lost photons is at least 99.9%, as mentioned independently. The experimental results reported in [13] indicate a very low level of detection probability, with an excitation rate of about 0.5 MHz as shown in its Figure 2c, whereas the maximal coincidence rate was about 175/s as shown in its Figure 3 [13]. Similarly, ref. [15] reports that "single photons were generated at a rate of 740 kHz, chosen to match the path length difference between the EOM fiber and the delay fiber. With this configuration, two independent, subsequently-generated photons can arrive at the beam splitter simultaneously". Yet, the maximum detected number of coincidences per second was only about 75, as displayed in its figures. Even recent experiments using optically nonlinear crystals for the parametric downconversion of photons report detection probabilities lower than 1%. Post-selection of data for a particular purpose is a typical procedure as "The raw data are sifted" [17]. All these bring to the fore the unavoidable amplification of spontaneously emitted photons as The propagation phase acquired by a photon (or photons) arises from dipole-photon interactions [8] through stimulated emission, and despite the small probability of absorption and spontaneous emission of a single photon, given the large concentration of molecules forming the dielectric device (see Avogadro's number), the scattering of single photons is certain.

The Glass Plate Beam Splitters and the HOM Dip
For a plate beam splitter sketched in Figure 1, the primary reflected and transmitted optical fields will lead to the transformation: where the subscripts of the reflection r ij and transmission coefficients t ij indicate the boundary interface, with the first subscript corresponding to the incoming direction of the photons. For a symmetric glass plate BS [3], the two incoming photons may be reflected with the same reflection coefficient, one at each of the external interfaces, without entering the BS. In this case, the two photons cannot cancel each other out, which would modify the statistical outcome. If one photon enters the BS and is reflected from the inside at the other interface, then the photon may be subject to quantum Rayleigh scattering and the time-synchronization is lost and, again, the statistical distribution would be modified.
These two options contradict the conventional interpretation of the HOM effect [3], claiming that the terms r 2 |1 c |1 d + t 2 |1 c |1 d disappear together for identical and synchronized photons launched into the input ports. These two terms have no timedependence in the mixed quantum state even though they occur at different times.
Consequently, the following description [18] is questionable and contradictory: "The phenomenon [HOM dip] is a beautiful manifestation of the interference of a quantum field . . . ". "If the photons were delayed from each other so that they did not arrive together at the beam splitter, then they would not exit together, and thus would not result in the two detectors at the output ports both registering a signal simultaneously. The change in the rate at which the detectors register in coincidence as the delay is varied gives the signature HOM dip, which provides a test of the nonclassical nature of the input light".
If two identical incoming photons meet inside the BS and cancel each other out or coalesce into a group of two photons, there is no role for the reflection and transmission coefficients in creating the HOM effect [3], which is also interpreted by means of the overall output distribution of measurements or the mixed quantum state in terms of the probability amplitudes.

The Cubic Prism Beam Splitter
For a cubic beam splitter made up of two butting prisms of different refractive indices as sketched in Figure 2, the field transformation becomes: The subscript 0 indicates free space, and subscripts 1 and 2 refer to the refractive index of each of the two prisms. In the case of the cubic prism beam splitter, any two external surfaces can form a resonant cavity through reflection or transmission at the diagonal interface. The output states will be affected by photons temporarily trapped inside the beam splitter, leading to the possibility of additional quantum states, as well as quantum Rayleigh coupling of photons [8].
The resonant cavity formed by two adjacent facets and the diagonal interface traps photons through quantum Rayleigh spontaneous emission or scattering. As a result, as more single photons enter the cavity sequentially and bounce back and forth in a Brownian motion, the number of excited electric dipoles increases, enabling a later photon to become amplified through quantum Rayleigh-stimulated emission, so that a small group of photons could propagate in a straight line.
It is claimed in [15] that with an initial state |Ψ i = |1 A |1 B , "Due to the indistinguishability of the photons, the detection of one photon in output port C or D at time t o projects the initial product state |Ψ i into the "which path" superposition state Ψ ± (t o ) = (|1 A |0 B ± |0 A |1 B )/ √ 2 of the remaining photon". In physical terms, the absorption of the first photon leads to the presence of another two photons in the pathwayentangled state. This is physically impossible at the level of one single measurement. Physically, describing time-dependent single quantum events by means of a time-independent quantum mixed state leads to counter-intuitive interpretations, which are rebutted under close scrutiny [8]. It was suggested that " . . . introducing an arbitrary differential phase ∆ϕ between the two components of |Ψ ± (t o ) . . . results in a phase-dependent wave function of the remaining single photon | Ψ ± (t o + τ) = (|1 A |0 B ± e i∆ϕ (τ) |0 A |1 B )/ √ 2 " [15] so "that the probability to detect the two photons in opposite output ports depends on the magnitude of the applied phase shift between the detection times", i.e., P coinc = sin 2 (∆ϕ/2)" [15]. However, this is physically impossible for one photon per radiation mode because of the quantum Rayleigh scattering [7], but it is possible with dynamic and coherent number states [8] that enable a smooth transition from the quantum to the classical regimes.

The Fiber-Optic Beam Splitter
The normalized longitudinal f ph and lateral f lat optical field profiles of a group of monochromatic photons were derived in [8] and have the form of a Wigner spectral component S (ω, t), which is a locally varying spectral amplitude [16], as opposed to a time-constant amplitude and phase of a Fourier spectrum, crossing a surface perpendicular to the wavevector of propagation. They are [8]: with ε the permittivity of the medium, z being the longitudinal coordinate given by the product of the speed of light c in the medium and the propagation time t. The radial distance from the peak of the field is denoted by ∆r. When a large number of photons propagate along the waveguide, the longitudinal field is multiplied by the number of photons, while the transverse field will fill the waveguide mode.
For an optical directional coupler such as the fiber-optic beam splitter, the evolution of the number of photons and their phases are presented in ref. [8], with the possibility of one waveguide capturing most of the photons, resulting in an asymmetric output. The coupling coefficient κ will have to take into consideration the temporarily discrete nature of the groups of photons by including the longitudinal field profile f ph next to the transverse spatial field f = f lat , of polarization e 1 or 2 , that is [8]: with indicating the reduced Planck constant, and z/t = v p is the phase velocity. Additionally, k o and n specify the free-space wavevector and the effective refractive index, respectively. The susceptibility of the medium is χ (1) . This spatio-temporal overlap is necessary for the quantum regime of discrete groups of photons. The phase-dependent coupling of photons [8] is critical in the operation of the optical waveguide beam splitters by creating an asymmetric output with the adjustable phase difference [19][20][21].
The same considerations will apply to integrated waveguide directional couplers used as beam splitters.

The Interference Filters
Experimental configurations for, apparently, two-photon quantum beats, e.g., [19][20][21], employ interference filters in order to control the coherence length of the photons.
The conventional interpretation of quantum optic experimental outcomes is based on global and mixed quantum states of one photon per radiation mode [1,2]. The corresponding wavefunctions are associated with the statistical distribution of an ensemble of measurements conducted under identical conditions: "A single-photon source coupling only to a single loss channel would emit a state given by |Ψ = √ P 0 |0 + √ P 1 |Ψ (1) with P 0 + P 1 = 1, and P 1 is the probability that a single-photon emission has happened . . . " [14]. The wave packet |Ψ (1) = dωψ(ω)â † (ω)|0 would be a superposition of Fourier spectral components ψ(ω) measured at the output, but each single measurement of the ensemble identifies only one spectral component at any given time [13]. The field profiles would be evaluated from a Fourier transform F [ψ(ω)] even though the spectral components are time-dependent and not simultaneously present.
The expectation value of the photonic optical field operatorâ would be proportional to Ψ â Ψ = √ P 0 P 1 ψ(ω)dω , implying that for an ensemble of single measurements delivering P 0 = 0 and P 1 = 1, the photonic field would vanish for a detection process of photon-dipole interaction. This contradicts the experimental results. Similarly, each photon is monochromatic, and space-and time-limited, and cannot be described by a spectral Fourier component, which is constant in time and space. As derived in ref. [8], an intrinsic photonic field distribution is carried by each interacting photon without any dependence on the measured statistical distribution of the ensemble of the mixed state. A pure state delivers one single measurement [22,23], whereas a mixed state describes the statistical distribution of an ensemble of measurements [23]. A monochromatic pulse composed of a number of identical photons will be represented by a time-varying spectrum S (ω, t) [16] and the spatial, physical distributions are given by Equation (12).
In addition, it is claimed that the coherence length of single photons reaching a photodetector at different times, one photon at a time, can be shaped by changing the time-independent superposition of Fourier components. The Fourier transform providing the photonic pulse shape and related coherence length would be controllable with a spectral filter.
From a physical perspective, a Fourier transform-or a superposition of spectral components-necessitates the simultaneous presence of the entire range of spectral components. However, this is not the case when only one photon at any given time crosses an interference filter. A single monochromatic photon propagating through a Fabry-Perot type filter, or a Bragg refractive index grating in a fiber, will be delayed randomly by repeated internal reflections and will acquire an integer multiple of a bias phase or time-delay. The higher the internal reflectivity of the cavity, the longer some photons may bounce back and forth inside the cavity, resulting in a spread-out group of photons or a "longer photon" output, which is interpreted as a longer coherence length. Such optical signals are best described by means of the mixed time-frequency (or Wigner-type) spectrum, e.g., [16], with the frequency amplitude itself being a function of time S (ω, t), specifying, in other words, a time-varying number of monochromatic photons being carried by time-varying photonic beam fronts. The time-stretching of the photonic group will be equivalent to pulse expansion for a narrower Fourier spectrum.
Thus, a group of photons simultaneously entering a resonant cavity of an interference filter will exit at different times as the higher the internal reflectivity, the longer the time that some photons will bounce back and forth inside the cavity. This process will cause the initially bunched group of photons to spread out in time and give rise to a longer coherence length for photon coincidence counting as in [20]. The wavefunction describing this output would take the following form at the spatial point of interaction r: where the times t m specify the existence of a group of n photons at location r. Recalling the non-Hermicity of the field operator [8], we find thatâ n = √ n e −i ϕ n − 1 , and the pure dynamic and coherent number state |Ψ n (ω, t) = ( |n (t) +|n(t) − 1 )/ √ 2 of a photonic beam front is monochromatic and time-dependent [8]. By contrast, the overall mixed state of the ensemble is multi-chromatic and time-independent (e.g., the bi-photon wavefunction).
The temporal profile of the optical field carried by a photon or any instantaneous photonic wave front should be determined from a pure quantum state wavefunction because it should be unaffected by the spectral distribution of an ensemble of measurements. However, for interference to take place, at least two coefficients c n (r, t) have to be non-zero in Equation (14).

Physical Aspects of Quantum Photonic Devices
The concepts of quantum interference of probability amplitudes and the collapse of an entangled wavefunction upon measurement are commonly invoked to explain experiments that are, allegedly, based on single and entangled photons. Yet, "Some of the experiments have been performed at light levels in the quantum regime, however, and this suggests strongly that the devices will work in the same way given single-photon sources and detectors" [24] (p. 264).
The approach based on a single and entangled photon per radiation mode was predicated on some form of quantum correlations that may exist between two photons that shared a common interaction (e.g., spontaneous parametric down-conversion or beam splitter interaction), along with the assumption that probability amplitudes of detecting a photon can be manipulated through a lack of knowledge about its propagation pathway [3]. However, recently, quantum-strong correlation of intensities [5] and polarization states [25] were measured with independent photons and explained in [8] and [11]. Since the quantum Rayleigh scattering prevents a single photon from propagating in a straight line inside a dielectric medium such as a beam splitter or an interference spectral filter, only groups of identical photons will provide the necessary time-synchronization for correlated measurements of coincident photons.
Equally, the intrinsic optical field profile of a group of photons is independent of the type of source that emitted them. It is the transient photons that give rise to instantaneous interference rather than the average number of photons.
Let us compare briefly the quantum and classical phase-sensitive resolutions obtained by means of correlated detections. The enhanced resolution of the 2-photon N00N state can be found in [26] (p. 86 and Figure 3.1). For a phase ϕ the modulation function is: m q (ϕ) = 0.5 (1 − cos2ϕ) = sin 2 ϕ indicating a reduced period of oscillation.
In the classical case, the intensity of an interference term 2 √ N 0 N m cosϕ between a number N 0 of unmodulated photons and a number N m of phase modulated photons is split into two equal branches and each branch is detected separately. Then, the photocurrents are multiplied to obtain for the classical modulation function Again, one obtains the same halving of the period but with the advantage of controlling the signal amplitude through the input number of photons. By increasing the number of splitting branches and multiplying the photocurrents, a reduction in the overall modulation function will become apparent, with the benefit of a higher rate of phasedependence. Therefore, a high phase-resolution is available without the complications of the quantum formalism.

Conclusions
Quantum-strong correlations of photonic detections between two spatially separated photodetectors have been, previously, experimentally proven with independent photons in published articles. The mathematical quantum interference of probability amplitudes for possible pathways corresponds to the "classical" distribution of optical amplitudes for multi-photon states, if only because a single photon is deflected by quantum Rayleigh scattering. Three types of beam splitters have been presented in the context of multi-photon states, and the quantum Rayleigh-stimulated coupling between two overlapping groups of identical photons have led to physically meaningful explanations for the experimentally demonstrated zero-count of coincident detections between two photodetectors, commonly known as the HOM dip.
The existence of the quantum Rayleigh (QR) scattering was well documented back in the 1970s in textbooks [27,28] and its absence from the theory of Quantum Optics developed since the early 1980s is still a puzzling question. A possible answer would be that the "miracles" of quantum optics would have needed explaining by other physical means, requiring a multi-disciplinary approach.
The quest for single-photon sources and photodetectors is unnecessary because for the required correlations of separated outputs, only identical devices operated in identical ways are needed. Even for phase modulated interference, the sensitivities of the classical case equal the performance of the theoretical quantum case, but without the complexities of the latter.

Conflicts of Interest:
The author declares no conflict of interest.

Appendix A. Quantum Optic Second-Order Parametric Equations of Motion
For a second-order nonlinearity or a χ (2) -based interaction, the frequencies are related by the equality: ω 1 + ω 2 = ω 3 , the pump being identified with ω 3 . Quantum mechanically, the Hamiltonian of interactionĤ int , which describes the exchange of photons by stimulated emission, has the expression: where the reduced Planck constant and the relevant susceptibility are, respectively, denoted byh and χ (2) q , the latter including the constant of proportionality relating the Hamiltonian expressed in terms of the electric field operators to the HamiltonianĤ int associated with photon annihilation and creation operators,â j andâ † j , of the ω j field (j = 1, 2, 3), as defined in [8].
The equations of motion for the expectation values are found from the Ehrenfest theorem presented in [8]. The composite wave function of the three interacting waves or photonic beams is: |Φ(r, t) = | Ψ n, 1 (r, t) | Ψ n, 2 (r, t) | Ψ n, 3 (r, t) (A2) Analogously to the linear wave interactions of [8], one obtains the following coupled wave equation of motion: ∂ ∂t After converting to the number of photons, s j = N j exp(ϕ j ), the equation of motion (A3) provides the rates of change in the number of photons N j and the phase ϕ j as follows: where the gain coefficient g 1 includes an overall coupling coefficient κ 1 that depends on the polarization states of the photons through the tensorial structure of the classical nonlinear susceptibility χ (2) . The phase difference (or the relative phase) between the three waves is θ 321 , β is the propagation constant and z/t = v p is the phase velocity. In (A4e), k 0;1 and n specify the free-space wavevector and the effective refractive index, respectively, while f identifies the field distribution. The condition for optimal amplification is found from Equations (A4a)-(A4c), for ∆ β = 0, to be θ op = ϕ 3 − ϕ 2 − ϕ 1 = −π/2 and is reached through the phase-pulling effect [12]. The mixing of three photons of different frequencies can lead to quadrature-dependent coupling of photons, with one pump photon splitting into a signal photon of the same quadrature state (p), and the idler photon coupled into the second quadrature state (q). As a consequence of the relative phase being −π/2 in the phase-dependent coefficient of Equation (A4b), i.e., ϕ p − ϕ p,s − ϕ q. i = −π/2, these quadrature waves will be amplified, whereas those shifted by π/2 will be de-amplified.