Next Article in Journal
Entropy Production in Quantum is Different
Next Article in Special Issue
Eigenvalues of Two-State Quantum Walks Induced by the Hadamard Walk
Previous Article in Journal
Non-Linear Dynamics Analysis of Protein Sequences. Application to CYP450
Previous Article in Special Issue
Can a Quantum Walk Tell Which Is Which?A Study of Quantum Walk-Based Graph Similarity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Directionally-Unbiased Unitary Optical Devices in Discrete-Time Quantum Walks

by
Shuto Osawa
1,*,
David S. Simon
1,2,* and
Alexander V. Sergienko
1,3,4,*
1
Department of Electrical and Computer Engineering, Boston University, 8 Saint Mary’s Street, Boston, MA 02215, USA
2
Department of Physics and Astronomy, Stonehill College, 320 Washington Street, Easton, MA 02357, USA
3
Department of Physics, Boston University, 590 Commonwealth Avenue, Boston, MA 02215, USA
4
Photonics Center, Boston University, 8 Saint Mary’s Street, Boston, MA 02215, USA
*
Authors to whom correspondence should be addressed.
Entropy 2019, 21(9), 853; https://doi.org/10.3390/e21090853
Submission received: 30 July 2019 / Revised: 22 August 2019 / Accepted: 29 August 2019 / Published: 31 August 2019
(This article belongs to the Special Issue Quantum Walks and Related Issues)

Abstract

:
The optical beam splitter is a widely-used device in photonics-based quantum information processing. Specifically, linear optical networks demand large numbers of beam splitters for unitary matrix realization. This requirement comes from the beam splitter property that a photon cannot go back out of the input ports, which we call “directionally-biased”. Because of this property, higher dimensional information processing tasks suffer from rapid device resource growth when beam splitters are used in a feed-forward manner. Directionally-unbiased linear-optical devices have been introduced recently to eliminate the directional bias, greatly reducing the numbers of required beam splitters when implementing complicated tasks. Analysis of some originally directional optical devices and basic principles of their conversion into directionally-unbiased systems form the base of this paper. Photonic quantum walk implementations are investigated as a main application of the use of directionally-unbiased systems. Several quantum walk procedures executed on graph networks constructed using directionally-unbiased nodes are discussed. A significant savings in hardware and other required resources when compared with traditional directionally-biased beam-splitter-based optical networks is demonstrated.

1. Introduction

The quantum approach to computing attracts public attention mainly because of its capability to execute some computational tasks faster when compared to classical computational devices [1,2]. Several physical platforms exist to realize quantum computation procedures. Linear optics has been one of the candidates because of its robustness against noise and the ease of quantum state manipulation at room temperature. The design of quantum computing gates with single photons has been proposed and is known as the Knill, Laflamme and Milburn (KLM) model [3]. This design makes use of linear-optical devices such as beam splitters and phase shifters. The quantum gate performance is executed probabilistically by the process of measuring auxiliary photons. While the KLM model has been used for gate-based quantum computation, other quantum-optical approaches to execute computational tasks have been developed. For example, quantum walks (QW) over optical networks of scattering centers have been considered as another promising tool in executing certain computational tasks [4,5,6,7,8]. The construction of such optical networks for quantum walks relies on the use of multiple beam splitters and phase shifters connected in a particular spatial graph pattern. A beam splitter is used as an elementary scattering center during the propagation, and many of them must be cascaded by connecting consecutively in order to form an extensive tree-like network [9,10]. Truly quantum mechanical information processing requires unitarity at every operation. The beam splitter in optics implements two-dimensional unitary transformations and can be seen as a probabilistic mixer of two spatial field modes.
The increase in dimensionality enables employing and manipulating more information, and this needs to be achieved in a coherent way. Optical networks are constructed to perform this task by constructing higher dimensional unitary matrices. It is known that higher dimensional unitary matrices can be decomposed using lower dimensional unitary matrices. By repeating this procedure, any complex unitary matrix can be eventually decomposed using only two-dimensional ones. The Reck decomposition model has been introduced to describe this procedure [11]. A symmetric version of the Reck model is often called the Clements model [12]. For instance, these two models have been used by researchers in designing and building experimental linear-optical networks for boson sampling purposes [13,14,15,16]. During the boson sampling process, photons propagate from one side of a complex nodal structure to the other side of the optical network, thus performing a computational task. Direct implementation of multimode optical device has been experimentally verified in integrated platforms [17,18,19,20]. Quantum walks over the network of quantum nodes represent another form of quantum information processing, as an alternative to the quantum gate model. QW can also perform certain computations more efficiently than classical algorithms [6,21,22,23,24,25,26,27,28]. Quantum walks in 1D and 2D systems have been experimentally demonstrated in optical systems [29,30,31,32,33,34,35,36,37,38,39,40].
The traditional quantum walk approach uses a coin operator and a shift operator to execute each elementary step. An alternative description of a quantum walk can be implemented using the scattering quantum walk, also known as the edge walk [41,42], which has been introduced to describe the quantum walk based on scattering at the nodes or vertices of a lattice on which the walk occurs. There is no need for a coin operator in this model. In order to execute some specific type of quantum walk, we need first to identify a network of scattering centers (a graph) on which the walk is performed. Many different special-purpose graphs can be formed using linear-optical devices in order to execute a particular computational procedure. Thanks to the Reck and Clements decomposition models, the majority of experimental demonstrations in this field, even some complex ones, could be realized using multiple directionally-biased two-dimensional optical devices such as beam splitters. However, the execution of such quantum walks calls for a large number of optical devices when the complexity and the required number of steps in the system increase. This is why quantum walks based on the use of directional devices demand a great deal of costly hardware real estate, which limits their scalability in the long run.
Recently, the original design of a directionally-unbiased linear-optical multiport was introduced [43]. This is a unitary coherent optical quantum information processing device that addresses two issues simultaneously: (i) it executes a higher dimensional unitary scattering process at every node of the network with fewer numbers of two-dimensional units for the device construction, and (ii) it scales down significantly the required amount of hardware resources by offering the possibility of reusing scattering units of the graph again and again. An array of such multiports can then form a graph upon which a photon can execute a quantum walk. In principle, the feature of full reversibility can be realized using special designs by incorporating commonly-used optical elements. This is referred to as “directional” or “directionally-biased” when a photon propagates only in one direction, meaning the input port and the output ports are never the same. This directionality could be circumvented in optics by placing mirrors so that a photon can leave the input port as well. This report will address multiple issues involved in designing, executing, testing, and applying both directional and directionally-unbiased devices. A higher dimensional quantum walk over a graph network based on the use of directionally-unbiased devices will be considered as an example of their practical applications.

2. Two-Dimensional Linear Optical Devices

Two-dimensional devices including interferometers are the main building blocks for any applications in classical and in quantum optics. These devices are unitary transformers that mix spatial optical modes without losses and realize the group of 2 × 2 unitary matrices denoted as U ( 2 ) . It has been shown that high-dimensional unitary matrices can be decomposed using U ( 2 ) matrices [44]. In order to have flexibility in quantum information processing, one needs to have some means of manipulating amplitude transition coefficients between the input and output fields. In principle, this could be achieved in two ways in optics: (i) by some kind of dynamic change in the input/output splitting ratio of a single beam splitter (BS) or (ii) by forming an interferometer with several beam splitters, thus offering tunability between output ports. In this section, we start with the basic properties of a beam splitter implementing the U ( 2 ) operation and discuss its features as a directionally-biased coupler. It will be followed by the consideration of integrated waveguided couplers and some well-known interferometers for implementing 2 × 2 transformations.

2.1. Lossless Optical Beam Splitter

A lossless beam splitter introduced in Figure 1 redirects incoming photons into two outgoing ports while maintaining energy conservation [45,46]. A BS can be represented using a 2 × 2 matrix, acting on two input and two output ports, denoted as E 1 , E 2 and E 3 , E 4 , as indicated in Figure 1. The transformation of the fields E 1 , E 2 is given by:
E 3 E 4 = T 13 R 23 R 14 T 24 E 1 E 2 ,
where T 13 , T 24 are the transmission from Port 1 to 3 and 2 to 4 and R 23 , R 14 are reflection for Port 2 to 3 and Port 1 to 4, respectively. The probability conservation relation between the input and output is:
| E 3 | 2 + | E 4 | 2 =   | E 1 | 2 + | E 2 | 2 .
By substituting Equation (1) in Equation (2),
| E 3 | 2 + | E 4 | 2 = ( | T 13 | 2 + | R 14 | 2 ) | E 1 | 2 + ( | T 24 | 2 + | R 23 | 2 ) | E 2 | 2 + T 13 R 23 * E 1 E 2 * + R 23 T 13 * E 2 E 1 * + T 24 R 14 * E 2 E 1 * + R 14 T 24 * E 1 E 2 * ,
and by comparing the result with Equation (2):
| T 13 | 2 + | R 14 | 2 =   | T 24 | 2 + | R 23 | 2 = 1 T 13 R 23 * + R 14 T 24 * = R 23 T 13 * + T 24 R 14 * = 0 .
Transmission and reflection coefficients T and R can be rewritten using amplitude and phase. Define T 13 = | T 13 | e i ϕ 13 , and so on, for all the transmission and reflection coefficients. Then, Equation (4) is reduced to:
| R 23 | | T 24 | = | R 14 | | T 13 | e i ( ϕ 14 + ϕ 23 ϕ 24 ϕ 13 ) .
In order to satisfy Equation (5), the phase values must be: ϕ 14 + ϕ 23 ϕ 24 ϕ 13 = ± π . This phase relation offers some flexibility in choosing the phase settings. Two different phase settings often appear in the literature for a beam splitter with a 50/50 power splitting ratio. When | R 23 | = | R 14 | = | R | = 1 2 ,   | T 13 |   =   | T 24 |   =   | T |   =   1 2 , one could choose ϕ 14 = ϕ 13 = ϕ 23 = 0 ,   ϕ 24 = π as an example. Other splitting ratios can be chosen as long as | T | 2 + | R | 2 = 1 is satisfied.
Example 1:
B S 1 = 1 2 1 1 1 1 .
Other phase settings could be ϕ 23 = ϕ 14 = ϕ R ,   ϕ 13 = ϕ 24 = ϕ T with | R 23 |   =   | R 14 |   =   | R | ,   | T 13 |   =   | T 24 |   =   | T | and substituting these in Equation (5).
| R | | T | = | R | | T | e 2 i ( ϕ T ϕ R ) ,
where ϕ R ϕ T = π 2 .
By choosing the phase settings ϕ T = 0 , ϕ R = π 2 , another example of the BS matrix can be produced.
Example 2:
B S 2 = 1 2 1 i i 1 .
Both examples are equivalent when appropriate phase shifters have been introduced before and after the beam splitter:
1 2 1 1 1 1 = 1 0 0 e i π 2 1 2 1 i i 1 1 0 0 e i π 2 .

2.2. Directionality of a Beam Splitter

A beam splitter is a symmetric device, meaning that any one of four ports can be used as an input, and its action is invariant under time reversal. At the same time, the device is not symmetric in the sense that the incoming photon cannot leave through the input port. We call this feature “directional-bias”; the choice of an input port biases the output to be in only two of the four possible output directions. This directional bias increases the required number of beam splitters when one enters the realm of higher dimensionality. In principle, this directional bias could be circumvented by placing external mirrors after the beam splitters so that they reverse the light propagation direction. This would allow the photon to leave through the input ports, and the system now becomes “directionally-unbiased”; all four output possibilities can still be realized, regardless of an input direction. Examples of ways to achieve this will be discussed in the coming sections.

2.3. 2 × 2 Integrated Directional Waveguide Coupler

A directional coupler is an integrated optics analog of a beam splitter. When two waveguides are brought close together, evanescent waves overlap and start coupling in the neighboring waveguide. Figure 2 illustrates a 2 × 2 integrated directional coupler and its cross-section. The coupling strength κ can be controlled by changing the distance between two waveguides. The BS and a directional coupler are both directionally-biased devices. The propagation of a photon through these devices can be described using a transfer matrix U. E o u t = U E i n , where E i n and E o u t are the input and output fields.
The transfer matrix of the directional coupler can be derived using coupled mode equations based on the Heisenberg equation [19,47]. The evolution in the z direction is given by:
i d A 1 d z = β A 1 + κ A 2 , i d A 2 d z = κ A 1 + β A 2 ,
where A j { j = 1 , 2 } are creation operators for a photon in the j th waveguide. β is a waveguide propagation constant, and κ is a coupling coefficient between the two waveguides.
Equation (10) can be rewritten in a matrix form:
d A 1 d z d A 2 d z = i β κ κ β A 1 A 2 .
We can solve for A 1 and A 2 by Equation (10) using differential equation solutions in the form of Equation (11) and finding eigenvalues and eigenvectors.
Eigenvalues with corresponding eigenvectors are given by:
λ 1 = β i κ i : 1 1 , λ 2 = β i + κ i : 1 1 ,
A 1 A 2 = c 1 e ( β z + κ z ) i 1 1 + c 2 e ( β z κ z ) i 1 1 ,
Initial conditions are given by: A 1 ( 0 ) = 1 , A 2 ( 0 ) = 0 , and c 1 = 1 2 , c 2 = 1 2 . The full transfer matrix can be reconstructed after solving also for the alternative initial condition: A 1 ( 0 ) = 0 , A 2 ( 0 ) = 1 .
U C o u p l e r = e β z i 2 e κ z i + e κ z i e κ z i e κ z i e κ z i e κ z i e κ z i + e κ z i = e β z i c o s ( κ z ) i s i n ( κ z ) i s i n ( κ z ) c o s ( κ z ) .

2.4. Interferometers as Two-Dimensional Devices

Interferometers are essential tools in quantum information processing and usually involve multiple beam splitters. The amplitude of each of the two outgoing modes can be modified by changing the relative phase between two paths. There are several major interferometer designs that offer 2 × 2 mode transformation. The Mach–Zehnder interferometer is a directionally-biased device that could be useful in realizing the Reck decomposition model, while the Michelson interferometer does not suffer from directional bias.

2.4.1. Mach–Zehnder Interferometer

The Mach–Zehnder interferometer shown in Figure 3 is a directionally-biased interferometer. Assume that each beam splitter has a 50/50 power splitting ratio between the two outgoing fields and the device is symmetric because the path length between two arms can be made identical. The beam splitter matrix U B S is applied twice, and the relative phase shift ϕ between the two modes by applying the matrix U p h a s e is introduced before the photon encounters the second beam splitter.
U M Z = U B S U p h a s e U B S .
By introducing expressions for the above matrices, one would obtain a specific formulation of the Mach–Zehnder transformation:
U M Z = 1 2 1 i i 1 e i ϕ 0 0 1 1 2 1 i i 1 = 1 2 e i ϕ 1 i ( e i ϕ + 1 ) i ( e i ϕ + 1 ) 1 e i ϕ = 1 2 e i ϕ 2 ( e i ϕ 2 e i ϕ 2 ) i e i ϕ 2 ( e i ϕ 2 + e i ϕ 2 ) i e i ϕ 2 ( e i ϕ 2 + e i ϕ 2 ) i e i ϕ 2 ( e i ϕ 2 e i ϕ 2 ) = e i ( ϕ 2 + π 2 ) s i n ( ϕ 2 ) c o s ( ϕ 2 ) c o s ( ϕ 2 ) s i n ( ϕ 2 ) .
The element-wise multiplication leads to the following input-output probability distribution:
P M Z = U M Z U M Z * = 1 2 1 c o s ϕ 1 + c o s ϕ 1 + c o s ϕ 1 c o s ϕ ,
where U * is a complex conjugate of U.
One can easily see that the Mach–Zehnder interferometer effectively serves as a tunable directionally-biased variable beam splitter, and this tunability plays a key role in higher dimensional interferometer-based optical networks.

2.4.2. Michelson Interferometer

The Michelson interferometer in Figure 4 is one example of the directionally-unbiased 2 × 2 device. Its layout could be used as an illustration of a general optical design principle that the directional bias within a 2 × 2 device could be circumvented by placing mirrors after the first beam splitter encounter and reversing directions of the optical flux.
A state of the input photon will be transformed by the first interaction with a BS.
E 3 E 4 = U B S E 1 E 2 ,
where U B S represents a linear transformation between the input fields and output fields of a beam splitter. In order to consider the transformation of fields from E 3 and E 4 to E 1 and E 2 , one must multiply the outcome with the inverse of matrix U B S from the left.
U B S 1 E 3 E 4 = U B S 1 U B S E 1 E 2 .
Since U B S is a unitary matrix U B S 1 = U B S , then:
E 1 E 2 = U B S E 3 E 4 .
This describes the transformation from E 3 and E 4 to E 1 and E 2 . This still represents a forward propagation; therefore, we must take a complex conjugate to reverse the propagation direction.
E 1 * E 2 * = U B S T E 3 * E 4 * .
This equation is read as a reverse propagation from E 3 and E 4 to E 1 and E 2 . Now, the 2 × 2 transformation of optical modes by the Michelson interferometer is described as:
U M i c h e l s o n = U B S T U P h a s e U B S .
The first transformation U B S describes the propagation from E 1 and E 2 to E 3 and E 4 , and then, phase shift U p h a s e introduces phase shifts between the two fields. The phases can be controlled by translating mirrors in the system as indicated in Figure 4. Finally, the reversed propagation and transformation from E 3 and E 4 to E 1 and E 2 is given by U B S T . These transformations complete the transformation of input fields by the Michelson interferometer. Such 2 × 2 directionally-unbiased devices based on the Michelson interferometer configuration could be used as elements for building higher dimensional interferometric systems. The Michelson interferometer is essentially the two-port version of the unbiased multiports introduced below, with U p h a s e providing the tunability.

3. Three- and Four-Dimensional Linear Optical Devices

The dimensionality of the 2 × 2 linear-optical device investigated in the previous section can be expanded to a more general situation covering a greater number of spatial modes. It has been shown in the past that one has to rely on using multiple 2 × 2 beam splitters in order to execute a high-dimensional transformation. This relationship is often called a Reck decomposition model (Reck model) [11]. The Reck model has been demonstrated experimentally [9]. There is also a symmetric directional alternative to Reck’s approach that is called Clements’ design [12,48]. This design can realize any unitary matrices and will be discussed in the 4 × 4 device section. In addition to Reck’s and Clements’ decomposition via multiple lower dimensional devices, a 3 × 3 directional transformation could be realized directly by exploiting a 3D optical integrated device in a waveguide configuration that is called an optical tritter [20,49,50]. Another decomposition model has been proposed as well [51]. This section examines these possible designs for 3 × 3 and 4 × 4 devices in detail. Four-dimensional devices are not just a simple extension of the three-dimensional devices. When the numbers of ports exceeds three, the distances between couplers are not identical. This means coupling strength would not be the same between couplers; therefore, it can change the final transfer matrix between the input fields and the output fields.

3.1. Reck Decomposition Design

It has been shown theoretically [11] that an arbitrary single N × N unitary matrix, U(N), can be decomposed into a succession of N ( N 1 ) 2 numbers of 2 × 2 mode mixing matrices. In order to understand the decomposition procedure, it is useful to understand the decomposition procedure for the 2 × 2 unitary matrix. Higher dimensional decomposition examples will be provided after the 2 × 2 example. An arbitrary unitary 2 × 2 matrix U(2) is defined as:
U ( 2 ) = A B C D ,
where A , B , C , D C . C is a set of complex numbers. It is always possible to find a unitary matrix T such that U(2) becomes diagonal after it is multiplied by the matrix T.
U ( 2 ) T = A 0 0 D ,
where A , D C .
The resulting diagonalized matrix will be turned into an identity matrix by multiplying it with an additional diagonal matrix P.
U ( 2 ) T P = 1 0 0 1 .
This procedure shows that any U(2) matrix can be transformed into an identity matrix. This result indicates that the inverse matrix ( T P ) 1 is the original U(2) we wanted. T and P are both unitary matrices; therefore, T 1 = T and P 1 = P where † is complex conjugate and transpose.
U ( 2 ) T P = I ( 2 ) U ( 2 ) = ( T P ) 1 = P T .
This procedure shows that a matrix U(2) is decomposed into matrices P and T . Such a diagonalization process can be applied in higher dimensions as well. In the 3 × 3 case, arbitrary U(3) matrices can be diagonalized using multiple 3 × 3 matrices with each matrix containing U(2) inside. T 3 , 1 , T 2 , 1 , and T 3 , 2 are the matrices containing U(2) inside. T 3 , 1 mixes spatial Modes 3 and 1; T 2 , 1 mixes spatial Modes 2 and 1; and T 3 , 2 mixes spatial Modes 3 and 2.
T 3 , 1 = U ( 2 ) 1 , 1 0 U ( 2 ) 1 , 2 0 1 0 U ( 2 ) 2 , 1 0 U ( 2 ) 2 , 2 ,   T 3 , 2 = 1 0 0 U ( 2 ) ,   T 2 , 1 = U ( 2 ) 0 0 1 ,
where U ( 2 ) 1 , 1 is the element of U(2) from the first row and the first column. The rest of three elements U ( 2 ) 1 , 2 , U ( 2 ) 2 , 1 , U ( 2 ) 2 , 2 follow the same rule. Our goal is to find a decomposition for an arbitrary matrix U(3). Assume that the U(3) matrix has the form of:
U ( 3 ) = A 1 B 1 C 1 D 1 E 1 F 1 G 1 H 1 I 1 .
As a first step, elements from the first row and the third column U ( 3 ) 1 , 3 and the third row and the first column U ( 3 ) 3 , 1 can be eliminated by multiplying a matrix T 3 , 2 .
U ( 3 ) T 3 , 2 = A 2 B 2 0 D 2 E 2 F 2 0 H 2 I 2 .
Repeat the elimination procedure for all the non-diagonal elements of U(3).
U ( 3 ) T 3 , 2 T 3 , 1 = A 3 B 3 0 D 3 E 3 0 0 0 I 3 ,
U ( 3 ) T 3 , 2 T 3 , 1 T 2 , 1 = A 4 0 0 0 E 4 0 0 0 I 4 .
U(3) is transformed into an identity matrix after multiplying by a diagonal matrix P. U(3) can be obtained by taking the inverse of ( T 3 , 2 T 3 , 1 T 2 , 1 P ) .
U ( 3 ) T 3 , 2 T 3 , 1 T 2 , 1 P = I ( 3 ) U ( 3 ) = ( T 3 , 2 T 3 , 1 T 2 , 1 P ) 1 = P T 2 , 1 T 3 , 1 T 3 , 2 .
Matrix elements A i through I i , i Z belong to C . Z is a set of integers. The U(3) matrix is decomposed into matrices P , T 2 , 1 , T 3 , 1 , and T 3 , 2 . This concludes that the knowledge of each individual beam splitter (or interferometer) in the system allows reconstructing a transfer matrix of the whole system. The same reconstruction process can be applied in the 4 × 4 case of U(4) decomposition.
U ( 4 ) T 4 , 3 T 4 , 2 T 4 , 1 T 3 , 2 T 3 , 1 T 2 , 1 P = I ( 4 ) U ( 4 ) = ( T 4 , 3 T 4 , 2 T 4 , 1 T 3 , 2 T 3 , 1 T 2 , 1 P ) 1 = P T 2 , 1 T 3 , 1 T 3 , 2 T 4 , 1 T 4 , 2 T 4 , 3 .
The matrices T i , j , i , j Z are 4 × 4 matrices, which contain U(2) matrices inside. The experimental setup in the case of 3 × 3 transformation is illustrated in Figure 5. The order of embedded U(2) matrices’ multiplication and their action is equivalent to the physical diagram outlined. The 4 × 4 case is given in Figure 6. Subfigures (a) and (b) in Figure 5 and Figure 6 are equivalent.

3.2. Clements Decomposition Design

The unitary matrix decomposition can also be realized in a slightly different configuration. The Clements design transforms the originally non-symmetric Reck configuration into a symmetric form [12], by which we mean that the situation is non-symmetric when photons in different input ports experience different numbers of beam splitters during their propagation and before exiting the unit. It would be helpful for any future consideration to introduce a simplified mesh representation for the systems outlined in Figure 7. For example, in the case of 4 × 4 transformation, its mesh decomposition via U(2) embedded matrices could be represented either by the original decomposition proposed by Reck. The crossing parts in the mesh designs mix spatial modes and consist of integrated couplers. The tunability of the power splitting ratio can be obtained either through the dynamical change of the coupling ratio between two waveguides in an integrated coupler or by forming an interferometer using two integrated couplers executed the same task. The graphical detail of the crossing parts in mesh design is indicated in Figure 8.
The mesh designs are illustrated in Figure 7. They are equivalent to the 4 × 4 Reck model and 4 × 4 Clements model. It is easy to note that a photon in the first path of the 4 × 4 Reck design illustrated in Figure 7a could encounter only one beam splitter, while the photon in the lowest path encounters at least three beam splitters prior to exiting the device. In the Clements symmetric design illustrated in Figure 7b, a photon in the first path and a photon in the last path encounter the same number of beam splitters. The loss tolerance of a quantum state becomes higher when a photon experiences the same number of beam splitter interactions [12]. The Clements designs are widely used instead of Reck decomposition, exactly because of its loss tolerance in quantum information processing.
In a similar way to the Reck model, a unitary matrix can be decomposed using multiple U(2)-based matrices. Unitary matrices’ realization based on the Reck model is decomposed by multiplying matrices from one side in succession. It is not necessary to multiply matrices from only one side to decompose the unitary matrix. The decomposition can be done by multiplying matrices from both sides. The U(3) case and U(4) case are given as an example.
U(3) case:
T 2 , 3 T 1 , 2 U ( 3 ) T 1 , 2 1 = P U ( 3 ) = T 1 , 2 1 T 2 , 3 1 P T 1 , 2 .
U(4) case:
T 3 , 4 T 2 , 3 U ( 4 ) T 1 , 2 1 T 3 , 4 1 T 2 , 3 1 T 1 , 2 1 = P U ( 4 ) = T 2 , 3 1 T 3 , 4 1 P T 1 , 2 T 2 , 3 T 3 , 4 T 1 , 2 .
The unitary matrix decomposition is possible, and this has been experimentally realized and demonstrated [12,52].

3.3. Integrated Optical Tritter and Quarter

Three- and four-dimensional directional linear optical devices will be introduced in this section. Integrated waveguide couplers [19,20,53,54] can be used to implement an optical tritter illustrated in Figure 9a and its cross-section shown in Figure 9b. The propagation dynamics of such a system can be described using the same formalism as in the previous case of the directional 2 × 2 coupler:
i d A 1 d z = β A 1 + κ A 2 + κ A 3 , i d A 2 d z = κ A 1 + β A 2 + κ A 3 , i d A 3 d z = κ A 1 + κ A 2 + β A 3 ,
where A j { j = 1 , 2 , 3 } are creation operators for a photon in the j th waveguide; β is a propagation constant in each waveguide; and κ is a coupling coefficient. It is assumed that all waveguides are identical and the distances between them are the same.
One could solve these equations using a matrix formalism in a similar way as the two-dimensional case:
d A 1 d z d A 2 d z d A 3 d z = i β κ κ κ β κ κ κ β A 1 A 2 A 3 .
One could find eigenvalues and corresponding eigenvectors:
λ 1 = β i + κ i : 1 1 0 , 1 0 1 , λ 2 = β i 2 κ i : 1 1 1 ,
A 1 A 2 A 3 = c 1 e β z i + κ z i 1 1 0 + c 2 e β z i + κ z i 1 0 1 + c 3 e β z i 2 κ z i 1 1 1 .
When the initial conditions are given by: A 1 ( 0 ) = 1 , A 2 ( 0 ) = 0 , A 3 ( 0 ) = 0 , then c 1 = c 2 = 1 3 , c 3 = 1 3 , where z is the propagation length. Other initial conditions are given as A 1 ( 0 ) = 0 , A 2 ( 0 ) = 1 , A 3 ( 0 ) = 0 , and A 1 ( 0 ) = 0 , A 2 ( 0 ) = 0 , A 3 ( 0 ) = 1 . After solving for all initial conditions, we can obtain a total transfer matrix for the system.
U I n t T r i t t e r = e β z i 3 2 e κ z i + e 2 κ z i e κ z i + e 2 κ z i e κ z i + e 2 κ z i e κ z i + e 2 κ z i 2 e κ z i + e 2 κ z i e κ z i + e 2 κ z i e κ z i + e 2 κ z i e κ z i + e 2 κ z i 2 e κ z i + e 2 κ z i .
Unlike two- and three-dimensional couplers where the distances between each pair of couplers are the same, for the four-dimensional coupler with its cross-section in Figure 10b, the coupling coefficients between diagonal coupling regions are different from those on the edges of the square. We again assume the coupling strength can be controlled.
i d A 1 d z = β A 1 + κ 1 A 2 + κ 2 A 3 + κ 1 A 4 , i d A 2 d z = κ 1 A 1 + β A 2 + κ 1 A 3 + κ 2 A 4 , i d A 3 d z = κ 2 A 1 + κ 1 A 2 + β A 3 + κ 1 A 4 , i d A 4 d z = κ 1 A 1 + κ 2 A 2 + κ 1 A 3 + β A 4 ,
where A j { j = 1 , 2 , 3 , 4 } are creation operators for a photon in the j th waveguide, β is a propagation constant, κ 1 is a coupling coefficient between two non-diagonal couplers, and κ 2 is a coupling coefficient for two diagonal couplers.
One could solve this combination of equations using matrix formalism
d A 1 d z d A 2 d z d A 3 d z d A 4 d z = i β κ 1 κ 2 κ 1 κ 1 β κ 1 κ 2 κ 2 κ 1 β κ 1 κ 1 κ 2 κ 1 β A 1 A 2 A 3 A 4 .
and by finding its eigenvalues and eigenvectors:
λ 1 = β i κ 2 i : 1 0 1 0 , 0 1 0 1 , λ 2 = β i 2 κ 1 i + κ 2 i : 1 1 1 1 , λ 3 = β i + 2 κ 1 i + κ 2 i : 1 1 1 1 ,
A 1 A 2 A 3 A 4 = c 1 e β z i κ 2 z i 1 0 1 0 + c 2 e β z i κ 2 z i 0 1 0 1 + c 3 e β z i 2 κ 1 z i + κ 2 z i 1 1 1 1 + c 4 e β z i + 2 κ 1 z i + κ 2 z i 1 1 1 1 .
When the initial conditions are given by: A 1 ( 0 ) = 1 , A 2 ( 0 ) = 0 , A 3 ( 0 ) = 0 , A 4 ( 0 ) = 0 , then c 1 = 1 2 , c 2 = 0 , c 3 = 1 4 , c 4 = 1 4 where z is the propagation length. After imposing the initial conditions and solving, we obtain a total transfer matrix for the system:
U I n t Q u a r t e r = e β z i 4 A B C B B C B A C B A B B A B C ,
where:
A = 2 e κ 2 z i + e ( 2 κ 1 + κ 2 ) z i + e ( 2 κ 1 + κ 2 ) z i , B = e ( 2 κ 1 + κ 2 ) z i + e ( 2 κ 1 + κ 2 ) z i , C = 2 e κ 2 z i + e ( 2 κ 1 + κ 2 ) z i + e ( 2 κ 1 + κ 2 ) z i .

4. Directionally-Unbiased Linear-Optical Designs

Directional devices were introduced in the previous sections. In this section, directionally-unbiased optical devices will be introduced. Some types of directionally-unbiased systems have existed in the literature for a long time [55,56], while new beam splitter-based designs have been introduced recently [43,57]. The new beam splitter-based design will be reviewed first, and then, a directionally-unbiased design based on an optical tritter will be introduced. These devices are the building blocks for quantum information processing applications, especially quantum walks. Three- and four-port designs will be reviewed in the following subsections.

4.1. Directionally-Unbiased Linear-Optical Three- and Four-Port Devices

Recently, several new designs for directionally-unbiased linear-optical multiports have been introduced, and they have applications in quantum simulations of Hamiltonians and topological phase simulations [43,58,59,60]. A three-port operation has been experimentally demonstrated using bulk optical devices [57]. The basic components of the multiport devices are beam splitters, mirrors, and phase shifters. The three-port operation is given in Figure 11. The entrance and exit ports are denoted as Port A, Port B, and Port C. A photon can enter any of the three ports, and the photon can leave any of the three ports. When a photon enters the system, the photon amplitude will be split at each beam splitter in the system. The final output photon amplitudes can be described by a transfer matrix constructed by adding all the possible paths the photon would take before it leaves the system. The coherence of the photon needs to be long enough to add all amplitudes coherently within the system. The four ports illustrated in Figure 12 will be introduced, as well as three ports in this subsection.
All input-output transfer elements for three ports could be considered as a coherent superposition of all possible paths inside the device:
A A = 1 4 e i ϕ C + 1 4 e i ϕ B i 1 8 e i ( ϕ B + ϕ C ) i 1 8 e i ( ϕ B + ϕ C ) + 1 16 e i ( ϕ A + ϕ B + ϕ C ) + 1 16 e i ( ϕ A + ϕ B + ϕ C ) 1 16 e i ( ϕ B + ϕ C + ϕ B ) 1 16 e i ( ϕ C + ϕ B + ϕ C ) 1 16 e i ( ϕ B + ϕ A + ϕ B ) 1 16 e i ( ϕ C + ϕ A + ϕ C ) + ,
A B = i 1 2 1 4 e i ϕ C i 1 8 e i ( ϕ A + ϕ B ) + i 1 8 e i ( ϕ A + ϕ C ) + i 1 8 e i ( ϕ B + ϕ C ) + 1 16 e i ( ϕ A + ϕ B + ϕ C ) 1 16 e i ( ϕ A + ϕ B + ϕ C ) 1 16 e i ( ϕ A + ϕ B + ϕ C ) + 1 16 e i ( ϕ C + ϕ B + ϕ C ) + 1 16 e i ( ϕ C + ϕ A + ϕ C ) + ,
A C = i 1 2 1 4 e i ϕ B + i 1 8 e i ( ϕ A + ϕ B ) i 1 8 e i ( ϕ A + ϕ C ) + i 1 8 e i ( ϕ B + ϕ C ) + 1 16 e i ( ϕ A + ϕ B + ϕ C ) 1 16 e i ( ϕ A + ϕ B + ϕ C ) 1 16 e i ( ϕ A + ϕ B + ϕ C ) + 1 16 e i ( ϕ B + ϕ C + ϕ B ) + 1 16 e i ( ϕ B + ϕ A + ϕ B ) + .
A A represents a transfer amplitude for the input A back to the output A. A similar geometric sum can be constructed for elements describing the photon coming in at B and leaving through Ports A, B, and C and, similarly, for an input C and the output through A, B, and C. Using the values above, a transfer matrix for this system can be reconstructed:
U m u l t i p o r t = U A A U B A U C A U A B U B B U C B U A C U B C U C C .
Here, U A A represents a transition from Port A to Port A. All terms are coherently summed; therefore, these transition amplitudes are for the long-time limit. The rest of the elements describe other possible (input → output) transitions.
The dimensionality and the number of optical elements is increased in the case of a four-port device as illustrated in Figure 12. It needs to be noted that four ports are slightly different from the three ports because of the numbers of beam splitter encounters before the input photon leaves the system. For example, the shortest path for Port A to Port B would be A B with one beam splitter encounter. Similarly, the shortest path for Port A to Port C would be A B C or A D C with two beam splitter encounters. The probability amplitude is lower for the path A to C because the photon encounters one extra beam splitter. This path-dependent amplitude difference needs to be considered for higher dimensional multiport implementation.
A A = 1 4 e i ϕ B + 1 4 e i ϕ D 1 16 e i ( ϕ B + ϕ C + ϕ B ) + 1 16 e i ( ϕ B + ϕ C + ϕ D ) 1 16 e i ( ϕ D + ϕ C + ϕ D ) + 1 16 e i ( ϕ D + ϕ C + ϕ B ) 1 16 e i ( ϕ B + ϕ A + ϕ B ) 1 16 e i ( ϕ D + ϕ A + ϕ D ) + ,
A B = i 1 2 i 1 8 e i ( ϕ D + ϕ C ) i 1 8 e i ( ϕ B + ϕ A ) + i 1 8 e i ( ϕ B + ϕ C ) + i 1 8 e i ( ϕ D + ϕ A ) + ,
A C = 1 4 e i ϕ B 1 4 e i ϕ D 1 16 e i ( ϕ B + ϕ C + ϕ D ) + 1 16 e i ( ϕ B + ϕ A + ϕ B ) + 1 16 e i ( ϕ D + ϕ C + ϕ D ) 1 16 e i ( ϕ D + ϕ C + ϕ B ) + 1 16 e i ( ϕ B + ϕ A + ϕ B ) + 1 16 e i ( ϕ D + ϕ A + ϕ D ) + ,
A D = i 1 2 i 1 8 e i ( ϕ B + ϕ C ) i 1 8 e i ( ϕ D + ϕ C ) i 1 8 e i ( ϕ D + ϕ A ) + i 1 8 e i ( ϕ B + ϕ A ) + .
The final transfer matrix consists of 16 (input → output) transition amplitudes:
U m u l t i p o r t = U A A U B A U C A U D A U A B U B B U C B U D B U A C U B C U C C U D C U A D U B D U C D U D D .

4.2. Constructing Reversible Optical Tritter and Quarter

It has been shown earlier that traditional integrated optical 3 × 3 and 4 × 4 couplers are directionally-biased devices. However, placing mirrors at each of the output ports of the device helps to eliminate this directional bias and return the optical signal back to any of the input ports. This type of reversible design has been also introduced in the area of linear interferometric networks and has been sometimes called a generalized Michelson interferometer [61,62]. The reversibility is introduced by mirrors, while additional phase shifters can be introduced before the mirrors. Both 3 × 3 and 4 × 4 implementations are illustrated in Figure 13. The reversed tritter can be realized using the same design formalism as a generalized Michelson interferometer. The input photon state is transformed by a tritter matrix, phase shifters, and a transposed tritter matrix.
U R e v T r i t t e r = U I n t T r i t t e r T U p h a s e U I n t T r i t t e r ,
U p h a s e = e i ϕ A 0 0 0 e i ϕ B 0 0 0 e i ϕ C ,
U I n t T r i t t e r = e β z i 3 2 e κ z i + e 2 κ z i e κ z i + e 2 κ z i e κ z i + e 2 κ z i e κ z i + e 2 κ z i 2 e κ z i + e 2 κ z i e κ z i + e 2 κ z i e κ z i + e 2 κ z i e κ z i + e 2 κ z i 2 e κ z i + e 2 κ z i = e β z i 3 A B B B A B B B A ,
where A = 2 e κ z i + e 2 κ z i and B = e κ z i + e 2 κ z i . Putting all the pieces together, the final expression is:
U R e v T r i t t e r = e 2 β z i 9 e i ϕ A A A + e i ϕ B B B + e i ϕ C B B e i ϕ A A B + e i ϕ B B A + e i ϕ C B B e i ϕ A A B + e i ϕ B B B + e i ϕ C B A e i ϕ A B A + e i ϕ B A B + e i ϕ C B B e i ϕ A B B + e i ϕ B A A + e i ϕ C B B e i ϕ A B B + e i ϕ B A B + e i ϕ C B A e i ϕ A B A + e i ϕ B B B + e i ϕ C A B e i ϕ A B B + e i ϕ B B A + e i ϕ C A B e i ϕ A B B + e i ϕ B B B + e i ϕ C A A .
U p h a s e = e i ϕ A 0 0 0 0 e i ϕ B 0 0 0 0 e i ϕ C 0 0 0 0 e i ϕ D ,
where ϕ A , ϕ B , ϕ C and ϕ D are phase shifts introduced before the second device encounter.
U I n t Q u a r t e r = e β z i 4 A B C B B C B A C B A B B A B C ,
where A = 2 e κ 2 z i + e 2 κ 1 z i + κ 2 z i + e 2 κ 1 z i + κ 2 z i , B = e 2 κ 1 z i + κ 2 z i + e 2 κ 1 z i + κ 2 z i , and C = 2 e κ 2 z i + e 2 κ 1 z i + κ 2 z i + e 2 κ 1 z i + κ 2 z i . The reversible quarter matrix is then derived from the equation below.
U R e v Q u a r t e r = U I n t Q u a r t e r T U p h a s e U I n t Q u a r t e r .

5. Discrete-Time Quantum Walks

Up to this section, photonics-based linear optical devices have been investigated. A major motivation for focusing on the directionally-unbiased versions is in their potential for implementing quantum walk applications. Quantum walks are motivated from classical random walks. Quantum walks can support superposition states and interference in the system where interference is absent in classical random walks. There are two types of quantum walks, discrete-time quantum walks and continuous-time quantum walks. During the discrete-time quantum walk, the evolution operator is applied in a discrete time fashion while the operator application timing is irrelevant in the continuous case. We focus on the discrete case in this review. The simplest classical and quantum walk design would be a walk performed on a line. In the case of a classical random walk on a line with an unbiased two-dimensional coin, a walker can hop one step to the right or to the left with equal probability depending on the result of a coin-toss event. The walker walks on a line for certain steps, and the probability at a specific position can be obtained by repeating the process. By recording all the probability at each location on the line, a probability distribution associated with that coin is constructed. Classical random walks involve intermediate measurement, meaning the position of the walker is measured right after a coin-toss event. In contrast, the quantum approach to random walks preserves the coherence of all possible paths by not measuring an intermediate state of the walker and, as a consequence, enables the quantum interference of available probability amplitudes. It is known that the probability distribution spreads faster in quantum walks compared to classical random walks. Classical random walks are useful for many randomized algorithm implementations [63]. It is natural to consider that quantum walks could achieve better outcomes than classical random walk-based algorithms, and it is indeed possible to gain algorithmic speedup using the fact that quantum walks can spread faster than classical random walks. Several different algorithms have been developed through quantum walks, and some are faster than classical algorithms. Hitting time, graph traversal speed from a point to another point in a graph, on a hypercube [21,22], and a glued tree are known to be exponentially faster in the quantum case [23,24]. Element distinctness [25], triangle finding [26], matrix product verification [27], and group commutativity testing [28] have been also investigated. Flexible graph construction is necessary to perform quantum walk-based algorithms. Any graphs consist of vertices and edges, and these need to be prepared in an experimentally realizable way. This task can be achieved through linear-optical devices, which have several input and output ports as discussed in previous sections. An experimental quantum walk implementation has been demonstrated in optical systems using optical cavities [29,30], optical rings [31], time-bins [32], Michelson interferometers [33], optical network [34], beam displacers [35], orbital angular momentum manipulation [36,37,38], and optical refraction [64]. The majority of these implementations are based on directional-optical devices; therefore, their implementation costs would rapidly increase as the dimensionality of the quantum walk system becomes higher. This applies to spatially-multiplexed quantum walk systems as they need to use beam splitters in a feed-forward manner. Time-multiplexed quantum walks are also commonly used since they can be compact. However, it would be challenging to perform node-by-node amplitude tuning. Integrated waveguide-based systems can be made directionally-unbiased and have been experimentally demonstrated [39,65,66]. Directionally-unbiased linear-optical multiport-based quantum walk configurations, which can realize amplitude tunability while offering an implementation resource reduction, will be introduced in the following several subsections.

5.1. Coin Walk: Quantum Walk on Vertices

The traditional quantum walk is illustrated using a position Hilbert space H P and a “coin” Hilbert space H C . A quantum walker’s position is described by the amplitudes in a position space spanned by { | m , m Z } , and a coin space is spanned by a two-dimensional computational basis { | R ( 1 , 0 ) T , | L ( 0 , 1 ) T } . The Hilbert space of the system is given by H = H P H C . We define a coin operator C ^ and a shift operator S ^ acting on each Hilbert space. The shift operator translates a walker’s position from |m〉 to | m 1 or | m + 1 depending on the result of the coin operation.
S ^ | m | R = | m + 1 | R and S ^ | m | L = | m 1 | L .
We can deduce a linear operator S ^ .
S ^ = n = | m + 1 m | | R R | + n = | m 1 m | | L L | .
The walker’s direction of the walk is decided by the result of the coin operator. The walk consists of applying, at each step, the coin operator, then the shift operator. The combined operation is given by:
V ^ = S ^ · ( I ^ C ^ ) .
This V ^ is applied on an initial state multiple times to perform walks with multiple steps.
| ψ ( t = N ) = V ^ N | ψ ( t = 0 ) .
The Hadamard coin operator H 2 ^ can be used to demonstrate the quantum walk. The coin operator C ^ in Equation (63) is substituted by H 2 ^ .
H 2 ^ = 1 2 1 1 1 1 .
One cycle of the quantum walk is completed by applying H ^ I ^ followed by S ^ . This process is performed multiple times without making intermediate measurements. Final measurements are made after a certain number of time steps. The probability distribution generated using a quantum walk behaves differently than a classical random walk. The standard deviation of the classical random walk on a line with N step is known to have a size of N [4]; on the other hand, a quantum walk on a line has a standard deviation of the order of N. This indicates that the quantum walk spreads faster than the classical random walk and can result in large speed increases in searching applications.

5.2. Scattering Quantum Walk: Quantum Walk on Edges

A different picture of the quantum walk is provided by the scattering model introduced [41,42]. This discrete-time scattering-based quantum walk is also called an edge walk. Unlike the coin model, the interference occurs on edges instead of performing the walk only on vertices. Each vertex works as a scattering center in this model. An input photon amplitude and phase will be controlled by a transmission and reflection coefficient at the scattering center. This model starts with a photon in a state | m 1 , m , representing a photon propagating from a vertex location m 1 to m, hence describing a state on an edge.
In this scattering model, a Michelson interferometer can serve as a vertex with two edges. One-step propagation starting at a specific edge is given by:
U ^ M i c h e l s o n | m 1 , m 1 2 ( | m , m + 1 + i | m , m 1 ) . U ^ M i c h e l s o n | m + 1 , m 1 2 ( | m , m 1 + i | m , m + 1 ) .
We will introduce a simplified description for a single scattering center, but we present the full description first. The unitary transformation represented by the Michelson interferometer in Equation (66) and illustrated in Figure 14b can be rewritten using the matrix below:
U F u l l = A R A R A L A R B R A R B L A R A R A L A L A L B R A L B L A L A R B R A L B R B R B R B L B R A R B L A L B L B R B L B L B L = 1 2 0 i 1 0 i 0 0 1 1 0 0 i 0 1 i 0 .
U ^ F u l l | A R 1 2 ( | B R + i | A L ) , U ^ F u l l | B L 1 2 ( | A L + i | B R ) .
| A R , | A L , | B R , and | B L correspond to | m 1 , m , | m , m 1 , | m , m + 1 , and | m + 1 , m , respectively. The propagation direction needs to be distinguished when multiple scattering centers are connected, but a simplified version can be used for a single scattering center. We will use the simplified matrix for a single scattering element in the upcoming sections so that we can directly make a comparison to the coin operators:
U ^ S i m p l i f i e d = A A B A A B B B = 1 2 i 1 1 i
By repeating this unitary matrix transformation process from an initial state, a walk can be implemented on a line. This scattering model is unitarily equivalent to the coin walk [67,68]. To see the unitary equivalence between the two models, we define a unitary operator E ^ .
E ^ | m 1 , m = | m | R , E ^ | m + 1 , m = | m | L ,
where |R〉 and |L〉 are defined in the coin model section. Consider a state evolution by operators E ^ U ^ with an initial condition | m 1 , m .
U ^ | m 1 , m = 1 2 ( | m , m + 1 + i | m , m 1 ) , E ^ U ^ | m 1 , m = 1 2 ( | m + 1 | R + i | m 1 | L ) .
V ^ E ^ | m 1 , m = V ^ | m | R = S ^ | m 1 2 ( | R + i | L ) = 1 2 ( | m + 1 | R + i | m 1 | L ) .
The former represents the edge walk, and the latter represents the coin walk. The outcomes are the same when the evolution operators U ^ and V ^ are multiplied by the operator E ^ ; therefore, these two formalisms are unitarily equivalent. This unitary equivalence E ^ U ^ = V ^ E ^ can be also seen as U ^ = E ^ V ^ E ^ . U ^ n = E ^ V ^ n E ^ because of unitarity of the operator E ^ where n is an integer. This result can be extended to higher dimensional walks. We can find the same equivalence for an initial state | m + 1 , m .

5.3. Higher Dimensional Coin Operators and Scattering Vertices

Quantum walks can be extended to higher dimensions by changing the dimension of the operators in the system and attaching additional edges to each vertex. It is possible to introduce scattering centers with different scattering amplitude ratios between output modes using directionally-unbiased devices. We will introduce several different coin operators and corresponding scattering centers in this section. The relationship between the coin model and the scattering model is deduced using an additional unitary operator as discussed in the previous subsection. There are several quantum coin operators with specific characteristics. The Hadamard coin, an unbiased coin, is one example.
The four-dimensional real-valued Hadamard coin H 4 is given as an example. This matrix is obtained by taking tensor product of two two by two real Hadamard matrices H 2 .
H 4 = H 2 H 2 = 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 .
In addition to the Hadamard coin, there are two other major specific coins used in quantum information processing. The first coin is motivated by Grover’s search algorithm [2].
C d = 2 d 1 2 d 2 d 2 d 2 d 1 2 d 2 d 2 d 2 d 1 ,
where d is the size of the matrix.
Matrices for d = 3 and 4 are given.
C 3 = 1 3 1 2 2 2 1 2 2 2 1 and C 4 = 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 .
This coin is biased in amplitudes (except for d = 4), yet symmetric under permutations of matrix labels. Another coin is a discrete Fourier transform (DFT) coin; this coin is unbiased; however, it is not symmetric under permutations. The Fourier transform matrix is given by:
U F o u r i e r = 1 d 1 1 1 1 1 ω ω 2 ω d 1 1 ω 2 ω 4 ω 2 ( d 1 ) 1 ω ( d 1 ) ω 2 ( d 1 ) ω ( d 1 ) ( d 1 ) ,
where ω = e 2 π i d .
Three-dimensional and four-dimensional Fourier coins are given by:
U F o u r i e r = 1 3 1 1 1 1 ω 3 ω 3 2 1 ω 3 2 ω 3 ,
where ω 3 = e 2 π i 3 .
U F o u r i e r = 1 2 1 1 1 1 1 ω 4 ω 4 2 ω 4 3 1 ω 4 2 ω 4 4 ω 4 6 1 ω 4 3 ω 4 6 ω 4 9 = 1 2 1 1 1 1 1 i 1 i 1 1 1 1 1 i 1 i ,
where ω 4 = e 2 π i 4 .

5.4. Equivalence between Higher Dimensional Coin Walk and Scattering Quantum Walk

The coin walk and the scattering walk were introduced in the previous subsections, as well as higher dimensional coin operators. It is possible to give a unitary equivalence relation between the two walks in higher dimensions as well. Consider a quantum walk on a 2D rectangular lattice. The center of the grid is given by coordinate ( m , n ). A photon on one of the edges around that grid is defined as | m , m , n 1 , n . This state is read as a photon propagation from a vertex location ( m , n 1 ) to ( m , n ). The unitary operator for one propagation step of an edge is defined as:
U ^ | m 1 , m , n , n 1 2 ( | m , m 1 , n , n + | m , m , n , n + 1 | m , m + 1 , n , n + | m , m , n , n 1 ) , U ^ | m , m , n + 1 , n 1 2 ( | m , m 1 , n , n + | m , m , n , n + 1 + | m , m + 1 , n , n | m , m , n , n 1 ) , U ^ | m + 1 , m , n , n 1 2 ( | m , m 1 , n , n + | m , m , n , n + 1 + | m , m + 1 , n , n + | m , m , n , n 1 ) , U ^ | m , m , n 1 , n 1 2 ( | m , m 1 , n , n | m , m , n , n + 1 + | m , m + 1 , n , n + | m , m , n , n 1 ) .
The corresponding coin operator of the coin walk is given by:
C ^ = 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 .
The operator transforms an initial coin state into a superposition state:
C ^ | L = 1 2 ( | L + | U + | R + | D ) ,
where | L = ( 1 , 0 , 0 , 0 ) T ,   | U = ( 0 , 1 , 0 , 0 ) T ,   | R = ( 0 , 0 , 1 , 0 ) T ,   | D = ( 0 , 0 , 0 , 1 ) T . A new shift operator is defined as follows:
S ^ = m n ( | m , n + 1 m , n | | U U | + | m , n 1 m , n | | D D | + | m + 1 , n m , n | | R R | + | m 1 , n m , n | | L L | ) .
One step of the coin walk is given by:
V ^ = S ^ ( I ^ C ^ ) .
We wish to find equivalence between the edge walk and the coin walk by finding a unitary operator E ^ . Define an operator E ^ transforming edge states into vertex states.
E ^ | m 1 , m , n , n = | m , n | R E ^ | m , m , n + 1 , n = | m , n | D E ^ | m + 1 , m , n , n = | m , n | L E ^ | m , m , n 1 , n = | m , n | U .
Consider two cases for the coin-based walk and edge walk starting with an initial state | m , m , n 1 , n .
Unitary transformation of the coin-based walk:
E ^ | m , m , n 1 , n = | m , n | U .
V ^ E ^ | m , m , n 1 , n = S | m , n ( 1 2 ( | L | U + | R + | D ) = S 1 2 ( | m , n | L | m , n | U + | m , n | R + | m , n | D ) = 1 2 ( | m 1 , n | L | m , n + 1 | U + | m + 1 , n | R + | m , n 1 | D ) .
Unitary transformation of the edge-based walk:
U ^ | m , m , n 1 , n = 1 2 ( | m , m 1 , n , n | m , m , n , n + 1 + | m , m + 1 , n , n + | m , m , n , n 1 ) .
E ^ U ^ | m , m , n 1 , n = 1 2 ( | m 1 , n | L | m , n + 1 | U + | m + 1 , n | R + | m , n 1 | D ) .
V ^ E ^ and E ^ U ^ both transform the initial state into the same state. Therefore, the outcomes are equivalent, and the coin walk and the scattering walk are unitarily equivalent. We went through a specific equivalence, which is the quantum walk generalized equivalence between two models, as found elsewhere [41,67].

5.5. Examples of Multi-Dimensional Quantum Walks on Graphs

The two types of quantum walks, the coin quantum walk and the scattering (edge) quantum walk, are both performed on graphs. Graphs with nodes implemented by higher dimensional coins are applicable to algorithm development. For example, the Grover search algorithm, when implemented via quantum walks on certain graphs with a superposition initial state, demonstrates significant speedup over classical algorithms. Many quantum walk applications are based on undirected graphs, meaning a walker can travel forward and backward in the system. Directionally-unbiased linear-optical devices possess reversibility and therefore can implement such undirected graphs. It is shown that a spatial search performed on a 2D lattice is faster than similar classical algorithms [69,70,71,72,73,74]. To observe the spatial search on a 2D lattice, the scattering centers have transmission and reflection coefficients equal to the Grover coin setting. One node in a graph is “marked” by introducing a different matrix on one specific scattering center in the lattice. Localization occurs on edges around the marked scattering center when the superposition state is sent in the system as an initial state. The graph geometry can be configured using directionally-unbiased devices. The rectangular lattice illustrated in Figure 15a would require four-port devices. Similarly, the hexagonal lattice illustrated in Figure 15b would require three-port devices. These optical quantum walk implementations through these optical devices are advantageous because of their amplitude tunability. The graphs can be implemented with the same Grover coin matrix setting throughout the graph vertices initially, then a marked coin can be introduced by tuning one of the vertices into a different coin. A quantum walk search can find the marked point faster than any classical search algorithms.
Other graph structures can be considered using directionally-unbiased devices. Quantum walks on a glued tree have been investigated theoretically and experimentally [39,75]. Hitting time, the time required to reach one point to another point on a graph, is commonly used to evaluate propagation speed on a specific graph. A quantum walk on a glued tree with three nodes gives exponential speedup when the three-dimensional Grover coin is used at the nodes [76]. We can form hypercubes using unbiased multiports, and it has been shown that the quantum walk hitting time is shorter than the classical walk case [22]. A walker starts on the left side of the graph, and the walker tries to reach the other end in a short amount of time. As indicated in Figure 16a, the randomly-connected middle part in the glued tree complicates the path finding procedure to reach the other end. The speedup applies in the case of the hypercube as well. A four-dimensional hypercube is illustrated in Figure 16b as an example. Classical algorithms cannot perform this search efficiently. On the other hand, quantum walks can perform exponentially faster than any classical algorithms to find the other end of the path.

6. Specific Transfer Matrix Examples Using Reversible Linear-Optical Devices

We consider several specific experimental configurations for an efficient realization of quantum walks in higher dimensions using linear optical devices and exploiting the very important feature of optical reversibility. We look into specific phase values and corresponding transfer matrices using formalisms covered in previous sections. The focus here is on realization of the Fourier coin and the Grover coin using directionally-unbiased linear-optical multiports, and reversible optical tritter and quarter configurations.

6.1. The Fourier Coin Realization

A three-dimensional Fourier coin has the form (see Section 5.3):
U F o u r i e r = 1 3 1 1 1 1 ω 3 ω 3 2 1 ω 3 2 ω 3 ,
where ω 3 = e 2 π i 3 .
This matrix can be generated with a reversible tritter containing phase shifters at all three ports, and directionally-unbiased linear-optical three-ports can perform the same job as well. For a reversible system, an input photon experiences the same phase shifts twice from the same phase shifters. U p h a s e i n and U p h a s e o u t would take care of such phase shifts. Figure 17 is a reversible tritter with phase shifters. Multiport designs are introduced in Section 4.1.
The matrix generation is performed by multiplying U p h a s e i n , U d e v i c e , and U p h a s e o u t in sequence. U d e v i c e can be any reversible optical device. The Fourier coin using the optical tritter is realized by the following procedure.
U F o u r i e r = U p h a s e o u t U T r i t t e r U p h a s e i n .
U p h a s e i n = U p h a s e o u t = e i ϕ a 0 0 0 e i ϕ b 0 0 0 e i ϕ c ,
U F o u r i e r = e i π 3 0 0 0 e i π 3 0 0 0 e i π 3 i 3 e i 2 π 3 1 1 1 e i 2 π 3 1 1 1 e i 2 π 3 e i π 3 0 0 0 e i π 3 0 0 0 e i π 3 = i 3 1 1 1 1 e i 4 π 3 e i 2 π 3 1 e i 2 π 3 e i 4 π 3 ,
where ϕ a , ϕ b , ϕ c are the phase shifts from phase shifters at the entrance ports. The reversible optical tritter and quarter can realize the Fourier coin when phases are set at specific values:
( ϕ A , ϕ B , ϕ C , κ z , ϕ a , ϕ b , ϕ c ) = ( 10 π 9 , 10 π 9 , 10 π 9 , 10 π 9 , π 3 , π 3 , π 3 ) ,
where ϕ A , ϕ B , ϕ C , and κ z are phase values from mirror units of the reversible tritter and propagation distance of the coupling region.
Similarly, the unbiased three-port operation can realize the same matrix with the settings:
U F o u r i e r = U p h a s e o u t U T h r e e p o r t U p h a s e i n .
U F o u r i e r = e i π 3 0 0 0 e i π 3 0 0 0 e i π 3 1 3 e i 2 π 3 e i 2 π 3 1 1 1 e i 2 π 3 1 1 1 e i 2 π 3 e i π 3 0 0 0 e i π 3 0 0 0 e i π 3 = 1 3 e i 2 π 3 1 1 1 1 e i 4 π 3 e i 2 π 3 1 e i 2 π 3 e i 4 π 3 ,
when:
( ϕ A , ϕ B , ϕ C , ϕ a , ϕ b , ϕ c ) = ( π 6 , π 6 , π 6 , π 3 , π 3 , π 3 ) ,
where ϕ A , ϕ B , and ϕ C are the phase values for the directionally-unbiased three-port operation.
Four-dimensional matrices are generated using the same methods. The reversible quarter with phase shifters is given in Figure 18. Multiport designs are introduced in Section 4.1.
The four-dimensional Fourier coin has the form of:
U F o u r i e r = 1 2 1 1 1 1 1 ω 4 ω 4 2 ω 4 3 1 ω 4 2 ω 4 4 ω 4 6 1 ω 4 3 ω 4 6 ω 4 9 ,
where ω 4 = e 2 π i 4 .
U F o u r i e r = U p h a s e o u t U Q u a r t e r U p h a s e i n ,
U p h a s e i n = U p h a s e o u t = e i ϕ a 0 0 0 0 e i ϕ b 0 0 0 0 e i ϕ c 0 0 0 0 e i ϕ d .
The reversible quarter can realize the Fourier coin as well with the phase settings equal to:
( ϕ A , ϕ B , ϕ C , ϕ D , κ z 1 , κ z 2 , ϕ a , ϕ b , ϕ c , ϕ d ) = ( π , π 4 , π , 5 π 4 , 7 π 4 , 7 π 8 , π 4 , π 2 , π 4 , π 2 ) ,
where ϕ A , ϕ B , ϕ C , ϕ D , and κ z 1 , κ z 2 are phase values from the mirror units of the reversible tritter and propagation distance of the coupling region.
Four ports:
U F o u r i e r = U p h a s e o u t U f o u r p o r t U p h a s e i n ,
when:
( ϕ A , ϕ B , ϕ C , ϕ D , ϕ a , ϕ b , ϕ c , ϕ d ) = ( 0 , π 2 , 0 , π 2 , π 4 , π 4 , 3 π 4 , π 4 ) ,
where ϕ A , ϕ B , ϕ C , ϕ D are the phase settings of the four-port operation.

6.2. The Grover Coin Realization

The Grover coin is realized using reversible designs as well. The procedure is identical to the Fourier coin case. The three-dimensional Grover coin takes the form of:
C 3 = 1 3 1 2 2 2 1 2 2 2 1 .
This can be realized using a reversible tritter with phase settings equal to:
( ϕ A , ϕ B , ϕ C , κ z ) = ( 11 π 6 , 11 π 6 , 11 π 6 , 11 π 6 ) ,
or using an unbiased three-port with settings:
( ϕ A , ϕ B , ϕ C ) = ( 3 π 2 , 3 π 2 , 3 π 2 ) .
The four-dimensional Grover coin operator is given by:
C 4 = 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ,
which can be realized with a reversible quarter,
U G r o v e r = U p h a s e o u t U Q u a r t e r U p h a s e i n .
The phase settings for this Grover coin realization is done by:
( ϕ A , ϕ B , ϕ C , ϕ D , κ z 1 , κ z 2 ) = ( 0 , 0 , 0 , 0 , π 8 , π 8 ) .
Similarly for the four ports:
U G r o v e r = U p h a s e o u t U F o u r p o r t s U p h a s e i n ,
( ϕ A , ϕ B , ϕ C , ϕ D ) = ( 3 π 2 , 3 π 2 , 3 π 2 , 3 π 2 ) .

7. Comparison between Directional- and Directionally-Unbiased Devices

Directional devices (the Reck and Clements decomposition model) can produce any unitary matrices U(N). However, when reversibility is introduced in the system (reversible tritters and directionally-unbiased linear-optical multiports), it imposes symmetry or the self-transpose property U i , j = U j , i , where i, j are matrix indices. Hence, reversible designs only produce the subset of symmetric unitary matrices. As an example, Equation (111) is a unitary matrix, but it is not a self-transpose matrix.
U = 0 1 0 0 0 1 1 0 0 .
Nevertheless, this reversible design can produce important coins, such as Grover and Fourier coins, for quantum walks. The properties of each directional- and directionally-unbiased design are compared in Table 1. In this table, directionally-unbiased three-port and four-port operations are denoted as 3-port and 4-port; Reversible tritter and quarters are denoted as Rev Tritter and Rev Quarter; 3-port Reck, 4-port Reck, and 4-port Clements represent the directional three-port Reck model, the directional four-port Reck model, and the directional four-port Clements model, respectively. In the table, we consider several different conditions for comparing different optical devices: the numbers of beam splitters used to form a device, the coherence length requirement to generate the final unitary matrix, and dense unitary matrix generation. The Grover and the Fourier coin generation are considered as indicated in the Conditions column in Table 1. The number of beam splitters varies depending on the design. Directionally-unbiased linear-optical multiports require the fewest beam splitters among all the devices in the table. A directionally-unbiased N-port requires N beam splitters, while other devices with N input and N output ports would require N(N + 1) beam splitters. The multiport device demands a long coherence length because different photons traveling paths with different travel lengths need to be coherently summed. In the reversible designs, Reck decomposition, and the Clements decomposition model, the input photon does not need long coherence lengths because the devices consist of balanced interferometers. All the devices listed in the table can realize both the Grover and the Fourier coins.

8. Summary

Directional optical designs and directionally-unbiased linear-optical designs were reviewed and investigated closely. We studied the use of directionally-unbiased linear-optical designs in quantum walk applications. This directionally-unbiased system allowed us to generate reversible graphs with vertices having multiple edges. This flexibility in graph generation can be useful in quantum walk-based algorithmic speedup as briefly mentioned in the introduction part of Section 5. Previously introduced directionally-unbiased designs [43,57], reversible optical tritter and quarter, can work as scattering centers having three and four edges for quantum walk applications. These designs are advantageous because of the implementation cost reduction by removing directional bias in the system. We focused on Grover and Fourier coin implementations, which are important matrices in quantum information processing. Grover and Fourier matrices can be realized when all the phases are at proper settings. Any unitary matrices can be realized using directional devices; however, this is not the case for the directionally-unbiased design. Directionally-unbiased designs cannot realize non-self-transpose unitary matrices. We have focused on quantum walk and search applications, but similar designs can be applied in other applications. For example, Hamiltonian simulations and topological phase simulations are possible immediate applications of directionally-unbiased designs [58,59,60].

Author Contributions

Conceptualization, S.O., D.S.S. and A.V.S.; software, S.O.; writing—original draft preparation, S.O.; writing—review and editing, D.S.S. and A.V.S.; supervision, D.S.S. and A.V.S.; funding acquisition, A.V.S.

Funding

This research was funded by the National Science Foundation, Emerging Frontiers in Research and Innovation-Advancing Communication Quantum Information Research in Engineering (EFRI-ACQUIRE) grant number ECCS-1640968, the Air Force Office of Scientific Research (AFOSR) grant number FA9550-18-1-0056, and the Northrop Grumman NexGen.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Shor, P.W. Polynomial-time algorithms for prime factorization and discrete logarithms on a quantum computer. SIAM Rev. 1999, 41, 303–332. [Google Scholar] [CrossRef]
  2. Grover, L.K. A fast quantum mechanical algorithm for database search. In Proceedings of the Twenty-Eighth Annual ACM Symposium on Theory of Computing, Philadelphia, PA, USA, 22–24 May 1996; pp. 212–219. [Google Scholar]
  3. Knill, E.; Laflamme, R.; Milburn, G.J. A scheme for efficient quantum computation with linear optics. Nature 2001, 409, 46–52. [Google Scholar] [CrossRef] [PubMed]
  4. Kempe, J. Quantum random walks: An introductory overview. Contemp. Phys. 2003, 44, 307–327. [Google Scholar] [CrossRef] [Green Version]
  5. Venegas-Andraca, S.E. Quantum walks: A comprehensive review. Quant. Inf. Proc. 2012, 11, 1015–1106. [Google Scholar] [CrossRef]
  6. Aharonov, Y.; Davidovich, L.; Zagury, N. Quantum random walks. Phys. Rev. A 1993, 48, 1687–1690. [Google Scholar] [CrossRef] [PubMed]
  7. Portugal, R. Quantum Walks and Search Algorithms; Springer Science & Business Media: New York, NY, USA, 2013. [Google Scholar]
  8. Childs, A.M. Universal computation by quantum walk. Phys. Rev. Lett. 2009, 102, 180501. [Google Scholar] [CrossRef] [PubMed]
  9. Carolan, J.; Harrold, C.; Sparrow, C.; Martín-López, E.; Russell, N.J.; Silverstone, J.W.; Shadbolt, P.J.; Matsuda, N.; Oguma, M.; Itoh, M.; et al. Universal linear optics. Science 2015, 349, 711–716. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Wang, J.; Paesani, S.; Ding, Y.; Santagati, R.; Skrzypczyk, P.; Salavrakos, A.; Tura, J.; Augusiak, R.; Mančinska, L.; Bacco, D.; et al. Multidimensional quantum entanglement with large-scale integrated optics. Science 2018, 360, 285–291. [Google Scholar] [CrossRef] [Green Version]
  11. Reck, M.; Zeilinger, A.; Bernstein, H.J.; Bertani, P. Experimental realization of any discrete unitary operator. Phys. Rev. Lett. 1994, 73, 58–61. [Google Scholar] [CrossRef]
  12. Clements, W.R.; Humphreys, P.C.; Metcalf, B.J.; Kolthammer, W.S.; Walmsley, I.A. Optimal design for universal multiport interferometers. Optica 2016, 3, 1460–1465. [Google Scholar] [CrossRef]
  13. Spring, J.B.; Metcalf, B.J.; Humphreys, P.C.; Kolthammer, W.S.; Jin, X.M.; Barbieri, M.; Datta, A.; Thomas-Peter, N.; Langford, N.K.; Kundys, D.; et al. Boson sampling on a photonic chip. Science 2013, 339, 798–801. [Google Scholar] [CrossRef] [PubMed]
  14. Broome, M.A.; Fedrizzi, A.; Rahimi-Keshari, S.; Dove, J.; Aaronson, S.; Ralph, T.C.; White, A.G. Photonic boson sampling in a tunable circuit. Science 2013, 339, 794–798. [Google Scholar] [CrossRef] [PubMed]
  15. Tillmann, M.; Dakić, B.; Heilmann, R.; Nolte, S.; Szameit, A.; Walther, P. Experimental boson sampling. Nat. Photonics 2013, 7, 540. [Google Scholar] [CrossRef]
  16. Crespi, A.; Osellame, R.; Ramponi, R.; Brod, D.J.; Galvao, E.F.; Spagnolo, N.; Vitelli, C.; Maiorino, E.; Mataloni, P.; Sciarrino, F. Integrated multimode interferometers with arbitrary designs for photonic boson sampling. Nat. Photonics 2013, 7, 545. [Google Scholar] [CrossRef]
  17. Weihs, G.; Reck, M.; Weinfurter, H.; Zeilinger, A. All-fiber three-path Mach–Zehnder interferometer. Opt. Lett. 1996, 21, 302–304. [Google Scholar] [CrossRef]
  18. Peruzzo, A.; Laing, A.; Politi, A.; Rudolph, T.; O’brien, J.L. Multimode quantum interference of photons in multiport integrated devices. Nat. Commun. 2011, 2, 224. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Spagnolo, N.; Vitelli, C.; Aparo, L.; Mataloni, P.; Sciarrino, F.; Crespi, A.; Ramponi, R.; Osellame, R. Three-photon bosonic coalescence in an integrated tritter. Nat. Commun. 2013, 4, 1606. [Google Scholar] [CrossRef] [Green Version]
  20. Meany, T.; Delanty, M.; Gross, S.; Marshall, G.D.; Steel, M.; Withford, M.J. Non-classical interference in integrated 3D multiports. Opt. Express 2012, 20, 26895–26905. [Google Scholar] [CrossRef]
  21. Moore, C.; Russell, A. Quantum walks on the hypercube. In Randomization and Approximation Techniques in Computer Science; Springer: Berlin/Heidelberg, Germany, 2002; pp. 164–178. [Google Scholar]
  22. Krovi, H.; Brun, T.A. Hitting time for quantum walks on the hypercube. Phys. Rev. A 2006, 73, 032341. [Google Scholar] [CrossRef] [Green Version]
  23. Childs, A.M.; Farhi, E.; Gutmann, S. An example of the difference between quantum and classical random walks. Quant. Inf. Proc. 2002, 1, 35–43. [Google Scholar] [CrossRef]
  24. Childs, A.M.; Cleve, R.; Deotto, E. Exponential algorithmic speedup by a quantum walk. In Proceedings of the 35th ACM Symposium on Theory of Computing, San Diego, CA, USA, 9–11 June 2003; pp. 59–68. [Google Scholar]
  25. Ambainis, A. Quantum walk algorithm for element distinctness. SIAM J. Comput. 2007, 37, 210–239. [Google Scholar] [CrossRef]
  26. Magniez, F.; Santha, M.; Szegedy, M. Quantum algorithms for the triangle problem. SIAM J. Comput. 2007, 37, 413–424. [Google Scholar] [CrossRef]
  27. Buhrman, H.; Špalek, R. Quantum verification of matrix products. In Proceedings of the Seventeenth Annual ACM-SIAM Symposium on Discrete Algorithm, Miami, FL, USA, 22–26 January 2006; pp. 880–889. [Google Scholar]
  28. Magniez, F.; Nayak, A. Quantum complexity of testing group commutativity. Algorithmica 2007, 48, 221–232. [Google Scholar] [CrossRef]
  29. Bouwmeester, D.; Marzoli, I.; Karman, G.P.; Schleich, W.; Woerdman, J. Optical galton board. Phys. Rev. A 1999, 61, 013410. [Google Scholar] [CrossRef]
  30. Knight, P.L.; Roldán, E.; Sipe, J.E. Optical cavity implementations of the quantum walk. Opt. Commun. 2003, 227, 147–157. [Google Scholar] [CrossRef] [Green Version]
  31. Knight, P.L.; Roldán, E.; Sipe, J. Quantum walk on the line as an interference phenomenon. Phys. Rev. A 2003, 68, 020301. [Google Scholar] [CrossRef] [Green Version]
  32. Schreiber, A.; Cassemiro, K.N.; Potoček, V.; Gábris, A.; Mosley, P.J.; Andersson, E.; Jex, I.; Silberhorn, C. Photons walking the line: A quantum walk with adjustable coin operations. Phys. Rev. Lett. 2010, 104, 050502. [Google Scholar] [CrossRef] [PubMed]
  33. Pandey, D.; Satapathy, N.; Meena, M.; Ramachandran, H. Quantum walk of light in frequency space and its controlled dephasing. Phys. Rev. A 2011, 84, 042322. [Google Scholar] [CrossRef]
  34. Zhao, Z.; Du, J.; Li, H.; Yang, T.; Chen, Z.B.; Pan, J.W. Implement quantum random walks with linear optics elements. arXiv Preprint, 2002; arXiv:quant-ph/0212149. [Google Scholar]
  35. Broome, M.A.; Fedrizzi, A.; Lanyon, B.P.; Kassal, I.; Aspuru-Guzik, A.; White, A.G. Discrete single-photon quantum walks with tunable decoherence. Phys. Rev. Lett. 2010, 104, 153602. [Google Scholar] [CrossRef]
  36. Goyal, S.K.; Roux, F.S.; Forbes, A.; Konrad, T. Implementing quantum walks using orbital angular momentum of classical light. Phys. Rev. Lett. 2013, 110, 263602. [Google Scholar] [CrossRef]
  37. Zhang, P.; Ren, X.F.; Zou, X.B.; Liu, B.H.; Huang, Y.F.; Guo, G.C. Demonstration of one-dimensional quantum random walks using orbital angular momentum of photons. Phys. Rev. A 2007, 75, 052310. [Google Scholar] [CrossRef]
  38. Cardano, F.; Massa, F.; Qassim, H.; Karimi, E.; Slussarenko, S.; Paparo, D.; de Lisio, C.; Sciarrino, F.; Santamato, E.; Boyd, R.W.; et al. Quantum walks and wavepacket dynamics on a lattice with twisted photons. Sci. Adv. 2015, 1, e1500087. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Tang, H.; Di Franco, C.; Shi, Z.Y.; He, T.S.; Feng, Z.; Gao, J.; Sun, K.; Li, Z.M.; Jiao, Z.Q.; Wang, T.Y.; et al. Experimental quantum fast hitting on hexagonal graphs. Nat. Photonics 2018, 12, 754–758. [Google Scholar] [CrossRef] [Green Version]
  40. Schreiber, A.; Gábris, A.; Rohde, P.P.; Laiho, K.; Štefaňák, M.; Potoček, V.; Hamilton, C.; Jex, I.; Silberhorn, C. A 2D quantum walk simulation of two-particle dynamics. Science 2012, 336, 55–58. [Google Scholar] [CrossRef] [PubMed]
  41. Feldman, E.; Hillery, M. Scattering theory and discrete-time quantum walks. Phys. Lett. A 2004, 324, 277–281. [Google Scholar] [CrossRef] [Green Version]
  42. Feldman, E.; Hillery, M. Modifying quantum walks: A scattering theory approach. J. Phys. A 2007, 40, 11343. [Google Scholar] [CrossRef]
  43. Simon, D.S.; Fitzpatrick, C.A.; Sergienko, A.V. Group transformations and entangled-state quantum gates with directionally unbiased linear-optical multiports. Phys. Rev. A 2016, 93, 043845. [Google Scholar] [CrossRef] [Green Version]
  44. Murnaghan, F.D. The Unitary and Rotation Groups; Lectures on Applied Mathematics; Spartan Books: Washington, DC, USA, 1962; Volume 3. [Google Scholar]
  45. Loudon, R. The Quantum Theory of Light; Oxford University Press: Oxford, UK, 2000. [Google Scholar]
  46. Saleh, B.E.; Teich, M.C.; Saleh, B.E. Fundamentals of Photonics; Wiley: New York, NY, USA, 1991; Volume 22. [Google Scholar]
  47. Bromberg, Y.; Lahini, Y.; Morandotti, R.; Silberberg, Y. Quantum and classical correlations in waveguide lattices. Phys. Rev. Lett. 2009, 102, 253904. [Google Scholar] [CrossRef]
  48. Sansoni, L.; Sciarrino, F.; Vallone, G.; Mataloni, P.; Crespi, A.; Ramponi, R.; Osellame, R. Two-particle bosonic-fermionic quantum walk via integrated photonics. Phys. Rev. Lett. 2012, 108, 010502. [Google Scholar] [CrossRef]
  49. Kowalevicz, A.; Sharma, V.; Ippen, E.; Fujimoto, J.G.; Minoshima, K. Three-dimensional photonic devices fabricated in glass by use of a femtosecond laser oscillator. Opt. Lett. 2005, 30, 1060–1062. [Google Scholar] [CrossRef]
  50. Suzuki, K.; Sharma, V.; Fujimoto, J.G.; Ippen, E.P.; Nasu, Y. Characterization of symmetric [3 × 3] directional couplers fabricated by direct writing with a femtosecond laser oscillator. Opt. Express 2006, 14, 2335–2343. [Google Scholar] [CrossRef] [PubMed]
  51. De Guise, H.; Di Matteo, O.; Sánchez-Soto, L.L. Simple factorization of unitary transformations. Phys. Rev. A 2018, 97, 022328. [Google Scholar] [CrossRef]
  52. Metcalf, B.J.; Thomas-Peter, N.; Spring, J.B.; Kundys, D.; Broome, M.A.; Humphreys, P.C.; Jin, X.M.; Barbieri, M.; Kolthammer, W.S.; Gates, J.C.; et al. Multiphoton quantum interference in a multiport integrated photonic device. Nat. Commun. 2013, 4, 1356. [Google Scholar] [CrossRef] [PubMed]
  53. Spagnolo, N.; Vitelli, C.; Sansoni, L.; Maiorino, E.; Mataloni, P.; Sciarrino, F.; Brod, D.J.; Galvão, E.F.; Crespi, A.; Ramponi, R.; et al. General rules for bosonic bunching in multimode interferometers. Phys. Rev. Lett. 2013, 111, 130503. [Google Scholar] [CrossRef] [PubMed]
  54. Spagnolo, N.; Aparo, L.; Vitelli, C.; Crespi, A.; Ramponi, R.; Osellame, R.; Mataloni, P.; Sciarrino, F. Quantum interferometry with three-dimensional geometry. Sci. Rep. 2012, 2, 862. [Google Scholar] [CrossRef]
  55. Zhang, J.G.; Sharma, A.; Ni, Y.D.; Li, Z. Investigation into network architecture and modulation scheme for MIL-STD-1773 optical fiber data buses. Aircr. Eng. Aerosp. Technol. 2000, 72, 126–137. [Google Scholar] [CrossRef]
  56. Saleh, A.; Kogelnik, H. Reflective single-mode fiber-optic passive star couplers. J. Light. Technol. 1988, 6, 392–398. [Google Scholar] [CrossRef]
  57. Osawa, S.; Simon, D.S.; Sergienko, A.V. Experimental demonstration of a directionally-unbiased linear-optical multiport. Opt. Express 2018, 26, 27201–27211. [Google Scholar] [CrossRef]
  58. Simon, D.S.; Fitzpatrick, C.A.; Osawa, S.; Sergienko, A.V. Quantum simulation of discrete-time Hamiltonians using directionally unbiased linear optical multiports. Phys. Rev. A 2017, 95, 042109. [Google Scholar] [CrossRef] [Green Version]
  59. Simon, D.S.; Fitzpatrick, C.A.; Osawa, S.; Sergienko, A.V. Quantum simulation of topologically protected states using directionally unbiased linear-optical multiports. Phys. Rev. A 2017, 96, 013858. [Google Scholar] [CrossRef] [Green Version]
  60. Simon, D.S.; Osawa, S.; Sergienko, A.V. Joint entanglement of topology and polarization enables error-protected quantum registers. New J. Phys. 2018, 20, 093032. [Google Scholar] [CrossRef]
  61. Vance, R.; Barrow, R. General linear differential interferometers. J. Opt. Soc. Am. A 1995, 12, 346–353. [Google Scholar] [CrossRef]
  62. Schwelb, O. Generalized analysis for a class of linear interferometric networks. I. Analysis. IEEE Trans. Microw. Theory Tech. 1998, 46, 1399–1408. [Google Scholar] [CrossRef]
  63. Motwani, R.; Raghavan, P. Randomized algorithms. ACM Comput. Surv. 1996, 28, 33–37. [Google Scholar] [CrossRef]
  64. Francisco, D.; Iemmi, C.; Paz, J.; Ledesma, S. Simulating a quantum walk with classical optics. Phys. Rev. A 2006, 74, 052327. [Google Scholar] [CrossRef]
  65. Peruzzo, A.; Lobino, M.; Matthews, J.C.; Matsuda, N.; Politi, A.; Poulios, K.; Zhou, X.Q.; Lahini, Y.; Ismail, N.; Wörhoff, K.; et al. Quantum walks of correlated photons. Science 2010, 329, 1500–1503. [Google Scholar] [CrossRef]
  66. Tang, H.; Lin, X.F.; Feng, Z.; Chen, J.Y.; Gao, J.; Sun, K.; Wang, C.Y.; Lai, P.C.; Xu, X.Y.; Wang, Y.; et al. Experimental two-dimensional quantum walk on a photonic chip. Sci. Adv. 2018, 4, eaat3174. [Google Scholar] [CrossRef] [Green Version]
  67. Venancio, B.; Andrade, F.; da Luz, M. Unveiling and exemplifying the unitary equivalence of discrete time quantum walk models. J. Phys. A 2013, 46, 165302. [Google Scholar] [CrossRef] [Green Version]
  68. Hillery, M.; Bergou, J.; Feldman, E. Quantum walks based on an interferometric analogy. Phys. Rev. A 2003, 68, 032314. [Google Scholar] [CrossRef] [Green Version]
  69. Benioff, P. Space searches with a quantum robot. arXiv Preprint, 2000; arXiv:quant-ph/0003006. [Google Scholar]
  70. Aaronson, S.; Ambainis, A. Quantum search of spatial regions. In Proceedings of the 44th Annual IEEE Symposium on Foundations of Computer Science, Cambridge, MA, USA, 11–14 October 2003; pp. 200–209. [Google Scholar]
  71. Childs, A.M.; Goldstone, J. Spatial search by quantum walk. Phys. Rev. A 2004, 70, 022314. [Google Scholar] [CrossRef]
  72. Tulsi, A. Faster quantum-walk algorithm for the two-dimensional spatial search. Phys. Rev. A 2008, 78, 012310. [Google Scholar] [CrossRef] [Green Version]
  73. Abal, G.; Donangelo, R.; Marquezino, F.L.; Portugal, R. Spatial search on a honeycomb network. Math. Struct. Comp. Sci. 2010, 20, 999–1009. [Google Scholar] [CrossRef] [Green Version]
  74. Abal, G.; Donangelo, R.; Forets, M.; Portugal, R. Spatial quantum search in a triangular network. Math. Struct. Comp. Sci. 2012, 22, 521–531. [Google Scholar] [CrossRef] [Green Version]
  75. Kempe, J. Discrete quantum walks hit exponentially faster. Probab. Theory Relat. Fields 2005, 133, 215–235. [Google Scholar] [CrossRef]
  76. Tregenna, B.; Flanagan, W.; Maile, R.; Kendon, V. Controlling discrete quantum walks: Coins and initial states. New J. Phys. 2003, 5, 83. [Google Scholar] [CrossRef]
Figure 1. Description of a beam splitter. A beam splitter is a device with two input and two output ports. A photon can enter the port E 1 and/or E 2 and will leave the port E 3 and/or E 4 .
Figure 1. Description of a beam splitter. A beam splitter is a device with two input and two output ports. A photon can enter the port E 1 and/or E 2 and will leave the port E 3 and/or E 4 .
Entropy 21 00853 g001
Figure 2. (a) A 2 × 2 directional coupler. k 1 and k 2 represent input and output spatial modes. A photon can enter either k 1 or k 2 , and a coupler transforms the input state. Two waveguides are closely located to allow evanescent coupling at the cross-section in the figure. (b) The cross-section area in the directional coupler illustrated in (a). The coupling strength between two waveguides is κ , and the propagation constant in each waveguide is β .
Figure 2. (a) A 2 × 2 directional coupler. k 1 and k 2 represent input and output spatial modes. A photon can enter either k 1 or k 2 , and a coupler transforms the input state. Two waveguides are closely located to allow evanescent coupling at the cross-section in the figure. (b) The cross-section area in the directional coupler illustrated in (a). The coupling strength between two waveguides is κ , and the propagation constant in each waveguide is β .
Entropy 21 00853 g002
Figure 3. The Mach–Zehnder interferometer. A photon can enter either the port E 1 or E 2 and will be transformed by the first beam splitter. The photon can leave either through a superposition of E 3 and E 4 and be transformed by the second beam splitter. Finally, the photon leaves the device either through E 5 and/or E 6 modes.
Figure 3. The Mach–Zehnder interferometer. A photon can enter either the port E 1 or E 2 and will be transformed by the first beam splitter. The photon can leave either through a superposition of E 3 and E 4 and be transformed by the second beam splitter. Finally, the photon leaves the device either through E 5 and/or E 6 modes.
Entropy 21 00853 g003
Figure 4. The Michelson interferometer. A photon can enter either the port E 1 and/or E 2 . The photon interacts with a beam splitter, mirror, and again with a beamsplitter. The photon leaves through either E 1 and/or E 2 , which are the same as the input ports. The phase can be controlled by translating mirrors in the system.
Figure 4. The Michelson interferometer. A photon can enter either the port E 1 and/or E 2 . The photon interacts with a beam splitter, mirror, and again with a beamsplitter. The photon leaves through either E 1 and/or E 2 , which are the same as the input ports. The phase can be controlled by translating mirrors in the system.
Entropy 21 00853 g004
Figure 5. (a) The 3 × 3 Reck model realization using three beam splitters. A photon can enter either Port 1, 2, or 3 and the photon leave either Port 1’, 2’, and/or 3’. A beam splitter in the system can be substituted by a Mach–Zehnder interferometer if one wants to give amplitude tuning at each beam splitter encounter. The beam splitter requirement will be increased to six when amplitude tuning by interferometers is imposed. (b) The information flow decomposition of the bulk 3 × 3 setup using a set of 2 × 2 unitary matrices. This is equivalent to the physical setup in (a).
Figure 5. (a) The 3 × 3 Reck model realization using three beam splitters. A photon can enter either Port 1, 2, or 3 and the photon leave either Port 1’, 2’, and/or 3’. A beam splitter in the system can be substituted by a Mach–Zehnder interferometer if one wants to give amplitude tuning at each beam splitter encounter. The beam splitter requirement will be increased to six when amplitude tuning by interferometers is imposed. (b) The information flow decomposition of the bulk 3 × 3 setup using a set of 2 × 2 unitary matrices. This is equivalent to the physical setup in (a).
Entropy 21 00853 g005
Figure 6. (a) The 4 × 4 Reck model using six beam splitters (12 with tunability). A photon can enter either one of four input ports and can leave through any of the four output ports. (b) The information flow representation in the case of decomposing a bulk 4 × 4 setup using 2 × 2 unitary matrices illustrated in (a).
Figure 6. (a) The 4 × 4 Reck model using six beam splitters (12 with tunability). A photon can enter either one of four input ports and can leave through any of the four output ports. (b) The information flow representation in the case of decomposing a bulk 4 × 4 setup using 2 × 2 unitary matrices illustrated in (a).
Entropy 21 00853 g006
Figure 7. (a) 4 × 4 mesh of the Reck design. Input photons flow from the left to the right. Each line cross-section consists of the two-dimensional two-mode mixer illustrated in Figure 8. This is identical to the setup in Figure 6a. (b) 4 × 4 mesh of symmetric Clements design.
Figure 7. (a) 4 × 4 mesh of the Reck design. Input photons flow from the left to the right. Each line cross-section consists of the two-dimensional two-mode mixer illustrated in Figure 8. This is identical to the setup in Figure 6a. (b) 4 × 4 mesh of symmetric Clements design.
Entropy 21 00853 g007
Figure 8. An integrated coupler requires tunability to change a transfer matrix of the system. The tunability is acquired through an interferometer with phase shifters. θ is an external phase shift, and ϕ is an internal phase shift.
Figure 8. An integrated coupler requires tunability to change a transfer matrix of the system. The tunability is acquired through an interferometer with phase shifters. θ is an external phase shift, and ϕ is an internal phase shift.
Entropy 21 00853 g008
Figure 9. (a) Integrated optical tritter. k i { i = 1 , 2 , 3 } represent input and output spatial modes. (b) Coupling region of the integrated tritter. β is a waveguide propagation coefficient, and κ is a coupling coefficient between the couplers.
Figure 9. (a) Integrated optical tritter. k i { i = 1 , 2 , 3 } represent input and output spatial modes. (b) Coupling region of the integrated tritter. β is a waveguide propagation coefficient, and κ is a coupling coefficient between the couplers.
Entropy 21 00853 g009
Figure 10. (a) Integrated quarter. k i { i = 1 , 2 , 3 } represent input and output spatial modes. (b) Cross-section area of the integrated quarter. A waveguide propagation constant is β , and coupling coefficients are κ 1 for the neighboring couplers and κ 2 for the diagonal couplers.
Figure 10. (a) Integrated quarter. k i { i = 1 , 2 , 3 } represent input and output spatial modes. (b) Cross-section area of the integrated quarter. A waveguide propagation constant is β , and coupling coefficients are κ 1 for the neighboring couplers and κ 2 for the diagonal couplers.
Entropy 21 00853 g010
Figure 11. Directionally-unbiased linear-optical three-port operation. Input ports can be used as outputs in this system. For example, if a photon is inserted at Port A, then the photon would leave Ports A, B, and C. The beam splitter (BS) splits the incoming photon into two outgoing directions. Mirror units are necessary to reverse the propagation direction so that the directional system becomes directionally-unbiased. Mirror units consist of a phase shifter P and mirror M, as illustrated in the figure.
Figure 11. Directionally-unbiased linear-optical three-port operation. Input ports can be used as outputs in this system. For example, if a photon is inserted at Port A, then the photon would leave Ports A, B, and C. The beam splitter (BS) splits the incoming photon into two outgoing directions. Mirror units are necessary to reverse the propagation direction so that the directional system becomes directionally-unbiased. Mirror units consist of a phase shifter P and mirror M, as illustrated in the figure.
Entropy 21 00853 g011
Figure 12. Directionally-unbiased linear-optical four-port operation. It has the same configuration as Figure 11 with one extra input/output port.
Figure 12. Directionally-unbiased linear-optical four-port operation. It has the same configuration as Figure 11 with one extra input/output port.
Entropy 21 00853 g012
Figure 13. (a) Reversible integrated tritter. The photon propagation direction is reversed by placing mirrors at the end of the coupling region. (b) Reversible integrated quarter. It has the same configuration as (a) with an extra input and output port.
Figure 13. (a) Reversible integrated tritter. The photon propagation direction is reversed by placing mirrors at the end of the coupling region. (b) Reversible integrated quarter. It has the same configuration as (a) with an extra input and output port.
Entropy 21 00853 g013
Figure 14. (a) Directionally-insensitive description for a scattering center. A photon can enter Port A and then leave either Port A or B. (b) Directionally-sensitive description for a scattering center. If the photon is initially in the state A R , then the photon will have an amplitude in the A L direction and the B R direction. Photon in states A R and A L do not interact, so they need to be distinguished when a graph is formed based on scattering centers.
Figure 14. (a) Directionally-insensitive description for a scattering center. A photon can enter Port A and then leave either Port A or B. (b) Directionally-sensitive description for a scattering center. If the photon is initially in the state A R , then the photon will have an amplitude in the A L direction and the B R direction. Photon in states A R and A L do not interact, so they need to be distinguished when a graph is formed based on scattering centers.
Entropy 21 00853 g014
Figure 15. (a) Rectangular lattice. This structure has four edges for every vertex. (b) Hexagonal lattice. This structure has three edges for every vertex.
Figure 15. (a) Rectangular lattice. This structure has four edges for every vertex. (b) Hexagonal lattice. This structure has three edges for every vertex.
Entropy 21 00853 g015
Figure 16. (a) Glued tree. Two trees are glued in a random manner in the middle part of the graph. A photon’s hitting time from one red circle to the other red circle is shorter than the classical walk on this graph. (b) Four-dimensional hypercube. Every vertex has four edges. A photon’s hitting time from one red circle to the other red circle is shorter than the classical walk on this graph.
Figure 16. (a) Glued tree. Two trees are glued in a random manner in the middle part of the graph. A photon’s hitting time from one red circle to the other red circle is shorter than the classical walk on this graph. (b) Four-dimensional hypercube. Every vertex has four edges. A photon’s hitting time from one red circle to the other red circle is shorter than the classical walk on this graph.
Entropy 21 00853 g016
Figure 17. Reversible tritter with phase shifters. The yellow squares are phase shifters, and the green squares are mirrors.
Figure 17. Reversible tritter with phase shifters. The yellow squares are phase shifters, and the green squares are mirrors.
Entropy 21 00853 g017
Figure 18. Reversible quarter with phase shifters. The yellow squares are phase shifters, and the green squares are mirrors.
Figure 18. Reversible quarter with phase shifters. The yellow squares are phase shifters, and the green squares are mirrors.
Entropy 21 00853 g018
Table 1. Comparison between reversible and non-reversible designs. In the table, 3-port and 4-port represent directionally-unbiased linear-optical multiports; Rev Tritter and Rev Quarter represent reversible tritter and reversible quarter; 3-Reck, 4-Reck and 4-Clements represent the three-port and four-port Reck decomposition models and the four-port Clements decomposition model.
Table 1. Comparison between reversible and non-reversible designs. In the table, 3-port and 4-port represent directionally-unbiased linear-optical multiports; Rev Tritter and Rev Quarter represent reversible tritter and reversible quarter; 3-Reck, 4-Reck and 4-Clements represent the three-port and four-port Reck decomposition models and the four-port Clements decomposition model.
Conditions3-port4-portRev TritterRev Quarter3-Reck4-Reck4-Clements
# of Beam splitters34--122020
Coherence LengthLongLongShortShortShortShortShort
U Generation
Grover Coin
Fourier Coin

Share and Cite

MDPI and ACS Style

Osawa, S.; Simon, D.S.; Sergienko, A.V. Directionally-Unbiased Unitary Optical Devices in Discrete-Time Quantum Walks. Entropy 2019, 21, 853. https://doi.org/10.3390/e21090853

AMA Style

Osawa S, Simon DS, Sergienko AV. Directionally-Unbiased Unitary Optical Devices in Discrete-Time Quantum Walks. Entropy. 2019; 21(9):853. https://doi.org/10.3390/e21090853

Chicago/Turabian Style

Osawa, Shuto, David S. Simon, and Alexander V. Sergienko. 2019. "Directionally-Unbiased Unitary Optical Devices in Discrete-Time Quantum Walks" Entropy 21, no. 9: 853. https://doi.org/10.3390/e21090853

APA Style

Osawa, S., Simon, D. S., & Sergienko, A. V. (2019). Directionally-Unbiased Unitary Optical Devices in Discrete-Time Quantum Walks. Entropy, 21(9), 853. https://doi.org/10.3390/e21090853

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop