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
and
, respectively; and the spacing between adjacent masses equals
. The lengths of the segments are
and
, where
and
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
and
, the system becomes the Su–Schrieffer–Heeger (SSH) model [
12]. For large values of
and
, in the long wavelength limit, the system approaches a continuous superlattice.
The displacement of a mass,
, in segment 1 of unit cell
N is defined as
. Similarly, the displacement of masses in segment 2 of unit cell
N is labelled
, with
. The masses obey bulk equations of motion in segments 1 and 2 for
and
, namely:
Equation (1a,b) supports plane wave solutions of the form:
The wave number
with
j = 1, 2 is related to the angular frequency via the well-known dispersion relation of infinite harmonic chains:
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 , we consider Bloch wave solutions, i.e., we choose and , where q is the wave number. We remove the upper script N on and when the periodicity is explicitly accounted for in the plane wave term .
The equations of the motion of masses at interfaces I and II are as follows:
Inserting the null terms
into Equation (4a) and
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:
We also introduce the following conditions of the continuity of displacement at I and II:
With these conditions, Equation (5a,b) becomes:
In Equation (7a), we use
. 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:
where
,
,
, and
, 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,
Equation (8) possesses non-trivial solutions when the determinant of the dynamic matrix is equal to zero. This condition leads to the dispersion relation:
This dispersion relation can be reformulated as:
2.2. SSH Limit
In the SSH scenario (
),
and
. Using the dispersion relations given by Equation (3) and basic trigonometric relations, Equation (10) reduces to:
This leads to the usual SSH dispersion relation:
Equation (8) can be manipulated algebraically to obtain a relationship between
and
, namely:
The amplitudes and do not differentiate between positive or negative , since m = 0 and l = 0 in Equation (2a,b).
Employing Equation (12) in the form of
, we recover the complex SSH amplitudes [
14]:
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, ], when . It is worth noting that these amplitudes cannot be zero for any value of the wave number q.
2.4. Amplitudes of the Displacement Field
To find the amplitudes, we use the transfer matrix method. Equation (8) can be rewritten as:
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
and
. By inserting that relation into Equation (16a), we obtain the transfer matrix relating amplitudes in one unit cell to amplitudes in the neighboring cell:
The components of the transfer matrix are given by:
Here,
and
. The star stands for a complex conjugation.
Due to the periodicity of the system, solutions will take the form of Bloch waves:
Combining Equations (18) and (17) leads to the eigenvalue problem [
16]:
If
is an eigenvalue of the transfer matrix, then the amplitudes
and
are given by:
The dispersion relation of Equation (9) can be rewritten as
, so Equation (20b) becomes:
After a number of trigonometric and algebraic manipulations, we obtain:
Furthermore, the amplitude
is given by:
In the continuum limit, these amplitudes simplify to the known expressions [
10]:
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
, such that
, then at
the dispersion relation becomes
,
, and
. This condition corresponds to a Fabry–Pérot resonance of the second medium in the superlattice unit cell [
13].
Since
, the sign of
determines whether
vanishes or not. From Equation (16b),
and
, meaning that
and
. At
, the amplitudes of the Bloch modes may be finite, while at
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
, 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
in the following manner. Let us first find expressions for the eigenvalues of the transfer matrix. These eigenvalues are solutions of the equation
. The Eigen values are given by:
For vibrational modes to be Bloch waves, we require
. According to Equation (23),
= 0 when
. In that case, the eigenvalues reduce to:
or:
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 . This condition is met only for a wave number q located on one side of the Brillouin zone, where can be identified as a sine function. For q on the other side of the Brillouin zone, is identifiable to minus one times a sine function, which has the opposite sign of . In other words, does not span the same region of the unit circle in the complex plane as . Under this condition, Equation (20b) leads to on one side of the Brillouin zone but on the other side. and on the side of the Brillouin zone, where and . 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 (
),
, and since in the denominator of Equation (23) the quantity
, both terms cancel each other;
cannot vanish anymore at
The component of the transfer matrix
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
in one region of
q, the only way the displacement can take the form of a Bloch wave with plane wave factor
over the entire Brillouin zone (for all,
q):
is by setting the terms in parentheses to zero. For any value of
m or
l, this implies
and
.
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:
where:
In these equations, we have:
and
.
N is the number of unit cells in the superlattice. In Equation (29),
,
, and
are the impedances of the type 1, detection, and source layers, which are chosen to be equal for illustrative purposes; however, varying
and
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
and
, as well as
and
. 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
, 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
(that is,
= π, which corresponds to the condition
). The amplitude of the backward-propagating wave equals zero under this condition. The same condition is satisfied for the fourth band
, where
. 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,
, 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
.
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:
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 , , , and .
The stress is expressed as:
These definitions allow us to introduce the two-component vector for
j = 1, 2:
Here, we define
. In a stack of unit cells, using the conditions of the continuity of displacement and stress, one can relate
right after a unit cell to
right before a unit cell using the following relation:
where:
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):
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
,
,
, and
, namely:
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:
with:
After some algebraic manipulation, we find:
One can verify algebraically that Equation (37) is identical to Equation (19). We can also observe that the determinants of the unimodular matrices and are identical and given by . When solving the eigenvalues, ε of the matrix M gives a characteristic equation of = 0. The finite superlattice takes a propagative eigenmode when . By setting , 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 When dealing with the finite superlattice, the unconstrained matrix M is diagonalized with the eigenvalues and , 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 , the approximate Bloch wave will make the contribution of the backward-propagating wave (i.e., ) 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.