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 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 . In other words, a type-II Dirac point is induced below the sound cone, represented by the blue lines in the band structure.
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 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.
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 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.
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 , 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.
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.
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.
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 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.
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 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.
(a) and (b) Eigenvalues and 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.
(a) and (b) Eigenvalues and 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.
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 and of S2 are plotted in Figs. 5(a) and 5(b), exhibiting a unimodular behavior in the PT-exact phase and inverse in the PT-broken phase.30,52 At the CPA-lasing point, goes to the maximum, while 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 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.
CPA-lasing effect. (a) and (b) The eigenvalues and 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.
CPA-lasing effect. (a) and (b) The eigenvalues and 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.
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 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).
AUTHOR DECLARATIONS
Conflict of Interest
The authors have no conflicts to disclose.
Author Contributions
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).
DATA AVAILABILITY
The data that support the findings of this study are available from the corresponding author upon reasonable request.