One-Dimensional Mass-Spring Chains Supporting Elastic Waves with Non-Conventional Topology

Abstract: There are two classes of phononic structures that can support elastic waves with non-conventional topology, namely intrinsic and extrinsic systems. The non-conventional topology of elastic wave results from breaking time reversal symmetry (T-symmetry) of wave propagation. In extrinsic systems, energy is injected into the phononic structure to break T-symmetry. In intrinsic systems symmetry is broken through the medium microstructure that may lead to internal resonances. Mass-spring composite structures are introduced as metaphors for more complex phononic crystals with non-conventional topology. The elastic wave equation of motion of an intrinsic phononic structure composed of two coupled one-dimensional (1D) harmonic chains can be factored into a Dirac-like equation, leading to antisymmetric modes that have spinor character and therefore non-conventional topology in wave number space. The topology of the elastic waves can be further modified by subjecting phononic structures to externally-induced spatio-temporal modulation of their elastic properties. Such modulations can be actuated through photo-elastic effects, magneto-elastic effects, piezo-electric effects or external mechanical effects. We also uncover an analogy between a combined intrinsic-extrinsic systems composed of a simple one-dimensional harmonic chain coupled to a rigid substrate subjected to a spatio-temporal modulation of the side spring stiffness and the Dirac equation in the presence of an electromagnetic field. The modulation is shown to be able to tune the spinor part of the elastic wave function and therefore its topology. This analogy between classical mechanics and quantum phenomena offers new modalities for developing more complex functions of phononic crystals and acoustic metamaterials.


Introduction
A new frontier in wave propagation involves media that have broken time-reversal symmetry associated with non-conventional topology.Topological electronic [1], electromagnetic [2,3], and phononic crystals [4][5][6][7][8][9][10][11][12][13][14][15][16][17][18][19][20] all have demonstrated unusual topologically constrained properties.In phononic crystals and acoustic metamaterials, symmetry breaking is linked to constraints on the topological form of acoustic wave functions.For instance, in the context of topology, for the well-known driven damped oscillator, the amplitude of the wave function has properties isomorphic to the evolution of a field of parallel vectors tangent to a strip-like manifold and perpendicular to the length of the strip.The direction along the length of the strip represents frequency space and the strip has to exhibit a torsion (vectors in the vector field change orientation) at the oscillator resonant frequency as the amplitude changes sign as one crosses the resonance (i.e., the amplitude accumulates a π-phase shift).Dissipation aside, one of the most central elements to symmetry breaking and topology of elastic waves, is dispersion.A simple, linear one-dimensional (1D) harmonic monoatomic chain is a dispersive system, but one that obeys time-reversal symmetry and supports elastic waves with conventional Crystals 2016, 6, 44 2 of 16 topology in wave vector space.Perturbing the 1D harmonic chain through linear or nonlinear coupling may create resonant phonon modes that are dispersive, but whose amplitude may depend on the frequency and wave vector.In this case, the interplay between the coupling and dispersion of the system may lead to symmetry breaking conditions and therefore non-conventional elastic wave topology.There exists two-classes of phonon structures possessing non-conventional topology, namely intrinsic and extrinsic systems.Time-reversal symmetry in intrinsic systems [4][5][6][7][8][9][10][11][12][13][14] is broken through internal resonance or symmetry breaking structural features (e.g., chirality) and without addition of energy from the outside.Energy is added to extrinsic topological systems to break time reversal symmetry [15][16][17][18][19][20].A common example of an extrinsic approach is that of time-reversal symmetry breaking of acoustic waves by moving fluids [21][22][23][24][25][26][27][28][29].Recently, extrinsic topological phononic crystals have demonstrated the astonishing property of non-reciprocity and backscattering-immune edge states and bulk states establishing classical equivalents of topological electronic insulators.The non-conventional topology of elastic waves in an intrinsic topological phononic structure has been associated with the notion of duality in the quantum statistics of phonons (i.e., boson vs. fermion) [4,5].In the current paper, we illustrate the topological properties of elastic waves in the two classes of topological phononic structures.We first uncover the spinorial characteristics of elastic waves in crystals composed of connected masses and springs.An externally applied spatio-temporal modulation of the spring striffness can also be employed to break the symmetry of the system further [30].The modulation is shown to be able to tune the spinor part of the elastic wave function and therefore its topology.Sections 2 and 3 of this paper, address intrinsic and extrinsic topological phononic structures, respectively.Conclusions are drawn in Section 4.

Intrinsic Topological Phononic Structures
Let us consider a system composed of two coupled one-dimensional harmonic crystals as illustrated in Figure 1.
Crystals 2016, 6, 44 2 of 15 (1D) harmonic monoatomic chain is a dispersive system, but one that obeys time-reversal symmetry and supports elastic waves with conventional topology in wave vector space.Perturbing the 1D harmonic chain through linear or nonlinear coupling may create resonant phonon modes that are dispersive, but whose amplitude may depend on the frequency and wave vector.In this case, the interplay between the coupling and dispersion of the system may lead to symmetry breaking conditions and therefore non-conventional elastic wave topology.There exists two-classes of phonon structures possessing non-conventional topology, namely intrinsic and extrinsic systems.Time-reversal symmetry in intrinsic systems [4][5][6][7][8][9][10][11][12][13][14] is broken through internal resonance or symmetry breaking structural features (e.g., chirality) and without addition of energy from the outside.Energy is added to extrinsic topological systems to break time reversal symmetry [15][16][17][18][19][20].A common example of an extrinsic approach is that of time-reversal symmetry breaking of acoustic waves by moving fluids [21][22][23][24][25][26][27][28][29].Recently, extrinsic topological phononic crystals have demonstrated the astonishing property of non-reciprocity and backscattering-immune edge states and bulk states establishing classical equivalents of topological electronic insulators.The non-conventional topology of elastic waves in an intrinsic topological phononic structure has been associated with the notion of duality in the quantum statistics of phonons (i.e., boson vs. fermion) [4,5].In the current paper, we illustrate the topological properties of elastic waves in the two classes of topological phononic structures.We first uncover the spinorial characteristics of elastic waves in crystals composed of connected masses and springs.An externally applied spatio-temporal modulation of the spring striffness can also be employed to break the symmetry of the system further [30].The modulation is shown to be able to tune the spinor part of the elastic wave function and therefore its topology.Sections 2 and 3 of this paper, address intrinsic and extrinsic topological phononic structures, respectively.Conclusions are drawn in Section 4.

Intrinsic Topological Phononic Structures
Let us consider a system composed of two coupled one-dimensional harmonic crystals as illustrated in Figure 1.In absence of external forces, the equations describing the motion of atoms at location n in the two coupled 1-D harmonic crystals are given by: In absence of external forces, the equations describing the motion of atoms at location n in the two coupled 1-D harmonic crystals are given by: m B 2 u n Bt 2 " K 0 pu n`1 ´2u n `un´1 q `KI pv n ´un q " 0 (1a) here u and v represent the displacement in the upper and lower chains, respectively.These displacements can be visualized as being oriented along the chains.The side springs are illustrated for the sake of simplicity as vertical spring but physically they would couple masses between chains along the direction of the displacements.
In the long wavelength limit the discrete Lagrangian is expressed as a continuous second derivative of position.Taking M " m for the sake of simplicity and mathematical tractability, the equations of motion (1a,b) can be rewritten as: where I is the 2 ˆ2 identity matrix, D " ¸and U " ˜u v ¸is the displacement vector.We also have defined α 2 " K I m and β 2 " K 0 m .Equation (2) takes a form similar to the Klein-Gordon equation.Using an approach paralleling that of Dirac, one can factor Equation (2) into the following form: In Equation (3), we have introduced the 4 ˆ4 matrices: A " . C ¯Ψ " 0 in the plane wave form: ψ j " a j e ikx e iωt with j = 1,2,3,4.This gives the eigen value problem: where δ " α ?

2
. We find two dispersion relations: ω " ˘βk and ω " ˘bpβkq 2 `2α 2 .The first set of dispersion relations corresponds to branches that start at the origin k = 0 and relates to symmetric eigen modes.The second set of branches represents antisymmetric modes with a cut off frequency at k = 0 of α ?
2. Assuming that a 1 " a 2 " a F and that a 3 " a 4 " a B , for the symmetric waves characterized by the first set of dispersion relations, then the Equation (4) reduce to pω `βkq a B " 0 and pω ´βkq a F " 0 which are satisfied by plane waves of arbitrary amplitudes, a F and a B , propagating in the forward (F) or backward (B) directions, respectively.This is the conventional character of Boson-like phonons.We now look for the eigen vectors that correspond to the second set of dispersion relations.Let us use the positive eigen value as an illustrative example: ω " `bpβkq 2 `2α 2 .One of the degenerate solutions of the system of four linear Equation ( 4) is: where a 0 is some arbitrary constant.Note that the negative signs reflect the antisymmetry of the displacement.Other solutions can be found by considering the complete set of plane wave solutions ψ j " a j e ˘ikx e ˘iωt with j = 1,2,3,4 as well as the negative frequency eigen value.The key result is that the second dispersion curve in the band structure is associated with a wave function whose amplitude shows spinorial character (Equation ( 5)).In this case, the displacement of the two coupled harmonic chains are constrained and the direction of propagation of waves in the two-chain system are not independent of each other.For instance, at k = 0, the antisymmetric mode is represented by a standing wave which enforces a strict relation between the amplitude of a forward propagating wave and a backward propagating wave.This characteristic was shown [4,5] to be representative of Fermion-like behavior of phonons.As k Ñ 8 , ω Ñ `βk , the first two terms in Equation ( 5) go to zero and only one direction of propagation (backward) is supported by the medium (third and fourth terms in Equation ( 5)).This example illustrates the difference in topology of elastic waves corresponding to the lower and upper bands in the band structure of the two-chain system.The constraint on the amplitude of waves in the upper band imparts a nonconventional spinorial topology to the eigen modes which does exists for modes in the lower band.The topology of the upper band can be best visualized by taking the limit M Ñ 8 .In that case, the system of Figure 1 becomes a single harmonic chain grounded to a substrate as shown in Figure 2.
where is some arbitrary constant.Note that the negative signs reflect the antisymmetry of the displacement.Other solutions can be found by considering the complete set of plane wave solutions = ± ± with j = 1,2,3,4 as well as the negative frequency eigen value.The key result is that the second dispersion curve in the band structure is associated with a wave function whose amplitude shows spinorial character (Equation ( 5)).In this case, the displacement of the two coupled harmonic chains are constrained and the direction of propagation of waves in the two-chain system are not independent of each other.For instance, at k = 0, the antisymmetric mode is represented by a standing wave which enforces a strict relation between the amplitude of a forward propagating wave and a backward propagating wave.This characteristic was shown [4,5] to be representative of Fermion-like behavior of phonons.As → ∞, → + , the first two terms in Equation ( 5) go to zero and only one direction of propagation (backward) is supported by the medium (third and fourth terms in Equation ( 5)).This example illustrates the difference in topology of elastic waves corresponding to the lower and upper bands in the band structure of the two-chain system.The constraint on the amplitude of waves in the upper band imparts a nonconventional spinorial topology to the eigen modes which does exists for modes in the lower band.The topology of the upper band can be best visualized by taking the limit → ∞.In that case, the system of Figure 1 becomes a single harmonic chain grounded to a substrate as shown in Figure 2.
We now write our solutions in the form: and where and ̅ are two by one spinors.Inserting the In that limit, the displacement v in Equation (1a,b) is negligible.Equation ( 2) becomes the Klein-Gordon equation: `α2 u " 0 with α 2 " K I {m and β 2 " K 0 {m.This equation describes only the displacement field, u.Equation (3) can be written as the set of Dirac-like equations: where σ x and σ y are the 2 ˆ2 Pauli matrices: We now write our solutions in the form: Ψ k " Ψ pk, ω k q " c 0 ξ k pk, ω k q e p˘qiω k t e p˘qikx and Ψ k " Ψ pk, ω k q " c 0 ξ k pk, ω k q e p˘qiω k t e p˘qikx where ξ k and ξ k are two by one spinors.Inserting the various forms for these solutions in Equation (6a,b) lead to the same eigen values that we obtained before for the upper band of the two-chain system, namely by ω " ˘aα 2 `β2 k 2 .Again, let us note that the band structure has two branches corresponding to positive frequencies and negative frequencies.Negative frequencies can be visualized as representing waves that propagate in a direction opposite to that of waves with positive frequency.The spinor part of the solutions for the different plane waves is summarized in the Table 1 below.Negative and positive k correspond to waves propagating in opposite direction.1 can be used to identify the symmetry properties of Ψ and Ψ in the allowed space: k, ω.We find the following transformation rules: which lead to the combined transformation: We have defined as transformations that change the sign of the frequency and wave number, respectively.As one crosses the gap at the origin k = 0, the multiplicative factor "i" indicates that the wave function accumulated a phase of π 2 .Also, the Pauli operator σ x enables the transition from the space of solutions Ψ to the space of Ψ.We also note the orthogonality condition Ψσ x Ψ " 0. The topology of the spinorial wave functions that reflects their symmetry properties is illustrated in Figure 3. as transformations that change the sign of the frequency and wave number, respectively.As one crosses the gap at the origin k = 0, the multiplicative factor "i" indicates that the wave function accumulated a phase of .Also, the Pauli operator enables the transition from the space of solutions Ψ to the space of Ψ .We also note the orthogonality condition Ψ Ψ = 0.The topology of the spinorial wave functions that reflects their symmetry properties is illustrated in Figure 3.We now shed more light on the properties of the spinorial solutions by treating in Equation (6a,b) as a perturbation, .When = 0, Equation (6a) reduces to the two independent equations: We now shed more light on the properties of the spinorial solutions by treating α in Equation (6a,b) as a perturbation, ε.When α " 0, Equation (6a) reduces to the two independent equations: whose solutions correspond to plane waves propagating in the forward direction, ϕ p0q 1 , with dispersion relation ω `" βk and the backward direction, ϕ p0q 2 with dispersion relation ω ´" ´βk.We rewrite Equation ( 9) in the form: Seeking the particular solutions to first-order, ϕ p1q 1 and ϕ p1q 2 , which follow in frequency the driving terms on the right-side side of Equation (10a,b), we find their respective amplitude to be: Since ω `" ω ´at k = 0 only, the amplitude of the first-order perturbed forward wave, a 1 changes sign as k varies from ´8 to `8.A similar but opposite change of sign occurs for the backward perturbed amplitude.These changes of sign are therefore associated with changes in phase of π and ´π for the forward and backward waves as one crosses the origin k = 0.These phase changes (sign changes) are characteristic of that occurring at a resonance.The gap that would occur at k = 0 in the band structure of the harmonic chain grounded to a substrate via side springs, should we push the perturbation theory to higher orders, may therefore be visualized as resulting from a resonance of the forward waves driven by the backward propagating waves and vice versa.However, since to first-order the amplitudes given by Equation (11a,b) diverge at the only point of intersection between the dispersion relations of the forward and backward waves, we have to use analytic continuation to expand them into the complex plane: where η Ñ 0 continues the eigen values ω ´and ω `into the complex plane.
By simple inspection, we can see that at the origin, k = 0, both amplitudes are pure imaginary quantities and therefore exhibit a phase of ˘π 2 .This is expected as Equation (12a,b) are representative of the amplitude of a driven damped harmonic oscillator which also shows a phase of π 2 with respect to the driving frequency at resonance.It is instructive to calculate the Berry connection [31] for these perturbed amplitudes.The Berry connection determines the phase change of a wave as some parameter takes the wave function along a continuous path on the manifold that supports it.Since the Berry phase applies to continuous paths, we cannot use that concept to determine the phase change across the gap of our system, i.e., between the positive and negative frequency branches of the band structure.Therefore, we resort to calculating the Berry connection for the first-order perturbed solution which still remains continuous but may capture the interaction between directions of propagation.Our intention is to characterize the topology of the spinorial part of the wave function Equation (12a,b), by calculating the change in phase of the waves as one crosses k = 0.It is important first to normalize the spinor: ξ " ˜a1 a 2 ¸.This normalized spinor takes the form: r ξ " ˜r a 1 r a 2 ¸" ?
The Berry connection is given by A pkq " ´i r ξ ˚Br ξ Bk " ´ir a 1 ˚Br a 1 Bk ´ir a 2 ˚Br a 2 B k .After several analytical and algebraic manipulations, we obtain: In Equation ( 13), the contribution to the Berry connection of r a 1 and r a 2 are identical.Using the identities: lim ηÑ0 η x 2 `η2 " πδ pxq and δ paxq " 1 a δ pxq, leads to The contribution of each direction of propagation to the spinorial part of the wave function accumulates a π 2 phase shift as one crosses the origin k = 0.

Extrinsic Phononic Structure
In the previous section, we revealed the spinorial character of the elastic wave function in the two-chain system and one chain coupled to the ground.The purpose of this section is to investigate the behavior of field Ψ when the parameter α (spring stiffness, K I ) is subjected to a spatio-temporal modulation, i.e., α " α 0 `α1 2sin pKx `Ωtq where α 0 and α 1 are constants (see Figure 4).Here, K " 2π L where L is the period of the modulation.Ω is the frequency modulation and its sign determines the direction of propagation of the modulation.
The question arises as to the effect of such a modulation on the state of the fermion-like phonons.The periodicity of the modulated one-dimensional medium suggests that we should be seeking solutions of Equation (1) in the form of Bloch waves: Ψ px, tq " ř k ř g ψ pk, g, tq e ipk`gqx where x P r0, Ls.The wave number k is limited to the first Brillouin zone: " ´π L , π L ‰ and g " 2π L l with l being an integer.Choosing Equation (6a), we obtain the modulated Dirac-like equation in the Fourier domain: where k ˚" k `g.In the previous section, we revealed the spinorial character of the elastic wave function in the two-chain system and one chain coupled to the ground.The purpose of this section is to investigate the behavior of field Ψ when the parameter (spring stiffness, KI) is subjected to a spatio-temporal modulation, i.e., = + 2 sin + Ω where and are constants (see Figure 4).Here, = where is the period of the modulation.Ω is the frequency modulation and its sign determines the direction of propagation of the modulation.The question arises as to the effect of such a modulation on the state of the fermion-like phonons.The periodicity of the modulated one-dimensional medium suggests that we should be seeking solutions of Equation ( 1) in the form of Bloch waves: Ψ( , ) = ∑ ∑ ( , , ) ( )   where ∈ 0, .The wave number is limited to the first Brillouin zone: , and = with being an integer.Choosing Equation (6a), we obtain the modulated Dirac-like equation in the Fourier domain: where * = + .Consistent with Quantum Field Theory (QFT) approaches, we solve Equation (15) using perturbation theory and in particular multiple time scale perturbation theory [32] up to second-order.Nowadays, the multiple-time-scales perturbation theory (MTSPT) for differential Consistent with Quantum Field Theory (QFT) approaches, we solve Equation ( 15) using perturbation theory and in particular multiple time scale perturbation theory [32] up to second-order.Nowadays, the multiple-time-scales perturbation theory (MTSPT) for differential equations is a very popular method to approximate solutions of weakly nonlinear differential equations.Several implementations of this method were proposed in various fields of mathematics, mechanics and physics [33][34][35][36][37][38].Moreover, Khoo et al. [39] have shown that the MTSPT is a reliable theoretical tool for studying the lattice dynamics of an anharmonic crystal.More recently, Swinteck et al. [40] applied successfully the MTSPT, as described in Reference [19], for solving propagation equations in a quadratically nonlinear monoatomic chain of infinite extent.Consequently, we used the MTSPT for solving Equation (15).
The parameter α 1 is treated as a perturbation ε.The wave function is written as a second-order power series in ε, namely: ψ pk ˚, τ 0 , τ 1 , τ 2 q " ψ p0q pk ˚, τ 0 , τ 1 , τ 2 q `εψ p1q pk ˚, τ 0 , τ 1 , τ 2 q `ε2 ψ p2q pk ˚, τ 0 , τ 1 , τ 2 q here ψ pjq with j " 0, 1, 2 are wave functions expressed to zeroth, first and second-order.We have also replaced the single time variable, t, by three variables representing different time scales: τ 0 " t, τ 1 " εt, and τ 2 " ε 2 t " ε 2 τ 0 .We can subsequently decompose Equation ( 15) into equations to zeroth, first and second-order in ε.The zeroth-order equation consists of the Dirac-like equation in absence of modulation: As seen in Section 2, its solutions take the form ψ p0q pk ˚, τ 0 , τ 1 , τ 2 q " a p0q pk ˚, τ 1 , τ 2 q e iω 0 τ 0 with a p0q pk ˚, τ 1 , τ 2 q " parts of the solution, respectively.We have the usual eigen values: ω 2 0 " α 2 `β2 pk ˚q2 .Inserting the zeroth-order solution into Equation ( 15) expressed to first-order leads to The solutions of the first-order Dirac equation are the sum of solutions of the homogeneous equation and particular solutions.The homogenous solution is isomorphic to the zeroth-order solution, it will be corrected in a way similar to the zeroth-order solution as one accounts for higher and higher terms in the perturbation series.Under these conditions and to ensure that there are no secular terms (i.e., terms that grow with time and that are incompatible with the assumption that ψ p1q must be a correction to ψ p0q in the particular solution of Equation ( 17), the pre-factors of terms like e iω 0 τ 0 are forced to be zero [39].Subsequently, the derivative of the amplitudes a p0q pk ˚, τ 1 , τ 2 q with respect to τ 1 must vanish and these amplitudes only depend on τ 2 .The right-hand side of Equation ( 17) reduces to the second term only.The particular solution of that simplified equation contains frequency shifted terms given by: In the preceding relations, we have defined: ω 0 `Ω " ω 0 pk ˚`Kq `Ω and ω 0 ´Ω " ω 0 pk ˚´Kq `Ω.We have also omitted the time dependencies for the sake of compactness.

+
We note the asymmetry of these quantities.The terms G, G' and F diverge within the Brillouin zone of the modulated systems when the condition pω 0 pk ˚q `Ωq 2 ´ω2 0 pk ˚´Kq Ñ 0 is satisfied but not when pω 0 pk ˚q ´Ωq 2 ´ω2 0 pk ˚`Kq Ñ 0 .This asymmetry reflects a breaking of symmetry in wave number space due to the directionality of the modulation.
Equation (19a,b) impose second-order corrections onto the zeroth-order solution.We multiply the relations (18a,b) by ε 2 e iω 0 τ 0 to obtain them in terms of ψ p0q and subsequently recombine them with the zeroth-order Equation ( 3).This procedure reconstructs the perturbative series of Equation ( 15) in terms of ψ p0q only: with ˜∅k To obtain Equation ( 20), we have also used: . Equation (20) shows that Equation ( 15) describing the dynamics of elastic waves in a harmonic chain grounded to a substrate via side springs, whose stiffness is modulated in space and time, is to second-order isomorphic to Dirac equation in Fourier domain for a charged quasiparticle including an electromagnetic field.The quantity ∅ k ˚plays the role of the electrostatic potential and A k ˚the role of a scalar form of the vector potential.The parentheses ´B Bt ´i∅ k ˚¯and pik ˚´i A k ˚q are the Fourier transforms of the usual minimal substitution rule.α 0 `mk ˚is the dressed mass of the quasiparticle.The mechanical system provides a mechanism for exchange of energy between the main chain modes and the side springs.The side springs lead to the formation of a fermion-like quasiparticle while their modulation provides a field through which quasiparticles interact.The strength and nature of the interaction is controllable through the independent modulation parameters, α 1 ,Ω and K.
The mechanical system allows for the exploration of a large parameter space of scalar QFT as the functions ∅ k ˚and A k ˚can be varied by manipulating the spatio-temporal modulation of the side spring stiffness.We imagine, therefore, that this classical phononic system can be employed to examine the behavior of scalar QFT from weak to strong coupling regimes, as well as, at all intermediate couplings.Further, the capacity to separate the ratio of the effective potentials opens avenues for the experimental realization of scalar fields whose behavior could previously only have been theorized.
We now seek solutions of Equation (20).These are solutions of Equation ( 18) with spinorial part a p0q pk ˚, τ 2 q satisfying the second-order conditions given by Equation (19a,b).These conditions can be reformulated as Br a Bτ 2 " iMr a where the vector r a " We find the eigen values λ " The respective eigen vectors are r e " a 0 a ω 0 pk ˚q `βk The coefficients are given by C " # ´?ω 0 pk ˚q´βk ?ω . We note that although the quantities G, G' and F may diverge the ratio G´G1 2F remains finite.The directed spatio-temporal modulation impacts both the orbital part and the spinor part of the zeroth-order modes.The orbital part of the wave function is frequency shifted to ω 0 `ε2 λ and ω 0 `ε2 λ 1 .The quantities ε 2 λ t and ε 2 λ 1 t represent phase shifts analogous to those associated with the Aharonov-Bohm effect [41] resulting from electrostatic and vector potentials ∅ k ˚and A k ˚.Near the resonant condition: pω 0 pk ˚q `Ωq 2 ´ω2 0 pk ˚´Kq Ñ 0 , it is the eigen value λ Ñ 1 pω 0 pk ˚q`Ωq´ω 0 pk ˚´Kq which diverges.This divergence is indicative of the formation of a gap in the dispersion relation [31].The eigen value λ1 Ñ 1 pω 0 pk ˚q`Ωq`ω 0 pk ˚´Kq " 1 2pω 0 pk ˚q`Ωq does not diverge.Since the frequency shift, ε 2 λ 1 , is expected to be small compared to ω 0 pk ˚q, the orbital term e iλτ 2 Ñ 1 .Considering that the lowest frequency ω 0 is α 0 , this condition would occur for all k*.Therefore, the term e ipω 0 `ε2 λ 1 qτ 0 " e iω 0 τ 0 will essentially contribute to the band structure in a perturbative way similar to that of the uncorrected zeroth-order solution or homogeneous parts of the first or second-order equations.The spinorial part of the zeroth-order solution is also modified through the coupling between the orbital and "spin" part of the wave function as seen in the expressions for r eand r e 1 .This coupling suggests an approach for the manipulation of the "spin" part of the elastic wave function by exciting the medium using a spatio-temporal modulation.Again, these alterations can be achieved by manipulating independently the magnitude of the modulation, α 1 as well as the spatio-temporal characteristics Ω and K.
The perturbative approach used here is showing the capacity of a spatio-temporal modulation to control the "spin-orbit" characteristics of elastic modes in a manner analogous to electromagnetic waves enabling the manipulation of the spin state of electrons [42].However, the pertubative method is not able to give a complete picture of the effect of the modulation on the entire band structure of the elastic modes.For this, the vibrational properties of the mechanical system are also investigated numerically beyond perturbation theory.We calculate the phonon band structure of the modulated elastic Klein-Gordon equation since its eigen values are identical to those of the modulated Dirac-like equation.We use a one-dimensional chain that contains N = 2400 masses, m = 4.361 ˆ10 ´9 kg, with Born-Von Karman boundary conditions.The masses are equally spaced by h = 0.1 mm.The parameters K 0 = 0.018363 kg¨m 2 ¨s´2 and K I = 2295 kg¨s ´2.The spatial modulation has a period L = 100 h and an angular frequency Ω = 1.934 ˆ10 5 rad/s.We have also chosen the magnitude of the modulation: α 1 " 1 10 α 0 .The dynamics of the modulated system is amenable to the method of molecular dynamics (MD).The integration time step is dt = 1.624 ˆ10 -9 s.The dynamical trajectories generated by the MD simulation are analyzed within the framework of the Spectral Energy Density (SED) method [43] for generating the band structure.To ensure adequate sampling of the system's phase-space the SED calculations are averaged over 4 individual MD simulations, each simulation lasting 2 20 time steps and starting from randomly generated initial conditions.We report in Figure 5, the calculated band structure of the modulated system.magnitude of the modulation: = .The dynamics of the modulated system is amenable to the method of molecular dynamics (MD).The integration time step is dt = 1.624 × 10 -9 s.The dynamical trajectories generated by the MD simulation are analyzed within the framework of the Spectral Energy Density (SED) method [43] for generating the band structure.To ensure adequate sampling of the system's phase-space the SED calculations are averaged over 4 individual MD simulations, each simulation lasting 2 20 time steps and starting from randomly generated initial conditions.We report in Figure 5, the calculated band structure of the modulated system.Figure 5 retains the essential features of the unperturbed band structure but for frequency shifted Bloch modes ( * ) ± Ω and two band gaps in the positive half of the Brillouin zone.The periodicity of the band structure in wave number space is retained but its symmetry is broken.The The brighter branches correspond to the usual zeroth-order type wave (e iω 0 pk`gqτ 0 ).The fainter branches parallel to the brighter ones are characteristic of first-order waves (e ipω 0 pk`gq˘Ω qτ 0 ).

Wave number (xL)
Figure 5 retains the essential features of the unperturbed band structure but for frequency shifted Bloch modes ω 0 pk ˚q ˘Ω and two band gaps in the positive half of the Brillouin zone.The periodicity of the band structure in wave number space is retained but its symmetry is broken.The frequency shifted modes are illustrative of the first-order particular solutions.Second-order frequency shifted modes ω 0 pk ˚q ˘2Ω do not show in the figure due to their very weak amplitude.The two band gaps occur at the wave vector k gap defined by the condition `ω0 `kgap `g˘`Ω˘2 ´ω2 0 `kgap `g ´K˘" 0 for g = 0 and g = K.It is the band folding due to the spatial modulation which enables overlap and hybridization between the frequency-shifted Bloch modes and the original Bloch modes of the lattice without the time dependency of the spatial modulation.The hybridization opens gaps in a band structure that has lost its mirror symmetry about the origin of the Brillouin zone.Considering the first gap and following a path in k space, starting at k = 0 at the bottom of the lowest branch, the wave function transitions from a state corresponding to a zeroth-order type wave, with orbital part (e iω 0 pkqτ 0 ) and spinor part ˜aω 0 `βk a ω 0 ´βk ¸to a wave having the characteristics of the first-order wave with orbital part (e ipω 0 pk ˚q´Ω qτ 0 ) and spinor part Ω and K enables strategies for tuning the spinorial character of the elastic wave.The effect of these "spin-orbit" manipulations of the elastic system can be measured by examining the transmission of plane waves as shown in [4].
It is also instructive to consider the symmetry of the Dirac-like equations in the presence of a spatio-temporal modulation to best understand its effect on the spinorial character of the wave function.In the case of a modulation with a general phase ϕ, Equation (6a,b) take the overall form: Applying the joint T-symmetry and parity symmetry to Equation (21a) does not result in Equation (21b) for all phases ϕ but a few special values.The modulated Equations (21a,b) have lost the symmetry properties of the unmodulated Dirac equations Equation (6a,b).The gap that formed at k gap in Figure 4 is, therefore, not a Dirac point.The transformations T ω Ñ ω k Ñ ´k and T ω Ñ ´ω k Ñ k do not apply near k gap .The constraints imposed on the spinorial component of the elastic wave function may be released in the vicinity of that wavenumber.This constraint was associated with Fermion-like wave functions which have the character of quasistanding waves, i.e., composed of forward and backward waves with a very specific proportion of their respective amplitudes.The release of the Dirac constraint associated with the impossibility for the medium to support forward propagating waves (+k gap ) but only backward propagating waves (´k gap ), may lead again to Boson-like behavior with no restriction on the amplitude of the backward propagating waves.

Conclusions
We presented phononic structures composed of two coupled one-dimensional harmonic chains and one harmonic chained grounded to a substrate that exhibit intrinsic non-conventional topology.This topology is associated with wave functions that possess spinorial and orbital components.The spinorial character of the wave function imparts a fermion-like character to the phonons.This behavior is reflected in a constraint on the amplitude of forward and backward going waves.We also developed a scalar Quantum Field Theory that demonstrates the analogy between the one-dimensional elastic system subjected to a spatio-temporal modulation of its elastic properties and the one-dimensional Dirac equation including an electromagnetic field.The directional spatio-temporal modulation enables the tuning of the spinorial and orbital components of the wave function.Since the spatio-temporal characteristics of the modulation are independent of each other they offer exquisite means of controlling the spinor components of the elastic wave.Practical physical realization of the modulation of elastic medium stiffness could be achieved by exploiting a variety of non-contact approaches including the photo-elastic effect [44], the magneto-elastic effect [45], and contact approaches such as the piezoelectric effect [46] or externally induced mechanical deformations [47].The analogy between classical mechanical systems such as the ones demonstrated here and quantum and electromagnetic phenomena offer a new modality for developing more complex functions of phononic crystals and acoustic metamaterials.Such mechanical analogues of electromagnetic and quantum phenomena have a long history.For instance, Maxwell in his seminal paper "A dynamical theory of the electromagnetic field" [48] sought an elastic model of electrical and magnetic phenomena and electromagnetic waves.Other mechanical models of physical phenomena abound, including quantum mechanical behavior.For instance, the localization of ultrasound waves in two-dimensional [49] and three-dimensional [50] disordered phononic media serve as mechanical of Anderson localization of electrons.Tunneling of classical waves through phononic crystal barriers establishes a correspondence with its quantum counterpart [51,52].In addition to the mass spring systems reported in this paper, other physical systems that support elastic waves that can be described by Klein-Gordon-like equations include plates and phononic crystal plates [53], phononic crystals which support rotational elastic waves [5], granular phononic materials [54], elastic and sound wave guides with slowly varying cross sectional area [55].
Finally, we have addressed the topology of elastic wave functions and symmetry breaking in two linear 1-D mass-spring systems as well as in an extrinsic non-linear 1-D mass-spring system subjected to a spatio-temporal modulation of stiffness.More complex symmetry breaking conditions are expected to arise when considering intrinsic non-linear systems [56,57], such as for instance, phenomena associated with nonlinear resonances between linear and nonlinear systems.Indeed, coupled nonlinear and linear mechanical systems have recently received attention for their application in targeted energy transfer, whereby undesirable mechanical energy is directed irreversibly from a linear system to a nonlinear system [58].

Figure 1 .
Figure 1.Schematic illustration of the phononic structure composed of two coupled 1-D harmonic crystals.The atoms in the lower and upper 1-D harmonic crystals have mass m and M, respectively.The force constant of the springs of each 1-D harmonic crystal is taken to be the same, K0 .The force constant of the coupling springs is KI.The periodicity of the crystal is a.

Figure 1 .
Figure 1.Schematic illustration of the phononic structure composed of two coupled 1-D harmonic crystals.The atoms in the lower and upper 1-D harmonic crystals have mass m and M, respectively.The force constant of the springs of each 1-D harmonic crystal is taken to be the same, K 0 .The force constant of the coupling springs is K I .The periodicity of the crystal is a.

Figure 2 .
Figure 2. Schematic illustration of the harmonic chain grounded to a substrate via side springs (system in Figure 1 when → ∞).

Figure 2 .
Figure 2. Schematic illustration of the harmonic chain grounded to a substrate via side springs (system in Figure 1 when M Ñ 8 ).

Figure 3 .
Figure 3. Schematic representation of the manifold supporting Ψ and Ψ .The manifold exhibits a local quarter-turn twist around k = 0 of the square cross section of the manifold reflects the orthogonality of Ψ and Ψ .The clored arrows are parallel transported on the manifold along the direction of wave number.Their change in orientation is indicative of the phase change.

Figure 3 .
Figure 3. Schematic representation of the manifold supporting Ψ and Ψ.The manifold exhibits a local quarter-turn twist around k = 0 of the square cross section of the manifold reflects the orthogonality of Ψ and Ψ.The clored arrows are parallel transported on the manifold along the direction of wave number.Their change in orientation is indicative of the phase change.

Figure 4 .
Figure 4. Schematic illustration of the time evolution (top to bottom) of the harmonic chain grounded to a substrate with spatial modulation of the side spring stiffness illustrated as a "pink glow" of varying width when and a "blue glow" when .

Figure 4 .
Figure 4. Schematic illustration of the time evolution (top to bottom) of the harmonic chain grounded to a substrate with spatial modulation of the side spring stiffness illustrated as a "pink glow" of varying width when α ą α 0 and a "blue glow" when α ă α 0 .

Figure 5 .
Figure 5. Band structure of the mechanical model system of Figure 2 calculated using the Spectral Energy Density (SED) method.The band structure is reported as a contour plot of the natural logarithm of the SED (color bar) versus normalized frequency and reduced wave number.The frequency is normalized to the lowest value of unperturbed band, namely 1.507 × 10 5 Hz.The horizontal axis is extended to the right beyond the first Brillouin zone [−π,π] to highlight the asymmetry and therefore the modulation-induced symmetry breaking of the band structure.The brighter branches correspond to the usual zeroth-order type wave (

Figure 5 .
Figure 5. Band structure of the mechanical model system of Figure2calculated using the Spectral Energy Density (SED) method.The band structure is reported as a contour plot of the natural logarithm of the SED (color bar) versus normalized frequency and reduced wave number.The frequency is normalized to the lowest value of unperturbed band, namely 1.507 ˆ10 5 Hz.The horizontal axis is extended to the right beyond the first Brillouin zone [´π,π] to highlight the asymmetry and therefore the modulation-induced symmetry breaking of the band structure.The brighter branches correspond to the usual zeroth-order type wave (e iω 0 pk`gqτ 0 ).The fainter branches parallel to the brighter ones are characteristic of first-order waves (e ipω 0 pk`gq˘Ω qτ 0 ).

˜b1 1 b 1 2 ¸.
The control of the position of the gap through Crystals 2016, 6, 44 13 of 16