Basic Cells Special Features and Their Influence on Global Transport Properties of Long Periodic Structures

Abstract

In this contribution, we address quantum transport in long periodic arrays whose basic cells, localized potentials $U\left(x\right)$, display certain particular features. We investigate under which conditions these “local” special characteristics can influence the tunneling behavior through the full structure. As the building blocks, we consider two types of $U\left(x\right)$s: combinations of either Pöschl–Teller, $U_0/{cosh}^2\left[\alpha x\right]$, potentials (for which the reflection and transmission coefficients are known analytically) or Gaussian-shaped potentials. For the latter, we employ an improved potential slicing procedure using basic barriers, like rectangular, triangular and trapezoidal, to approximate $U\left(x\right)$ and thus obtain its scattering amplitudes. By means of a recently derived method, we discuss scattering along lattices composed of a number, N, of these $U\left(x\right)$s. We find that near-resonance energies of an isolated $U\left(x\right)$ do impact the corresponding energy bands in the limit of very large Ns, but only when the cell is spatially asymmetric. Then, there is a very narrow opening (defect or rip) in the system conduction quasi-band, corresponding to the energy of the $U\left(x\right)$ quasi-state. Also, for specific $U_0$’s of a single Pöschl–Teller well, one has 100% transmission for any incident $E>0$. For the $U\left(x\right)$ parameters rather close to such a condition, the associated array leads to a kind of “reflection comb” for large Ns; $|T_N{\left(k\right)|}^2$ is not close to one only at very specific values of k, when $|T_N{|}^2\approx 0$. Finally, the examples here—illustrating how the anomalous transport comportment in finite but long lattices can be inherited from certain singular aspects of the $U\left(x\right)$s—are briefly discussed in the context of known effects in the literature, notably for lattices with asymmetric cells.

1. Introduction

Regular and repetitive structures are common in the natural world [1,2,3,4], arising from uncountable phenomena [5]. So, often the associated patterns are used to probe, and thus to understand, a vast number of distinct processes [6,7]. Potentially, their organization and related general—in the sense of recurrent—symmetries should determine usual traits and properties [6,8], resulting in rather universal features. In fact, the global characteristics of energy bands in arbitrary crystals, certain cyclical attributes of the natural elements (even allowing the construction of a periodic table) and the replicating trends in biology [8] are all good examples of such an ubiquity.

Fitting perfectly well into the above description is the relevant problem of wave (either elastic, electromagnetic or matter) scattering in periodic lattices [9,10,11]. Nonetheless, there are situations in which the dynamics in determined orderly networks can display unusual comportment. In some instances, this is due to uniquely engineered systems [12,13,14]. One illustration, which is a bit controversial, is the emergence of classical-like chaotic behavior due to quantum scattering in Fibonacci lattices (for a historical account and review of the original literature, see [15]). Another relates to particular (quasi-transparent) scattering states, which are connected to the non-trivial zeros of the Riemann zeta function (and so, with the Riemann’s hypothesis) in the so-called logarithmic chains; refer to [16] and refs. therein.

But anomalous or unusual effects do not appear only when waves are propagating in a too-tailored structure. They can also take place in typical lattices, as nicely discussed and exemplified in refs. [17,18,19,20]. However, there may be a caveat. The mentioned effects might be (theoretically) unveiled and explained only if one considers large enough lattices, when multiple scattering [21] becomes fundamental and is treated exactly or at least in the proper orders of approximation [22,23,24].

Based on a well-established Green’s function approach [25] (for a review, see [26]), very recently a method has been developed [27] to calculate the scattering properties of arbitrarily long lattices formed by N equally spaced building blocks (or cells). The framework fully includes multiple scattering processes, and the final exact expressions are simply dependent on the basic cell reflection, R, and transmission, T, coefficients. Therefore, the results are exact if R and T are exact. Moreover, for N extremely large—e.g., numerically N has been analyzed up to ${10}^{10}$ in [27]—one recovers most of the essential features of the energy bands corresponding to the infinite lattice case.

Employing the above-mentioned protocol, in this contribution we shall discuss the transmission properties across long regular lattices assuming special building blocks. We suppose that symmetric and asymmetric basic cells display particular characteristics. Then, we examine if these “local” features can introduce additional trends in the transport behavior of the whole structure. In general, we find that the near-resonance energies of the cell influence some aspects of the band’s shape, but only when the cells are asymmetric. In this case, we observe a narrow gap opening, a rip, in the band corresponding to the energy of the quasi-state of the cell. This effect, nevertheless, becomes noticeable only for relatively large Ns, i.e., large structures.

To construct such special cells, we consider two instances, both relying on localized Us. The difference is that, in the first, for Pöschl–Teller (second, for Gaussian) potentials, R and T are (are not) known analytically. For the latter, we observe that an usual procedure is to approximate a continuous localized $U\left(x\right)$ by a collection of $\mathcal{N}$ very narrow barriers, often rectangular ones (see Figure 1a). The shape of $U\left(x\right)$ is then reproduced by properly setting the barrier heights. Thus, one can employ distinct techniques, such as the transfer matrix method [28], to compute (in general numerically) $T_U$ and $R_U$. But not uncommonly, $\mathcal{N}$ must be large, especially when $U\left(x\right)$ is a peak-shaped function. We show that the usage of basic triangular and trapezoidal (besides rectangular) barriers, whose scattering coefficients are analytically accessible, can considerably reduce $\mathcal{N}$. Such basic potentials form a composed building block (CBB), allowing good analytic approximations for $T_U$ and $R_U$. Finally, once we have the cell scattering amplitudes, we can study arrays of N equally spaced $U\left(x\right)$s via the method in [27].

The paper is organized as follows. We present a very brief overview of the general approach in [27] and its basic equations in Section 2. We exemplify how to approximate a continuous localized potential by a collection of barriers in Section 3. We present some properties of long lattices in Section 4. We explore anomalous behavior in the transmission probability originated from asymmetric cells in Section 5. Final remarks and conclusions are drawn in Section 6. More technical necessary results are provided in the Appendices.

2. Brief Review on the Approach Used to Study Transport in 1D Periodic Structures

Consider an arbitrary system composed of $\mathcal{M}$ non-overlapping compact support potentials, $U_m$. So, each $U_m$ ($m=1,\cdots ,\mathcal{M}$) is non-null only within a finite spatial interval. The resulting structure is illustrated in Figure 1b. Further, assume that the scattering coefficients of $U_m$ are given by $T_m$ and $R_m^{(\pm )}$. Here, these coefficients are written up to the phases usually associated with either the width (transmission) or the end positions of $U_m$ (for details, see Appendix A). The superscript $(\pm )$ identifies the incident direction for a wave incoming from $x\to \mp \infty$. Obviously, if the potential is symmetric, $R_m^{(+)}=R_m^{(-)}=R_m$.

For the separation between potentials $m=\mathcal{M}-1$ and $m=\mathcal{M}$ being denoted as $d_{\mathcal{M}}$, Figure 1b, the reflection, $R_{1,\mathcal{M}}^{(\pm )}$, and transmission, $T_{1,\mathcal{M}}$, coefficients for the whole array formed by $U_1,U_2,\dots ,U_{\mathcal{M}}$ can be calculated from that formed by $U_1,U_2,\dots ,U_{\mathcal{M}-1}$ through the recurrence relations

(1)$\begin{array}{ccc}\hfill R_{1,\mathcal{M}}^{(+)}& =& R_{1,\mathcal{M}-1}^{(+)}+{\displaystyle \frac{R_{\mathcal{M}}^{(+)}{T}_{1,\mathcal{M}-1}^{2}exp\left[2ik{d}_{\mathcal{M}}\right]}{1-R_{1,\mathcal{M}-1}^{(-)}{R}_{\mathcal{M}}^{(+)}exp\left[2ik{d}_{\mathcal{M}}\right]}},\hfill \\ \hfill R_{1,\mathcal{M}}^{(-)}& =& R_{\mathcal{M}}^{(-)}+{\displaystyle \frac{R_{1,\mathcal{M}-1}^{(-)}{T}_{\mathcal{M}}^{2}exp\left[2ik{d}_{\mathcal{M}}\right]}{1-R_{1,\mathcal{M}-1}^{(-)}{R}_{\mathcal{M}}^{(+)}exp\left[2ik{d}_{\mathcal{M}}\right]}},\hfill \\ \hfill T_{1,\mathcal{M}}& =& {\displaystyle \frac{T_{1,\mathcal{M}-1}^{(+)}{T}_{\mathcal{M}}^{(+)}exp\left[ik{d}_{\mathcal{M}}\right]}{1-R_{1,\mathcal{M}-1}^{(-)}{R}_{\mathcal{M}}^{(+)}exp\left[2ik{d}_{\mathcal{M}}\right]}}.\hfill \end{array}$

Such expressions follow directly from the exact Green’s function method in [25].

Now, suppose a periodic finite lattice with N identical localized potentials, $U\left(x\right)$, where any two neighbor Us are separated by L. In this case, we can simplify the notation by writing $T_{1,\mathcal{M}}=T_N$ and $R_{1,\mathcal{M}}^{(\pm )}=R_N^{(\pm )}$ and the scattering amplitudes of U by $R^{(\pm )}$ and T. Moreover, for any $\mathcal{M}$ we have $d_{\mathcal{M}}=L$. Then, it has been shown [27] that the above recurrence relations lead to (using the fact that always $|R^{(+)}{|}^2=|R^{(-)}{|}^2={\left|R\right|}^2$)

(2)$|T_N{|}^2=|1-{\displaystyle \frac{C_N}{C_{N-1}}}|{\displaystyle \frac{{\left|T\right|}^2}{{\left|R\right|}^2}},{|R_N^{(\pm )}|}^2=|1-{\displaystyle \frac{C_N}{C_0}}{|}^2{\displaystyle \frac{1}{{\left|R\right|}^2}},$

where

(3)$C_N={\displaystyle \frac{(1+C_0{\mathsf{\Gamma }}_{-}){({\mathsf{\Gamma }}_{-}/{\mathsf{\Gamma }}_{+})}^N-(1+C_0{\mathsf{\Gamma }}_{+})}{(C_0-{\mathsf{\Gamma }}_{+}){({\mathsf{\Gamma }}_{-}/{\mathsf{\Gamma }}_{+})}^N-(C_0-{\mathsf{\Gamma }}_{-})}}$

and

(4)$\gamma ={\displaystyle \frac{1+(T^2-R^{(+)}{R}^{(-)})exp\left[2ikL\right]}{iTexp\left[ikL\right]}},{\mathsf{\Gamma }}_{\pm }={\displaystyle \frac{\gamma \pm \sqrt{{\gamma }^2+4}}{2}}.$

The above expressions are exact and have no limitations regarding the values of N. So, the transmission and reflection probabilities for extremely long arrays are easy to obtain from these ready-to-use formula.

Approximating a Rapidly Decaying Potential Through Compact Support Potentials

Suppose a continuous quickly decaying potential, $U\left(x\right)$, i.e., it rapidly tends to zero as $\left|x\right|$ increases. In principle, we can approximate $U\left(x\right)$ as a block of juxtaposed $\mathcal{N}$ barriers, as illustrated in Figure 1a. Then, from our present framework, the scattering coefficients of $U\left(x\right)$, $T_U$ and $R_U^{(\pm )}$ can be computed from $T_{1,\mathcal{M}}$ and $R_{1,\mathcal{M}}^{(\pm )}$ in Equation (1), with $\mathcal{M}=\mathcal{N}$ and all ds set to zero.

We should remark that, in the potential slicing technique, one often assumes a relatively large set of rectangular barriers to approximate $U\left(x\right)$ and then employs the transfer matrix method for the calculations (for a didactic review, see, e.g., Ref. [29]). However, this scheme is considerably simplified and computationally speed up if we allow a more diverse set of barriers, namely, not only rectangular (r) but also triangular (t) and trapezoidal ($\tilde{\mathrm{t}}$) shapes, all presenting exact analytical scattering coefficients; see Appendix B. This drastically reduces $\mathcal{N}$, also making Equation (1) simple to handle.

3. Certain Special $\mathit{U}\mathit{\left(}\mathit{x}\mathit{\right)}$s and Their CBBs Description

For the sake of nomenclature, we will refer to the collection of $\mathcal{N}$ localized potentials closely describing a certain U as the composed building block (CBB) representing $U\left(x\right)$. So, the goal is to have $T_U$ and $R_U^{(\pm )}$ well approximated by the corresponding $T_{CBB}=T_{1,\mathcal{N}}$ and $R_{CBB}^{(\pm )}=R_{1,\mathcal{N}}^{(\pm )}$. Further, unless otherwise explicitly mentioned, we will suppose $\mathrm{m}=\hslash =1$, with $\mathrm{m}$ being the particle mass. Accordingly, the energy is given by $E=k^2{\hslash }^2/2\mathrm{m}=k^2/2$.

3.1. The Gaussian Potential

First, we discuss the Gaussian potential

(5)$U_G\left(x\right)=U_{0}exp\left[-{\displaystyle \frac{{(x-\mu )}^2}{2{\sigma }^2}}\right],$

for which there are no known analytic expressions for the reflection and transmission amplitudes. We approximate $U_G$ in four different ways, considering four sets of CBBs: triangle–triangle (I), triangle–rectangle–triangle (II), triangle–trapezoid–rectangle–trapezoid–triangle (III) and triangle–trapezoidal–trapezoid–rectangle–trapezoid–trapezoid–triangle (IV). The constructions are depicted in Figure 2a. To estimate how similar the CBB structures are to the original Gaussian potential, we can compute the difference in the areas. For $U_G$ with the parameters $\mu =0$, $\sigma =\sqrt{2}/2$, $U_0=1$ and for CBBs with the parameters given in Table 1, the areas of CBB-I, CBB-II, CBB-III and CBB-IV differ from that of $U_G$, respectively, by 8.11%, 6.48%, 3.50% and 2.56%. The parameters in Table 1 (and in the next table) were manually adjusted (in trial-and-error attempts), seeking to minimize the areas’ differences. Of course, automatized numerical procedures would improve the agreement even more. We intend to implement algorithms to do so in a future contribution.

As previously mentioned, we obtain $T_{CBB-I,II,III,IV}$ from the recurrence relations in Equation (1) and from the basic shapes of individual Rs and Ts given in Appendix B (the resulting final analytic expressions are a bit lengthy and are thus not shown here, but they are easy to derive and handle from any algebraic software). In Figure 2b, we compare the transmission probabilities for these CBBs with that of the actual Gaussian calculated via the numerically accurate Wronskian method of Ref. [30,31]. From a direct visual inspection, one realizes that the agreement is very good for CBB-III and CBB-IV and fairly good for CBB-I and CBB-II. As an extra test, we have integrated the curves ${\left|T\left(E\right)\right|}^2$ in the energy region shown in Figure 2b. The difference between the results for the Gaussian and the distinct CBB-I,II,III,IV are, respectively, 2.59%, 1.59%, 0.43%, 0.36%. These small values, especially for the latter two, quantitatively confirm the observed, rather fine, concordance.

3.2. The Pöschl–Teller Potential

Next, we address the Pöschl–Teller potential,

(6)$U_{PT}\left(x\right)={\displaystyle \frac{U_0}{{cosh}^2\left[\alpha (x-c)\right]}},$

which recently has found many applications in condensed matter physics and quantum chemistry [32,33,34,35,36]. The exact quantum scattering amplitudes for $U_{PT}$ are presented in Appendix C. For a full analysis of $U_{PT}$ in terms of a Green’s function approach—the key technique used in [27] to derive the expressions in Section 2—refer to Ref. [37].

Previously, seven elementary barriers were enough to yield a very good approximation for the Gaussian potential (CBB-IV in Figure 2a). So, we again use this same number to construct the CBB for the Pöschl–Teller potential. The resulting CBB, of concrete configuration t$\tilde{\mathrm{t}}$$\tilde{\mathrm{t}}$r$\tilde{\mathrm{t}}$$\tilde{\mathrm{t}}$t, is depicted in the inset of Figure 3. Setting $c=0$, $\alpha =1$ and $U_0=2$, the corresponding parameters for the CBB’s first t$\tilde{\mathrm{t}}$$\tilde{\mathrm{t}}$r basic shapes are (by symmetry, those for the last three barriers, namely, $\tilde{\mathrm{t}}$$\tilde{\mathrm{t}}$t, are akin) $w=1.30,0.70,1.15,0.3$, $U_a=0,0.11,0.47$, $U_b=0.11,0.47,2$ and $\tilde{U}=2$. We notice that the areas of this CBB and of the actual $U_{PT}\left(x\right)$ differ by only 2.92%. Figure 3 compares the transmission probability for the CBB and for the $U_{PT}$ (this latter obtained from the exact expression in Equation (A7)). As expected, the curves are very similar, e.g., when integrated, differing by only 0.36%.

4. Simple Finite Periodic Lattices

We shall now consider “simple” finite periodic lattices. By simple we mean arrays of N equally spaced cells, $U\left(x\right)$s. For these cells, we assume a single symmetric $U_G$ or $U_{PT}$. Also, the distance between two successive Us are such that the potentials have a fairly negligible overlap (for the curious scattering properties of two highly superposed Pöschl–Teller potentials, see [37]).

4.1. The Gaussian Case

For the Gaussian arrays, we set a distance, $\overline{\mu }$, between the cells. As the transmission and reflection amplitudes of $U_G$, we shall consider those of the corresponding CBB-IV. Hence, for $\overline{w}$ denoting the base length of the whole CBB-IV structure, as illustrated in Figure 4a, we define the lattice parameter as $L=\overline{\mu }-\overline{w}$.

Naturally, we should check if by using CBB-IV, instead of the actual $R_G$ and $T_G$, we can still obtain a reliable description for the full lattice transport properties. So, as a verification, we compare $|T_N{\left(k\right)|}^2$ from the expression in Eqution (2), setting $R=R_{CBB-IV}$ and $T=T_{CBB-IV}$, with numerical schemes. For $N=2$ and $N=3$, the traditional slicing method (based on the transfer matrix procedure) employing $\mathcal{N}=88$ rectangular barriers for each Gaussian has been implemented in [38], leading to very accurate results. Figure 4b,c show plots for the transmission probability as a function of k. As one can see, the agreement between the two approaches is very good, but in our case computationally rather inexpensive. In this way, in the following we analyze a lattice with $N={10}^4$ Gaussians, a configuration that, to the best of our knowledge, has not been addressed previously in the literature for this type of cell.

We further remark that while the present analysis might be feasible using other frameworks, such as the transfer or scattering matrix approaches (see, e.g., [39,40,41]), the current method is computationally more efficient. This fact stems from Equation (3), where the number of cells only enters as an exponent of certain quantities, instead of representing the number of matrices to be multiplied (for more details, see [27]). Also, although we treat the Gaussian potential as a CBB, we emphasize that the general expressions for the transport along the full lattice are exact and do not rely on iterative procedures.

For our finite lattices of $N={10}^4$ equal Gaussian barriers, we discuss four distinct configurations. The corresponding cells display the following parameters pairs $(U_0,\sigma )$; see Equation (5): (4, 0.1), (2, 0.2), (1, 0.4), (0.66, 0.6). Then, the cells of the four cases have the same area. The CBB-IVs used to model the Gaussian barriers are depicted in Figure 5a, with the parameters listed in Table 2. The distances between successive barriers, $\overline{\mu }$, are chosen, such that $L=1$; see Figure 5b. The resulting transmission probabilities as a function of k for the distinct arrays are shown in Figure 5c–f. The associated individual $|T_{CBB-IV}{\left(k\right)|}^2$ is also shown. Since the areas of the different Gaussians are equal and always $L=1$, effectively we have a direct relation between $U_0$ and $\sigma$, and thus we can assume $U_0$ as the only free parameter. Hence, in Figure 5c–f the variation of $U_0$ is simply equivalent to re-scaling the ks, explaining why the quasi-band structures in the four instances are similar, only with distinct widths. One way to picture this re-scaling is to draw a parallel with a tight-binding model. In such a case, the “potential strength” (independent on area) is proportional to the effective hopping integral, which, in turn, is proportional to the band width [42,43]. So, as $\sigma$ decreases, the band widths increase accordingly. Further, notice that for an array of Gaussians the function $|T_{CBB-IV}{\left(k\right)|}^2$ acts as the envelope of the centers of the allowed quasi-bands [44]. This complies with similar features of lattices formed by delta and triangular barriers [27].

We lastly remark that certain works addressing semiconductor superlattices assume Gaussian profiles [45,46,47,48,49]. For instance, for the Gaussians being approximated by rectangular barriers, a photovoltaic device in which a Gaussian superlattice is inserted into a GaAs solar cell has been investigated [50]. Similar constructions have been considered in the analysis of the electronic properties of graphene-based superlattices [51] as well as of phonon tunneling in semiconductor heterostructures [52]. Our present approach would be a very valuable tool for all these problems, permitting one to treat a much larger number of Gaussian potentials, eventually leading to a more realistic description of the aforementioned devices.

4.2. The Pöschl–Teller Case

Differently from Section 4.1, for $U_{PT}$ the scattering coefficients are known analytically; see Appendix C. However, there is a caveat. The Pöschl–Teller is a rapidly decaying but not a compact support potential. Therefore, to be able to use the results of Section 2, a few simple modifications are necessary. In fact, how to proceed in such a context has been fully discussed in Ref. [37]. The first change is to introduce extra phases for $R^{(\pm )}$ and T, and the second is to consider a k-dependent L, both to be implemented in the auxiliary Equation (4).

In particular, the Pöschl–Teller is one of the potentials investigated in detail in [37]. So, here we just use the expressions derived in [37]. For $L\to L\left(k\right)$ in Equation (4), we have $U_{PT}\left(x\right)=U_0/{cosh}^2\left[\alpha (x-c)\right]$ (with $ϵ=U_0/E=2U_0/k^2$), and

(7)$L\left(k\right)={\displaystyle \frac{1}{\alpha }}\left\{\begin{array}{cc}2I\left(\alpha c\right)-ln\left[\right|ϵ-1\left|\right],\hfill & ϵ\ne 1,\hfill \\ 2ln[cosh[\alpha c\left]\right],\hfill & ϵ=1,\hfill \end{array}\right.$

where

(8)$I\left(\xi \right)=ln\left[sinh\left[\xi \right]+\sqrt{{cosh}^2\left[\xi \right]-ϵ}\right]+{\displaystyle \frac{\sqrt{ϵ}}{4}}ln\left[{\left({\displaystyle \frac{\sqrt{ϵ}sinh\left[\xi \right]-\sqrt{{cosh}^2\left[\xi \right]-ϵ}}{\sqrt{ϵ}sinh\left[\xi \right]+\sqrt{{cosh}^2\left[\xi \right]-ϵ}}}\right)}^2\right].$

The above $L\left(k\right)$ implies that the width between the centers of two consecutive $U=U_{PT}$ in the lattice is $2c$ and is thus akin to the distance $\overline{\mu }$ depicted in Figure 4a for the Gaussians. The extra phases multiplying the transmission and reflection coefficients in Equation (4) are given in Appendix C.

Interestingly, when $U_0<0$ and $1+8m|U_0|/{\hslash }^2{\alpha }^2={(2n+1)}^2$ for $n=0,1,2,\dots$, there is full transmission across the Pöschl–Teller well for any incident $k>0$. For example, for $U_0=-1$ and $\alpha =1$, Figure 6a shows that the transmission probability is always one. Likewise, one would obtain $100\%$ transmission along a whole array of these wells, regardless their number, N. Figure 6a also displays $|T_{PT}{|}^2$ for single wells with the same $\alpha$, but for depths that are slightly different from $-1$, namely, $U_0=-0.99$ (1% shallower), $U_0=-1.089$ (10% deeper) and $U_0=-0.891$ (10% shallower). As expected, the further the $|U_0|$s are from the mentioned special values, the larger the necessary onset for k to make $|T_{PT}{\left(k\right)|}^2\approx 1$.

For the three cases with $U_0\ne -1$ in Figure 6a, we suppose finite periodic lattices with $N={10}^4$ and the separation between two successive wells equal to $2c=14$ (inset of Figure 6a). The resulting $|T_N{|}^2$ are shown in Figure 6b–d. We point out that as $U_0$ approaches $-1$, the forbidden quasi-bands tend to become comparatively narrower.

To better appreciate this last behavior, we present in Figure 7 the plot of $|T_N{|}^2$ versus k for an extremely long lattice of $N={10}^7$ Pöschl–Teller wells with $U_0=-0.999$, hence with a depth much closer to the reflectionless case of $U_0=-1$ than the examples in Figure 6b–d. The other parameter values are those of Figure 6. Notice that the emerging forbidden quasi-bands are anomalously narrow (essentially spikes), a kind of reflection comb. This pattern continues for larger ks, although we have not analyzed the eventual threshold $k_{rc}$ for which the comb-like comportment is lost (this will be the subject of a future contribution). In Figure 7, the average separation between instances of $|T_N|\approx 0$ is $\Delta k=0.23$.

Additionally, we mention that, even for hugely large but finite Ns, typically the allowed quasi-bands are not characterized by perfect transmission. As nicely shown in [44], one can still have fluctuations from $|T_N{|}^2=1$. This is indeed seen in Figure 6 as the “dark regions” in the plots of $|T_N{|}^2$. On the other hand, for the allowed quasi-bands in Figure 7, we practically do not see a departure from $|T_N{|}^2\approx 1$. Thence, at least in this regard the array in Figure 7 already exhibits features of an actual infinite lattice.

5. Anomaly in the Lattice Transmission Induced by Asymmetric Basic Cells

The recurrence relations in Equation (1) allow one to construct a multitude of distinct symmetrical and asymmetrical building blocks (SBB and ABB). This should lead to a great diversity of scattering processes taking place in the associated arrays, consequently impacting their global transport features. In this last section, we shall exemplify such a type of phenomenology, contrasting symmetric and asymmetric cells. So, we consider lattices having building blocks formed by two barriers, left and right, whose centers are a distance, d, apart. Consequently, the BBs are spatially symmetric (asymmetric) if these two barriers are equal (different). We notice that the transport properties of lattices having, as a cell, a single asymmetric potential, like trapezoidal and triangular barriers, have already been discussed elsewhere (for an overview, see, e.g., [27]). However, as far as we know, the particular effect we characterize here is not directly comparable with the ones in these other studies.

We start by supposing an ABB composed of two Gaussians of parameters $U_0=4,\sigma =0.1$ and $U_0=0.66,\sigma =0.6$, moreover, with $d=1$. The calculations are then approximated by our previous CBB-IVs (the comparison between the double Gaussians and the double CBB-IVs, as well as the parameters for the latter, are presented, respectively, in the inset of Figure 8a and in Table 2). The transmission probability through the ABB, denoted as $|T_{N=1}{|}^2$, is shown in Figure 8a. We observe that around $k=0.9$ there is a peak where $|T_1{|}^2=1$. In other words, $k_{res}\approx 0.9$ is a quasi-state resonance, allowing a 100% transmission through the ABB structure.

For the lattices formed by this ABB, we set the distance between the cells (constituted by two CBB-IVs) to be 1. In Figure 8b–f we depict the resulting $|T_N{|}^2$s versus k for N equal to 5, 10, 20, 50 and ${10}^4$. As expected, as N increases we observe the formation of quasi-band structures for the associated systems. Nonetheless, another effect also arises. The transmission probability for $k\approx k_{res}$ gradually shifts away from one, eventually vanishing for large Ns. Figure 8g represents a blow-up of Figure 8f in a k interval about $k_{res}$. It highlights the emergence of a type of defect or rip in the allowed quasi-band corresponding to $k_{res}$. Hence, instead of permitting a full propagation along the array when the incident wave has $k\approx k_{res}$, the original quasi-bond state of the ABB cell turns into a trapped state of the finite lattice—provided N is large enough. Although not explicitly displayed here (but see below), we report that, for Gaussian SBBs, we have not observed the above-mentioned rips in the lattice bands around the cell $k_{res}$s.

To investigate if such a phenomenon is not only due to the particularities of Gaussians, we also suppose the Pöschl–Teller potential. To simplify the analysis, we set the same $\alpha =1$ and $c=5$ for both cell barriers (recall that the centers are $2c=10$ apart), defining $U_{0,1}$ (left) and $U_{0,2}$ (right) as their heights. Obviously, for SBB (ABB) we have $U_{0,1}=U_{0,2}$ ($U_{0,1}\ne U_{0,2}$). The calculations for $|T_N{|}^2$ follow the procedure in the previous section. We only mention that to obtain $R^{(\pm )}$ and T for the Pöschl–Teller SBB and ABB—from the recurrence relations in Equation (4)—we set $L\to L_1\left(k\right)/2+L_2\left(k\right)/2$ (refer to Figure 9a), with $L_j$ given by Equation (7). Figure 9b,c display $|T_1{\left(k\right)|}^2$ for SBBs with $U_{0,1}=U_{0,2}=0.9$ and $U_{0,1}=U_{0,2}=0.8$. Resonance $k_{res}$s are observed for both SBBs. The transmission probabilities for the related lattices with $N={10}^4$ cells are also shown. Note that no defects (or rips) are identified in the allowed quasi-bands corresponding to these $k_{res}$s.

The situation is distinct for ABBs. In Figure 9d, we have similar plots, but now for a basic cell with $U_{0,1}=0.9$ and $U_{0,2}=0.8$. As in Figure 8, the defects do appear for the quasi-bands corresponding to the ks for which $|T_1{|}^2\approx 1$. For example, for $k\approx 1.27,1.55,1.83$ it reads $|T_1{|}^2\approx 0.97,0.99,1.00$ and we clearly see a kind of defect in the quasi-bands related to such ks. For an ABB with $U_{0,1}=0.9$ and $U_{0,2}=0.6$, the results are shown in Figure 9e, displaying a somehow similar phenomenology. However, for $k\approx 1.25,1.53,1.82$, then $|T_1{|}^2=0.82,0.95,0.99$. So, for this ABB the transmission probability local peaks tend to be smaller. As a consequence, the defects are more like a very narrow forbidden quasi-band than a rip.

The present analysis is also repeated for cells formed by two rectangular barriers in Appendix D. We find exactly the same type of behavior.

A natural question is why the defects arise only when the double barrier cells are asymmetric. This can be explained by considering the multiple scattering processes taking place within each cell and in between the different cells. First, the sole cell resonances, $k_{res}$s, are a direct consequence of the multiple scattering occurring between its left and right barriers. But if the cells are symmetric, once placed in an array, the full structure in fact can be viewed as a set of $2N$ equal barriers, and in a sense the identity of the individual double barrier cells is lost. Therefore, any emergent undulatory process is a consequence of the multiple scattering taking place for all the localized $2N$ potentials, each contributing at the same footing. Even if the distances between two cell barriers, d, and two consecutive cells, L, are distinct, this only introduces two length scales (so two characteristic phases), but the single barriers scattering coefficients are all the same along the array. Conversely, for asymmetric cells, effectively we have two distinct multiple scattering processes: that in each cell and that among the distinct cells along the whole lattice. Thus, in general, the $k_{res}$s are not simply washed out by organizing the isolated cells into a network. As a result, the corresponding quasi-state resonances must influence the full structure comportment, and by translational invariance symmetry yield lattice trapping states such that $|T_N\left(k_{res}\right){|}^2\approx 0$.

As a parallel, the Su–Schrieffer–Heeger (SSH) model [42,53,54], widely employed in the study of polymers such as polyacetylene, exhibits a similar gap opening behavior. This tight-binding model depicts a single spinless electron on an 1D lattice with two (distinct) site unit cells, simulating a lattice distorted by the Peierls instability [55,56]. Thus, the electrons have only one degree of freedom, hopping between sites. The gap opening is then caused by the difference between the interior and exterior hopping potentials, with the former (latter) referring to the dynamics within the unit cell (connecting adjacent unit cells). Hence, phenomenologically the two hopping potentials in the SSH correspond to the two distinct multiple scattering mechanisms unveiled in our work.

Finally, we shall mention that anomalous behavior related to the transport properties of a lattice often arises from asymmetrical configurations. For instance, numerical analysis of 1D photonic arrays has revealed that asymmetric defects significantly affects beam propagation [57]. Also, when considering non-periodic 1D lattices with momentum conservation, an asymmetric potential can induce normal heat conduction due to a finite-size effect, with such an anomaly emerging in the thermodynamic limit [58]. Analogous results have been observed in three-dimensional networks. As an example, we cite the shape memory of alloy cellular lattices. Micro-structural imperfections within the cells lead to spatial asymmetries, which significantly influence the mechanical response of the entire lattice [59].

6. Final Remarks and Conclusions

In the present work, we have analyzed the influence of the local features of certain basic cells on the transmission properties of finite, but long, periodic arrays. So, based on a framework developed in [27], we have developed a simple scheme to generate distinct building blocks (the cells) with different structures and to calculate their resulting reflection, $R^{(\pm )}$, and transmission, T, coefficients. It is worth mentioning that for the many distinct potentials (in particular those addressed here), which could be used to construct the cells, the proposed approach tends to be computationally more efficient than the potential slicing technique (based on rectangular barriers and the transfer matrix method). As case studies, we have discussed combinations of Gaussian and of Pöschl–Teller potentials, obtaining very good numerical results.

Importantly, we have unveiled certain anomalous behavior for transport in structures induced by their cells’ rather special characteristics. For lattices composed of Pöschl–Teller wells, when the isolated potentials have their parameter values close to the conditions of 100% transmission (a remarkable property of $U_{PT}$), the transport profile along the array resembles an unusual "reflection comb". In other words, $|T_N{|}^2$ is not one only for very determined values of k, when it practically vanishes.

A second anomalous effect was observed in periodic very long structures formed by double barrier cells. For asymmetric shapes, the allowed quasi-bands corresponding to the $k_{res}$s display a kind of defect, like a rip, which becomes narrower and narrower for $|T_N\left(k_{res}\right){|}^2$ closer and closer to 1. Nevertheless, such a phenomenon does not take place if the cells are spatially symmetric. The proper physical reasons for this have been discussed.

Since the segmentation of different continuous localized potentials into rectangular barriers is a common procedure to treat contexts like tunneling through gate oxides [60,61] and quantum dots [62], our protocol here could be an important new tool to deal with these same types of problems. Moreover, well-tailored CBBs could mimic relevant continuous potentials, e.g., approximating Lennard–Jones-like and Morse-like Us, which are often employed to describe inter-atomic effective interactions [63].

We finally mention that other possible applications for the present general method relate to looking for anomalous behavior in very promising systems (aimed at optical and thermal devices) like 1D photonic crystals [64,65,66] and low-dimensional solid structures such as nanotubes and nanowires [67].

Author Contributions

Conceptualization, M.G.E.d.L.; methodology, L.R.N.O.; software, L.R.N.O.; validation, L.R.N.O.; formal analysis, L.R.N.O.; investigation, L.R.N.O. and M.G.E.d.L.; writing—original draft preparation, L.R.N.O.; writing—review and editing, L.R.N.O. and M.G.E.d.L.; supervision, M.G.E.d.L.; project administration, M.G.E.d.L.; funding acquisition, M.G.E.d.L. All authors have read and agreed to the published version of the manuscript.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

The data supporting this study’s findings are available from the corresponding author upon request.

Acknowledgments

We would like to thank L. S. Schulman and M. Koehler for fruitful and elucidating discussions.

Conflicts of Interest

The authors declare no conflicts of interest.

Abbreviations

The following abbreviations are used in this manuscript:

CBB	Composed building block
WM	Wronskian method
TM	Transfer matrix method
SBB	Symmetrical building block
ABB	Asymmetrical building block

Footnotes

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.

Figures and Tables

Figure 1. (a) Schematics of a continuous localized potential approximated by a set of [Forumla omitted. See PDF.] elementary (here rectangular) barriers. The darker gray potential indicates the n-th barrier, with its reflection and transmission coefficients shown. (b) An array of [Forumla omitted. See PDF.] arbitrary localized, compact support potentials.

Figure 2. (a) Distinct CBBs that approximate a Gaussian potential: CBB-tt (I), CBB-trt (II), CBB-t[Forumla omitted. See PDF.]r[Forumla omitted. See PDF.]t (III) and CBB-t[Forumla omitted. See PDF.][Forumla omitted. See PDF.]r[Forumla omitted. See PDF.][Forumla omitted. See PDF.]t (IV). Here, CBB-abc… means that the CBB is formed, from left to right, by the sequence a, b, c, … of juxtaposed barriers; moreover, t, [Forumla omitted. See PDF.] and r stand, respectively, for triangular, trapezoidal and rectangular shapes. The parameters are those in Table 1. (b) As a function of incident energy, the [Forumla omitted. See PDF.]s of CBBs III and IV are compared with the Gaussian potential actual transmission probability, calculated with the Wronskian method (WM) in Ref. [30,31]. The agreement is very good for both CBBs. In the inset is the same type of comparison, but for CBBs I and II.

Figure 3. The exact probability transmission versus k for the Pöschl–Teller potential compared with that of the CBB (whose composition is shown in the inset). For the parameter values, see main text.

Figure 4. (a) Schematics of two successive Gaussians (dashed curves) belonging to an array of N localized barriers. As indicated, each Gaussian is approximated by a proper CBB-IV. The full transmission probability versus k for the case of (b) [Forumla omitted. See PDF.] and (c) [Forumla omitted. See PDF.], calculated from the present approach using the CBB-IV and from the transfer matrix method (TM) in [38] (the corresponding curves have been digitalized directly from Ref. [38]). Each Gaussian reads [Forumla omitted. See PDF.]; in addition, [Forumla omitted. See PDF.]. For the first four BBs of the CBB-IV, namely t[Forumla omitted. See PDF.][Forumla omitted. See PDF.]r, the parameters are [Forumla omitted. See PDF.], [Forumla omitted. See PDF.], [Forumla omitted. See PDF.] and [Forumla omitted. See PDF.]. To align with Ref. [38], specifically for this example, we set [Forumla omitted. See PDF.].

Figure 5. (a) The four CBB-IVs that model the continuous Gaussian potentials, forming arrays of [Forumla omitted. See PDF.] barriers. (b) Illustration of two successive Gaussians composing a lattice in the cases of [Forumla omitted. See PDF.] and [Forumla omitted. See PDF.]. The parameters are always chosen such that [Forumla omitted. See PDF.]. Transmission probabilities as a function of k for the arrays with (c) [Forumla omitted. See PDF.], (d) [Forumla omitted. See PDF.], (e) [Forumla omitted. See PDF.], (f) [Forumla omitted. See PDF.]. The dashed curves represent the transmission probability for the corresponding single CBB-IV.

Figure 6. (a) Transmission probability as a function of k for a single Pöschl–Teller well with [Forumla omitted. See PDF.] and [Forumla omitted. See PDF.] equal to [Forumla omitted. See PDF.], [Forumla omitted. See PDF.], [Forumla omitted. See PDF.] and [Forumla omitted. See PDF.]. For the three latter cases, being the basic cells of finite periodic lattices with [Forumla omitted. See PDF.] and [Forumla omitted. See PDF.] (which is the separation between the centers of two successive wells; inset in (a)), the corresponding [Forumla omitted. See PDF.] versus k plots are presented in (b–d).

Figure 7. Similar to Figure 6b, but for [Forumla omitted. See PDF.] and [Forumla omitted. See PDF.]. The extremely narrow (basically spikes) forbidden quasi-bands occur around the wavenumbers [Forumla omitted. See PDF.], [Forumla omitted. See PDF.], [Forumla omitted. See PDF.], [Forumla omitted. See PDF.], [Forumla omitted. See PDF.], [Forumla omitted. See PDF.], [Forumla omitted. See PDF.], [Forumla omitted. See PDF.], [Forumla omitted. See PDF.]. Apart from for the spikes, there are almost no fluctuations from [Forumla omitted. See PDF.].

Figure 8. (a) For the Gaussian ABB discussed in the main text (and approximated by two CBB-IVs, inset) of parameters [Forumla omitted. See PDF.] (left, the red graphic in the inset) and [Forumla omitted. See PDF.] (right, the blue graphic in the inset) barriers, the resulting transmission probability, [Forumla omitted. See PDF.], shown as a function of k. For the N finite periodic lattices, the corresponding [Forumla omitted. See PDF.] are displayed in (b) [Forumla omitted. See PDF.], (c) [Forumla omitted. See PDF.], (d) [Forumla omitted. See PDF.], (e) [Forumla omitted. See PDF.], (f) [Forumla omitted. See PDF.]. In all cases, the distance between successive cells (formed by two CBB-IVs) are equal to one. A blow up of (f) in a particular k interval emcopassing [Forumla omitted. See PDF.] is shown in (g).

Figure 9. (a) Schematics of a cell formed by two Pöschl–Teller barriers, where the distance between them as a function of k is [Forumla omitted. See PDF.]. The transmission probability of a single cell, [Forumla omitted. See PDF.], and of a lattice with [Forumla omitted. See PDF.] cells, [Forumla omitted. See PDF.], are shown for the cases of (b) [Forumla omitted. See PDF.], (c) [Forumla omitted. See PDF.], (d) [Forumla omitted. See PDF.], [Forumla omitted. See PDF.], (e) [Forumla omitted. See PDF.], [Forumla omitted. See PDF.]. In all cases, [Forumla omitted. See PDF.] and [Forumla omitted. See PDF.]. The transmission probabilities for the isolated left and right barriers are also displayed.

Table 1

Parameters employed to construct the CBBs, Figure 2a, that approximate the Gaussian barrier in Equation (5), with $\mu =0$, $\sigma =\sqrt{2}/2$, $U_0=1$. The CBB-II,III,IV have a central rectangular barrier, with $w=0.4$ and $\tilde{U}=1$. The triangular (t) and trapezoidal ($\tilde{\mathrm{t}}$) shapes to the left and right sides of the CBB centers are specular images of one another. So, only the values of the left Basic Shapes are indicated.

CBB	Basic Shape	w	${\mathit{U}}_{\mathit{a}}$	${\mathit{U}}_{\mathit{b}}$
I	t	$1.5$	0	$1.15$
II	t	$1.3$	0	1
III	t	$0.9$	0	$0.15$
	$\tilde{\mathrm{t}}$	$1.1$	$0.15$	1
IV	t	$0.6$	0	$0.07$
	$\tilde{\mathrm{t}}$	$0.5$	$0.07$	$0.3$
	$\tilde{\mathrm{t}}$	$0.9$	$0.3$	1

Table 2

Parameters of the CBB-IVs in Figure 5a, displayed in the same way as in Table 1. For the rectangular Basic Shape $\tilde{U}=U_0$ and $w=0.05,0.1,0.2,0.3$, respectively, for $\sigma =0.1,0.2,0.4,0.6$.

Gaussian (${\mathit{U}}_{\mathbf{0}}\mathbf{,}\mathit{\sigma }$)	Basic Shape	w	${\mathit{U}}_{\mathit{a}}$	${\mathit{U}}_{\mathit{b}}$
$(4,0.1)$	t	$0.08$	0	$0.3$
	$\tilde{\mathrm{t}}$	$0.07$	$0.3$	$1.3$
	$\tilde{\mathrm{t}}$	$0.125$	$1.3$	4
$(2,0.2)$	t	$0.16$	0	$0.15$
	$\tilde{\mathrm{t}}$	$0.14$	$0.15$	$0.65$
	$\tilde{\mathrm{t}}$	$0.25$	$0.65$	2
$(1,0.4)$	t	$0.35$	0	$0.09$
	$\tilde{\mathrm{t}}$	$0.25$	$0.09$	$0.33$
	$\tilde{\mathrm{t}}$	$0.5$	$0.33$	1
$(0.66,0.6)$	t	$0.5$	0	$0.05$
	$\tilde{\mathrm{t}}$	$0.3$	$0.05$	$0.17$
	$\tilde{\mathrm{t}}$	$0.85$	$0.17$	$0.66$

Appendix A. General R and T Coefficients

Assume a compact support potential, U, which is non-zero only in an interval $a<x<b$ and whose usual quantum reflection and transmission coefficients are denoted by r and t. Considering a plane wave of wavenumber k, incoming from $x<a$ (+) or $x>b$ (−), we can define R and T as $r^{(\pm )}\left(k\right)=exp[\pm 2ik{c}^{(\mp )}]R^{(\pm )}\left(k\right)$ and $t^{(\pm )}\left(k\right)=exp[-ik(b-a)]T^{(\pm )}\left(k\right)$, with $c^{(+)}=b$ and $c^{(-)}=a$. In other words, the scattering amplitudes for a compact support potential can be written as a product of a complex function (R and T) and a phase related to the onset points of U [25]. Note that ${\left|r\right|}^2={\left|R\right|}^2$ and ${\left|t\right|}^2={\left|T\right|}^2$. In addition, for problems with time-reversal symmetry, $T^{(\pm )}=T$. Also, if U is symmetric, $R^{(\pm )}=R$.

When addressing a multiscattering problem, the exact Green’s function is expressed as a sum over all possible scattering paths, where each term is a product of R or T by an exponential involving the classic action. Since the paths already take into account the phases $exp[\pm 2ik{c}^{(\mp )}]$ and $exp[-ik(b-a\left)\right]$, one should use R and T, instead of r and t. This explains the usage of $R^{(\pm )}$ and T in Section 2.

Appendix B. Rectangular, Trapezoidal and Triangular Barriers

Here we summarize the quantum amplitudes for rectangular, trapezoidal and triangular barriers, which are used to construct the CBBs discussed in the present work.

First, consider the symmetric rectangular barrier of width w and height $\tilde{U}$. Its scattering amplitudes are easily found in textbooks. They are given by (A1)$$R_{rect}=-{\displaystyle \frac{i({\kappa }^2+k^2)sinh\left[\kappa w\right]}{\mathcal{A}\left(k\right)}},T_{rect}={\displaystyle \frac{2\kappa k}{\mathcal{A}\left(k\right)}},$$ where $\kappa =\sqrt{2\tilde{U}-k^2}$ and $\mathcal{A}\left(k\right)=2\kappa kcosh\left[\kappa w\right]+i({\kappa }^2-k^2)sinh\left[\kappa w\right]$.

Suppose now the trapezoidal potential (A2)$$\begin{array}{ccc}\hfill U\left(x\right)& =& {\displaystyle \frac{(U_b-U_a)x+(U_{a}b-U_{b}a)}{(b-a)}},a\le x\le b,\hfill \end{array}$$ and $U\left(x\right)=0$ otherwise. It is asymmetric since $U\left(a\right)=U_a$ and $U\left(b\right)=U_b$. From Equation (A2), we can also obtain a triangular barrier by setting either $U_a$ or $U_b$ to zero. Its reflection and transmission coefficients are analytically calculated (see, for example, [27]), so that (A3)$$R_{trap}^{(\pm )}\left(k\right)={\displaystyle \frac{{\mathcal{B}}^{(\pm )}\left(k\right)}{\mathcal{C}\left(k\right)}},T_{trap}\left(k\right)=-{\displaystyle \frac{2\eta \left(k\right)}{\pi \mathcal{C}\left(k\right)}},$$ with (A4)$$\alpha =2{\displaystyle \frac{(U_a-U_b)}{(b-a)}},\beta =2{\displaystyle \frac{(U_{a}b-U_{b}a)}{(b-a)}},y_x=\left[-x+{\displaystyle \frac{(\beta -k^2)}{\alpha }}\right]{\alpha }^{1/3},\eta \left(k\right)={\displaystyle \frac{i}{k}}{\alpha }^{{\displaystyle \frac{1}{3}}},$$ and (A5)$$\begin{array}{ccc}\hfill {\mathcal{B}}^{(\pm )}\left(k\right)& =& [\mathrm{Ai}\left(y_a\right)\mp \eta \left(k\right){\mathrm{Ai}}^{\prime }\left(y_a\right)][\mathrm{Bi}\left(y_b\right)\mp \eta \left(k\right){\mathrm{Bi}}^{\prime }\left(y_b\right)]\hfill \\ & & -[\mathrm{Ai}\left(y_b\right)\mp \eta \left(k\right){\mathrm{Ai}}^{\prime }\left(y_b\right)][\mathrm{Bi}\left(y_a\right)\mp \eta \left(k\right){\mathrm{Bi}}^{\prime }\left(y_a\right)],\hfill \\ \hfill \mathcal{C}\left(k\right)& =& [\mathrm{Ai}\left(y_a\right)+\eta \left(k\right){\mathrm{Ai}}^{\prime }\left(y_a\right)][\mathrm{Bi}\left(y_b\right)-\eta \left(k\right){\mathrm{Bi}}^{\prime }\left(y_b\right)]\hfill \\ & & -[\mathrm{Ai}\left(y_b\right)-\eta \left(k\right){\mathrm{Ai}}^{\prime }\left(y_b\right)][\mathrm{Bi}\left(y_a\right)+\eta \left(k\right){\mathrm{Bi}}^{\prime }\left(y_a\right)].\hfill \end{array}$$ Above, $\mathrm{Ai}\left(y\right)$, ${\mathrm{Ai}}^{\prime }\left(\mathrm{y}\right)$ ($\mathrm{Bi}\left(y\right)$, ${\mathrm{Bi}}^{\prime }\left(\mathrm{y}\right)$) represent, respectively, the Airy function of the first (second) kind and its derivative with respect to y.

Appendix C. Scattering Coefficients for the U_PT(x) = U₀/cosh²[α x] Potential

For the $U_{PT}$ potential of Equation (6), in Equation (4) one uses $R^{(\pm )}=R_{PT}exp\left[ik\varphi \left(K\right)\right]$ and $T=T_{PT}exp\left[ik\varphi \left(K\right)\right]$ for (see [37] and the references therein) (A6)$$\begin{array}{ccc}\hfill T_{PT}& =& {\displaystyle \frac{k}{i\alpha }}{\displaystyle \frac{\mathsf{\Gamma }(-ik/\alpha -s)\mathsf{\Gamma }(-ik/\alpha +s+1)}{\mathsf{\Gamma }(-ik/\alpha +1)\mathsf{\Gamma }(-ik/\alpha +1)}},\hfill \\ \hfill R_{PT}& =& {\displaystyle \frac{\mathsf{\Gamma }(-ik/\alpha -s)\mathsf{\Gamma }(-ik/\alpha +s+1)\mathsf{\Gamma }(ik/\alpha )}{\mathsf{\Gamma }(-s)\mathsf{\Gamma }(s+1)\mathsf{\Gamma }(-ik/\alpha )}}.\hfill \end{array}$$

Here, $\mathsf{\Gamma }(·)$ is the gamma function, $ϵ=U_0/E=2U_0/k^2$ and $s=(-1+\sqrt{1-8U_0/{\alpha }^2})/2$. A phase $\varphi \left(k\right)$ multiplying the reflection and transmission amplitudes of the basic cells is necessary [37] when the potential is not of compact support. For the Pöschl–Teller potential, it reads (A7)$$\varphi \left(k\right)={\displaystyle \frac{1}{\alpha }}\left(2\sqrt{ϵ}ln[\sqrt{ϵ}+1]+(1-\sqrt{ϵ})ln\left[\right|ϵ-1\left|\right]\right).$$

Appendix D. Double Rectangular Barrier Cells and Anomalous Behavior

Consider the building block formed by two rectangular barriers, as illustrated in Figure A1. If the widths $w_1$ and $w_2$ and heights ${\tilde{U}}_1$ and ${\tilde{U}}_2$ are the same (different), we have an SBB (ABB). The reflection and transmission coefficients for a single rectangular potential are given in Appendix B. Hence, the coefficients for the double barrier cell are simply determined from the recurrence relations in Equation (1).

Assume $w_1=w_2=1$, ${\tilde{U}}_1={\tilde{U}}_2=2$ and an unitary separation between any two successive cells for a lattice with $N={10}^4$. For this case, in Figure A1b we show the transmission probabilities as a function of k for $N=1$ and $N={10}^4$. Note that, for this symmetric cell, the $k_{res}$s do not influence the lattice quasi-band structure. On the other hand, for an ABB cell with the exact same parameters, unless for ${\tilde{U}}_2=2.2$, the profiles of $|T_1{\left(k\right)|}^2$ and $|T_{{10}^4}{\left(k\right)|}^2$ are depicted in Figure A1c. Now, we clearly observe rip-like defects (around the cell $k_{res}$s) in the lattice-allowed quasi-bands. These corroborate the results in Section 5

Figure A1. (a) Illustration of a building block formed by two rectangular barriers, with widths [Forumla omitted. See PDF.] and heights [Forumla omitted. See PDF.]. Graph of the transmission probability in terms of the wavenumber k, when we consider (b) a symmetrical and (c) an asymmetrical cell. In (b), [Forumla omitted. See PDF.], and in (c), [Forumla omitted. See PDF.]. The remain parameters are [Forumla omitted. See PDF.], [Forumla omitted. See PDF.] in both figures. The width between any two consecutive barriers is unitary.

References

1. Rosen, J. Symmetry Rules: How Science and Nature Are Founded on Symmetry; 1st ed. Springer: Berlin/Heidelberg, Germany, 2008.

2. Bensoussan, A.; Lions, J.-L.; Papanicolaou, G. Asymptotic Analysis for Periodic Structures; 2nd ed. American Mathematical Society: Providence, RI, USA, 2011.

3. Guest, S.D.; Fowler, P.W.; Power, S.C. Rigidity of periodic and symmetric structures in nature and engineering. Philos. Trans. R. Soc. A; 2014; 372, 20130358. [DOI: https://dx.doi.org/10.1098/rsta.2013.0358] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/24379434]

4. Nefyodov, E.I.; Smolskiy, S.M. Periodic processes and structures in nature, science, and engineering. Electromagnetic Fields and Waves. Textbooks in Telecommunication Engineering; Springer: Cham, Switzerland, 2019.

5. Bran, A.M.; Stadler, P.F.; Jost, J.; Restrepo, G. The six stages of the convergence of the periodic system to its final structure. Commun. Chem.; 2023; 6, 87. [DOI: https://dx.doi.org/10.1038/s42004-023-00883-9] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/37130929]

6. Leitherer, A.; Ziletti, A.; Ghiringhelli, L.M. Robust recognition and exploratory analysis of crystal structures via Bayesian deep learning. Nat. Commun.; 2021; 12, 6234. [DOI: https://dx.doi.org/10.1038/s41467-021-26511-5] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34716341]

7. Barkat, O. Study and Simulation of the Characteristics of Periodic Structures; Our Knowledge Publishing: Delhi, India, 2021.

8. Johnston, I.G.; Dingle, K.; Greenbury, S.F.; Camargo, C.Q.; Doye, J.P.; Ahnert, S.E.; Louis, A.A. Symmetry and simplicity spontaneously emerge from the algorithmic nature of evolution. Proc. Natl. Acad. Sci. USA; 2022; 119, e2113883119. [DOI: https://dx.doi.org/10.1073/pnas.2113883119] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35275794]

9. Newton, R.G. Scattering Theory of Waves and Particles; 2nd ed. Springer: Berlin/Heidelberg, Germany, 2013.

10. Ishimaru, A. Electromagnetic Wave Propagation, Radiation, and Scattering: From Fundamentals to Applications; 1st ed. John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2017.

11. Touhei, T. Theory of Elastic Wave Propagation and Its Application to Scattering Problems; 1st ed. CRC Press: Boca Raton, FL, USA, 2024.

12. Jia, Z.; Liu, F.; Jiang, X.; Wang, L. Engineering lattice metamaterials for extreme property, programmability, and multifunctionality. J. Appl. Phys.; 2020; 127, 150901. [DOI: https://dx.doi.org/10.1063/5.0004724]

13. Schnitzer, O.; Craster, R.V. Bloch waves in an arbitrary two-dimensional lattice of subwavelength Dirichlet scatterers. SIAM J. Appl. Math.; 2017; 77, pp. 2119-2135. [DOI: https://dx.doi.org/10.1137/16M107222X]

14. Ojeda-Aciego, M.; Outrata, J. Concept lattices and their applications. Int. J. Gen. Syst.; 2015; 45, pp. 55-56. [DOI: https://dx.doi.org/10.1080/03081079.2015.1072921]

15. da Luz, M.G.; Anteneodo, C. Nonlinear dynamics in meso and nano scales: Fundamental aspects and applications. Philos. Trans. R. Soc. A; 2011; 369, pp. 245-259. [DOI: https://dx.doi.org/10.1098/rsta.2010.0301]

16. Citrin, D.S. Quasitransparent states in the logarithmic chain and nontrivial zeros of the Riemann zeta function. Phys. Rev. B; 2024; 110, L081406. [DOI: https://dx.doi.org/10.1103/PhysRevB.110.L081406]

17. Ptitsyna, N.; Shipman, S.P. Guided modes and anomalous scattering by a periodic lattice. Proceedings of the 2008 12th International Conference on Mathematical Methods in Electromagnetic Theory; Odessa, Ukraine, 29 June–2 July 2008.

18. Buonocore, S.; Sen, M.; Semperlotti, F. Occurrence of anomalous diffusion and non-local response in highly-scattering acoustic periodic media. New J. Phys.; 2019; 21, 033011. [DOI: https://dx.doi.org/10.1088/1367-2630/aafb7d]

19. Cutolo, A.; Palumbo, S.; Carotenuto, A.R.; Sacco, E.; Fraldi, M. A class of periodic lattices for tuning elastic instabilities. Extrem. Mech. Lett.; 2022; 55, 101839. [DOI: https://dx.doi.org/10.1016/j.eml.2022.101839]

20. Krasnok, A.; Baranov, D.; Li, H.; Miri, M.A.; Monticone, F.; Alú, A. Anomalies in light scattering. Adv. Opt. Photonics; 2019; 11, pp. 892-951. [DOI: https://dx.doi.org/10.1364/AOP.11.000892]

21. Martin, P.A. Multiple Scattering: Interaction of Time-Harmonic Waves with N Obstacles; 1st ed. Cambridge University Press: Cambridge, UK, 2006.

22. Harrington, D.R. Multiple scattering, the Glauber approximation, and the off-shell eikonal approximation. Phys. Rev.; 1969; 184, 1745. [DOI: https://dx.doi.org/10.1103/PhysRev.184.1745]

23. Stamnes, K.; Li, W.; Stamnes, S.; Hu, Y.; Zhou, Y.; Chen, N.; Fan, Y.; Hamre, B.; Lu, X.; Huang, Y. et al. Laser light propagation in a turbid medium: Solution including multiple scattering effects. Eur. J. Phys. D; 2023; 77, 110. [DOI: https://dx.doi.org/10.1140/epjd/s10053-023-00694-6]

24. Glauber, R.J.; Osland, P. Asymptotic Diffraction Theory and Nuclear Scattering; 1st ed. Cambridge University Press: Cambridge, UK, 2019.

25. da Luz, M.G.E.; Heller, E.J.; Cheng, B.K. Exact form of Green’s functions for segmented potentials. J. Phys. A; 1998; 31, 2975. [DOI: https://dx.doi.org/10.1088/0305-4470/31/13/007]

26. Andrade, F.M.; Schmidt, A.G.M.; Vicentini, E.; Cheng, B.K.; da Luz, M.G.E. Green’s function approach for quantum graphs: An overview. Phys. Rep.; 2016; 647, pp. 1-46. [DOI: https://dx.doi.org/10.1016/j.physrep.2016.07.001]

27. Oliveira, L.R.N.; da Luz, M.G.E. Quantum scattering in one-dimensional periodic structures: A Green’s function approach solved through continued fractions. Phys. Rev. B; 2024; 110, 054303. [DOI: https://dx.doi.org/10.1103/PhysRevB.110.054303]

28. Kalotas, T.M.; Lee, A.R. A new approach to one-dimensional scattering. Am. J. Phys.; 1991; 59, 48. [DOI: https://dx.doi.org/10.1119/1.16705]

29. Abdolkader, T.M.; Shaker, A.; Alahmadi, A.N.M. Numerical simulation of tunneling through arbitrary potential barriers applied on MIM and MIIM rectenna diodes. Eur. J. Phys.; 2018; 39, 045402. [DOI: https://dx.doi.org/10.1088/1361-6404/aab5cf]

30. Fernández, F.M. Quantum Gaussian wells and barriers. Am. J. Phys.; 2011; 79, 752. [DOI: https://dx.doi.org/10.1119/1.3574505]

31. Fernández, F.M. Wronskian method for one-dimensional quantum scattering. Am. J. Phys.; 2011; 79, 877. [DOI: https://dx.doi.org/10.1119/1.3596393]

32. Cevik, D.; Gadella, M.; Kuru, Ş.E.N.G.Ü.L.; Negro, J. Resonances and antibound states for the Pöschl–Teller potential: Ladder operators and SUSY partners. Phys. Lett. A; 2016; 380, pp. 1600-1609. [DOI: https://dx.doi.org/10.1016/j.physleta.2016.03.003]

33. Shizgal, B.D. Pseudospectral method of solution of the Schrödinger equation with non classical polynomials: The Morse and Pöschl–Teller (SUSY) potentials. Comput. Theor. Chem.; 2016; 1084, pp. 51-58. [DOI: https://dx.doi.org/10.1016/j.comptc.2016.03.002]

34. Al Sakkaf, L.; Al Khawaja, U. Bound-states spectrum of the nonlinear Schrödinger equation with Pöschl-Teller and square-potential wells. Phys. Rev. E; 2022; 106, 024206. [DOI: https://dx.doi.org/10.1103/PhysRevE.106.024206] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/36109915]

35. Eyube, E.S.; Yusuf, I.; Omugbe, E.; Makasson, C.R.; Onate, C.A.; Mohammed, B.D.; Balami, B.Y.; Tahir, A.M. Energy spectrum and magnetic susceptibility of the improved Pöschl-Teller potential. Phys. B Condens. Matter.; 2024; 416483. [DOI: https://dx.doi.org/10.1016/j.physb.2024.416483]

36. Ahmed, F.; Bouzenada, A.; Moreira, A.R. Effects of Pöschl-Teller potential on approximate ℓ ≠ 0-states solution in topological defect geometry and Shannon entropy. Phys. Scr.; 2024; 99, 075411. [DOI: https://dx.doi.org/10.1088/1402-4896/ad56df]

37. da Luz, M.G.E.; Cheng, B.K.; Beims, M.W. Asymptotic Green functions: A generalized semiclassical approach for scattering by multiple barrier potentials. J. Phys. A; 2001; 34, 5041. [DOI: https://dx.doi.org/10.1088/0305-4470/34/24/303]

38. Manga Rao, V.S.C.; Prasad, M.D.; Gupta, S.D. Self-similarity in nonclassical transmission through a Fibonacci sequence of Gaussian barriers. J. Mod. Phys. C; 2005; 16, 327. [DOI: https://dx.doi.org/10.1142/S0129183105007133]

39. Sánchez-Soto, L.L.; Monzón, J.J.; Barriuso, A.G.; Cariñena, J.F. The transfer matrix: A geometrical perspective. Phys. Rep.; 2012; 513, pp. 191-227. [DOI: https://dx.doi.org/10.1016/j.physrep.2011.10.002]

40. Pereyra, P. The transfer matrix method and the theory of finite periodic systems. From heterostructures to superlattices. Phys. Status Solidi B; 2022; 259, 2100405. [DOI: https://dx.doi.org/10.1002/pssb.202100405]

41. Rumpf, R.C. Improved formulation of scattering matrices for semi-analytical methods that is consistent with convention. Prog. Electromagn. Res. B; 2011; 35, pp. 241-261. [DOI: https://dx.doi.org/10.2528/PIERB11083107]

42. Asbóth, J.K.; Oroszlány, L.; Pályi, A. A Short Course on Topological Insulators; 1st ed. Springer: Cham, Switzerland, 2016.

43. Majidi, M.A. Electron hopping integral renormalization due to anharmonic phonons. J. Phys. Conf. Ser.; 2018; 1011, 012080. [DOI: https://dx.doi.org/10.1088/1742-6596/1011/1/012080]

44. Sprung, D.W.L.; Wu, H.; Martorell, J. Scattering by a finite periodic potential. Am. J. Phys.; 1993; 61, pp. 1118-1123. [DOI: https://dx.doi.org/10.1119/1.17306]

45. Tung, H.H.; Lee, C.P. A novel energy filter using semiconductor superlattices and its application to tunneling time calculations. IEEE J. Quantum Electron.; 1996; 32, 2122. [DOI: https://dx.doi.org/10.1109/3.544758]

46. Gómez, I.; Domínguez-Adame, F.; Diez, E.; Bellani, V. Electron transport across a Gaussian superlattice. J. Appl. Phys.; 1999; 85, 3916. [DOI: https://dx.doi.org/10.1063/1.369764]

47. Diez, E.; Gómez, I.; Domínguez-Adame, F.; Hey, R.; Bellani, V.; Parravicini, G.B. Gaussian semiconductor superlattices. Phys. E; 2000; 7, 832. [DOI: https://dx.doi.org/10.1016/S1386-9477(00)00071-0]

48. Lara, G.A.; Orellana, P. Resonant tunneling through Gaussian superlattices. Rev. Mex. Fis.; 2002; 48, 40.

49. Silba-Vélez, M.; Pérez-Álvarez, R.; Contreras-Solorio, D.A. Transmission and escape in finite superlattices with Gaussian modulation. Rev. Mex. Fis.; 2015; 61, 132.

50. Cabrera, C.I.; Contreras-Solorio, D.A.; Hernandez, L.; Enciso, A.; Rimada, J.C. Gaussian superlattice in GaAs/GaInNAs solar cells. Rev. Mex. Fis.; 2017; 63, 223.

51. Zhang, Y.P.; Yin, Y.H.; Lü, H.H.; Zhang, H.Y. Electronic band gap and transport in graphene superlattice with a Gaussian profile potential voltage. Chin. Phys. B; 2013; 23, 027202. [DOI: https://dx.doi.org/10.1088/1674-1056/23/2/027202]

52. Villegas, D.; de León-Pérez, F.; Pérez-Álvarez, R. Gaussian superlattice for phonons. Microelectron. J.; 2005; 36, 411. [DOI: https://dx.doi.org/10.1016/j.mejo.2005.02.033]

53. Su, W.P.; Schrieffer, J.R.; Heeger, A.J. Solitons in polyacetylene. Phys. Rev. Lett.; 1979; 42, 1698. [DOI: https://dx.doi.org/10.1103/PhysRevLett.42.1698]

54. Su, W.P.; Schrieffer, J.R.; Heeger, A.J. Erratum: Soliton excitations in polyacetylene. Phys. Rev. B; 1983; 28, 1138. [DOI: https://dx.doi.org/10.1103/PhysRevB.28.1138]

55. Batra, N.; Sheet, G. Understanding basic concepts of topological insulators through Su–Schrieffer–Heeger (SSH) model. arXiv; 2019; arXiv: 1906.08435

56. Meier, E.J.; An, F.A.; Gadway, B. Observation of the topological soliton state in the Su–Schrieffer–Heeger model. Nat. Commun.; 2016; 7, 13986. [DOI: https://dx.doi.org/10.1038/ncomms13986]

57. Jovanović, S.; Krasić, M.S. Asymmetric defects in one-dimensional photonic lattices. Laser Phys.; 2021; 31, 023001. [DOI: https://dx.doi.org/10.1088/1555-6611/abd8d6]

58. Liao, W.; Li, N. Energy and momentum diffusion in one-dimensional periodic and asymmetric nonlinear lattices with momentum conservation. Phys. Rev. E; 2019; 99, 062125. [DOI: https://dx.doi.org/10.1103/PhysRevE.99.062125]

59. Ravari, M.K.; Esfahani, S.N.; Andani, M.T.; Kadkhodaei, M.; Ghaei, A.; Karaca, H.; Elahinia, M. On the effects of geometry, defects, and material asymmetry on the mechanical response of shape memory alloy cellular lattice structures. Smart Mater. Struct.; 2016; 25, 025008. [DOI: https://dx.doi.org/10.1088/0964-1726/25/2/025008]

60. Fukuda, M.; Mizubayashi, W.; Kohno, A.; Miyazaki, S.; Hirose, M. Analysis of tunnel current through ultrathin gate oxides. Jpn. J. Appl. Phys.; 1998; 37, L1534. [DOI: https://dx.doi.org/10.1143/JJAP.37.L1534]

61. Städele, M.; Tuttle, B.R.; Hess, K. Tunneling through ultrathin SiO₂ gate oxides from microscopic models. J. Appl. Phys.; 2001; 89, 348. [DOI: https://dx.doi.org/10.1063/1.1330764]

62. Fukuda, M.; Nakagawa, K.; Miyazaki, S.; Hirose, M. Resonant tunneling through a self-assembled Si quantum dot. Appl. Phys. Lett.; 1997; 70, 2291. [DOI: https://dx.doi.org/10.1063/1.118816]

63. Ramírez, C.; González, F.H.; Galván, C.G. Solving Schrödinger equation with scattering matrices. Bound states of Lennard-Jones potential. J. Phys. Soc. Jpn.; 2019; 88, 094002. [DOI: https://dx.doi.org/10.7566/JPSJ.88.094002]

64. Lambropoulos, K.; Simserides, C. Spectral and transmission properties of periodic 1D tight-binding lattices with a generic unit cell: An analysis within the transfer matrix approach. J. Phys. Commun.; 2018; 2, 035013. [DOI: https://dx.doi.org/10.1088/2399-6528/aab065]

65. Kumar, V.; Singh, K.S.; Ojha, S. Band structure, reflection properties and abnormal behaviour of one-dimensional plasma photonic crystals. Prog. Electromagn. Res. M; 2009; 9, pp. 227-241. [DOI: https://dx.doi.org/10.2528/PIERM09101701]

66. Hardhienata, H.; Aziz, A.I.; Rahmawati, D.; Alatas, H. Transmission characteristics of a 1D photonic crystal sandwiched by two graphene layers. J. Phys. Conf. Ser.; 2018; 1057, 012003. [DOI: https://dx.doi.org/10.1088/1742-6596/1057/1/012003]

67. De Vita, F.; Dematteis, G.; Mazzilli, R.; Proment, D.; Lvov, Y.V.; Onorato, M. Anomalous conduction in one-dimensional particle lattices: Wave-turbulence approach. Phys. Rev. E; 2022; 106, 034110. [DOI: https://dx.doi.org/10.1103/PhysRevE.106.034110]