Patent ID: 12253719

DETAILED DESCRIPTION

The following description of exemplary embodiments of the invention is presented for the purpose of illustration and description. It is not intended to be exhaustive or to limit the claimed invention to the precise form described, and persons skilled in the art will appreciate that many modifications and variations are possible. The embodiments have been chosen and described in order to best explain the principles of the invention and its practical applications to thereby enable others skilled in the art to best make and use the invention in various embodiments and with various modifications as are suited to the particular use contemplated.

GMZIs and MZIs can be used as the building blocks for many photonic switching networks. Certain embodiments described herein relate to a GMZI that uses Mach-Zehnder lattice (MZL) subcomponents to suppress the effect of phase-shifter loss in all switch configurations. (MZLs are interferometric devices which can be used for optical switching with digital-like behavior, see e.g. J. Van Campenhout, et al. Optics Express 17, 23793 (2009); J. Van Campenhout, et al. Optics Express 19, 11568 (2011)). Examples of GMZIs and optical circuits incorporating GMZIs are described below with reference toFIGS.11A-18B.

Certain embodiments described herein provide a “Hadamard-type” GMZI. This type of GMZI is an N-mode device which can enact N switching configurations (where N is any power of two). Each of these configurations performs a permutation operation on the light at the input ports, and it is possible to route light from any one input port to any chosen target output port. The configurations are selected using one layer of active phase shifter devices which are set to values 0 or π. Thus, a Hadamard-type GMZI can be used as an N×N optical switch.

For purposes of the present description, the phase-shifter loss model is taken to be transmittance T(φ)=T(π)|φ|πwhere φ is a controllable phase shift between 0 and 2π. This loss model is relevant to phase-shifter devices based on carrier depletion and injection, as well as devices which suffer from unwanted polarization rotation. The phase shifter error model is assumed to be Gaussian-distributed random phase error with standard deviation scaling ∝φ.

Some embodiments can suppress loss and phase error, with additional interferometric stages achieving lower loss and more “digital: behavior for the N×N switch. Switches of this kind can be used in a variety of applications in classical and quantum photonics which depend on high-speed optical switching where it is important to minimize switch losses.

1. Phase Shifters Using Mach Zehnder Lattice

FIG.1Ashows a simplified circuit diagram of a Hadamard-type GMZI100according to some embodiments. GMZI100has a set of N input ports102and a set of N output ports104, where N=2afor integer a. (In other words, N is a power of 2.) Input ports102are coupled to inputs of a first Hadamard-type passive interferometer110, and output ports105are coupled to outputs of a second Hadamard-type passive interferometer112. (Example implementations of Hadamard-type passive interferometers are described below.) Each output port103of first Hadamard-type interferometer110is coupled to a corresponding input port104of Hadamard-type passive interferometer112by a phase shifter106. A control logic circuit (not shown), which can be a classical digital logic circuit, can generate control signals for phase shifters106. By selecting an appropriate set of control signals for phase shifters106, GMZI100can be operated to output, at output ports10adesired permutation of a group of photons received on input ports102at output ports105.

According to some embodiments, each phase shifter106can be implemented using a coupled pair of Mach Zehnder lattices (MZL).FIG.1Bshows a simplified circuit diagram of an MZL-based phase shifter circuit106according to some embodiments. Phase shifter106includes two input ports126and128and two output ports130and132. Input port126can be coupled to an output port of Hadamard-type GMZI110ofFIG.1A, and output port130can be coupled to an input port of Hadamard-type GMZI112ofFIG.1A. Input ports128and output ports132can be coupled to vacuum. (As used herein, a “port” can be a section of waveguide to which other waveguides can be coupled.)

Phase shifter circuit106can include a first (1×2) MZL122and a second (2×1) MZL124, each of which includes a pair of waveguides defining first and second optical paths. MZL122includes waveguides160-1and160-2, and MZL124includes waveguides160-3and160-4. Each of MZL122,124can have a number (S) of stages140. Each stage140can include an active phase shifter142disposed along one of the waveguides160-i. Active phase shifters142can be implemented using a variety of techniques, including electro-optic materials for which the refractive index changes as a function of an applied electric field. (Additional examples of phase shifter implementations are described below.)

Each active phase shifter142can apply a selectable phase shift in a range from 0 to π (in radians). Each stage140can also include a dummy phase shifter144on the other optical path. Dummy phase shifter144can apply a constant (e.g., zero) phase shift. Beam splitter regions (also referred to as “directional coupler regions”)150can be included on the optical paths between the pairs of phase shifters142,144, coupling successive stages together. Each directional coupler150can have a transmittance (Ti, also referred to as “transmissivity,” and the transmittance for different directional couplers can be equal or different. For instance, transmittance can be tapered such that it is higher for directional couplers150at one end of MZL122and lower at the other end of MZL122. Such tapering can be in either direction (increasing transmissivity at successively-encountered directional couplers or decreasing transmissivity at successively-encountered directional couplers). In some embodiments, directional coupler regions150can be formed by arranging lengths of waveguides160-3and160-4(or waveguides160-1and160-2) in proximity to each other to provide a desired transmittance; examples are described below.

MZL122and MZL124are coupled together by a pair of phase shifters154,156. Phase shifter154couples waveguide160-2to waveguide160-3, and phase shifter156couples waveguide160-1to160-4. Phase shifters154,156can apply fixed phase shifts, which can be different from each other. For example phase shifter154can apply a phase shift of π while phase shifter156applies a phase shift of 0; more generally the phase shifts applied by phase shifters154,156can differ by π radians. As shown, the coupling can be such that waveguide160-1in MZL122, which passes through active phase shifters142-1through142-S is coupled to waveguide160-4in MZL124, which passes through dummy phase shifters144-(S+1) through144-2S. Conversely, waveguide160-2in MZL122, which passes through dummy phase shifters144-1through144-S is coupled to waveguide160-3in MZL124, which passes through active phase shifters142-(S+1) through142-2S. In this arrangement, a photon entering port126(or port128) that reflects at each directional coupler would pass through S active phase shifters142and S dummy phase shifters144before exiting port132(or port134).

For operation as a phase shifter, one input port126(or128) of phase shifter circuit106can be coupled to receive a photon from an upstream component or device (e.g., from one of the output ports of Hadamard-type passive interferometer110shown inFIG.1A) while the other input port128(or126) is coupled to vacuum Similarly, one output port130(or132) of phase shifter circuit106can be coupled to a downstream component or device (e.g., to one of the input ports of Hadamard-type passive interferometer112shown above). Coupling to vacuum can be achieved, e.g., by terminating a waveguide at a dead end; the termination can include photon absorbing materials to avoid contamination with unwanted light. Each active phase shifter142-1through142-S in MZL122can be set to apply a phase shift of p while each active phase shifter142-1through142-S in MZL122can be set to apply a phase shift φπ. In some embodiments, φπand φ0are chosen such that |φπ−φ0|=π. In this configuration, MZL122can act as a 1×2 switch while MZL124acts as a 2×1 switch. For a photon entering at port126, one of phase shifters154,156can be selected (by controlling the phase shift φπof active phase shifters142-1through142-S), thereby achieving a phase swing of π radians, and active phase shifters142-(S+1) through142-2S can be used to propagate the photon from the selected one of phase shifters154,156to output port130(by controlling the phase shift φ0of active phase shifters142-(S+1) through142-2S).

Phase shifter circuit106can achieve low loss phase shifting of a photon. For instance, when all the phase shifters142are off, the light crosses from one mode into the other due to the combined transmittance of the directional couplers; however when the phase shifters in one arm are set to a suitable value φ the light is confined to the opposite arm. (In general φ can take a range of values which varies according to the precise design.) Because the loss of the phase shifters is assumed to increase with p, the output loss is suppressed, and increasing the number S of MZL stages can strengthen the effect. In some embodiments, tapering the transmittance values for the directional couplers can also enable some additional optimization.

FIG.2shows a modified MZI circuit200that provides a point of comparison for the performance of phase shifter circuit106ofFIG.1B. MZI circuit200is configured to provide a phase shift between a fixed input port202and a fixed output port204. In some GMZI implementations, MZI circuit200can be used instead of a single phase-shifter in order to remove setting-dependence for loss and phase error. Setting phase shifters206,208to (φ, 0) or (0, φ) achieves a phase swing of π (for any value of φ). The transmittance is sin2(φ/2), and setting φ to π maximizes extinction ratio.

FIGS.3A and3Bshow graphs310,320illustrating performance of modified MZI circuit200ofFIG.2as determined from modeling of circuit behavior. Graph310shows output loss as a function of φ, and graph320shows phase variability as a function of φ for different percentage values of phase variation (Δφ). It is assumed that phase-shifter loss Lπ=100 mdB. As shown inFIG.3A, the settings-averaged output loss is ˜Lπ/2 when the phase variability is small (i.e., around 1% of π), but increases significantly with increasing phase-shifter variability. The output phase variability can be shown to be half the value for the component phase shifter which is switched on.

FIG.4shows a graph400illustrating performance of phase shifter circuit106ofFIG.1Baccording to some embodiments, as determined from modeling of circuit behavior. In this example, each MZL122,124has S=3 stages, which corresponds to six interferometric stages across MZLs122and124. Graph400shows output loss as a function of 6φ. As shown, output losses in this example are 27 mdB, 36 mdB, and 65 mdB for 1%, 5% and 10% phase-shifter imprecision respectively. These values are approximately halved compared to modified MZI circuit200(which achieves 51 mdB, 77 mdB, and 157 mdB respectively).

In addition, the spread, or variability, of the output phase shift of phase shifter106in this example is improved by about a factor of 3 (for all three curves) compared to modified MZI circuit200.FIG.5shows a graph500of output phase error (100×Δ(output phase)/π) as a function of 6φ according to some embodiments.

For a large GMZI, the output loss due to reduced extinction ratio arising specifically from the active phase-shifters can be approximately −10 log10[1−(pc×π)2] where pc is is the percentage error on the phase. Hence, the factor of 3 improvement in phase-shifter variability translates into a factor of 10 improvement in loss associated with the extinction ratio. For instance, in some embodiments a 10% error for the component phase shifters can give a110mdB additional loss for the output of the GMZI using the MZI scheme but only a 12 mdB loss using an MZL-based scheme of the kind shown inFIG.1.

It is noted that phase shifter106has an increased total length of phase shifters across all interferometric stages as compared to modified MZI circuit200. For one implementation of phase shifter106with S=3, the total length of phase shifters is almost 3π. The increased total length may introduce additional phase-shifter loss that scales with the total length.

2. Mach Zehnder Interferometer Implementations

The foregoing description makes references to Mach Zehnder interferometers (MZI). In some embodiments, an MZI can be used to couple photons propagating in the same time bin in different waveguides. This section provides additional examples of MZI behavior and implementations, in part to illustrate techniques that can be used to construct GMZIs and/or phase shifters such as phase shifter106. For purposes of the present description, a photon (or absence thereof) in a waveguide in a given time bin is referred to as a “mode”; an “occupied” mode has a photon, while an “unoccupied” or “vacuum” mode does not. In some embodiments, mode occupancy can be used to encode information, including quantum information (e.g., qubits or qudits); however, the present disclosure is not limited to any particular use-case.

FIG.6Ashows a schematic diagram610(also referred to as a circuit diagram or circuit notation) for coupling of two modes. The modes are drawn as horizontal lines612,614, and the mode coupler616is indicated by a vertical line that is terminated with nodes (solid dots) to identify the modes being coupled. In the more specific language of linear quantum optics, the mode coupler616shown inFIG.6Arepresents a 50/50 beam splitter that implements a transfer matrix:

T=12⁢(1ii1),(1)
where T defines the linear map for the photon creation operators on two modes. (In certain contexts, transfer matrix T can be understood as implementing a first-order imaginary Hadamard transform.) By convention the first column of the transfer matrix corresponds to creation operators on the top mode (referred to herein as mode1, labeled as horizontal line612), and the second column corresponds to creation operators on the second mode (referred to herein as mode6, labeled as horizontal line614), and so on if the system includes more than two modes. More explicitly, the mapping can be written as:

(a1†a2†)input↦12⁢(1-i-i1)⁢(a1†a2†)output,(2)
where subscripts on the creation operators indicate the mode that is operated on, the subscripts input and output identify the form of the creation operators before and after the beam splitter, respectively and where:
ai|ni,nj>=√{square root over (ni)}|ni−1,nj>
aj|ni,nj>=√{square root over (nj)}|ni,nj−1>
aj†ni,nj)=√{square root over (nj+1)}|ni,nj+1>  (3)
For example, the application of the mode coupler shown inFIG.6Aleads to the following mappings:

a1input†↦12⁢(a1output†-i⁢a2output†)⁢a2input†↦12⁢(-i⁢a1output†+a2output†)(4)
Thus, the action of the mode coupler described by Eq. (1) is to take the input states |10>, |01>, and |11> to

❘"\[LeftBracketingBar]"10〉↦❘"\[LeftBracketingBar]"10〉-i❘"\[RightBracketingBar]"⁢01〉2(5)❘"\[LeftBracketingBar]"01〉↦-i⁢❘"\[LeftBracketingBar]"10〉+❘"\[RightBracketingBar]"⁢01〉2❘"\[LeftBracketingBar]"11〉↦-i2⁢(❘"\[LeftBracketingBar]"20〉+❘"\[RightBracketingBar]"⁢02〉)

FIG.6Bshows a physical implementation of a mode coupling that implements the transfer matrix T of Eq. (1) for two photonic modes in accordance with some embodiments. In this example, the mode coupling is implemented using a waveguide beam splitter600, also sometimes referred to as a directional coupler or mode coupler. Waveguide beam splitter600can be realized by bringing two waveguides602,604into close enough proximity that the evanescent field of one waveguide can couple into the other. By adjusting the separation d between waveguides602,604and/or the length l of the coupling region, different couplings between modes can be obtained. In this manner, a waveguide beam splitter600can be configured to have a desired transmissivity. For example, the beam splitter can be engineered to have a transmissivity equal to 0.5 (i.e., a 50/50 beam splitter for implementing the specific form of the transfer matrix T introduced above). If other transfer matrices are desired, the reflectivity (or the transmissivity) can be engineered to be greater than 0.6, greater than 0.7, greater than 0.8, or greater than 0.9 without departing from the scope of the present disclosure.

In addition to mode coupling, some unitary transforms may involve phase shifts applied to one or more modes. In some photonic implementations, variable phase-shifters can be implemented in integrated circuits, providing control over the relative phases of the state of a photon spread over multiple modes. Examples of transfer matrices that define such a phase shifts are given by (for applying a +i and −i phase shift to the second mode, respectively):

s=(100i)⁢s†=(100-i)(6)
For silica-on-silicon materials some embodiments implement variable phase-shifters using thermo-optical switches. The thermo-optical switches use resistive elements fabricated on the surface of the chip, that via the thermo-optical effect can provide a change of the refractive index n by raising the temperature of the waveguide by an amount of the order of 10−5K. One of skill in the art with access to the present disclosure will understand that any effect that changes the refractive index of a portion of the waveguide can be used to generate a variable, electrically tunable, phase shift. For example, some embodiments use beam splitters based on any material that supports an electro-optic effect, so-called χ2and χ3 materials such as lithium niobite, BBO, KTP, and the like and even doped semiconductors such as silicon, germanium, and the like.

Beam-splitters with variable transmissivity and arbitrary phase relationships between output modes can also be achieved by combining directional couplers and variable phase-shifters in a Mach-Zehnder Interferometer (MZI) configuration700, e.g., as shown inFIG.7A. Complete control over the relative phase and amplitude of the two modes702a,702bin dual rail encoding can be achieved by varying the phases imparted by phase shifters706a,706b, and706cand the length and proximity of coupling regions704aand704b.FIG.7Bshows a slightly simpler example of a MZI710that allows for a variable transmissivity between modes702a,702bby varying the phase imparted by the phase shifter706.FIGS.7A and7Bare examples of how one could implement a mode coupler in a physical device, but any type of mode coupler/beam splitter can be used without departing from the scope of the present disclosure.

In some embodiments, beam splitters and phase shifters can be employed in combination to implement a variety of transfer matrices. For example,FIG.8Ashows, in a schematic form similar to that ofFIG.6A, a mode coupler800implementing the following transfer matrix:

Tr=12⁢(111-1).(7)
Thus, mode coupler800applies the following mappings:

❘"\[LeftBracketingBar]"10〉↦❘"\[LeftBracketingBar]"10〉+❘"\[RightBracketingBar]"⁢01〉2⁢❘"\[LeftBracketingBar]"01〉↦❘"\[LeftBracketingBar]"10〉-❘"\[RightBracketingBar]"⁢01〉2⁢❘"\[LeftBracketingBar]"11〉↦12⁢(❘"\[LeftBracketingBar]"20〉+❘"\[RightBracketingBar]"⁢02〉).(8)
The transfer matrix Trof Eq. (7) is related to the transfer matrix T of Eq. (1) by a phase shift on the second mode. This is schematically illustrated inFIG.8Aby the closed node807where mode coupler816couples to the first mode (line612) and open node808where mode coupler816couples to the second mode (line614). More specifically, Tr=sTs, and, as shown at the right-hand side ofFIG.8A, mode coupler816can be implemented using mode coupler616(as described above), with a preceding and following phase shift (denoted by open squares818a,818b). Thus, the transfer matrix Trcan be implemented by the physical beam splitter shown inFIG.8B, where the open triangles represent +i phase shifters.

Similarly, networks of mode couplers and phase shifters can be used to implement couplings among more than two modes. For example,FIG.9shows a four-mode coupling scheme that implements a “spreader,” or “mode-information erasure,” transformation on four modes, i.e., it takes a photon in any one of the input modes and delocalizes the photon amongst each of the four output modes such that the photon has equal probability of being detected in any one of the four output modes. (The well-known Hadamard transformation is one example of a spreader transformation.) As inFIG.6A, the horizontal lines912-915correspond to modes, and the mode coupling is indicated by a vertical line916with nodes (dots) to identify the modes being coupled. In this case, four modes are coupled. Circuit notation902is an equivalent representation to circuit diagram904, which is a network of first-order mode couplings. More generally, where a higher-order mode coupling can be implemented as a network of first-order mode couplings, a circuit notation similar to notation902(with an appropriate number of modes) may be used.

FIG.10illustrates an example optical device1000that can implement the four-mode mode-spreading transform shown schematically inFIG.9in accordance with some embodiments. Optical device1000includes a first set of optical waveguides1001,1003formed in a first layer of material (represented by solid lines inFIG.10) and a second set of optical waveguides1005,1007formed in a second layer of material that is distinct and separate from the first layer of material (represented by dashed lines inFIG.10). The second layer of material and the first layer of material are located at different heights on a substrate. One of ordinary skill will appreciate that an interferometer such as that shown inFIG.10could be implemented in a single layer if appropriate low loss waveguide crossing were employed.

At least one optical waveguide1001,1003of the first set of optical waveguides is coupled with an optical waveguide1005,1007of the second set of optical waveguides with any type of suitable optical coupler, e.g., the directional couplers described herein (e.g., the optical couplers shown inFIGS.6B,7A,7B). For example, the optical device shown inFIG.10includes four optical couplers1018,1020,1022, and1024. Each optical coupler can have a coupling region in which two waveguides propagate in parallel. Although the two waveguides are illustrated inFIG.10as being offset from each other in the coupling region, the two waveguides may be positioned directly above and below each other in the coupling region without offset. In some embodiments, one or more of the optical couplers1018,1020,1022, and1024are configured to have a coupling efficiency of approximately 50% between the two waveguides (e.g., a coupling efficiency between 49% and 51%, a coupling efficiency between 49.9% and 50.1%, a coupling efficiency between 49.99% and 50.01%, and a coupling efficiency of 50%, etc.). For example, the length of the two waveguides, the refractive indices of the two waveguides, the widths and heights of the two waveguides, the refractive index of the material located between two waveguides, and the distance between the two waveguides are selected to provide the coupling efficiency of 50% between the two waveguides. This allows the optical coupler to operate like a 50/50 beam splitter.

In addition, the optical device shown inFIG.10can include two inter-layer optical couplers1014and1016. Optical coupler1014allows transfer of light propagating in a waveguide on the first layer of material to a waveguide on the second layer of material, and optical coupler1016allows transfer of light propagating in a waveguide on the second layer of material to a waveguide on the first layer of material. The optical couplers1014and1016allow optical waveguides located in at least two different layers to be used in a multi-channel optical coupler, which, in turn, enables a compact multi-channel optical coupler.

Furthermore, the optical device shown inFIG.10includes a non-coupling waveguide crossing region1026. In some implementations, the two waveguides (1003and1005in this example) cross each other without having a parallel coupling region present at the crossing in the non-coupling waveguide crossing region1026(e.g., the waveguides can be two straight waveguides that cross each other at a nearly 90-degree angle).

Those skilled in the art will understand that the foregoing examples are illustrative and that photonic circuits using beam splitters and/or phase shifters can be used to implement many different transfer matrices, including transfer matrices for real and imaginary Hadamard transforms of any order, discrete Fourier transforms, and the like. One class of photonic circuits, referred to herein as “spreader” or “mode-information erasure (MIE)” circuits, has the property that if the input is a single photon localized in one input mode, the circuit delocalizes the photon amongst each of a number of output modes such that the photon has equal probability of being detected in any one of the output modes. Examples of spreader or MIE circuits include circuits implementing Hadamard transfer matrices. (It is to be understood that spreader or MIE circuits may receive an input that is not a single photon localized in one input mode, and the behavior of the circuit in such cases depends on the particular transfer matrix implemented.) In other instances, photonic circuits can implement other transfer matrices, including transfer matrices that, for a single photon in one input mode, provide unequal probability of detecting the photon in different output modes.

3. GMZI Implementations

In some embodiments, fast and low-loss optical switch networks can enable scalable quantum information processing using photonic qubits. For example, such networks can be employed within a linear-optical quantum computing (LOQC) system, since many such systems relies on non-deterministic processes of single-photon generation, entanglement generation and fusion measurements, and they also have important applications for quantum communications, such as enabling all-photonic quantum repeaters. However, the present disclosure is not limited to any particular use-case.

Advantageously, one or more embodiments disclosed herein provide for low loss, fast, and minimally-decohering photonic switch networks. Some embodiments provide for switch networks having a minimization of depth and count and are particularly suited for implementations that include active phase shifters, which are historically the largest contributors to the size and amount of noise in switch networks. Examples of switch networks will now be described. Such networks can be used, for instance, in any of the embodiments described above.

Components that can be used in photonic platforms include waveguides, directional couplers, passive and active (fast) phase shifters, crossings, single-photon detectors and heralded single-photon sources (HSPSs). S witch networks can be categorized according to their primary function as follows. N-to-1 (M) muxes (also referred to as N×1 muxes) map one (or multiple M) inputs to designated output ports. The inputs are commonly assumed to be probabilistic and of the same type, although more complicated assumptions apply in some problems. For example, a N-to-4 photon mux extracts groups of four photons from N HSPSs. Sometimes it is necessary to carefully distinguish the number of output (input) ports from the number of principal target outputs (inputs). Most commonly, the excess ports must be populated with the vacuum state, and the switch network is required to access specific distributions (“patterns”) of the outputs (inputs) across the ports. We refer to switch networks as permutation networks when their primary purpose is to rearrange (subsets of) inputs, where the inputs should generally be regarded as inequivalent. Furthermore, switch networks are also classified on the basis of the photonic degree of freedom distinguishing their inputs. Schemes based on space and time are the most common, but the use of frequency, orbital angular momentum, and combinations of multiple degrees of freedom has also been proposed.

In some embodiments, Mach-Zehnder Interferometers (MZIs) may be used which are networks that implement identity or swap operations on two inputs. Two possible realizations of this type of circuit are shown inFIGS.11A and111B.FIGS.11A and11Bshow building blocks of composite switch networks.FIGS.11A and11Bshow 2-to-2 MZIs that implement identity or swap operations on the inputs. The circuits consist of two directional couplers with an active phase shifter1102(gray) on one or both arms between them. The push-pull configuration shown inFIG.11Aalso has a fixed passive −π/2 phase shift1104(white) on one arm and selects between the two operations by setting the top or bottom active phase to −π/2. The configuration shown inFIG.11Buses a 0 or −π active phase to select the operation. Many switch network architectures are built by connecting multiple MZIs to form various topologies.

The Generalized Mach-Zehnder Interferometer (GMZI) is an extension of an MZI with N>2 inputs and M≥1 outputs, shown inFIG.11C. This configuration allows a set of permutations to be performed on the inputs, as discussed in further detail below, making this device a powerful block for the construction of composite N-to-1 and N-to-M switch networks.FIG.11Cshows a N-to-M GMZI made of two passive balanced splitter networks1112,1114(white) and a layer of N active phase shifters1102(gray). Varying the settings of the active phases selects specific permutations of the N inputs and routes them to M>1 output ports.

There are a number of spatial mux schemes that select one of multiple inputs from distinct locations in space. For example, a N-to-1 GMZI can be used as a mux, since it allows routing of any input to a single output port. The advantages of this scheme are its low constant active phase shifter depth (1) and count (N). However, the total propagation distance and the number of waveguide crossings increase rapidly with N. This downside of the monolithic GMZI structure is obviated by constructing composite switch networks of 2-to-1 MZIs, at the cost of increasing the component depth and count. Two examples of N-to-1 schemes of this kind include the “log-tree” and “chain”, both of which can be built with no crossings.

FIGS.12A and12Bshow spatial N-to-1 muxes, with inputs at N spatially-distinct locations (ports).FIG.12Ashows a log-tree mux1202(N=8 example). 2-to-1 MZIs1204form a tree structure with 2(2┌log2(N)┐−1) active phase shifters arranged in ┌log2(N)┐ layers.FIG.12Bshows a chain mux1210(N=4 example). (N−1) MZIs1204are connected through one output and input to form a line. The active phase shifter count is the same as for the log-tree, but the depth varies between 1 and (N−1).

In a “log-tree”, the MZIs form a converging symmetric tree of degree 2, where the chosen input is routed from one of the leaves to the root, as shown inFIG.12A. An asymmetric variant of this scheme, known as a “chain”, includes MZIs cascaded to form a linear topology in which each block selects either the output of the previous block or the new input, as shown inFIG.12B. The depth of the network traversed by the output depends on the chosen input, which can worsen the interference of resources from different chains, due to imbalanced losses and errors. The switching logic of this scheme presents an interesting advantage: while being very simple and entirely local to each individual MZI, it minimizes the amount of error on by selecting the input available closest to the output. Analysis of these three schemes in the context of single photon multiplexing shows that all three architectures require components with performance well beyond the state-of-the-art to achieve a multiplexing efficiency high enough for use in LOQC.

In temporal multiplexing, resources can be input at the same spatial location but different times, and the aim is to produce an output in a specific time bin. This requires networks with fewer components, but the output time bins become longer. There are two main kinds of temporal schemes: designs with storage devices, such as cavities or fiber loops, and designs based on networks of delays The former simply consist of a storage device1302and a single 2×2 switch network1304used to choose whether to store or output each input, as shown inFIG.13A. This can be thought of as the temporal version of a chain mux, and it presents the same advantage in terms of switching logic. The log-tree also has a temporal equivalent known as a “binary-division delay network”. This scheme consists of a series of MZIs1304with delays of different lengths between them, as illustrated inFIG.13B.

FIGS.13A and13Bshow N-to-1 temporal muxes, with inputs in N distinct time bins.FIG.13Ashows a storage loop scheme (time chain). A 2×2 MZI receives one resource per time bin T and routes it to a storage device (a delay line here) or discards it. After N time bins, the chosen input is output. The number of active phase shifters in the path of the chosen input varies between 1 and N.FIG.13Bshows a binary delay network (time log-tree). The scheme comprises a series of ┌log2(N)┐+1 MZIs with delays of lengths 2nT between them, where T is the duration of a time bin at the input and n=0, . . . ┌log2(N)┐−1. The active phase shifter depth scales as with the number of input time bins as ┌log2(N)┐.

The topologies described above can be generalized by replacing each MZI with a GMZI with n inputs, as shown inFIGS.14A-14D. This introduces a trade-off between the active phase shifter depth and count, which decreases with n, and the number of waveguide crossings and propagation distance within each block, which increases with n. In addition, this modification turns temporal schemes into hybrid networks, where multiple spatially distinct resources are input in each time bin. The trade-offs introduced by the parameter n can be exploited to optimize the structure of these schemes for different regimes of physical error rates.

FIGS.14A-14Dshow examples of generalized N-to-1 composite multiplexing networks, obtained by replacing the MZI sub-blocks with n×1 GMZIs1404.FIG.14Ashows a generalized spatial log-tree1402(n=3 example with some first layer GMZIs omitted for simplicity). The degree of the tree is n and its depth is ┌lognN┐.FIG.14Bshows a generalized spatial chain1412. Each stage after the first takes n−1 new inputs, so that the depth of the network varies between 1 and [(N−1)/(n−1)].FIG.14Cshows a generalized delay network (time log-tree)1422. The GMZIs enclose ┌lognN┐ layers of n−1 delays with lengths ni, . . . (n−1)ni, where i=0, . . . , ┌lognN┐−1 is the index of the layer of delays. The number of active phase shifters on a path across the scheme is ┌lognN┐+1.FIG.14Dshows a generalized storage loop scheme1432. n−1 inputs enter the GMZI in every time bin. After ┌N/(n−1)┐ time bins, the GMZI outputs the chosen input.

In applications such as LOQC, which rely on the interference of multiplexed resources, multiplexing is used to produce synchronized outputs. The schemes described so far achieve this by having a single predetermined output spatio-temporal bin. However, when large output probabilities are needed this leads to a large of resources, which can be understood as follows. The number of available resources for a network of size N follows a binomial distribution with average valueN=Np, where p is the probability of an input being populated. The probability of a network successfully producing an output is then pmux=1−(1−p)N. For the typical situation with large N and small p values, the binomial distribution is well approximated by a Poissonian distribution, and so pmux≅1−e−Np. It follows that the average number of inputs scales as Np=−ln(1−pmux), and so the number of available resources that are not used grows rapidly as pmuxapproaches 1. An alternative approach that leads to major efficiency improvements is relative multiplexing. Rather than routing resources to single pre-allocated outputs, this technique uses spatial or temporal log-tree networks to synchronize selected inputs in variable space-time locations, chosen depending on the resources available at any particular instant.

N-to-M schemes in the literature are generally based on the spatial degree of freedom. The simplest of these is a GMZI with more than one output, which has the appealing feature of a single layer of N active phase shifters. However, it only gives access to N permutations, and therefore to limited combinations of inputs. Consequently, the N×M GMZI is more useful when used as a permutation network or as a building block for larger schemes. More flexible routing is achieved by using smaller networks to build composite topologies, known as “switch fabrics”. However, the component depth and count and the size of the crossing networks of these schemes tend to be large, and these downsides trade against each other, making the networks impractical for use in the field of quantum applications.

As an example, Spanke's tree network1500, shown inFIG.15A, allows arbitrary rerouting of the inputs with a constant active switch depth of 2, at the cost of a large number of active phase shifters and waveguide crossings. However, the number of active phase shifters and waveguide crossings scales as O(NM). On the other hand, the scheme1520shown inFIG.15Bavoids large crossing networks, but has an active phase shifter count O(NM) and depth that varies between 1 and M, resulting in variable error rates on the outputs.

FIGS.15A and15Bshow examples of N-to-M switch networks.FIG.15Ashows a Spanke network1500. Two layers of interconnected GMZIs allow arbitrary routing of N inputs to M outputs. The fixed active phase shifter depth of 2 makes this scheme interesting, but the scaling of the number of active phase shifters and crossings scaling as (NM) poses challenges for large sizes.FIG.15Bshows a concatenated GMZI1520. This scheme consists of M concatenated GMZIs with progressively fewer outputs. No complex crossing networks are required between its building blocks, but the O(NM) active phase shifter count and variable depth up to M limit the maximum feasible network size.

For quantum applications, where low error rates are required, N-to-M muxes need to be simplified to reduce the number of active phase shifters, both in total and along the path to the output, as well as the complexity of the crossing networks. The routing algorithms associated with these networks also need to be simplified, to avoid the need for unfeasibly long delays for the inputs. The complexity of the logic is largely determined by its generality, so restricting the operation of the networks to specific tasks is helpful to reduce processing times. These provide guiding principles for the design of additional schemes.

A general switch network implements a set of unitary transfer matrices Uk, where each unitary routes light between a subset of input and output ports. If Ukroutes light from port t to port s, then its sth row and tth column must be zero apart from |Us,t|=1, and similarly for other pairings of input and output ports. The aim of this section is to elucidate the sets of routing operations that are achievable using the simplest form of a many-mode switching network, which is to say one corresponding to transfer matrices Uk=WDkV†, where the unitary matrices W, V†describe passive interferometers, and the Dkform a set of diagonal phase matrices. The phase matrices are implemented physically using a single layer of fast phase shifters acting on every mode, and for simplicity, we will write D in terms of a phase vector d, Ds,t=dsδs,t. The discussion below provides a comprehensive treatment of these switch networks and presents several new constructions.

An important class of switch networks is obtained by considering sets of permutation matrices {Uk=WDkV†}. By adding the fixed passive network corresponding to e.g. U1−1(so, the inverse of an arbitrary permutation from that set), we obtain a new set {UkU1−1}{WD′kW†} of pairwise commuting permutation matrices. So it makes sense to restrict the discussion to the case where the {Uk} are commuting. Switch networks of this type were introduced above as “Generalized Mach-Zehnder interferometers” (GMZIs). Here we need a more precise definition for GMZIs, and we will define them as switch networks having the following specific properties:(i) {Uk=WDkW†} is a set of transfer matrices corresponding to commuting permutations of N modes. The entries of Dkare given by roots of unity (up to an overall global phase factor eiϕkwhich can be chosen at will).(ii) The GMZI switch setting Dkroutes light from input port 1 to output port k.

From these properties it is straightforward to prove that the GMZI must have exactly N settings, and that for any choice of input and output port, there is exactly one setting which routes light between the ports.

From a mathematical standpoint, the set of operations implemented by a GMZI on N modes forms an abelian group of order N. This fact is very helpful here as it allows us to characterize the entire family of GMZIs defined by (i), (ii) using well-known results from group theory (namely the basis theorem for finite abelian groups). In particular, for any GMZI, {Uk} must be isomorphic to a direct sum of cyclic groups, where the order of each of the cyclic groups is a power of a prime number.

To be more concrete, we define groups of commuting permutations([n1, n2, . . . , nr]) generated by matrices C(n1)⊗I(n2)⊗I(n3). . . , I(n1)⊗C(n2)⊗I(n3). . . , I(n1)⊗I(n2)⊗C(n3). . . , where (C(n))i,j, =δi,(j+1 mod n)is a cyclic permutation matrix of size n, and I(nl)is the nl×nlidentity matrix, and ⊗ is the Kronecker product on matrices (The Kronecker product here acts at the level of linear-optical transfer matrices and should not be confused with tensor product operations on quantum state spaces), and the group operation is matrix multiplication. Then, any GMZI on N modes, satisfying properties (i), (ii) above, must implement a set of permutation operations which corresponds to one of the possibilities for([n1, n2, . . . , nr]) with N=Πl=1rnl(up to fixed mode permutations at the input and output).

The different types of GMZIs of fixed size can now be determined using the fact that([n1, n2]) and([n1n2]) are isomorphic if and only if n1and n2are coprime. For example, for N=8, we can identify three fundamentally different types of GMZI:(i)([2,2,2]), permutations are generated by Pauli matrices X⊗I(2)⊗I(2), I(2)⊗X⊗I(2), I(2)⊗I(2)⊗X.(ii) {([4,2])}, permutations are generated by matrices

C(4)⊗I(2)⁢where⁢C(4)=(1111),and⁢I(4)⊗X.(iii)([8]), permutations are generated by matrix

C(8)=(11111111).

We refer to GMZIs implementing([2,2, . . . , 2]), i.e. permutations of the form of swaps on subsets of modes, as “Hadamard-type” GMZIs due the type of passive interferometer which is used (explained below). Similarly, we refer to GMZIs implementing([N]) as “discrete-Fourier-transform (DFT)-type”.

The discussion above characterizes the routing power of linear-optical circuits using one-layer of fast phase shifters in the switch network. In particular, a GMZI on N modes is limited to N routing operations, which is obviously small compared to the N! possible mode rearrangement operations. However, the possibility of implementing different sets of permutation operations is exploited by some of designs for spatial and temporal muxes which are discussed herein. Strictly speaking the limitation to N operations originates in property (ii) above—i.e. the ability to route light from any input port to any output port. More general constructions using a single stage of active phase shifts can be trivially obtained by acting with separate GMZIs on subsets of modes. The resulting transfer matrices are given by the direct sum of the individual GMZIs' transfer matrices. For example, using three MZIs in parallel results in a switch network on 6 modes, allowing 8 different settings. Such a construction can implement abelian groups of permutations of maximum order, which are given in J. M. Burns and B. Goldsmith, Bull. London Math. Soc. 21, 70 (1989), with the number of operations scaling to good approximation as ˜3N/3.

We now turn to linear-optical circuits that can implement the GMZIs defined above. In particular, a circuit that can implement the routing operations([n1, n2, . . . , nr]) on N=Πl=1rnlmodes must enact transfer matrices of the form,
Pk=(C(n1))ki⊗(C(ni))k2⊗ . . . ⊗(C(nr))kr,
with settings vector k where 0≤kl<nlwith l=1, . . . , r. This can be achieved using a circuit with transfer matrices WDkW†as follows:

W=W(n1)⊗W(n2)⊗…⊗W(nr)⁢with⁢(W(nl))s,t=eı2π⁢st/nlnl,
where the W(nl)are DFT matrices; the kth setting of the fast phase shifters is given by
Dk=Dk1(n1)⊗Dk2(n2)⊗ . . . ⊗(C(nr),with (dk(n))s=e−12πks/nfor Dk(n).

One route to constructing practical interferometers for W and W†is to reduce them to networks of beam-splitter and phase-shifter components using generic unitary decompositions from M. Reck et al., Phys. Ref. Lett. 73, 58 (1994), or W. R. Clements et al., Optica 3, 1460 (2016). These decompositions have optical depth (number of optical elements encountered on the longest path through the interferometer) scaling as 2N−3 and N respectively. This means that the transmittance along the longest path will scale with an exponent which is proportional to the size parameter N—which presents a severe experimental limitation for scaling to large GMZI sizes.

GMZI networks—having a lot of special structure—allow for specific decompositions of the type given by equation1600shown inFIG.16, where the matrices S, correspond to crossing networks which reorder modes within the interferometer. Since the subexpressions of the form I(N/nl)⊗V(nl)correspond to repeated blocks of modes interfering according to unitary V(nl), the equation for WinFIG.16can be seen to describe stages of local interference separated by crossing networks. Note also that since the bracketed expressions in the decomposition commute there is some freedom in the configuration of the crossing networks, and some of them can be treated as relabelings of modes rather than physical circuit elements.FIG.17Aillustrates the construction of a Hadamard-type GMZI using the decomposition, as well as simplification which is possible when the GMZI is used as a N-to-1 mux.

FIGS.17A and17Bshow Hadamard-type GMZI constructions: (i) inFIG.17A, illustration of a linear-optical circuit1710for a GMZI on N=16 modes, for which the fast phase shifters are set to configurations of 0 and π to select one of 16 operations from([2,2,2,2]); (ii) inFIG.17B, possible simplification of the circuit1720when only one output port is required—as is the case when the GMZI is used as a N-to-1 mux. The passive interferometers are constructed following the decomposition of W with stages of interference using 50:50 beam-splitters or directional couplers on pairs of adjacent modes, separated by crossings networks. Note that the phases in the physical interferometer generally differ from the constructions given in the main text, and this implies minor modifications for the transfer matrices and phase-shifter settings. An equivalent circuit notation1740is shown inFIG.17B.

For more general GMZI types, we note that the unitary matrices V(nl)can be decomposed into elementary beam-splitter and phase-shifter operations using the generic decomposition methods mentioned above. Alternatively, since the V(nl)are assumed to be discrete Fourier transforms, they can be recursively decomposed into smaller discrete Fourier transforms acting on sets of local modes Inl/(nl′)⊗V(nl), Inl/(nl)⊗V/(nl′)(for any sizes satisfying nl=nl′× nl″) together with crossings networks and additional phase shifts.

One more subtle feature of the GMZI constructions that was remarked on above is that the matrices Dkfor the GMZIs are determined up to a setting-dependent global phase factor eiϕk. In principle these global phases can be freely set over a range [0,2π) (provided the active phase shifters themselves are configured with sufficient phase range). For an application such as single-photon multiplexing, the global phase factors have no role in the operation of the switch network. However, they can be useful if the switch network is applied to only some part of the input states (e.g. single rails from dual-rail qubits) or if it is incorporated in larger interferometers. In these cases, additional functionality can be absorbed into the operation of the switch network without adding extra layers of switching.

This idea is very useful for LOQC, where it is often desirable to multiplex some circuit which generates entangled states, whilst also applying internal adaptive corrections to its output. An example of this occurs when multiplexing Bell states from a standard BSG circuit. This circuit produces a Bell state across four modes with probability 3/16, but the Bell states do not conform to dual-rail qubit encoding (i.e. with qubits allocated to fixed pairs of modes) in a third of cases. Although this problem can be addressed using an additional MZI at the mux output to perform an optional mode-swap operation, a more elegant solution is presented inFIGS.18A and18B.

FIGS.18A and18Bshow examples of larger GMZI to implement adaptive swaps of rails while multiplexing Bell states generated with n2standard BSGs.FIG.18Ashows sending the two rails that might need to be swapped (circled at1804) through a single GMZI1806of size N=n1n2(n1=n2=2 in this diagram) allows multiplexing and permutation operations to be combined while avoiding the need for an additional switching stage.FIG.18Bshows that the modular structure of the GMZI1806can be exploited to apply portions of the circuit at different locations and to optimize the physical implementation. In this example, the network which incorporates the swap operation can be decomposed into two 2-to-1 GMZIs1808with extra directional couplers1810applied at the output of the BSGs and between the two output rails.

In this approach, a mux on n2copies of the BSG implements multiplexing and swap operations, using a size N=n1n2GMZI on n1=2 inner rails from each BSG, and regular n2-to-1 multiplexing for the outer rails. The ability to permute the rails increases the success probability for generating a dual-rail encoded Bell state from 1/8 to 3/16, and thereby decreases the amount of multiplexing needed to reach any particular target output probability by a factor of ˜1.55.

More generally, the transfer matrices associated with a GMZI that implements the routing operations([n1, n2]) are

P(k1,k2)=(C(n1))k1⊗(C(n2))k2=(C(n1)⊗I(n2))k1⁢(I(n1)⊗C(n2))k2.
This can be interpreted as n1separate copies of n2-to-1 GMZIs (second term) with an additional set of permutations of the n1outputs also available (first term). So, permutations of n1rails can be implemented while multiplexing each one n2times by sending all N=n1n2inputs through a single larger GMZI rather than smaller separate ones. The key advantage of this method is that the depth and total number of active phase shifters do not change (1and N respectively).

Using a larger GMZI comes at the cost of increasing the optical depth of the circuit, particularly in terms of waveguide crossings. As seen from the expression of W above, the passive interferometers in a GMZI can be decomposed into smaller networks connected by layers of crossings. This modular structure can be exploited to distribute parts of the circuit across different locations and avoid large on-chip crossing networks. In the BSG example, the implementation shown inFIG.18Bhighlights how the first layer of crossings can be realized in a different way, e.g. using long distance phase-stable optical routing, to mitigate the impact of the largest crossing network in the interferometer.

The discussion so far presented a large family of GMZIs and explained their key properties, taking an approach focused on achievable sets of permutations which is different to earlier works. As well as N-to-1 muxing (potentially with extra functionality as explained above, these GMZIs have assorted applications as building blocks for spatial and temporal muxes. Alternative constructions of GMZIs are also possible, and it is valuable to explore them with a view to minimizing practical requirements on fast phase shifters. However, it is not feasible to exhaust all possible GMZI designs, as some properties for Hadamard matrices are not known. Instead we will highlight some specific new constructions with useful properties.

One observation is that phase swing requirements (where the swing is defined per phase shifter as the difference between the maximum and minimum phase shifts across all GMZI settings) can sometimes be reduced by introducing fixed phase-shift offsets. For some of the constructions above, the phase shifter settings correspond to complete sets of roots of unity, and the phase swing is π for Hadamard interferometers and >π for the other GMZI types. Table 1 shows examples of reduced swing for GMZI sizes N=2,3,4 including examples of GMZIs with reduced phase swing using fixed phase-shift offsets. It is assumed that all the fast phase shifter components are identical and access the same range of phase shifts (which is minimized). Note that the use of offsets necessitates modification of the GMZI transfer matrices by additional phase factors corresponding to setting-dependent “global” phases at the output.

TABLE 1PhaseGMZI typeoffsetsCommentHadamard(−3π/2, 0)Swing reduced from π to π/2, coincidingN = 2with MZI variant in FIG. 11A.DFT(−4π/3, 0, 0)Swing reduced from 4π/3 to 2π/3.N = 3Hadamard(−π, 0, 0, 0)Swing unchanged at π, but for each settingN = 4only one phase shifter is set to π and theothers to 0.

To find some more subtle constructions, we can consider general constraints on GMZIs implementing transfer matrices Uk=WDkV†on N modes, which are required to act minimally as N-to-1 muxes. It is straightforward to prove a lemma stating that (a), V in this case must be proportional to a complex Hadamard matrix (i.e. V must satisfy |Vs,t|=1/√{square root over (N)} as well as being unitary), and (b) the phase vectors dkmust be orthogonal. A simple consequence of this result is that it is never possible to construct any GMZI for which the phase-shifter swing is less than π/2 (since it is never possible to achieve 0 for the real part of (dk, dk′>). Similarly, when the phase-shifter values are restricted to {0, π/2} it is not possible to find more than 2 orthogonal vectors dkfor any even value of N (and never more than 1 for odd values of N), which is to say that it is not possible to do better than a 2-to-1 mux.

As another application of this lemma, one can look for sets of orthonormal phase vectors {dk} and construct a GMZI which uses these as phase settings for a N-to-1 mux, by choosing V to have row vectors vk=dk, and any unitary W with first row vector w1=(1,1, . . . , 1)/√{square root over (N)}. An interesting and non-trivial example of such a set of phase vectors is given in Table 2. More specifically the able below shows examples of six orthogonal phase vectors with a subset d1, . . . , d4having a reduced phase swing of 2π/3 (compared to 4π/3 for the entire set). A N=6 GMZI constructed using these settings can implement a 4-to-1 mux which has phase swing of only 2π/3 (by restricting to the first four phase-shifter settings). Furthermore, it is easily seen that this example is not related to the constructions above since the only possibility would be the GMZI implementing([6])≅([3,2]), for which individual phase settings range on six values (compared to three in Table 2).

TABLE 2Settings for a N = 6 GMZI acting as a 6-to-1 muxd1= (1, 1, 1, e−21π/3, e−21π/3, e−21π/3)/√{square root over (6)}d2= (1, e−21π/3, e−21π/3, 1, e−21π/3, 1)/√{square root over (6)}d3= (e−21π/3, 1, e−21π/3, e−21π/3, 1, 1)/√{square root over (6)}d4= (e−21π/3, e−21π/3, 1, 1, 1, e−21π/3)/√{square root over (6)}d5= (1, e−21π/3, e−41π/3, e−21π/3, 1, e−41π/3)/√{square root over (6)}d6= (e−21π/3, 1, e−41π/3, 1, e−21π/3, e−41π/3)/√{square root over (6)}

Finally, we turn to a new way of using GMZIs when phase settings are modified from those connecting single input and output ports. Taking Hadamard-type GMZIs with transfer matrices Uk=WDkW†on N modes, consider first when the phase vector dk′for Dk′is modified so that −π phases are set to a (common) value −ϕ, while the 0 phases are unchanged. In this case Uk′is modified to

U~k′(ϕ)=e-i⁢ϕ/2[cos⁡(ϕ2)⁢I(N)+i⁢sin⁡(ϕ2)⁢Uk′].

This unitary maps a single photon incident at one input port to a superposition across the mode at the input and the output under the permutation Uk, with weighting controlled by the value of ϕ. Further modification of the phase settings can achieve mappings from one input to arbitrary pairs of output ports—suppose it is desired to map from input port p1to output ports q1and q2, then this can be implemented by finding the (unique) settings k1, k2with U=WDk1(2)W†: pq1(2), and choosing phase vector

d˜=e-i⁢ϕ/2[cos⁡(ϕ2)⁢dk+i⁢sin⁡(ϕ2)⁢dk′].
The transfer matrix for the GMZI is then

U~(ϕ)=e-i⁢ϕ2[cos⁡(ϕ2)⁢Uk+i⁢sin⁡(ϕ2)⁢Uk′],
where the individual phase settings are taken from the set {0, −ϕ, −π, −π, −ϕ}. Note that a second input port p2is also mapped to the pair q1and q2, where UkUk′:p1p2. We call a GMZI used according to the equation above for Ũ(ϕ) a switchable pairwise coupler and it can be useful in spatial and temporal muxes (with the proviso that paired ports receive the vacuum state to avoid contamination of the intended input).

4. Additional Embodiments

While the invention has been described with reference to specific embodiments, those skilled in the art will appreciate that variations and modifications are possible. For instance, the number of stages in an MZL-based phase shifter can be varied as desired, and different combinations of phase shifts can be applied. Phase shifters of the kind described herein can be used in a Hadamard-type GMZI (as described above) or in other GMZI circuits or other circuits.

Further, embodiments described above include references to specific materials and structures (e.g., optical fibers), but other materials and structures capable of producing, propagating, and operating on photons can be substituted.

Although not explicitly shown in the figures, it should be understood that active switching elements can be operated under control of classical control logic and/or classical decision logic circuits. Such circuits can be implemented on-chip with the waveguides, beam splitters, detectors and/or and other photonic circuit components or off-chip as desired. Any of the classical logic circuits described herein can be implemented using a microprocessor, microcontroller, field programmable gate array (FPGA), application-specific integrated circuit (ASIC) or any other digital logic circuitry. In some embodiments, some or all of the classical logic circuits can be implemented in a classical computer system.

It should be understood that all numerical values used herein are for purposes of illustration and may be varied. In some instances ranges are specified to provide a sense of scale, but numerical values outside a disclosed range are not precluded.

It should also be understood that all diagrams herein are intended as schematic. Unless specifically indicated otherwise, the drawings are not intended to imply any particular physical arrangement of the elements shown therein, or that all elements shown are necessary. Those skilled in the art with access to this disclosure will understand that elements shown in drawings or otherwise described in this disclosure can be modified or omitted and that other elements not shown or described can be added. The terms “upstream” and “downstream” are used herein in reference to the direction of photon propagation along an optical path such as an optical fiber or other waveguide and are not intended to imply any particular physical arrangement of waveguides.

A recitation of “a”, “an” or “the” is intended to mean “one or more” unless specifically indicated to the contrary. Reference to “one” of a particular component, feature, or other element is not intended to preclude additional co-existing instances of that component, feature, or other element, unless specifically indicated to the contrary. The use of “or” is intended to mean an “inclusive or,” and not an “exclusive or” unless specifically indicated to the contrary.

This disclosure provides a description of the claimed invention with reference to specific embodiments. Those skilled in the art with access to this disclosure will appreciate that the embodiments are not exhaustive of the scope of the claimed invention, which extends to all variations, modifications, and equivalents.