Next Article in Journal
Asymptotic Solutions of Steady Lamellar Eutectic Growth in Directional Solidification for Small Tangent Values of the Contact Angles
Next Article in Special Issue
Immunity to Backscattering of Bulk Waves in Topological Acoustic Superlattices
Previous Article in Journal
Piezoelectric Properties and Thermal Stability of Pb(Yb1/2Nb1/2)O3-BiScO3-PbTiO3 Ternary Ceramics
Previous Article in Special Issue
Phononic Band Structure by Calculating Effective Parameters of One-Dimensional Metamaterials
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Pseudo-Spin Polarized One-Way Elastic Wave Eigenstates in One-Dimensional Phononic Superlattices

by
Pierre A. Deymier
1,*,†,
Keith Runge
1,†,
Alexander Khanikaev
2,3,† and
Andrea Alù
3,4,†
1
Department of Materials Science and Engineering, University of Arizona, Tucson, AZ 85721, USA
2
Department of Electrical Engineering, The City College of New York, New York, NY 10031, USA
3
Physics Program, Graduate Center, City University of New York, New York, NY 10016, USA
4
Photonics Initiative, Advanced Science Research Center, City University of New York, New York, NY 10031, USA
*
Author to whom correspondence should be addressed.
Current affiliation: New Frontiers of Sound Science and Technology Center, The University of Arizona, Tucson, AZ 85721, USA.
Crystals 2024, 14(1), 92; https://doi.org/10.3390/cryst14010092
Submission received: 15 December 2023 / Revised: 11 January 2024 / Accepted: 13 January 2024 / Published: 19 January 2024
(This article belongs to the Special Issue Advances and Application of Phononic Crystals)

Abstract

:
We investigate a one-dimensional discrete binary elastic superlattice bridging continuous models of superlattices that showcase a one-way propagation character, as well as the discrete elastic Su–Schrieffer–Heeger model, which does not exhibit this character. By considering Bloch wave solutions of the superlattice wave equation, we demonstrate conditions supporting elastic eigenmodes that do not satisfy the translational invariance of Bloch waves over the entire Brillouin zone, unless their amplitude vanishes for a certain wave number. These modes are characterized by a pseudo-spin and occur only on one side of the Brillouin zone for a given spin, leading to spin-selective one-way wave propagation. We demonstrate how these features result from the interplay of the translational invariance of Bloch waves, pseudo-spins, and a Fabry–Pérot resonance condition in the superlattice unit cell.

1. Introduction

Breaking reciprocity is one principle for forming acoustic or elastic waves traveling only in one direction (that is, one-way wave propagation). A variety of mechanisms can be used to break reciprocity, which range from symmetry breaking to exploiting non-linearity and spatiotemporal modulations [1,2]. For example, breaking time-reversal symmetry by using gyroscopic inertial effects in lattice structures can lead to the topologically protected one-way propagation of elastic waves localized at the edges of the lattice [3]. The dependence on the amplitude of the dynamic behavior of the non-linear waves may also lead to non-reciprocity [4]. Moreover, the modulation in space and time of the physical properties of elastic media introduces a bias that breaks time-reversal symmetry and subsequently reciprocity, leading to the one-way propagation of bulk elastic waves [5,6,7].
While breaking reciprocity is one route to one-way propagation [8,9], other topologically protected systems that preserve reciprocity have been demonstrated. For example, static one-dimensional superlattices that do not break time-reversal symmetry and obey linear elasticity have been shown to support robust, topologically protected elastic eigenwaves with non-zero amplitude in the forward propagation direction but with zero amplitude in the opposite direction [10]. These one-way-propagating, topologically protected eigenwaves occur along elastic bands with a non-zero Berry phase. Due to time-reversal symmetry, these eigenwaves possess a pseudo-spin that breaks mirror symmetry, and an eigenwave with opposite handedness is expected to be supported by the system in the reverse propagation direction. Hence, differing from non-reciprocal one-way systems, these responses do not generally isolate but can be used to efficiently direct signals towards desired directions and are robust against back-reflections, as long as disorder or defects do not couple the two pseudo-spin hands.
Here, we develop a model for one-dimensional discrete phononic superlattices, which enables us to bridge continuum superlattices that support reciprocal one-way eigenwaves and the elastic analogue of the well-known one-dimensional Su–Schrieffer–Heeger (SSH) model [11,12], which does not generally support topologically protected, one-way-propagating bulk waves. The discrete superlattice model is investigated using the transfer matrix method to calculate its band structure and corresponding eigenmodes. We show that the transfer matrix exhibits, under specific conditions, eigenvalues that do not span the complete unit circle in the complex plane. Subsequently, such eigenmodes do not satisfy the translational invariance of the elastic modes in the form of Bloch waves over the entire Brillouin zone, unless their amplitude vanishes for a certain wave number. This condition occurs on one side of the Brillouin zone and not the other, leading to one-way wave propagation. In the case of the continuous limit of a finite superlattice, we show that reciprocity is obeyed, as expected, and the one-way nature is associated with an inherent handedness of the eigenmodes (pseudo-spin), which breaks the mirror symmetry. One-way propagation is identified as the combined effect of translational invariance, pseudo-spin, and Fabry–Pérot resonances [13] in one of the constitutive media of the superlattice. These results offer an original view of unidirectional propagation, which can be probed via exciting one-way modes by spin-selective sources in 1D periodic systems.

2. Model and Method

2.1. Model System and Dynamic Equations of Motion

The binary superlattice system under study is formed by periodically repeating unit cells composed of two segments of different one-dimensional mass–spring harmonic chains (see Figure 1). Each mass is identical and equals 1; the stiffness levels of the springs in segments 1 and 2 equal K 1 and K 2 , respectively; and the spacing between adjacent masses equals a . The lengths of the segments are d 1 = n a and d 2 = p a , where n and p are integers. Segments 1 and 2 are separated by interfaces labelled I and II. In this model, we assume that the stiffness of the springs (bond strength) in a segment remains the same up to the interfaces. When n = 1 and p = 1 , the system becomes the Su–Schrieffer–Heeger (SSH) model [12]. For large values of n and p , in the long wavelength limit, the system approaches a continuous superlattice.
The displacement of a mass, m [ 0 , n 1 ] , in segment 1 of unit cell N is defined as u A , m N . Similarly, the displacement of masses in segment 2 of unit cell N is labelled u B , l N , with l [ 0 , p 1 ] . The masses obey bulk equations of motion in segments 1 and 2 for m [ 1 , n 2 ] and l [ 1 , p 2 ] , namely:
d 2 u A , m N d t 2 = K 1 u A , m + 1 N u A , m N K 1 u A , m N u A , m 1 N
d 2 u B , l N d t 2 = K 2 u B , l + 1 N u B , l N K 2 u B , l N u B , l 1 N
Equation (1a,b) supports plane wave solutions of the form:
u A , m N = A + N e i k 1 m a + A N e i k 1 m a e i ω t
u B , l N = B + N e i k 2 l a + B N e i k 2 l a e i ω t
The wave number k j with j = 1, 2 is related to the angular frequency via the well-known dispersion relation of infinite harmonic chains:
ω 2 = 4 K i sin k j a 2 2
The behavior of the superlattice will, therefore, be dependent on the overlap of the bulk bands (density of states) of the materials forming segments 1 and 2.
Since the system is periodic with period L = d 1 + d 2 , we consider Bloch wave solutions, i.e., we choose A ± N = e i q N L A ± and B ± N = e i q N L B ± , where q is the wave number. We remove the upper script N on A ± and B ± when the periodicity is explicitly accounted for in the plane wave term e i q N L .
The equations of the motion of masses at interfaces I and II are as follows:
d 2 u A , 0 N d t 2 = K 1 u A , 1 N u A , 0 N K 2 u A , 0 N u B , p 1 N 1
d 2 u B , 0 N d t 2 = K 2 u B , 1 N u B , 0 N K 1 u B , 0 N u A , n 1 N
Inserting the null terms K 1 u A , 0 N u A , p 1 N 1 K 1 u A , 0 N u A , p 1 N 1 into Equation (4a) and K 2 u B , 0 N u B , 1 N K 2 u B , 0 N u B , 1 N into Equation (4b), we obtain, by virtue of the bulk equations of motion, the conditions of the continuity of forces at interfaces I and II:
K 1 u A , 0 N u A , p 1 N 1 = K 2 u A , 0 N u B , p 1 N 1
K 2 u B , 0 N u B , 1 N = K 1 u B , 0 N u A , n 1 N
We also introduce the following conditions of the continuity of displacement at I and II:
u A , 0 N = u B , p N 1
u B , 0 N = u A , n N
With these conditions, Equation (5a,b) becomes:
K 1 u A , 0 N u A , 1 N = K 2 u B , p N 1 u B , p 1 N 1
K 2 u B , 0 N u B , 1 N = K 1 u A , n N u A , n 1 N
In Equation (7a), we use u A , p 1 N 1 = u A , 1 N . Inserting the general solutions given by Equation (2a,b) and the Bloch wave form into Equations (6a,b) and (7a,b) yields the following set of dynamic equations:
α 1 β 1 1 1 f α 1 δ ( 1 ) f β 1 δ + ( 1 ) δ ( 2 ) δ + ( 2 ) 1 1 e i q L α 2 e i q L β 2 f δ ( 1 ) f δ + ( 1 ) e i q L α 2 δ ( 2 ) e i q L β 2 δ + ( 2 ) A + A B + B = 0
where f = K 1 K 2 , α 1 = 1 β 1 = e i k 1 n a , α 2 = 1 β 2 = e i k 2 p a , and δ ± ( j ) = 1 e ± i k j a , with j = 1, 2. We also enforce the periodicity of the system, requiring that the solutions in Equation (2a,b) take the form of Bloch waves; that is, B ± N 1 = B ± N e i q L .
Equation (8) possesses non-trivial solutions when the determinant of the dynamic matrix is equal to zero. This condition leads to the dispersion relation:
f δ + ( 2 ) δ ( 2 ) δ + ( 1 ) δ ( 1 ) 2 cos q L α 1 β 2 + β 1 α 2 δ + ( 2 ) f δ + ( 1 ) δ ( 2 ) f δ ( 1 ) + α 1 α 2 + β 1 β 2 δ + ( 2 ) f δ ( 1 ) δ ( 2 ) f δ + ( 1 ) = 0 .
This dispersion relation can be reformulated as:
cos q L = cos k 1 d 1 cos k 2 d 2 + 1 2 1 f 4 sin k 2 a 2 2 sin k 1 a sin k 2 a + f 4 sin k 1 a 2 2 sin k 1 a sin k 2 a + 4 sin k 1 a 2 2 sin k 2 a 2 2 sin k 1 a sin k 2 a sin k 1 d 1 sin k 2 d 2

2.2. SSH Limit

In the SSH scenario ( n = p = 1 ), sin k 1 d 1 sin k 2 d 2 = sin k 1 a sin k 2 a and cos k 1 d 1 cos k 2 d 2 = cos k 1 a cos k 2 a . Using the dispersion relations given by Equation (3) and basic trigonometric relations, Equation (10) reduces to:
ω 4 2 K 1 + K 2 ω 2 + 2 K 1 K 2 1 cos q L = 0
This leads to the usual SSH dispersion relation:
ω 2 = K 1 + K 2 ± K 1 2 + K 2 2 + 2 K 1 K 2 cos q L
Equation (8) can be manipulated algebraically to obtain a relationship between A = A + + A and B = B + + B , namely:
A ω 2 + K 1 + K 2 = K 1 + K 2 e i q L B
The amplitudes A = A + + A and B = B + + B do not differentiate between positive or negative k j , since m = 0 and l = 0 in Equation (2a,b).
Employing Equation (12) in the form of ω 2 + K 1 + K 2 2 = K 1 + K 2 e i q L K 1 + K 2 e + i q L , we recover the complex SSH amplitudes [14]:
A B K 1 + K 2 e i q L K 1 + K 2 e + i q L
The dispersion relation (Equation (12)) exhibits band inversion when varying the spring stiffness. The complex amplitudes (Equation (14)) support a non-conventional topology over the Brillouin zone, q [ π L ,   π L ], when K 1 < K 2 . It is worth noting that these amplitudes cannot be zero for any value of the wave number q.

2.3. Continuum Limit

In the long wavelength limit, δ ± ( j ) i k j a . For large values of n and p, the dispersion relation Equations (9) or (10) reduces to the known expression [15]:
cos q L = cos k 1 d 1 cos k 2 d 2 1 2 1 F + F sin k 1 d 1 sin k 2 d 2
In Equation (15), we define F = k 1 K 1 k 2 K 2 . The solutions for the amplitudes are addressed separately in the next subsection.

2.4. Amplitudes of the Displacement Field

To find the amplitudes, we use the transfer matrix method. Equation (8) can be rewritten as:
1 1 f δ ( 1 ) f δ + ( 1 ) A + N A N = α 2 β 2 α 2 δ ( 2 ) β 2 δ + ( 2 ) B + N 1 B N 1
α 1 β 1 f α 1 δ ( 1 ) f β 1 δ + ( 1 ) A + N A N = 1 1 δ ( 2 ) δ + ( 2 ) B + N B N
Since the second relation, Equation (16b), applies to the unit cell N-1, as well as cell N, it can be inverted to obtain a relation between B + N 1 B N 1 and A + N 1 A N 1 . By inserting that relation into Equation (16a), we obtain the transfer matrix relating amplitudes in one unit cell to amplitudes in the neighboring cell:
A + N A N = t 11 t 12 t 21 t 22 A + N 1 A N 1
The components of the transfer matrix are given by:
t 11 = 1 f δ + ( 2 ) δ ( 2 ) δ + ( 1 ) δ ( 1 ) α 1 f δ + ( 1 ) δ ( 2 ) δ + ( 2 ) f δ ( 1 ) α 2 + f δ + ( 1 ) δ + ( 2 ) δ 2 + f δ ( 1 ) β 2
t 12 = 1 f δ + ( 2 ) δ ( 2 ) δ + ( 1 ) δ ( 1 ) β 1 f δ + ( 1 ) δ ( 2 ) δ + ( 2 ) f δ + ( 1 ) α 2 + f δ + ( 1 ) δ + ( 2 ) δ 2 + f δ + ( 1 ) β 2
Here, t 21 = t 12 * and t 22 = t 11 * . The star stands for a complex conjugation.
Due to the periodicity of the system, solutions will take the form of Bloch waves:
A + N A N = e i q L A + N 1 A N 1
Combining Equations (18) and (17) leads to the eigenvalue problem [16]:
t 11 t 12 t 21 t 22 e i q L 1 0 0 1 A + A = 0
If e i q L is an eigenvalue of the transfer matrix, then the amplitudes A + and A are given by:
A + = t 12
A = t 11 e i q L
The dispersion relation of Equation (9) can be rewritten as 2 cos q L = t 11 + t 22 , so Equation (20b) becomes:
A = 1 2 t 11 t 22 i sin q L
After a number of trigonometric and algebraic manipulations, we obtain:
A = i sin k 1 d 1 cos k 2 d 2 + 1 δ + ( 2 ) δ ( 2 ) δ + ( 1 ) δ ( 1 ) 2 f δ + ( 1 ) δ ( 1 ) + 1 f δ + 2 δ 2 + δ + ( 1 ) + δ ( 1 ) δ + ( 2 ) + δ ( 2 ) cos k 1 d 1 sin k 2 d 2 sin q L
Furthermore, the amplitude A + N is given by:
A + = 2 i f δ + ( 2 ) δ ( 2 ) δ + ( 1 ) δ ( 1 ) e i k 1 d 1 f δ + ( 1 ) δ + ( 2 ) f δ + ( 1 ) δ ( 2 ) sin k 2 d 2
In the continuum limit, these amplitudes simplify to the known expressions [10]:
A + = 1 2 F 1 F sin k 1 d 1 sin k 2 d 2 + i 2 F 1 F cos k 1 d 1 sin k 2 d 2
A = i sin k 1 d 1 cos k 2 d 2 + 1 2 F + 1 F cos k 1 d 1 sin k 2 d 2 sin q L

2.5. Conditions for One-Way Propagation in Infinite Superlattices

For the continuous system, we have shown [17] that if a band contains the frequency point ω 0 , such that sin k 2 d 2 = 0 , then at ω 0 ( q 0 ) the dispersion relation becomes cos q 0 L = cos k 1 d 1 + k 2 d 2 , A + = 0 , and A = i sin ( k 1 d 1 + k 2 d 2 ) sin q 0 L . This condition corresponds to a Fabry–Pérot resonance of the second medium in the superlattice unit cell [13].
Since k 1 d 1 + k 2 d 2 > 0 , the sign of q 0 determines whether A vanishes or not. From Equation (16b), A + = 0 and A = 0 , meaning that B + = 0 and B = 0 . At ω 0 ( + q 0 ) , the amplitudes of the Bloch modes may be finite, while at ω 0 ( q 0 ) the Bloch modes equal zero. This behavior was shown to be associated with the non-conventional topology of the Bloch modes of the continuous superlattice [18]. The Berry connection undergoes a π jump at q 0 , where the amplitudes vanish. The dispersion branches supporting zero-amplitude modes in the band structure of the superlattice exhibit a π Berry phase [18].
This argument also extends to a discrete superlattice with p > 1 in the following manner. Let us first find expressions for the eigenvalues of the transfer matrix. These eigenvalues are solutions of the equation λ 2 λ t 11 + t 11 * + t 11 t 11 * t 12 t 12 * = 0 . The Eigen values are given by:
λ = 1 2 t 11 + t 11 * ± 1 2 t 11 t 11 * 2 + 4 t 12 t 12 *
For vibrational modes to be Bloch waves, we require λ e ± i q L . According to Equation (23), A + = t 12   = 0 when sin k 2 d 2 = 0 . In that case, the eigenvalues reduce to:
λ = 1 2 t 11 + t 11 * ± 1 2 t 11 t 11 *
or:
λ = R e a l ( t 11 )   ± i I m a g ( t 11 )
where Real and Imag stand for the real and imaginary parts.
The condition for the discrete superlattice to support Bloch waves takes the form of R e a l t 11 ± I m a g t 11 = cos q L ± i sin q L . This condition is met only for a wave number q located on one side of the Brillouin zone, where I m a g t 11 can be identified as a sine function. For q on the other side of the Brillouin zone, I m a g t 11 is identifiable to minus one times a sine function, which has the opposite sign of sin q L . In other words, λ does not span the same region of the unit circle in the complex plane as e ± i q L . Under this condition, Equation (20b) leads to A = 0 on one side of the Brillouin zone but A 0 on the other side. B + = 0 and B = 0 on the side of the Brillouin zone, where A + = 0 and A = 0 . This result indicates that the two counterpropagating eigenmodes have different modal structures with defined handedness levels, which can be interpreted as a pseudo-spin associated with the propagation direction. In the present scenario, only one pseudo-spin state is allowed to propagate in the lattice in a specific direction, associated with spin momentum locking.
It is interesting to note that for the SSH model ( n = p = 1 ), sin k 2 d 2 = sin k 2 a , and since in the denominator of Equation (23) the quantity δ + ( 2 ) δ ( 2 ) sin k 2 a , both terms cancel each other; A + cannot vanish anymore at sin k 2 ( d 2 = a ) = 0 . The component of the transfer matrix t 12 0 for the SSH system. Therefore, the SSH model does not lead to zero-amplitude Bloch modes, as was noted at the end of Section 2.2. This is also true for a superlattice with a segment of harmonic lattice 1 extending beyond one lattice spacing but with a segment of harmonic lattice 2 limited to one lattice spacing (n > 1 and p = 1). The reverse (n = 1 and p > 1) is not true, thereby illustrating the importance of the choice of origin in the topology (Berry or Zak phase) of vibrational modes in superlattices [14,18].
From a physical point of view, if λ e ± i q L in one region of q, the only way the displacement can take the form of a Bloch wave with plane wave factor e i q N L over the entire Brillouin zone (for all, q):
u A , m N = e i q N L A + e i k 1 m a + A e i k 1 m a e i ω t  
u B , l N = e i q N L B + e i k 2 l a + B e i k 2 l a e i ω t
is by setting the terms in parentheses to zero. For any value of m or l, this implies A ± = 0 and B ± = 0 .

2.6. Reciprocity Condition

Reciprocity implies that the signal received by a detector emitted by a vibration source remains the same upon the interchange of the source and the receiver. Since the considered system is linear, time-invariant, and bias-free, we expect reciprocity to be satisfied. In order to investigate reciprocity in the one-way propagation scenario, we focus on the continuum limit of the superlattice. We consider a finite continuous superlattice sandwiched between a source layer “S” and a detection layer “D” with impedances that may differ from those of materials 1 and 2. We calculate the transmission coefficient of the finite superlattice using transfer matrices, as is achieved in [19]. The transmission coefficient is given by:
T = 4 Z D Z S Z D Z 1 b Z 1 Z S c 2 + d + Z D Z S a 2
where:
a = λ μ 2 sin N θ sin θ + cos N θ
b = σ sin N θ sin θ
c = ζ sin N θ sin θ
d = λ μ 2 sin N θ sin θ + cos N θ
In these equations, we have:
λ = cos k 1 d 1 cos k 2 d 2 F sin k 1 d 1 sin k 2 d 2
μ = cos k 1 d 1 cos k 2 d 2 1 F sin k 1 d 1 sin k 2 d 2
σ = sin k 1 d 1 cos k 2 d 2 + F cos k 1 d 1 sin k 2 d 2
ζ = sin k 1 d 1 cos k 2 d 2 1 F cos k 1 d 1 sin k 2 d 2
and cos θ = λ + μ 2 . N is the number of unit cells in the superlattice. In Equation (29), Z 1 , Z D , and Z S are the impedances of the type 1, detection, and source layers, which are chosen to be equal for illustrative purposes; however, varying Z D and Z S will change the results quantitively but not qualitatively. Here, instead of interchanging the source and detector, we swap media 1 and 2 by interchanging the indices 1 and 2 in Equation (29). In the long wavelength limit, we have k 1 = ω c 1 and k 2 = ω c 2 , as well as c 1 = a K 1 and c 2 = a K 2 . In Figure 2, we present the transmission coefficient as a function of the frequency for two superlattices.
The two transmission coefficients are identical for all frequencies, confirming that reciprocity is achieved. In Figure 2, the gaps correspond to the usual Bragg scattering. These gaps arise when λ + μ 2 > 1 , which corresponds to non-propagative modes.
Each band is composed of N resonances. We also note that the lower envelopes of the transmission coefficient of the 2nd and 4th transmission bands are asymmetric, while the 1st, 3rd, and 5th bands, which have conventional topologies, exhibit a more symmetrical behavior. We note in Figure 3a that the lower envelope of the transmission coefficient of the second band approaches one near ω c 1 d 1 = 2.618 (that is, k 2 d 2 = 1.2 × 2.618 = π, which corresponds to the condition sin k 2 d 2 = 0 ). The amplitude of the backward-propagating wave equals zero under this condition. The same condition is satisfied for the fourth band ω c 1 d 1 = 5.236 , where k 2 d 2 = 2 π . This behavior arises for the even-numbered bands in the transmission spectrum. These bands are known to correspond to dispersion bands of the infinite superlattice with a non-conventional topology, for which the Berry or Zak phase is equal to π. The odd-numbered bands are associated with conventional topologies and zero Berry or Zak phases. This type of asymmetry in the transmission of topologically non-conventional bands has been observed theoretically for electronic waves in finite semiconductor superlattices [20] and electromagnetic superlattices [18].
In Figure 3, we can clearly see that the asymmetry in the dispersion of the transmission of the even-numbered bands is associated with the Fabry–Pérot resonance condition of medium 2 in a single unit cell, sin k 2 d 2 = 0 , where the transmission coefficient is equal to one. This is exactly the condition for which, in an infinite superlattice, spin-selective one-way propagation occurs, as imposed by the requirement that the waves take the form of Bloch waves, i.e., the waves satisfy translational invariance. Note that for a single unit cell, gaps exist that again correspond to the condition λ + μ 2 > 1 .
To shed more light on the relationship between the behaviors of the finite system and the infinite system, below we give the transfer matrix across one unit cell derived in [19]. The displacement field in the two media 1 and 2 of a unit cell is expressed as a general ansatz formed as a linear combination of transmitted and reflected waves:
U 1 ( x ) = A + e i k 1 x + A e i k 1 x e i ω t
U 2 ( x ) = B + e i k 2 x + B e i k 2 x e i ω t
Here, x is the coordinate along the direction perpendicular to the layers forming the superlattice. There is no specific constraint on the form of the amplitudes A + , A , B + , and B .
The stress is expressed as:
S 1 ( x ) = i ω Z 1 A + e i k 1 x A e i k 1 x e i ω t  
S 2 ( x ) = i ω Z 2 B + e i k 2 x B e i k 2 x e i ω t
These definitions allow us to introduce the two-component vector for j = 1, 2:
W j x = U j ( x ) S j ( x ) = e i k j x e i k j x i ω Z j e i k j x i ω Z j e i k j x A +   o r   B + A o r   B
Here, we define H j x = e i k j x e i k j x i ω Z j e i k j x i ω Z j e i k j x . In a stack of unit cells, using the conditions of the continuity of displacement and stress, one can relate W 1 x = d 1 + d 2 right after a unit cell to W 1 x = 0 right before a unit cell using the following relation:
W 1 x = L = d 1 + d 2 = M W 1 x = 0
where:
M = λ 1 Z 1 ω σ Z 1 ω ζ μ
The symbol in Equation (34) was defined earlier. We remark that there is no constraint imposed on the ansatz of Equation (30a,b). For the periodic infinite superlattice and a unit cell labelled, N, we must use a continuous version of the Bloch wave ansatz given by Equation (28a,b):
U 1 N ( x ) = e i q N L A + p e i k 1 ( x N L ) + A p e i k 1 ( x N L ) e i ω t     for   N L < x < N L + d 1
U 2 N ( x ) = e i q N L B + p e i k 1 ( x N L d 1 ) + B p e i k 1 ( x N L d 1 ) e i ω t   for   N L + d 1 < x < N + 1 L
This ansatz achieves translational periodicity, as it represents Bloch waves. In Equation (35a,b), the upper script “p” stands for periodic. This is necessary as the periodicity redefines the amplitude parameters. Comparing Equation (30a,b) and Equation (35a,b), we can establish constraints that would be imposed by translational invariance on the parameters A + , A , B + , and B , namely:
A + A N = e i k 1 N L 0 0 e i k 1 N L e i q N L A + p A p
B + B N = e i k 1 ( N L + d 1 ) 0 0 e i k 1 ( N L + d 1 ) e i q N L B + p B p
These are actually unitary transformations, which rotate the representation amplitudes of the periodic system into the representation of the amplitudes in Equation (30a,b).
By using Equation (32) and inserting Equation (36a,b) into Equation (33), i.e., N = 1 for the left hand side of the equal sign and N = 0 for the right hand side, one obtains the following relation:
t 11 t 12 t 21 t 22 e i q L 1 0 0 1 A + p A p = 0 .
with:
T = t 11 t 12 t 21 t 22 = H 1 1 ( 0 ) M H 1 ( 0 )
After some algebraic manipulation, we find:
t 11 = λ + μ 2 + i σ ζ 2 ,   t 12 = λ μ 2 + i σ + ζ 2 ,   and   t 21 = t 12 *   and   t 22 = t 11 * .
One can verify algebraically that Equation (37) is identical to Equation (19). We can also observe that the determinants of the unimodular matrices M and T are identical and given by λ μ σ ζ = 1 . When solving the eigenvalues, ε of the matrix M gives a characteristic equation of ε 2 λ + μ + 1 = 0. The finite superlattice takes a propagative eigenmode when λ + μ 2 1 . By setting cos θ = λ + μ 2 , the characteristic equation takes the form of the dispersion relation of the infinite superlattice (Equation (10)) in the limits of the long wavelength, provided one identifies θ with qL. Equations (37) and (19) achieve translational invariance and lead to one-way propagation (zero amplitudes for q < 0 and non-zero amplitudes for q > 0) when sin k 2 d 2 = 0 . When dealing with the finite superlattice, the unconstrained matrix M is diagonalized with the eigenvalues ε + = e i θ and ε = e i θ , which account for forward and backward propagation and are associated with the different pseudo-spins of the eigenvectors. The diagonalized matrix is then brought to the power N to calculate the transmission coefficient of Equation (29) for the finite lattice with N unit cells. The wave function of the finite superlattice will try to approach the Bloch wave solutions. Under the condition sin k 2 d 2 = 0 , the approximate Bloch wave will make the contribution of the backward-propagating wave (i.e., e i θ ) small compared to the forward wave, thereby leading to high values of the transmission coefficient. One-way propagation in the infinite superlattice arises from the combined effects of the constraint of the translational invariance of the wave function and of a Fabry–Pérot resonance condition in the superlattice unit cell.

3. Conclusions

One-dimensional discrete phononic superlattices composed of alternating segments of two different harmonic crystals exhibit Bloch modes with specified handedness, with a finite amplitude for one direction and zero amplitude for the opposite direction. By using the transfer matrix method, i.e., solving the scattering of acoustic waves by the superlattice, we showed that the one-way transport nature of these eigenwaves arises from the periodicity of the system, even when the eigenvalues of the transfer matrix do not span the complete unit circle in the complex plane. The one-way propagation behavior of the discrete superlattice approaches that of a continuous superlattice in the long wavelength limits. A necessary condition for one-way propagation in our superlattice model is that the second harmonic crystal segment possesses a spatial extent that exceeds the lattice parameter. This condition enables Fabry–Pérot-type resonances. In the limits of a superlattice composed of single lattice parameter segments, the model reduces to the SSH system, which does not exhibit one-way propagation. In this case, the eigenvalues of the transfer matrix span the complete unit circle in the complex plane and allow Bloch modes with non-zero amplitude to exist over the entire Brillouin zone. By comparing an infinite superlattice and a finite superlattice in the long wavelength limits, we also showed that reciprocity is always achieved. One-way propagation in the infinite superlattice results from the combined effects of translational periodicity, the non-trivial handedness of the associated eigenvector associated with a pseudo-spin, and the Fabry–Pérot resonances of the second medium in the unit cell of the binary superlattice. Phononic superlattices supporting topologically protected acoustic or elastic waves that avoid back-reflections may provide attractive solutions for designing low-loss devices such as bulk or surface acoustic wave devices for telecommunication systems [21]. Another application of current interest is the use of one-way propagation for the management of thermal phonons. Specifically, recent advances in this context have been reported for phonon focusing [22], thermal switching [23], and thermal barriers [24]. Furthermore, binary superlattices that can be manufactured using thin-film deposition methods [25] or colloid-assembled nanocrystal superlattices [26] can serve as platforms for phonon engineering. These offer viable platforms for the experimental investigation of the effects discussed herein.

Author Contributions

Conceptualization, P.A.D. and K.R.; methodology, P.A.D., K.R., A.K. and A.A..; formal analysis, P.A.D., K.R., A.K. and A.A..; writing—original draft preparation, P.A.D. and K.R.; writing—review and editing, A.K. and A.A.; funding acquisition, P.A.D., K.R., A.K. and A.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the US National Science Foundation (NSF) grant #2242925 through the Science and Technology Center New Frontiers of Sound (NewFoS).

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Rasmussen, C.; Quan, L.; Alù, A. Acoustic nonreciprocity. J. Appl. Phys. 2021, 129, 210903. [Google Scholar] [CrossRef]
  2. Nassar, H.; Yousefzadeh, B.; Fleury, R.; Ruzzene, M.; Alù, A.; Daraio, C.; Norris, A.N.; Huang, G.; Haberman, M.R. Nonreciprocity in acoustic and elastic materials. Nat. Rev. Mater. 2020, 5, 667. [Google Scholar] [CrossRef]
  3. Wang, P.; Lu, L.; Bertoldi, K. Topological phononic crystals with one-way elastic edge waves. Phys. Rev. Lett. 2015, 115, 104302. [Google Scholar] [CrossRef]
  4. Nadkarni, N.; Arrieta, A.F.; Chong, C.; Kochmann, D.M.; Daraio, C. Unidirectional transition waves in bistable lattices. Phys. Rev. Lett. 2016, 116, 244501. [Google Scholar] [CrossRef] [PubMed]
  5. Swinteck, N.; Matsuo, S.; Runge, K.; Vasseur, J.O.; Lucas, P.; Deymier, P.A. Bulk elastic waves with unidirectional backscattering-immune topological states in a time-dependent superlattice. J. Appl. Phys. 2015, 118, 063103. [Google Scholar] [CrossRef]
  6. Nassar, H.; Chen, H.; Norris, A.N.; Haberman, M.R.; Huang, G.L. Non-reciprocal wave propagation in modulated elastic metamaterials. Proc. R. Soc. A 2017, 473, 20170188. [Google Scholar] [CrossRef]
  7. Nassar, H.; Xu, X.C.; Norris, A.N.; Huang, G.L. Modulated phononic crystals: Non-reciprocal wave propagation and Willis materials. J. Mech. Phys. Solids 2017, 101, 10–29. [Google Scholar] [CrossRef]
  8. Khanikaev, A.B.; Fleury, R.; Mousavi, S.H.; Alù, A. Topologically robust sound propagation in an angular-momentum-biased graphene-like resonator lattice. Nat. Commun. 2015, 6, 8260. [Google Scholar] [CrossRef]
  9. Fleury, R.; Khanikaev, A.B.; Alù, A. Floquet topological insulators for sound. Nat. Commun. 2016, 7, 11744. [Google Scholar] [CrossRef]
  10. Deymier, P.A.; Runge, K.; Vasseur, J.O. Geometric phase and topology of elastic oscillations and vibrations in model systems: Harmonic oscillator and superlattice. AIP Adv. 2016, 6, 121801. [Google Scholar] [CrossRef]
  11. Su, W.P.; Schrieffer, J.R.; Heeger, A.J. Solitons in Polyacetylene. Phys. Rev. Lett. 1979, 42, 1698. [Google Scholar] [CrossRef]
  12. Huang, H.; Chen, J.; Huo, S. Recent advances in topological elastic metamaterials. J. Phys. Condens. Matter 2021, 33, 503002. [Google Scholar] [CrossRef] [PubMed]
  13. Pérot, A.; Fabry, C. On the application of interference phenomena to the solution of various problems of spectroscopy and metrology. Astrophys. J. 1899, 9, 87. [Google Scholar] [CrossRef]
  14. Hasan, M.A.; Calderin, L.; Runge, K.; Deymier, P.A. Spectral Analysis of Amplitudes and Phases of Lattice Vibrations: Topological Applications. J. Acoust. Soc. Am. 2019, 146, 748. [Google Scholar] [CrossRef] [PubMed]
  15. Camley, R.E.; Djafari-Rouhani, B.; Dobrzynski, L.; Maradudin, A.A. Transverse elastic waves in periodically layered infinite and semi-infinite media. Phys. Rev. B 1983, 27, 7318. [Google Scholar] [CrossRef]
  16. Yariv, A.; Yeh, P. Optical Waves in Crystals: Propagation and Control of Laser Radiation; Wiley: New York, NY, USA, 1984. [Google Scholar]
  17. Deymier, P.A.; Runge, K. Sound Topology, Duality, Coherence and Wave-Mixing: An Introduction to the Emerging New Science of Sound; Springer Series in Solid-State Sciences; Springer: Berlin/Heidelberg, Germany, 2017; Volume 188. [Google Scholar]
  18. Xiao, M.; Zhang, Z.Q.; Chan, C.T. Surface impedance and bulk band geometric phases in one-dimensional systems. Phys. Rev. X 2014, 4, 021017. [Google Scholar] [CrossRef]
  19. Mizuno, S.; Tamura, S.-I. Theory of acoustic-phonon transmission in finite-size superlattice systems. Phys. Rev. B 1992, 45, 734. [Google Scholar] [CrossRef]
  20. Luo, X.; Shi, J.; Zhang, Y.; Niu, Z.; Miao, D.; Mi, H.; Huang, W. Filtering electrons by mode coupling in finite semiconductor superlattices. Sci. Rep. 2022, 12, 7502. [Google Scholar] [CrossRef]
  21. Deymier, P.A.; Runge, K. One-way propagation of topologically non-conventional bulk transverse elastic waves in infinite and finite superlattices: Application to low-loss acoustic wave devices. Appl. Phys. Lett. 2023, 123, 012202. [Google Scholar] [CrossRef]
  22. Anufriev, R.; Ramiere, A.; Maire, J.; Nomura, M. Heat guiding and focusing using ballistic phonon transport in phononic nanostructures. Nat. Commun. 2017, 8, 15505. [Google Scholar] [CrossRef]
  23. Ishibe, T.; Kaneko, T.; Uematsu, Y.; Sato-Akaba, H.; Komura, M.; Iyoda, T.; Nakamura, Y. Tunable Thermal Switch via Order-Order Transition in Liquid Crystalline Block Copolymer. Nano Lett. 2022, 22, 6105. [Google Scholar] [CrossRef]
  24. Ahn, C.; Fong, S.W.; Kim, Y.; Lee, S.; Sood, A.; Neumann, C.M.; Asheghi, M.; Goodson, K.E.; Pop, E.; Wong, H.-S.P. Energy-Efficient Phase-Change Memory with Graphene as a Thermal Barrier. Nano Lett. 2015, 15, 6809. [Google Scholar] [CrossRef]
  25. Kothari, K.; Malhotra, A.; Maldovan, M. Cross-plane heat conduction in III-V semiconductor superlattices. J. Phys. Condens. Matter 2019, 31, 345301. [Google Scholar] [CrossRef]
  26. Yazdani, N.; Jansen, M.; Bozyigit, D.; Lin, W.M.M.; Volk, S.; Yarema, O.; Yarema, M.; Juranyi, F.; Huber, S.D.; Wood, V. Nanocrystal superlattices as phonon-engineered solids and acoustic metamaterials. Nat. Commun. 2019, 10, 4236. [Google Scholar] [CrossRef]
Figure 1. Schematic representation of a one-dimensional discrete superlattice. A periodically repeating unit cell N is composed of equally spaced identical unit masses coupled through linear springs with stiffnesses of K 1 and K 2 . See text for more details.
Figure 1. Schematic representation of a one-dimensional discrete superlattice. A periodically repeating unit cell N is composed of equally spaced identical unit masses coupled through linear springs with stiffnesses of K 1 and K 2 . See text for more details.
Crystals 14 00092 g001
Figure 2. Transmission coefficient T, corresponding to the first 6 bands of two finite superlattices, with N = 20 as a function of the frequency for (a) Z D = Z S = Z 1 . The two superlattices are related by an interchange of layer indices 1 and 2, which is equivalent to interchanging the source S and detector D. We have also chosen k 2 d 2 = 1.2 k 1 d 1 . Each band supports 20 resonances, corresponding to the number of periods of the finite superlattice, N.
Figure 2. Transmission coefficient T, corresponding to the first 6 bands of two finite superlattices, with N = 20 as a function of the frequency for (a) Z D = Z S = Z 1 . The two superlattices are related by an interchange of layer indices 1 and 2, which is equivalent to interchanging the source S and detector D. We have also chosen k 2 d 2 = 1.2 k 1 d 1 . Each band supports 20 resonances, corresponding to the number of periods of the finite superlattice, N.
Crystals 14 00092 g002
Figure 3. The same scheme as in Figure 1. The black solid line is for a single unit cell N = 1, sandwiched between the source and the detector. The arrows mark the location where sin k 2 d 2 = 0 .
Figure 3. The same scheme as in Figure 1. The black solid line is for a single unit cell N = 1, sandwiched between the source and the detector. The arrows mark the location where sin k 2 d 2 = 0 .
Crystals 14 00092 g003
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Deymier, P.A.; Runge, K.; Khanikaev, A.; Alù, A. Pseudo-Spin Polarized One-Way Elastic Wave Eigenstates in One-Dimensional Phononic Superlattices. Crystals 2024, 14, 92. https://doi.org/10.3390/cryst14010092

AMA Style

Deymier PA, Runge K, Khanikaev A, Alù A. Pseudo-Spin Polarized One-Way Elastic Wave Eigenstates in One-Dimensional Phononic Superlattices. Crystals. 2024; 14(1):92. https://doi.org/10.3390/cryst14010092

Chicago/Turabian Style

Deymier, Pierre A., Keith Runge, Alexander Khanikaev, and Andrea Alù. 2024. "Pseudo-Spin Polarized One-Way Elastic Wave Eigenstates in One-Dimensional Phononic Superlattices" Crystals 14, no. 1: 92. https://doi.org/10.3390/cryst14010092

APA Style

Deymier, P. A., Runge, K., Khanikaev, A., & Alù, A. (2024). Pseudo-Spin Polarized One-Way Elastic Wave Eigenstates in One-Dimensional Phononic Superlattices. Crystals, 14(1), 92. https://doi.org/10.3390/cryst14010092

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop