Topological Qubits as Carriers of Quantum Information in Optics

Winding number is a topologically significant quantity that has found valuable applications in various areas of mathematical physics. Here, topological qubits are shown capable of formation from winding number superpositions and so of being used in the communication of quantum information in linear optical systems, the most common realm for quantum communication. In particular, it is shown that winding number qubits appear in several aspects of such systems, including quantum electromagnetic states of spin, momentum, orbital angular momentum, polarization of beams of particles propagating in free-space, optical fiber, beam splitters, and optical multiports.


Introduction
Topological quantities, such as winding numbers, Chern numbers, and Euler characteristics [1,2], depend on global properties of a system and are unchanged by smooth, local deformations of the space on which they are defined. They have become prominent in areas of physics ranging from particle physics and condensed matter physics to quantum communications and optics. Part of the reason for this surge of interest is their stability against disturbances: since they are integer-valued, rather than continuous functions of the system parameters, small continuous deformations of those parameters often leave them unchanged. In particular, the idea of information processing using quantum bits (qubits) that are encoded into topologically stable quantities has become an active research area, due to the possibilty of constructing quantum computing systems with reduced need for error correction. The most prominent candidate for such topological qubits are anyons [3,4], which are unfortunately challenging to produce and manipulate experimentally. Here we point out that topologically-based qubits appear naturally in optical systems and that these optical topological qubits are often relatively easy to produce, manipulate, and detect, allowing the possibility of their use in topological quantum information processing.
Consider a continuous closed curve s(t) confined to a plane. This plane could be in position or momentum space, for example, or in an abstract Hilbert space. Most simply considered, the winding number ν [5] of s(t) is the net number of counterclockwise turns traced out about the origin as parameter t runs from 0 to 1. The winding number is therefore well described using the polar coordinates s(t) = {r(t), φ(t)} in the plane. The winding number ν ∈ Z, can be given as In optical contexts, the winding number might count the number of rotations of a polarization vector, the amount of orbital angular momentum present, or the number of times a phase rotates as the photon momentum circles the Brillioun zone of a spatially periodic system.
Although this characterization can be generalized to include the winding number relative to any point p in the plane, the winding number about the origin is considered here without loss of generality; see Figure 1. Winding number is a topologically significant quantity, in that it is invariant under smooth deformations of the curve which don't involve crossing the origin. When the origin is crossed, ν changes by a discrete integer amount. A wide range of physical systems whose descriptions have interesting topological aspects have been considered in optics, solid state physics, elementary particle physics, and elsewhere in the study of nature; a number of these can be considered via the winding number or its various generalizations, such as the Chern number [6]. In the quantum theoretical context in particular, quantum states with various values of winding number can be identified. Among the quantum applications is the characterization of distinct topological phases of materials and the topologically protected edge or boundary states between them [7][8][9][10][11][12]. Long known in momentum space in solid state systems, such topological states have recently been studied in photonic quantum walk systems as well [13][14][15][16]. Figure 1. The winding number ν ∈ Z of a continuous curve about some point p counts the number of times the curve encircles p. Here, p ia taken to be the origin, and the curve parametrized as s(t) : R → R 2 in the plane can be given a polar coordinate description via functions {r(t), φ(t)}, with parameter t (0 ≤ t ≤ 1). Specifically, s(t) = r(t) · (cos φ(t), sin φ(t)), with φ(1) = φ(0) + 2πν.
In applications, the fact that the homotopy class of a curve does not depend on the the radial part of the parametrized curve s(t) is of great value; for example, it means that physical quantities depending only on the winding number are very stable against the effects of small perturbations to the system. As a consequence of the lack of radial dependence, one can simply consider homotopy classes of maps from the circle to itself: any such curve is continuously deformable to one of maps of the circle, S 1 → S 1 : S → S n with multiplication being that on the complex unit circle. (Recall that the i-th homotopy group π(S n ) captures the continuous mapping of the i-dimensional sphere S i into the n-dimensional sphere S n .) This does not distinguish two maps that can be continuously deformed from one to the other; only equivalence classes of mappings are pertinent. The homotopy groups of spheres algebraically describe the spheres as topological spaces with an operation of addition defined on their equivalence classes, giving rise to an Abelian group structure. The homotopy classes of maps from a circle to a topological space form the first homotopy group, which is some subgroup of the group of integers Z under addition; the winding number of a curve in the complex plane C (the set of points x + iy) labels its homotopy class [1,17].
Here, we provide a unified overview of the appearance of topological qubits arising in several areas of linear quantum optics, including quantum electromagnetic states (EM) of polarization, spin, momentum, orbital angular momentum (OAM) [18][19][20]. In particular, we show how winding number qubits, the natural description of which is in the so-called Bloch sphere S 2 (cf., e.g., [21], Section 1.2), appear in a sequence of three progressively more useful examples.

Quantum Momentum-Space Winding Number
Consider a complex-valued electromagnetic field, representing an optical beam. Recall that the winding number ν indicates the different equivalence classes in the homotopy group. Each equivalence class consists of the set of closed curves in space for which a point with the parametrized phase circles the origin in the complex plane ν times, as the angle in physical space circles the beam axis once. From the quantum perspective, the EM field can be viewed as possessing excited states, corresponding to discrete numbers of particles, called photons. The spin angular momentum of the photons gives rise to linear polarization states, which can become correlated with a winding number. Superpositions of polarization states then lead to superposition of winding number, or in other words, to winding number qubits.
One context in which such polarization-based winding number qubits can be formed was demonstrated already in the 1980's [22]. For example, consider light traveling along an optical fibre wound helically into a coil characterized by a helix pitch angle θ relative to the local axis of the fiber ( Figure 2). The photons propagate with momentum k(τ) parallel to the fiber axis along an optical path parametrized by arclength τ along the fiber. These photons have a general spin state | k, σ such that where eigenvalue σ characterizes the photon helicity. Superpositions of photon helicity states lead to linearly polarized light. This polarized light can then be sent, for instance, into a single-mode isotropic fiber that is helically wound; cf., [23]. The momentum k is of fixed magnitude, but its direction rotates as the fiber is traversed; accordingly k takes a path on the intersection of a plane with the momentum 3-sphere, or in other words, on a circle in momentum space ( Figure 3a). In this case the number of windings about the axis in momentum space also corresponds to the windings in physical space.   In the process, the quantum state gains a non-trivial geometrical Berry phase, γ(C), which depends on the winding number ν(C), the area Ω(C) subtended by the curve C relative to the origin, and the photon helicity σ.
One can consider a non-trivial geometrical phase known as the Berry phase [24][25][26][27], which arises as a system is taken around a closed loop in parameter space. In his original paper on the subject [24], Berry showed that this phase takes the form where C is the closed path traced out on the Bloch sphere by the quantum state of the photon as it winds along the curve, and where Ω(C) is the solid angle subtended by the curve C relative to the origin of a relevant Hilbert space. This proportionality between geometric phase and areas or solid angles is generic behavior for geometric phases. In the case of interest here, the parameter being varied adiabatically is the direction of the momentum vector as the light winds around the helix, and the polarization state will trace out a curve C on the Bloch sphere as in Figure 3b. The Bloch sphere is often used to describe qubit states in a two-state system. In the current case, the basis qubits |0 and |1 are represented by vertically and horizontally polarized optical states, respectively. These lie at the poles of the Bloch sphere (which is often called the Poincaré sphere in the optical context), with diagonal linear polarizations, circular polarizations, and other equal superpositions of the horizontal and vertical states all lying on the equator. If C is a circle subtending a cone with the origin as apex, and θ is the helix pitch angle relative to the local axis of the fiber, then geometric arguments [22] show that the solid angle subtended by C is given by The geometric phase, Equation (3), is therefore proportional to the winding number. Consider a physical situation in which laser photons enter a non-birefringingent material in such a uniformly helically-wound fiber in the linear polarization superposition state Ignoring any overall dynamical phase factors, the kinetics of the situation dictate that the exit polarization state of the photons leaving the fiber at the other end will be where exp(±iγ) are the phases picked up during parallel transport of the optical state along the fiber. The physical interpretation of the Berry phase is that it measures the rotation of the linear polarization direction [23]. Using a fixed reference beam for comparison, γ can be measured interferometrically, so that one can find the winding number for fixed pitch angle θ: Accordingly, for example, for γ = π and a pitch angle of θ = π 3 , one has achieved winding number ν = 1. A coherent superposition can then be constructed via an optical-fiber Mach-Zehnder interferometer (cf., e.g., [21], Section 1.6) with an entry beamsplitter of variable transmittance that provides one path with such a twisted fiber and the other with an untwisted but otherwise identical fiber and an exit beamsplitter on which the two paths are recombined after an appropriate corresponding relative phase shift between the two arms to produce superpositions of the state of Equation (5), which has winding number 0, and that of Equation (6), that is, a variable-amplitude α i (α ∈ C, i ∈ {0, 1}, |α 0 | 2 + |α 1 | 2 = 1) superposition of winding number 0 and winding number 1: that is, a generic winding number qubit state on the surface of the Bloch sphere.
The topological aspects of this qubit are in some sense trivial; although it is produced by different windings in the coils, the state that leaves the coil is simply a superposition of two computational basis states with a relative phase shift. There is in no sense any topological protection of the output, despite the topological nature of its generation. In the next section, we consider a second example that partially remedies this deficiency.

Winding Number and Quantum OAM
In the study of quantum electromagnetic states of nonzero orbital angular momentum, it is most common to consider those arising in paraxial laser light beams, that is, beams propagating nearly parallel to a given axis. It is convenient to consider these states in the basis formed by the Laguerre-Gauss (LG) mode functions [28] of the field; these can be seen to depend on the winding number associated with the angle about the propagation axis φ; see Figure 4. In this context, the wave has a phase proportional to the azimuthal angle, ψ ∼ e ilφ . The wavefunction is required to be single-valued, which forces parameter l to be integer-valued. It is in fact a conserved quantum number for the system: it characterizes the orbital angular momentum (OAM) eigenvalue of the corresponding field eigenstate and plays the role of winding number, l ≡ ν (see below). z φ Figure 4. An optical wavefront with nonzero OAM. The phase factor of the wave, e ilφ , depends on the azimuthal angle φ, giving the surfaces of constant phase a characteristic corkscrew shape that clearly exhibits winding about the point of the plane penetrated by axis of propagation, z. The Poynting vector, S is perpendicular to the wavefront at all points, rotating as the wave propagates along the z-axis.
Recall that, in general, the EM field possesses energy, momentum, and angular momentum, the last being associated with either spin or spatial degrees of freedom [29]. To see directly the connection to winding number, first note that, in free space, the Poynting vector S = c 2 P of the EM field, indicating the direction and magnitude of the field momentum P is the vector cross-product of the electric ( E) and magnetic ( B) field intensities. It has been theoretically demonstrated [30] and experimentally observed [31] that LG modes of a laser beam carry an OAM for quantum states of linearly polarized light that are (in the paraxial case) independent of the corresponding spin-angular momentum (polarization) of the beam photons. Whereas the polarization is a result of the intrinsic spin-1 nature of the light, the OAM is a result of spatial structure imposed on the shape of its wavefronts. It has also been shown experimentally that the OAM carried by LG modes is conserved in spontaneous parametric down conversion (SPDC), with the OAM of the incoming pump beam determining the combined OAM of the two outgoing daughter modes [32].
In general, fields of laser beams can be described by the product of a Hermite polynomial and a spatial Gaussian, that is, they are described via the Hermite-Gaussian (HG) basis. However, it is convenient and most common to consider the states under consideration in this case in the basis formed by the Laguerre-Gauss (LG) mode functions L α p (x) of the field. If the phase front of the propagating field is helical, the Poynting vector of the beam possesses a non-trivial azimuthal component corresponding to the OAM circulating about the beam axis (along the z direction), tracing out a continuous closed curve in the x-y plane and forming an optical vortex on the axis with annular intensity profiles with radii increasing with the magnitude |l| of the winding number, the modes being describable in that plane.
For the EM field in the paraxial case, only the region close to the propagation axis need be considered [33][34][35][36][37][38][39]; the LG laser mode functions of a paraxial cylindrical beam have the helical phase factor exp (−ilφ) .
This phase becomes undefined on the beam axis, so the axis corresponds to a singular vortex line. Applying the angular momentum operator about the propagation axis,L z = ih ∂ ∂φ , it is clear that the discrete integer value l, often called the topological charge, determines the quantized angular momentum eigenvalue where l ∈ Z with the phase front turning counter-clockwise for l > 0 and clockwise for l < 0. From Equation (10), the phase lφ winds around the origin l times each time the spatial azimuthal angle φ circles the axis once; thus, l is the winding number of the phase function about the phase singularity on the beam axis. Such OAM states are described as solutions of the paraxial Helmholtz equation-i.e., in the approximation of a slowly varying envelope A(r) of the wave U( r) = A( r) exp(ikz) in the z direction, where ∂A ∂z kA and ∂ 2 A ∂z 2 k 2 A and slow means the variation of A( r) is slow within the distance λ = 2π/k of one wavelength along z; cf. e.g., [40], p. 51. As a specific example of such a field, the Laguerre-Gauss (LG) mode functions of the EM field LG l p have winding number l associated with the azimuthal dependence of the phase φ as in Equation (10); the linear increase of phase tilts the wavefronts, producing a corkscrew-shaped wavefront that rotates about the axis like a screw as it propagates (Figure 4). The value of l determines how often the phase completes a full cycle of 2π as one winds around the spatial axis once ( Figure 5). l = 0 l = 2 l = 1 Figure 5. Plots of intensity (top row) and phase (bottom row) in the transverse plane for several LG modes with different l values, but with p = 0 in each case. For l = 0, the phase is constant in each plane; for l = 0, it increases linearly with angle as the axis is circled. The amplitude vanishes on the axis for the l = 0 cases, but is a maximum there for l = 0.
The LG beam field state with orbital angular momentum lh with winding number l and with p radial nodes can be written explicitly as [41] where E 0 is a constant and w(z) = w 0 1 + z z R is the beam radius at distance z; z R = πw 2 0 λ is called Rayleigh range. The arctan(z/z R ) term is the well-known Gouy phase that routinely appears in beam optics. Integer p indexes the radial properties of the mode, with a central dark spot, u l p surrounded by p dark nodal rings; see Figure 6. LG modes with different p values. In addition to the intensity zero on the axis, there are p nodal rings on which the amplitude vanishes (the dark blue rings). The phase is discontinuous in the radial direction at each nodal ring, while increasing linearly in the azimuthal direction in each case.
Creation and annihilation operators can be defined that connect quantum states with different p values, and a differential operator may be defined for which these states serve as eigenstates [42][43][44]. Examples of the intensity and phase of these beams are given in Figure 5 at p = 0 for different values of l; plots of intensity for different p at fixed l are shown in Figure 6. Making use of orbital angular momentum beamsplitters [45], one can form qubit states of angular momentum and thereby of winding number [46]. In particular, one can make use of devices (analogous to polarization beam splitters) that select the optical path based on OAM, with one path for each of the distinguishable OAM states [45].
A coherent superposition of two OAM states can be constructed via an OAM analog of the Mach-Zehnder interferometer (cf., e.g., [21], Section 1.6) with an entry OAM beamsplitter of variable transmittance on which the two paths are recombined after an appropriate corresponding relative phase shift between the two arms to produce superpositions of the state of Equation (5), which has winding number 0, and that of Equation (6), that is, a variable-amplitude α i (α ∈ C, i ∈ {0, 1}, |α 0 | 2 + |α 1 | 2 = 1) superposition of winding number 0 and winding number 1: that is, again, a generic qubit state on the surface of the Bloch sphere. This OAM-based qubit avoids the main deficiency of the helical fiber example of the previous section. Each term in the outgoing superposition state has a well-defined winding number about a central vortex singularity. This winding number can be easily measured by standard methods, and the total winding number exhibits robust topological protection. There is still one drawback, however: in free space propagation, turbulence can cause the central vortex to break up into a set of daughter vortices, each of lower winding number. Although the total topological charge of all of the vortices together is conserved, the daughter vortices can become widely separated under propagation, so that any finite-size detector will register a reduced winding number [47,48]. Thus, the effective winding number can be reduced by external disturbances. In the next section, we consider another physical implementation of optical topological qubits that does not suffer from this problem

Multiports and Winding Number
Recall that quantum optical beams of photons can be divided by beam splitters (BS)-producing superposition states of paths for photons in ordinary space-and by their polarized versions (PBS), which selectively direct them in beams according to polarization characteristics. Ordinary beam splitters can be combined and/or modified to form more complex linear optical systems which elicit specific additional beam properties including winding number.
Standard beam splitters produce two output beams from one input beam, with no backward propagation back out to the input direction. Basic multiport generalizations of standard beamsplitters allow for more than two output beams, but only allow uni-directional net movement of photons-i.e., these quantum systems preclude reversals of beam direction inside themselves. Recently, a generalized multiport (MBS) has been proposed [49] and experimentally demonstrated [50] that allows a reversal of beam direction. These directionally unbiased multiports can be viewed as optical scattering vertices and thus can be combined into networks that give rise to scattering-based quantum walks on graph structures [51][52][53].
In such unbiased multiports, the beamsplitter input/output ports can be attached at their internal vertices, as shown in the inset of Figure 7 left. Adding a phase shifter allows the output of the multiport to be altered, so that different phase shifts at the vertices influence the interference between photon paths; when sufficiently small (see [49]), such a unit is described by an n × n unitary transition matrix U whose rows and columns correspond to the input and output states at the ports-rows and columns correspond to the ports A, B, C.  Figure 7. The directionally-unbiased three-port (left) and four-port (middle). The rectangles after the beam splitters represent a vertex mirror unit consisting of a mirror and a phase-shifter. The distance between each beam splitter and the adjacent mirror unit is half the distance between two adjacent beam splitter. Right figure: more schematic representations of the three-and four-port units as optical scattering vertices that connect to three or four external legs for use in optical graphs or networks.
When the internal phase shifts at all three mirror units are the same, one haŝ where θ is the total phase shift at each mirror unit (comprising the mirror and phase plate).
Consider the optical units shown on the left in Figure 8 [51-53] each of which corresponds to an implementation of the abstract diamond graph on the right. Each diamond graph unit consists of a pair of three-point scattering vertices connected at two edges, with an additional phase shift on one connecting edge. Each of the optical units contains a pair of three-ports; each graph edge corresponds to an optical path and each multiport provides a physical implementation of an abstract three-point scattering vertex. Let the multiports be very small relative to the distance between them, so they can be treated as approximately pointlike in scale; each unit cell contains four multiports and two phase shifters. The phase shifts φ a and φ b in the two diamonds may be different from each other. In the configuration of units shown in Figure 9, the phase plates rotate the polarization by 90 • . A time unit T is defined as the time to go from one diamond graph to the next. Photons present in the areas labeled a and b, bounded by the phase plates, represent two substates at each lattice site. When φ a = φ b , the two triangle graphs inside each cell will have different transmittance; take the transmission probabilities corresponding to the two graphs T a = |t a | 2 and T b = |t b | 2 then. The polarization state flips upon each passage through a phase plate allows the a and b states to be easily distinguished. Photons can be introduced to and released from this system via optical switches and circulators [18].  Figure 8. The diamond graph (left), consisting of two three-edge scattering vertices and a phase shift, has been used theoretically to implement scattering-based quantum walks [51][52][53]. Two directionally-unbiased three-ports connected to each other on two edges, with a phase shift between them provides a means to experimentally implement this graph (right).
One lattice site A and B subsites Figure 9. Diamond graphs, made from two directionally unbiased three-ports, as shown in Figure 8; each cell (indicated by the dashed curves) is made from four three-ports. The two diamond graphs at each site may have different internal phase shifts, φ a and φ b . The red rectangles are phase plates that rotate the polarization by 90 • , so photons in the right half of the cell have vertical polarization (state b), while those in the left half have horizontal polarization (state a). The size of the diamond graphs is exaggerated for clarity-they should be small compared to the distance separating them.
This system is directly analogous to a one-dimensional lattices of atoms, such as a polyacetylene molecule, in which topological effects are well-known. In such systems, there is an alternation of two distinct types of unit cells, leading to a periodic modulation of hopping amplitudes for electrons moving through the structure. As a result of this periodic modulation in space, the energy levels wrap around the Brillouin zone in a topologically nontrivial manner, leading to a nonzero winding number as the quasi-momentum traverses a full Brillouin zone. Here, similar effects occur in an optical system by having unit cells (pairs of diamond graphs) whose internal phase shifts alternate periodically. The transition amplitudes for photons hopping through the system then become periodically modulated in space. As in the condensed matter case, this leads to energy bands that wind nontrivially around the Brillouin zone. The transition amplitudes depend on the phase shifts in the diamond graphs, which can be arranged to be easily controlled (by varying path lengths or by using electro-optically controlled phase plates); as a result, the topological behavior of the system can also be readily controlled by the experimenter, in contrast to the situation in condensed matter systems.
Consider a discrete-time view of such a system. Its behavior can be described by a quantum HamiltonianĤ and the corresponding discrete-time evolution matrixÛ = e −iĤT that takes the system forward one time-step T. (Here, the units are chosen such thath = 1.) There is then one photon collision with a diamond graph per unit time. For initial photon quantum state |ψ(0) , the state at time t = nT is |ψ(nT) =Û n |ψ(0) .
One can then define a quasi-momentum k on a one-dimensional periodic lattice made from a sequence of repeating unit cells, labeled by an integer, m; the quasi-momentum k and the position variable m is dimensionless and discrete. A quantum Hamiltonian expression for a system behaving in this way isĤ = |t a | N ∑ m=1 (|m, b m, a| + |m, a m, b|) (16) where the diamond graph transmission amplitudes for phase shifts are [51][52][53]: where j = a, b.
At each fixed lattice site m or each fixed k, this Hamiltonian can be written in terms of the identity matrix and the Pauli matrices [21].Ĥ describing the dynamics in a two dimensional "internal" subspace labeled by the two substates present at each lattice site. The two energy levels are then separated by a momentum-dependent gap. Without loss of generality, an appropriate redefinition of states allows the quantum state amplitudes to be made real: Transitions in which the photon "reflects off" the diamond graph, returning to the same subcell add a constant background term to the Hamiltonian and can be ignored for current purposes. The remaining terms are those that take a photon from one subcell to an adjacent subcell in a single time step of the corresponding quantum walk.
The nonvanishing coefficients of the Pauli matrices in Equation (18) are now As k goes from 0 to 2π, d then traces out paths labeled by their winding numbers ν about the origin. These winding numbers are functions of photon hopping amplitudes: ν(v, w) = ν(|t a |, |t b |). Because t a and t b vary only weakly with k, the path traced out by d(k) encloses the origin and, so, has nonzero winding number if |t b | > |t a |. Different values of the phase shifts allow solutions with both zero and nonzero winding numbers to occur: Evaluation of Equations (17), (21) and (22) for a range of φ a and φ b values readily shows that varying these phases causes the path traced out by d to shift horizontally and change radius, leading to transitions between winding numbers 0 and 1 due to the lack of enclosure to enclosure of the origin in this space as shown in Figure 10. This system is topologically equivalent to the well-known Su-Schreiffer-Heeger (SSH) system [54] and reduces exactly to the pure SSH model when the transmission amplitudes are independent of k. Different phase values in the two diamond graphs therefore lead to different topological phases, distinguished by their distinct winding numbers. By attaching two chains of these graphs with different winding numbers on each chain there arise localized, topologically-protected states at their mutual boundaries [7][8][9]18]. With appropriate choices of phase shifts φ a and φ b , a persistent probability of finding the photon at the boundary between the two path topologies appearing in the complex system described in Section 4 can be achieved; a localized quantum optical state protected from noise influences by winding-number [8,9,18] arises when two chains of directionally unbiased three-ports, the two chains having different winding numbers [18,20].
In addition, it has been shown [55,56] that at the boundaries between regions of different winding number there is strong reflection of the wave function. By appropriate choice of parameters transmission between regions can be made arbitrarily small; thus winding-number bit flip errors can be made negligible. Furthermore, by entangling the winding number with a more easily measured variable, such as polarization, the value of the winding number bit can be read in a straightforward manner by measurement of the variable entangled to it [20].
Such winding number-polarization entangled states can be achieved in practice by making the phase shifts in the diamond graphs polarization-dependent [20]. In this way, if a photon is prepared in a polarization superposition state |ψ 0 = α 0 | ↑ + α 1 | → (23) and placed into the chain, it will give rise to the winding number superposition state |ψ = α 0 |ν = 0, ↑ + α 1 |ν = 1, → , where the 0 and 1 represent winding number values [20]. Interpreting these winding numbers as bit values, the photons now form a winding number qubit, directly analogous to those of Equations (8) and (13), and spanning the full set of quantum states in a two-dimensional complex Hilbert space. The topological stability of the winding number state suppresses bit flips, but since each winding number is unambiguously associated with a given value of a more easily measured variable (polarization), the bit value can easily be read out by polarization measurements.
The multiport system produces winding number without the disadvantages of the previous examples: the qubits can be readily produced, manipulated, and detected. Examining Figure 7, it can be seen that there are no ports open to the vacuum at any of the beam splitters; as a result, no vacuum noise is introduced, in contrast to other quantum optical approaches. Further, since the qubit values in this system depend only on the discrete topological class of the energy band, small perturbations of any kind generally introduce no errors; thus normal levels of noise in the system have little effect, with robust topological protection of the qubits.

Conclusions
It has been shown that winding number, a topologically significant quantity, can be used as a coding parameter for the carrying of quantum information in linear optical systems, the most common realm for quantum communication. In particular, it was shown that winding number can be coded in several physical quantities of such systems, including quantum electromagnetic states of spin, momentum, orbital angular momentum, polarization of beams of particles propagating in free-space, optical fiber, beam splitters, and optical multiports. By taking two winding number values as bit values, the photons considered are seen as carriers of winding number qubits, because bit pairs have become qubit pairs by being encoded into a pair of superposable quantum states in a two-dimensional complex Hilbert space. Making use of the topologically-induced robustness of winding number qubits opens new possibilities for building error-protection into information processing systems.
Author Contributions: All authors contributed equally to the research and writing of this manuscript.