Theory of Quantum Path Computing with Fourier Optics and Future   Applications for Quantum Supremacy, Neural Networks and Nonlinear   Schr\"odinger Equations

Burhan Gulbahar

arXiv:1908.02274·quant-ph·July 13, 2020

Theory of Quantum Path Computing with Fourier Optics and Future Applications for Quantum Supremacy, Neural Networks and Nonlinear Schr\"odinger Equations

Burhan Gulbahar

PDF

TL;DR

This paper introduces a novel photonic quantum path computing framework using Fourier optics, enabling scalable quantum problem solving without photon interactions, with applications in quantum supremacy, neural networks, and nonlinear Schrödinger equations.

Contribution

It provides a theoretical model for photonic QPC using Fourier optical setups, extending its problem solving capabilities and discussing future quantum computing applications.

Findings

01

QPC can sample over 2^100 paths for quantum supremacy.

02

Theoretical modeling of QPC with Gaussian and Hermite-Gaussian sources.

03

Potential applications include quantum neural networks and nonlinear Schrödinger solutions.

Abstract

The scalability, error correction and practical problem solving are important challenges for quantum computing (QC) as more emphasized by quantum supremacy (QS) experiments. Quantum path computing (QPC), recently introduced for linear optic based QCs (LOQCs) as an unconventional design, targets to obtain scalability and practical problem solving. It samples the intensity from the interference of exponentially increasing number of propagation paths obtained in multi-plane diffraction (MPD) of classical particle sources. QPC exploits MPD based quantum temporal correlations of the paths and freely entangled projections a<t different time instants, for the first time, with the classical light source and intensity measurement while not requiring photon interactions or single photon sources and receivers. In this article, photonic QPC is defined, theoretically modeled and numerically analyzed…

Tables3

Table 1. Table 1: Iteration parameters for FPI modeling of MPD with the kernels K H O subscript 𝐾 𝐻 𝑂 K_{HO} and K L C T ( a , b , c , d ) superscript subscript 𝐾 𝐿 𝐶 𝑇 𝑎 𝑏 𝑐 𝑑 K_{LCT}^{(a,b,c,d)} for Gaussian sources

	Formula for $K_{H O}$ based MPD ( $ω t_{j, j + 1} / π \notin ℤ$ )	Formula for $K_{L C T}^{(a, b, c, d)}$ based MPD ( $b_{j, j + 1} \neq 0$ )
$Ψ_{0} (x_{0}) = exp (- x_{0}^{2} / (2 σ_{0}^{2}))$ $/$ $\sqrt{σ_{0} \sqrt{π}}$ $&$ $Ψ_{1} (x_{1}) = χ_{0} exp (A_{0} x_{1}^{2} + ı B_{0} x_{1}^{2})$
$A_{0}$	$- {\hat{m}}_{0}^{2} σ_{0}^{2} / (2 \cos^{2} (α_{0}) {\hat{m}}_{0}^{2} σ_{0}^{4} + 2 λ_{0}^{2})$	$- 2 π^{2} σ_{0}^{2} / (4 π^{2} a_{01}^{2} σ_{0}^{4} + b_{01}^{2})$
$B_{0}$	$\cos (α_{0}) {\hat{m}}_{0} (λ_{0}^{2} - \sin^{2} (α_{0}) {\hat{m}}_{0}^{2} σ_{0}^{4}) / (2 λ_{0} (\cos^{2} (α_{0}) {\hat{m}}_{0}^{2} σ_{0}^{4} + λ_{0}^{2}))$	$π d_{01} / b_{01} - 4 π^{3} a_{01} σ_{0}^{4} / (b_{01} (4 π^{2} a_{01}^{2} σ_{0}^{4} + b_{01}^{2}))$
$χ_{0}$	$π^{- 1 / 4} \sqrt{{\hat{m}}_{0} σ_{0} / (\cos (α_{0}) {\hat{m}}_{0} σ_{0}^{2} + ı λ_{0})}$	$\exp (- ı π / 4) \sqrt{2 \sqrt{π} σ_{0} / (b_{01} - 2 ı π a_{01} σ_{0}^{2})}$

Table 2. Table 2: Polynomial expressions in ( 18 - 21 ) for the case of N = 3 𝑁 3 N=3 and Gaussian source with σ 0 subscript 𝜎 0 \sigma_{0} ( b j , j + 1 ≠ 0 subscript 𝑏 𝑗 𝑗 1 0 b_{j,j+1}\neq 0 for j ∈ [ 0 , 2 ] 𝑗 0 2 j\in[0,2] )

$p o l_{1}$	$π (- 2 ı b_{01} q_{7} π β_{2}^{2} - 2 π (2 π (a_{01} q_{7} - b_{12} q_{12}) β_{2}^{2} + b_{12} b_{23} ı q_{19}) σ_{0}^{2} + b_{01} b_{12} b_{23} q_{11})$
$p o l_{2}$	$2 π (4 β_{1}^{2} π^{2} (a_{01} q_{7} - b_{12} q_{12}) β_{2}^{2} - a_{01} b_{01} b_{12}^{2} b_{23} + 2 b_{12} ı π q_{17}) σ_{0}^{2} + b_{01} (4 β_{1}^{2} ı π^{2} q_{7} β_{2}^{2} + b_{01} b_{12}^{2} b_{23} (- ı) - 2 b_{12} q_{23} π)$
$p o l_{3}$ , $p o l_{4}$	$- 2 b_{01} b_{12} b_{23} π (2 a_{01} π σ_{0}^{2} + b_{01} ı)$ , $π (b_{01} (b_{01} b_{12} q_{12} - 2 ı β_{1}^{2} q_{7} π) - 2 π (2 π (a_{01} q_{7} - b_{12} q_{12}) β_{1}^{2} + a_{01} b_{01} b_{12} ı q_{12}) σ_{0}^{2})$
$p o l_{5}$	$- 4 β_{2}^{2} b_{01}^{2} b_{12}^{3} b_{23} π^{2} q_{20} (b_{12}^{2} b_{01}^{4} + 4 β_{1}^{2} π^{2} q_{13} b_{01}^{2} + 4 π^{2} σ_{0}^{2} (2 β_{1}^{2} q_{14} b_{01}^{2} + (4 π^{2} q_{28} β_{1}^{2} + a_{01}^{2} b_{01}^{2} b_{12}^{2}) σ_{0}^{2}))$
$p o l_{6}$	$(4 b_{01}^{2} π^{2} q_{11}^{2} β_{1}^{4} + b_{01}^{4} b_{12}^{2} + 4 π^{2} σ_{0}^{2} (2 β_{1}^{2} b_{01}^{2} b_{12}^{2} + (4 π^{2} q_{19}^{2} β_{1}^{4} + a_{01}^{2} b_{01}^{2} b_{12}^{2}) σ_{0}^{2}))$ $\times ((16 β_{1}^{4} π^{4} q_{7}^{2} β_{2}^{4} + b_{01}^{2} b_{12}^{4} b_{23}^{2} + 4 b_{12}^{2} π^{2} q_{8}) b_{01}^{2} + 4 π^{2} σ_{0}^{2} (2 β_{1}^{2} b_{01}^{2} q_{9} b_{12}^{2} + q_{22} σ_{0}^{2}))$
$p o l_{7}$ , $p o l_{10}$	$16 β_{1}^{2} β_{2}^{4} b_{01} b_{12} π^{4} (4 a_{01} π^{2} q_{19} σ_{0}^{4} + b_{01}^{2} q_{11}) q_{29}$ , $ı β_{1}^{2} β_{2}^{2} b_{01} b_{12} σ_{0} (b_{01} - 2 ı a_{01} π σ_{0}^{2}) (b_{01} (b_{01} b_{12} - 2 ı β_{1}^{2} q_{11} π) - 2 π q_{18} σ_{0}^{2})$
$p o l_{8}$	$- 8 β_{1}^{2} β_{2}^{2} b_{01} b_{12}^{2} b_{23} π^{3} (4 a_{01} π^{2} q_{19} σ_{0}^{4} + b_{01}^{2} q_{11}) (b_{12}^{2} b_{01}^{4} + 4 β_{1}^{2} π^{2} q_{13} b_{01}^{2} + 4 π^{2} σ_{0}^{2} (2 β_{1}^{2} q_{14} b_{01}^{2} + (4 π^{2} q_{28} β_{1}^{2} + a_{01}^{2} b_{01}^{2} b_{12}^{2}) σ_{0}^{2}))$
$p o l_{9}$ , $p o l_{12}$	$- 8 β_{2}^{4} b_{01}^{2} b_{12}^{2} π^{3} q_{20} q_{29}$ , $2 β_{2}^{2} b_{12}^{2} π^{2} (- 4 β_{1}^{2} q_{13} π^{2} b_{01}^{2} - 8 β_{1}^{2} π^{2} q_{14} σ_{0}^{2} b_{01}^{2} - 4 (4 π^{4} q_{28} β_{1}^{2} + a_{01}^{2} b_{01}^{2} b_{12}^{2} π^{2}) σ_{0}^{4} - b_{01}^{4} b_{12}^{2})$
$p o l_{11}$	$(2 a_{01} π σ_{0}^{2} + b_{01} ı) (2 π q_{15} σ_{0}^{2} + b_{01} (2 π q_{11} β_{1}^{2} + b_{01} b_{12} ı))$ $\times (b_{01} q_{16} - 2 π (4 β_{1}^{2} π^{2} (a_{01} q_{7} - b_{12} q_{12}) β_{2}^{2} - a_{01} b_{01} b_{12}^{2} b_{23} + 2 b_{12} ı π q_{17}) σ_{0}^{2})$
$p o l_{13}$	$(16 β_{1}^{4} π^{4} q_{7}^{2} β_{2}^{4} + b_{01}^{2} b_{12}^{4} b_{23}^{2} + 4 b_{12}^{2} π^{2} q_{8}) b_{01}^{2} + 4 π^{2} σ_{0}^{2} (2 β_{1}^{2} b_{01}^{2} q_{9} b_{12}^{2} + q_{22} σ_{0}^{2})$
$p o l_{14}$	$4 π^{3} σ_{0}^{2} (σ_{0}^{2} (a_{01}^{2} b_{01}^{2} b_{12}^{4} b_{23}^{2} d_{23} + 16 π^{4} β_{1}^{4} β_{2}^{4} q_{4} (a_{01} q_{7} - b_{12} q_{12}) + 4 π^{2} b_{12}^{2} q_{1}) + 2 β_{1}^{2} b_{01}^{2} b_{12}^{2} (4 π^{2} β_{2}^{2} q_{6} + b_{12}^{2} b_{23}^{2} d_{23}))$ $+ π b_{01}^{2} (16 π^{4} β_{1}^{4} β_{2}^{4} q_{5} q_{7} + b_{01}^{2} b_{12}^{4} b_{23}^{2} d_{23} + 4 π^{2} b_{12}^{2} q_{2})$
The functions $q_{j}$ for $j \in [1, 30]$ utilized while defining the polynomials are defined as follows: $q_{1} \equiv b_{12}^{2} b_{23}^{2} d_{23} β_{1}^{4} - 2 a_{01} b_{12} b_{23}^{2} q_{11} d_{23} β_{1}^{4} + a_{01}^{2} q_{2}$ , $q_{2} \equiv a_{23}^{2} b_{01}^{2} b_{12}^{2} d_{23} β_{2}^{4} - a_{23} b_{01}^{2} b_{12} q_{27} β_{2}^{4} + b_{23} q_{3}$ , $q_{3} \equiv b_{23} d_{23} q_{11}^{2} β_{1}^{4} + 2 β_{2}^{2} b_{01}^{2} b_{23} d_{23} β_{1}^{2} + β_{2}^{4} b_{01}^{2} d_{12} q_{26}$ , $q_{4} \equiv b_{12} (b_{12} - q_{12} d_{23}) - a_{01} b_{12} q_{11} + a_{01} d_{23} q_{7}$ , $q_{5} \equiv - b_{01} b_{23} d_{23} + a_{12} b_{01} q_{25} + b_{12} d_{01} q_{25}$ , $q_{6} \equiv (β_{1}^{2} b_{23}^{2} + β_{2}^{2} q_{12}^{2}) d_{23} - β_{2}^{2} b_{12} q_{12}$ , $q_{7} \equiv - b_{01} b_{23} + a_{12} b_{01} q_{12} + b_{12} d_{01} q_{12}$ , $q_{8} \equiv a_{23}^{2} b_{01}^{2} b_{12}^{2} β_{2}^{4} + 2 a_{23} b_{01}^{2} b_{12} b_{23} d_{12} β_{2}^{4} + b_{23}^{2} q_{24}$ , $q_{9} \equiv 4 π^{2} (β_{1}^{2} b_{23}^{2} + β_{2}^{2} q_{12}^{2}) β_{2}^{2} + b_{12}^{2} b_{23}^{2}$ , $q_{10} \equiv {(b_{12} q_{12} - a_{01} q_{7})}^{2}$ , $q_{11} \equiv a_{12} b_{01} + b_{12} d_{01}$ , $q_{12} \equiv a_{23} b_{12} + b_{23} d_{12}$ , $q_{13} \equiv b_{12}^{2} d_{01}^{2} β_{1}^{2} + 2 a_{12} b_{01} b_{12} d_{01} β_{1}^{2} + (a_{12}^{2} β_{1}^{2} + β_{2}^{2}) b_{01}^{2}$ , $q_{14} \equiv 2 β_{1}^{2} π^{2} β_{2}^{2} + b_{12}^{2}$ , $q_{15} \equiv a_{01} b_{01} b_{12} - 2 ı β_{1}^{2} q_{19} π$ , $q_{16} \equiv - 4 ı β_{1}^{2} q_{7} π^{2} β_{2}^{2} + b_{01} b_{12}^{2} b_{23} ı + 2 b_{12} π q_{23}$ , $q_{17} \equiv a_{01} q_{23} - β_{1}^{2} b_{12} b_{23}$ , $q_{18} \equiv 2 π q_{19} β_{1}^{2} + a_{01} b_{01} b_{12} ı$ , $q_{19} \equiv a_{01} a_{12} b_{01} - b_{12} + a_{01} b_{12} d_{01}$ , $q_{20} \equiv b_{01}^{2} + 4 π^{2} σ_{0}^{2} (β_{1}^{2} + a_{01}^{2} σ_{0}^{2})$ , $q_{21} \equiv b_{12}^{2} b_{23}^{2} β_{1}^{4} - 2 a_{01} b_{12} b_{23}^{2} q_{11} β_{1}^{4} + a_{01}^{2} q_{8}$ , $q_{22} \equiv 16 β_{1}^{4} π^{4} q_{10} β_{2}^{4} + a_{01}^{2} b_{01}^{2} b_{12}^{4} b_{23}^{2} + 4 b_{12}^{2} π^{2} q_{21}$ , $q_{23} \equiv b_{23} q_{11} β_{1}^{2} + a_{23} β_{2}^{2} b_{01} b_{12} + β_{2}^{2} b_{01} b_{23} d_{12}$ , $q_{24} \equiv a_{12}^{2} b_{01}^{2} β_{1}^{4} + b_{12}^{2} d_{01}^{2} β_{1}^{4} + 2 a_{12} b_{01} b_{12} d_{01} β_{1}^{4} + 2 β_{2}^{2} b_{01}^{2} β_{1}^{2} + β_{2}^{4} b_{01}^{2} d_{12}^{2}$ , $q_{25} \equiv b_{23} d_{12} d_{23} + b_{12} (a_{23} d_{23} - 1)$ , $q_{26} \equiv b_{23} d_{12} d_{23} - b_{12}$ , $q_{27} \equiv b_{12} - 2 b_{23} d_{12} d_{23}$ , $q_{28} \equiv q_{13} a_{01}^{2} - 2 β_{1}^{2} b_{12} q_{11} a_{01} + β_{1}^{2} b_{12}^{2}$ , $q_{29} \equiv 4 b_{01}^{2} π^{2} q_{11} q_{7} β_{1}^{4} + b_{01}^{4} b_{12}^{2} q_{12} + 4 π^{2} σ_{0}^{2} (2 β_{1}^{2} b_{01}^{2} q_{12} b_{12}^{2} + q_{30} σ_{0}^{2})$ and finally $q_{30} \equiv 4 π^{2} q_{19} (a_{01} q_{7} - b_{12} q_{12}) β_{1}^{4} + a_{01}^{2} b_{01}^{2} b_{12}^{2} q_{12}$ .

Table 3. Table 3: Iteration parameters for FPI modeling of MPD with the kernel K L C T ( a , b , c , d ) superscript subscript 𝐾 𝐿 𝐶 𝑇 𝑎 𝑏 𝑐 𝑑 K_{LCT}^{(a,b,c,d)} for Hermite-Gaussian sources ( b j , j + 1 ≠ 0 subscript 𝑏 𝑗 𝑗 1 0 b_{j,j+1}\neq 0 )

$Ψ_{0} (x_{0}) = (2^{1 / 4} / \sqrt{W_{0} 2^{l} l!}) exp (- π x_{0}^{2} / W_{0}^{2}) H_{l} (\sqrt{2 π} x_{0} / W_{0})$ and $Ψ_{1} (x_{1}) = χ_{01} \exp (u_{01} x_{1}^{2}) H_{l} (g_{01} x_{1})$
Symbol	Formula	Symbol	Formula	Symbol	Formula
$g_{01}$	$\sqrt{\frac{2 π W_{0}^{2}}{a_{01}^{2} W_{0}^{4} + b_{01}^{2}}}$	$u_{01}$	$- \frac{π W_{0}^{2}}{a_{01}^{2} W_{0}^{4} + b_{01}^{2}} + ı (\frac{π d_{01}}{b_{01}} - \frac{(π W_{0}^{2}) (a_{01} W_{0}^{2})}{b_{01} (a_{01}^{2} W_{0}^{4} + b_{01}^{2})})$	$χ_{01}$	$\frac{2^{1 / 4} \sqrt{W_{0}}}{\sqrt{2^{l} l!}} \sqrt{\frac{a_{01} W_{0}^{2} - ı b_{01}}{a_{01}^{2} W_{0}^{4} + b_{01}^{2}}} {(\frac{a_{01} W_{0}^{2} - ı b_{01}}{\sqrt{a_{01}^{2} W_{0}^{4} + b_{01}^{2}}})}^{l}$

Equations122

 \frac{δ q ( x , t )}{δ t} + \frac{δ ^{2} q ( x , t )}{δ x ^{2}} + 2 κ ∣ q (x, t) ∣^{2} q (x, t) = 0

 \frac{δ q ( x , t )}{δ t} + \frac{δ ^{2} q ( x , t )}{δ x ^{2}} + 2 κ ∣ q (x, t) ∣^{2} q (x, t) = 0

\displaystyle\begin{split}f_{BB}[k]\,\equiv\,I[k]\equiv\Bigg{|}\sum_{s_{1}=1}^{K_{1}}\ldots\sum_{s_{N-1}=1}^{K_{N-1}}\,e^{(\widetilde{A}_{\vec{s}}\,+\,\imath\,\widetilde{B}_{\vec{s}})(k\,T_{s})^{2}}\,\widetilde{\Upsilon}_{\vec{s}}\,e^{\vec{x}_{N-1,\vec{s}}^{T}\,\mathbf{\widetilde{H}}_{\vec{s}}\,\vec{x}_{N-1,\vec{s}}}\,e^{(\vec{\widetilde{h}}_{\vec{s}}^{T}\,\vec{x}_{N-1,\vec{s}})\,k\,T_{s}}\Bigg{|}^{2}&\end{split}

\displaystyle\begin{split}f_{BB}[k]\,\equiv\,I[k]\equiv\Bigg{|}\sum_{s_{1}=1}^{K_{1}}\ldots\sum_{s_{N-1}=1}^{K_{N-1}}\,e^{(\widetilde{A}_{\vec{s}}\,+\,\imath\,\widetilde{B}_{\vec{s}})(k\,T_{s})^{2}}\,\widetilde{\Upsilon}_{\vec{s}}\,e^{\vec{x}_{N-1,\vec{s}}^{T}\,\mathbf{\widetilde{H}}_{\vec{s}}\,\vec{x}_{N-1,\vec{s}}}\,e^{(\vec{\widetilde{h}}_{\vec{s}}^{T}\,\vec{x}_{N-1,\vec{s}})\,k\,T_{s}}\Bigg{|}^{2}&\end{split}

n \sum π_{n} [P_{N - 1, s_{n, N - 1}}] ⊙ [P_{N - 2, s_{n, N - 2}}] ⊙ ... ⊙ [P_{1, s_{n, 1}}] ⊙ [ρ_{0}]

n \sum π_{n} [P_{N - 1, s_{n, N - 1}}] ⊙ [P_{N - 2, s_{n, N - 2}}] ⊙ ... ⊙ [P_{1, s_{n, 1}}] ⊙ [ρ_{0}]

\ket{\Psi_{{\color[rgb]{0,0,0}3}}}\equiv a_{uu}\ket{U_{A}U_{B}}+\,a_{lu}\ket{L_{A}U_{B}}+\,a_{ll}\ket{L_{A}L_{B}}+\,a_{ul}\ket{U_{A}L_{B}}

\ket{\Psi_{{\color[rgb]{0,0,0}3}}}\equiv a_{uu}\ket{U_{A}U_{B}}+\,a_{lu}\ket{L_{A}U_{B}}+\,a_{ll}\ket{L_{A}L_{B}}+\,a_{ul}\ket{U_{A}L_{B}}

\Psi_{{\color[rgb]{0,0,0}3}}(x)=\Psi_{{\color[rgb]{0,0,0}3},1}(x)\,+\,\Psi_{{\color[rgb]{0,0,0}3},2}(x)\,+\ \Psi_{{\color[rgb]{0,0,0}3},3}(x)\,+\,\Psi_{{\color[rgb]{0,0,0}3},4}(x)

\Psi_{{\color[rgb]{0,0,0}3}}(x)=\Psi_{{\color[rgb]{0,0,0}3},1}(x)\,+\,\Psi_{{\color[rgb]{0,0,0}3},2}(x)\,+\ \Psi_{{\color[rgb]{0,0,0}3},3}(x)\,+\,\Psi_{{\color[rgb]{0,0,0}3},4}(x)

K_{L C T}^{(a, b, c, d)} (x_{1}, x_{0}) \approx \frac{1}{b} e^{\frac{-  π}{4}} e^{\frac{ π}{b} (d x_{1}^{2} - 2 x_{1} x_{0} + a x_{0}^{2})}

K_{L C T}^{(a, b, c, d)} (x_{1}, x_{0}) \approx \frac{1}{b} e^{\frac{-  π}{4}} e^{\frac{ π}{b} (d x_{1}^{2} - 2 x_{1} x_{0} + a x_{0}^{2})}

M_{L C T} = [a c b d]

M_{L C T} = [a c b d]

\mbox F r ees p a ce p r o p a g a t i o n : M_{F S}

\mbox F r ees p a ce p r o p a g a t i o n : M_{F S}

\mbox F o u r i er t r an s f or m : M_{F T}

\displaystyle\mbox{Fractional Fourier transform of order $\alpha$: }\mathbf{M}_{FRFT}

\displaystyle\mbox{Thin lens of focal length $f$: }\mathbf{M}_{L}

\mbox S c a l in g : M_{S}

M_{L C T} \equiv [10 \frac{2 π ℏ τ _{b, j}^{⋆}}{m _{λ}} 1] [1 - \frac{1}{λ f _{j, j + 1}} 01] [10 \frac{2 π ℏ τ _{a, j}^{⋆}}{m _{λ}} 1]

M_{L C T} \equiv [10 \frac{2 π ℏ τ _{b, j}^{⋆}}{m _{λ}} 1] [1 - \frac{1}{λ f _{j, j + 1}} 01] [10 \frac{2 π ℏ τ _{a, j}^{⋆}}{m _{λ}} 1]

M_{H O} = [cos (ω t) - \frac{m _{λ} s i n ( ω t )}{2 π ℏ t _{01}} \frac{2 π ℏ t _{01} s i n ( ω t )}{m _{λ}} cos (ω t)]

M_{H O} = [cos (ω t) - \frac{m _{λ} s i n ( ω t )}{2 π ℏ t _{01}} \frac{2 π ℏ t _{01} s i n ( ω t )}{m _{λ}} cos (ω t)]

\displaystyle\begin{split}\Psi_{N,n}^{G}(x_{N})\,=\,&\chi_{0}\,\bigg{(}\prod_{j=1}^{N-1}\chi_{j,n}\bigg{)}\,e^{(\widetilde{A}_{N-1,n}\,+\,\imath\,\widetilde{B}_{N-1,n})\,x_{N}^{2}}\,e^{(C_{N-1,n}\,+\,\imath\,D_{N-1,n})\,x_{N}}\end{split}

\displaystyle\begin{split}\Psi_{N,n}^{G}(x_{N})\,=\,&\chi_{0}\,\bigg{(}\prod_{j=1}^{N-1}\chi_{j,n}\bigg{)}\,e^{(\widetilde{A}_{N-1,n}\,+\,\imath\,\widetilde{B}_{N-1,n})\,x_{N}^{2}}\,e^{(C_{N-1,n}\,+\,\imath\,D_{N-1,n})\,x_{N}}\end{split}

\displaystyle\begin{split}\Psi_{N}^{G}(x_{N})\,\equiv\,{\color[rgb]{0,0,0}\sum_{n=0}^{N_{p}-1}\Psi_{N,n}^{G}(x_{N})\,=\,}\sum_{n=0}^{N_{p}-1}\widetilde{\Upsilon}_{N,n}^{G}\,e^{\vec{x}_{N-1,n}^{T}\,\mathbf{\widetilde{H}}_{N-1,n}^{HO/G}\,\vec{x}_{N-1,n}}\,e^{(\widetilde{A}_{N-1,n}\,+\,\imath\,\widetilde{B}_{N-1,n})\,x_{N}^{2}}\,e^{(\vec{\widetilde{h}}_{N-1,n}^{T}\vec{x}_{N-1,n})\,x_{N}}&\end{split}

\displaystyle\begin{split}\Psi_{N}^{G}(x_{N})\,\equiv\,{\color[rgb]{0,0,0}\sum_{n=0}^{N_{p}-1}\Psi_{N,n}^{G}(x_{N})\,=\,}\sum_{n=0}^{N_{p}-1}\widetilde{\Upsilon}_{N,n}^{G}\,e^{\vec{x}_{N-1,n}^{T}\,\mathbf{\widetilde{H}}_{N-1,n}^{HO/G}\,\vec{x}_{N-1,n}}\,e^{(\widetilde{A}_{N-1,n}\,+\,\imath\,\widetilde{B}_{N-1,n})\,x_{N}^{2}}\,e^{(\vec{\widetilde{h}}_{N-1,n}^{T}\vec{x}_{N-1,n})\,x_{N}}&\end{split}

Ψ_{3}^{G} (x_{3}) = n = 0 \sum 5 Υ_{3}^{G} e^{x_{2, n}^{T} H_{2}^{L C T / G} x_{2, n}} e^{(A_{2} +  B_{2}) x_{3}^{2}} e^{(h_{2}^{T} x_{2, n}) x_{3}}

Ψ_{3}^{G} (x_{3}) = n = 0 \sum 5 Υ_{3}^{G} e^{x_{2, n}^{T} H_{2}^{L C T / G} x_{2, n}} e^{(A_{2} +  B_{2}) x_{3}^{2}} e^{(h_{2}^{T} x_{2, n}) x_{3}}

H_{2}^{L C T / G} = \frac{p o l _{1} ( β _{2} , σ _{0} )}{p o l _{2} ( β _{1} , β _{2} , σ _{0} )} \frac{p o l _{3} ( σ _{0} )}{ p o l _{2} ( β _{1} , β _{2} , σ _{0} )} 0 \vspace 0.1 in \frac{p o l _{4} ( β _{1} , σ _{0} )}{p o l _{2} ( β _{1} , β _{2} , σ _{0} )}

H_{2}^{L C T / G} = \frac{p o l _{1} ( β _{2} , σ _{0} )}{p o l _{2} ( β _{1} , β _{2} , σ _{0} )} \frac{p o l _{3} ( σ _{0} )}{ p o l _{2} ( β _{1} , β _{2} , σ _{0} )} 0 \vspace 0.1 in \frac{p o l _{4} ( β _{1} , σ _{0} )}{p o l _{2} ( β _{1} , β _{2} , σ _{0} )}

h_{2}^{T} = [\frac{p o l _{5} ( β _{1} , β _{2} , σ _{0} )}{p o l _{6} ( β _{1} , β _{2} , σ _{0} )} \frac{p o l _{7} ( β _{1} , β _{2} , σ _{0} )}{p o l _{6} ( β _{1} , β _{2} , σ _{0} )}] +  [\frac{p o l _{8} ( β _{1} , β _{2} , σ _{0} )}{p o l _{6} ( β _{1} , β _{2} , σ _{0} )} \frac{p o l _{9} ( β _{1} , β _{2} , σ _{0} )}{p o l _{6} ( β _{1} , β _{2} , σ _{0} )}]

h_{2}^{T} = [\frac{p o l _{5} ( β _{1} , β _{2} , σ _{0} )}{p o l _{6} ( β _{1} , β _{2} , σ _{0} )} \frac{p o l _{7} ( β _{1} , β _{2} , σ _{0} )}{p o l _{6} ( β _{1} , β _{2} , σ _{0} )}] +  [\frac{p o l _{8} ( β _{1} , β _{2} , σ _{0} )}{p o l _{6} ( β _{1} , β _{2} , σ _{0} )} \frac{p o l _{9} ( β _{1} , β _{2} , σ _{0} )}{p o l _{6} ( β _{1} , β _{2} , σ _{0} )}]

Υ_{3}^{G} = - 2 (- 1)^{3/4} 2 π^{5/4} \frac{p o l _{10} ( β _{1} , β _{2} , σ _{0} )}{p o l _{11} ( β _{1} , β _{2} , σ _{0} )}

Υ_{3}^{G} = - 2 (- 1)^{3/4} 2 π^{5/4} \frac{p o l _{10} ( β _{1} , β _{2} , σ _{0} )}{p o l _{11} ( β _{1} , β _{2} , σ _{0} )}

A_{2} +  B_{2} = \frac{p o l _{12} ( β _{1} , β _{2} , σ _{0} )}{p o l _{13} ( β _{1} , β _{2} , σ _{0} )} +  \frac{p o l _{14} ( β _{1} , β _{2} , σ _{0} )}{b _{23} p o l _{13} ( β _{1} , β _{2} , σ _{0} )}

A_{2} +  B_{2} = \frac{p o l _{12} ( β _{1} , β _{2} , σ _{0} )}{p o l _{13} ( β _{1} , β _{2} , σ _{0} )} +  \frac{p o l _{14} ( β _{1} , β _{2} , σ _{0} )}{b _{23} p o l _{13} ( β _{1} , β _{2} , σ _{0} )}

\displaystyle\begin{split}\Psi_{N,n}^{HG}(x_{N})\,=\,\chi_{01}\bigg{(}\prod_{j=1}^{N-1}\chi_{j,j+1,n}\bigg{)}\,e^{\widetilde{u}_{N-1,N,n}\,x_{N}^{2}}\,e^{\,v_{N-1,N,n}\,x_{N}}\,H_{l}(\widetilde{g}_{N-1,N,n}\,x_{N}\,+\,h_{N-1,N,n})&\end{split}

\displaystyle\begin{split}\Psi_{N,n}^{HG}(x_{N})\,=\,\chi_{01}\bigg{(}\prod_{j=1}^{N-1}\chi_{j,j+1,n}\bigg{)}\,e^{\widetilde{u}_{N-1,N,n}\,x_{N}^{2}}\,e^{\,v_{N-1,N,n}\,x_{N}}\,H_{l}(\widetilde{g}_{N-1,N,n}\,x_{N}\,+\,h_{N-1,N,n})&\end{split}

\displaystyle\begin{split}\Psi_{N,n}^{HG}(x_{N})\,=\,\chi_{01}\,\bigg{(}\prod_{j=1}^{N-1}\widetilde{\chi}_{a,j,j+1,n}\bigg{)}\,e^{\frac{-\,\imath\,\pi\,(N-2)}{4}}\,e^{\vec{x}_{N-1,n}^{T}\,\widetilde{\mathbf{H}}_{N-1,n}^{LCT/HG}\,\vec{x}_{N-1,n}}\hskip 61.42993pt&\\ \times\,e^{(\vec{\widetilde{\gamma}}_{N-1,n}^{T}\,\vec{x}_{N-1,n})\,x_{N}}\,e^{\widetilde{u}_{N-1,N,n}\,x_{N}^{2}}\,H_{l}(\widetilde{g}_{N-1,N,n}\,x_{N}\,+\,\vec{\widetilde{\eta}}^{T}_{N-1,n}\,\vec{x}_{N-1,n})&\end{split}

\displaystyle\begin{split}\Psi_{N,n}^{HG}(x_{N})\,=\,\chi_{01}\,\bigg{(}\prod_{j=1}^{N-1}\widetilde{\chi}_{a,j,j+1,n}\bigg{)}\,e^{\frac{-\,\imath\,\pi\,(N-2)}{4}}\,e^{\vec{x}_{N-1,n}^{T}\,\widetilde{\mathbf{H}}_{N-1,n}^{LCT/HG}\,\vec{x}_{N-1,n}}\hskip 61.42993pt&\\ \times\,e^{(\vec{\widetilde{\gamma}}_{N-1,n}^{T}\,\vec{x}_{N-1,n})\,x_{N}}\,e^{\widetilde{u}_{N-1,N,n}\,x_{N}^{2}}\,H_{l}(\widetilde{g}_{N-1,N,n}\,x_{N}\,+\,\vec{\widetilde{\eta}}^{T}_{N-1,n}\,\vec{x}_{N-1,n})&\end{split}

Ψ_{N}^{H G} (x_{N}) =

Ψ_{N}^{H G} (x_{N}) =

=

Υ_{N}^{H G} e^{u_{N - 1, N} x_{N}^{2}} n = 0 \sum N_{p} - 1 e^{x_{N - 1, n}^{T} H_{N - 1}^{L C T / H G} x_{N - 1, n}} e^{(γ_{N - 1}^{T} x_{N - 1, n}) x_{N}} H_{l} (g_{N - 1, N} x_{N} + η_{N - 1}^{T} x_{N - 1, n})

Υ_{N}^{H G} e^{u_{N - 1, N} x_{N}^{2}} n = 0 \sum N_{p} - 1 e^{x_{N - 1, n}^{T} H_{N - 1}^{L C T / H G} x_{N - 1, n}} e^{(γ_{N - 1}^{T} x_{N - 1, n}) x_{N}} H_{l} (g_{N - 1, N} x_{N} + η_{N - 1}^{T} x_{N - 1, n})

\displaystyle\begin{split}\Psi_{N}^{HG,U}(x_{N})\,=\,\Upsilon_{{\color[rgb]{0,0,0}N}}^{HG}\,e^{u_{N-1,N}\,x_{N}^{2}}\,\sum_{n=0}^{N_{p}-1}\,e^{(\vec{\gamma}_{N-1}^{T}\,\vec{x}_{N-1,n})\,x_{N}}\,e^{\vec{x}_{N-1,n}^{T}\,\mathbf{H}_{N-1}^{LCT/HG}\,\vec{x}_{N-1,n}}\,(H^{\star}(x_{N})\,+2\,\vec{\eta}^{T}_{N-1}\,\vec{x}_{N-1,n})^{l}&\end{split}

\displaystyle\begin{split}\Psi_{N}^{HG,U}(x_{N})\,=\,\Upsilon_{{\color[rgb]{0,0,0}N}}^{HG}\,e^{u_{N-1,N}\,x_{N}^{2}}\,\sum_{n=0}^{N_{p}-1}\,e^{(\vec{\gamma}_{N-1}^{T}\,\vec{x}_{N-1,n})\,x_{N}}\,e^{\vec{x}_{N-1,n}^{T}\,\mathbf{H}_{N-1}^{LCT/HG}\,\vec{x}_{N-1,n}}\,(H^{\star}(x_{N})\,+2\,\vec{\eta}^{T}_{N-1}\,\vec{x}_{N-1,n})^{l}&\end{split}

\Psi_{p,j}(p_{j})=\frac{1}{\sqrt{2\,\pi\,\hbar}}\int dx_{j}\,\Psi_{j}(x_{j})\exp\big(-\imath\,x_{j}\,p_{\color[rgb]{0,0,0}j}\,/\,\hbar\big{missing})

\Psi_{p,j}(p_{j})=\frac{1}{\sqrt{2\,\pi\,\hbar}}\int dx_{j}\,\Psi_{j}(x_{j})\exp\big(-\imath\,x_{j}\,p_{\color[rgb]{0,0,0}j}\,/\,\hbar\big{missing})

W_{j} (x_{j}, p_{j}) = \frac{1}{π ℏ} \int d y Ψ_{j} (x_{j} - y) Ψ_{j}^{*} (x_{j} + y) e^{\frac{ 2 p _{j} y}{ℏ}}

W_{j} (x_{j}, p_{j}) = \frac{1}{π ℏ} \int d y Ψ_{j} (x_{j} - y) Ψ_{j}^{*} (x_{j} + y) e^{\frac{ 2 p _{j} y}{ℏ}}

P_{E} (j) \equiv \int d x_{j} ∣ Ψ_{j} (x_{j}) ∣^{2}

P_{E} (j) \equiv \int d x_{j} ∣ Ψ_{j} (x_{j}) ∣^{2}

N_{p a t h} =

N_{p a t h} =

N_{p a t h} =

N_{p a t h} =

=

Peer Reviews

No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.

Videos

No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.

Full text

11affiliationtext: Department of Electrical and Electronics Engineering, Ozyegin University, 34794 Istanbul, Turkey**affiliationtext: [email protected]

The fully edited version of the published article is available in https://www.nature.com/articles/s41598-020-67364-0 with DOI https://doi.org/10.1038/s41598-020-67364-0. This is the author accepted copy of the original article.

@mathmargin

Theory of Quantum Path Computing with Fourier Optics and Future Applications for Quantum Supremacy, Neural Networks and Nonlinear Schrödinger Equations

Burhan Gulbahar

Abstract

The scalability, error correction and practical problem solving are important challenges for quantum computing (QC) as more emphasized by quantum supremacy (QS) experiments. Quantum path computing (QPC), recently introduced for linear optic based QCs (LOQCs) as an unconventional design, targets to obtain scalability and practical problem solving. It samples the intensity from the interference of exponentially increasing number of propagation paths obtained in multi-plane diffraction (MPD) of classical particle sources. QPC exploits MPD based quantum temporal correlations of the paths and freely entangled projections at different time instants, for the first time, with the classical light source and intensity measurement while not requiring photon interactions or single photon sources and receivers. In this article, photonic QPC is defined, theoretically modeled and numerically analyzed for arbitrary Fourier optical or quadratic phase set-ups while utilizing both Gaussian and Hermite-Gaussian source laser modes. Problem solving capabilities already including partial sum of Riemann theta functions are extended. Important future applications, implementation challenges and open issues such as universal computation and quantum circuit implementations determining the scope of QC capabilities are discussed. The applications include QS experiments reaching more than $2^{100}$ Feynman paths, quantum neuron implementations and solutions of nonlinear Schrödinger equation.

The scalability of quantum resources including qubits and quantum gates, improved error correction capabilities and practical problem solving ability are the most important challenges for modern quantum computing (QC). Recent quantum supremacy experiments (QS) of Google as a success milestone for the human history of computing emphasize the importance of these properties in their set-up with $53$ qubits and $20$ cycles reaching Hilbert space size of $\approx 2^{66}$ Feynman paths making it significantly difficult to classically calculate their result [1, 2]. QS experiments target to show the computational capability of QCs such that feasible computations obtained with QCs require significant resources to perform with classical computers [3]. In the global and highly competitive race including technology giants, a wide variety of but quite complex hardware architectures are used. For example, Google, IBM and Rigetti Computing use superconducting circuits while Microsoft using topological anions generated by frozen nanowires, Ion-Q using ion traps at room temperature and D-Wave using quantum annealing technology [4]. On the other hand, both the challenges of scalability and practical problem solving capability continue to exist such that QCs that can show significant advantages in practical problem solving compared to conventional computers require much more resources such as thousands of logical qubits and hundreds of thousands physical qubits [5]. Therefore, building QC system architectures which are more tolerant to noise and decoherence combined with capabilities of error correction, practical problem solving, low hardware complexity and resource scalability is significantly important for near-term advantages of QC.

Linear optic (LO) based QCs (LOQCs) have a special place among the existing QC architectures due to easy manipulation of photons, unique features of photons not interacting with the environment in terms of decoherence, working at room temperature and maturity in classical optics for centuries [7, 8, 9]. For example, boson sampling [10] as a candidate for LOQC based QS promises experimental implementation of QS in the future while the recent experimental achievements in Ref. 11 improve the complexity of the solved problems gradually reaching Hilbert space size of $\approx 2^{48}$ . However, the existing LOQC candidates for QS have still the fundamental challenges of scalability of the resources, e.g., the requirement of single photon sources and receivers, and practical problem solving capability, e.g., matrix permanents in boson sampling. Furthermore, multi-photon entanglement resources and quantum circuits are challenging to create due to the difficulty in the interaction of photons with each other. Multi-plane diffraction (MPD) based QC system denoted by quantum path computing (QPC) as shown in Fig. 1 is recently proposed in Ref. 6 as one of the simple LOQC architectures targeting to realize scalability and practical problem solving capabilities. Sampling from the interference of exponentially increasing number of propagation paths and the freely entangled projections at different time instants [12] are utilized for QC purposes, for the first time, by exploiting coherent and classical sources such as standard laser sources. Exponentially increasing number of Feynman paths with respect to the given amount of slits and diffraction planes makes the classical calculation of the interference output significantly hard [6]. The unique form of temporal correlation freely available among the exponentially increasing number of Feynman paths in the MPD set-up is denoted as quantum path entanglement (QPE) in Ref. 12 as a novel resource to exploit for QC based on the coherence and superposition of the classical light source [13, 14]. In fact, a quantum mechanical propagator for photons with the form of the classical Fresnel diffraction integral is verified in Ref. 15 such that the classical intensity of the field is proportional to the probability density of photon detection for the position observable transversal to the propagation and in the limit for large number of quanta. Therefore, FO based set-up exploits classical light source and its intensity measurement while exploiting quantum temporal correlations in a unique MPD design. QPC promises a significant alternative to cope with the fundamental challenges of scalability of multi-photon entanglement resources and the complex requirements of single photon sources and detection mechanisms observed in conventional linear optical QC systems [9, 16, 17, 18]. It has important all-in-one advantages combining the utilization of the classical sources, i.e., either fermion or boson, hardware simplicity based on diffraction slits and detection with conventional photon counting intensity detectors without requiring simultaneous detection of multiple photons in multiple registers [10].

Besides that, another QC architecture related to multi-slit structures is denoted as duality computer (DC) which exploits duality parallelism for performing different gate operations on the sub-wave functions through sub-waves corresponding to each slit [19, 20, 21, 22] while utilized in various machine learning [23] and photonic chip applications [24]. The quantum wave divider (QWD) divides wave functions into sub-waves. In addition, sub-sub-waves can also be obtained in a multi-level QWD while quantum wave combiner (QWC) is utilized to combine the waves after performing operations on each path. It is firstly required as an open issue to explicitly model the theoretical computational complexity and QC capabilities of a multi-level version of DC with respect to the targeted QWD/QWC configurations before comparing with MPD design. For example, combined operations of uniform QWD/QWC in a DC leaving the state unchanged is not comparable with MPD modifying the state with diffractions as time evolves. Single level complexity discussions in the DC related literature include the requirement of an extra qudit [20, 21] for simulating DC device with an ordinary quantum computer which is not comparable with the tensor product structure of MPD requiring multiple qudits. The division of each wave into a sub-wave and then into a sub-sub-wave could require an exponentially increasing number of slit resources making it challenging as discussed next for Ref. 25. On the other hand, DC based architectures have a long history and maturity along two decades with many capabilities and applications including practical database search without requiring extra resources [19] verified with quantum circuit modeling in Ref. 20. DC utilizes linear combination of unitaries (LCU) for QC while becoming one of the five major techniques for designing quantum algorithms [26], in multi-party secure computation [27] and HHL quantum algorithm [28].

Furthermore, simple optical setups exploiting wave-particle duality and interferometers have the cost of exponential complexity of resources either in time, space or energy domains to achieve QC advantages as discussed in Ref. 6. For example, wave particle (WP) computer [29] exploits full optical interconnections of an $N\times N$ input signal array with an $N\times N$ output signal creating $N^{4}$ channels with a tensor product. WP computer also utilizes a filter array between the input and output to increase the number of connections in an additive manner with respect to the connections in each inter-planar region. Such architectures, including also Ref. 25 for a slit based modeling, provide advantages of parallelism compared with classical models without exploiting temporal correlations of quantum histories and their tensor product structure [12]. They utilize the tensor product only for a single inter-planar propagation, i.e., a single measurement plane directly detecting propagation from the input array. The rich set of two and three dimensional alternative optical interconnection architectures and opto-electronic computing are discussed in detail in Ref. 30 by also including multi-stage interconnection topology. Analog Fourier optics (FO) and its digital equivalent, i.e., digital FO architectures composed of smart pixel arrays of two-dimensional electronic processing units connected with optical interconnections, exploit speed and parallelism advantages of the optical design [31]. Furthermore, programmable directed logic networks are discussed in Ref. 32 by emphasizing the energy efficiency of optical architectures.

On the other hand, QPC formulation is performed for electron based set-up in Ref. 6 while theoretical studies modeling QPE in Ref. 12 and classical optical communications in Ref. 33 formulate free space propagation (FSP) of light. They do not generalize to arbitrary set-ups of FO, i.e., first order centered optical or quadratic-phase systems including arbitrary sections of free space, thin lenses, graded index media and spatial filters [34] and mathematically characterized as linear canonical transforms (LCTs)[35]. LCTs are linear integral transforms including the Fresnel and fractional Fourier transform (FRFT), scaling, chirp multiplication and some other operations as special cases while being equivalent to spatial distribution of light in phase-space optics for quadratic-phase systems [34]. Besides that, previous MPD studies utilize Gaussian sources without extending to Hermite-Gaussian (HG) beams compatible with the standard laser sources within the paraxial approximation [36]. Photonic QPC formulation is not available while important applications of QPC other than the partial sum of Riemann theta function (RTF) and period finding presented in Ref. 6 are not discussed and theoretically analyzed yet.

Diffractive and phase space optics are also getting attention in quantum technologies with periodic single plane diffraction for implementing quantum logic gates using quantum Talbot effect [62], for testing D-dimensional (qudit) Bell inequality with free space entangled quantum carpets [63] and for the evaluation of entanglement over the entire transverse field distribution of the photons [64] while without any discussion regarding the MPD based advantages. Proposed theoretical modeling and system design of photonic QPC with widely available optical components, e.g., thin lenses, free space and diffraction planes as a form of spatial filtering, provide a unique opportunity to exploit conventional FO for QC. The large amount of theoretical and experimental maturity in FO since the last century is combined with MPD based system design to realize scalable and low complexity QC systems with important capabilities and global resources for efficient implementation and development.

In this article, QPC set-up is defined, theoretically modeled and numerically analyzed for FO with arbitrary LCTs between diffraction planes. QPC system exploiting diffraction in an unconventional manner maintains photonic advantages including decoherence and noise while avoids the need to interact with multiple photons by eliminating many problems encountered in multi-photon entanglement and circuit implementations. Furthermore, the quantum nature of FO is discussed based on the experimental [37, 15] and theoretical [38] studies verifying the validity of Fresnel diffraction formulation for quantum optical propagation. Classical monochromatic light sources of both Gaussian and Hermite-Gaussian (HG) beams are utilized compatible with the standard laser sources within the paraxial approximation [36]. LCT based design which provides more flexibility is numerically compared with FSP in terms of improvement on the detection efficiency and the interference complexity defined with the magnitudes of the interfering paths and negative volume of Wigner distribution function [6, 39].

Important future applications of photonic QPC are, for the first time, introduced and theoretically modeled in an introductory and brief manner. These include the feasibility of QS experiments compared with alternative technologies, adapting certified random number generation protocols for the photonic QPC architecture [41, 42, 43, 44], quantum neural network (QNN) implementations and making the solutions of nonlinear Schrödinger equation (NLSE) easier. The detailed modeling and utilization of photonic QPC for these applications are presented as open issues.

The potential of QS experiments with photonic QPC is presented in this article to reach more than $2^{100}$ Feynman paths in a scalable set-up with several tens of diffraction planes while requiring experimental implementations for better modeling and verifying the scalability for large scale QPC set-ups. A feasible method is proposed to exponentially increase the number of Feynman paths with the cost of linearly increasing number of planes and slits allowing to obtain significantly large Hilbert space. However, it is an open issue to verify QS capability both complexity theoretically and experimentally based on the promising results in Ref. 6 and the modeling in this article such as by performing analogous modeling and experiments in Refs. 10, 11 and 40 achieved for Boson sampling. Moreover, QPC with Gaussian sources results in unique mathematical forms of wave functions on the sensor plane in (16) to be exploited for the solutions of the partial sum of RTF [45, 46, 47, 48], period finding [52] and Diophantine approximation [53] similar to the algorithms and methods in Ref. 6 but with much more design flexibility due to LCTs, diversity of the tools and maturity in the science of FO. HG sources result in different forms in (25) and (27) while closely related to the standard RTF form and requiring future studies to exploit for the solutions of numerical problems in various scientific disciplines. On the other hand, open issues and challenges for FO based QPC design are discussed to determine the scope of the proposed design for QC purposes, e.g., universal quantum computation capability, implementations of quantum circuit gates and basic search algorithms such as Grover search.

Neural networks (NNs) exploiting the quantum advantages, i.e., QNNs, improve the capabilities of classical NNs with quantum interference and superposition for deep learning applications [54] in various disciplines [55]. On the other hand, linear and unitary framework of quantum mechanics results in the challenges of implementing non-linear and dissipative dynamics of classical neural networks [55]. The state-of-the-art neuron implementations utilize various methods to introduce non-linearity including quantum measurements [55, 56]. The quantum interference among the exponentially increasing number of paths and the entanglement denoted as QPE in Ref. 12 are promising for designing and practically implementing novel design of QNNs. QPC set-up has inherently nonlinear formulation with respect to slit positions to encode the input and it operates on the quantum superposition of the inputs. Besides that, implementations of diffractive NNs utilizing single-layer [57, 58] and all-optical multi-layer diffractive architectures [59] do not exploit interference among the paths or quantum domain advantages. Photonic QPC succeeds to combine the implementations of QC and QNNs with the same hardware design of MPD as a uniquely valuable unconventional hardware architecture.

NLSE solution is very important in the analysis and performance measurement of fiber optic cables [60, 61]. It is also necessary for the solution of nonlinear Fourier transform (NLFT) which is a transformation that finds a wide range of applications with increasing importance [60]. NLSE and NLFT play a similar role for nonlinear and integration equations compared with the role of the FT in linear systems. NLFT is a transformation system for expressing the signal in the time plane by using nonlinear periodic waves or solitons [60]. It is also referred to as scattering transform. NLSE is expressed as follows [60, 61]:

[TABLE]

where $q(x,t)$ is the solution wave function that provides the periodic boundary condition ( $q(x+l,t)=q(x,t)$ and period $l>0$ ) and $\kappa$ is some variable. In this article, the speed up in NLSE solution is conjectured by exploiting RTF summations in QPC.

The remainder of the paper is organized as follows. We firstly define and theoretically model photonic QPC and its extension for FO followed by the discussion of the performance based on Wigner distribution function. Then, future applications including QS, quantum neuron implementation and solution of NLSE are introduced, theoretically modeled and the challenges are discussed. Numerical analysis for photonic MPD is provided and then open issues for realizing photonic QPC are presented.

Results

Quantum Path Computing with Optical Multi-plane Diffraction and Coherent Light Sources

MPD set-up introduced in Ref. 6 as shown in Fig. 1 is extended to optical implementations for QC by using coherent laser sources and conventional photodetectors. The set-up is composed of $N-1$ diffraction planes with $K_{j}$ slits on each plane for $j\in[1,N-1]$ and a single sensor plane indexed with $N$ while the central position of a slit is given by $X_{j,i}$ for $i\in[1,K_{j}]$ as shown in Fig. 1(a). Each slit is assumed to apply a spatial filtering of ${\color[rgb]{0,0,0}G(X_{j,i},\widetilde{\beta}_{j,i},x_{j})\,\equiv\,}\exp\big(-(x_{j}-X_{j,i})^{2}\,/\,(2\,\widetilde{\beta}_{j,i}^{2})\big{missing})$ , i.e., slit mask function, where $\widetilde{\beta}_{j,i}$ determines the slit width. The wave function on $j$ th plane is denoted with $\Psi_{j}(x_{j})$ which is the wave form after diffraction through the previous planes, i.e., with the indices $k\in[1,j-1]$ , while before diffraction through the slits on $j$ th plane. There is an exponentially increasing number of propagation paths through the slits until to the final sensor plane, i.e., $N_{p}\equiv\prod_{j=1}^{N-1}K_{j}$ , while $n$ th path includes the diffraction through a single slit on each plane with the corresponding wave function $\Psi_{j,n}(x_{j})$ on $j$ th plane. Assume that $n$ th path passes through the slit indexed with $s_{n,j}$ on $j$ th plane and we define the path vectors $\vec{x}_{N-1,n}^{T}\equiv\begin{bmatrix}X_{1,s_{n,1}}\,\,X_{2,s_{n,2}}\,\,\ldots\,\,X_{N-1,s_{n,N-1}}\end{bmatrix}$ and $\vec{x}_{N-1,\vec{s}}^{T}\equiv\begin{bmatrix}X_{1,s_{1}}\,\,X_{2,s_{2}}\,\,\ldots\,\,X_{N-1,s_{N-1}}\end{bmatrix}$ where $\vec{s}\equiv\begin{bmatrix}s_{1}\,\,s_{2}\,\,\ldots\,\,s_{N-1}\end{bmatrix}$ and $(.)^{T}$ is the transpose operation. The mapping between the path index $n$ and slit index $s_{n,j}$ for the path is defined with the function $n=f_{s2n}(\vec{s})$ where $s_{n,j}$ is predefined for each $n$ . Furthermore, $\vec{0}_{k}$ is the column vector of length $k$ with all zeros and $\mathbf{0}_{k}$ is the square matrix of all zeros with the size $k\times k$ . Similarly, rectangular matrices are shown with $\mathbf{0}_{k,l}$ . In the rest of the article, a parameter $B$ depending only on $\vec{\beta}_{N-1,n}\equiv[\widetilde{\beta}_{1,s_{n,1}}\,\,...\,\,\widetilde{\beta}_{N-1,s_{n,N-1}}]$ but not on $\vec{x}_{N-1,n}$ is denoted with $\widetilde{B}_{j,n}$ on each $j$ th plane including $\widetilde{.}$ over the symbol. Therefore, if the slits are chosen with the same $\widetilde{\beta}_{j,{s_{n,j}}}=\beta_{j}$ specific to each plane, then $\widetilde{B}_{j,n}$ becomes independent of $n$ and is converted to the notation $B_{j}$ .

MPD set-up shown in Fig. 1 is utilized for QC denoted by QPC in Ref. 6 by sampling the interference of exponentially increasing number of interfering paths. The capability to theoretically characterize QPC with quantum FO provides future applications for both QC and quantum information theory by exploiting energy efficient combination of optical elements. Intensity sampling on the sensor plane ( $I[k]$ ) for the MPD set-up in Refs. 6, 12 generates a black-box (BB) function $f_{BB}[k]$ with the following special form to utilize in solutions of important and classically hard number theoretical problems:

[TABLE]

where $k\in\mathbb{Z}$ , $T_{s}\in\mathbb{R^{+}}$ is a sampling interval, $\widetilde{A}_{\vec{s}}\in\mathbb{R^{-}}$ , $\widetilde{B}_{\vec{s}}\in\mathbb{R^{+}}$ and $\widetilde{\Upsilon}_{\vec{s}}\,\in\mathbb{C}$ . The complex valued matrix $\mathbf{\widetilde{H}}_{\vec{s}}\equiv\mathbf{\widetilde{H}}_{R,\vec{s}}\,+\,\imath\,\mathbf{\widetilde{H}}_{I,\vec{s}}$ and the vector $\vec{\widetilde{h}}_{\vec{s}}\equiv\vec{\widetilde{c}}_{\vec{s}}+\,\imath\,\vec{\widetilde{d}}_{\vec{s}}$ have the values depending on $\widetilde{\beta}_{j,i}$ for $j\in[1,N-1]$ and $i\in[1,K_{j}]$ corresponding to the specific selection of slits in the path $\vec{s}$ , inter-plane durations for the particle propagation, particle mass $m$ (for electron based set-ups in Refs. 6, 12), beam width $\sigma_{0}$ of the Gaussian source wave packet and Planck’s constant $\hbar$ . The calculation of (2) in an efficient manner is significantly hard while two different methods utilizing (2) for practical problems are introduced. Solution for the partial sum of RTF or multi-dimensional theta function is the first application with importance in number theory and geometry [45, 46, 47, 48]. The second method utilizes MPD with the phase term $\vec{\widetilde{d}}_{\vec{s}}^{T}\,\vec{x}_{N-1,\vec{s}}$ in $\exp\big(\vec{\widetilde{h}}_{\vec{s}}^{T}\,\vec{x}_{N-1,\vec{s}}\,k\,T_{s}\big{missing})$ for period finding [52] and the solution of specific instances of SDA problems [53].

The basic unit of QC systems, i.e., the qubit, is defined on a two-state system where discretized degrees of freedom (DoF) of photons including path, transverse-spatial modes and time/frequency bins are exploited to create high-dimensional entanglement [49]. For example, multi-slit structures are already utilized to define spatial qudits by projecting the wave function into the transverse position and momentum Hilbert spaces through slits and characterizing their properties using their propagation in free space [50, 51]. The qubit states in Ref. 51 are expressed in the basis $\ket{l},\ket{r}$ representing the photons passing through either the left or the right slit. However, entangled multiple photons, e.g., photon pair $A$ and $B$ , are conventionally generated through spontaneous parametric down-conversion (SPDC) to realize multi-photon entangled state, e.g., $\big{(}1\,/\,\sqrt{2}\big{)}\,\big{(}\ket{l_{A}}\ket{r_{B}}\,+\,\ket{r_{A}}\ket{l_{B}}\big{)}$ . The fundamental difference of MPD based qudits from multi-photon slit based entangled spatial qudits is the utilization of the tensor product structure for each single photon in time domain rather than spatially among multiple photons obtained through SPDC [6]. The projection events through the slits of consecutive planes are freely entangled at two different time instants denoted as QPE with the detailed modeling presented in Ref. 12 based on consistent histories and entangled histories frameworks. The presented free entanglement in time domain provides an important advantage exploiting directly the classical light sources and not requiring the difficult coupling of multiple photons. The concept of free entanglement is introduced for boson sampling exploiting boson statistics of a number of indistinguishable bosons while they still require multiple photons, and generation and detection mechanisms for single photons [10].

History state in MPD is composed of diffraction events as follows [12]:

[TABLE]

where $\mathbf{P}_{j,s_{n,j}}$ represents the projection operator through the slit indexed with $s_{n,j}$ , $\pi_{n}$ as [math] or $1$ determines a compound set of trajectories, $\odot$ denotes tensor product operation and $\left[\rho_{0}\right]$ denotes the initial state. The analogy of MPD based multiple qubits with the general two qubit state of two photons is represented as shown in Fig. 2 in the basis of $\ket{U}$ and $\ket{L}$ for the upper and lower slits, respectively. The general state for the projection through two diffraction planes indexed with $A$ and $B$ is represented as follows:

[TABLE]

where the amplitudes are denoted by $a_{ij}$ , and $i$ and $j$ denote the projection through upper or lower slits. The $A$ and $B$ in the MPD set-up denote the indices of planes for the projection of a single photon at different time instants rather than the indices of two photons as in the entangled state of two spatial qubits of two photons. There are four different projection history states where the wave function whose intensity to be measured on the final detection plane, i.e., $\Psi_{{\color[rgb]{0,0,0}3}}(x)$ , is described as the interference of four different wave function histories corresponding to each trajectory, i.e., $\Psi_{{\color[rgb]{0,0,0}3},j}(x)$ for $j\in[1,4]$ :

[TABLE]

$\Psi_{{\color[rgb]{0,0,0}3,}1}(x)$ corresponds to $a_{uu}\ket{U_{A}U_{B}}$ and the other components are defined as shown in Fig. 2. Each component depends in a complex manner on the slit geometries as modeled by the RTF. QPC applications of MPD based high dimensional entangled states do not include any measurement by closing or opening slits but a final measurement on the detector plane obtaining the complicated interference pattern of exponentially many Feynman paths [6].

QPC based on FO promises expanding the set of solvable problems both with LCT based general system design and also the sources including HG beams. Furthermore, a discussion is included to utilize non-Gaussian slits with the proposed mathematical modeling in the Open Issues and Discussion section. Propagation through Fourier optical systems based on Fresnel diffraction is modeled emphasizing the quantum nature of Fresnel diffraction and FO in the Methods section. Next, MPD modeling is proposed for Fresnel diffraction and arbitrary LCT based optical systems by utilizing the proposed kernels.

Quantum Path Computing with Fourier Optical Systems

The set-up in Fig. 1 defined with FSP and electron based theoretical formulation is extended to optical systems of LCT as shown in Fig. 3. The kernel of one dimensional (1D) quadratic-phase system or LCT converting the input signal $f(x_{0})$ to the output $\int_{-\infty}^{\infty}dx_{0}\,K_{LCT}^{(a,b,c,d)}(x_{1},x_{0})\,f(x_{0})$ is represented as follows:

[TABLE]

In matrix notation, it is shown with the following unit-determinant matrix:

[TABLE]

where $a\,d-b\,c=1$ and the matrix for the composition of two consecutive systems represented by $\mathbf{M}_{1}$ and $\mathbf{M}_{2}$ is calculated by the multiplication $\mathbf{M}_{2}\,\mathbf{M}_{1}$ . The kernel matrices $K_{FS}(x_{1},x_{0})$ and $K_{HO}(x_{1},x_{0})$ denoting FSP kernel in phase-space optics [34] and the kernel based on quantum mechanical harmonic oscillator (HO) modeling of the evolution of light wave function [15] (in analogy with FRFT implementation), respectively, are defined in the Methods section while discussing quantum FO. Some simple examples of LCT matrices for propagation of length $L_{j,j+1}$ between $j$ th and $(j+1)$ th planes are given as follows [34]:

[TABLE]

where $\mathbf{M}_{S}$ scales with $\Psi_{j+1}(x_{j+1})\equiv(1\,/\,\sqrt{a_{j,j+1}})\,\Psi_{j}(x_{j+1}\,/\,a_{j,j+1})$ , the kernel for $\mathbf{M}_{L}$ is $e^{-\imath\,\pi\,x_{1}^{2}\,/\,(\lambda\,f)}$ and $m_{\lambda}\equiv 2\,\pi\,\hbar\,/\,(\lambda\,c)$ is defined in the Methods section after discussing (47).

The varying forms of wave functions on the measurement plane extending (2) are modeled which are promising to be utilized in QC applications. It is presented next such that obtained forms are similar to (2) while having higher flexibility of system design. The theoretical modeling of BB functions for quantum HO based or FRFT based light propagation modeling with Gaussian sources is presented next with the wave function in (16). Gaussian source case is also extended to arbitrary LCTs. Similarly, the wave functions for arbitrary LCTs with HG sources are presented in (25) and (27) next.

An arbitrary LCT with the matrix parameters $\{a_{j,j+1},b_{j,j+1},c_{j,j+1},d_{j,j+1}\}$ is implemented in phase-space optics by consecutive applications of FSP of length $L_{a,j,j+1}$ , then thin lens of focal length $f_{j,j+1}$ , and another FSP of length $L_{b,j,j+1}$ [35]. LCT matrix $\mathbf{M}_{LCT}$ is calculated as follows:

[TABLE]

where $\tau^{\star}_{a,j}\equiv L_{a,j,j+1}\,/\,c$ and $\tau^{\star}_{b,j}\equiv L_{b,j,j+1}\,/\,c$ , and the middle matrix is for the effect of thin lens [34]. FRFT with scaling is a special case of LCT as discussed in the Methods section. Therefore, FSP, FRFT and arbitrary LCT based QPC set-ups are implemented with the universal configuration in Fig. 3.

QPC with Fresnel Diffraction and FRFT by using Gaussian Sources

Firstly, two special cases of LCTs are considered, i.e., FSP of light and propagation modeled with FRFTs denoting graded-index media propagation as the solution of the quantum HO in Ref. 15. Furthermore, we assume that the source distribution has a Gaussian form of $\Psi_{0}(x_{0})=\mbox{exp}\big{(}-\,x_{0}^{2}\,/\,(2\,\sigma_{0}^{2})\big{)}\,/\,\sqrt{\sigma_{0}\,\sqrt{\pi}}$ while HG waveforms as eigenfunctions of FRFTs [34] are considered for the general case of LCTs in the next section. It is assumed that the optical system between the planes results in the kernels $K_{FS}(x_{1},x_{0})$ and $K_{HO}(x_{1},x_{0})$ defined in (45) and (47) based on Fresnel diffraction integral for free space and quantum HO solution [15], respectively. The definition and the derivation of HO based kernel with the following kernel matrix for the propagation duration of $t_{01}$ are detailed in the Methods section while we are discussing the quantum mechanical modeling of FO, i.e., denoting with quantum FO:

[TABLE]

The important observation is that iterative integration with $K_{HO}(x_{1},x_{0})$ results in the final intensity distribution of MPD with the same form of $K_{FS}(x_{1},x_{0})$ while with different algorithms for calculating the iteration parameters as shown in Table 1 in the Methods section. The kernel $K_{FS}(x_{1},x_{0})$ has the same form with $K_{m,FS}(x_{1},x_{0})$ used for QPC modeling in Ref. 6 by replacing the electron mass $m$ with the photon equivalent mass $m_{\lambda}$ . Therefore, the same formulations are utilized for the cases of FS and HO solutions while modeling the sampled wave function on the sensor plane with iterations and the resulting structure of problem solving capabilities.

The wave function for $n$ th path on the sensor plane for the general case of non-uniform slit widths is given by the following by using the iterative formulation:

[TABLE]

It is further simplified by extraction of $\vec{x}_{N-1,n}$ dependent parts and summing the contributions from each path as follows:

[TABLE]

where $\widetilde{\Upsilon}_{N,n}^{G}\equiv\chi_{0}\,\big{(}\prod_{j=1}^{N-1}\sqrt{\widetilde{\xi}_{j,n}}\big{)}$ , and the complex vector $\vec{\widetilde{h}}_{N-1,n}$ and the matrix $\mathbf{\widetilde{H}}_{N-1,n}^{HO/G}$ are defined in the Methods section for the HO case with simplified formulation compared with the case for electron based FSP set-up in Ref. 6. The corresponding iteration parameters are given in Table 1 in the Methods section.

We have not included the effects of special forms of $K_{HO}(x_{j+1},x_{j})$ with $\omega\,t_{j,j+1}\,=\,k\,\pi$ for $k\in\mathbb{Z}$ corresponding to integer multiples of FRFT order $2$ since the result is $\Psi_{j+1}(x_{j+1})\equiv\Psi_{j}(\pm x_{j+1})$ (inserting $\pm x_{j+1}$ into $\Psi_{j}(x_{j})$ ) [34]. This case can be simply realized by assuming that spatial filtering operations of the slits on $j$ th and $(j+1)$ th planes are combined on a single plane by also noting that whether the wave function is reversed or not. For example, multiple inter-plane propagation intervals can result in multiple reversals with the overall effect of the identity and combined spatial filtering of Gaussian slits.

QPC with Arbitrary Linear Canonical Transforming Optical Systems

Gaussian Sources: The resulting final intensity of MPD propagation for the case of $K_{LCT}^{(a,b,c,d)}$ ( $b_{j,j+1}\neq 0$ ) with Gaussian sources has the same form with $K_{HO}$ while with different algorithms for calculating iteration parameters in Table 1 in the Methods section and replacing $\widetilde{\mathbf{H}}_{N-1,n}^{HO/G}$ with $\widetilde{\mathbf{H}}_{N-1,n}^{LCT/G}$ . Therefore, all the derivations utilized for $K_{HO}$ including the explicit forms of the wave function are applicable. We have not included $K_{LCT}^{(a_{j,j+1},\,b_{j,j+1},\,c_{j,j+1},\,d_{j,j+1})}$ with $b_{j,j+1}=0$ for simplicity. Two simple cases are scaling and chirp multiplication with $a_{j,j+1}=d_{j,j+1}=1$ resulting in $\Psi_{j+1}(x_{j+1})\,\equiv\,\exp(\imath\,\pi\,c\,x_{j+1}^{2})\,\Psi_{j}(x_{j+1})$ [34]. These cases further improve the flexibility of the LCT system to realize the desired transformation on the wave function.

A simple example is presented with explicit expressions as follows for $K_{1}=3$ , $K_{2}=2$ , $N=3$ , $b_{j,j+1}\neq 0$ for $j\in[0,2]$ , fixed slit width parameters $\beta_{1}$ and $\beta_{2}$ for simplicity and $\vec{x}_{2,n}^{T}\equiv\begin{bmatrix}X_{1,s_{n,1}}\,\,X_{2,s_{n,2}}\end{bmatrix}$ :

[TABLE]

where the following parameters are derived by using Table 1 in the Methods section:

[TABLE]

where it is observed that the parameters are some rational complex polynomials of $\beta_{1}$ , $\beta_{2}$ and $\sigma_{0}$ with the coefficients in terms of polynomial expressions of $a_{j,j+1}$ , $b_{j,j+1}$ and $d_{j,j+1}$ for $j\in[0,2]$ with varying orders of ( $\beta_{1}$ , $\beta_{2}$ , $\sigma_{0}$ ): reaching ([math], $2$ , $2$ ) for $pol_{1}$ ; ( $2$ , $2$ , $2$ ) for $pol_{2}$ ; ([math], [math], $2$ ) for $pol_{3}$ ; ( $2$ , [math], $2$ ) for $pol_{4}$ ; ( $6$ , $4$ , $8$ ) for $pol_{5}$ , $pol_{7}$ , $pol_{8}$ and $pol_{9}$ ; ( $8$ , $4$ , $8$ ) for $pol_{6}$ , ( $4$ , $2$ , $5$ ) for $pol_{10}$ ; ( $4$ , $2$ , $6$ ) for $pol_{11}$ ; ( $4$ , $4$ , $4$ ) for $pol_{12}$ , $pol_{13}$ and $pol_{14}$ .

The highly complicated expressions for the polynomials are explicitly shown in Table 2 in the Methods section such that they are obtained by using the iteration method in Table 1 in the same section. It is possible by using the explicit expressions directly to perform various gedanken experiments and computational complexity analysis with any number of slits and LCTs.

A simple numerical example is presented as shown in Fig. 4. The scaling property of thin lens is utilized in the Numerical Results section to improve the intensity of the diffraction on the final detection plane. For example, a simple Gaussian source beamwidth of $\sigma_{0}=20$ ( $\mu$ m) and $\lambda=650$ (nm) shown in Fig. 4(b) is scaled by shifting the position of the lens of focal length $60$ (mm) inside the interval of fixed length of $L_{01}=L_{a,01}+L_{b,01}=t_{01}\times c$ with $t_{01}=1$ (ns) as shown in Fig. 4(a). The shift is modeled with the ratio $r_{L}=L_{a,01}\,/\,L_{01}$ . It is observed in Fig. 4(c) that the intensity of the wave function can be focused with respect to the positions of the slits on the first plane.

Hermite-Gaussian Sources: If the source is chosen as the standard HG waveform of $\Psi_{0}(x_{0})\,=\,(2^{1/4}\,/\,\sqrt{W_{0}\,2^{l}\,l!})\,\mbox{exp}\big{(}-\,\pi\,x_{0}^{2}\,/\,W_{0}^{2}\big{)}\,H_{l}(\sqrt{2\,\pi}\,x_{0}\,/\,W_{0})$ for $K_{LCT}^{(a_{j,j+1},\,b_{j,j+1},\,c_{j,j+1},\,d_{j,j+1})}$ ( $b_{j,j+1}\neq 0$ ) where $H_{l}(x)\,\equiv\,(-1)^{l}\,e^{x^{2}}\,d^{l}e^{-x^{2}}\,/\,dx^{l}$ is the $l$ th order Hermite polynomial [34, 65], then $\Psi_{N,n}(x_{N})$ for $n$ th path is obtained as follows by using the integral equality of HG functions in the Methods section:

[TABLE]

where the parameters $\chi_{01}$ , $\chi_{j,j+1,n}$ , $\widetilde{u}_{j,j+1,n}$ , $v_{j,j+1,n}$ , $\widetilde{g}_{j,j+1,n}$ and $h_{j,j+1,n}$ obtained in an iterative manner for $j\in[1,N-1]$ are calculated with simple algebra for $n$ th path and shown in Table Matrix formulation for LCT system with Hermite-Gaussian sources in the Methods section. Simple algebraic manipulations of (22) to extract $\vec{x}_{N-1,n}$ dependent parts result in the following simplification:

[TABLE]

where $\widetilde{\chi}_{a,j,j+1,n}$ for $j\in[1,N-1]$ is defined in Table Matrix formulation for LCT system with Hermite-Gaussian sources, the vectors $\vec{\widetilde{\gamma}}_{N-1,n}$ and $\vec{\widetilde{\eta}}_{N-1,n}$ , and the matrix $\widetilde{\mathbf{H}}_{N-1,n}^{LCT/HG}$ are defined in the Methods section. It is observed in (23) that each different path results in a different shift on Hermite polynomial determined with $\vec{\widetilde{\eta}}^{T}_{N-1,n}\,\vec{x}_{N-1,n}$ even for the uniform $\beta_{j}$ for each slit on $j$ th plane. As a result, the final wave function on the sensor plane denoted with $\Psi_{N}^{HG}(x_{N})$ for the general case of non-uniform slit widths defined with $\widetilde{\beta}_{j,n}$ for $j\in[1,N-1]$ and $n\in[0,N_{p}-1]$ is given by the following:

[TABLE]

where $\widetilde{\Upsilon}_{N,n}^{HG}\,\equiv\,\chi_{01}\,\big{(}\prod_{j=1}^{N-1}\widetilde{\chi}_{a,j,j+1,n}\big{)}$ $e^{-\,\imath\,\pi\,(N-2)\,/\,4}$ and with the similarity to the form in (16) for the HO solution except multiplicative Hermite polynomial for each $n$ th path. The complexity of calculating the Gaussian form in (16) is classically hard as thoroughly discussed in Ref. 6 which requires to compute a special form of partial sum of RTF while the complex vector $\vec{\widetilde{h}}_{N-1,n}$ and the matrix $\mathbf{\widetilde{H}}_{N-1,n}^{HO/G}$ varying for each path making it much harder compared with the computation of conventional partial sum of RTF. Therefore, the complexity characterization of computing $\Psi_{N}^{HG}(x_{N})$ is an open issue while it is expected to be significantly hard since each summation term depends on path index $n$ with varying vector and matrix parameters while also including a product term of Hermite polynomial for each path making it harder.

If the uniform slit width case is chosen and the path independent variables are denoted with $\widetilde{\Upsilon}_{N,n}^{HG}=\Upsilon_{N}^{HG}$ , $\widetilde{\chi}_{a,j,j+1,n}=\chi_{a,j,j+1}$ , $\widetilde{\mathbf{H}}_{N-1,n}^{LCT/HG}=\mathbf{H}_{N-1}^{LCT/HG}$ , $\vec{\widetilde{\gamma}}_{N-1,n}=\vec{\gamma}_{N-1}$ , $\widetilde{u}_{N-1,N,n}=u_{N-1,N}$ , $\widetilde{g}_{N-1,N,n}=g_{N-1,N}$ , $\vec{\widetilde{\eta}}_{N-1,n}=\vec{\eta}_{N-1}$ , then (25) is transformed into the following:

[TABLE]

It is further simplified as follows by using the useful identity $H_{l}(x\,+\,y)=(H+2\,y)^{l}$ in Ref. 66 where $H^{k}\equiv H_{k}(x)$ :

[TABLE]

where $\big{(}H^{\star}(x_{N})\big{)}^{k}\,\equiv\,H_{k}(g_{N-1,N}\,x_{N})$ . The computational complexity of calculating $\Psi_{N}^{HG,U}(x_{N})$ is similarly expected to be significantly hard since the mathematical form is more complicated compared with the partial sum of RTF.

The set-up parameters including the slits, lenses and inter-plane distances are required to be tuned in order to obtain the desired vectors $\vec{\widetilde{h}}_{N-1,n}$ , $\vec{\widetilde{\gamma}}_{N-1,n}$ , $\vec{\widetilde{\eta}}_{N-1,n}$ and matrices $\widetilde{\mathbf{H}}_{N-1,n}^{HO/G}$ , $\widetilde{\mathbf{H}}_{N-1,n}^{LCT/G}$ and $\widetilde{\mathbf{H}}_{N-1,n}^{LCT/HG}$ in (16), (25) and (27) for the targeted number theoretical problems. Next, Wigner distribution is defined where its negative volume is regarded as an indicator of non-classicality.

Wigner Distribution, Negativity and Path Magnitudes

The momentum domain wave function $\Psi_{p,j}(p_{j})$ is defined as Fourier transform of spatial representation of wave function $\Psi_{j}(x_{j})$ on $j$ th plane as follows:

[TABLE]

The distribution of energy through space-momentum phase-space is described by Wigner distribution function defined as follows [39, 6]:

[TABLE]

The negative volume of Wigner function defined in Ref. 39 and utilized in Ref. 6 to describe the increasing non-classicality or time-domain entanglement resources in Ref. 12 is described as $V_{j}\equiv$ $\big{(}\int\int|W_{j}(x_{j},p_{j})|\,dx_{j}\,dp_{j}\,-\,1\big{)}\,/\,2$ . On the other hand, the probability of the particle to be detected on $j$ th plane, i.e., to be diffracted through $(j-1)$ th plane, is computed as follows:

[TABLE]

In this article, the increasing interference complexity and non-classical nature of MPD based time-domain entanglement resources are assumed to be characterized by utilizing $V_{j}$ and by observing the magnitudes of the interfering paths defined as $P_{E,n}(j)\equiv\int dx_{j}\,|\Psi_{j,n}(x_{j})|^{2}$ for each $n$ th path. Therefore, more paths with large magnitudes and $V_{j}$ emphasize increasing interference and non-classicality. Characterizing the correlation between the distribution of path magnitudes and $V_{j}$ is an open issue since the behavior is highly set-up specific as observed in the Numerical Results section. On the other hand, path magnitudes throughout the whole plane may not reflect their localized characteristics such as effecting some sample locations more compared with the others. Therefore, it is an open issue to characterize the interference complexity in terms of the intensity distribution of the paths while the path magnitudes are taken as a reference for simplicity in this article. Next, the potential future applications of QPC architecture based on FO are presented.

Future Applications

Applications for Quantum Supremacy and Certified Random Number Generation

The scalable structure of QPC with coherent light sources and low complexity FO promise the large scale implementation of QS experiments. Recently, some experiments are performed with 53 qubits with superconductor based architectures of Google as a milestone for the human history of computational capabilities [1, 2]. In this article, we propose that a similar experiment could be formulated in the QPC set-up as the problem of finding the distribution of light intensity on the photodetector array plane by randomly generated slit positions and widths analogical to random circuit sampling [43, 67]. Fig. 5(a) shows the system architecture that is similar to the randomly generated quantum gates. The aim is to perform a complexity analysis of a randomly generated QPC architecture. The total number of Feynman paths is expressed as follows for $L$ diffraction planes:

[TABLE]

Suppose that energy decreases by $1/s$ ( $s>1$ ) to a total of $1/s^{L}$ and the number of significant paths decreases by $1/r$ ( $r>1$ ) leading to a decrease in the total number of effective paths by $1/r^{L}$ approximating the final intensity distribution. Therefore, if we define the number of paths that can be realized for unit source energy with $\widetilde{N}_{path}$ , then the following is obtained:

[TABLE]

If we measure the Hilbert space size created by the total number of paths by defining the virtual qubit number and assuming $s\equiv 2^{s^{*}}$ , $r\equiv 2^{r^{*}}$ , $m\equiv 2^{m^{*}}{\color[rgb]{0,0,0}2^{2}}$ (assuming the minimum of $m=8$ with $m^{*}=1$ ), $k\equiv 2^{k^{*}}$ , then the following definition is obtained:

[TABLE]

The basic expression determining the size of the Hilbert space is denoted with the gain $G=m^{*}+k^{*}-s^{*}-r^{*}$ . It provides the cumulative effect of increasing number of paths due to the linearly increasing number of the slits with the coefficient $m$ and the initial number of the slits $k$ combined with the decreasing number of paths due to the inter-plane attenuation coefficients $s$ and $r$ for the effective number of the paths. As shown in Figure 5(b) for $L=N-1$ , even at very low gain rates, e.g., $G=-5$ with $m^{*}=1$ , that is, where the spreading energy drops very quickly and the number of significant paths is too low, the number of virtual qubits reaches hundreds. Furthermore, $m^{*}+k^{*}$ , which can be designed flexibly in a multi-slit architecture, is adapted against the low gain. In addition, even with $N=10$ planes, Hilbert space size of approximately reaching hundreds of virtual qubits is obtained for the case of high gain $G$ . For example, assume the worst case situation such that diffracted photon forms a large amplitude path by diffracting through a locally limited number of slits on the next plane denoted by the parameter $\widetilde{r}$ . In other words, the slit locations distant apart on consecutive planes will not form a large amplitude path and the number of effective paths increases as the multiplication by $\widetilde{r}^{j}$ after diffracting through $j$ consecutive planes (by assuming there is enough number of closely spaced slits on the next plane). This can be adjusted by increasing the inter-plane distance such that each diffracted beam will expand to a larger area on the consecutive plane. Moreover, assume that the slits are placed close enough to keep the probability of the diffraction through the next plane roughly constant, i.e., $1\,/\,s$ , as observed in MPD simulation studies in Ref. 6. Then, if the condition $s<\widetilde{r}$ is satisfied, the number of effective paths will increase with the multiplication by $(\widetilde{r}\,/\,s)^{j}$ through $j$ consecutive planes. For example, effective number of paths in Ref. 6 is observed to be increasing with $(\widetilde{r}\,/\,s)^{3}>70$ even by removing many effective paths (Fig. 7(b) in Ref. 6). In other words, assuming $\widetilde{r}\,/\,s\approx 2^{2}$ allows to reach $2^{100}$ effective Feynman paths with $N\approx 51$ planes. Therefore, without requiring extensive simulations, it is clear that the number of paths increases exponentially with specially adjusted set-up parameters of inter-plane distance, slit widths and distributions.

Each path will form a unique contribution to the overall intensity. There is no apparent way of calculating the exact final intensity other than identifying and summing the contribution of each path. Then, increasing the complexity of MPD set-up makes it harder to calculate the contribution of each path until reaching to the QS scale. It is possible to compare roughly with Google QS experiment where the computational complexity requires the calculation of $4^{31}\times 2^{4}=2^{66}$ different Feynman paths for $53$ qubits and $20$ cycles (Table XI in Ref. 1) as shown in Fig. 5(b). Compared to the $q_{path}=66$ , the proposed QPC system architecture suggests to perform future QS experiments with a simple system structure.

The open issues include the rigorous characterization of the computational complexity of sampling from MPD exploring the relations among the paths in terms of magnitude and distribution. Furthermore, the inter-plane gain $G$ is required to be both theoretically modeled and experimentally measured for random and large diffraction architectures. Another open issue is to analyze the modeling of the sampling problem of QPC with universal quantum circuits and to determine the computational complexity class, e.g., the relations with $BQP$ [10] and complexity theoretical fundamentals of QS experiments [40]. Experimental implementation requires slit design and manufacturing, sensitive photon detection due to the attenuation after large number of planes and spatially coherent light sources covering all the paths reaching to the detector plane. The number of slits, i.e., determined by the parameters $k$ , $m$ and $N$ , is limited by the capability to realize significant number of small width slits, e.g., in micrometer scale, on an appropriate planar surface such as by patterning metallic slits on glass substrate [68]. The beam width and inter-planar distances should be adapted for spatial coherence of the light diffracting through planes [12]. However, the linear modeling of the architecture in Fig. 5(a) and the expansion of light beam through propagation allow to realize a feasible architecture in future experimental implementations.

On the other hand, achieving QS experiments allows to adapt certified random number generation protocols for the QPC architecture [41, 42, 43, 44]. Although there are recent high speed, e.g., on the orders of several Gbit/s, random number generation protocols working in a local manner and exploiting the sampling of interference based intensity fluctuations of laser pulses such as Refs. 69 and 70, the idea of randomness extraction from QS experiments proposed by Aaronson [41] in a way allowing to download from a remote and trusted public source is new. The user interacts with a remote QC and makes it to generate random bits without any trust to the QC itself. Similar to Aaronson’s protocol, it is possible to firstly collect random numbers from a trusted computer. Then, using these numbers, the widths and planar distributions of the slits are determined to have a random diffraction set-up by assuming that the mechanical structure of the device can be modified remotely. Then, intensity distribution in the photodetector array is measured. Therefore, both the random structure of MPD set-up and interference of the exponential number of paths result in a very difficult measurement output intensity to efficiently calculate with classical computers. However, it is an open issue whether it is possible to utilize the proposed remote QC device based on QPC similar to the Aaronson’s protocol which realizes sampling from the n-qubit output of the quantum circuit and performs Heavy Output Generation (HOG) tests. On the other hand, QPC does not allow to sample the probability of a single path but the interference of exponentially many number of paths. Therefore, it is a challenge and open issue to adapt the interference sampling operation in QPC for a similar complexity theoretical proof of randomness generation.

Neuromorphic applications with quantum neuron implementations

In a classical artificial neuron implementation, the output is represented as $O=f(\sum_{i=1}w_{i}\,x_{i}+b)$ where $f(.)$ is the nonlinear activation function, $w_{i}$ is the $i$ th synaptic weight, $x_{i}$ is the $i$ th input and $b$ is the bias as shown in Fig. 6(a). QPC based neuron has inherent nonlinearity with the form of $f(\vec{x}^{T}\,\mathbf{H}\,\vec{x}\,+\,\vec{x}^{T}\,\vec{h}\,+\,b)$ between the input slit position vector $\vec{x}$ and output (O) as shown in Fig. 6(b) based on the MPD formulation for Gaussian sources in (16). The quadratic weighting $\vec{x}^{T}\,\mathbf{H}\,\vec{x}$ , linear weighting $\vec{x}^{T}\,\vec{h}$ and bias $b$ are fed into the nonlinear function $f(.)\equiv exp(.)$ for the LCT set-up where $\mathbf{H}\equiv\mathbf{\widetilde{H}}_{N-1,n}^{LCT/G}$ , $\vec{h}\equiv\vec{\widetilde{h}}_{N-1,n}\,x_{N}$ , the bias $b\equiv\log(\widetilde{\Upsilon}_{N,n}^{G})\,+\,(\widetilde{A}_{N-1,n}\,+\,\imath\,\widetilde{B}_{N-1,n})\,x_{N}^{2}$ , $\vec{x}\equiv\vec{x}_{N-1,n}$ and $x\equiv x_{N}$ is the measurement position.

The example output through the slits $Y_{j}$ for $j\in[1,2]$ in Fig. 6(b) depends on interfering quantum superposition of input combinations as follows while assuming path independent forms of the variables for simplicity, i.e., $\mathbf{H}\equiv\mathbf{H}_{N-1}^{LCT/G}$ , $\vec{h}\equiv\vec{h}_{N-1}\,x_{N}$ and $b\equiv\log(\Upsilon_{N}^{G})\,+\,(A_{N-1}\,+\,\imath\,B_{N-1})\,x_{N}^{2}$ :

[TABLE]

where $G(y,\beta_{N},x_{N})$ is the slit mask function depending on the output slit ${\color[rgb]{0,0,0}y\,=\,}Y_{i}$ for $i\in[1,2]$ of the $N$ th plane and $\beta_{N}$ is the fixed slit width with path-independent assumption for simplicity. The parameters are possible to depend on each path with variable slit masks. Exponentially large number of synaptic chains (paths) through slit inputs, their quantum interference and simplicity to sample the intensity output $|O|^{2}$ provide significant opportunities to exploit for quantum advantages.

The challenges include designing the quantum neuron based on the slit positions as inputs while changing the weight in a controllable manner. Besides that, extensive simulation studies are required to practically observe the quantum advantages for various problems. The positions of the slits are required to be modified dynamically with special designs. In addition, large scale QNN implementations both in simulations and experiments are required to observe the performances in various problems, e.g., pattern recognition or machine learning for very large problem sizes.

Solution of nonlinear Schrödinger equation

Finite-band solutions of NLSE in (1) are expressed with RTF as follows [60, 61]:

[TABLE]

where Riemann spectrum is $(\mathbf{Y},\vec{k},\vec{\omega},\vec{\delta}^{-},\vec{\delta}^{+})$ , Riemann period matrix $\mathbf{Y}$ [71] is calculated using $\vec{k}$ and $\vec{\omega}$ together with nonlinear spectrum data (Appendix to Section 24 in Ref. 48) and the partial sum of RTF denoted as $\Theta_{M}$ converging to $\Theta$ for $M\rightarrow\infty$ is defined as follows:

[TABLE]

where $\mathbf{{\color[rgb]{0,0,0}\Gamma}}$ is a complex matrix, $\vec{y}$ is a complex vector and $\vec{a}^{T}\equiv[a_{1}\,\,a_{2}\,\,\ldots a_{N-1}]$ . As shown in Ref. 6, if $j\in[1,N-1]$ and $a_{j}\in S_{M}\equiv[-M,M]$ , we select $X_{j,i}\in S_{M}\Delta x_{j}$ and also if the slit widths are kept constant for each plane, then $A_{\vec{s}}$ , $B_{\vec{s}}$ , $\Upsilon_{\vec{s}}$ , $\mathbf{H_{\vec{s}}}$ , $\vec{c}_{\vec{s}}$ and $\vec{d}_{\vec{s}}$ values become path independent, i.e., $A_{N-1}$ , $B_{N-1}$ , $\Upsilon_{N}$ , $\mathbf{H}_{N-1}$ , $\vec{c}_{N-1}$ and $\vec{d}_{N-1}$ while the superscript $(.)^{LCT/G}$ removed for simplicity, and (2) is transformed as follows for $x\equiv k\,T_{s}$ :

[TABLE]

where ${\color[rgb]{0,0,0}\mathbf{\widehat{\Gamma}}}\equiv-\,\mathbf{D}\,\mathbf{H}_{\color[rgb]{0,0,0}N-1}\,\mathbf{D}\,/\,\pi$ and the diagonal matrix $\mathbf{D}$ is formed of the elements $\{\Delta x_{1},\,\Delta x_{2},\ldots,\Delta x_{N-1}\}$ . In other words, we can achieve the absolute value of the partial sum of particular RTF by using the measurement result on the sensor plane. Although the calculation of RTF function is classically quite difficult [6, 48], it has important applications in areas including geometry, arithmetic and number theory [47], nonlinear spectral theory for ocean and water sciences [48], cryptography and the solution modeling of NLSEs [60, 61].

Measurements in the QPC system allow to obtain information about $|q(x,t)|$ . The most important challenge for utilizing QPC in the solutions of NLSE is to determine the set equation parameters $(\mathbf{Y},\vec{k},\vec{\omega},\ \vec{\delta}^{\color[rgb]{0,0,0}\pm})$ which can be implemented with specific QPC design. Theoretical modeling and extensive simulations are required to examine all the practical sample parameters and systems in the literature where NLSE solutions are achieved with RTF based solution.

Next, numerical simulations are achieved to analyze the effects of FO based components such as lenses and HG sources on the intensity distribution obtained with QPC. Simulation studies for large scale implementations of QPC for the future applications are open issues.

Numerical Results

MPD set-up with two diffraction planes and single sensor plane is numerically analyzed for both Gaussian and HG sources with beam width and waist sizes of $\sigma_{0}=20\,\mu$ m and $W_{0}=200\,\mu$ m, respectively. The set-up is shown in Fig. 7 with $N=3$ planes. The source waveforms are shown in Figs. 8(a) and (b), respectively. HG order is set to $l=10$ with highly oscillatory and negative initial $V_{0}$ of $1.076$ . The wavelength of the light is chosen in the red spectrum of $\lambda=650$ nm while the low cost laser sources are commercially available in a wide spread manner.

Two different set-ups composed of LCT and FSP systems as shown in Figs. 7(a) and (b), respectively, are compared where the LCT system includes thin lenses between the diffraction planes while not included in the FSP system. The kernel based on HO in (47) giving FRFT as a special case of general LCT formulation is not numerically analyzed since LCT based system is a more general version while various combinations including FRFT systems are applicable with the formulation in Tables 1 and Matrix formulation for LCT system with Hermite-Gaussian sources in the Methods section. Therefore, two different set-ups with the kernels $K_{FS}$ in (48) and $K_{LCT}^{a,b,c,d}$ in (6) for the inter-plane propagation are compared for the same design of set-up in terms of the properties of the slits and planes. Inter-plane distance vector is given by $\vec{L}^{T}\equiv[31.5\,\,\,30.75\,\,\,0.9]$ (cm). The distances of the first plane to the first lens and the second plane to the second lens are denoted with the vector $\vec{L}_{a}^{T}\equiv[L_{a,01}\,\,\,L_{a,12}]$ where both the distances are set to $21$ cm while the lenses of the focal length $\vec{f}^{T}\equiv[f_{01}\,\,f_{12}]=[63\,\,63]$ mm focus the light intensity to more compact areas on the consecutive planes compared with FSP. It is assumed that the propagation between the second and third planes includes only FSP without any thin lens to simplify the set-up. $K_{1}=11$ and $K_{2}=27$ slits are used on the first and second planes, respectively. The slit positions and widths on the first plane, as shown in Figs. 8(c) and (d), respectively, are adapted to the maximum intensity locations of HG source propagation on the first plane in LCT system while the ones on the second plane are chosen uniformly with the separation of $40\,\mu$ m and the width of $\widetilde{\beta}_{2,j}\equiv\beta_{2}=8\,\mu$ m. $K_{LCT}^{a_{01},\,b_{01},\,c_{01},\,d_{01}}$ and $K_{LCT}^{a_{12},\,b_{12},\,c_{12},\,d_{12}}$ are calculated by using (13).

FSP has less control over the propagation of light compared with LCT based FO. FSP spreads the light without any tuning to the slit positions by reducing the probability of the photon to reach to the consecutive planes after diffraction. Therefore, in numerical analysis, LCT is shown to improve the probability of photon detection on the sensor plane ( $P_{E}$ ) and also the negative volume of Wigner function compared with FSP. The vectors of $P_{E}(j)$ and $V_{j}$ composed of the values on the first, second and third (sensor) planes for Gaussian sources are denoted with $\vec{P}_{E}^{G,FSP}$ and $\vec{V}^{G,FSP}$ , respectively, for FSP while with $\vec{P}_{E}^{G,LCT}$ and $\vec{V}^{G,LCT}$ for LCT. The cases with HG sources are denoted with the superscript of $HG$ . It is an open issue to adapt LCT parameters with respect to any given set-up including inter-plane distances, slit locations and widths in a way to maximize the interference and the probability of the photon reaching to the sensor plane.

Hermite-Gaussian Sources

The waveforms on the three planes in spatial domain are shown in Figs. 9(a), (b), (c) and 9(e), (f), (g) for FSP and LCT cases, respectively. It is observed that LCT focuses the light better on the slit locations while FSP reduces $P_{E}$ significantly. $\vec{P}_{E}^{HG,LCT}=[1\,\,0.515\,\,0.202]$ is much more improved compared with $\vec{P}_{E}^{HG,FSP}=[1\,\,0.042\,\,0.0062]$ . The magnitudes of the interfering paths are shown in Figs. 9(d) and (h), for $j=2$ and $j=3$ , respectively, while Wigner distributions on the second and third planes scaled with $\hbar$ are shown in Figs. 9(i) and (j) for FSP and, (k) and (l) for LCT. It is observed that LCT provides significantly larger path magnitudes while $\vec{V}^{HG,LCT}=[1.076\,\,1.47\,\,2.17]$ is also improved compared with $\vec{V}^{HG,FSP}=[1.076\,\,1.28\,\,1.25]$ . Observe that HG source has already negative Wigner volume of $1.076$ which is much further improved by LCT set-up compared with FSP.

Gaussian Sources

The waveforms in spatial domain for Gaussian sources are shown in Figs. 10(a), (b), (c) and 10(e), (f), (g) for FSP and LCT cases, respectively. $\vec{P}_{E}^{G,LCT}=[1\,\,0.35\,\,0.124]$ and $\vec{P}_{E}^{G,FSP}=[1\,\,0.077\,\,0.0161]$ values are obtained where the magnitudes of the interfering paths are shown in Figs. 10(d) and (h), for $j=2$ and $j=3$ , respectively. Similar to the HG sources, LCT improves diffraction probabilities significantly compared with FSP. Wigner distributions scaled with $\hbar$ are shown in Figs. 10(i), (j), (k) and (l) having different characteristics compared with HG case in Figs. 9(i), (j), (k) and (l). It is similarly observed that LCT provides significantly larger path magnitudes. $\vec{V}^{G,LCT}=[0\,\,0.842\,\,1.426]$ and $\vec{V}^{G,FSP}=[0\,\,1.21\,\,0.93]$ are obtained with increasing interference complexity through diffraction on consecutive planes in LCT case while starting with purely classical Gaussian source of zero negative Wigner volume. $V_{2}$ of $0.842$ for LCT is smaller than $1.21$ for FSP on the second plane. This is due to the both the specific set-up parameters and more diverse distribution of the path magnitudes in LCT after diffraction from the first plane as shown in Fig. 10(d). It becomes more difficult on the third plane to correlate the distribution of the path magnitudes shown in Fig. 10(h) with $V_{N}=V_{3}$ shown in Figs. 10(j) and (l). In other words, complexity behaves differently compared with the transmission probability while requiring simultaneous maximization depending on specific set-up as an open issue as discussed in the Results section while presenting Wigner distribution analysis.

Open Issues and Discussion

There are some open issues to best exploit photonic QPC method based on MPD and FO. Mathematical formulation correlating specific set-up parameters to path magnitude distribution and negative volume of Wigner function is an open issue. Iterative formulation of the vectors $\vec{\widetilde{h}}_{N-1,n}$ , $\vec{\widetilde{\gamma}}_{N-1,n}$ , $\vec{\widetilde{\eta}}_{N-1,n}$ and the matrices $\widetilde{\mathbf{H}}_{N-1,n}^{HO/G}$ , $\widetilde{\mathbf{H}}_{N-1,n}^{LCT/G}$ and $\widetilde{\mathbf{H}}_{N-1,n}^{LCT/HG}$ in (16), (25) and (27) are complicated as shown in the Methods section with complicated parameters and iterations in Tables 1 and Matrix formulation for LCT system with Hermite-Gaussian sources. Therefore, adapting the physical set-up parameters to the desired form of partial sum of RTF for the target number theoretical problem, and characterizing the path distributions and the negative volume of Wigner function explicitly are important open issues.

In the proposed formulation, the qubits are obtained through the tensor product structure of projections at different time instants on the contrary with the spatial encoding and entanglement of multiple photons. However, the histories of photon trajectories are not formulated for realizing conventional quantum gate implementations. Implementations of the quantum circuit gates are required to obtain universal QC architectures. Therefore, proposed QC formulation is limited to utilization of the interference of exponentially increasing number of Feynman paths based on the superposition and coherence properties of light source. Implementations of quantum circuits and fundamental search algorithms such as Grover search are future works to clarify the potential future scope of QPC based computing architectures in terms of universal QC capabilities. On the other hand, FO based QPC implementation has two main advantages resembling Boson sampling advantages in a different context [10]: (a) not utilizing multiple photons as qubits getting rid of the coupling disadvantages while exploiting single photon trajectories and (b) utilizing the free entanglement, for the first time, of the classical light obtained through freely available temporal correlations among the projections at different time instants. Boson sampling compared with QPC utilizes still multiple indistinguishable photons (but not as qubits) while requires single photon generation and detection to exploit the free entanglement among indistinguishable photons through multi-mode interferometer with the regarding boson statistics.

Realizing perfectly Gaussian slits compared with the conventional rectangular apertures is an important open issue for matching the experimental results with the proposed theoretical model. However, any slit structure can be represented as a composition of Gaussian slits by using the method defined in Ref. 72 and applied successfully in optical diffraction theory and experiments [73, 74]. The one dimensional slit mask function ${\color[rgb]{0,0,0}\widehat{G}}(x)$ is represented as follows:

[TABLE]

where $a_{i}$ and $\beta_{i}$ are found with optimization based on the experimental measurement results while increasing $K$ provides more accurate results. If the perfect Gaussian slits are replaced with the superposition in (40), then the summations in (16), (25) and (27) should be made for each $\beta_{i}$ of the single slit. The functional form with partial sum of RTF should be calculated and summed for each combination of $\beta_{i}$ through all the slits. Therefore, non-Gaussian slits can possibly realize the solutions of much harder computational complexity problems as an open issue.

There are some factors effecting the degree of compatibility between the theory and practice. These include imperfection in optical set-up, e.g., finite size lens effects, planar thickness, characterization of slit functions, sources and detector efficiency. The theoretical model should be extended including all the set-up parameters having diverging effects on the final intensity distribution. Similarly, the effects of exotic paths, i.e., trajectories between the slits on the same plane, should be included in the mathematical model as thoroughly discussed in Refs. 38, 75, 6. All these considerations potentially lead to unavoidable errors requiring quantum error correction studies adapted to QPC architectures [76]. Another open issue is related to the utilization of the measurements on all the sensor planes for computational purposes not only the final sensor plane since they include diffraction through previous planes. Theoretical models are required to exploit the sensor measurement results.

Methods

Quantum Fourier Optics

In scalar diffraction theory, the first Rayleigh-Sommerfeld formula of the Huygens-Fresnel principle for the propagation of light on planar surfaces is described as follows by using the Green’s theorem [77]:

[TABLE]

where $U_{I}(P_{1})$ is the wave amplitude at the point $P_{1}$ , $U$ is the distribution on the planar screen where diffraction occurs, $\Sigma$ denotes the integration over the slit including its multiplicative effects on the wave amplitude, $G_{-}\equiv exp(\imath\,k\,r_{01})\,/\,r_{01}-exp(\imath\,k\,\widetilde{r}_{01})\,/\,\widetilde{r}_{01}$ is the Green’s function vanishing on the diffraction surface for the first type of solution of Rayleigh-Sommerfeld formula, $r_{01}\equiv|\vec{r}_{0}-\vec{r}_{1}|$ and $k\equiv 2\,\pi\,/\,\lambda$ for the monochromatic light source of wavelength $\lambda$ as shown in Fig. 11 [77]. Assuming that $r_{01}\gg\lambda$ , the following approximation holds in rectangular coordinates:

[TABLE]

where the kernel $K_{FS}(\vec{r}_{1},\vec{r}_{0})$ for FSP is defined as follows:

[TABLE]

where $r_{01}=\sqrt{z^{2}\,+\,(x_{1}-x_{0})^{2}\,+\,(y_{0}-y_{0})^{2}}$ . The kernel for Fresnel diffraction integral is obtained by further approximation of $r_{01}$ in the near-field for large $z$ resulting in (45). This expression is the convolution integral conventionally used in phase-space optics for FSP [34].

Recently, scalar diffraction theory and Fresnel diffraction integral are discussed in Ref. 15 to be validly representing the evolution of light wave function modeled with the Hamiltonian of the quantized electromagnetic field $H=(\hat{p}^{2}+\omega^{2}\,\hat{q}^{2})\,/\,2$ as the Feynman’s path integral (FPI) solution of the quantum mechanical HO [78]. Fresnel diffraction nature of the propagation is verified with experimental photon counting studies for single photons. The wave function amplitude of light field in one dimension on a plane $\Psi(x_{0})$ is modeled to propagate into the amplitude $\Psi(x_{1})$ on another plane (Eq. 16 in Ref. 15 transformed into a simpler form) with the following formulation:

[TABLE]

where the kernel based on HO is the following:

[TABLE]

where $c$ is the velocity of light, $\omega\equiv 2\,\pi\,c\,/\,\lambda$ , $\omega\,t\neq n\,\pi$ for $n\in\mathbb{Z}$ , $t_{01}$ is the propagation duration between the planes and $m_{\lambda}\equiv\hbar\,k\,/\,c$ is the defined equivalent mass of photon propagation. In addition, the approximated FSP kernel in (45) is simply converted to the following in 1D system:

[TABLE]

The kernel for massive particles with the mass $m$ such as an electron is expressed as follows [78, 6]:

[TABLE]

In other words, the formulation based on phase-space optics for photon and electron propagation wave amplitudes have the similar form in (48) and (49) except an overall phase factor. The form in (48) is utilized in Ref. 33 for defining QSM while targeting only classical communications.

On the other hand, both the kernels $K_{HO}(x_{1},x_{0})$ and $K_{FS}(x_{1},x_{0})$ are special cases of LCTs defined for quadratic-phase optics [34]. As a class of linear integral transforms, they include as special cases the Fresnel transform and FRFT, simple scaling, chirp multiplication and some other operations. Spatial distribution of light in phase-space optics for the class denoted by quadratic-phase systems is mathematically equivalent to LCTs (Chapters 3 and 8 in Ref. 34). These optical systems include arbitrary combinations of the sections of free space in the Fresnel approximation, thin lenses and sections of quadratic graded-index media. In Ref. 15, FRFT nature of the kernel $K_{HO}(x_{1},x_{0})$ is shown both theoretically and experimentally while emphasizing the applicability of all the properties of Fourier analysis to quantum optics. In this article, propagation of the wave function is extended to the general case of LCTs providing flexibility to utilize arbitrary optical set-ups by enlarging the functional structures and number theoretical problems exploited in QPC. Furthermore, a better control is obtained for the energy flow of the light through the slits.

The kernel matrices for $K_{HO}(x_{1},x_{0})$ and $K_{FS}(x_{1},x_{0})$ are given as follows:

[TABLE]

$\mathbf{M}_{HO}$ has the same form with the propagation of light in quadratic graded-index media of having the refractive index distribution of $n^{2}(x)=n_{0}^{2}(1-(x\,/\,\chi)^{2})$ where $n_{0}$ and $\chi$ are the medium parameters. The parameter matrix of the propagation through the quadratic graded-index medium of length $d_{gri}$ is given by the following (Section 8.3.3 in Ref. 34):

[TABLE]

where $\alpha=d_{gri}\,/\,\chi$ . There is a FRFT relation between scaled versions of the input $\hat{f}(x)$ and output $\hat{g}(x)$ with FRFT order $\alpha$ as $\hat{g}(x)=e^{-\imath\,d_{gri}\,/\,(2\,\chi)}\lambda_{\chi}^{-1\,/4}f_{a}(x\,/\,\sqrt{\lambda_{\chi}})$ where $f(x)\equiv\lambda_{\chi}^{1\,/4}\,\hat{f}(x\,\sqrt{\lambda_{\chi}})$ and $f_{a}(x)$ denotes the $a$ th order FRFT of $f(x)$ . FRFT operation of order $\alpha$ is represented with the parameter matrix of $a=d=\cos(\alpha)$ and $b=\sin(\alpha)$ . As a result, $\mathbf{M}_{HO}$ represents a FRFT relation between the input and output scaled with the parameter $\sqrt{\lambda_{\chi}}$ where the parameters are $\alpha=\omega\,t$ and $\lambda_{\chi}\equiv 2\,\pi\,\hbar\,t_{01}\,/\,m_{\lambda}$ while as a special case of LCTs.

Matrix formulation for HO/LCT System with Gaussian sources

The following formulation is valid for both HO and LCT based design with Gaussian sources where the corresponding iteration parameters are defined in Table 1. The elements in the vector $\vec{\widetilde{h}}_{N-1,n}=\vec{\widetilde{c}}_{N-1,n}\,+\,\imath\,\vec{\widetilde{d}}_{N-1,n}$ are defined as follows:

[TABLE]

where $\vec{\widetilde{v}}_{k,j,n}$ for $k\in[0,j-1]$ is given as follows:

[TABLE]

Here, the matrix multiplication $\prod_{i=1}^{k}\mathbf{U}_{i}$ denotes $\mathbf{U}_{1}\,\mathbf{U}_{2}\ldots\mathbf{U}_{k}$ for any matrix $\mathbf{U}_{i}$ for $i\in[1,k]$ and $\widetilde{p}_{4,j,n}$ , $\widetilde{p}_{5,j,n}$ , $\widetilde{\zeta}_{j,c,n}$ and $\widetilde{\zeta}_{j,d,n}$ for $j\in[1,N-1]$ are defined in Table 1. Assume that $\mbox{diag}\{\vec{y}_{1},\ldots,\vec{y}_{K}\}$ and $\mbox{diag}\{\mathbf{y}_{1},\ldots,\mathbf{y}_{K}\}$ define the operators creating block diagonal matrices by putting the vectors $\vec{y}_{j}$ and the matrices $\mathbf{y}_{j}$ for $j\in[1,K]$ , respectively, (all the vectors or the matrices having the same dimensions) to the main diagonal and making zero the remaining elements. The matrix $\mathbf{\widetilde{H}}_{N-1,n}^{HO/G}$ is more simplified as follows compared with the more complicated form achieved for electron based FSP in Ref. 6:

[TABLE]

where the diagonal matrices are defined as follows:

[TABLE]

$2\times 2$ block $\widetilde{\mathbf{K}}_{j,b,n}$ and $1\times 2$ vector $\widetilde{\vec{k}}_{j,c,n}^{T}$ for $j\in[2,N-1]$ are defined as follows:

[TABLE]

$\widetilde{\mathbf{V}}_{N-1,n}$ is a lower triangular block matrix defined as follows:

[TABLE]

Expanding $\mathbf{\widetilde{H}}_{N-1,n}^{HO/G}$ in terms of real and imaginary parts is achieved by finding the real and imaginary parts of $\widetilde{p}_{1,j,n}$ for $j\in[1,N-1]$ and $\widetilde{p}_{3,j,n}$ for $j\in[2,N-1]$ , and $\widetilde{\mathbf{K}}_{j,b,n}$ and $\widetilde{\vec{k}}_{j,c,n}$ for $j\in[2,N-1]$ since $\widetilde{\mathbf{V}}_{N-1,n}$ is a real matrix. This is easily achieved by using the explicit forms of $\widetilde{p}_{1,j,n}$ and $\widetilde{p}_{3,j,n}$ in Table 1. Some variables and constants used in Table 1 (not defined in the table) are the following: $\imath\equiv\sqrt{-1}$ , $x_{j,n}$ denotes $X_{j,s_{n,j}}$ , $\lambda_{0}\,=\,\hbar\,t_{0,1}$ , $\alpha_{j}\,=\,\omega\,t_{j,j+1}$ for $j\in[0,N-1]$ , $\widehat{m}_{j}=m_{\lambda}\,/\,\sin(\omega\,t_{j,j+1})$ and $\widetilde{\lambda}_{j,n}^{*}$ as the conjugate of $\widetilde{\lambda}_{j,n}$ for real values of $(a_{j,j+1},b_{j,j+1},c_{j,j+1},d_{j,j+1})$ .

Polynomials for the example in (18-21) are presented in Table 2 for the simple case of $N=3$ and Gaussian source. It is possible by using the explicit modeling to make various gedanken experiments and to perform complexity theoretical calculations.

Matrix formulation for LCT system with Hermite-Gaussian sources

Iteration parameters utilized in (22-26) are presented in Table Matrix formulation for LCT system with Hermite-Gaussian sources. Some variables and constants used in Table Matrix formulation for LCT system with Hermite-Gaussian sources (not defined previously or for Table 1) are the following: $\widetilde{\tau}_{a,j,n}\,\equiv\,b_{j,j+1}\,\widetilde{u}_{j-1,j,n}\,+\,\imath\,\pi\,a_{j,j+1}$ for $j\in[2,N-1]$ , $\tau_{a,1}\,\equiv\,b_{12}\,{\color[rgb]{0,0,0}u}_{01}\,+\,\imath\,\pi\,a_{12}$ and $\widetilde{\Gamma}_{1,n}\,=\,2\,\widetilde{\beta}_{1,n}^{2}\,b_{1,2}\,{\color[rgb]{0,0,0}g}_{0,1}^{2}\,+\,\widetilde{\tau}_{1,n}$ . The formulation in (22) is obtained by using the integral equality for Hermite polynomials (Section 16.5 in Ref. 65) in an iterative manner along the planes:

[TABLE]

The parameters $\vec{\widetilde{\eta}}_{N-1,n}$ and $\widetilde{\gamma}_{j,j+1,n}$ in $\vec{\widetilde{\gamma}}_{N-1,n}^{T}\equiv\big{[}\widetilde{\gamma}_{12,n}\,\,\widetilde{\gamma}_{23,n}\,\,\ldots\,\,\widetilde{\gamma}_{N-1,N,n}\big{]}$ utilized in (23) are defined as follows by using the iterations in Table Matrix formulation for LCT system with Hermite-Gaussian sources:

[TABLE]

where the following are defined:

[TABLE]

where $\vec{\widetilde{\Xi}}_{N-2,n}^{T}$ is defined for $N>2$ and $\widetilde{\gamma}_{j,j+1,n}=\prescript{N-1}{}{\widetilde{\gamma}_{j,j+1,n}^{\star}}$ while $\prescript{l}{}{\widetilde{\gamma}_{j,j+1,n}^{\star}}$ and $\widetilde{\varepsilon}_{j,j+1,n}$ are defined as follows:

[TABLE]

where $\widetilde{\gamma}_{12,n}\equiv\widetilde{v}_{a,12,n}$ and $\widetilde{\varepsilon}_{12,n}\equiv\widetilde{h}_{a,12,n}$ for $N=2$ , and $\prescript{1}{}{\widetilde{\gamma}_{1,2,n}^{\star}}\equiv\widetilde{v}_{a,12,n}$ . Finally, $\widetilde{\mathbf{H}}_{N-1,n}^{LCT/HG}$ becomes equal to the following by using the iterations in Table Matrix formulation for LCT system with Hermite-Gaussian sources :

[TABLE]

where $\widetilde{\mathbf{H}}_{k,n}^{LCT/HG,\chi}$ is the following:

[TABLE]

Bibliography78

The reference list from the paper itself. Each links out to its DOI / PubMed record.

1[1] Arute, F. et al. Quantum supremacy using a programmable superconducting processor. \Journal Title Nature 574 , 505–510 (2019).
2[2] Cho, A. Google claims quantum computing milestone. https://science.sciencemag.org/content/365/6460/1364 (2019).
3[3] Preskill, J. Quantum computing in the NISQ era and beyond. \Journal Title Quantum 2 , 79 (2018).
4[4] Gomes, L. Quantum computing: Both here and not here. \Journal Title IEEE Spectrum 42–47 (2018).
5[5] Childs, A. M., Maslov, D., Nam, Y., Ross, N. J. & Su, Y. Toward the first quantum simulation with quantum speedup. \Journal Title Proceedings of the National Academy of Sciences 115 , 9456–9461 (2018).
6[6] Gulbahar, B. Quantum path computing: computing architecture with propagation paths in multiple plane diffraction of classical sources of fermion and boson particles. \Journal Title Quantum Information Processing 18 , 167 (2019).
7[7] Kok, P. et al. Linear optical quantum computing with photonic qubits. \Journal Title Reviews of Modern Physics 79 , 135 (2007).
8[8] Li, Y., Humphreys, P. C., Mendoza, G. J. & Benjamin, S. C. Resource costs for fault-tolerant linear optical quantum computing. \Journal Title Physical Review X 5 , 041007 (2015).