The pseudo-magnetic field, an artificial synthetic gauge field, has attracted intense research interest in the classical wave system. The strong pseudo-magnetic field is realized in a two-dimensional photonic crystal (PhC) by introducing the uniaxial linear gradient deformation. The emergence of the pseudo-magnetic field leads to the quantization of Landau levels. The quantum-Hall-like edge states between adjacent Landau levels are observed in our designed experimental implementation. The combination of two reversed gradient PhCs gives rise to the spatially nonuniform pseudo-magnetic field. The propagation of the large-area edge state and the interesting phenomenon of the snake state induced by the nonuniform pseudo-magnetic field is experimentally demonstrated in a PhC heterostructure. This provides a good platform to manipulate the transport of electromagnetic waves and to design useful devices for information processing. |
1.IntroductionMany interesting phenomena emerge under the strong magnetic field, for example, the Landau level (LL) and quantum Hall effect.1,2 By artificially introducing the deformation of the graphene lattice by the strain, the Dirac points of the graphene generate a shift in the momentum space, resulting in an effective gauge field of vector potential.3–8 The behavior of electrons in this system is similar to that in a genuine magnetic field; therefore, the synthetic gauge field is called a pseudo-magnetic field (PMF). The PMF leads to quantization of LLs and quantum-Hall-like effects and possesses opposite signs at two inequivalent valleys, so that the time-reversal symmetry is not broken in a system, which provides a new method to control the motion of electrons. Unlike electronic systems, photons in a photonic crystal (PhC) are uncharged and do not directly couple to an external magnetic field; therefore, various exotic phenomena related to magnetic fields are not accessible in a PhC system (for example, the Aharonov–Bohm effect and Landau quantization). Many works have been studied on how to construct PMFs by using artificial metamaterial in various classical wave systems.9–23 Similar to the strained graphene, for example, the PMF can be generated by the uniaxial or triaxial deformation in sonic crystals10–13 and PhCs.14–21 Different from the perpendicular PMF generated by the shift of Dirac points, many works have been discovered that show the various in-plane PMFs can be obtained by spatially varying mass terms. This synthetic gauge field has been utilized to trap the photonic modes in the photonic Dirac cavity guide,24,25 guide the topological modes in the Dirac waveguide,26,27 and particularly realize the chiral LLs.28 Graphene applied a nonuniform perpendicular magnetic or PMF possesses peculiar properties, such as snake state29–35 and quantum confinement.36 The snake state describes the motion of the trajectory of electrons at the boundary where the magnetic field changes its direction. By analogy to this phenomenon, the snake state has been studied in an elastic wave system.37 However, the nonuniform PMF perpendicular PMF generated by the shift of degenerated Dirac points in the momentum space has not been experimentally demonstrated in a PhC system. In our work, we theoretically and experimentally realize the uniform and nonuniform PMFs in PhC heterostructures when the uniaxial linear gradient deformation is applied and the lattice symmetry is decreased. The edge state and the snake state in the gap between adjacent LLs are directly observed in the measured distributions of electric fields. 2.Photonic LLs and Edge StatesWe construct a two-dimensional PhC heterostructure in a spatially nonuniform PMF. The PMFs applied perpendicular to the PhC are constant within each of the three regions, where along the direction is positive, and gives the ratio of two PMFs on the upper and lower regions of a PhC heterostructure. We define the synthetic PMF as when a PhC has 21 layers in the direction. As shown in Fig. 1(a), the upper and lower experimental samples both have 11 layers corresponding to the PMF strength and the parameter . The upper and lower regions of the PhC sample satisfy the mirror symmetry [see the red and blue boxes in Fig. 1(a)]. The uniform PMF in our designed structure is realized in a gradient PhC with a linearly uniaxial deformation along the direction, which uses the same method as Ref. 12. As shown in Fig. 1(b), we first consider a triangular lattice PhC composed of metallic cylinders embedded in the air, in which the lattice constant is and the radius of the cylinder is 4 mm. In the calculated TM ( along the axis) band, degenerated Dirac cones appear at and points [see Fig. 1(c)]. When the filling ratio of the unit cell maintains unchanged and the metallic cylinders are deformed into an elliptical shape, which means the symmetry of the PhC is reduced from to , the Dirac cones shift along the direction in the Brillouin zone, but do not open the bandgap. The two axes of the metallic ellipse are represented by and . The shift of the Dirac cones denoted by has a linear relation with [see Fig. 1(d)] and gives rise to the vector potential .12 As we apply a linear gradient deformation in the direction and translational invariance along the direction, the of th row of the PhC structure is modulated linearly as where and are the length of elliptical pillars at the first and last layers of a gradient PhC, and are fixed as 3.3 and 5.3 mm, respectively, and the number of layers is modulated. Based on the relation between and PMF , the linear variation of causes the PMF along and axes for and valleys in PhCs. Without breaking the time-inversion symmetry, this method can produce similar phenomena caused by external real-magnetic fields. As shown in Fig. 1(e), the Dirac cones at the and valleys become discrete energy level plateaus, where is the LL index, is the photonic analog of the cyclotron frequency, and is the energy level of the Dirac point. The gradient PhC has layers of metallic pillars in the direction and obtains a uniform PMF in Fig. 1(e). The LLs are dispersive in this designed PhC structure, and flatter energy levels can be obtained by tuning the shift of Dirac cones. Therefore, the strong PMF in the gradient PhC leads to the quantization of photonic LLs.We then study the nonuniform PMF in a PhC heterostructure. As displayed in the schematic of Fig. 2(a), the upper and lower parts of the heterostructure have the same 21 layers but the reversed gradient; hence the PMFs of the two parts are antiparallel and the defined ratio at valley. The LLs are dispersive in this structure. Away from the and points, the LLs split into two bands along one direction in the momentum space due to the nonuniform PMF. The bands between adjacent LLs correspond to propagating edge states located at the middle region of the PhC heterostructure. Since the time-reversal symmetry is preserved, the edge states at distinct valleys propagate in opposite directions. The distributions of eigen-electric field at selected points for LL and edge states are shown in the right panels of Fig. 2(a). For the parallel and asymmetric case of corresponding to the combined 21- and 41-layer PhCs, LLs split into two bands along two directions, which become more complicated, and additional bands emerge between adjacent LLs, as shown in Fig. 2(b). We add three layers of metallic cylinders in the middle region as transition layers (PMF is zero) displayed in Fig. 2(c). Here, the number of layers is . The various LLs are deformed, and additional bands indicated by yellow regions correspond to the interesting electromagnetic transport, the so-called “snake state.” The parallel PMFs in the two PhCs cannot cause the snake states as displayed in Fig. 2(d). The PMF in the PhC leads to the transport phenomenon of edge states, which is similar to the quantum Hall edge state. The energy bands between the different LLs represent the dispersion of the edge states, as shown in Figs. 1(e) and 1(f). The localized edge state at the upper boundary gradually evolves into bulk states distributed in the middle of the structure ( LL) and finally becomes the edge state at the lower boundary near the valley. Here, we experimentally observe the propagation of the edge state. As shown in Fig. 3(a), we place an excited source marked by a red star in the left upper corner of a gradient PhC indicated by the red region. The experimental sample consists of 11 layers of aluminum pillars, and other parameters are the same with the simulations in Fig. 1. The right- (left-) propagating edge state are selectively excited on the top boundary with a positive (negative) group velocity in the simulations and experiments, as shown in Figs. 3(b) and 3(d), corresponding to 10.72 and 10.7 GHz, respectively, between and LLs. The chiral edge state between and LLs spatially distributes at the bottom boundary, shown in Fig. 3(c), which has a large-area distribution of the electric field, in contrast to the localization of edge state on the top side. Owing to the limitation of the scanning field in experiments, the measured area of the electric field is a part of the PhC sample, indicated by a white rectangle in the simulated distribution. An aluminum strip regarded as a perfect electric conductor at the microwave region is placed at the top boundary to prevent the edge state from scattering outside, and other boundaries are clad with absorbing materials. To describe the confinement of the edge state on the top boundary, the parameter is defined as38 where is the whole area of the PhCs, denoted by the green dotted box is the edge region as shown in Fig. 3(b), and is the electric field intensity. As shown in Fig. 3(e), the large imply well-localized microwave on the top boundary. The simulated and measured transmission spectra between input and output ports are displayed in Figs. 3(f) and 3(g). The frequency region with high transmission marked by yellow is consistent with the edge state dispersion between and LLs in Fig. 1(e). At LL with 10.78 GHz, the bulk modes are excited, in which electromagnetic waves spread into the central part.3.Transport of Interface States and Snake StatesIn the nonuniform PMFs constructed by two reversed gradient PhCs, the interface state distributes in the middle domain with a large area displayed in Fig. 2(a), in contrast to the edge states, which are spatially tightly confined on the top boundary of a gradient PhC with the uniform PMF in Fig. 3(b). Hence, it can carry a large amount of energy to design devices, for example, the photonic and acoustic large-area waveguides.39–42 To demonstrate the transport of the large-area interface state, we implement the experiments. The experimental sample is displayed in Fig. 3(a). The excited source is placed at the left of the PhCs marked by a blue star. The four boundaries of the experimental sample are all wrapped with absorbing material. The simulated and experimentally measured electric fields are shown in Figs. 4(a) and 4(b), which indicate that the interface state distributes in the middle domain corresponding to the dispersion between and LLs. The normalized electric field intensity at (560) mm along the direction is shown in Figs. 4(c) and 4(d). The electric field concentrates and distributes uniformly in the middle domain. At the interface of two domains with opposite magnetic fields, the cyclotron motions of electrons are driven by Lorenz force and present snake-like orbits. By analogy to this interesting transport, the snake state can be achieved in a PhC system. We design a structure composed of two reversed gradient PhCs. The experimental sample is shown in Fig. 1(a), in which the upper and lower parts have 11 layers. Fewer layers lead to a large PMF when the variation of the vector potential is fixed, which means the elliptical pillars at the top and bottom remain unchanged. The PMFs of two gradient PhCs have equal magnitudes but opposite directions. There are three transition layers of the circle pillars placed in the middle of the sample (yellow region), which have zero PMF. The larger transition layers lead to the wider frequency range of the snake state. The excited source is placed at the left of the PhC, as indicated by a red star in Fig. 1(a). The simulated and experimental distributions of the snake states are shown in Figs. 5(a) and 5(b). A disorder is introduced by randomly moving the position of three PhC pillars located at the bend of the wavy trajectory of the electromagnetic waves. The measured electric field demonstrates the propagation of the snake state is unaffected, as shown in Fig. 5(c). However, the transport of the snake state in PMF is fragile when the disorder or defect is introduced on the propagating path, which distinguishes it from the edge state in the quantum spin Hall effect.43,44 The intensity of electric fields at three positions along the direction, as shown in Fig. 5(e), experimentally demonstrates the bending property of the snake state, which agrees with the simulation in Fig. 5(d). 4.ConclusionThe PMF is achieved in a two-dimensional PhC when the lattice remains unchanged and metallic pillars are uniaxially deformed along the direction. A linear gradient PhC has a spatially uniform and strong PMF, which leads to the quantization of LLs and quantum Hall-like edge state confined on a boundary. The nonuniform PMF in a heterostructure composed of two reversed gradient PhCs has a large-area interface state in the middle domain. By adding a transition region with a zero PMF, the snake state arises from the deformed LLs. We implement the experiments to demonstrate the transport of the edge state and snake state. Our proposed structures have unique advantages in obtaining nonuniform PMFs, such as easy implementation, which provides a good platform to observe the manipulation of electromagnetic waves caused by the PMFs. Code and Data AvailabilityData underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request. AcknowledgmentsThis work was supported by the Graduate Innovation Program of China University of Mining and Technology (Grant No. 2023WLJCRCZL273), the Fundamental Research Funds for the Central Universities (Grant No. 2023ZDYQ11003), the China Postdoctoral Science Foundation (Grant No. 2023M743784), the State Key Laboratory of Millimeter Waves (Grant No. K202407), the Basic Research Program of Xuzhou (Grant No. KC22016), the Key Academic Discipline Project of China University of Mining and Technology (Grant No. 2022WLXK06), and the National Natural Science Foundation of China (Grant No. 12274315). ReferencesK. V. Klitzing, G. Dorda and M. Pepper,
“New method for high-accuracy determination of the fine-structure constant based on quantized Hall resistance,”
Phys. Rev. Lett., 45 494 https://doi.org/10.1103/PhysRevLett.45.494 PRLTAO 0031-9007
(1980).
Google Scholar
D. J. Thouless et al.,
“Quantized Hall conductance in a two-dimensional periodic potential,”
Phys. Rev. Lett., 49 405 https://doi.org/10.1103/PhysRevLett.49.405 PRLTAO 0031-9007
(1982).
Google Scholar
F. Guinea, M. I. Katsnelson and A. K. Geim,
“Energy gaps and a zero-field quantum Hall effect in graphene by strain engineering,”
Nat. Phys., 6 30 https://doi.org/10.1038/nphys1420 NPAHAX 1745-2473
(2010).
Google Scholar
N. Levy et al.,
“Strain-induced pseudo–magnetic fields greater than 300 Tesla in graphene nanobubbles,”
Science, 329 544 https://doi.org/10.1126/science.1191700 SCIEAS 0036-8075
(2010).
Google Scholar
G. Salerno et al.,
“How to directly observe Landau levels in driven-dissipative strained honeycomb lattices,”
2D Mater., 2 034015 https://doi.org/10.1088/2053-1583/2/3/034015
(2015).
Google Scholar
G. Salerno et al.,
“Propagating edge states in strained honeycomb lattices,”
Phys. Rev. B, 95 245418 https://doi.org/10.1103/PhysRevB.95.245418
(2017).
Google Scholar
É. Lantagne-Hurtubise, X. X. Zhang and M. Franz,
“Dispersive Landau levels and valley currents in strained graphene nanoribbons,”
Phys. Rev. B, 101 085423 https://doi.org/10.1103/PhysRevB.101.085423
(2020).
Google Scholar
B. L. Wu et al.,
“Transport property of inhomogeneous strained graphene,”
Chin. Phys. B, 30 030504 https://doi.org/10.1088/1674-1056/abe3e3 1674-1056
(2021).
Google Scholar
H. Abbaszadeh et al.,
“Sonic Landau levels and synthetic gauge fields in mechanical metamaterials,”
Phys. Rev. Lett., 119 195502 https://doi.org/10.1103/PhysRevLett.119.195502 PRLTAO 0031-9007
(2017).
Google Scholar
C. Brendel et al.,
“Pseudomagnetic fields for sound at the nanoscale,”
Proc. Natl. Acad. Sci. U. S. A., 114 3390 https://doi.org/10.1073/pnas.1615503114
(2017).
Google Scholar
Z. Yang et al.,
“Strain-induced gauge field and Landau levels in acoustic structures,”
Phys. Rev. Lett., 118 194301 https://doi.org/10.1103/PhysRevLett.118.194301 PRLTAO 0031-9007
(2017).
Google Scholar
X. Wen et al.,
“Acoustic Landau quantization and quantum-Hall-like edge states,”
Nat. Phys., 15 352
–356 https://doi.org/10.1038/s41567-019-0446-3 NPAHAX 1745-2473
(2019).
Google Scholar
J. Luo et al.,
“Pseudomagnetic fields and Landau levels for out-of-plane elastic waves in gradient snowflake-shaped crystals,”
Phys. Lett. A, 383 125974 https://doi.org/10.1016/j.physleta.2019.125974 PYLAAG 0375-9601
(2019)).
Google Scholar
M. C. Rechtsman et al.,
“Strain-induced pseudomagnetic field and photonic Landau levels in dielectric structures,”
Nat. Photonics, 7 153
–158 https://doi.org/10.1038/nphoton.2012.302 NPAHBY 1749-4885
(2013).
Google Scholar
H. Schomerus and N. Y. Halpern,
“Parity anomaly and Landau-level lasing in strained photonic honeycomb lattices,”
Phys. Rev. Lett., 110 013903 https://doi.org/10.1103/PhysRevLett.110.013903 PRLTAO 0031-9007
(2013).
Google Scholar
O. Jamadi et al.,
“Direct observation of photonic Landau levels and helical edge states in strained honeycomb lattices,”
Light Sci. Appl., 9 144 https://doi.org/10.1038/s41377-020-00377-6
(2020).
Google Scholar
M. Bellec et al.,
“Observation of supersymmetric pseudo-Landau levels in strained microwave graphene,”
Light Sci. Appl., 9 146 https://doi.org/10.1038/s41377-020-00351-2
(2020).
Google Scholar
C. R Mann, S. A. R. Horsley and E. Mariani,
“Tunable pseudo-magnetic fields for polaritons in strained metasurfaces,”
Nat. Photonics, 14 669
–674 https://doi.org/10.1038/s41566-020-0688-8 NPAHBY 1749-4885
(2020).
Google Scholar
J. Guglielmon, M. C. Rechtsman and M. I. Weinstein,
“Landau levels in strained two-dimensional photonic crystals,”
Phys. Rev. A, 103 013505 https://doi.org/10.1103/PhysRevA.103.013505
(2021).
Google Scholar
Z. T. Huang et al.,
“Strain-tunable synthetic gauge fields in topological photonic graphene,”
(2021). Google Scholar
F. S. Deng et al.,
“Valley-dependent beams controlled by pseudomagnetic field in distorted photonic graphene,”
Opt. Lett., 40 3380
–3383 https://doi.org/10.1364/OL.40.003380 OPLEDP 0146-9592
(2015).
Google Scholar
W. Wang et al.,
“Moiré fringe induced gauge field in photonics,”
Phys. Rev. Lett., 125 203901 https://doi.org/10.1103/PhysRevLett.125.203901 PRLTAO 0031-9007
(2020).
Google Scholar
H. T. Teo et al.,
“Pseudomagnetic suppression of non-Hermitian skin effect,”
(2023). Google Scholar
R. Barczyk et al.,
“Interplay of leakage radiation and protection in topological photonic crystal cavities,”
Laser Photonics Rev., 16 200071 https://doi.org/10.1002/lpor.202200071
(2022).
Google Scholar
K. Chen et al.,
“Photonic Dirac cavities with spatially varying mass term,”
Sci. Adv., 9 eabq4243 https://doi.org/10.1126/sciadv.abq4243 STAMCV 1468-6996
(2023).
Google Scholar
S. Kiriushechkina et al.,
“Spin-dependent properties of optical modes guided by adiabatic trapping potentials in photonic Dirac metasurfaces,”
Nat. Nanotechnol., 18 875
–881 https://doi.org/10.1038/s41565-023-01380-9 NNAABX 1748-3387
(2023).
Google Scholar
A. Vakulenko et al.,
“Adiabatic topological photonic interfaces,”
Nat. Commun., 14 4629 https://doi.org/10.1038/s41467-023-40238-5 NCAOBW 2041-1723
(2023).
Google Scholar
H. Jia et al.,
“Experimental realization of chiral Landau levels in two-dimensional Dirac cone systems with inhomogeneous effective mass,”
Light Sci. Appl., 12 165 https://doi.org/10.1038/s41377-023-01209-z
(2023).
Google Scholar
L. Oroszlány et al.,
“Theory of snake states in graphene,”
Phys. Rev. B, 77 081403 https://doi.org/10.1103/PhysRevB.77.081403
(2008).
Google Scholar
T. K. Ghosh et al.,
“Conductance quantization and snake states in graphene magnetic waveguides,”
Phys. Rev. B, 77 081404 https://doi.org/10.1103/PhysRevB.77.081404
(2008).
Google Scholar
J. R. Williams and C. M. Marcus,
“Snake states along graphene p-n junctions,”
Phys. Rev. Lett., 107 046602 https://doi.org/10.1103/PhysRevLett.107.046602 PRLTAO 0031-9007
(2011).
Google Scholar
T. Taychatanapat et al.,
“Conductance oscillations induced by ballistic snake states in a graphene heterojunction,”
Nat. Commun., 6 6093 https://doi.org/10.1038/ncomms7093 NCAOBW 2041-1723
(2015).
Google Scholar
L. Cohnitz et al.,
“Interaction-induced conductance from zero modes in a clean magnetic graphene waveguide,”
Phys. Rev. B, 92 085422 https://doi.org/10.1103/PhysRevB.92.085422
(2015).
Google Scholar
P. Rickhaus et al.,
“Snake trajectories in ultraclean graphene p–n junctions,”
Nat. Commun., 6 6470 https://doi.org/10.1038/ncomms7470 NCAOBW 2041-1723
(2015).
Google Scholar
C. Zuo et al.,
“Reverse strain-induced snake states in graphene nanoribbons,”
Phys. Rev. B, 105 195420 https://doi.org/10.1103/PhysRevB.105.195420
(2022).
Google Scholar
Y. N. Ren et al.,
“Magnetic field-tunable valley-contrasting pseudomagnetic confinement in graphene,”
Phys. Rev. Lett., 129 076802 https://doi.org/10.1103/PhysRevLett.129.076802 PRLTAO 0031-9007
(2022).
Google Scholar
M. Yan et al.,
“Pseudomagnetic fields enabled manipulation of on-chip elastic waves,”
Phys. Rev. Lett., 127 136401 https://doi.org/10.1103/PhysRevLett.127.136401 PRLTAO 0031-9007
(2021).
Google Scholar
B. Yang et al.,
“Details of the topological state transition induced by gradually increased disorder in photonic Chern insulators,”
Opt. Express, 28 31487
–31498 https://doi.org/10.1364/OE.405820 OPEXFF 1094-4087
(2020).
Google Scholar
M. Wang et al.,
“Topological one-way large-area waveguide states in magnetic photonic crystals,”
Phys. Rev. Lett., 126 067401 https://doi.org/10.1103/PhysRevLett.126.067401
(2021).
Google Scholar
M. Wang et al.,
“Valley-locked waveguide transport in acoustic heterostructures,”
Nat. Commun., 11 3000 https://doi.org/10.1038/s41467-020-16843-z NCAOBW 2041-1723
(2020).
Google Scholar
J. Q. Wang et al.,
“Extended topological valley-locked surface acoustic waves,”
Nat. Commun., 13 1324 https://doi.org/10.1038/s41467-022-29019-8 NCAOBW 2041-1723
(2022).
Google Scholar
Q. L. Chen et al.,
“Photonic topological valley-locked waveguides,”
ACS Photonics, 8 1400
–1406 https://doi.org/10.1021/acsphotonics.1c00029
(2021).
Google Scholar
M. Z. Hasan and C. L. Kane,
“Colloquium: topological insulators,”
Rev. Mod. Phys., 82 3045 https://doi.org/10.1103/RevModPhys.82.3045 RMPHAT 0034-6861
(2010).
Google Scholar
B. K. Alexander and S. Gennady,
“Two-dimensional topological photonics,”
Nat. Photonics, 11 763
–773 https://doi.org/10.1038/s41566-017-0048-5 NPAHBY 1749-4885
(2017).
Google Scholar
BiographyBin Yang is currently a master’s student in the School of Materials and Physics at China University of Mining and Technology. His research focuses on photonic crystals and topological insulators under the guidance of Associate Professor Yuting Yang. He strives to advance the theory and explore novel optical phenomena in these fields. Yuting Yang received her PhD from Soochow University, Suzhou, China, in 2019. During 2018 to 2019, she worked as a visiting PhD student at Nanyang Technological University. She is currently an associate professor in China University of Mining and Technology. Her research interest focuses on topological insulators, photonic crystals, and non-Hermitian topological electric circuits. |