Theory of Quantum Path Entanglement and Interference with Multiplane Diffraction of Classical Light Sources

Quantum history states were recently formulated by extending the consistent histories approach of Griffiths to the entangled superposition of evolution paths and were then experimented with Greenberger–Horne–Zeilinger states. Tensor product structure of history-dependent correlations was also recently exploited as a quantum computing resource in simple linear optical setups performing multiplane diffraction (MPD) of fermionic and bosonic particles with remarkable promises. This significantly motivates the definition of quantum histories of MPD as entanglement resources with the inherent capability of generating an exponentially increasing number of Feynman paths through diffraction planes in a scalable manner and experimental low complexity combining the utilization of coherent light sources and photon-counting detection. In this article, quantum temporal correlation and interference among MPD paths are denoted with quantum path entanglement (QPE) and interference (QPI), respectively, as novel quantum resources. Operator theory modeling of QPE and counterintuitive properties of QPI are presented by combining history-based formulations with Feynman’s path integral approach. Leggett–Garg inequality as temporal analog of Bell’s inequality is violated for MPD with all signaling constraints in the ambiguous form recently formulated by Emary. The proposed theory for MPD-based histories is highly promising for exploiting QPE and QPI as important resources for quantum computation and communications in future architectures.


Introduction
Quantum temporal correlations are analyzed with diverse methods by utilizing histories or trajectories of evolving quantum systems with more recent emphasis on mathematical formulation of the entangled superposition of quantum histories in Reference [1], i.e., denoted with the entangled histories framework. These varying methods include Feynman's path integral (FPI) formalism [2] as the most fundamental of all inherently including histories, consistent histories approach defined by Griffiths [3][4][5], and the recently formulated entangled histories framework [1] and two-state vector formalism [6,7] while all emphasizing correlations in time as standard quantum mechanical (QM) formalisms without violating Copenhagen interpretations. Multiplane diffraction (MPD) design as a simple linear optical system was recently proposed for quantum computing (QC) [8,9] and for modulator design in classical optical communications [10] by exploiting the tensor product structure of quantum temporal correlations as quantum resources while utilizing only the classical light sources and conventional photon-counting intensity detection. The MPD architecture generates interference of an exponentially increasing number of propagation trajectories along the diffraction events through multiple slits on the consecutive planes. The simplicity of source and detection in MPD setup combined with the highly important promise of the utilization of the tensor product structure of the temporal correlations as quantum resources motivates the definition and study of quantum trajectories or histories in MPD as novel quantum resources. These new resources denoted as quantum path entanglement (QPE) and quantum path interference (QPI) are defined and theoretically modeled in this article in terms of the temporal correlations and interference among the trajectories, respectively, to be exploited for future quantum computing and communications systems.
In this article, MPD design is, for the first time, proposed for defining QPE and QPI as novel quantum resources. Operator theory modeling for MPD-based resources is presented by combining the consistent histories approach of Griffiths [1,[3][4][5] and the entangled histories framework in Reference [1] with the FPI approach as the inherent structure of MPD creating Feynman paths. MPD creates quantum propagation paths through individual slits in a superposition in which the linear combinations result in evolving quantum history states. It has low experimental complexity with classical light sources and conventional photon-counting detection for near-future experimental verification. The theory of QPE and QPI based on MPD proposed in this article provides a set of tools to explore new structures composed of the correlations and interference among the paths for future applications in quantum computing and communications and provides QM foundational studies based on quantum histories.
The concept of the entangled histories is defined in References [1,11] as the quantum history which cannot be described as a definite sequence of states in time. There is a superposition of multiple timelines of sequences of events. In this article, we follow similar terminology and denote the temporal correlation among the quantum propagation paths unique to the MPD design with QPE, i.e., emphasizing the entanglement among the path histories similar to References [1,11]. Tensor product structure among the temporal correlations of multiple time instants is utilized as a novel resource for computing in References [8,9] and for communications in Reference [10] in an analogical manner to the multiparticle spatial correlations of the conventional quantum entanglement resources. MPD provides a simple system design inherently including such states having correlations among the paths denoted with QPE. A concrete example of a history state in MPD composed of diffraction events through N planes is defined as follows: ∑ n π n P N,s n,N P N−1,s n,N−1 ... P 1,s n,1 where P j,s n,j is the projection operator for diffraction through the slit indexed with s n,j on jth plane and for nth trajectory, π n as 0 or 1 allows to choose a compound set of trajectories, denotes tensor product operation, and [ρ 0 ] denotes the initial state. The quantum state of the light after diffraction through consecutive N planes includes a superposition of different trajectories through the slits. Experiments for entangled histories has just been, for the first time, performed in Reference [11] by using the polarization states of a single photon and by creating Greenberger-Horne-Zeilinger (GHZ)-type states. MPD-based design compared with complex single photon setup allows the classicality of light sources and simple intensity detection (or photon counting) as a significantly low complexity tool to study quantum histories and QM foundations with near-future experiments. MPD utilizes simple and widely available coherent sources such as Gaussian wave packets of standard laser output conventionally denoted as classical light.
In this article, an important property of MPD-based QPE is, for the first time, presented: Leggett-Garg Inequality (LGI) violations as the temporal analog of Bell's inequality. One of the fundamental tools to analyze quantum temporal correlations of a system is to check the violations of LGIs [12]. LGIs, as proposed by Leggett and Garg in 1985, check a system in terms of the fundamental principles of macroscopic realism (MR) and noninvasive measurability (NIM) such that the systems obeying these rules satisfy the intuition about the classical macroscopic world [13]. QM systems violate LGIs such that MR principles implying the existence of a preexisting value of a macroscopic system and the NIM principle implying the measurement of the value without disturbing the system are both invalidated [14,15].
LGI violations [12][13][14]16,17] are utilized for various purposes such as testing temporal correlations of a single system as an indicator of the quantumness and analyzing QC systems, e.g., Grover's algorithm violating temporal Bell inequality [18]. The simple LGI inequality with three-time formulation violated with various QM setups is defined as follows: C 01 + C 12 − C 02 ≤ 1 (2) where C ij ≡ Q i Q j is the expected value of the multiplication of the dichotomic observables Q i as the measurement outcomes at time t i . The left-hand side is maximally violated by QM systems with the value of 3 / 2. The violation analysis of LGIs is, for the first time, performed for MPD by utilizing the recently proposed ambiguous form by Emary in Reference [16] with the precautions regarding the signaling-in-time (SIT) problem in order to convince a macrorealist about the noninvasive nature of measurements, i.e., to prevent signaling forward in time with measurements. This is achieved by inferring event probabilities from ambiguous measurements rather than direct measurements and by modifying the fundamental inequality in Equation (2)  LGIs. However, MPD utilizes photon-counting intensity detection with a significantly low experimental complexity. It is also simpler compared with the LGI violating architectures utilizing single-photon sources and Mach-Zehnder interferometers [11,20,21]. Besides that, light sources not fully coherent in terms of spatial and temporal dimensions are theoretically modeled while the violation of LGI and QPI are numerically analyzed for specific MPD setup geometry satisfying coherence of light under Gaussian source beam assumptions. On the other hand, LGIs are interpreted in a quantum contextual framework in Reference [22], where the contextuality implies the impossibility to consider a quantum measurement as revealing a preexisting property independent of the set of measurements. It is also analyzed in relation with consistent histories approach in Reference [23]. Furthermore, nonlocality and contextuality are presented as important quantum resources [24]. Therefore, the relation of the proposed QPE and QPI resources with quantum contextuality is an open issue to be explored.
The other important property of quantum histories is the interference among them denoted by QPI. To the best of the author's knowledge, theoretical modeling of the interference among quantum history states leading to a counterintuitive observation to be easily verified experimentally has not been previously formulated. Implementation of the theoretically modeled QPI setup will significantly improve our understanding about QM fundamentals regarding time. QPI is the temporal analogue of the spatial interference obtained in Young's double-slit setup. Destructive and constructive interferences among the paths are observed in the time domain for the QPI case. A special case is modeled such that decreasing the number of photons to diffract through a plane by removing a Feynman path results in an increase in the number of photons diffracting through the next plane due to the interference between two quantum trajectories. This is proposed, for the first time, as a counterintuitive nature of the interference among the quantum histories.
The novel contributions of the article are summarized as follows: • The paper is organized as follows. We firstly define MPD setup with diffractive projection and measurement operators in Sections 2.1 and 2.2. It is followed by the history state modeling of QPE in Section 2.3. Then, we present theoretical modeling of the violation of LGI in Section 2.4, followed by QPI scenario in Section 2.5. Then, numerical analysis is presented in Section 2.6. We provide the conclusions and discuss future applications of QPE and QPI based on MPD setup in Section 3. Finally, the methods utilized for theoretical modeling are presented in Section 4.

MPD Setup for Quantum Temporal Correlations
MPD setup is formed from N − 1 diffraction planes of multiple slits in front of a classical light source and the measurement of interference pattern with N sensor planes, i.e., both diffraction and sensing on the same plane, as shown in Figure 1a. It is also possible to locally count the diffracted photons with the measurement planes inserted between the diffraction planes as discussed in Section 2.6. The utilized light source is assumed to be coherent as the closest analog of a classical light field emphasizing the absence of nonclassical states of light such as single photon generation, squeezed light, or multiple particles of entangled photons [25]. The standard laser output is almost perfectly a coherent state corresponding to the fundamental transverse modes of light field distribution producing Gaussian beams. This coherent Gaussian wave function keeps the position and momentum uncertainties stationary as emphasized by Glauber [26]. It is an eigenstate of the annihilation operatorâ for the harmonic oscillator, i.e.,â |α = α |α , represented as follows in the complete orthonormal basis of the number states |n of the single mode oscillator [26]: where its representation in the position basis gives the Gaussian form. Therefore, the source is assumed to have normalized Gaussian wave function Each plane is assumed to be capable of performing measurement with photodetectors for counting the number of photons hitting the detector area. Therefore, a plane either allows projective diffraction of light through slits denoted by the operator symbol P or performs measurement denoted by M on its sensor array positions where there are no slits. Gaussian slits are utilized with FPI modeling for simplicity [2,8] as mathematically described in Equation (7) in the next subsection. Light is assumed to perform free space propagation between consecutive planes. The plane with the index j has S j slits, where the central positions and widths of slits are denoted by X j,i and W j,i , respectively, and j ∈ [1, N − 1] and i ∈ [1, S j ]. The widths of the slits are assumed to be the same on each plane but not constrained among different planes. Distance between the ith and jth planes is denoted by L i,j , where the distance from the light transmitter source to the first plane is given by L 0,1 . Light is assumed to have propagation in the z-axis with the velocity given by c, while quantum superposition interference is observed in the x-axis as a one-dimensional model which can be easily extended to two dimensions (2D) [8]. Interplane distances and durations are denoted by the vectors L T = [L 0,1 ... L N−1,N ] and t T = [t 0,1 ... t N−1,N ] ≡ L T / c, respectively, where transpose is denoted by (.) T . The value t T is accurate with the assumption L j−1,j W j,i , X j,i for j ∈ [1, N − 1] and i ∈ [1, S j ] such that QM effects are emphasized in the x-axis. Nonrelativistic modeling is assumed. We do not consider the effects of environment dephasing or decohering of the interference pattern for double-slit setups [27,28]. Furthermore, minor effects of exotic paths [29] on the numerical results are ignored as discussed in Reference [8] without affecting the main modeling.
Free-particle evolution kernel for the optical propagation paths between time-position values (t j , x j ) and (t j+1 , x j+1 ) is defined as follows [9,10] with the same form for electron propagation [2,8]: where ∆t = t j+1 − t j , ∆x = x j+1 − x j , m ≡h k / c is the virtual mass term for the photon with the wave number k = 2 π / λ, and λ is the wavelength of the light. The validity of Fresnel diffraction formulation for quantum optical propagation is verified based on recent experimental [30] and theoretical [31] studies, while Fourier optics [32] extension of MPD is recently proposed in Reference [9]. Therefore, the Fresnel diffraction integral for free space proposed in Equation (4) and its consecutive application with FPI formalism are theoretically valid and highly reliable for the simple design of MPD. The proposed theoretical model significantly promises to be verified with near-future experiments due to the simplicity of the setup. Then, the propagated wave function |Ψ j = ∞ −∞ dx j |x j Ψ j (x j ) on the jth plane becomes as follows by utilizing Equation (4) consecutively in FPIs [8][9][10]: where ψ j,n (x j ) is the contribution for each nth propagation path through the slits on the overall superposition, and the definitions of the notations n and N j are explained next while Υ j = χ 0 ∏ j−1 l=1 √ ξ l ; the constants A j−1 , B j−1 , χ 0 , and ξ l for l ∈ [1, j − 1]; H j−1 = H R,j−1 + ı H I,j−1 ; and the vectors − → c j−1 and − → d j−1 depending on the group of {h, m, σ 0 , t l,l+1 , and β l } for l ≤ j − 1 are explicitly defined in Reference [8]. Explicit forms of the parameters required for double-and triple-plane setups are provided in Section 4 while formulating LGIs and QPI, respectively, in the following discussions. The total number of paths just before diffraction on the jth plane is calculated by N j = ∏ j−1 l=1 S l , while the set of slit positions for the path indexed with n ∈ [0, N j − 1] is denoted by − → x n ≡ [X 1,s n,1 X 2,s n,2 ... X j−1,s n,j−1 ] T while each nth path is indexed by the set of diffracted slits as the following: Path n ≡ {s n,1 , s n,2 , ... s n,j−1 ; s n,l ∈ [1, S l ], l ∈ [1, j − 1]} (6) where the specific slit on lth plane for nth trajectory is indexed with s n,l . The same symbol of the position vector − → x is used for both the dimensions N and j. The size of the vector is inferred from the index of the current plane analyzed throughout the text. The position on the jth plane is denoted by x j . In Equation (5) for j = N, each path reaching the Nth plane is indexed by n for n ∈ [0, N p − 1] as shown in Figure 1b for a simple example of N p = 4, where total number of paths is given by multiplying the number of slits on each plane as N p ≡ ∏ N−1 j=1 S j . The vector − → x n ≡ [X 1,s n,1 X 2,s n,2 ... X N−1,s n,N−1 ] T denotes the set of slit positions ordered with respect to the plane indices for nth path for the case of N planes. Next, diffraction and measurement operators are theoretically defined by emphasizing the operator algebra of multiplane evolution.

Diffractive Projection and Measurement Operators
Projection operator denotes the light to be in the Gaussian slit in a coarse-grained sense [8,33] as follows: where Total diffraction through all slits of the jth plane has the Measurement operators are redefined due to the proposed Gaussian slit design such that trace preserving equality is satisfied, i.e., M † j M j + P † j P j = I, where I is the identity operator and (.) † or (.) H denotes Hermitian or conjugate transpose operation. It is assumed that wave function at time t = t 0 evolves to |Ψ j and |Ψ + j for just before and just after diffraction on the jth plane at t − j and t + j , respectively. The state of the light at t + j has experienced either M j or P j . The measurement operator on the jth plane is defined as the following: Therefore, if we define the measurement operator in FPI formalism as multiplication of the wave function with m j (x j ) reducing the probability to measure the light while approaching the slit center, then the following is obtained by using Equation (8): There are two different types of detection mechanisms in MPD design denoted by Rec 1 and Rec 2 . In Rec 1 , all of the planes for j ∈ [1, N] have detectors measuring the incident light and Rec 1 is the model proposed in this article forming a complete set of diffractive projection P j and measurement M j operators until the final detector plane N for j ∈ [1, N − 1] and j ∈ [1, N]. In this article, Rec 1 modeling is utilized to model history-based time evolution of the light. An example is shown in Figure 1b, where there is a total of seven different sets of consecutive events forming a complete set of histories. The proposed setup is modeled compatible with the consistent histories approach defined in Reference [3] or the entangled histories framework in Reference [1]. On the other hand, the receiver type with the sensors only on the final plane is denoted by Rec 2 . In Rec 2 , i.e., the modeling utilized in Reference [8] for QC, only the final intensity distribution or interference pattern on the detector plane is measured. There is either no detection at the time t + N or the light is detected on the final detector plane with the index N. An operator denoting no detection is defined as M o to form a complete set for Rec 2 ; then, Next, consistent histories approach is applied for MPD setup.

History State Modeling of QPE
Following the definition of consistent histories [3][4][5] and entangled histories [1], a history state is defined for MPD based on the set of projections M j and P j,i on each jth plane for j ∈ [1, N] and i ∈ [1, S j ]. History Hilbert space is defined as follows: where H j denotes the set of projections on planes and denotes tensor product operation. Hilbert space until t + j includes both projections P j and M l on the planes with the indices l ≤ j since the light is detected at some plane until t j or still diffracting through the jth plane. A general history state with QPE composed of superposition of trajectories is denoted as follows based on the notation (similar to bra-ket but with different notations of (.| and |.) for the histories corresponding to .| and |. , respectively) in Reference [1]: where |Ψ j ) is some history state between times t 0 and t j for t j > t 0 , the projector O n (t j ) denotes either of M l or P l,i for l ≤ j and i ∈ [1, S l ], and π n as 0 or 1 is some permutation choosing a compound set of histories indexed by n. Observe that t j includes measurements M l for l ≤ j as possible events such that the state does not change after measurement. It also includes events with zero probability such as jth plane projection at times not equal to t j . Some examples for N = 4 are as follows: The state |Ψ a 4 ) shows that the light is detected on the first plane at t 1 while not changing at consecutive time states, i.e., without diffracting even from the first plane. In |Ψ b 4 ), the light is diffracted from the first slit of the first plane at t 1 , then is diffracted from the fourth slit of the second plane at t 2 and from the second slit of the third plane at t 3 , and is finally measured on the fourth plane. The third example |Ψ c 4 ) is a state with zero probability due to the orthogonality of the operators on different planes. A simple example for three planes with two slits is shown in Figure 1b with seven different history states while N p = 4 of them reach the final detector plane as consecutively diffracted trajectories. History Hilbert space summing to the identity denoted by I H as the family based upon an initial state and neglecting the histories with zero probability is described as follows [3]: where M j α denotes α ≡ N + 1 − j consecutive measurements of M j on the same plane. This includes all the possible history states and evolution for the light until t = t N starting from t 0 . A chain operator is presented in Reference [1] to define the inner product between history states which maps a history state to an operator. The chain operator provides history states with positive semi-definite inner products. This operator is inherently defined in the MPD system as the free-particle evolution kernel K(x 1 , t 1 ; x 0 , t 0 ). Assume that the free-particle evolution operator with the notation U j+1,j acts as the bridging operator connecting projections at times t j and t j+1 . Then, chain operator denoted by χ t j+1 ,t j for the time duration (t j , t j+1 ) is defined as follows: where O n (t j+1 ) and O n (t j ) in 4 = denote the cases which are not presented in the first three definitions. I is the identity operator equalizing the consecutive measurements on the same plane, i.e., M l j = M j , for any integer l. Furthermore, it bridges dynamically not possible history states which have zero probability to occur as discussed in Reference [3]. These include consecutive measurements on different planes such as M j+1 M j , future projection or measurements at a previous time such as M j P j+1 , or consecutive sets of the same projector P j at future times such as P j P j , where free-space propagation in the z-axis prevents this. Then, the compound history state mapped or affected by the chain operator is defined as follows: where V j+1,j denotes either U j+1,j or I. Besides that, MPD allows to model and explore varying kinds of superposition of history states and QPEs similar to the specific entangled states discussed in Reference [1] resembling the temporal counterpart of Bell states. For example, entangled history states of the GHZ type is experimentally tested in Reference [11]. It is an open issue to utilize MPD to generate and test such states with important implications and applications based on QPE. Next, probability amplitudes of histories are modeled.

Event Probabilities
The probabilities characterize the statistical properties of the measurement of classical light. It is assumed that the probability is proportional to the square of the wave function with Born's postulate. It is calculated by integrating the number of photons on the detector area at a specific position for a time interval T enough to obtain the statistical properties [30,31]. The normalized probability is easy to calculate by measuring the number of photons for each event by forming a histogram and then by dividing the number of photons for the specific event to the total number of source photon counts. The number of photons at a particular position x is frequently denoted with the integral element |Ψ(x)| 2 dx, while |Ψ(x)| 2 is denoted as the intensity of the light at the particular position. The probability for the particular history state is found with the positive semi-definite inner product defined as follows: where tr{.} is the trace operation. Assume that two specific elementary history states corresponding to specific diffraction paths indexed with n and n ∈ [0, N p − 1] composing the superposition wave function in Equation (5) are denoted by |ψ N, n ) and |φ N, n ), respectively. These paths include only the diffraction projections at the planes with the indices j ∈ [1, N − 1] denoted by P j,s n,j and P j,s n,j , respectively. If the initial state ρ 0 = |Ψ 0 Ψ 0 | and M N = I are included, then the weight of an elementary diffraction history denoted by the inner product W n ≡ (ψ N, n |ψ N, n ) in Reference [3] becomes the following: where the trace is realized with respect to the position, tr{ρ 2 0 } = 1 is utilized, and ψ N, n (x N ) in position basis of the Nth plane is calculated by putting j = N and n = n in the defined wave function ψ j,n (x j ) in Equation (5). Similarly, inner product between history states is defined as follows: The probability for the light to be diffracted through the ith slit on the jth plane with the projection P j,i is denoted by Prob P j,i . Similarly, probability to be measured on the jth plane with measurement projection M j is denoted by Prob M j . Prob P j,i is calculated by using the weight of the compound history Ω N,{j,i} including the targeted event P j,i as follows: where Ω N,{j,i} is defined as follows: ...
where elementary diffraction history states include diffraction events P l,i l for l < j and i l ∈ [1, S l ] until the jth plane and diffraction event P j,i on the jth plane at t j , and where the events O n (t j+1 ) to [O n (t N )] denote any dynamically possible projector at the times between t j+1 and t N . Probability for the events after diffraction will not have any effect on diffraction probability through P j,i , and those projections are discarded. Then, it is easily calculated by using Equations (5) and (7) and with Ψ j | P † j,i P j,i |Ψ j as follows: Similarly, diffraction through one of several slits in a superposition of s slits on the jth plane is given by the following expression: It is important to emphasize the practical meaning of the probabilities of the light diffractions or measurements on the plane. In practice, the probabilities are proportional to the number of photons for each event, e.g., the number of photons passing through a particular slit integrated over a long measurement time for calculating projection probabilities or the number of photons detected on the specific area of the planes for the measurement projection. Histogram-based modeling for counting the photons for all planes and the slits provides the normalized overall probability for each event summing to a total of unity. Photon or particle counting with classical light is already achieved in various studies characterizing the exotic properties of the paths [31] or Fresnel diffraction properties [30].
Quantumness and temporal correlations for the MPD system design are analyzed by explicitly providing theoretical formulation of LGIs next.

Modeling of the Violation of LGI in MPD
LGIs test the temporal correlations by measuring at different times in analogy with spatial Bell's inequalities for the entanglement between spatially separated systems [12,13]. Three-time correlation-based inequality is defined in Equation (2) as K = C 01 + C 12 − C 02 ≤ 1, where the bound is violated quantum mechanically with dichotomic systems, i.e., Q j = ±1 for j ∈ [0, 2], reaching the bound 3/2 for a two-level system with the maximum LGI violation of 1 / 2. C j 1 ,j 2 ≡ Q j 1 Q j 2 is the expected value of the multiplication of the dichotomic observables, which is equal to C j 1 , is the probability for the measurement of Q j 1 and Q j 2 at times t j 1 and t j 2 , respectively, and t 2 > t 1 > t 0 . Noninvasiveness or non-disturbing structure of the measurements should be clearly satisfied in order to reduce the "clumsiness loophole" [15,16], i.e., experimental limitations and disturbance of the clumsy measurements making it difficult to convince a macrorealist. In Reference [16], ambiguous measurements are utilized to revise Equation (2) by including the effect of signaling. In this article, the same formulation is extended for the MPD setup exploiting simple architecture of slits.
The correlation and entanglement in time are tested with the two-plane setup, where each plane includes triple slits as shown in Figure 2a. It is assumed that the light diffracting through the first plane is taken into account while calculating probability amplitudes, i.e., utilizing negative measurement techniques. For example, if the measured state is set to P 1,1 , then the second and third slits are closed, forcing the light to diffract through only the first slit setting the measurement result. Furthermore, denote p 1 (i 1,1 ) ≡ Prob P 1,i 1,1 and p 1 ({i 1,1 , i 1,2 }) ≡ Prob P 1,i s , where i s = {i 1,1 , i 1,2 } for i 1,1 , i 1,2 ∈ [1, 3] and i 1,1 = i 1,2 . The probability p 1 ({i 1,1 , i 1,2 }) corresponds to the measurement result for P 1,i 1,1 ∪ P 1,i 1,2 being projected in one of the slits with the indices i 1,1 and i 1,2 on the first plane. Similarly, p 1 ({1, 2, 3}) denotes the overall projection on superposition in all three slits. On the other hand, assume that p 1,2 ({i 1,1 , i 1,2 },î 2 ) denotes the probability for the history: where Oˆi 2 (t 2 ) is one of [P 2,1 ], [P 2,2 ], [P 2,3 ], or [M 2 ] denoted byî 2 = 1, 2, 3, and 4, respectively. Similar to Equations (26) and (27), p 1,2 ({i 1,1 , i 1,2 },î 2 ) forî 2 ∈ [1, 3] is found as follows: where the elementary wave function is found with Equation (5) by using i as the path index as follows: where 3], χ 0 , ξ 1 , A 1 , B 1 , H R,1 , H I,1 , c 1 , and d 1 are defined in Section 4. Similarly, p 1,2 ({i 1,1 , i 1,2 }, 4) is defined as follows: where i 1,1 , i 1,2 ∈ [1, 3] and i 1,1 = i 1,2 . p 1,2 (i 1,1 ,î 2 ), and p 1,2 ({1, 2, 3},î 2 ) denote the probabilities for 4]. The same formulation is valid also for the second plane for p 2 (i 2,1 ), p 2 ({i 2,1 , i 2,2 }), and p 2 ({1, 2, 3}) for i 2,1 , i 2,2 ∈ [1, 3] and i 2,1 = i 2,2 . Observe that, at time t 2 , it is assumed that M 2 is also included in calculations providing a complete set I = M 2 + ∑ 3 i=1 P 2,i . Negative measurement methodology for the first plane is utilized such that the light only diffracting through the first plane is utilized in calculating probabilities. Therefore, all the probability calculations based on Equations (26), (27), (29), and (31) are normalized by , p 1,2 (i 1,1 ,î 2 ), p 1,2 ({i 1,1 , i 1,2 },î 2 ), and p 1,2 ({1, 2, 3},î 2 ) are assumed to be normalized through the rest of the article. The normalized operator is defined as P Assume that an ambiguous measurement set of three projections composed by O A The setups for ambiguous measurements are shown in Figure 2b-d, respectively. In addition, an assignment of dichotomic indices for the measurement results is designed denoted by 4]. Q 0 ≡ 1 denotes initial condition [ρ 0 ] with unity probability. These dichotomic indices can be assigned arbitrarily while they are chosen in Section 2.6 based on the maximization of LGI violation by comparing all the possible assignment combinations. Then, utilizing a similar architecture to the ambiguous LGI, i.e., Equation (14) in Reference [16], a conversion matrix D is defined inferring the probability p 1 (i 1,1 ) from the ambiguous measurements with p 1 (i 1,1 ) is the element at the i 1,1 th row and the i A 1 th column of the conversion matrix D, and p 1 (i 1,1 ) denotes the inferred probability such that a macrorealist will not observe any problem. Similarly, p 1,2 (i 1,1 ,î 2 ) becomes the following: where p A 1,2 (i A 1 ,î 2 ) denotes the probability for the history Oˆi 4], while it is found in Equations (29) and (31) by normalizing as follows: where i A ) / 2 since the following probability relation holds: Therefore, D 11 = 0.5, D 12 = 0.5, and D 13 = −0.5. Similarly, D 21 = D 23 = D 32 = D 33 = 0.5 and D 22 = D 31 = −0.5. Then, a macrorealist is convinced that the inferred probabilities are utilized for the calculations of C 01 , C 12 , and C 02 with Equation (2) by replacing p 1,2 (i 1,1 ,î 2 ) with p 1,2 (i 1,1 ,î 2 ) and by properly defining the degree of signaling level between the first and second planes for the ambiguous measurements increasing the required LGI bound. Then, following the similar methodology in Reference [16] (Equations (5) and (14) in Reference [16]), K = C 01 + C 12 − C 02 with Q 0 = 1 is easily transformed into the combination of the standard LGI term free of the invasive measurement and a signaling term as follows by firstly replacing the measured probabilities with the inferred ones and then by inserting into K: where the first term is the standard LGI definition with inferred probabilities and the second term includes inferred signaling terms ∆ S (î 2 ) between the first and second planes for the measurement Oˆi 2 (t 2 ) for eachî 2 . It is modeled as a signaling quantifier showing the influence of the measurement at time t 1 to the measurement at time t 2 and is defined by utilizing ambiguous measurements as follows: where ∑ 3 i 1,1 =1 D i 1,1 ,i A 1 = 1 / 2. Therefore, the no-signaling-in-time (NSIT) condition for the definition with ambiguous measurement of ∆ S (î 2 ) = 0 is expected to convince a macrorealist about the reliability of the measurement setup. Then, the violation of LGI free of the invasive measurement becomes the following: |∆ S (î 2 )| (38) where the summation at the right side of the inequality, i.e., K V − 1, shows the invasiveness of the measurements and the signaling level. Therefore, measured values p A 1,2 (i A 1 ,î 2 ) are utilized to check the violation compatible with respect to the objections of a macrorealist. If Equations (29), (31), and (33) are inserted into Equations (35), (36), and (38), then the following is obtained: ; and the variables f T , G 1 , and G 2 are defined in Section 4. Next, quantum interference among the paths, i.e., QPI, is defined for the MPD setup.

Modeling of QPI
Double-slit interference gives a clear indication of quantumness showing wave-particle duality and spatial interference as emphasized by Feynman. MPD setup presents the complementary phenomenon of the temporal interference among the paths which cannot be explained in any classical way showing that paths interfere in time, destructively and constructively decreasing and increasing the probability of the consecutive events, respectively. A gedanken experiment shown in Figure 3 is designed with three planes. The target is to analyze interference effects of opening both slits on the first plane in terms of the probability of the light to diffract through first (PL-1), second (PL-2), and third (PL-3) planes. History states at times t 1 , t 2 , and t 3 with three types of projections indexed by the superscripts a, b, and c are defined with the setups shown in Figure 3a-c, respectively, as follows: where superposition event at time t 1 is defined as P s ≡ [P 1,1 ] + [P 1,2 ] and the event probabilities are defined as follows: It is observed that p 1 ({1, 2}) = p 1 (1) + p 1 (2) while interference exists among consecutive planes. The targeted scenario for the relation between |Ψ a 3 ), i.e., the superposition of |Ψ b 3 ) and |Ψ c 3 ), and |Ψ c 3 ) is as follows: The superposition of |Ψ b 3 ) and |Ψ c 3 ) on PL-1 increases the probability for the light to diffract at time t 1 in 1 >.
In 2 >, the superposition constructively interferes to increase also the probability to diffract through the slit on PL-2 at time t 2 . However, they destructively interfere in 3 <, decreasing the probability to diffract through the slit on PL-3 at time t 3 . Assuming starting with the setup in Figure 3c, with the second slit (the one with X 1,2 > 0) open on PL-1, if the first slit is additionally opened as shown in Figure 3a, then the probability for the light to diffract through PL-1 and PL-2 increases while decreasing the probability to diffract through PL-3. The counterintuitive observation based on a classical logic with balls passing through the slits is described as follows. We open the second slit on PL-1 and observe that the total number of balls passed through them for a statistical experiment increases. It becomes more probable to pass through PL-1 with two slits in a classically logical manner. Furthermore, the probability to pass through the single slit on PL-2 or the the total number of balls passing through PL-2 somehow increases. However, we observe that the probability to pass through the single slit on the consecutive PL-3 counterintuitively decreases in spite of the fact that more balls are coming from the second plane. This is complementary to the conventional spatial interference extensively studied in double-slit interference experiments. QPI modeling requires the calculation of ψ 3,j (x 3 ) in addition to ψ 2,j (x 3 ) for j ∈ [0, 1] as modeled in Equation (30). The explicit parameters required for the calculation of both ψ 2,j (x 3 ) and ψ 3,j (x 3 ) are presented in Section 4.

Numerical Results
Two different calculation scenarios are performed denoted by Sim 1 and Sim 2 as shown in Table 1. Sim 1 calculates violation of LGI while Sim 2 shows an example of interference in time for the numerical analysis of QPI. Physical parameters are monochromatic laser source wavelength λ = 650 (nm) as a widely available resource, light velocity of 3 × 10 8 (m/s) in the z-direction, andh = 1.05 × 10 −34 (J × s) as Planck's constant. The wavelength allows another degree of freedom to be adapted based on experimental design requirements or the targeted system design. The layouts used in the simulations Sim 1 and Sim 2 are shown in Figure 4a,b, respectively. Furthermore, illustrative measurement setups for practically counting the number of photons compared with the emitted photons in unit time are presented in Figure 4c,d in order to calculate the probabilities p 1 ({1, 2}) and p 1,2 ({1, 3}, 2), respectively.

Violation of LGI
There are two planes (PL-1 and PL-2) as shown in Figures 2 and 4a. The slit positions on PL-1 and PL-2 are set to D s + [−∆x 0 ∆x] × β 1 and [−∆x 0 ∆x] × β 2 , respectively, where the positions on PL-1 are shifted with varying D s and the ratio of inter-slit distances to the slit widths is fixed in both planes. ∆x is chosen larger than seven to realize independence of Gaussian slits, i.e., exp((X 1,i − X 1,i+1 ) 2 / 2 β 2 1 ) 1. The distance between the planes is set to L such that the time duration is t 01 = t 12 ≡ L / c. Gaussian source parameter σ 0 is varied between 10 (µm) and 800 (µm), compatible with standard laser resources including fiber lasers allowing smaller diameters reaching tens of micrometers.
The shift of the slits on PL-1 results in varying levels of temporal correlation for the diffraction through the slits on PL-2.
LGI violation (K A − K V ) is analyzed for varying D s , t 01 = t 02 , β 1 , β 2 , and σ 0 . In Figure 5a, it is shown with the signaling level (K V − 1) for varying D s for t 01 = 0.2 (ns), t 12 = 0.1 (ns), ∆x = 7, β 1 = 15 (µm), β 2 = 30 (µm), and σ 0 = 130 (µm). The maximum violation is analyzed for different values of Q 1,i 1,1 and Q 2,î 2 for i 1,1 ∈ [1, 3] andî 2 ∈ [1, 4], respectively, and the signs maximizing the violation are chosen for each D s shift. In Figure 5b, distributions of the sign assignments maximizing the violation are shown. Different setups realized with varying shift on PL-1 result in different optimized sign assignments for maximum violation. Furthermore, violation decreases as the interplane slit distance increases, i.e., decreasing to zero with K A = K V ≈ 1. It is observed that LGI is violated significantly, reaching close to 0.3 for the specific setup shown in Figure 5a. The signaling ∑ 4 k=1 |∆ S (k)| is close to zero, as shown in Figure 5a for the marked region between the violation peaks. In the next simulations, the signaling is shown to decrease approximately to zero for varying setup values. It is shown as a proof of concept in Figure 6a that the amount of signaling is K V − 1 ≈ 2.3 × 10 −3 for a violation of K A − K V ≈ 0.2122 while satisfying NIM and signaling-in-time-related assumptions discussed in Reference [16] and utilizing a NIM-free violation. In Figure 6, the effects of varying ∆x, σ 0 , and t 01 = t 12 on the violation of LGI are shown for varying β 1 and β 2 pairs. D s and the signs of Q 1,i 1,1 and Q 2,î 2 for i 1,1 ∈ [1, 3] andî 2 ∈ [1, 4], respectively, are chosen to maximize the violation for each pair and specific σ 0 value. β 1 Table 1 which can be further improved by increasing the range of values and the resolution in simulations such as for the values of β 1 , β 2 , and D s . However, the provided simple parameter set shows LGI violation under no-signaling conditions as a proof of concept. More specific observations provide more information about the nature of LGI violations for MPD setup as discussed next. It is observed in Figure 6a that violation becomes smaller as the relative distance between slits compared with the slit width parameter β increases from ∆x = 7 to ∆x = 11 for t 01 = t 12 = 0.1 (ns). Furthermore, σ 0 values maximizing the violation and the functional behavior with respect to varying σ 0 are approximately the same as ∆x changes for fixed t 01 = t 12 . The range of violation reaches ≈ 0.4 and ≈ 0.23 for ∆x = 7 and ∆x = 11, respectively, while decreases as ∆x increases. The signaling term is approximately zero for the ∆x = 11 case with the maximum violation amplitude of K A − K V ≈ 0.2122 as emphasized previously. It is an open issue to design MPD setups maximizing the violation, e.g., similar to the boundary value of 0.5 observed in conventional QM setups or 0.464, as calculated in Reference [16] with inverted measurements for a three-level quantum optical system. Optimum slit width maximizing violation for each σ 0 decreases as ∆x increases, as shown in Figure 6b. In other words, as ∆x increases, β 1 and β 2 are getting smaller to stabilize ∆x × β 1 and ∆x × β 2 for better violation. As a result, increasing the relative inter-slit distance results in more classical behavior while decreasing violation. On the other hand, increasing σ 0 further does not have any effect on the optimum β 1 , β 2 , and D s as the source behaves as a plane wave for the specific setup. Violation amplitudes with respect to different β 1 and β 2 pairs for the maximizing σ 0 values (extracted from Figure 6a) of 30 (µm) and 150 (µm) for ∆x = 7 and ∆x = 11, respectively, are shown in Figure 6c,e, respectively, for t 01 = t 12 = 0.1 (ns). There is a decrease in both the range of β 1 and β 2 values and the maximum violation for ∆x = 11. On the other hand, increasing interplane distance two times, i.e., making t 01 = t 12 = 0.2 (ns), is observed not to change the maximum violation regime of 0.4 while increasing both the value of σ 0 for the maximum violation to 230 (µm), as shown in Figure 6a, and the (β 1 , β 2 ) values giving the maximum violation amplitudes, as shown in Figure 6d. Increasing interplane distance improves the spread of the wave function on the consecutive plane while requires larger widths of source beam and slits in order to have similar violation amplitudes.
The results in Figure 6 assume fully coherent source both temporally and spatially. The realistic modeling of the temporal and spatial coherence of the light source is presented in Section 4.2. In our simulations, the total duration that light propagates is smaller than 1 ns, i.e., corresponding to ≈ 10 −9 × 3 × 10 8 = 30 (cm), which is much smaller than the temporal coherence time of the conventional single-mode lasers, i.e., ∆t > 10 −6 (s) for single-mode fiber lasers with ∆ f of a few KHz [34]. On the other hand, the analysis of the spatial coherence is realized by defining the setup diameters D 1 and D 2 on the first and second planes, respectively, for the areas covering the slits. Then, these are compared with the spatial coherence diameter defined in Equation (58) depending on the duration of the propagation t from the source to the diffraction plane and on the standard deviation of the Gaussian source denoted by σ 0 . The detailed modeling and discussion are provided in Section 4.2, where the spatial coherence diameters D c and diffraction setup diameters D 1 and D 2 are described and it is targeted that D c is larger than both D 1 and D 2 as described in Equations (59) and (60). Then, D c (t 01 , σ 0 ) corresponds to the propagation from the source to the first plane and D c (t 12 , β 1 ) is for the propagation from the first plane to the second plane. It is assumed that the Gaussian slit modifies the diffracted wave with β 1 as the new source standard deviation while leaving the analysis with respect to the parameters of the wave function in Equation (5), e.g., A 1 , as an open issue, as discussed in Section 4.2. Assume that the analysis providing the minimum signaling with ∆x = 11 and t 01 = t 12 = 0.1 (ns) is targeted for coherence. In Figure 6f, the comparisons of D c (t 01 , σ 0 ) vs. D 1 and D c (t 12 , β 1 ) vs. D 2 are shown. It is observed that the spatial coherence diameter covers the simulation parameters for both the peak and no-signaling cases. The LGI violation curve for ∆x = 11 in Figure 6a is plotted again in Figure 6g by also emphasizing the violation amplitudes for the simulation parameters in Figure 6f. It is an open issue to design MPD setups compatible with the coherence properties of the practical sources while satisfying quantum properties including the violation of LGI.
Besides that, similar to the analysis of spatial and temporal coherence properties of the violation of LGI, the parameter set resulting in the maximum constructive and destructive interferences on PL-2 and PL-3, respectively, is analyzed for compatibility with the coherence of practical sources. The inequalities defined in Equations (61)-(63) in Section 4.2 are targeted to be satisfied where the spatial coherence diameter D c and diffraction setup diameters D 1 , D 2 , and D 3 are described. The inequalities are calculated for the defined group of σ 0 , β 1 , β 2 , β 3 , t 01 , t 12 , and t 13 , where X 3,1 is given in Figure 7c and X 2,1 ∈ [140, 170] (µm) is targeted. It is found that the practical coherence properties are satisfied as shown in Figure 7f comparing D c (t 12 , β 1 ) with D 2 and D c (t 23 , β 2 ) with D 3 where D c (t 01 , σ 0 ) ≈ 607 (µm) > D 1 ≈ 270 (µm).

Discussion and Conclusions
In this article, two novel resources for quantum technologies, i.e., QPE and QPI, are introduced based on tensor product structure of quantum history states for the simple linear optical setup of MPD. Operator theory modeling is presented by combining conventional history-based approaches of Griffith [3][4][5] and entangled histories framework of Reference [1] with FPI modeling. The inherent Feynman path generation mechanism of MPD setup is exploited for realizing quantum trajectories and for studying quantum temporal correlations. Following the similar terminology of entangled quantum trajectories in previous formulations, the temporal correlation among the quantum propagation paths of MPD design is denoted with QPE as a novel quantum resource. The state of the light after diffraction through consecutive planes is represented as a superposition of different trajectories through the slits as the main definition of QPE. Its two fundamental properties are theoretically and numerically analyzed: violation of LGI as the temporal analog of Bell's inequality with the ambiguous form and no-signaling recently proposed by Emary [16] and QPI as a counterintuitive phenomenon defining the quantum interference of histories.
LGIs are violated reaching > 0.2 with a NIM-free formulation complementing the recent experimental implementation of the violation of LGI in Reference [19] utilizing linear polarization degree of freedom of the classical light with a photon-counting simple MPD setup of classical light. Besides that, it is an open issue to further design novel structures allowing the implementation of specific QPE states in analogy with GHZ states experimentally implemented for entangled trajectories [11].
Furthermore, QPI is numerically analyzed for a scenario where the decrease in the number of photons diffracting through a plane counterintuitively results in an increase in the number of photons diffracting through the next plane due to interference between two quantum history or trajectory states. MPD design introducing QPE and QPI is providing a test bed to understand the nature of temporal correlations improving our understanding of quantum mechanics and to design novel structures exploiting history-based formulation for practical purposes of quantum technologies. The simplicity of both the setup and proposed theoretical modeling allows varying kinds of gedanken experiments for implementing paradoxes emphasizing the quantum nature with counterintuitive scenarios.
MPD has significantly low experimental complexity, relying on light sources of finite spatial and temporal coherence combined with conventional intensity detection or photon counting by promising near-future implementations. The experimental implementation is further simplified based on the mature science of Fourier optics with recent proposals in Reference [9]. Experimental and simulation studies for single-plane diffraction architectures verify the proposed theoretical modeling based on Fresnel diffraction of Fourier optics [30] and FPIs [31,35,36]. Therefore, proposed MPD theory simply extends previous modeling to multiple planes with FPI approach. One of the experimental challenges is the design of Gaussian slit, while any slit mask can be represented in terms of Gaussian basis as discussed in Reference [9]. Furthermore, based on the results in this article, any slit mask is potentially expected to have unique MPD setup parameters resulting in both the violation of LGI and QPI as a future work both to theoretically design and to experimentally verify. Moreover, these experimental implementations require only the calculation of total photon counts on the total area of the specific slits. This simplifies the slit design by allowing various geometries and the experimental setup regarding the detectors related to the size, position, precision, and noise [37]. The overall number of photon counts is satisfied more easily compared with detailed representations of the interference pattern in spatial basis.
QPE as an entanglement resource promises future applications for computing [8,9], communications [10], and varying quantum technologies while enriching the conventional resources of quantum entanglement based on correlations among multiple spatial quantum units. In addition, MPD-based formulation of QPE and QPI allows future studies emphasizing the importance of time in QM fundamentals such as regarding entanglement in time [7,38,39], quantum cosmology, and gravity [40][41][42].

Parameters for FPI Modeling of the Violation of LGI
LGI and wave functions are modeled with FPI modeling [8], and the resulting parameters to be utilized in Equations (5), (30), (39), and (40) are provided in Table 2. The functions utilized in Equations (39) and (40) are defined as follows, where the variables k i for i ∈ [1,11] defined in Table 2 depend on the system setup parameters t 01 , t 12 , β 1 , β 2 , and σ 0 , where f 1 (x, l) ≡ e k 10 x 2 + l T e 4 , f 2 (y) ≡ −2 1 T 3 g y − ∑ 3 i=1 e 1 (X 1,i , y), f 1,2 (x, y, l) ≡ l T g y + e 1 (x, y), Besides that, the explicit parameters required for the calculation of ψ 3,i (x 3 ) for i ∈ [0, 1], i.e., two different paths through two different slits on the first plane and the single slit both on the second and the third planes, are shown in Table 3 based on the iterative modeling in Reference [8]. Table 2. Parameters for modeling LGI and path integrals ψ 2,i (x 2 ) for i ∈ [0, 2] : three paths .

Formula
Formula Formula

Temporal and Spatial Coherence of the Light Sources
The coherence of the light field is characterized both temporally and spatially [43]. Its temporal coherence length is inversely proportional to the the full width half maximum (FWHM) of the emission peak in the wavelength spectrum as follows [43,44]: where n r is the refractive index of the medium and ∆λ is the wavelength domain representation of the spectral width. The interference fringes in a double-slit experiment performed with this source are observed if the path length difference between different paths through the slits is smaller than ∆l as described in Reference [43]. Spatial coherence depends on the size of the light field source where it is experimentally characterized with double-slit interference experiment as clearly described in Reference [43], as shown in Figure 8a. The source is assumed to have a side length ∆s, and ∆θ is the angle from the slits separated with D c to the center of the source. The distance between the source and slit planes is L. Then, spatial coherence distance D c satisfies the following relation for the interference fringes to be observed: In Reference [45], more detailed calculation of the spatial coherence diameter of the 2D Gaussian intensity source I(x 0 , y 0 ) = I 0 exp − (x 2 0 + y 2 0 ) / (2 σ 2 0 ) is achieved by using van Cittert-Zernike theorem [46]. It is observed that an incoherent source of uniform intensity with circular diameter ∆s results in D c ∝ λ L / (π ∆s) while the coherence diameter of an incoherent source with Gaussian-distributed intensity and the standard deviation of a propagating coherent Gaussian source are both some multiple of λ L / (π σ 0 ). Therefore, in this article, it is assumed that D c is given by the beam width of the propagating Gaussian wave, which is approximated as 2 √ 2 σ D for the far field, as shown in Figure 8b, where the intensity drops to 1/e 2 at D c / 2. The accuracy of the estimation is further improved to include near field due to the diverse parameter ranges of the LGI and time domain interference setups as follows: where σ D ≡ 1 / √ −2 A 0 for A 0 < 0 since the free-space propagating Gaussian wave function on the detector plane Ψ D (x) is proportional to exp(A 0 x 2 + ı B 0 x 2 ), A 0 = −m 2 σ 2 0 / (2h 2 L 2 / c 2 + 2 m 2 σ 4 0 ), and t ≡ L / c based on trivial application of FPI kernel in Equation (4). The slit positions on the first plane should be inside the coherence area for the proposed numerical simulations to be more compatible with the future experimental implementations.
In numerical simulations of Sim 1 for the violation of LGI with two planes of triple slits, the position of the furthest edge of a slit on PL-1 is found by D s + ∆x β 1 + √ 2 β 1 for D s > 0 by assuming that the Gaussian slit has a width of 2 √ 2 β 1 similar to the Gaussian source formulation. Then, for the specific set of (σ 0 , D s , β 1 , ∆x) shown in Figure 6a,b, the simulation results are reliable if the slits are in the spatial coherence area formulated as follows: D c (t 01 , σ 0 ) ≥ D 1 ≡ 2 D s + (∆x + where D 1 corresponds to the diameter of the area on PL-1, where the slits reside as shown in Figure 8c. (c) (d) Figure 8. (a) The conventional modeling for the spatial coherence of light sources based on double-slit diffraction [43], where ∆θ ∆s ≤ λ is required for the fringes to be observed determining the spatial coherence diameter (D c ); (b) free-space propagation of Gaussian beam, where D c is approximated as the 1/e 2 intensity beamwidth of 2 √ 2 σ 0 with the standard deviation of σ D . The descriptions of the calculation of the setup diameters on the planes to include the slits are denoted by D j for j ∈ [1, 3] with respect to the location and the standard deviation of the source on the previous plane (σ 0 for the first plane and β j−1 for the jth plane) for (c) LGI violation numerical analysis Sim 1 with two planes of triple slits on each plane and (d) interference in time scenario Sim 2 with three planes.
Similarly, the maximum difference between the positions of the slits on PL-1 and PL-2, i.e., the upper slit on PL-1 and the lower slit on PL-2, is calculated by using D 2 ≡ 2 D s + (∆x + √ 2) β 1 + 2 (∆x + √ 2) β 2 . It is assumed that the projected light through the slit on the (j − 1)th plane has spatial coherence larger than D c (t j−1,j , β j−1 ) by assuming that the slit masking modifies the beam width to ≈ 2 √ 2 β j−1 . More accurate analysis left as an open issue requires formulation of the coherence for the wave function ψ j,n (x j ) for each nth path in Equation (5) by also considering the parameter A j for calculating the beam width on the jth plane instead of utilizing −m 2 β 2 j−1 / (2h 2 L 2 / c 2 + 2 m 2 β 4 j−1 ) calculated by replacing σ 0 with β j−1 in the expression of A 0 . Then, the validity of the simulation requires the following inequality between PL-1 and PL-2: D c (t 12 , β 1 ) ≥ D 2 (60) Similar formulations for Sim 2 with three planes as shown in Figure 8d for the time domain interference result in the following requirements for the results in Figure 7 to be compatible with the sources practically not fully coherent: Funding: This research was funded by Ozyegin University Research Grant.

Conflicts of Interest:
The author declares no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.