Compact Localized States in Engineered Flat-Band PT\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathscr{P}}{\mathscr{T}}$$\end{document} Metamaterials

The conditions leading to flat dispersionless frequency bands in truly one-dimensional parity-time (PT\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathscr{P}}{\mathscr{T}}$$\end{document}) symmetric metamaterials comprised of split‐ring resonators (SRRs) arranged in a binary pattern are obtained analytically. In this paradigmatic system, in which the SRRs are coupled through both electric and magnetic dipole-dipole forces, flat-bands may arise from tailoring its natural parameters (such as, e.g., the coupling coefficients between SRRs) and not from geometrical effects. For sets of parameters which values are tailored to flatten the upper band of the spectrum, the solution of the corresponding quadratic eigenvalue problem reveals the existence of compact, two-site localized eigenmodes. Numerical simulations confirm the existence and the dynamic stability of such modes, which can be formed through the evolution of single-site initial excitations without disorder or nonlinearity.

where λ E and λ M (λ ′ E and λ ′ M ) are respectively the electric and magnetic coupling coefficients between SRRs with center-to-center distance d (d′), γ is the gain/loss coefficient (γ > 0), and the overdots denote derivation with respect to the normalized temporal variable τ ω = t LC , with ω = − LC LC 1 . By substituting the plane wave solution into Eqs (1 and 2), where κ is the normalized wavevector and Ω is the frequency in units of ω LC , and requesting nontrivial solutions for the resulting stationary problem, we obtain In the exact PT phase, Eq. (4) gives a gapped spectrum with two frequency bands separated by a gap. The FB condition is obtained by requesting κ Ω = κ d d ( )/ 0 2 for any κ in the first Brillouin zone. After tedious calculations, we get www.nature.com/scientificreports www.nature.com/scientificreports/ The values of ′, λ E , λ M , and γ, for which the expression in the first squared brackets in Eq. (6) equals to zero, provide a physically acceptable parameter set which flattens the upper band of the spectrum over the whole Brillouin zone. For such a parameter set, the calculated frequency spectrum (from Eq. (4)) contains a completely flat, isolated upper band ( Fig. 2(a)) with zero group velocity v g for any κ in the first Brillouin zone (inset).
The electric and magnetic coupling coefficients between single SRRs depend crucially on their mutual position. In the model adopted here, the SRRs assume planar geometry; even in this case, however, the coupling coefficients depend strongly on the mutual orientation of the gaps. For the configuration of Fig. 1, these coefficients have been calculated accurately and plotted in Fig. 2(a) of ref. 15 using the approach presented in ref. 50 , where the effect of retardation has been also taken into account. In these works, single SRRs having typical dimensions have been considered, while the authors have verified that the resonance frequencies of pairs of SRRs coupled with the calculated coefficients match the resonances found by direct numerical simulations using commercial software (CST Microwave Studio). This indicates that the calculated coupling coefficients can quantitatively describe the near-field interaction between SRRs. According to ref. 15 , the coupling coefficients used in Fig. 2(a) correspond to distances d∼2 and d′∼3.7 between neighboring SRRs in units of their radii, that seems in principle feasible. On the contrary, for parameter sets not satisfying Eq. (6), non-flat bands with κ− dependent group velocities v g such as those shown in Fig. 2(b) are typically obtained. The band-boundaries (BBs), i.e., the extremal frequencies in each band, are shown as a function of the gain/loss coefficient γ in Fig. 2(c). A critical value of that coefficint, γ = γ c , separates the exact or unbroken from the broken PT phase. For γ < γ c (exact or unbroken PT phase), the BBs of the upper band (red and black dotted curves) practically coincide for a substantial interval of γ indicating zero bandwidth (i.e., a FB), while the width of the lower band (limited by the blue and green dotted curves) remains almost constant up to γ = γ c . The flatness of the upper band for low γ can be seen more clearly in the inset. In Fig. 2(d), the BBs are plotted as a function of ′  ; the width of the upper band (limited by the red and black solid curves) consecutively decreases, goes through zero at a critical value of ′  indicated by the arrow, and then increases with increasing ′  (see also the inset). At that critical value of ′  the two bands coincide with those shown in Fig. 2(a). www.nature.com/scientificreports www.nature.com/scientificreports/ The Quadratic Eigenvalue Problem. Eqs (1 and 2) can be written in matrix form as is an N− dimensional vector, with N being the total number of SRRs, and the N × N matrices M and K are given by 51 with Ω being the eigenvalue and Q e the corresponding N− dimensional eigenvector. It can be solved by standard eigenproblem solvers after its linearization by the classical augmentation procedure 52,53 , which transforms square matrices of order N to 2N. Then, Eq. (10) is reduced to a standard eigenvalue problem (SEP) In what follows, γ < γ c so that the system is in the exact (unbroken) PT phase, and therefore all its eigenvalues are real. A few selected eigenmodes Q e are shown in Fig. 3. The eigenmodes from the non-flat (lower) band are all extended, as that shown in Fig. 3(a). The other three eigenmodes in Fig. 3 belong to the flat (upper) band, and therefore they correspond to the same frequency eigenvalue Ω κ = Ω FB . Figure 3(b-d) show a partially localized, a highly localized, and a compact two-site localized eigenmode, respectively. The insets in Fig. 3(d) enlarge the localization region. The existence of compact localized eigenvectors in a class of FB models has been also demonstrated in ref. 25 . In that work, a model with at least one FB was considered, whose eigenvectors in the Bloch representation may be mixed to obtain highly localized FB eigenvectors, due to macroscopic degeneracy. While there is not any general theorem which states that among all these combinations there will be compact localized eigenvectors, such eigenvectors do exist in several FB models. Compact localized states were constructed, which are actually exact FB localized eigenstates, and then classified according to the number of unit cells they occupy. Such compact localized states occupying one unit cell form a complete and orthogonal basis that makes possible to detangle them from the rest of the lattice. The inversion of the detangling procedure provides the most general FB generator having localized eigenstates that occupy one unit cell. Such models include cross-stich and diamond chains, both quasi-one dimensional at the single-site level.
Here, a purely one-dimensional system is considered in the form of a PT symmetric SRR chain arranged in a binary pattern, whose geometry does not in general provides a FB. However, that system does possess a FB when the coupling coefficients and the gain/loss factor satisfy the condition Eq. (6). In that case, the resulting QEP is directly solved after being transformed into a SEP, demonstrating the existence of a compact localized FB eigenstate as that shown in Fig. 3(d). In the light of these findings, it would be expected that the PTMM in the exact PT phase supports compact localized excitations which in general could be expressed as linear superpositions of a small number of eigenmodes.
The existence of localized states or modes in discrete lattices has been investigated intensively in the past. It has now been established that it may be due to quenched disorder (random lattices), a phenomenon based on wave interference which is known as Anderson localization 54 . In the presence of randomly distributed impurities in a metal, for example, different paths taken by an electron can interfere destructively, leading to localization of its wavefunction. The concept of Anderson localization is applicable to, and has actually been observed experimentally in a variety of physical systems. Localization may also appear in nonlinear lattices which are perfectly periodic; this effect is known as intrinsic localization, which leads to localized states of the discrete breather type 55 . Discrete breather excitations can be created in nonlinear lattices from an initially extended state through the standard modulational instability mechanism. Such localized states have been also experimentally observed in a variety of physical systems. Apart from the case of quenched disorder and nonlinearity, localized states can also arise in tight-binding (nearest-neighbor) models from particular lattice geometries in which destructive interference leads to the emergence of FBs. FB models typically possess localized eigenstates, although there is no general theorem to guarantee their existence. Since in the model considered here, neither disorder nor nonlinearity is www.nature.com/scientificreports www.nature.com/scientificreports/ present, the apearance of compact localized eigenstates for the PTMM is only due to the existence of the FB in its frequency spectrum.
n n n where typically A = 0.1. For N e SRRs at each end of the PTMM, gain has been replaced by equal amount of loss, as if the PT metamaterial were embedded into a lossy metamaterial. It is found empirically that this is the most effective way to stabilize localized modes in PTMMs, even in the presence of nonlinearity 11 (in which case the localized modes are of the discrete breather type). The need for introducing these lossy parts at each end of the PTMM comes from the fact that the initial condition Eq. (12) is clearly not an exact eigenmode of Eqs (1 and 2). However, after long time-integration, that type of initial condition may relax to an exact eigenmode with a lower energy than the initial one. During this process, some of the initial energy spreads towards the ends of the PTMM, where it is dissipated. Thus, these lossy parts help the excess energy to go smoothly away during the transient phase of integration, and allow the formation of a stable compact localized FB state with constant energy. This gradual energy removal is clearly observed in Fig. 4(a) discussed below. The total energy, the (energetic) participation number, and the second moment respectively, where ε = ∑ n n n n is the center of energy with ε = E E / n n tot being the energy density, are shown in Fig. 4 as functions of τ. The four curves are obtained for parameter sets which satisfy the FB condition; the coupling coefficients have very similar values, while γ increases from 0.03 to 0.12 in steps of 0.03. In all cases, a steady FB state is reached at the end of the integration time; however, the transient period is longer and E tot is higher for higher γ. In the insets of Fig. 4(a-c), the quantities E tot , P, and m 2 , respectively, are plotted for the first 25,000 time www.nature.com/scientificreports www.nature.com/scientificreports/ units (t.u.). The energy E tot strongly fluctuates, but it remains on average constant until τ∼12,000 t.u. At the same time, P increases while fluctuating strongly and reaches P = 6, while m 2 increases ∝τ 2 indicating ballistic spreading of the initial state. Around τ∼12,000 t.u., the energy which spreads out from the initial state has reached the lossy ends of the PTMM where it is dissipated. After that time, E tot , P, and m 2 decrease gradually until they saturate to constant values with vanishingly small fluctuations, indicating that a steady state has been reached. The constancy of m 2 , in particular, indicates that energy spreading has been stopped. Part of the initial excitation remains localized on two SRRs, as indicated by the value of P = 2 which measures the number of the energetically strongest excited sites.
Inspection of the q n (τ) in Fig. 5 reveals that the excitation is localized to only two, neighboring SRRs. The charges stored in the capacitors of these SRRs oscillate in anti-phase with frequency Ω FB ; charge oscillations in the rest of the SRRs seems to have vanishingly small amplitude (see also Fig. 6(a)). Note that if the FB condition Eq. (6) is not satisfied, the initial excitations disperse rapidly in the lattice without leaving behind any trace of localization. Compact localized excitations, which tails decay as a stretched exponential or superexponential, often appear in discrete systems with nonlinear dispersion 56 . Here, the compact localized states appear solely due to the FB in the absence of nonlinearity or disorder, and hence there are significant differences. In Fig. 6(b), the quantity is plotted as a function of n; the three curves correspond to different integration times, τ 0 . In obtaining the results in Fig. 6(b), convergence to a steady localized FB state was accelerated by continually eliminating the energy spreading away from the localization region during half of the integration time, i.e., for τ 0 /2 time units. During this time, we set = = (localization region). Note that for the chosen initial condition, the compact two-site FB localized state is generated at the sites = + n N /2 1 and = + n N /2 2, i.e., in the middle of the localization region. For τ τ > /2 0 the integration proceeds without elimination of the energy for τ 0 /2 more time units. The almost horizontal segments between < <   n 16 6 3 and < <   n 66 112 (i.e., those parts of the lattice which do not belong either to the lossy ends nor to the localization region) correpond to the tails of the localized FB states. With increasing τ 0 , the tails become more and more negative until they saturate at ∼ − y 35 n (the limit of double precision arithmetics) at τ 0 = 2 × 10 6 . In the inset, the y n profile at τ 0 = 3 × 10 6 t.u. (black curve) is fitted by  www.nature.com/scientificreports www.nature.com/scientificreports/ b = −35.684 and c = −0.109955. Using the procedure of energy elimination, the relaxation time towards the formation of a compact localized FB state is reduced considerably. For example, after 3 × 10 6 time units of time-integration, the quantities E tot , P, and m 2 take respectively the values 2.855 × 10 −3 , 2.005, and 2.826 × 10 −1 . In order to reach these values without energy elimination, the integration time needed is respectively ∼9.13 × 10 6 , ∼2.82 × 10 7 , and ∼4.62 × 10 7 time units. Note that without energy elimination there are still significant fluctuations around the above values for E tot , P, and m 2 . Thus, the quantities E tot , P, and m 2 do not seem to relax at the same rate when they are calculated with and without elimination of energy. However, as far as the localization region is concerned, a comparison of the relaxation times for .  P 2 005, i.e., 3 × 10 6 and ∼2.82 × 10 7 , indicates a difference of an order of magnitude. Note that .  P 2 005 indicates that there are only two strongly excited sites in the PTMM.

Conclusion
By tailoring the model parameters of a 1D PTMM comprising SRRs arranged in a binary pattern, an isolated and completely FB may appear in its frequency spectrum. The analytical condition, which those parameters have to satisfy in order to flatten the upper band, is obtained. The solution of the QEP reveals the existence of compact, localized FB eigenmodes, in which most of the energy is concentrated in two neighboring SRRs which are separated by distance d′. This is consistent with earlier results on nonlinear binary PTMMs, for which the fundamental discrete breathers are two-site ones (and not single-site ones) 11 . The formation of compact, two-site localized FB states from single-site initial excitations is numerically confirmed. Since the FB is isolated, these FB states could be continued into compact breather-like (nonlinear) excitations. Note that the possibility of flattening one of the bands of the spectrum is solely due to model parameter engineering, and not to any geometrical effects. In this aspect, the existence of two types of coupling between SRRs, i.e., electric and magnetic coupling, is crucial. Note also that similar results would have been obtained for γ = 0, in which case only three parameters, i.e., λ E , λ M , and ′, should be matched to satisfy the corresponding FB condition. However, as it is demonstrated here, the band-flattening capability is not harmed by the PT symmetry, as long as the PTMM is in the exact PT phase.   Fig. 5. (b) The function y n after integration for τ 0 = 1.5 × 10 6 (black), 2 × 10 6 (red), 3 × 10 6 (green), time units. Inset: Fitting of y n with an appropriate function (see text).