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.

1. 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: where is the projection operator for diffraction through the slit indexed with on jth plane and for nth trajectory, as 0 or 1 allows to choose a compound set of trajectories, ⊙ denotes tensor product operation, and 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: where is the expected value of the multiplication of the dichotomic observables as the measurement outcomes at time . The left-hand side is maximally violated by QM systems with the value of . 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) by including a signaling term and by providing a NIM-free bound as described in detail in the Results section. The violation of LGI with no-signaling assumption reaching >0.2, i.e., left-hand side of >1.2, is numerically obtained for three-time formulation of LGI in MPD setup. The optimization study to maximize it to the calculated bounds [16] is left as an open issue. Besides that, a novel system design, i.e., MPD, violating LGIs with classical light sources is proposed in this article, complementing the recent experimental result in Reference [19] utilizing linear polarization degree of freedom of the classical light to violate 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:

introduction and operator theory modeling of two novel quantum resources, i.e., QPE and QPI, denoting temporal correlations and the interference among quantum trajectories, respectively, in MPD while utilizing the tensor product structure for future quantum computing and communication architectures and foundational QM studies;

operator theory modeling of MPD-based resources QPE and QPI by combining history-based previous formulations of quantum histories [1,3,4,5] with FPI formalism;

theoretical modeling and numerical analysis of MPD setup for the violation of LGI, with the ambiguous and no-signaling forms recently proposed by Emary in Reference [16], reaching of correlation amplitude numerically obtained for three-time formulation while leaving the maximization of the violation to the boundary levels as an open issue;

a novel setup, i.e., MPD, violating the ambiguous form of LGI with classical light sources complementing the recent experiment utilizing linear polarization degree of freedom of the classical light [19] while MPD setup with remarkably low complexity design utilizing classical light sources and photon-counting intensity detection;

theoretical modeling and numerical analysis of counterintuitive properties and examples of the interference among MPD-based Feynman paths denoted as QPI promising to be easily verified experimentally in future studies;

the modeling and numerical analysis of the coherence properties of the light sources in terms of spatial and temporal dimensions while discussing design issues for MPD setup with coherent light sources; and

discussion for future applications of QPE and QPI as quantum resources and experimental implementations.

The paper is organized as follows. We firstly define MPD setup with diffractive projection and measurement operators in Section 2.1 and Section 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.

2. Results

2.1. MPD Setup for Quantum Temporal Correlations

MPD setup is formed from 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 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 with the standard deviation term .

Figure 1

(a) System model of the free propagating light with velocity c in the z-direction and MPD through N planes, where jth plane includes slits at positions for and interplane distance of . (b) Example of three plane diffractions () with two slits for the first and second planes showing all the possible seven types of histories composed of diffractions or projections , , , and through slits and measurements , , and on the planes. There are paths detected on the third plane.

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 or performs measurement denoted by 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 slits, where the central positions and widths of slits are denoted by and , respectively, and and . 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 , where the distance from the light transmitter source to the first plane is given by . 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 and , respectively, where transpose is denoted by . The value is accurate with the assumption for and 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 and is defined as follows [9,10] with the same form for electron propagation [2,8]: where , , is the virtual mass term for the photon with the wave number , 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 on the jth plane becomes as follows by utilizing Equation (4) consecutively in FPIs [8,9,10]: where is the contribution for each nth propagation path through the slits on the overall superposition, and the definitions of the notations n and are explained next while ; the constants , , , and for ; ; and the vectors and depending on the group of , m, , , and for 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 , while the set of slit positions for the path indexed with is denoted by while each nth path is indexed by the set of diffracted slits as the following: where the specific slit on lth plane for nth trajectory is indexed with . The same symbol of the position vector 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 . In Equation (5) for , each path reaching the Nth plane is indexed by n for as shown in Figure 1b for a simple example of , where total number of paths is given by multiplying the number of slits on each plane as . The vector 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.

2.2. 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 is the slit projection function and the effective slit width is , i.e., leading to a drop in the intensity, where and . Projectors are mutually exclusive with high accuracy such that slit distances are chosen large enough to satisfy for . Total diffraction through all slits of the jth plane has the operator . Measurement operators are redefined due to the proposed Gaussian slit design such that trace preserving equality is satisfied, i.e., , where is the identity operator and or denotes Hermitian or conjugate transpose operation. It is assumed that wave function at time evolves to and for just before and just after diffraction on the jth plane at and , respectively. The state of the light at has experienced either or . 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 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 and . In , all of the planes for have detectors measuring the incident light and is the model proposed in this article forming a complete set of diffractive projection and measurement operators until the final detector plane N for and . In this article, 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 . In , 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 or the light is detected on the final detector plane with the index N. An operator denoting no detection is defined as to form a complete set for ; then, . Next, consistent histories approach is applied for MPD setup.

2.3. 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 and on each jth plane for and . History Hilbert space is defined as follows: where denotes the set of projections on planes and ⊙ denotes tensor product operation. Hilbert space until includes both projections and on the planes with the indices since the light is detected at some plane until 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 is some history state between times and for , the projector denotes either of or for and , and as 0 or 1 is some permutation choosing a compound set of histories indexed by n. Observe that includes measurements for 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 . Some examples for are as follows: The state shows that the light is detected on the first plane at while not changing at consecutive time states, i.e., without diffracting even from the first plane. In , the light is diffracted from the first slit of the first plane at , then is diffracted from the fourth slit of the second plane at and from the second slit of the third plane at , and is finally measured on the fourth plane. The third example 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 of them reach the final detector plane as consecutively diffracted trajectories. History Hilbert space summing to the identity denoted by as the family based upon an initial state and neglecting the histories with zero probability is described as follows [3]: where denotes consecutive measurements of on the same plane. This includes all the possible history states and evolution for the light until starting from . 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 . Assume that the free-particle evolution operator with the notation acts as the bridging operator connecting projections at times and . Then, chain operator denoted by for the time duration (, ) is defined as follows: where and in denote the cases which are not presented in the first three definitions. is the identity operator equalizing the consecutive measurements on the same plane, i.e., , 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 , future projection or measurements at a previous time such as , or consecutive sets of the same projector at future times such as , 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 denotes either or .

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 , while 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 is the trace operation. Assume that two specific elementary history states corresponding to specific diffraction paths indexed with and composing the superposition wave function in Equation (5) are denoted by and , respectively. These paths include only the diffraction projections at the planes with the indices denoted by and , respectively. If the initial state and are included, then the weight of an elementary diffraction history denoted by the inner product in Reference [3] becomes the following: where the trace is realized with respect to the position, is utilized, and in position basis of the Nth plane is calculated by putting and in the defined wave function 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 is denoted by . Similarly, probability to be measured on the jth plane with measurement projection is denoted by . is calculated by using the weight of the compound history including the targeted event as follows: where is defined as follows: where elementary diffraction history states include diffraction events for and until the jth plane and diffraction event on the jth plane at , and where the events to denote any dynamically possible projector at the times between and . Probability for the events after diffraction will not have any effect on diffraction probability through , and those projections are discarded. Then, it is easily calculated by using Equations (5) and (7) and with as follows: is calculated with . Similarly, diffraction through one of several slits in a superposition of s slits on the jth plane is given by the following expression: where and for , .

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.

2.4. 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 , where the bound is violated quantum mechanically with dichotomic systems, i.e., for , reaching the bound for a two-level system with the maximum LGI violation of . is the expected value of the multiplication of the dichotomic observables, which is equal to , where is the probability for the measurement of and at times and , respectively, and . 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 , then the second and third slits are closed, forcing the light to diffract through only the first slit setting the measurement result. Furthermore, denote and , where for and . The probability corresponds to the measurement result for being projected in one of the slits with the indices and on the first plane. Similarly, denotes the overall projection on superposition in all three slits. On the other hand, assume that denotes the probability for the history: where is one of , , , or denoted by , and 4, respectively. Similar to Equations (26) and (27), for is found as follows: where the elementary wave function is found with Equation (5) by using i as the path index as follows: where , , , , , , , , , , , and are defined in Section 4. Similarly, is defined as follows: where and . , and denote the probabilities for and , where and . The same formulation is valid also for the second plane for , , and for and . Observe that, at time , it is assumed that is also included in calculations providing a complete set . 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 . The probabilities denoted by , , for , , , and are assumed to be normalized through the rest of the article. The normalized operator is defined as for .

Figure 2

(a) The violation of Leggett–Garg Inequality (LGI) with the setup of two planes with triple slits where the event set at time is , , and and, at time , are , , , and and ambiguous measurement setups by closing (b) the third, (c) the second, and (d) the first slits on the first plane.

Assume that an ambiguous measurement set of three projections composed by , , and is defined. 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 and for and , respectively, where and . denotes initial condition 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 is defined inferring the probability from the ambiguous measurements with , where denotes the probability for the history for , is the element at the th row and the th column of the conversion matrix , and denotes the inferred probability such that a macrorealist will not observe any problem. Similarly, becomes the following: where denotes the probability for the history , where and , while it is found in Equations (29) and (31) by normalizing as follows: where and correspond to the the event , i.e., equals , , and for , 2, and 3, respectively. For example, for the proposed setup, since the following probability relation holds: Therefore, , , and . Similarly, and . Then, a macrorealist is convinced that the inferred probabilities are utilized for the calculations of , , and with Equation (2) by replacing with 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]), with 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 between the first and second planes for the measurement for each . It is modeled as a signaling quantifier showing the influence of the measurement at time to the measurement at time and is defined by utilizing ambiguous measurements as follows: where . Therefore, the no-signaling-in-time (NSIT) condition for the definition with ambiguous measurement of 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: where the summation at the right side of the inequality, i.e., , shows the invasiveness of the measurements and the signaling level. Therefore, measured values 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: where ; ; ; the functions , , , and ; and the variables , , and are defined in Section 4. Next, quantum interference among the paths, i.e., QPI, is defined for the MPD setup.

2.5. 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 , , and 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 is defined as and the event probabilities are defined as follows: It is observed that while interference exists among consecutive planes. The targeted scenario for the relation between , i.e., the superposition of and , and is as follows: The superposition of and on PL-1 increases the probability for the light to diffract at time in . In , the superposition constructively interferes to increase also the probability to diffract through the slit on PL-2 at time . However, they destructively interfere in , decreasing the probability to diffract through the slit on PL-3 at time . Assuming starting with the setup in Figure 3c, with the second slit (the one with ) 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.

Figure 3

Setup for constructive and destructive interferences in time for the probabilities to diffract through each plane showing the history states (a) as the superposition of and , (b) , and (c) . The targeted scenario with classically counterintuitive nature where a specific example of interference pattern (represented as the number of lambs denoting the number of photons for a practical counting experiment) for the cases of (d) two slits on PL-1 both open and (e) only the second slit open. The operation of closing the first slit decreases the number of photons diffracted through PL-2 while counterintuitively increases the number of photons through PL-3 since we classically expect a decrease. This scenario shows the interference of histories at two different time instants for PL-2 and PL-3 with firstly constructive and then destructive effects, respectively.

QPI modeling requires the calculation of in addition to for as modeled in Equation (30). The explicit parameters required for the calculation of both and are presented in Section 4.

2.6. Numerical Results

Two different calculation scenarios are performed denoted by and as shown in Table 1. calculates violation of LGI while shows an example of interference in time for the numerical analysis of QPI. Physical parameters are monochromatic laser source wavelength (nm) as a widely available resource, light velocity of (m/s) in the z-direction, and (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 and 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 and , respectively.

Table 1

Simulation and system parameters.

ID	Property	Value	ID	Property	Value
Sim_1	X→1T	D_s+−Δx0Δx×β_1	Sim_2	X→1T	−44×β_1
	X→2T	−Δx0Δx×β_2		X_2,1 (μm)	0,500
	Δx; D_s	{7,11}; 0,3000 (μm)		X_3,1 (μm)	−600,800
	t_01=t_12 (ns)	{0.1,0.2}		t_01,t_12,t_23 (ns)	0.5, 0.2, 0.1
	β_1, β_2 (μm)	1,50, 1,100		β_1, β_2, β_3 (μm)	25, 35, 45
	σ_0 (μm)	10,800		σ_0 (μm)	200

Figure 4

The layouts used in (a) and (b) , where for , the fixed values of the parameters are (m), (ns), (ns), (ns), (m), (m), and (m) in addition to the fixed values of the slit positions on the first plane. The practical measurement setups to be utilized in future experiments are illustrated for the probabilities (c) and (d) . The measurement planes count the detected number of photons compared with the number of photons emitted by the source in unit time.

2.6.1. Violation of LGI

There are two planes (PL-1 and PL-2) as shown in Figure 2 and Figure 4a. The slit positions on PL-1 and PL-2 are set to and , respectively, where the positions on PL-1 are shifted with varying and the ratio of inter-slit distances to the slit widths is fixed in both planes. is chosen larger than seven to realize independence of Gaussian slits, i.e., . The distance between the planes is set to L such that the time duration is . Gaussian source parameter 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 () is analyzed for varying , , , , and . In Figure 5a, it is shown with the signaling level () for varying for (ns), (ns), , (m), (m), and (m). The maximum violation is analyzed for different values of and for and , respectively, and the signs maximizing the violation are chosen for each 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 . It is observed that LGI is violated significantly, reaching close to for the specific setup shown in Figure 5a. The signaling 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 for a violation of while satisfying NIM and signaling-in-time-related assumptions discussed in Reference [16] and utilizing a NIM-free violation.

Figure 5

(a) LGI violation () and signaling () for varying , where (ns), (ns), , (m), (m), and (m) and (b) the corresponding dichotomic sign assignments for ambiguous measurements maximizing the violation for each .

Figure 6

(a) Maximum LGI violation () and the corresponding amount of signaling () for varying , , and and (b) the corresponding values of , , and maximizing the violation for each assuming fully coherent sources. Maximum violation for varying pairs for fully coherent sources where (c) and (ns) at the maximizing (m), (d) and (ns) at (m), and (e) and (ns) at (m). It is observed that there is a large set of slit pairs and beam width resulting in LGI violation reaching for and for , respectively, while there are local peaks for pairs for all cases. Increasing values expands the pairs for similar values of violations. (f) The comparison of the spatial coherence diameters with the diffraction setup diameters and for the first and second planes, respectively, where the targeted case is and (ns), i.e., analyzed as the red curve in Figure 6a, and (g) the corresponding LGI violation curve plotted again by emphasizing the coherence including the peak points.

In Figure 6, the effects of varying , , and on the violation of LGI are shown for varying and pairs. and the signs of and for and , respectively, are chosen to maximize the violation for each pair and specific value. and are chosen in the sets (m) and (m), respectively, while choosing the maximum violation pairs. Similarly, is chosen in the interval of (m) with the resolution of (m). There are important various observations for the specific simulation constraints in Table 1 which can be further improved by increasing the range of values and the resolution in simulations such as for the values of , , and . 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 to for (ns). Furthermore, values maximizing the violation and the functional behavior with respect to varying are approximately the same as changes for fixed . The range of violation reaches and for and , respectively, while decreases as increases. The signaling term is approximately zero for the case with the maximum violation amplitude of as emphasized previously. It is an open issue to design MPD setups maximizing the violation, e.g., similar to the boundary value of observed in conventional QM setups or , as calculated in Reference [16] with inverted measurements for a three-level quantum optical system. Optimum slit width maximizing violation for each decreases as increases, as shown in Figure 6b. In other words, as increases, and are getting smaller to stabilize and 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 further does not have any effect on the optimum , , and as the source behaves as a plane wave for the specific setup. Violation amplitudes with respect to different and pairs for the maximizing values (extracted from Figure 6a) of 30 (m) and 150 (m) for and , respectively, are shown in Figure 6c,e, respectively, for (ns). There is a decrease in both the range of and values and the maximum violation for . On the other hand, increasing interplane distance two times, i.e., making (ns), is observed not to change the maximum violation regime of while increasing both the value of for the maximum violation to 230 (m), as shown in Figure 6a, and the (, ) 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 (cm), which is much smaller than the temporal coherence time of the conventional single-mode lasers, i.e., (s) for single-mode fiber lasers with of a few KHz [34]. On the other hand, the analysis of the spatial coherence is realized by defining the setup diameters and 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 . The detailed modeling and discussion are provided in Section 4.2, where the spatial coherence diameters and diffraction setup diameters and are described and it is targeted that is larger than both and as described in Equations (59) and (60). Then, corresponds to the propagation from the source to the first plane and is for the propagation from the first plane to the second plane. It is assumed that the Gaussian slit modifies the diffracted wave with as the new source standard deviation while leaving the analysis with respect to the parameters of the wave function in Equation (5), e.g., , as an open issue, as discussed in Section 4.2. Assume that the analysis providing the minimum signaling with and (ns) is targeted for coherence. In Figure 6f, the comparisons of vs. and vs. 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 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.

2.6.2. QPI Analysis

The three-plane setup (PL-1, PL-2, and PL-3) shown in Figure 3 is numerically analyzed for fixed values of (m), (m), (m), (m), (ns), (ns), (ns), and . Sampling value of (m) is utilized in the analysis. Constructive and destructive interferences of two different paths at times and , respectively, are performed by designing the slit positions with respect to spatial constructive and destructive interferences on PL-2 and PL-3, respectively. In Figure 7a, single slit position for PL-2, i.e., , is chosen on the constructive interference regions, where is larger than due to the superposition. Then, for each constructive slit position, the destructive interference regions on PL-3 such that is smaller than are searched while the magnitudes and the corresponding are shown in Figure 7b,c, respectively. It is observed that m maximizes the destructive interference while the corresponding wave function amplitude on PL-3 is shown in Figure 7d, showing the the maximum destructive interference at m. Then, for each pair, probability amplitudes of the histories are shown in Figure 7e with the marked areas where constructive interference occurs on PL-2 but with the destructive interference obtained on PL-3. The conditions with counterintuitive nature in Equations (53)–(55) for interference in time are satisfied. The top two pairs of curves satisfy and due to the constructive interference while the bottom pair of the curves satisfies . In other words, the probability for the light to diffract through the second plane is decreased after blocking the first slit on PL-1 which counterintuitively results in an increase for the diffraction probability through the third plane. As a result, the proposed design of the setup and utilization of spatial interference result in interference of the quantum paths in time for the projection histories with classically counterintuitive probabilistic results.

Figure 7

(a) compared with and for diffraction through the layer PL-2, (b) for varying on PL-3 such that destructive interference is maximized for each with respect to while m maximizes the destructive interference, (c) maximizing the destructive interference for varying , (d) the comparison of and on PL-3 for specific m showing the destructive interference maximized with m, and (e) the marked regions satisfy the counterintuitive scenario in (53)–(55) for varying with the corresponding pair in Figure 7c. Constructive and destructive interferences are observed for diffraction through PL-2 and PL-3, respectively, with different kinds of correlation of the paths at different times as a proof-of-concept numerical simulation of quantum path interference (QPI) in time between the two paths. (f) The comparison of setup diameters on the second and third planes, i.e., and , respectively, with the spatial coherence diameters and , respectively, in the targeted range of (m) in Figure 7e.

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 and diffraction setup diameters , , and are described. The inequalities are calculated for the defined group of , , , , , , and , where is given in Figure 7c and (m) is targeted. It is found that the practical coherence properties are satisfied as shown in Figure 7f comparing with and with where (m) (m).

3. 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 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].

4. Methods

4.1. 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 for defined in Table 2 depend on the system setup parameters , , , , and , where , , , , , , , , , , , , , , , , , and . Besides that, the explicit parameters required for the calculation of for , 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 for three paths).

	Formula		Formula		Formula
k_1	−ℏmt_12(a_t,σ+ℏ^2t_01t_12)2k_11	k_8	−141β12+1β12+d_t,σ	ξ_1	β12mϑ_tβ12m(mσ02+ıb_t)+ıℏt_1,2ϑ_t
k_2	ℏm^3σ02t_12β12+d_t,σ/k_11	k_9	121β12−1β12+d_t,σ	A_1	−β12m^2ℏ^2t0,12+m^2σ02Ξ_1(2α_t,σ,β)
k_3	−β12m^2(a_t,σ+ℏ^2t_01t_12)/k_11	k_10	−1/(β12+d_t,σ)	B_1	β14m^3ℏt_0,1+mℏt_1,2(ℏ^2t0,12+m^2Ξ12)(2α_t,σ,β)
k_4	β12m^4σ02(β12+d_t,σ)/k_11	k_11	β14m^2m^2σ02Ξ_2+b_t+β12m^2β22a_t,σ+2c_t,σ+ℏ^2t122a_t,σ	H_R,1	−m^2β12(b_t+m^2σ04)+c_t,σ(2α_t,σ,β)
k_5	−m^2β12(m^2σ02Ξ_2+b_t)+c_t,σk_11	a_t,σb_tc_t,σ	ℏ^2t012+m^2σ04ℏ^2(t_01+t_12)^2ℏ^2σ02t122	H_I,1	mℏt_1,2a_t,σ+ℏ^2t_0,1t_1,2(2α_t,σ,β)
k_6	−m^22β12(m^2σ02Ξ_2+b_t)/(4k_11)+m^2(−β22a_t,σ−2c_t,σ)/(4k_11)	α_t,σ,β	β14m^2b_t+m^2σ04+2β12m^2c_t,σ+ℏ^2t1,22a_t,σ	c_1	β12m^2a_t,σ+ℏ^2t_0,1t_1,2/α_t,σ,β
k_7	β22m^2a_t,σ/(2k_11)	χ_0	π^−1/4mσ_0mσ02+ıℏt_0,1	d_1	−mℏt_1,2ℏ^2t0,12+m^2σ02Ξ_1α_t,σ,β

Table 3

Parameters for modeling the path integrals of QPI for two paths).

	Formula		Formula	j∈[1,2]	Formula
H_2	ν_2,2ζ_1,c+ıζ_1,d^2+ν_1,10ν_3,2ζ_1,c+ıζ_1,dν_1,2	ν_2,2	−β22ℏt_2,32ıς_2	ν_1,j	−2ℏt_j,j+1(A_j−1+ıB_j−1)+ım2ıς_j
		ν_3,2	−ℏt_2,3ıς_2	ζ_j	4B_j−1βj4ℏmt_j,j+1+βj4m^2+ℏ^2tj,j+12ϱ_j
c→_2ν_4,2ζ_1,c+ν_5,2ζ_1,dζ_2,cd→_2ν_4,2ζ_1,d−ν_5,2ζ_1,cζ_2,d		ν_4,2	β22ζ_2,c	ζ_j,c	(2B_j−1ℏmt_j,j+1βj2+βj2m^2)/ζ_j
		ν_5,2	−2ℏt_2,3A_2m	ζ_j,d	ℏmt_j,j+12A_j−1βj2−1/ζ_j
A_0	−m^2σ02/(2ℏ^2t0,12+2m^2σ04)	A_2	β22m^22A_1β22−12ζ_2	ϱ_jξ_j	4βj4Aj−12+Bj−12−4A_j−1βj2+1βj2m/ς_j
B_0	ℏmt_0,1/(2ℏ^2t0,12+2m^2σ04)	B_2	2B_1β24m^2+ℏmt_2,3ϱ_22ζ_2	ς_j	ℏt_j,j+12βj2(B_j−1−ıA_j−1)+ı+βj2m

4.2. 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 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 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 , and is the angle from the slits separated with to the center of the source. The distance between the source and slit planes is L. Then, spatial coherence distance 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 is achieved by using van Cittert–Zernike theorem [46]. It is observed that an incoherent source of uniform intensity with circular diameter results in 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 . Therefore, in this article, it is assumed that is given by the beam width of the propagating Gaussian wave, which is approximated as for the far field, as shown in Figure 8b, where the intensity drops to at . 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 for since the free-space propagating Gaussian wave function on the detector plane is proportional to , , and 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.

Figure 8

(a) The conventional modeling for the spatial coherence of light sources based on double-slit diffraction [43], where is required for the fringes to be observed determining the spatial coherence diameter (); (b) free-space propagation of Gaussian beam, where is approximated as the intensity beamwidth of with the standard deviation of . The descriptions of the calculation of the setup diameters on the planes to include the slits are denoted by for with respect to the location and the standard deviation of the source on the previous plane ( for the first plane and for the jth plane) for (c) LGI violation numerical analysis with two planes of triple slits on each plane and (d) interference in time scenario with three planes.

In numerical simulations of 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 for by assuming that the Gaussian slit has a width of similar to the Gaussian source formulation. Then, for the specific set of shown in Figure 6a,b, the simulation results are reliable if the slits are in the spatial coherence area formulated as follows: where corresponds to the diameter of the area on PL-1, where the slits reside as shown in Figure 8c.

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 . It is assumed that the projected light through the slit on the th plane has spatial coherence larger than by assuming that the slit masking modifies the beam width to . More accurate analysis left as an open issue requires formulation of the coherence for the wave function for each nth path in Equation (5) by also considering the parameter for calculating the beam width on the jth plane instead of utilizing calculated by replacing with in the expression of . Then, the validity of the simulation requires the following inequality between PL-1 and PL-2: Similar formulations for 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: