We explore the distinctive properties associated with a type-II Dirac point in a simply structured phononic crystal with a lattice deformation. This type-II Dirac point emerges at the Brillouin zone boundary, resulting from the lifting of two degenerate bands and featuring a conical-like Fermi surface in the equi-frequency curve. A practical implementation of such a phononic crystal is achieved with LEGO bricks. Upon introducing a periodic parity-time (PT) symmetric non-Hermitian perturbation, the phononic crystal undergoes a transition from PT-symmetric phase to PT-broken phase, causing the deformation of type-II Dirac point into an oval of exceptional points in the band structure. Based on the eigenmodes of the type-II Dirac point, a k p perturbation theory can be used to characterize these systems before and after the phase transition. Using a scattering matrix, we analyze the symmetric and broken phases and demonstrate that broadband unidirectional transparency and a coherent perfect absorber and laser can be realized with such a phononic crystal slab.

When the effective mass of a particle approaches zero, a linear dispersion emerges from the solution of the Dirac equation.1 Over the past decade, extensive research has been conducted on the Dirac conical dispersion in electromagnetic,2–4 acoustic,5–8 and mechanical wave systems.9,10 Numerous interesting phenomena, including zero-reflective-index properties11–18 and topological transitions,2–8,19–21 are directly linked to the distinctive features of Dirac cones. When certain non-Hermitian perturbations are introduced, the Dirac point can be deformed and spawn a ring of exceptional points (EPs),22–29 which is associated with unique opportunities such as unidirectional transparency,30,31 super-prism effect,25 and single-mode lasers.32,33

Dirac points can be categorized into two classes: type-I Dirac points exhibit a point-like Fermi surface and vanishing density of states (DOS) at the frequency of the Dirac point, while type-II Dirac points feature strongly titled dispersion, conical-like Fermi surfaces, and non-vanishing DOS at the frequency of Dirac point.34–38 Although less prevalent than its higher-dimensional counterpart, the type-II Weyl point,39–46 the manifestation of a type-II Dirac point, serves to enhance comprehension and facilitate research in the realm of topological semimetals with strongly tilted dispersion.

In this study, we present a straightforward realization of a type-II Dirac point in a two-dimensional phononic crystal (PC). In contrast to previous implementations of type-II linear dispersions, which are predominantly reliant on the intricate design of unit cell structures,34–46 our approach leverages lattice symmetry and band folding. We construct a rectangular lattice of unit cells by deforming a supercell consisting of two copies of square unit cells. This configuration induces a type-II Dirac point at the Brillouin zone boundary, accompanied by a conical-like Fermi surface in the equi-frequency curve (EFC). The PC is experimentally realized using simple LEGO bricks. When PT-symmetric non-Hermitian periodicity is introduced, the type-II Dirac point is deformed into a titled oval of EPs. Meanwhile, the response of the PC to the incident waves from the left and right sides becomes disparate and the scattering matrix of the PC undergoes a phase transition. In contrast to one-dimensional PT-symmetric systems, which yield unidirectional transparency at a single frequency, our phase transition results in two curves of EPs in the eigenvalues of the scattering matrix, corresponding to broadband unidirectional transparency. Within the titled oval of EPs, the PC with PT-symmetric non-Hermitian periodicity exhibits characteristics of both a coherent perfect absorber (CPA) and a laser oscillator when there is an appropriate phase difference in the waves incident from different sides.

We start with a two-dimensional PC depicted in the inset at the upper left corner of Fig. 1(a). The PC has a supercell comprising two solid rods, which are treated as sound-hard bodies in the air. The lattice constant is 2a in the x direction and a in the y direction. Both rods have a square cross section with side length L = a/3. We calculate the band structure of this supercell along the boundaries of the rectangular irreducible Brillouin zone, using the finite-element package COMSOL Multiphysics, and plot the results in Fig. 1(a). The dimensionless frequency, 2πc/a, is adopted, where c is the speed of sound in air. The band structure in Fig. 1(a) can be interpreted as the folded band structure of the square unit cell along the kx = π/(2a) line, resulting in a doubly degenerate line along the XM boundary. We then modify the supercell by moving the two square rods closer until they touch, as shown in the inset of Fig. 1(b). The original double degeneracy along the XM boundary is lifted, and the bands split into two linear branches with distinct slopes, intersecting at k x , k y = π 2 a , 0.505 π a . In other words, a type-II Dirac point is induced below the sound cone, represented by the blue lines in the band structure.

FIG. 1.

Induction of the type-II Dirac point by deforming a two-dimensional phononic crystal. (a) Band structure of a supercell, which is comprised of two sound-hard square rods (blue) in air (gray) (inset). The blue line represents the sound cone. (b) When the two square rods are shifted together, a type-II Dirac point (black dot) is induced at the XM boundary of the Brillouin zone. (c) and (d) Band structure near the type-II Dirac point along the x-direction and y-direction, respectively. The red circles are the results from numerical simulation using COMSOL, and the black curves are obtained from k p analysis. (e) A zoom-in plot of the oblique conical dispersion surfaces near the type-II Dirac point. (f) Equi-frequency curves at different frequencies near the midpoint of the XM boundary.

FIG. 1.

Induction of the type-II Dirac point by deforming a two-dimensional phononic crystal. (a) Band structure of a supercell, which is comprised of two sound-hard square rods (blue) in air (gray) (inset). The blue line represents the sound cone. (b) When the two square rods are shifted together, a type-II Dirac point (black dot) is induced at the XM boundary of the Brillouin zone. (c) and (d) Band structure near the type-II Dirac point along the x-direction and y-direction, respectively. The red circles are the results from numerical simulation using COMSOL, and the black curves are obtained from k p analysis. (e) A zoom-in plot of the oblique conical dispersion surfaces near the type-II Dirac point. (f) Equi-frequency curves at different frequencies near the midpoint of the XM boundary.

Close modal
To analytically investigate the titled linear dispersion in the vicinity of the type-II Dirac point shown in Fig. 1(b), we adopt a k p method based on perturbation, which gives the following 2 × 2 effective Hamiltonian:5 
(1)
where Δ k = k k 0 , k 0 is the wave vector at the Dirac point, Ψ n k 0 ( r ) is the Bloch wave function of the nth band at the Dirac point, ρ r ( r ) is the normalized mass density, Ω is the area of the unit cell, and p l j (l, j = 1 or 2) is the coefficient of linear term in k, which takes the following values calculated from the Bloch wave function obtained from COMSOL:
(2)
Subsequently, we can obtain the slopes of the linear dispersions in x and y directions as
(3)
where s x = p 11 x 2 + p 12 x 2 , s y 1 = p 11 y 2 + p 12 y 2 , and s y 2 = p 12 y 2 + p 22 y 2 . We plot the linear dispersions using the slopes from Eq. (3) with black lines in Figs. 1(c) and 1(d), and found that they, indeed, overlap with the band dispersions near the type-II Dirac point calculated from COMSOL.

Enlarged views of the dispersion surface and the EFCs at different frequencies near the type-II Dirac point are plotted in Figs. 1(e) and 1(f), respectively. A tilted conical dispersion is clearly observed. Above or below the frequency of the Dirac point ω D = 0.33 2 π c / a , the EFC shows two untouched quadratic branches at the Brillouin zone boundary. At the frequency ωD, the EFC exhibits a cross at the Dirac point with conical-like Fermi surfaces. Such EFCs distinguish the type-II Dirac point from type-I Dirac points, which have a point-like EFC at ωD.

Due to the simple structure of the PC, we can construct it using LEGO bricks. As shown in Fig. 2(a), we assemble the LEGO bricks on a thin LEGO baseplate, featuring a periodic square array of fixing points for the bricks. Each unit cell measures 48 × 24 mm2, with two 1 × 1 LEGO bricks affixed at the center to form a scatterer with a rectangular cross section of 16 × 8 mm2 and a height of 9.6 mm. The small cylindrical studs on the baseplate and the top of the bricks do not impact our experimental results due to their small size. Since the type-II Dirac point is located outside the sound cone, the associated bands can manifest as acoustic surface waves near the PC–air interface (supplementary material). This proves advantageous as the real-space field distribution can be conveniently measured by raster-scanning the top of the PC with a microphone, as shown in Fig. 2(b). Subsequently, a two-dimensional Fourier transform is applied to convert real-space fields to k-space to obtain the EFC. In our experiment, a loudspeaker with a small opening port is placed at the edge of the PC as the sound source. In Fig. 2(c), two arcs are observed in the EFC at the frequency of Dirac point 4800 Hz. The brighter arc corresponds to the forward propagation wave, and the dimmer arc arises from the backscattering wave components in the PC. These arcs represent the different branches of phononic bands in the EFC. The numerically calculated EFC at the frequency of the Dirac point is plotted in Fig. 2(d), showing excellent agreement with the experimental results. If we position the source on the other edge of the PC, the resulting EFC will be the mirror image of the original EFC. Together, they constitute the complete EFC at the Dirac point’s frequency, with the arc crossings representing the type-II Dirac points.

FIG. 2.

Experimental realization of type-II Dirac point. (a) The PC is composed by LEGO bricks. (b) A photograph of the experimental setup. The top cladding plate is removed to show the PC. LEGO bricks of different colors are identical in their properties. The equi-frequency curve at the frequency of type-II dirac point from (c) the Fourier transform of scanned acoustic field and (d) the numerical results in COMSOL are plotted.

FIG. 2.

Experimental realization of type-II Dirac point. (a) The PC is composed by LEGO bricks. (b) A photograph of the experimental setup. The top cladding plate is removed to show the PC. LEGO bricks of different colors are identical in their properties. The equi-frequency curve at the frequency of type-II dirac point from (c) the Fourier transform of scanned acoustic field and (d) the numerical results in COMSOL are plotted.

Close modal
A Dirac point in a Hermitian system can be deformed into a ring of EPs under certain non-Hermitian perturbations.24–26 An intriguing question arises: what happens to a type-II Dirac point when non-Hermiticity is taken into account? As illustrated in the inset of Fig. 3(a), we replace the right and left parts of air in the unit cell with gain (colored blue) and loss (colored red) materials, characterized by the speed of sound v g = 343 10 i m / s and v l = 343 + 10 i m / s , respectively. The mass density of both gain and loss materials remains 1.21 kg/m3. The resulting band diagram along the XM boundary is presented in Fig. 3(a). The black solid line represents the real part of the band structure, depicting a titled degenerate region, while the blue dashed line delineates a split region in the imaginary part of the band structure. The non-Hermitian perturbation deforms the type-II Dirac point into a titled oval of EPs (supplementary material), marking the interface between PT-exact phase and PT-broken phase of the band diagram. At the XM boundary of Brillouin zone, there exist two EPs, which are labeled as EP1 and EP2 at different frequencies ω 1 = 0.313 2 π c / a and ω 2 = 0.347 2 π c / a . Near the original type-II Dirac point, the effective k p Hamiltonian with non-Hermitian perturbations is expressed as (supplementary material)
(4)
where β i ( r ) is the imaginary part of the relative compressibility. The associated eigenstates of H coalesce at EPs. In Fig. 3(b), we demonstrate the consistency between the dispersion calculated from the Hamiltonian with non-Hermitian perturbation and the finite-element simulation near the original type-II Dirac point. The semi-minor and semi-major axes of the oval of EPs are along the kx and ky directions, respectively. This is clearly shown in Fig. 3(c), illustrating the real part of the dispersion surface near the original type-II Dirac point. To provide a comprehensive perspective of the EPs oval, the corresponding imaginary part of the dispersion surface is presented in Fig. 3(d), complementing the information conveyed in Fig. 3(c). The yellow and green lines shown in Fig. 3(c), intersecting the EPs oval at EP1 and EP2 in Fig. 3(a), represent eigenstates that can perfectly couple with an externally incident plane wave from the left or right side in the air. Such coupling is facilitated by the surface impedance matching (supplementary material).
FIG. 3.

Oval of EPs when PT-symmetric non-Hermitian periodicity is considered in the PC. (a) The real part (black solid line) and imaginary part (blue dashed line) of the band structure at the XM boundary of the Brillouin zone. The red dashed lines are the original type-II dirac conical dispersions. (b) The real part (solid lines and red circles) and imaginary part (dashed lines and blue circles) of the band structure near the oblique exceptional oval, calculated from simulations (circles) or k p analysis (lines). (c) The real part of the dispersion surface in the vicinity of the original Dirac point. The yellow and green lines are the eigenstates that can couple perfectly with incident plane waves in the air from the left and right sides, respectively. (d) The imaginary part of the dispersion surface in the vicinity of the original Dirac point.

FIG. 3.

Oval of EPs when PT-symmetric non-Hermitian periodicity is considered in the PC. (a) The real part (black solid line) and imaginary part (blue dashed line) of the band structure at the XM boundary of the Brillouin zone. The red dashed lines are the original type-II dirac conical dispersions. (b) The real part (solid lines and red circles) and imaginary part (dashed lines and blue circles) of the band structure near the oblique exceptional oval, calculated from simulations (circles) or k p analysis (lines). (c) The real part of the dispersion surface in the vicinity of the original Dirac point. The yellow and green lines are the eigenstates that can couple perfectly with incident plane waves in the air from the left and right sides, respectively. (d) The imaginary part of the dispersion surface in the vicinity of the original Dirac point.

Close modal
As a result of complex interference inside PT-symmetric systems, there is an interesting phenomenon called “unidirectional transparency,” which has attracted broad interest.30,31 Here, we consider a PC slab consisting of 14 layers of unit cells. Near the oval of EPs in the band diagram, the response of the PC is scanned as a function of both frequency and incident angle. The analytic calculations are based on the following scattering matrix:30,31
(5)
where t is the transmission coefficient and rl and rg are the reflection coefficients from loss side and gain side of the PC slab, respectively.31, rl and rg are unequal in magnitude but can differ in phase by only 0 or π.30 The eigenvalues of S1 are given by
(6)
and are plotted in Figs. 4(a) and 4(b), respectively. When rgrl < 0 (rgrl > 0), both eigenvalues are complex (real). When rg = 0 or rl = 0, EPs of S1 occurs, with only one eigenvalue λ1 = λ2 = t and one eigenvector. Interestingly, the trajectory of EPs (dashed line) corresponds to the trajectory of eigenstates in the dispersion surface whose surface impedance is perfectly matched with air (supplementary material),47,48 implying broadband unidirectional transparencies. The open circles shown in Figs. 4(a) and 4(b) correspond to PT-broken phase in Fig. 3(a). The two end points, i.e., points A and B in Fig. 4(a), correspond to EP1 and EP2, which are also EPs of S1. We note that there are stripes between the trajectories of EPs in Fig. 4(a) at which λ1 is nearly 1. They are not EPs. These stripes are related to the finite size effect of PCs and the definition of scattering matrix. When the scattering matrix is defined as S1, its eigenvalues will pass through 1 repeatedly.30 At these points, there exist two orthogonal eigenvectors, while the eigenvectors at EPs coalesce. Figure 4(c) illustrates the acoustic pressure field distributions when a plane wave with frequency ω 1 = 0.313 2 π c / a is incident at an angle θ1 = 45.1° corresponding to the EP1 in Fig. 3(a). When the plane wave is incident from the loss side, the system achieves total transmission with vanishing reflection. Conversely, when the incident wave comes from the gain side, the transmission remains at unity, but the reflection is significantly amplified (about 298%). In Fig. 4(d), the acoustic pressure field distributions are presented for incident waves at the frequency ω 2 = 0.347 2 π c / a and incident angle θ2 = 54.5°, associated with EP2 in Fig. 3(a). In contrast to the previous scenario, while the transmission still reaches unity irrespective of whether the wave is incident from the gain or loss side, the reflection is zero when the incident wave comes from the gain side but high (about 295%) when the incident wave comes from the loss side. In Figs. 4(e) and 4(f), we present the acoustic pressure field distribution corresponding to the excitation of incident waves at specific points marked by an open square (point C) and open diamond (point D) in Fig. 4(a). The frequency and angle of incidence are ω = 0.344 2 π c / a and θ = 44.3° for point C and ω = 0.321 2 π c / a and θ = 56.6° for point D. In both cases, unity transmission is consistently achieved, while the reflection is zero when the incident wave comes from the loss side for point C and the gain side for point D. 51% and 64% reflection are observed if the wave incident from gain (point C) and loss (point D) side, respectively. These findings highlight the presence of broadband unidirectional transparency at the EP curves of the S matrix.
FIG. 4.

(a) and (b) Eigenvalues λ 1 and λ 2 of scattering matrix S1 near the PT-broken phase. Open circles correspond to the PT-broken phase shown in Fig. 3(a). (c)–(f) The pressure field distributions under the excitation of plane wave incident with the profile associated with A–D points shown in (a). Unidirectional transparency is found in (c) and (e) if incident from the loss side and in (d) and (f) if incident from the gain side.

FIG. 4.

(a) and (b) Eigenvalues λ 1 and λ 2 of scattering matrix S1 near the PT-broken phase. Open circles correspond to the PT-broken phase shown in Fig. 3(a). (c)–(f) The pressure field distributions under the excitation of plane wave incident with the profile associated with A–D points shown in (a). Unidirectional transparency is found in (c) and (e) if incident from the loss side and in (d) and (f) if incident from the gain side.

Close modal
In the PT-broken phase, the PC can behave both as a laser oscillator and as a CPA, which is the so-called “CPA-laser.”49–51 This CPA-lasing solution can be predicted by another version of the scattering matrix:
(7)

The difference between S1 and S2 is the way the vector of the output field is written.30,52 For the same PC slab discussed earlier, the eigenvalues λ 1 and λ 2 of S2 are plotted in Figs. 5(a) and 5(b), exhibiting a unimodular ( λ 1 = λ 2 = 1 ) behavior in the PT-exact phase and inverse ( λ 1 = 1 / λ 2 ) in the PT-broken phase.30,52 At the CPA-lasing point, λ 2 goes to the maximum, while λ 1 goes to zero. Figure 5(c) exhibits the simulated results of such a typical CPA-laser. The incident wave comes from both sides, with frequency ω = 0.332 2 π c / a and incident angle θ = 50.42°. When the phase difference of the incoming waves is Δφ = φL − φG = 0.5π, the coherent interference of incident waves confines the acoustic pressure field in the loss regions of the PC, resulting in perfect absorption. On the contrary, when the phase difference of the incoming waves is Δφ = 1.5π, the acoustic pressure field is confined in the gain regions of the PC, and the outgoing waves are strongly amplified. The frequency and incident angle of CPA-lasing point correspond to the position of the original Dirac point.

FIG. 5.

CPA-lasing effect. (a) and (b) The eigenvalues λ 1 and λ 2 of the scattering matrix S2. (c) At the CPA-lasing point, the PC can behave as a perfect absorber (left) or a laser (right) with appropriate phase difference in the incident wave from different sides. Colors represent the amplitude of the acoustic pressure field.

FIG. 5.

CPA-lasing effect. (a) and (b) The eigenvalues λ 1 and λ 2 of the scattering matrix S2. (c) At the CPA-lasing point, the PC can behave as a perfect absorber (left) or a laser (right) with appropriate phase difference in the incident wave from different sides. Colors represent the amplitude of the acoustic pressure field.

Close modal

In summary, we have theoretically designed and experimentally demonstrated the existence of an acoustic type-II Dirac point in a two-dimensional PC. Through the implementation of a rectangular supercell, a band-folding mechanism was employed to generate doubly degenerate bands along the Brillouin zone boundary. By deforming the lattice, the band degeneracy is lifted, resulting in the emergence of a strongly tilted Dirac point. Notably, when introducing a PT-symmetric non-Hermitian perturbation, the band dispersion exhibits a titled oval of EPs, indicating the transition from the PT-symmetric phase to the PT-broken phase. The band dispersion of the PC can be accurately characterized by the k p method. Interestingly, there exists broadband unidirectional transparency along two lines on the dispersion surfaces where the eigenstates can perfectly couple with externally incident plane waves in the air. In the PT-broken phase, the PC can function as both a laser oscillator and a CPA under specific conditions of phase difference between the left and right sides. Our proposed approach establishes a versatile platform for future investigations on the characteristics and applications of type-II Dirac points in photonic and phononic crystals.

See the supplementary material for the band structure near the type-II Dirac point in the phononic crystal, analysis of parity-time symmetric non-Hermitian perturbations in our phononic crystals, acoustic pressure field distributions of eigenstates in the PT-exact phase and the PT-broken phase, surface impedance of eigenstates associated with EPs of the scattering matrix, reverse direction of unidirectional transparency, and amplification near the CPA-lasing point.

This work was partially supported by the King Abdullah University of Science and Technology (KAUST) Office of Sponsored Research (OSR) under Award No. OSR-2020-CRG9-4374 and KAUST Baseline Research Fund No. BAS/1/1626-01-01. C.X. also acknowledges the support from the Jiangsu Specially Appointed Professor Program (Grant No. 164080H00244). G.M. was supported by the Hong Kong Research Grants Council (Grant Nos. RFS2223-2S01, 12302420, and 12301822). J.M. was supported by Guangdong Basic and Applied Basic Research Foundation (Grant No. 2021A1515010322).

The authors have no conflicts to disclose.

Changqing Xu: Conceptualization (equal); Data curation (equal); Formal analysis (equal); Funding acquisition (equal); Investigation (equal); Methodology (equal); Writing – original draft (lead); Writing – review & editing (equal). Jun Mei: Conceptualization (equal); Formal analysis (equal); Funding acquisition (lead); Investigation (equal); Methodology (equal); Project administration (equal); Writing – review & editing (supporting). Guancong Ma: Conceptualization (equal); Formal analysis (equal); Funding acquisition (lead); Investigation (equal); Methodology (equal); Project administration (equal); Supervision (equal); Writing – review & editing (supporting). Ying Wu: Conceptualization (equal); Formal analysis (equal); Funding acquisition (lead); Investigation (equal); Methodology (equal); Project administration (equal); Supervision (lead); Writing – review & editing (lead).

The data that support the findings of this study are available from the corresponding author upon reasonable request.

1.
P. A. M.
Dirac
, “
The quantum theory of the electron
,”
Proc. R. Soc. A
117
,
610
(
1928
).
2.
S. A.
Skirlo
,
L.
Lu
, and
M.
Soljačić
, “
Multimode one-way waveguides of large Chern numbers
,”
Phys. Rev. Lett.
113
,
113904
(
2014
).
3.
L.
Lu
,
J. D.
Joannopoulos
, and
M.
Soljačić
, “
Topological photonics
,”
Nat. Photonics
8
,
821
829
(
2014
).
4.
A.
Slobozhanyuk
,
S. H.
Mousavi
,
X.
Ni
,
D.
Smirnova
,
Y. S.
Kivshar
, and
A. B.
Khanikaev
, “
Three-dimensional all-dielectric photonic topological insulator
,”
Nat. Photonics
11
,
130
136
(
2017
).
5.
J.
Mei
,
Y.
Wu
,
C. T.
Chan
, and
Z.-Q.
Zhang
, “
First-principles study of dirac and dirac-like cones in phononic and photonic crystals
,”
Phys. Rev. B
86
,
035141
(
2012
).
6.
S.-Y.
Yu
,
X.-C.
Sun
,
X.
Ni
,
Q.
Wang
,
X.-J.
Yan
,
C.
He
,
X.-P.
Liu
,
L.
Feng
,
M.-H.
Lu
, and
Y.-F.
Chen
, “
Surface phononic graphene
,”
Nat. Mater.
15
,
1243
1247
(
2016
).
7.
Z.
Yang
,
F.
Gao
,
X.
Shi
,
X.
Lin
,
Z.
Gao
,
Y.
Chong
, and
B.
Zhang
, “
Topological acoustics
,”
Phys. Rev. Lett.
114
,
114301
(
2015
).
8.
J.
Lu
,
C.
Qiu
,
L.
Ye
,
X.
Fan
,
M.
Ke
,
F.
Zhang
, and
Z.
Liu
, “
Observation of topological valley transport of sound in sonic crystals
,”
Nature Phys.
13
,
369
374
(
2017
).
9.
D.
Torrent
,
D.
Mayou
, and
J.
Sánchez-Dehesa
, “
Elastic analog of graphene: Dirac cones and edge states for flexural waves in thin plates
,”
Phys. Rev. B
87
,
115143
(
2013
).
10.
J.
Vila
,
R. K.
Pal
, and
M.
Ruzzene
, “
Observation of topological valley modes in an elastic hexagonal lattice
,”
Phys. Rev. B
96
,
134307
(
2017
).
11.
X.
Huang
,
Y.
Lai
,
Z. H.
Hang
,
H.
Zheng
, and
C. T.
Chan
, “
Dirac cones induced by accidental degeneracy in photonic crystals and zero-refractive-index materials
,”
Nat. Mater.
10
,
582
586
(
2011
).
12.
P.
Moitra
,
Y.
Yang
,
Z.
Anderson
,
I. I.
Kravchenko
,
D. P.
Briggs
, and
J.
Valentine
, “
Realization of an all-dielectric zero-index optical metamaterial
,”
Nat. Photonics
7
,
791
795
(
2013
).
13.
M.
Dubois
,
C.
Shi
,
X.
Zhu
,
Y.
Wang
, and
X.
Zhang
, “
Observation of acoustic dirac-like cone and double zero refractive index
,”
Nat. Commun.
8
,
14871
(
2017
).
14.
Y.
Li
,
S.
Kita
,
P.
Muñoz
,
O.
Reshef
,
D. I.
Vulis
,
M.
Yin
,
M.
Lončar
, and
E.
Mazur
, “
On-chip zero-index metamaterials
,”
Nat. Photonics
9
,
738
742
(
2015
).
15.
C.
Xu
and
Y.
Lai
, “
Configurable dirac-like conical dispersions in complex photonic crystals
,”
Phys. Rev. B
95
,
045124
(
2017
).
16.
C.
Xu
,
H.
Chu
,
J.
Luo
,
Z. H.
Hang
,
Y.
Wu
, and
Y.
Lai
, “
Three-dimensional electromagnetic void space
,”
Phys. Rev. Lett.
127
,
123902
(
2021
).
17.
M.
Minkov
,
I. A. D.
Williamson
,
M.
Xiao
, and
S.
Fan
, “
Zero-index bound states in the continuum
,”
Phys. Rev. Lett.
121
,
263901
(
2018
).
18.
C.
Xu
,
G.
Ma
,
Z. G.
Chen
,
J.
Luo
,
J.
Shi
,
Y.
Lai
, and
Y.
Wu
, “
Three-dimensional acoustic double-zero-index medium with a fourfold degenerate dirac-like point
,”
Phys. Rev. Lett.
124
,
074501
(
2020
).
19.
L.-H.
Wu
and
X.
Hu
, “
Scheme for achieving a topological photonic crystal by using dielectric material
,”
Phys. Rev. Lett.
114
,
223901
(
2015
).
20.
C.
Xu
,
Z.-G.
Chen
,
G.
Zhang
,
G.
Ma
, and
Y.
Wu
, “
Multi-dimensional wave steering with higher-order topological phononic crystal
,”
Sci. Bull.
66
,
1740
1745
(
2021
).
21.
S.
Yves
,
R.
Fleury
,
T.
Berthelot
,
M.
Fink
,
F.
Lemoult
, and
G.
Lerosey
, “
Crystalline metamaterials for topological properties at subwavelength scales
,”
Nat. Commun.
8
,
16023
(
2017
).
22.
K. G.
Makris
,
R.
El-Ganainy
,
D. N.
Christodoulides
, and
Z. H.
Musslimani
, “
Beam dynamics in PT symmetric optical lattices
,”
Phys. Rev. Lett.
100
,
103904
(
2008
).
23.
A.
Szameit
,
M. C.
Rechtsman
,
O.
Bahat-Treidel
, and
M.
Segev
, “
PT-symmetry in honeycomb photonic lattices
,”
Phys. Rev. A
84
,
021806(R)
(
2011
).
24.
B.
Zhen
,
C. W.
Hsu
,
Y.
Igarashi
,
L.
Lu
,
I.
Kaminer
,
A.
Pick
,
S.-L.
Chua
,
J. D.
Joannopoulos
, and
M.
Soljačić
, “
Spawning rings of exceptional points out of dirac cones
,”
Front. Opt.
525
,
354
358
(
2015
).
25.
A.
Cerjan
,
A.
Raman
, and
S.
Fan
, “
Exceptional contours and band structure design in parity-time symmetric photonic crystals
,”
Front. Opt.
116
,
203902
(
2016
).
26.
H.
Zhou
,
C.
Peng
,
Y.
Yoon
,
C. W.
Hsu
,
K. A.
Nelson
,
L.
Fu
,
J. D.
Joannopoulos
,
M.
Soljačić
, and
B.
Zhen
, “
Observation of bulk Fermi arc and polarization half charge from paired exceptional points
,”
Science
359
,
1009
1012
(
2018
).
27.
K.
Ding
,
C.
Fang
, and
G.
Ma
, “
Non-Hermitian topology and exceptional-point geometries
,”
Nat. Rev. Phys.
4
,
745
760
(
2022
).
28.
K.
Ding
,
G.
Ma
,
M.
Xiao
,
Z. Q.
Zhang
, and
C. T.
Chan
, “
Emergence, coalescence, and topological properties of multiple exceptional points and their experimental realization
,”
Phys. Rev. X
6
,
021007
(
2016
).
29.
K.
Ding
,
G.
Ma
,
Z. Q.
Zhang
, and
C. T.
Chan
, “
Experimental demonstration of an anisotropic exceptional point
,”
Phys. Rev. Lett.
121
,
085702
(
2018
).
30.
L.
Ge
,
Y. D.
Chong
, and
A. D.
Stone
, “
Conservation relations and anisotropic transmission resonances in one-dimensional PT-symmetric photonic heterostructures
,”
Phys. Rev. A
85
,
023802
(
2012
).
31.
C.
Shi
,
M.
Dubois
,
Y.
Chen
,
L.
Cheng
,
H.
Ramezani
,
Y.
Wang
, and
X.
Zhang
, “
Accessing the exceptional points of parity-time symmetric acoustics
,”
Nat. Commun.
7
,
11110
(
2016
).
32.
J.
Zhang
,
B.
Peng
,
Ş. K.
Özdemir
,
K.
Pichler
,
D. O.
Krimer
,
G.
Zhao
,
F.
Nori
,
Y.
Liu
,
S.
Rotter
, and
L.
Yang
, “
A phonon laser operating at an exceptional point
,”
Nat. Photonics
12
,
479
484
(
2018
).
33.
L.
Feng
,
Z. J.
Wong
,
R.-M.
Ma
,
Y.
Wang
, and
X.
Zhang
, “
Single-mode laser by parity-time symmetry breaking
,”
Science
346
,
972
975
(
2014
).
34.
H.-X.
Wang
,
Y.
Chen
,
Z. H.
Hang
,
H.-Y.
Kee
, and
J.-H.
Jiang
, “
Type-II dirac photons
,”
npj Quantum Mater.
2
,
54
(
2017
).
35.
C.
Hu
,
Z.
Li
,
R.
Tong
,
X.
Wu
,
Z.
Xia
,
L.
Wang
,
S.
Li
,
Y.
Huang
,
S.
Wang
,
B.
Hou
,
C. T.
Chan
, and
W.
Wen
, “
Type-II dirac photons at metasurfaces
,”
Phys. Rev. Lett.
121
,
024301
(
2018
).
36.
X.
Wu
,
X.
Li
,
R.-Y.
Zhang
,
X.
Xiang
,
J.
Tian
,
Y.
Huang
,
S.
Wang
,
B.
Hou
,
C. T.
Chan
, and
W.
Wen
, “
Deterministic scheme for two-dimensional type-II dirac points and experimental realization in acoustics
,”
Phys. Rev. Lett.
124
,
075501
(
2020
).
37.
G. G.
Pyrialakos
,
N. S.
Nye
,
N. V.
Kantartzis
, and
D. N.
Christodoulides
, “
Emergence of type-II dirac points in graphyne like photonic lattices
,”
Phys. Rev. Lett.
119
,
113901
(
2017
).
38.
M. B.
de Paz
,
A.
González-Tudela
, and
P. A.
Huidobro
, “
Manipulating generalized dirac cones in subwavelength dipolar arrays
,”
Phys. Rev. A
106
,
033505
(
2022
).
39.
Z.
Yang
and
B.
Zhang
, “
Acoustic type-II Weyl nodes from stacking dimerized chains
,”
Phys. Rev. Lett.
117
,
224301
(
2016
).
40.
A. A.
Soluyanov
,
D.
Gresch
,
Z.
Wang
,
Q.
Wu
,
M.
Troyer
,
X.
Dai
, and
B. A.
Bernevig
, “
Type-II Weyl semimetals
,”
Nature
527
,
495
498
(
2015
).
41.
B.
Xie
,
H.
Liu
,
H.
Cheng
,
Z.
Liu
,
S.
Chen
, and
J.
Tian
, “
Experimental realization of type-II Weyl points and Fermi arcs in phononic crystal
,”
Phys. Rev. Lett.
122
,
104302
(
2019
).
42.
B.
Yang
,
Q.
Guo
,
B.
Tremain
,
L. E.
Barr
,
W.
Gao
,
H.
Liu
,
B.
Béri
,
Y.
Xiang
,
D.
Fan
,
A. P.
Hibbins
, and
S.
Zhang
, “
Direct observation of topological surface-state arcs in photonic metamaterials
,”
Nat. Commun.
8
,
97
(
2017
).
43.
J.
Noh
,
S.
Huang
,
D.
Leykam
,
Y. D.
Chong
,
K. P.
Chen
, and
M. C.
Rechtsman
, “
Experimental observation of optical Weyl points and Fermi arc-like surface states
,”
Nat. Phys.
13
,
611
617
(
2017
).
44.
W.-J.
Chen
,
M.
Xiao
, and
C. T.
Chan
, “
Photonic crystals possessing multiple Weyl points and the experimental observation of robust surface states
,”
Nat. Commun.
7
,
13038
(
2016
).
45.
X.
Huang
,
W.
Deng
,
F.
Li
,
J.
Lu
, and
Z.
Liu
, “
Ideal type-II weyl phase and topological transition in phononic crystals
,”
Phys. Rev. Lett.
124
,
206802
(
2020
).
46.
X.
Wu
,
H.
Fan
,
T.
Liu
,
Z.
Gu
,
R.-Y.
Zhang
,
J.
Zhu
, and
X.
Zhang
, “
Topological phononics arising from fluid-solid interactions
,”
Nat. Commun.
13
,
6120
(
2022
).
47.
C.
Liu
,
J.
Luo
, and
Y.
Lai
, “
Acoustic metamaterials with broadband and wide-angle impedance matching
,”
Phys. Rev. Materials
2
,
045201
(
2018
).
48.
C.
Liu
,
C.
Ma
,
X.
Li
,
J.
Luo
,
N. X.
Fang
, and
Y.
Lai
, “
Wide-angle broadband nonreflecting acoustic metamaterial fence
,”
Phys. Rev. Appl.
13
,
054012
(
2020
).
49.
S.
Longhi
, “
PT-symmetric laser absorber
,”
Phys. Rev. A
82
,
031801
(
2010
).
50.
Y. D.
Chong
,
L.
Ge
, and
A. D.
Stone
, “
Coherent perfect absorbers: Time-reversed lasers
,”
Phys. Rev. Lett.
106
,
093902
(
2011
).
51.
Z. J.
Wong
,
Y.-L.
Xu
,
J.
Kim
,
K.
O’Brien
,
Y.
Wang
,
L.
Feng
, and
X.
Zhang
, “
Lasing and anti-lasing in a single cavity
,”
Nat. Photonics
10
,
796
801
(
2016
).
52.
A.
Novitsky
,
D.
Lyakhov
,
D.
Michels
,
A. A.
Pavlov
,
A. S.
Shalin
, and
D. V.
Novitsky
, “
Unambiguous scattering matrix for non-Hermitian systems
,”
Phys. Rev. A
101
,
043834
(
2020
).

Supplementary Material