Pseudospin-1 Physics of Photonic Crystals

We review some recent progress in the exploration of pseudospin-1 physics using dielectric photonic crystals (PCs). We show some physical implications of the PCs exhibiting an accidental degeneracy induced conical dispersion at the Γ point, such as the realization of zero refractive index medium and the zero Berry phase of a loop around the nodal point. The photonic states of such PCs near the Dirac-like point can be described by an effective spin-orbit Hamiltonian of pseudospin-1. The wave propagation in the positive, negative, and zero index media can be unified within a framework of pseudospin-1 description. A scale change in PCs results in a rigid band shift of the Dirac-like cone, allowing for the manipulation of waves in pseudospin-1 systems in much the same way as applying a gate voltage in pseudospin-1/2 graphene. The transport of waves in pseudospin-1 systems exhibits many interesting phenomena, including super Klein tunneling, robust supercollimation, and unconventional Anderson localization. The transport properties of pseudospin-1 systems are distinct from their counterparts in pseudospin-1/2 systems, which will also be presented for comparison.


Introduction
Graphene has become a fertile platform to explore phenomena related to Dirac particles predicted in fundamental physics and to realize peculiar physical phenomena in materials science [1][2][3][4][5][6][7][8][9][10][11][12][13][14][15]. These interesting phenomena include unconventional half-integer quantum Hall effect in graphene when subjected to magnetic fields [4][5][6], Klein paradox that a classically forbidden region is transparent for Dirac electrons [7,15], weak antilocalization due to the destructive interference between two counter-propagating backscattering waves [8,9], Zitterbewegung of Dirac electrons in the presence of confining potentials [10,11,15], and supercollimation of electron wave packets in graphene subjected to onedimensional (1D) disordered potentials [14]. The low-energy quasiparticles in graphene can be described by a massless Dirac equation and the wave function can be expressed in the form of a two-component spinor [15][16][17] due to the two sublattice degrees of freedom in graphene. As such, graphene is called a "pseudospin-1/2" system. We note that the spin here is not the intrinsic spin of electrons, but a pseudospin referring to the spatial degree of freedom of the wave function. Such pseudospin-1/2 systems can also be found in other systems [18][19][20][21][22][23][24], such as topological insulators in which the surface state dispersion exhibits a Dirac cone [18,19], and photonic and phononic crystals in triangular or honeycomb lattices in which Dirac cones are found at the corners of the Brillouin zone [20][21][22][23].
With the rapid progress of experimental techniques, systems with higher pseudospin values, such as pseudospin-1, have been constructed using various artificial materials via a fine tuning of system parameters. A pseudospin-1 system is characterized by two linear bands meeting and intersecting with an additional flat band at a Dirac-like point. Such systems have attracted quite a bit of attention. The most typical lattice structure to realize pseudospin-1 systems is the Lieb lattice. Its sublattice symmetry protects the existence of a flat band. Pseudospin-1 systems have been realized experimentally in different artificially constructed systems. For example, it is found that the above Dirac-like dispersion can be realized at the corner of the Brillouin zone by placing bosonic cold atoms into an optical Lieb lattice [25][26][27]. In addition, a photonic Lieb lattice formed by an array of optical waveguides can also exhibit the Diraclike dispersion [28][29][30][31]. More recently, two groups presented two different methods for realizing an electronic Lieb lattice 2 Research through atom manipulation with a low-temperature scanning tunneling microscope (STM) [32,33]. Drost et al. removed atoms from a chlorine layer placed on top of a Cu(100) crystal surface, leaving the desired Lieb lattice formed by the atomic vacancies [32]. In contrast, Slot et al. achieved the Lieb lattice not by vacancies but through adding carbon monoxide (CO) molecules to the top of a Cu(111) surface [33]. Theoretical works also predict that such a Dirac-like cone can be found in artificial crystals of ultracold atoms in Dice (T 3 ) or stacked triangular lattices [34][35][36][37][38][39][40], and some electronic materials, such as transition-metal oxide SrTiO 3 /SrIrO 3 /SrTiO 3 trilayer heterostructures [41], blue phosphorene oxide [42], and graphene-In 2 Te 2 bilayer [43]. In those systems, the pseudospin-1 character emerges from the embedded sublattice degrees of freedom. In addition, some two-dimensional (2D) photonic crystals (PCs) are found to carry a Dirac-like cone at k=0 induced by the accidental degeneracy of monopole and dipole excitations [21,22,[44][45][46][47][48][49][50], which combine to give three degrees of freedom. For such PCs, calculations using effective medium theory show that both effective permittivity and permeability reach zero at the Dirac-like point frequency [22,[44][45][46][47][48][49][50]. Above/below the Dirac-like point frequency, the effective parameters are positive/negative. Different pseudospin number leads to different boundary conditions and makes a marked difference on the transport properties of the pseudospin-1 systems. Many interesting transport phenomena, which are different from those in pseudospin-1/2 systems, have been predicted in pseudospin-1 systems. For example, in the presence of 1D potential barrier, there is a so-called super Klein tunneling effect for pseudospin-1 systems [26,27,35,47,48], that is unity transmission for all incident angles when the incident energy is half of the potential barrier. In a pseudospin-1 superlattice formed by a Kronig-Penney type of photonic potential, an electromagnetic (EM) wave packet can propagate in the superlattice without any distortion of shape [47]. We call such phenomenon "supercollimation". When pseudospin-1 systems are subjected to 1D disordered potentials, two unconventional localization behaviors are found for obliquely incident waves. One is nonuniversal critical behavior for which the critical exponent of the localization length depends strongly on the type of disorder [51]. Another is the existence of a minimum localization length at some critical disorder strength beyond which the waves become less localized [52].
In the presence of a circular potential barrier, a perfect caustic phenomenon can occur for large scatterer size when the incident energy equals half of the barrier due to the super-Klein tunneling effect [53]. Furthermore, in the low-energy regime, a superscattering phenomenon occurs for an arbitrarily weak scatterer, i.e., extraordinarily strong scattering characterized by an unusually large cross section [54]. It was also reported that when a spatially uniform electric field is suddenly applied to an electronic pseudospin-1 system, the resulting current can be enhanced by the flat band in both the linear and nonlinear response regimes, compared with that in the pseudospin-1/2 system [55].
Unlike artificially ultracold atom and electronic Lieb lattices, PCs do not require extremely low temperature and atomic level manipulation. The photonic pseudospin-1 system was first experimentally demonstrated in the microwave frequency regime [44] and later in the optical regime [46]. The notion of pseudospin-1 was subsequently extended to aperiodic systems, as conical dispersions at k=0 can also be realized in some photonic quasicrystals as well as the effective zero refractive index [50].
Since pseudospin-1 systems were proposed to realize using artificially constructed ultracold-atom systems around nine years ago, many new interesting phenomena have been discovered in the past few years as mentioned before. While recent experiments demonstrate the viability of realizing pseudospin-1 systems with ultracold-atom or electronic Lieb lattices, it still remains challenging to achieve a long enough coherence length, which makes it difficult to have complex junctions in such systems to observe the theoretically predicted novel behaviors experimentally. Thus, in the review, we give a comprehensive review of the theoretical predictions and experimental progress on the pseudospin-1 systems with a focus on the newly discovered phenomena in PC-based systems. Due to the ease of fabrication and absence of complex interactions between photons, these phenomena may finally be tested in systems formed by PCs. In Section 2, we introduce the relation between the pseudospin-1 dielectric PCs and the zero refractive index effective media, and the experimental realizations of such PCs. In Section 3, we show the calculation of the Berry phase for a loop encircling the Dirac-like point. In Section 4, we relate the Dirac-like cone dispersion of PCs to the spin-orbit Hamiltonian description of pseudospin-1 based on the zero refractive index effective medium description, and we show the photonic counterpart of the gate voltage in graphene. The 1D transport phenomena in pseudospin-1 systems are reviewed in Section 5, including Klein tunneling, supercollimation and unconventional Anderson localization behaviors in disordered 1D potentials for oblique incident angles. A brief summary is given in Section 6.

Effective Medium Description of Pseudospin-1 PCs
It is well known that the lattice symmetry of the 2D PCs with a triangular lattice guarantees the existence of Dirac cones at the Brillouin zone corners with a pseudospin-1/2 character [20][21][22][23]. However, effective medium theory is only good for small k vectors (long wavelength) and is not applicable in general for Dirac cones located at the Brillouin zone boundary. At the Brillouin zone center, the lattice symmetry alone gives quadratic dispersions only. It is found that for some 2D dielectric PCs with C 4 or C 3 symmetries, e.g., square or triangular lattices, the monopole and dipole excitations can become triply degenerate by fine tuning of the radii or permittivity of cylinders so that a Dirac-like cone appears at k=0, in which the conical dispersions intersect with an additional flat band at the Dirac-like point [22,44,45]. Such degeneracy is not given by the lattice symmetry. Instead, it occurs only for some specific values of the radius and permittivity of the dielectric rods composed of the PCs and hence the degeneracy is accidental. As an example, in Figure 1(a), we show the three eigenfrequency surfaces near  k=0 calculated for the transverse electric (TE) modes for a 2D square lattice PC with dielectric cylinders in air. Here, TE modes refer to the EM modes with the electric field along the axis of cylinders. The inset of Figure 1(a) shows the unit cell of the PC. The permittivity and radii of the dielectric cylinders are taken as = 12.5 and = 0.2 ( is the lattice constant). It can be clearly seen that near k=0, a Dirac cone (blue) meets and intersects with a nearly flat band (green) at the Dirac-like point frequency. The modes in the Dirac cone are transverse EM modes with a finite group velocity, while those modes in the flat band are longitudinal EM modes with the locally averaged magnetic field vector aligned with the k-vector and having nearly zero group velocity [44,45]. The optical properties of such 2D PCs exhibiting a Diraclike cone near k=0 is the same as an effective medium with zero refractive index at the Dirac-like point frequency. In Figure 1(b), we show the effective permittivity ( ) and permeability ( ) near the Dirac-like point frequency of the PC calculated by the effective medium theory [56]. We note that both (red solid circles) and (blue open circles) reach zero at the Dirac-like point frequency = 1.0826 ( / ) (c is the speed of light in vacuum), and their dispersions are linear with a frequency change given by = − , giving rise to a nearly constant effective impedance as shown in the inset of Figure 1(b). As we will show later, these features are very important when we map the Diraclike cone dispersion at k=0 onto a spin-orbit Hamiltonian of pseudospin 1 based on the effective medium description. It should be pointed out that the effective material constitutive parameters of 2D PCs are anisotropic. If we assume that the PCs are arranged periodically in the xy plane and the rods are along z-direction, the effective permittivity and permeability can be described by diagonal tensors ← → = diag( , , ) and ← → = diag( , , ), respectively, with = and = for triangular and square lattices [57]. For TE modes, only , , and enter the wave propagation problem. Thus, the effective parameters and shown in Figure 1(b) are actually = and = = . Compared with a single zero material ( = 0 or = 0 but not both), the double-zero material here has a finite impedance, which is a desirable feature in eliminating strong reflections from the interface.
The PC with = = 0 was first realized experimentally in the microwave regime by arranging alumina rods ( = 3.75mm and = 8.8) in a square lattice ( = 17mm) in air [44]. Using this PC structure, some unique properties arising from zero refractive index were demonstrated experimentally, such as cloaking and lensing [44]. The PC at optical frequencies was implemented by a stack of square-cross-section silicon rods ( = 13.7, width=260 nm, and = 600nm) embedded in SiO 2 ( = 2.25) [46]. The experimentally observed angular selectivity of transmission and directive spontaneous emission from quantum dots placed within the material provide direct evidence that the PC has indeed a near-zero refractive index near the Diraclike point frequency [46]. We note here that since the PC is all dielectric, it is much less absorptive than conventional metal-based metamaterials in which the metallic components induce large ohmic loss at high frequencies. Periodicity is not a necessary condition to realize the Dirac-like cone at k=0. The Dirac-like cone dispersion and zero refractive index were also realized experimentally in a 2D 12-fold photonic quasicrystal constructed by placing dielectric rods at the 2D dodecagonal lattice sites generated by the square-triangle tiling model [50].

Berry Phase of the Dirac-Like Cone
It is well known that the Dirac cone of graphene and its photonic counterpart give rise to a Berry phase of [15,58]. It is interesting to know if the additional flat band affects the Berry phase of the Dirac-like cone. In general, one can calculate numerically the eigenmodes of the PC to obtain the Berry phase. For square or triangular lattice PCs composed of cylindrical rods, one can conveniently use the 4 Research multiple scattering theory (MST) to obtain the eigenmodes [44,45]. For the Dirac-like cones studied in this review, it is sufficient to consider the monopolar ( = 0) and dipolar ( = ±1) terms only [44]. Here m denotes the angular momentum number of the eigenmode. The MST equations can be expressed as [44] ( where 0 , 1 , and 1 are real, and and k are the magnitude and angle of k in the polar coordinate [44,45]. By letting the determinant of the 3 × 3 matrix in (1) equal to zero, one can obtain three solutions: , which corresponds to the flat band. When = ±1, the eigenvectors are |Φ ,k ⟩ = , corresponding to the two linear bands [45]. Using the three eigenvectors, one can obtain the Berry phase for each band of the Dirac-like cone; that is, = ∮ ⟨Φ ,k | ∇ k Φ ,k ⟩ ⋅ k = 0 [45], which is different from that of the Dirac cone in grapheme [15].

Spin-Orbit Hamiltonian Near k = 0 from Maxwell's Equations.
We start from the Maxwell's equations and apply effective medium theory to demonstrate that the Diraclike conical dispersion of PCs near k = 0 can be related to an effective spin-orbit Hamiltonian of pseudospin 1. We will assume that bands forming the Dirac cone are derived from monopole and dipole excitations and under those circumstances, we can safely apply effective medium theory to describe EM wave propagation in 2D PCs near a Dirac-like point at the Brillouin zone center. The waves are assumed to be TE modes propagating in the effective medium (xy plane) of 2D PCs shown in Section 2, with electric field along zdirection and magnetic field in the xy plane. As we discussed in Section 2, for the effective material constitutive tensors of PCs, we only need to take into consideration , , and . We note that all three parameters are functions of frequency . By letting = and = = , the Maxwell's equations can be expressed in the following matrix form [47]: wherẽis̃= (− − , , − ) for transverse modes and̃= ( − , , − − ) with being spatially independent for longitudinal modes, respectively. We note that̃are different from the eigenvectors of PCs in Section 3, where the three components correspond to the mode amplitudes of Mie scattering with the angular momentum number m=-1, 0, and 1, respectively. As shown in Section 2, for PCs with a conical dispersion at k = 0, the effective permittivity and permeability approach linearly to zero at the Dirac-like point frequency , and we have ≅ ( − )̃and ≅ ( − )̃in the neighborhood of [44,47]. The system parameters̃= ( / )| = and̃= ( / )| = must be positive definite in order for the energy density to be nonnegative [44,47].
Through the Fourier transform of̃in k-space, i.e., for longitudinal modes, with 2 ) ∫ − k⋅r r. Using the above linear approximations of and , (2) can be transformed to the following equation in k-space where = − , and ( ) is the projection of the 2D wavevector k on the x (y) axis. We note that, for longitudinal modes with any nonzero k, the component (k) is zero. Using a similarity transformation with = diag(2√̃/̃, √ 2, 2√̃/̃) and̃( k) = (k) , we can rewrite (3) as Research 5 As can be seen in Figure 1(b), the effective constitutive parameters of the PCs are negative ( , < 0) below the Diraclike point frequency , while they are positive ( , > 0) above . Right at , the effective refractive index is zero ( = = 0). The sign change of in (4) is directly related to the sign change of and in (2) when the frequency crosses over . Thus, the positive and negative index behaviors of the system are treated by (4) in a unified framework. Equation (4) can be further rewritten as by using the following spin-1 matrices, i.e., S =̂+̂with Here V = 1/√̃̃is the group velocity of the conical dispersion [47]. Equation (4) or (5) represents a spin-orbit interaction of pseudospin 1. Three normalized eigenvectors can be obtained for (4) or (5): Here k is the angle between the wavevector k and the positive x-direction. The three bands of the Dirac-like cone can be clearly seen from the corresponding eigenvalues = V |k| ( = 0, ±1). Here = ±1 correspond to the upper ( = +1) and lower ( = −1) conical bands, while = 0 corresponds to the flat band. The values also represent three pseudospin states in which the pseudospin is either parallel ( = +1), perpendicular ( = 0), or antiparallel ( = −1) to the wavevector k. The chiralities of the upper and lower cones are opposite. From (7), we can obtain that the Berry phase equals 0 for a loop around the Dirac-like cone, according to the formula = ∮⟨ (k) |∇ k | (k) ⟩ ⋅ k. This is consistent with the results obtained from the MST theory for a realistic PC structure shown in Section 3. For comparison, the corresponding Berry phase is for pseudospin-1/2 systems (e.g., graphene). This shows immediately that while both pseudospin-1 and pseudospin-1/2 systems are characterized by conical dispersions, the physics involving wave scattering are quite different. Next, we will show that the normalized eigenvector (7) corresponds to a transverse mode in a positive medium ( , > 0), while (k) −1 corresponds to one in a negative medium ( , < 0). For the TE modes under investigations, the electric field (k) in the normalized eigenvectors (k) ( = ±1) is taken as 1. For arbitrary values of (k) , we have (k) = (k) (k) ±1 . From the Maxwell's equations, the magnetic field H (k) can be expressed as Using the first-order approximation of , i.e., wherek = k/|k| is a unit vector indicating the direction of wavevector k. We see that = +1 describes a transverse mode in positive media for which the three vectors E (k) , H (k) , andk form a right-handed triad, while = −1 describes one in negative media for which they form a left-handed triad. In terms of electromagnetic fields, the eigenvectors are rewritten as ) . (8) When = 0, the eigenvector (k) corresponds to a longitudinal mode. In terms of electromagnetic fields, it can be expressed as [47], ) .
It should be mentioned that in our effective medium description the group velocity of the flat band is zero. However, in realistic PC structures, for k near the Dirac-like point, the group velocity can be small and negligible. When k is far away from the Dirac-like point, the effective medium description fails so that the band is not flat any more, i.e., it has finite group velocity and the band modes are not longitudinal modes either. Due to the orthogonality between the transverse and longitudinal modes, these flat band modes near the Dirac-like point are almost impossible to be excited by incidence of a propagating transverse mode in experiments. Thus, we ignore these longitudinal modes when we study the transport properties of the transverse waves in Section 5. To detect these flat band modes experimentally, one may put a point source inside the 2D PCs. In this case, a localized excitation appears in the vicinity of the source, which moves with the position of source. For ultracold atoms in a Lieb lattice, the atoms can be transferred coherently into the flat band by actively engineering the population and phase on each lattice site [25]. For an array of optical waveguides with a Lieb lattice configuration, the flat band mode can be excited by a structured excitation of the three sublattices [28][29][30][31].

Photonic Counterpart of Gate Voltage.
It is well known that electron transport in graphene can be manipulated by applying a gate voltage to form complex junctions, with which many interesting transport phenomena have been observed experimentally or predicted theoretically, such as Klein tunneling [7,15] and supercollimation of electron wave packets [12][13][14]. One might be interested to know whether we can find the photonic counterpart of gate voltage in graphene which can shift the Dirac-like cones rigidly up or down in frequency with the group velocity essentially unchanged, and whether the application of such photonic potentials in PCs will result in any new physical phenomena. It is also well known in classical wave physics that one possible way to manipulate the propagation of waves is to make the structural variation in the system to form a certain structural pattern.
We will see that, for dielectric PCs, a particular kind of structural change is equivalent to applying a gate voltage in electrons. Such structurally induced photonic potential will play a significant role in the dynamics of pseudospin-1 wave packets.
Here we show that such a frequency upshift or downshift of the Dirac-like cone can be easily achieved by changing the length scale of a dielectric PC [47]. Suppose that a dielectric structure with a spatial configuration of dielectric (r) is scaled uniformly by a factor to achieve a new configuration (r) = (r/ ), the mode frequency and wavevector k of the new configuration can be obtained by rescaling the mode frequency and wavevector k of the old one through the relations = / and k = k/ , according to the scaling properties of the Maxwell's equations [59]. We assume that the materials consisting of the structure have frequencyindependent permittivity in the operational frequency range.
where and indicate the respective group velocities before and after the scaling. The change of length scales of dielectric PCs plays the role of the gate voltage in graphene, shifting rigidly the conical dispersions of PCs in frequency. As a photonic analog of electron potential, the length-scaling induced Dirac-like point frequency shift (Δ ) is effectively a photonic gate potential with = Δ = − . If the local length scales of the PCs are modulated along one direction, for example, the x-direction, we can construct a 1D photonic potential ( ) = ( ) − . The total Hamiltonian now becomes I is a 3 × 3 identity matrix. From (11), pseudospin-1 waves propagating along the x-direction do not experience backscattering from the 1D photonic potential, which leads to one-way transport phenomenon. From the EM wave point of view, this is a result of impedance matching because the impedance is a scale invariant constant near the Diraclike point frequency, as shown analytically in [47] and numerically in the inset of Figure 1(b). To see how pseudospin-1 waves are scattered by 1D potential, we need to consider the boundary conditions at the interface. Assuming that = ( 1 , 2 , 3 ) is a general solution of (11), we can obtain three boundary conditions at an interface = 0 by integrating the wave equation = over a small interval ∈ [ 0 − , 0 + ] along the xdirection and taking the limit → 0 [27,35,47], If the wave function takes the form of (8) or (9), we directly obtain that (12)- (14) are equivalent to the continuity of , , and , respectively, i.e., the boundary conditions required by TE waves. We note here that only the first two boundary conditions, i.e., Eqs. (12) and (13), are independent because the continuity of 2 implies the continuity of . It is consistent with EM wave theory in which, for time-harmonic fields, the continuity of the tangential E and H components across interfaces implies the continuity of normal B and D components [60].

Klein Tunneling of Pseudospin-1 EM Waves.
We now consider the scattering of pseudospin-1 EM waves travelling in the xy plane by a 1D square photonic potential barrier. The scattering process is schematically depicted in Figure 2(a). The wave functions in the three regions of Figure 2(a) can be written in terms of the eigenvectors in (7). The wave function in region I can be expressed as with = sgn( ), = arctan( 0 / 0 ), 0 = |q 0 | cos , 0 = |q 0 | sin , and |q 0 | = | |/V . In region II, we have  [47].
with = + arctan( 0 / 1 ), = sgn( − 0 ), and 1 = −√( 0 − 2 ) 2 /(V 2 ) − 2 0 . In region III, we have A transfer matrix M can be defined by the relation, where 0 =̃0 0 0 , 0 =̃0 − 0 0 , 1 =̃1 0 ( 0 + ) and 1 =̃1 − 0 ( 0 + ) . Here = 0 and = 0 + are the left and right boundaries of the square barrier, respectively [see Figure 2(a)]. The transfer matrix can be obtained from the boundary conditions shown in (12) and (13) as with the elements, We can obtain the transmission through the barrier from the transfer matrix M, For normal incidence ( = 0), it is obvious that ( = 0) = 1 for any value of . Such perfect transmission phenomenon at normal incidence is a manifestation of Klein tunneling effect in pseudospin-1 systems. For pseudospin-1 EM waves, the perfect transmission can be interpreted as a consequence of impedance matching at normal incidence. At a finite incident angle, for values of 1 satisfying the relation 1 = ( ∈ integers), we will have ( ) = 1 due to Fabry-Perot resonances. More interestingly, when = 0 /2, we have 1 = − 0 and = − , and then obtain ( ) = 1 for any incident angle. Such all-angle transparency is called "super Klein tunneling", which can only happen for pseudospin-1 systems.
The super Klein tunneling effect was verified numerically using a realistic PC sandwich structure [47]. It can also be understood from the standpoint of impedance matching in EM wave theory. When = 0 /2, the permittivities and permeabilities in the background and potential barrier are equal in magnitude but opposite in sign [47], thus the impedance is matched for all incident angles, which in turn leads to all-angle unity transmission. In EM wave theory, it is known that two adjacent layers, in case that they have equal thickness and their permittivities and permeabilities are equal in magnitude but opposite in sign, can "optically cancel" each other in space, and are called "complementary materials" in the metamaterial literature [61]. Interestingly, one may take the pseudospin-1 photonic system as a candidate to realize "complementary materials". We note that "complementary materials" were extensively studied in the field of metamaterials which typically require the use of metallic inclusions to achieve the effective negative medium in order to cancel the optical path in a positive medium. However, the huge ohmic loss of metals at high frequencies makes it not so promising to realize them at optical frequencies.
Pseudospin-1 photonic system is hence a good candidate since the PCs are all dielectric so that they have low loss even at high frequencies. Recently, experimental realization of complementary materials has been demonstrated in the microwave regime using PCs composed of alumina ( = 8.1) rods in air [62,63]. It looks promising to push the working frequency to the optical regime in the near future.
We will show below how pseudospin-1 EM waves are scattered when they meet a 1D potential barrier in the case that deviates from 0 /2. To study the scattering properties in realistic PC systems, a photonic potential barrier can be constructed using a PC sandwich structure. As shown in Figure 2 Figure 3. The perfect transmission for the case of normal incidence ( = 0) is a manifestation of "Klein tunneling". At finite incident angles, we find some other transmission peaks with ( ) = 1 at angles satisfying the relation 1 = ( ∈ integers) as a consequence of Fabry-Perot resonances. When the incident angle is larger than the critical angle = sin −1 (( 0 − )/ ) ≅ 41 0 , total reflection occurs and no transmission is found. The result of ( ) from (21) for = 0.0436 / 2 is also shown in Figure 3 (blue solid line). Excellent agreement is found between the prediction from (21) and the numerical result for the realistic structure. For comparison, we plot the transmission in Figure 3 for Dirac electrons in graphene (green dotted line) through a  potential barrier with the same width = 40 2 . We take the energy and potential height 0 as /ℏV = /V and 0 /ℏV = /V , respectively, with V being the Fermi velocity of Dirac electrons. Note that Dirac electrons and pseudospin-1 EM waves have transmission peaks at the same angles due to Fabry-Perot effect. However, Dirac electrons with a pseudospin of 1/2 in graphene have much deeper transmission dips, which indicates that electrons experience much stronger anisotropic scattering by the 1D potential. It is found that the above difference is due to their different boundary conditions [47].

Supercollimation of Wave Packets in a Superlattice of
Photonic Crystals. The supercollimation here refers to a transport phenomenon that a wave packet can propagate a long distance in a 2D system while preserving its shape. Such a phenomenon was predicted theoretically in certain graphene superlattices consisting of two alternating layers with equal thickness due to anisotropic renormalization of the group velocity [13]. The occurrence of this phenomenon for pseudospin-1/2 system requires a specific condition on the product of the barrier height difference and layer thickness. However, we will see that it is more robust for pseudospin-1 EM waves as collimation can occur without any condition [47]. We consider here a pseudospin-1 superlattice formed by a Kronig-Penney type of photonic potential along the x-direction [see Figure 4(a)]. The thickness of each layer is taken to be d, the potential height is 0 , and lattice constant is L. For one unit cell of the periodic structure, the transfer matrix is the product of two parts, where M( ) is the part for the barrier layer as shown in (19) and P( − ) describes the part of transfer matrix in the background layer, Using (22) in conjunction with the Bloch theorem, we obtain the following dispersion relation for the superlattice [47], , is the Bloch wavevector, and is the wavevector component along y direction. When = 0 /2, it is found from (24) that the equifrequency curve is a straight line along the axis with = 0. It is clear that the group velocity in the y direction is zero for this equifrequency line. This result holds for all nonzero values of 0 . We plot the dispersion relation of (24) in Figure 4(b) for a superlattice with = 15 2 , = 30 2 , and 0 = 2 /15, where 2 is the lattice constant of PC2, and 0 is the potential shift between PC1 and PC2 (see Figure 2). A wedge structure can be clearly seen in Figure 4(b), exhibiting a dispersionless behavior along the direction even when ̸ = 0 /2. A wedge equation, = V | | + 0 /2 ( = ±1), can be found when | | << 0 /2V by expanding (24) in the vicinity of = 0 /2. It can be seen clearly from this equation that within the wedge-like dispersion structure, the 1D periodic modulation of photonic potential reduces the group velocity in y direction to zero. On the other side, it makes no impact on the one in x-direction. Similar to the case in some special graphene (pseudospin-1/2) superlattice [12,13], such wedge structure will lead to the supercollimation of pseudospin-1 wave packets in a superlattice of PCs due to the strong anisotropic renormalization of group velocity mentioned above, which means that a wave packet constructed in the frequency range of the wedge-like dispersion can be guided to propagate without distortion of shape along the potential modulation direction of the superlattice, independent of its initial direction of motion. It is worth noting here that for graphene, as a pseudospin-1/2 system, the wedge structure can only be achieved under certain special conditions that the periodic potentials 0 and the layer thickness d should satisfy the relation 0 = 2 ℏV / [12][13][14]. However, for a pseudospin-1 system, such a structure can be found for all nonzero values of 0 if each layer of the superlattice has the same thickness, i.e., = 2 [47].

Figures 4(c)-4(e)
show numerically the supercollimation phenomenon for a pseudospin-1 system in a Gaussian wave packet propagation simulation. In the simulation, we take the potential height as the photonic potential shift between PC1 and PC2 shown in Figure 2(b), i.e., 0 = 2 /15, and the width = 15 2 . The superlattice is connected to a lead which is a PC with a photonic potential 0 /2, and then a Gaussian wave packet is sent from the lead towards the superlattice. The initial wave packet in the form = 0 exp[−|r − r c | 2 / 2 0 + k c ⋅ (r − r c )] is shown in Figure 4(c). Here r c is the initial center of wave packet at a distance away from the left end of the superlattice,  Figure 4(e)]. When the superlattice is absent, the wave packet propagates along the direction of initial center wave vector marked by the red arrow and spreads sideway rapidly. However, when the superlattice is present, the Gaussian wave packet always propagates along the potential modulation direction ( -direction), irrespective of the initial direction of motion. Furthermore, it is found that at the incident angle = 45 0 , the wave packet is stretched and tilted, which can be attributed to the strong dependence of the reflections on the value of around = 45 0 when the wave packet enters the superlattice from the lead. The wave packet will remain undistorted after it enters the superlattice [47].
In conventional 1D disordered materials, all states must become localized due to the coherent backscattering effect [77][78][79][80][81]. However, for pseudospin-1 systems, waves propagating in the normal direction experience no backscattering from a disordered 1D potential and only acquire a random phase in the spatial wave function [47]. Such one-way transport behavior was first discovered in pseudospin-1/2 systems [15,73,82,83]. For pseudospin-1 EM waves considered in this review, the absence of backscattering in normal incidence can be understood from the classical wave impedance matching between any two neighboring layers, as we discussed in Section 4.2. Anderson localization can hence occur only for obliquely incident waves.
The systems considered here are 1D pseudospin-1 stacks of N uniform layers. Each layer is assumed to have a thickness d, but has a step-wise-constant random potential V(x). Here we take the Dirac-like point of the background medium as the zero of the photonic potential, i.e., V = 0. We consider the propagation of a plane wave through the multilayers at an oblique incident angle ( ̸ = 0) with a reduced frequency . The wave equation of pseudospin-1 EM waves in a 1D random potential can be written in a Hamiltonian form as (11), For convenience, we use the normalized frequency = /V and random potential ( ) = ( )/V in the following. The normalized potential in each layer fluctuates independently, with a uniform distribution in the interval [− , ], where is the random strength of the normalized potential. The localization length , defined as the reciprocal of the Lyapunov exponent , can be obtained from the following relation, where is the transmission coefficient and ⟨⟩ denotes ensemble averaging.
In Figure 5(a), we show the dependence of the localization length on the random strength at a particular normalized frequency = 0.02 for three different incident angles. Results are calculated by the transfer matrix method (TMM) and averaged over an ensemble of 4000 configurations. The sample size N is taken to be five times that of the localization length. We see that at small randomness, decays with the randomness according to a general form ∝ −2 as expected. However, as increases, drops suddenly to a minimum at a critical random strength = , regardless of incident angle, and rises immediately afterward. The sudden change of localization behavior near indicates the occurrence of some sharp transition between the two regions: < and > . In Figure 6, we plot the localization length as a function of incident angle in the two regions. We find that the asymptotic -dependence of indeed changes from ∝ sin −4 ( < ) to ∝ sin −2 ( > ) when moves across the critical random strength = . The V-shape turnaround of localization length occurring at a critical random strength and the change of the asymptotic behaviors in the -dependence are properties unique to pseudospin-1 systems and not found in any ordinary disordered systems. In ordinary materials, stronger disorder causes enhanced backscattering, and therefore a monotonic decrease of localization length with increasing disorder. In Figure 5(b), we show a similar plot as Figure 5(a) except that we now fix the incident angle at sin = 0.3 and plot as a function of random strength at three different values of . It is clearly seen that the same localization characteristics are found in three different regions of , i.e., < , ≈ , and > .     , where → is now a 2D Pauli vector in xy plane and E is the incident energy of pseudospin-1/2 quasiparticles. The localization length is plotted against the random strength for different incident angles and energies in Figures 5(c) and 5(d). Note here that for pseudospin-1/2 systems, we define the normalized energy and random potential ( ) as = /ℏV and ( ) = ( )/ℏV , respectively. It is found that there is also a localization length minimum in pseudospin-1/2 systems after which the localization length increases again, but the sharp change in found in pseudospin-1 systems is absent in pseudospin-1/2 transport. Instead, the localization length crosses over smoothly from a decreasing behavior at small random strength to an increasing one at large random strength. The minimum occurs at the random strength around a few . Furthermore, as can be seen in Figure 6, the -dependence of localization length in both decreasing and increasing regions shows an asymptotic behavior of ∝ sin −2 , in sharp contrast to pseudospin-1 behavior.
The above anomalous localization behaviors can be understood from the following analytic analysis. It turns out that (25) can be transformed into a scalar equation using a new coordinate u [52]:  where = ∫ 0 ( − ( )) and ( ) = ( ). In (27), the oblique angle appears in the scattering term. In the case of normal incidence, i.e., = 0, the scattering term vanishes, and (27) then describes wave propagation in a homogeneous medium and contains two general solutions, Ψ ∝ ± = exp[± ∫ 0 ( − ( )) ]. We can see that the random phase accumulation due to 1D random potential during the oneway transport is now absorbed in the new coordinate u. Similarly, a scalar wave equation for pseudospin-1/2 systems can be constructed according to pseudospin-1/2 Hamiltonian equation, Here is the ith interface of the layered structure in the u coordinate and v is the normalized random potential in ith layer. We can qualitatively understand from the scattering terms in (27) and (28) why the -dependence of in the two pseudospin systems behaves differently at small . For ordinary disordered materials, the localization length in disordered 1D systems has the same order of magnitude as the mean free path, and the latter is inversely proportional to the square of the scattering strength [79]. For small values of , the 2 dependence of the effective scattering potential in (27) leads to a ∝ −4 (or sin −4 ) behavior, whereas in addition to the 2 scattering term, (28) has another scattering term Ψ ∑ +1 =1 ( − ) located at all N+1 interfaces. This interface term dominates at small and lead to a ∝ 2 (or sin −2 ) behavior. The sudden drop of near = for pseudospin-1 systems can be understood from the diverging behavior of the scattering terms in (27) when | − ( )| < | | in some layers. In this case the waves become evanescent inside those layers. When exceeds the critical value = , the probability of having evanescent waves goes down as increases. Meanwhile, the scattering potentials in the propagating layers are weakened in general, as can be seen from the scattering terms in (27). Thus, the localization length goes up with increasing . However, such a sudden drop of localization length is smeared out by the additional interface scattering terms in (28) so that a smooth crossover from the localization length decreasing behavior to an increasing one is found for pseudospin-1/2 systems.
Analytical solutions obtained from the surface Green function (SGF) method [52] show that for pseudospin-1 systems, when > , the evanescent waves emerge and contribute to a sin 2 term, which dominates over the sin 4 term from propagating waves, resulting in the transition of the -dependence shown in Figure 6. For pseudospin-1/2 systems, both evanescent and propagating waves give a sin 2 dependence in , leading to the same asymptoticdependence for all s [52].
The 2D transport properties of pseudospin-1 waves also have unique features. The scattering of pseudospin-1 waves from a circularly symmetric potential barrier has been studied [53,54]. For the case of 2D PCs, such a potential barrier can be realized by changing the length scale in a selected circular domain. Many interesting scattering properties of pseudospin-1 waves from such a circular potential barrier have been predicted theoretically, including perfect caustics, revival resonance, isotropic transport, and superscattering [53,54].

Conclusions
We have reviewed the realization of pseudospin-1 physics in classical wave systems using 2D dielectric photonic crystals possessing Dirac-like dispersions at k=0. The physics of those systems near the Dirac-like point frequency can be described by an effective pseudospin-1 Hamiltonian which has three degrees of freedom, representing the projection of the real space wave function to monopole and dipole excitations. A photonic potential which can shift the Diraclike cone rigidly can be realized simply by changing the length scale of the dielectric PCs. Due to the unique Dirac-like conical dispersion, waves in pseudospin-1 systems interact with external potentials in unusual ways, leading to many exotic transport properties. We have also reviewed a number of unusual pseudospin-1 wave transport phenomena, including super Klein tunneling, supercollimation, and unconventional Anderson localization. Without the requirements of extremely low temperature and atom manipulation, the fabrication of photonic systems is much easier than that for ultracold-atom and electronic systems. With the recent rapid progress of modern microelectronic technologies, one can now achieve long enough coherent length to form complex junctions, especially in the microwave regime. Also, the absence of complex interactions between photons facilitates the advances of theories and makes the experimental observations and characterizations much easier. With the rapid development of the field of pseudospin-1 physics, there emerge lots of interesting phenomena, such as the superscattering of pseudospin-1 waves from a circular potential barrier [54], perfect caustics inside a circular barrier [53], current enhancement induced by the flat band in nonequilibrium transport of pseudospin-1 particles [55], and nonuniversal long-wavelength critical behavior of Anderson localization [51]. We expect that photonic crystals can be a good platform to test experimentally these theoretical predictions of pseudospin-1 physics. The demonstration of the above features of pseudospin-1 systems may provide us novel ways in controlling the propagation of pseudospin-1 particles, leading to possible applications in optical circuits and optical imaging beyond the diffraction limit.

Conflicts of Interest
The authors declare that they have no conflicts of interest.