Abstract
The dynamical equations of motion for a discrete, one-dimensional harmonic chain with side restoring forces are analogous to the relativistic Klein–Gordon equation. Dirac factorization of the discrete Klein–Gordon equation introduces two equations with time reversal (T) and parity (P) symmetry-breaking conditions. The Dirac-factored equations enable the exploration of the properties of the solutions of the dynamic equations under PT symmetry-breaking conditions. The spinor solutions of the Dirac factored equations describe two types of acoustic waves: one with a conventional topology (Berry phase equal to 0) and the other with a non-conventional topology (Berry phase of π). In the latter case, the acoustic wave is isomorphic to the quantum spin of an electron, also known as an “acoustic pseudospin,” which requires a closed path, corresponding to two Brillouin zones (BZs), to restore the original spinor. We also investigate the topology of evanescent waves supported by the Dirac-factored equations. The interface between topologically conventional and non-conventional chains exhibits topological surface states. The Dirac-factored equations of motion of the one-dimensional harmonic chain with side springs can serve as a model for the investigation of the properties of acoustic topological insulators.
1 Introduction
A topological insulator cannot be adiabatically transformed into an ordinary insulator without passing through an intermediate conducting state (). While the bulk system is insulating, the surface can support conduction that is topologically protected, meaning that the surface states are insensitive to local perturbations (). Topological insulators can exhibit quantum mechanical properties such as the quantum spin Hall (QSH) effect () and the anomalous quantum Hall (QH) effect, which occurs in the absence of an external magnetic field due to the breaking of time-reversal symmetry (). Acoustic analogs of the QSH effect have been implemented (a) by tuning an accidental double Dirac cone in graphene-like lattices (; ) and (b) via a BZ folding mechanism (; ; ; ). QH-related effects have been realized in acoustics by arranging circulating flows into periodic settings to form an acoustic lattice that breaks time reversal symmetry (; ; ).
Phononic structures can support elastic waves with non-conventional topology by breaking symmetry (). Dirac factorization of the wave equation reveals potential topological properties that may result from symmetry breaking (which might be brought about by structural or external perturbations). For instance, the equations of motion of two coupled one-dimensional harmonic systems () can be factored in the long wavelength limit as a product of two Dirac-like equations, each of which breaks parity symmetry and time reversal symmetry. Propagative wave solutions to each Dirac equation can satisfy two possible dispersion relations, giving rise to symmetric and anti-symmetric eigenmodes. The former exhibit the conventional character of Boson-like phonons, while the latter exhibit Fermion-like behavior of phonons (; ). The topological properties of evanescent waves in product parity-time symmetry have also been investigated near exceptional points ().
The wave equation for a two-dimensional plate coupled to a rigid substrate can also be subject to Dirac factorization (). These factors are analogous to the long-wavelength limit of the Qi, Wu, and Zhang (QWZ) model of the anomalous quantum Hall effect (). The Dirac factorization reveals waves with spin-like degrees of freedom that have a gapped band structure, which is similar to the spin Hall effect. demonstrate a method inspired by the Dirac factorization of the Klein–Gordon equation to establish a connection between topological mechanical modes and the topological band theory in electronic systems. This leads to the prediction of new topological bulk mechanical phases with distinct boundary modes. Topological phonons can also be classified using local symmetries () by adapting the classification of non-interacting electron systems to mechanical systems.
Here, we study solutions of the Dirac factorizations of the discrete one-dimensional harmonic chain with side restoring forces and investigate the appearance of edge modes at the interface of conventional and non-conventional topologies. In Section 2, we introduce the Dirac factorization of the discrete Klein–Gordon equation, and in Section 3, we find the dispersion relation and amplitude vectors corresponding to propagative wave solutions of the Dirac equations. Section 4 addresses evanescent wave solutions and their dispersion relation. In Sections 5 and 6, we compute the respective Berry phases of the amplitude vectors of the propagative and evanescent waves. Section 7 addresses the existence of edge modes at the interface between two topologically different semi-infinite media that obey the acoustic Dirac equation. Finally, in Section 8, we summarize and draw conclusions.
2 Model system and equation of motion for the discrete harmonic chain
We consider the model of the one-dimensional harmonic chain illustrated in Figure 1. We assume that the chain lies along the -axis. The chain is composed of identical masses, , interacting with their neighbors via linear spring forces with a spring constant ; each of the masses is connected to a rigid substrate through harmonic springs with a spring constant . The coordinates of the mass at rest are , where is the spacing between adjacent masses at rest.
FIGURE 1
We denote the displacement of the mass by . Newton’s equation for the mass is
Taking and in Equation 1, we can rewrite the above equation as
The quantity can be identified as a discrete second derivative with respect to position. Thus, Equation 2 can be thought of as a discrete version of the Klein–Gordon equation It can be shown that the continuous Klein–Gordon equation (when interpreted as an operator acting on vectors with two components) can be factorized as a product of Dirac equations:—where , and are the Pauli matrices and is the 2 × 2 identity matrix ().
Now, let us provide a physical interpretation of amplitude vectors with two components. We note that if , the equation corresponding to the first operator in the square bracket of Equation 3 becomes . Using a plane wave solution with , this equation reduces to the system . We obtain two solutions for the angular velocity of the plane wave . These correspond to plane waves propagating in the positive and negative directions. In this case, the components and are now independent of each other and of the wave number. The amplitude of the plane wave propagating in the positive direction is independent of that of the wave propagating in the opposite direction. When , the plane wave solutions of the equation obey the system of linear equations with the dispersion relation . The amplitude components and are given by . The components of the × amplitude vector are not independent of each other. This indicates that the directions of propagation are not independent of each other anymore; it is parameter that couples those directions. The wave function has the character of quasi-standing waves, which are composed of forward and backward waves with a very specific proportion of their respective amplitudes and .
In order to achieve a factorization of the discrete Klein–Gordon equation as a product of Dirac-like equations, we will interpret the operator as acting on vectors with four components. Here, is the forward difference operator defined by , and is the backward difference operator defined by . These difference operators are linear (note that ). The operators and can also be interpreted as operators on vectors where the difference operations are carried out component-wise.
Let and We can reformulate Equation 2 in terms of a product of two operators:
It can be verified that this product is the same as (with being the 4 × 4 identity matrix in this equation). Note that the matrices in the above product are tensor products of 2 × 2 matrices and hence are 4 × 4 matrices. Thus, the operators in the factorization are acting on vectors with four components, and the product of the operators is a four-dimensional version of the discrete Klein–Gordon equation:—where .
Now the solution to Equation 4 is a vector. We have seen that the long wavelength limit (continuous limit) of the operators in the square brackets of Equation 3 leads to solutions with two components, corresponding to the mutually dependent amplitudes of the forward and backward waves. The four components of the discrete system reflect the fact that the forward difference and backward difference operators act differently on the forward and backward amplitudes of a quasi-standing wave. However, this solution form does not correspond to any straightforward physical interpretation.
The tensor products of the matrices appearing in the above equation are as follows:
Taking , , and , the Dirac factorization of the discrete Klein–Gordon equation becomes—where
3 Eigenvectors and dispersion relation
3.1 Dispersion relation
We will now find propagative solutions to the Dirac equations in Equation 5. Let us consider an ansatz taking the form of a plane wave with a amplitude vector, :
We have the matrices
. Then, with the ansatz in Equation 6, the Dirac equations become the following system of equations:
In matrix form, Equation 7 becomes
The determinant of this 4 × 4 matrix isand there exist non-zero solutions to Equation 8 if the determinant is zero. Equation 9 gives us
Equation 10 gives us the dispersion relations
The dispersion relation is illustrated in Figure 2.
FIGURE 2
3.2 Eigenvectors
We now solve Equation 8 for the components of the amplitude vector , , , and , satisfying the dispersion relations given by Equation 11. Let us redefine such that , , and . For the sake of simplifying the notation, we also define . With this notation, Equation 8 becomes the system of four linear equations:
To find solutions to the system of equations given by Equations 12 a–d, we hypothesize that
Recognizing that , , and , the normalized amplitude eigenvector is obtained in the form
4 Evanescent waves
In order to find non-propagative solutions of the Dirac equations (Equation 5), we consider the ansatz
We follow the procedure used in Section 3 to find the dispersion relation for waves of the form given by Equation 15:
When , then , and when = 0; that is, or equivalently . We denote this value of by . The points are illustrated schematically in Figure 3.
FIGURE 3
The normalized amplitude eigenvector is of the form—where with .
5 Berry phase for propagative waves
5.1 Continuous contribution to the Berry phase
The contribution of the continuous part of the function, to the Berry connection () is given by
Note that since , and therefore, Equation 18 gives us for all for which is continuous. So, the continuous contribution to the Berry connection is zero.
5.2 Contribution of discontinuities in the eigenvectors to the Berry connection
The eigenvector given by Equation 14 contains components in the form of square roots— and —which are square roots of complex numbers.
The square root of a complex number is given by the formula
Here, and . Near the origin on both sides (positive and negative) of the Brillouin zone: , so . We therefore have .
The phase of in the vicinity of the origin of the first Brillouin zone satisfies .
If —that is —then and The quantity remains continuous at the origin of the Brillouin zone; the same is true for the quantity Therefore, in this case, there is no discontinuity in the complex amplitude, which leads to the Berry phase being equal to zero, as discussed in in Section 5.1.
If , then —that is, and On the positive side of the origin of the Brillouin zone , and . On the negative side of the origin of the Brillouin zone , and . The quantity undergoes a π phase discontinuity at the origin of the Brillouin zone.
Considering the component , we still have but . We still have a discontinuity when , but then On the positive side of the origin of the Brillouin zone , and . On the negative side of the origin of the Brillouin zone , and . The quantity undergoes a -π phase discontinuity at the origin of the Brillouin zone.
In the case of to calculate the discontinuity contribution of the components at to the Berry phase, consider two amplitude vectors on both sides of the origin of the Brillouin zone:
Their inner product gives Therefore, the change in geometric phase as the amplitude vector is crossing the origin of the first Brillouin zone is . Now, the Berry phase is therefore the sum of the contributions from the continuous and the discontinuous parts of the amplitude vector, so the Berry phase is equal to . Note that this phase is independent of the specific value of as long as we consider the Dirac-factored equation with .
The Berry phase is a topological invariant of the system. Dirac-factored equations describe two types of acoustic waves: one with a conventional topology (Berry phase equal to 0) and the other one with an unconventional topology (Berry phase of π). In the latter case, the acoustic waves are isomorphic to the quantum spin of an electron, which requires a closed path, corresponding to two Brillouin zones to recover the original eigen amplitude vector. This is an example of an acoustic pseudospin.
6 Berry phase for evanescent waves
We now compute the Berry connection of the unit amplitude vector given by Equation 17 along the loop given by the dispersion relation of the evanescent waves with . Since is a real number, we have the restriction , which is equivalent to —that is, and
We first compute the Berry connection along the top and the bottom halves of the loop, when and
Case 1 with or .In this case, the unit amplitude vector is . Note that the quantities inside the square root are non-negative because of the restriction on , so the unit amplitude vector is a continuous function of Therefore, the Berry connection is given byWe note that since , which leads to Therefore, the Berry connection when is zero along the top and bottom halves of the loop.
Case 2 with or .In this case, the unit amplitude vector is . Again, the quantities inside the square root are non-negative. The Berry connection is thus given by A similar calculation to that in case gives us along the top and bottom halves of the loop.
Case 3The Berry connection near the points where the loop intersects the -axis.We wish to determine whether near , , the unit amplitude vector at with positive , is parallel or anti-parallel to , the unit amplitude vector at with negative . We calculate the dot productAs , note that . So in that limit, the dispersion relation given by Equation 16 gives us Therefore, when we have as When we have as So, when , as . Therefore, and are parallel. When , as So, and are anti-parallel.Similarly, when if , and are anti-parallel, and if , and are parallel.Therefore, along the closed loop over the evanescent mode, the Berry phase amounts to for both
7 Interface modes
We now consider a system composed of two semi-infinite chains described by the acoustic Dirac equation, but differing only in the value of the parameter . Such a system may be realized by considering the mass spring system in Figure 1 (which is governed by the Klein–Gordon equation) () and by considering evanescent waves corresponding to the Dirac equation for masses with a negative label and Dirac equation for masses with a positive label (since solutions to the Dirac equations are solutions to the Klein–Gordon equation). The consideration of such waves is mathematical. Figure 4 describes the set up.
FIGURE 4
At the interface, the Dirac equations are
The Equations 19 a,b can be expanded as
We now seek solutions to these equations that take the form of evanescent waves decaying on both sides of the interface. These are solutions of the form
Equations 21 a,b give are solutions to the bulk acoustic Dirac equation of infinite chains, each with its respective values of . At the interface, we consider
The existence of solutions in the forms given by Equations 22a,b, which satisfy Equations 20a,b, implies the existence of interface modes between the two chains with different topologies.
Inserting Equations 22a,b into Equations 20a,b gives us a system of eight linear equations in the eight unknowns , , , , , , , and :
Equation 23 has solutions if the matrix is invertible—that is, if its determinant is not equal to zero. By defining ,, and , we calculate the determinant of that matrix to be
Note that this determinant is independent of k’. To illustrate, let us set and plot .
The is non-zero for the majority of frequencies corresponding to evanescent waves, except for the two frequencies . There exist solutions to Equation 24 within this range of real frequencies. The interface between the chains with trivial and nontrivial topologies supports localized interfacial modes. These are topological interfacial modes. In Figure 5 there are two values of for which det() = 0. This implies that there are no real frequency solutions, although there may be solutions with complex frequencies. These waves may decay as a function of time.
FIGURE 5
8 Conclusion
We have here demonstrated that the Dirac factorization of the equations of motion of a mass and spring model exposes the potential for topological insulator behavior in acoustic systems. In particular, for a 1-d harmonic mass and spring chain attached elastically to a rigid substrate, the equations of motion give rise to a discrete version of the Klein–Gordon equation that can be factorized into Dirac equations with broken time-reversal and parity symmetry. Propagative and evanescent wave solutions of the Dirac factored equations were obtained. For propagative modes, the Berry phase for the Dirac equation with the + sign was found to be zero (conventional topology), and that with the – sign was found to be (non-conventional topology). In contrast, for the evanescent mode, the Berry phase was for both the + and – equations. Using the distinction between topologies for the propagative waves, we demonstrated the existence of a topologically protected interface mode between conventional and non-conventional topologies.
Dirac factorization of classical wave equations exposes the possibility of topological insulators arising from broken symmetries. It reveals the possibilities offered by symmetry breaking in terms of the direction of wave propagation. However, additional physical conditions or mechanisms are needed to break T- or P-symmetry and realize one-way propagating waves in physical systems.
Statements
Data availability statement
The original contributions presented in the study are included in the article/supplementary material; further inquiries can be directed to the corresponding author.
Author contributions
AB: Writing – original draft, Writing – review and editing. KR: Writing – original draft, Writing – review and editing. PD: Writing – original draft, Writing – review and editing.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This work was supported by the Science and Technology Center New Frontiers of Sound (NewFoS) through NSF cooperative agreement # 2242925.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1
BerryM. V. (1984). Quantal phase factors accompanying adiabatic changes. Proc. R. Soc. Lond. A39245–57. 10.1098/rspa.1984.0023
2
CalderinL.HasanM. A.JenkinsN. G.LataT.LucasP.RungeK.et al (2019). Experimental demonstration of coherent superpositions in an ultrasonic pseudospin. Sci. Rep.9, 14156. 10.1038/s41598-019-50366-y
3
ChangC.-Z.LiuC.-X.MacDonaldA. H. (2023). Colloquium: Quantum anomalous hall effect. Rev. Mod. Phys.95 (1), 011002. 10.1103/RevModPhys.95.011002
4
ChenZ.HeH.LiH.LiM.KouJ.LuY.et al (2024). Observation of parity-time symmetry for evanescent waves. Commun. Phys.7, 339. 10.1038/s42005-024-01816-1
5
DengY.GeH.TianY.LuM.JingY. (2017). Observation of zone folding induced acoustic topological insulators and the role of spin-mixing defects. Phys. Rev. B96, 184305. 10.1103/PhysRevB.96.184305
6
DeymierP.RungeK. (2016). One-dimensional mass-spring chains supporting elastic waves with non-conventional topology. Crystals6 (4), 44. 10.3390/cryst6040044
7
DeymierP. A.RungeK. (2022). Revealing topological attributes of stiff plates by dirac factorization of their 2D elastic wave equation. Appl. Phys. Lett.120 (8), 081701. 10.1063/5.0086559
8
DeymierP. A.RungeK.SwinteckN.MuralidharanK. (2014). Rotational modes in a phononic crystal with fermion-like behavior. J. Appl. Phys.115, 163510. 10.1063/1.4872142
9
DeymierP. A.RungeK.SwinteckN.MuralidharanK. (2015). Torsional topology and fermion-like behavior of elastic waves in phononic structures. Comptes Rendus Mécanique343, 700–711. 10.1016/j.crme.2015.07.003
10
HasanM. Z.KaneC. L. (2010). Colloquium: topological insulators. Rev. Mod. Phys.82 (4), 3045–3067. 10.1103/RevModPhys.82.3045
11
HeC.NiX.GeH.SunX.-C.ChenY.-B.LuM.-H.et al (2016). Acoustic topological insulator and robust one-way sound transport. Nat. Phys.12, 1124–1129. 10.1038/nphys3867
12
KaneC. L.LubenskyT. C. (2014). Topological boundary modes in isostatic lattices. Nat. Phys.10, 39–45. 10.1038/nphys2835
13
KhanikaevA. B.FleuryR.MousaviS. H.AlùA. (2015). Topologically robust sound propagation in an angular-momentum-biased graphene-like resonator lattice. Nat. Commun.6, 8260. 10.1038/ncomms9260
14
MeiJ.ChenZ.WuY. (2016). Pseudo-time-reversal symmetry and topological edge states in two-dimensional acoustic crystals. Sci. Rep.6, 32752. 10.1038/srep32752
15
NiX.HeC.SunX.-C.LiuX.-p.LuM.-H.FengL.et al (2015). Topologically protected one-way edge mode in networks of acoustic resonators with circulating air flow. New J. Phys.17, 053016. 10.1088/1367-2630/17/5/053016
16
QiX.-L.ZhangS.-C. (2010). The quantum spin hall effect and topological insulators. Phys. Today63 (1), 33–38. 10.1063/1.3293411
17
QiX.-L.ZhangS.-C. (2011). Topological insulators and superconductors. Rev. Mod. Phys.83 (4), 1057–1110. 10.1103/RevModPhys.83.1057
18
QiX.-L.WuY.-S.ZhangS.-C. (2006). Topological quantization of the spin hall effect in two-dimensional paramagnetic semiconductors. Phys. Rev. B74, 085308. 10.1103/PhysRevB.74.085308
19
SüsstrunkR.HuberS. D. (2016). Classification of topological phonons in linear mechanical metamaterials. Proc. Natl. Acad. Sci. U.S.A.113 (33), E4767–E4775. 10.1073/pnas.1605462113
20
XiaB.-Z.LiuT.-T.HuangG.-L.DaiH.-Q.JiaoJ.-R.ZangX.-G.et al (2017). Topological phononic insulator with robust pseudospin-dependent transport. Phys. Rev. B96, 094106. 10.1103/PhysRevB.96.094106
21
XueH.YangY.ZhangB. (2022). Topological acoustics. Nat. Rev. Mater7, 974–990. 10.1038/s41578-022-00465-6
22
YangZ.GaoF.ShiX.LinX.GaoZ.ChongY.et al (2015). Topological acoustics. Phys. Rev. Lett.114, 114301. 10.1103/PhysRevLett.114.114301
23
YvesS.FleuryR.LemoultF.FinkM.LeroseyG. (2017). Topological acoustic polaritons: robust sound manipulation at the subwavelength scale. New J. Phys.19, 075003. 10.1088/1367-2630/aa66f8
24
ZhangZ.WeiQ.ChengY.ZhangT.WuD.LiuX. (2017). Topological creation of Acoustic pseudospin multipoles in a flow-free symmetry-broken metamaterial lattice. Phys. Rev. Lett.118, 084303. 10.1103/PhysRevLett.118.084303
Summary
Keywords
topological insulator, Berry phase, Dirac equation, Klein-Gordon equation, interface mode
Citation
Basu A, Runge K and Deymier PA (2025) The acoustic Dirac equation as a model of topological insulators. Front. Acoust. 3:1615210. doi: 10.3389/facou.2025.1615210
Received
20 April 2025
Accepted
29 August 2025
Published
29 September 2025
Volume
3 - 2025
Edited by
Michael J. Leamy, Georgia Institute of Technology, United States
Updates
Copyright
© 2025 Basu, Runge and Deymier.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Abhirup Basu, abhirup.basu56@gmail.com
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.