Simulation of On-Chip Broadband Photon Spin Router Base on Nondiffracting Surface Plasmon Beam Launching

: The development of a photonic device based on a non-diffracting surface plasmon polariton (SPP) beam can effectively improve the anti-interference ability. Furthermore, an easily adjustable on-chip routing device is highly desirable and extremely important in practical optical communication applications. However, no non-diffracting SPP-beam-based spin routing devices with high tunability in multiple degrees of freedom have been reported. In this study, we theoretically designed a simple micro-nano structure to realize a highly adjustable non-diffracting SPP-beam-based spin router using Finite-Difference Time-Domain (FDTD) simulation. The simulation results show that the structure enables spin-controlled nondiffracting SPP-beam directional launching. The launching direction of the nondiffracting SPP beam can be dynamically rotated counterclockwise or clockwise by changing the incident angle. Hence, the routing SPP beam can be coupled to different output waveguides to provide dynamic tunability. Moreover, this device shows good broadband response ability. This work may motivate the design and fabrication of future practical photon routing devices.


Introduction
Due to the strong light-electron coupling at the interface of a metal and a medium, surface plasmon polaritons (SPPs) have enticing properties, such as the subwavelength confinement and surface enhancement of the optical field. Due to the benefits provided by these features, SPPs are widely used in designing on-chip photonics devices, such as optical switches [1], optical logical gates [2,3], and other optical parameter detecting devices [4,5], which will be important in future information processing and all optical communication. In these planar-integrated nanophotonics devices, on-chip manipulation and control of SPPs is of crucial importance. Therefore, a large number of studies have been conducted regarding the control of SPPs on metal surfaces with well-designed nanostructures [6][7][8] by manipulating optical parameters such as wavelength [9], phase [10][11][12], and polarization state [13][14][15] of the incident light. An optical router is an important component in communication systems. When designing optical routing devices for signal transmission, in general, the aim is to uniquely connect a certain launching path of the excited optical signal with a certain optical state of the incident signal. Therefore, the light parameters, such as wavelength and polarization state, which determine the SPPs' directional launching, are the key functionality for designing SPP-based on-chip photon routing devices. In 2013, Jiao Lin et al. reported that the polarization-sensitive plasmonic apertures having the shape of an ear of wheat in a gold film enabled polarization-controlled tunable SPP directional launching [16]. Since then, controllable SPP directional launching [17][18][19][20] and on-chip spin routing devices [21][22][23] have received significant research attention. However, most

Principle and Simulation Results
The 3D structure diagram of our device is shown in Figure 1a and the 2D structure diagram is shown in Figure 1b. To realize the polarization-controlled SPP launching, the basic structure of the device consists of two parallel columns of a nano-rectangular slot with width W and length L, which are etched on a gold film. The distance between the center points of the two adjacent slots on different columns is S 2 + (D/2) 2 , where S is the perpendicular distance of the two columns and D is the distance between the center points of the two adjacent slots on one column. When D is shorter than the effective wavelength of the SPPs (λ spp ), the columns will generate a quasi-plane SPP wave. For the case of plane wave incidence, only the component polarized perpendicular to the long side of the nano-rectangular slot can significantly excite the SPPs. Therefore, the initial phase of the SPPs generated by a single rectangle slot depends on its orientation angle α and the handedness of the incident light. This is the so-called Pancharatnam-Berry (PB) phase [29], which is given by φ(α) = σ ± α − sgn(α) · π 2 , where σ ± is the spin quantum number, i.e., σ + = 1 represents LCP and σ − = −1 represents RCP; sgn() is the sign function; and α ∈ (−π, π) is the polar angle of the orientation direction of the nano-rectangular slot in the local coordinate system (x , y ) (see Figure 1b). Thus, under LCP incidence and slot parameters of α 1 = π 4 and α 2 = −π 4 , for the SPP waves propagating to the left side of the columns, the phase difference between the SPPs exited on column 1 and column 2 is ∆φ le f t = k spp · S − π 2 , whereas for the SPP waves propagating to the right side of the columns, the phase difference between the SPPs exited on column 1 and column 2 is ∆φ right = k spp · S + π 2 , where k spp represents the wavevector of the SPP, and k spp = 2π λ spp . If we design S = λ spp 4 , the phase difference will be ∆φ le f t = 0 and ∆φ right = π. This means that the SPPs wave will generate constructive interference on the left side of the columns and destructive interference on the right side. Similarly, when the incident light is RCP, the opposite happens. In this way, we can realize different directions of SPP launching under different circular polarization excitations. If the columns have a slant angle θ with the y-axis, the propagation direction of the SPP plane wave will have the same angle with the x-axis. To realize nondiffracting SPP-beam excitation, an aperture structure consisting of an upper subarray with slant angle θ and a lower subarray with the slant angle −θ was designed, as shown in the inset of Figure 2a. The structure can generate a cosine-Gauss SPP beam [30], which is a kind of novel nondiffracting beam. This is demonstrated by the FDTD Appl. Sci. 2021, 11, 10643 3 of 13 simulation result given in Figure 2a, which shows a nondiffracting SPP beam launching along the negative x-axis under LCP incidence. This is because, when the incident light is LCP, the SPP beams excited on the upper and lower subarray launch toward the left side and converge on the negative x-axis. Then, the constructive interference of the two converging SPP beams generates the nondiffracting SPP beam. In the simulation, the incident light is a plane wave with LCP, the light wavelength (λ 0 ) is 633 nm, and the dielectric constant of the gold is −11.8351 + 1.24102i. Hence, the wavelength of the SPPs is λ spp = λ 0 (ε m + ε d )/(ε m ε d ) = 605.67nm, where ε m is the real part of the permittivity of metal and ε d is the permittivity of the dielectric medium (i.e., air). The mesh precision of the FDTD model is 10 nm, the boundary condition is a perfect match layer (PML), and the simulation parameters are h = 120 nm, W = 40 nm, L = 200 nm, S = 150 nm and D = 300 nm. The value of the slant angle θ is 10 • and every column has 13 nano-rectangular slots. Obviously, in this aperture structure (the inset of Figure 2a), when the incident light is RCP, the SPP beams excited on the upper and lower subarray will launch toward the right side and cannot converge, and hence cannot generate a nondiffracting SPP beam in this case. Therefore, in order to launch a nondiffracting SPP beam in different directions under different circular polarization incidences, a symmetric pair of subarrays should be etched on the gold film, shown as the inset of Figure 2b. Figure 2b shows the simulated SPP intensity distribution result of this structure under LCP incidence. The result shows that, although the participation of the symmetric pair results in some side lobes, it will not affect the non-diffraction characteristic of the main lobe of the SPP beam. As this structure (the inset of Figure 2b) is mirror symmetric, it is easy to know that it will launch a nondiffracting SPP beam along the positive x-axis under RCP incidence. Hence, we can realize a nondiffracting SPP-beam launching-based spin router. In the simulation of Figure 2b, the distance (δ, see the inset of Figure 2b) between the left and right symmetric structures is δ = 5 µm and the other parameters are that same as those used in Figure 2a.
We then chose the transversal electric field intensity distributions at x = ±13 µm, shown in Figure 2c, to observe the intensity distribution of the nondiffracting SPP beam. It can be found that the transversal electric field intensity distributions are similar to the cosine-Gauss-like beam [30]. We found that the intensity distribution is slightly asymmetrical with the x-axis. This is because the subarrays are not symmetrical with the x-axis.
The key point of a nondiffracting SPP beam in improving the anti-interference ability of photonic devices is that when there is an obstacle on the propagation path, the SPP beam can effectively bypass the obstacle and continue to propagate. This self-healing ability of our device is demonstrated by the FDTD simulation result of Figure 2d, in which a gold cylinder with a radius of 150 nm is set on the propagation path at the point of coordinate (−13 µm, 0) as an obstacle. Obviously, although the obstacle causes the SPPs to diverge, the SPP beam will return to being nondiffracting behind the obstacle. Thus, the proposed structure can improve the anti-interference ability of the photonic device. In addition, the tunability of a device is very important in practical applications. Hence, realizing tunable nondiffracting SPP-beam launching is expected in our device. For this, we investigated the influence of the incident angle on the SPP launching direction. If the light is incident at an inclined angle γ with the z-axis in the y-z plane, it carries an in-plane wave vector k in = k 0 · sin(γ), where k 0 is the free-space wave vector of the incident light. This extra wave vector will change the SPP-beam launching direction [31,32].
participation of the symmetric pair results in some side lobes, it will not affect the nondiffraction characteristic of the main lobe of the SPP beam. As this structure (the inset of Figure 2b) is mirror symmetric, it is easy to know that it will launch a nondiffracting SPP beam along the positive x-axis under RCP incidence. Hence, we can realize a nondiffracting SPP-beam launching-based spin router. In the simulation of Figure 2b, the distance ( δ , see the inset of Figure 2b) between the left and right symmetric structures is = 5µm and the other parameters are that same as those used in Figure 2a.  We then chose the transversal electric field intensity distributions at = ±13 μm, shown in Figure 2c, to observe the intensity distribution of the nondiffracting SPP beam. It can be found that the transversal electric field intensity distributions are similar to the cosine-Gauss-like beam [30]. We found that the intensity distribution is slightly asymmetrical with the x-axis. This is because the subarrays are not symmetrical with the x-axis. The key point of a nondiffracting SPP beam in improving the anti-interference ability of photonic devices is that when there is an obstacle on the propagation path, the SPP beam can effectively bypass the obstacle and continue to propagate. This self-healing ability of our device is demonstrated by the FDTD simulation result of Figure 2d, in which a gold cylinder with a radius of 150 nm is set on the propagation path at the point of coordinate (−13 µm, 0) as an obstacle. Obviously, although the obstacle causes the SPPs to diverge, the SPP beam will return to being nondiffracting behind the obstacle. Thus, the proposed structure can improve the anti-interference ability of the photonic device. In addition, the tunability of a device is very important in practical applications. Hence, realizing tunable nondiffracting SPP-beam launching is expected in our device. For this, we investigated the influence of the incident angle on the SPP launching direction. If the light is incident at an inclined angle γ with the z-axis in the y-z plane, it carries an inplane wave vector light. This extra wave vector will change the SPP-beam launching direction [31,32]. Figure 3 illustrates how kin influences the nondiffracting SPP-beam launching direction when the incident light has an inclined angle with the z-axis in the y-z plane. In the figure, the line segments OP and OQ, respectively, are the equiphase lines of the SPP  Figure 3 illustrates how k in influences the nondiffracting SPP-beam launching direction when the incident light has an inclined angle with the z-axis in the y-z plane. In the figure, the line segments OP and OQ, respectively, are the equiphase lines of the SPP beams excited by the upper and lower nano-rectangular slot arrays when γ = 0; the line segments OP and OQ , respectively, are the equiphase lines of the SPP beams excited by the upper and lower nano-rectangular slot arrays when γ = 0; OF is the nondiffracting SPP-beam launching direction when γ = 0; and OF is the nondiffracting SPP-beam launching direction when γ = 0. It can be easily seen that under γ = 0, the equiphase lines of the SPP beams excited by the upper and lower nano-rectangular slot arrays (OP and OQ ) are rotated counterclockwise with respect to the equiphase lines under γ = 0 by β, which is given by: Appl. Sci. 2021, 11, x FOR PEER REVIEW 5 of 13 = ⋅ .
(1) In this case, the main lobe of the nondiffracting SPP beam will be located on the red dashed line of Figure 3, which has an angle β with the x-axis. Figure 4a-d shows the results of modulating the SPP-beam launching direction by changing γ. We can observe that when the inclined angle is 5° and −5°, the main lobe of the nondiffracting SPP beam is inclined without destroying its nondiffracting property, and the inclined angle of the nondiffracting SPP beam is approximately 5°, which shows good agreement with the calculation result of the equation. We then further extracted the transversal electric field intensity distributions at = ±13μm when the incident light is LCP and RCP, respectively. The results are shown in  In this case, the main lobe of the nondiffracting SPP beam will be located on the red dashed line of Figure 3, which has an angle β with the x-axis. Figure 4a-d shows the results of modulating the SPP-beam launching direction by changing γ. We can observe that when the inclined angle is 5 • and −5 • , the main lobe of the nondiffracting SPP beam is inclined without destroying its nondiffracting property, and the inclined angle of the nondiffracting SPP beam is approximately 5 • , which shows good agreement with the calculation result of the equation.  In this case, the main lobe of the nondiffracting SPP beam will be located on the red dashed line of Figure 3, which has an angle β with the x-axis. Figure 4a-d shows the results of modulating the SPP-beam launching direction by changing γ. We can observe that when the inclined angle is 5° and −5°, the main lobe of the nondiffracting SPP beam is inclined without destroying its nondiffracting property, and the inclined angle of the nondiffracting SPP beam is approximately 5°, which shows good agreement with the calculation result of the equation. We then further extracted the transversal electric field intensity distributions at = ±13μm when the incident light is LCP and RCP, respectively. The results are shown in  We then further extracted the transversal electric field intensity distributions at x = ±13 µm when the incident light is LCP and RCP, respectively. The results are shown in Figure 5a. The figure shows that, despite the change in the inclined angle of the incident light, the intensity distribution curves still approximate the cosine-Gauss distribution and have good consistency for the incident light with different polarization states. However, because the aperture array is not symmetrical with the x-axis, both the intensity and transverse displacement of the main lobe maximum strength point under γ = 5 • are slightly larger than those under γ = −5 • . The self-healing ability of the main lobe of the SPP beam is also demonstrated by setting a gold cylinder with a radius of 150 nm at the point of coordinate (13 µm, 0.94 µm) under RCP incidence and inclined angle γ = 5 • (see Figure 5b). Hence, it is demonstrated that, in practical applications, the SPP launching directions can be effectively regulated by controlling γ without destroying its nondiffracting property.
light, the intensity distribution curves still approximate the cosine-Gauss distribution and have good consistency for the incident light with different polarization states. However, because the aperture array is not symmetrical with the x-axis, both the intensity and transverse displacement of the main lobe maximum strength point under = 5°are slightly larger than those under = −5°. The self-healing ability of the main lobe of the SPP beam is also demonstrated by setting a gold cylinder with a radius of 150 nm at the point of coordinate (13 μm, 0.94 μm) under RCP incidence and inclined angle = 5° (see Figure 5b). Hence, it is demonstrated that, in practical applications, the SPP launching directions can be effectively regulated by controlling γ without destroying its nondiffracting property.

Multi-Nondiffracting SPP-Beam Launching
The above results of the intensity distribution indicate that some sidelobes exist beside the main lobe of the nondiffracting SPP beam. It is obvious that these sidelobes are induced by the interference of the SPP electric field generated by different subarrays located on different sides of the y-axis. By properly controlling the distance between the two subarrays on different sides of the y-axis and the length of the subarrays, we found that the overlapping region of the SPP beams generated by different subarrays located on different sides of the y-axis can be enlarged and the sidelobes can be converted into a new nondiffracting SPP beam. Therefore, this provides an opportunity to realize multichannel nondiffracting SPP-beam directional launching. To demonstrate this, by reducing the distance between the subarrays on different sides of the y-axis and extending their length appropriately, we realized directional launching of three nondiffracting SPP beams. The FDTD simulation results are shown in Figure 6. In the simulation, the distance between

Multi-Nondiffracting SPP-Beam Launching
The above results of the intensity distribution indicate that some sidelobes exist beside the main lobe of the nondiffracting SPP beam. It is obvious that these sidelobes are induced by the interference of the SPP electric field generated by different subarrays located on different sides of the y-axis. By properly controlling the distance between the two subarrays on different sides of the y-axis and the length of the subarrays, we found that the overlapping region of the SPP beams generated by different subarrays located on different sides of the y-axis can be enlarged and the sidelobes can be converted into a new nondiffracting SPP beam. Therefore, this provides an opportunity to realize multichannel nondiffracting SPP-beam directional launching. To demonstrate this, by reducing the distance between the subarrays on different sides of the y-axis and extending their length appropriately, we realized directional launching of three nondiffracting SPP beams. The FDTD simulation results are shown in Figure 6. In the simulation, the distance between the two subarrays on the different sides of the y-axis is δ = 1 µm, every column has 25 nano-rectangular slots, and the other parameters are the same as those used in Figure 2a. It can be seen that there are three nondiffracting SPP beams when the incident light is LCP and RCP. Because the distance between the two subarrays on the different sides of the y-axis does not affect the propagating direction of the SPP wave, the nondiffracting SPP-beam launching direction remains unchanged when adjusting the distance (see Figure 6a). We then further extracted the transversal electric field intensity distributions at x = ±11 µm, which are shown in Figure 6b. Figure 6c-f shows the results of modulating the three SPP-beam launching directions by changing γ. We can observe that when the incidence inclined angle is 5 • and −5 • , the inclined angle of the three nondiffracting SPP beams is approximately 5 • . the two subarrays on the different sides of the y-axis is δ = 1 μm, every column nano-rectangular slots, and the other parameters are the same as those used in Fig  It can be seen that there are three nondiffracting SPP beams when the incident light and RCP. Because the distance between the two subarrays on the different sides o axis does not affect the propagating direction of the SPP wave, the nondiffractin beam launching direction remains unchanged when adjusting the distance (see Fig  We then further extracted the transversal electric field intensity distributions ±11 μm, which are shown in Figure 6b. Figure 6c-f shows the results of modulat three SPP-beam launching directions by changing γ . We can observe that wh incidence inclined angle is 5° and −5°, the inclined angle of the three nondiffracti beams is approximately 5°. The figure shows that the nondiffracting property is maintained when we r the three SPP-beam launching directions by controlling γ . The incident of RC under an inclined angle = 5° was selected as an example to detect the self-healing of the three nondiffracting SPP beams. We placed a gold cylinder obstacle with a of 150 nm at the points (10 μm, 4.96 μm), (11 μm, 0.90 μm), and (10 μm, −3.33 simulate the SPP intensity distribution. The FDTD simulation results are shown in 7. It can be seen that all three directional launching SPP beams have self-healing after the obstacle. The figure shows that the nondiffracting property is maintained when we regulate the three SPP-beam launching directions by controlling γ. The incident of RCP light under an inclined angle γ = 5 • was selected as an example to detect the self-healing ability of the three nondiffracting SPP beams. We placed a gold cylinder obstacle with a radius of 150 nm at the points (10 µm, 4.96 µm), (11 µm, 0.90 µm), and (10 µm, −3.33 µm) to simulate the SPP intensity distribution. The FDTD simulation results are shown in Figure 7. It can be seen that all three directional launching SPP beams have self-healing ability after the obstacle. Considering that the broadband response of a device is expected in practical applications, we further simulated the response of our multichannel nondiffracting SPPbeam launching device in a broadband wavelength range. We chose the two points of (−11 µm, 0) and (11 µm, 0) on the two main SPP intensity lobes under different circular polarization incidences as observation points to observe the directional launching ability of this device, and used the contrast ratio of the two observation points' SPP intensity as the evaluation criterion. Figure 8a shows the SPP intensity contrast ratio of the two points with LCP incident light in the wavelength range from 500 to 1000 nm. It can be seen that the SPP intensity contrast ratio is larger than 10 dB between 565 and 875 nm. The SPP intensity distributions on the metal/air interface when the incident wavelengths are 565 and 875 nm are shown in Figure 8 b,c, respectively. It can be seen that, although the perpendicular distance S = 150 nm of the two columns designed for the wavelength 633 nm cannot perfectly match the at these two incident wavelengths, the device can still realize SPP directional launching.

Simulation of Practical Spin Routing Device
In this section, we demonstrate an application example of routing optical signals to single or multiple subwavelength waveguides using nondiffracting SPP-beam launching. For this, we added a PMMA film with 250 nm thickness to the original metal/dielectric interface. In addition, we chose the communication wavelength window at 980 nm as the working wavelength to design the geometric parameters. First, we changed the structure parameters of the single nano-slot to realize the highest SPP excitation efficiency at an incident wavelength of 980 nm. The optimized width and length of the single nano-slot was chosen to be 40 and 220 nm, respectively. In the optimization, we set h = 120 nm and Considering that the broadband response of a device is expected in practical applications, we further simulated the response of our multichannel nondiffracting SPP-beam launching device in a broadband wavelength range. We chose the two points of (−11 µm, 0) and (11 µm, 0) on the two main SPP intensity lobes under different circular polarization incidences as observation points to observe the directional launching ability of this device, and used the contrast ratio of the two observation points' SPP intensity as the evaluation criterion. Figure 8a shows the SPP intensity contrast ratio of the two points with LCP incident light in the wavelength range from 500 to 1000 nm. It can be seen that the SPP intensity contrast ratio is larger than 10 dB between 565 and 875 nm. The SPP intensity distributions on the metal/air interface when the incident wavelengths are 565 and 875 nm are shown in Figure 8b,c, respectively. It can be seen that, although the perpendicular distance S = 150 nm of the two columns designed for the wavelength 633 nm cannot perfectly match the λ spp 4 at these two incident wavelengths, the device can still realize SPP directional launching. Considering that the broadband response of a device is expected in practical applications, we further simulated the response of our multichannel nondiffracting SPPbeam launching device in a broadband wavelength range. We chose the two points of (−11 µm, 0) and (11 µm, 0) on the two main SPP intensity lobes under different circular polarization incidences as observation points to observe the directional launching ability of this device, and used the contrast ratio of the two observation points' SPP intensity as the evaluation criterion. Figure 8a shows the SPP intensity contrast ratio of the two points with LCP incident light in the wavelength range from 500 to 1000 nm. It can be seen that the SPP intensity contrast ratio is larger than 10 dB between 565 and 875 nm. The SPP intensity distributions on the metal/air interface when the incident wavelengths are 565 and 875 nm are shown in Figure 8 b,c, respectively. It can be seen that, although the perpendicular distance S = 150 nm of the two columns designed for the wavelength 633 nm cannot perfectly match the at these two incident wavelengths, the device can still realize SPP directional launching.

Simulation of Practical Spin Routing Device
In this section, we demonstrate an application example of routing optical signals to single or multiple subwavelength waveguides using nondiffracting SPP-beam launching. For this, we added a PMMA film with 250 nm thickness to the original metal/dielectric interface. In addition, we chose the communication wavelength window at 980 nm as the working wavelength to design the geometric parameters. First, we changed the structure parameters of the single nano-slot to realize the highest SPP excitation efficiency at an incident wavelength of 980 nm. The optimized width and length of the single nano-slot was chosen to be 40 and 220 nm, respectively. In the optimization, we set h = 120 nm and

Simulation of Practical Spin Routing Device
In this section, we demonstrate an application example of routing optical signals to single or multiple subwavelength waveguides using nondiffracting SPP-beam launching. For this, we added a PMMA film with 250 nm thickness to the original metal/dielectric interface. In addition, we chose the communication wavelength window at 980 nm as the working wavelength to design the geometric parameters. First, we changed the structure parameters of the single nano-slot to realize the highest SPP excitation efficiency at an incident wavelength of 980 nm. The optimized width and length of the single nano-slot was chosen to be 40 and 220 nm, respectively. In the optimization, we set h = 120 nm and the refractive index of the PMMA to be 1.49. The simulated relationship between the excitation efficiency and the incident wavelength under these parameters is shown in Figure 9. In the simulation, the incident wavelength ranged from 0.7 to 1.2 mm and the electric field polarization direction is perpendicular to the long side of the nano-slot. The SPP intensity at the point O located on the center point of the nano-slot edge was chosen to evaluate the excitation efficiency. The simulation result shows that the structural parameters can meet our requirement; thus, we used the nano-slot with these geometric parameters to design practical spin-routing devices.
Appl. Sci. 2021, 11, x FOR PEER REVIEW 9 of 13 the refractive index of the PMMA to be 1.49. The simulated relationship between the excitation efficiency and the incident wavelength under these parameters is shown in Figure 9. In the simulation, the incident wavelength ranged from 0.7 to 1.2 mm and the electric field polarization direction is perpendicular to the long side of the nano-slot. The SPP intensity at the point O located on the center point of the nano-slot edge was chosen to evaluate the excitation efficiency. The simulation result shows that the structural parameters can meet our requirement; thus, we used the nano-slot with these geometric parameters to design practical spin-routing devices. For incident light of 980 nm, the dielectric constant of the gold is −39.938 + 2.75855i, so the wavelength of the SPPs is 967.65 nm. We hence set = 242 nm and = 400 nm in our design. To connect the routing device to the second-order on-chip photonic functional element, we designed strip waveguides to transmit the SPP signal. The detailed device arrangement is shown in the inset of Figure 10a. The brown area presents the PMMA and the orange area presents the gold film. The strip waveguides with width of 150 nm and thickness of 250 nm were fabricated on the gold film. In order to couple the nondiffracting SPP beam to different waveguides under the different incidence angles and provide high coupling efficiency between the strip waveguides and the nondiffracting SPPs beam, the strip waveguides have the same slant angle as the nondiffracting SPP beams' propagation directions. Therefore, we can realize on-chip nondiffracting signal routing by controlling the incidence inclined angle and the spin state. We demonstrated this by FDTD simulation. In the simulation, we designed three strip PMMA waveguides which can couple the nondiffracting signal with the inclined incidence angles of −7.5°, 0°, and 7.5°; the distance between the two subarrays on the different sides of y-axis is δ = 5 μm; and every column has 12 nano-rectangular slots. The lengths of the three strip PMMA waveguides are all 5 µm in the horizontal direction. Due to the significantly different transmission losses of the SPPs between metal/air and metal/PMMA interfaces, the nondiffracting SPP signal can only enable stable transmission in the PMMA strip waveguides, which can be seen in Figure 10a-c. Therefore, there is high coupling efficiency between the strip PMMA waveguides and the nondiffracting SPP signals. Moreover, the nondiffracting SPPS beams excited by different inclined angles will couple to different PMMA strip waveguides. In this way, we realized a dynamic tunable on-chip photonic routing device. For incident light of 980 nm, the dielectric constant of the gold is −39.938 + 2.75855i, so the wavelength of the SPPs is 967.65 nm. We hence set S = 242 nm and D = 400 nm in our design. To connect the routing device to the second-order on-chip photonic functional element, we designed strip waveguides to transmit the SPP signal. The detailed device arrangement is shown in the inset of Figure 10a. The brown area presents the PMMA and the orange area presents the gold film. The strip waveguides with width of 150 nm and thickness of 250 nm were fabricated on the gold film. In order to couple the nondiffracting SPP beam to different waveguides under the different incidence angles and provide high coupling efficiency between the strip waveguides and the nondiffracting SPPs beam, the strip waveguides have the same slant angle as the nondiffracting SPP beams' propagation directions. Therefore, we can realize on-chip nondiffracting signal routing by controlling the incidence inclined angle and the spin state. We demonstrated this by FDTD simulation. In the simulation, we designed three strip PMMA waveguides which can couple the nondiffracting signal with the inclined incidence angles of −7.5 • , 0 • , and 7.5 • ; the distance between the two subarrays on the different sides of y-axis is δ = 5 µm; and every column has 12 nano-rectangular slots. The lengths of the three strip PMMA waveguides are all 5 µm in the horizontal direction. Due to the significantly different transmission losses of the SPPs between metal/air and metal/PMMA interfaces, the nondiffracting SPP signal can only enable stable transmission in the PMMA strip waveguides, which can be seen in Figure 10a-c. Therefore, there is high coupling efficiency between the strip PMMA waveguides and the nondiffracting SPP signals. Moreover, the nondiffracting SPPS beams excited by different inclined angles will couple to different PMMA strip waveguides. In this way, we realized a dynamic tunable on-chip photonic routing device. Considering that the wavelength will affect the spin sensitivity of the device, which may eventually lead to a failed identification of the input polarization state, we chose the center points of the two horizontal waveguides at the two sides of the aperture structure, i.e., the points (−17.5 µm, 0.09 µm) and (17.5 µm, 0.09 µm), to observe the spin sensitivity of this device in the wavelength range of 700 to 1200 nm. We used the SPP intensity contrast ratio of the two points as the evaluation criterion. Figure 10d shows the SPP intensity contrast ratio of the two points under LCP incident light in the wavelength range of 700 to 1200 nm, in which the red and black lines represent the SPP intensity at the points (−17.5 µm, 0.09 µm) and (17.5 µm, 0.09 µm) respectively, and the blue line represents the SPP intensity contrast ratio of the two points. It can be seen that the SPP intensity at the point (17.5 µm, 0.09 µm) is always very low, and the SPP intensity at the point (−17.5 µm, 0.09 µm) has a peak with a half height width of about 120 nm (950 to 1070 nm). In this wavelength range, the SPP intensity contrast ratio is larger than 8 dB. This means that the designed device has a wide working bandwidth.
Similarly, we designed a three-channel routing device, which can be used for spincontrolled multicasting. The structure of the strip PMMA waveguides is illustrated in the inset of Figure 11a. Figure 11a shows the coupling of three nondiffracting SPP beams to the three PMMA waveguides when the inclined incidence angle is 0°. In the simulation, the distance between the two subarrays on the different sides of the y-axis is δ = 1 μm, every column has 20 nano-rectangular slots, and the other parameters are the same as those used in Figures 9 and 10. It can be seen that the signal transmitted by the middle PMMA strip waveguide is the strongest because the intensity of the nondiffracting SPP beam of the main lobe is larger than that of the side lobes. All the three strip PMMA waveguides can receive nondiffracting SPP signals correctly and have high coupling efficiency. We also tested the working bandwidth of this multichannel photon router. Similar to Figure 10d, we used the SPP intensity contrast ratio of the center points of one waveguide and its symmetrical waveguide at the other side of the aperture structure as the evaluation criterion. The relationship between the intensity contrast ratios of the three pairs of waveguides and the incident wavelength are shown in Figure 11b-d, respectively. It can be seen that at the center point of the left-hand side middle waveguide, i.e., (−15.5 Considering that the wavelength will affect the spin sensitivity of the device, which may eventually lead to a failed identification of the input polarization state, we chose the center points of the two horizontal waveguides at the two sides of the aperture structure, i.e., the points (−17.5 µm, 0.09 µm) and (17.5 µm, 0.09 µm), to observe the spin sensitivity of this device in the wavelength range of 700 to 1200 nm. We used the SPP intensity contrast ratio of the two points as the evaluation criterion. Figure 10d shows the SPP intensity contrast ratio of the two points under LCP incident light in the wavelength range of 700 to 1200 nm, in which the red and black lines represent the SPP intensity at the points (−17.5 µm, 0.09 µm) and (17.5 µm, 0.09 µm) respectively, and the blue line represents the SPP intensity contrast ratio of the two points. It can be seen that the SPP intensity at the point (17.5 µm, 0.09 µm) is always very low, and the SPP intensity at the point (−17.5 µm, 0.09 µm) has a peak with a half height width of about 120 nm (950 to 1070 nm). In this wavelength range, the SPP intensity contrast ratio is larger than 8 dB. This means that the designed device has a wide working bandwidth.
Similarly, we designed a three-channel routing device, which can be used for spincontrolled multicasting. The structure of the strip PMMA waveguides is illustrated in the inset of Figure 11a. Figure 11a shows the coupling of three nondiffracting SPP beams to the three PMMA waveguides when the inclined incidence angle is 0 • . In the simulation, the distance between the two subarrays on the different sides of the y-axis is δ = 1 µm, every column has 20 nano-rectangular slots, and the other parameters are the same as those used in Figures 9 and 10. It can be seen that the signal transmitted by the middle PMMA strip waveguide is the strongest because the intensity of the nondiffracting SPP beam of the main lobe is larger than that of the side lobes. All the three strip PMMA waveguides can receive nondiffracting SPP signals correctly and have high coupling efficiency. We also tested the working bandwidth of this multichannel photon router. Similar to Figure 10d, we used the SPP intensity contrast ratio of the center points of one waveguide and its symmetrical waveguide at the other side of the aperture structure as the evaluation criterion. The relationship between the intensity contrast ratios of the three pairs of waveguides and the incident wavelength are shown in Figure 11b-d, respectively. It can be seen that at the center point of the left-hand side middle waveguide, i.e., (−15.5 µm, 0.09 µm), the SPP intensity has a peak with a half height width of about 140 nm (935 to 1075 nm), within which the SPP intensity contrast ratio of the point (−15.5 µm, 0.09 µm) and its symmetrical point (15.5 µm, 0.09 µm) is larger than 9 dB. At the center points of the left-hand side upper and lower waveguides, i.e., (−15.5 µm, 4.6 µm) and (−15.5 µm, −4.5 µm), the SPP intensity has a peak with a half height width of about 55 nm, within which the SPP intensity contrast ratio of these points and their symmetrical points on the right-hand side is larger than 6.5 dB. 2021, 11, x FOR PEER REVIEW 11 of 13 µm, 0.09 µm), the SPP intensity has a peak with a half height width of about 140 nm (935 to 1075 nm), within which the SPP intensity contrast ratio of the point (−15.5 µm, 0.09 µm) and its symmetrical point (15.5 µm, 0.09 µm) is larger than 9 dB. At the center points of the left-hand side upper and lower waveguides, i.e., (−15.5 µm, 4.6 µm) and (−15.5 µm, −4.5 µm), the SPP intensity has a peak with a half height width of about 55 nm, within which the SPP intensity contrast ratio of these points and their symmetrical points on the right-hand side is larger than 6.5 dB.

Conclusions
In conclusion, a novel photon spin routing device with a simple structure was proposed. It has good anti-interference ability and dynamic tunability because it is based on nondiffracting surface plasmon beam launching and the launching direction of the nondiffracting SPP beam can be dynamically rotated to couple to different output waveguides by changing the incident angle. For the single non-diffracting SPP-beam directional launching device, the SPP beam can reconstruct itself after passing through an obstacle, and the transmission distance is larger than 20 mm. Moreover, the device also has a good broadband response. The simulation results show that the SPP intensity contrast ratio is larger than 8 dB between 950 and 1070 nm. For the multi-nondiffracting SPP-beam launching device, the non-diffracting transmission distance is larger than 15 mm. The SPP intensity contrast ratio of the main lobe is larger than 9 dB within 140 nm (935 to 1075 nm), and the SPP intensity contrast ratio of the side lobe is larger than 6.5 dB within 50 nm (930 to 980 nm) broadband. These results may be useful for practical on-chip photon spin routing and spin-controlled multicasting.

Conclusions
In conclusion, a novel photon spin routing device with a simple structure was proposed. It has good anti-interference ability and dynamic tunability because it is based on nondiffracting surface plasmon beam launching and the launching direction of the nondiffracting SPP beam can be dynamically rotated to couple to different output waveguides by changing the incident angle. For the single non-diffracting SPP-beam directional launching device, the SPP beam can reconstruct itself after passing through an obstacle, and the transmission distance is larger than 20 mm. Moreover, the device also has a good broadband response. The simulation results show that the SPP intensity contrast ratio is larger than 8 dB between 950 and 1070 nm. For the multi-nondiffracting SPP-beam launching device, the non-diffracting transmission distance is larger than 15 mm. The SPP intensity contrast ratio of the main lobe is larger than 9 dB within 140 nm (935 to 1075 nm), and the SPP intensity contrast ratio of the side lobe is larger than 6.5 dB within 50 nm (930 to 980 nm) broadband. These results may be useful for practical on-chip photon spin routing and spin-controlled multicasting.