Topologically Enhanced Harmonic Generation in a Nonlinear Transmission Line Metamaterial

Topologically Enhanced Harmonic Generation in a Nonlinear Transmission Line Metamaterial

You Wang Division of Physics and Applied Physics, School of Physical and Mathematical Sciences,
Nanyang Technological University, Singapore 637371, Singapore
   Li-Jun Lang Division of Physics and Applied Physics, School of Physical and Mathematical Sciences,
Nanyang Technological University, Singapore 637371, Singapore
   Ching Hua Lee Institute of High Performance Computing, A*STAR, Singapore 138632, Singapore Department of Physics, National University of Singapore, Singapore 117551, Singapore    Baile Zhang Division of Physics and Applied Physics, School of Physical and Mathematical Sciences,
Nanyang Technological University, Singapore 637371, Singapore
Centre for Disruptive Photonic Technologies, Nanyang Technological University, Singapore 637371, Singapore
   Y. D. Chong Division of Physics and Applied Physics, School of Physical and Mathematical Sciences,
Nanyang Technological University, Singapore 637371, Singapore
Centre for Disruptive Photonic Technologies, Nanyang Technological University, Singapore 637371, Singapore
July 26, 2019

Nonlinear transmission lines (NLTLs) are nonlinear electronic circuits commonly used for parametric amplification and pulse generation Cullen (1958); Tien (1958); Landauer (1960a, b). It has previously been shown that harmonic generation can be enhanced, and shock waves suppressed, in so-called “left-handed” NLTLs Lai et al. (2004); Kozyrev and van der Weide (2005, 2008); Powell et al. (2009), a manifestation of the unique properties of left-handed media Vesalago (1968). Here, we demonstrate that harmonic generation in a left-handed NLTL can be greatly increased by the presence of a topological edge state. Our NLTL is a nonlinear analogue of the Su-Schrieffer-Heeger (SSH) lattice Su et al. (1979). Recent studies of nonlinear SSH circuits have investigated the solitonic and self-focusing behaviors of modes at the fundamental harmonic Hadad et al. (2016); Goren et al. (2018); Hadad et al. (2018). We find, however, that frequency-mixing processes in an SSH NLTL have important effects that have previously been neglected. The presence of a topological edge mode at the first harmonic can produce strong higher-harmonic signals that propagate into the lattice, acting as an effectively nonlocal cross-phase nonlinearity. We observe maximum third-harmonic signal intensities that are 5 times that of a comparable left-handed NLTL of a conventional design, and a 250-fold intensity contrast between the topologically nontrivial and trivial lattice configurations. Our work may have applications for compact electronic frequency generators, as well as for advancing the fundamental understanding of the effects of nonlinearities on topological states.

Topological edge states—robust bound states that are guaranteed to exist at the boundary between media with “topologically incompatible” bandstructures—were first discovered in condensed matter physics Bernevig and Hughes (2013). Recently, however, electronic LC circuits have emerged as a highly promising method of realizing these remarkable phenomena Jia et al. (2015); Albert et al. (2015); Lee et al. (); Imhof et al. (); Goren et al. (2018); Zhu et al. (2018); Hadad et al. (2018); Serra-Garcia et al. (). Compared to other classical platforms like photonics Wang et al. (2009); Lu et al. (2014); Khanikaev and Shvets (2017); Ozawa et al. (), acoustics Yang et al. (2015); Fleury et al. (2016); He et al. (2016), and mechanical lattices Huber (2016); Lee et al. (2018), which have also been used to realize topologically nontrivial bandstructures and topological edge states, electronic circuits have several compelling advantages: extreme ease of experimental analysis; the ability to fabricate complicated structures via printed circuit board (PCB) technology; and the intriguing prospect of introducing nonlinear and/or amplifying circuit elements to easily study how topological edge states behave in novel physical regimes. Notably, circuits have been used to study the Su-Schrieffer-Heeger (SSH) chain (the simplest one-dimensional topologically-nontrivial lattice) Lee et al. (), nonlinear SSH chains supporting solitonic edge states Hadad et al. (2018), two-dimensional topological insulator lattices Jia et al. (2015), and the corner states of high-order topological insulators Imhof et al. (); Serra-Garcia et al. ().

One of the most interesting questions raised by the emergence of topologically nontrivial classical lattices is how topological edge states interact with nonlinear media. Previous studies have focused particularly on nonlinearity-induced local self-interactions in the fundamental harmonic, which can give rise to solitons with anomalous plateau-like decay profiles in nonlinear SSH chains Hadad et al. (2016, 2018), or chiral solitons in two-dimensional lattices Lumer et al. (2013); Ablowitz et al. (2014); Leykam and Chong (2016); Lumer et al. (2016); Kartashov and Skryabin (2016); Gulevich et al. (2017). It has also been suggested that topological edge states in nonlinear lattices could be used for robust traveling-wave parametric amplification Peano et al. (2016), optical isolation Zhou et al. (2017), and other applications Rosenthal et al. (2018); Chacón et al. (); Kruk et al. (2017, 2018).

Figure 1: (a) Schematic of a left-handed transmission line circuit with alternating capacitances: linear capacitors with capacitance , and back-to-back varactors with nonlinear capacitance . The capacitances act like hoppings in a nonlinear Su-Schrieffer-Heeger (SSH) model. An input voltage is applied at points or to probe the topologically trivial or nontrivial lattice. (b) Photograph of the printed circuit board. (c) Capacitance versus bias voltage. Dots are calculated from varactor manufacturer data, and the solid curve is the fit based on Eq. (5). (d)–(e) Calculated eigenfrequencies of a finite closed linear circuit, versus the capacitance ratio . The characteristic frequency is , and the lattice has 40 sites. Two cases are shown: (d) -type capacitors at the edge, for which the gap is topologically nontrivial; (e) -type capacitors at the edge, for which the gap is trivial. The red dotted curves indicate the band-edge frequencies and . The blue dashes indicate the operating regime of the nonlinear circuit, with in the linear (low-voltage) limit and effectively increasing with voltage amplitude.

We report on the implementation of a nonlinear SSH chain, based on a left-handed NLTL, in which the topological edge state induces highly efficient harmonic generation. The first harmonic mode is localized to the lattice edge, similar to a linear topological edge state, whereas the higher-harmonic waves propagate into the lattice bulk, with voltage amplitudes reaching over an order of magnitude larger than the first harmonic signal. The intensity of the generated third harmonic signal has a maximum of times that of the input first-harmonic signal, compared to for a comparable conventional left-handed NLTL without a topological edge state. Moreover, the third-harmonic intensity is 250 times larger than in a “trivial” circuit (which has equivalent parameters but lacks a topological edge state in the linear limit) with the same input parameters. Although previous studies have emphasized the role of local self-interactions (including in a previous demonstration of a nonlinear SSH circuit based on weakly-coupled LC resonators Hadad et al. (2018)), the higher-harmonic signals in our circuit play a decisive role in the nonlinear modulation acting on the first-harmonic mode; they have the effect of driving the entire lattice, not just the edge, deeper into the nontrivial regime. This is due to the fact that the left-handed NLTL has an unbounded dispersion curve supporting traveling-wave higher-harmonic modes.

Circuit design—The transmission line circuit, shown schematically in Fig. 1(a), contains inductors of inductance and capacitors of alternating (dimerized) capacitances and . We will shortly treat the case where the capacitors are nonlinear (the and elements are always linear). First, consider the linear limit where is a constant. We define the characteristic angular frequency , and the capacitance ratio


The case of corresponds to a standard (non-dimerized) left-handed transmission line Lai et al. (2004); Kozyrev and van der Weide (2005, 2008); Powell et al. (2009). This type of transmission line is characterized by having sites separated by capactitors, and connected to ground by inductors, rather than vice versa.

Let us treat the points adjacent to the capacitors as lattice sites, indexed by an integer , and close the circuit by grounding the edges [the left edge is the site labelled A in Fig. 1(a)]. Using Kirchhoff’s laws, we can show that a mode with angular frequency satisfies SM ()


where denotes the complex voltage on site . The matrix has the form of the SSH Hamiltonian:


Thus, the eigenfrequency modes of the circuit have a one-to-one correspondence with the SSH eigenstates.

The band diagram for the linear closed circuit is shown in Fig. 1(d). Note the lack of an upper cutoff frequency, a characteristic of left-handed transmission lines Kozyrev and van der Weide (2008). There is a bandgap for . For , the bandgap contains edge states, which are zero-eigenvalue eigenstates of that can be characterized via a topological invariant derived from the Zak phase Bernevig and Hughes (2013). Note, however, that the edge state’s angular frequency is


not zero. Conversely, for , there is a finite bandgap below , which is topologically trivial and contains no edge states. If we swap the two types of capacitors (so that the -type capacitors are the ones at the edge), then the bandgap is trivial and the bandgap nontrivial, as shown in Fig. 1(e).

Next, consider a nonlinear circuit with each capacitor consisting of a pair of back-to-back varactors. The nonlinear capacitance decreases with the magnitude of the bias voltage (the voltage between the end-points of the capacitor), as shown in Fig. 1(c). For theoretical analyses, it is convenient to model this nonlinearity by


where , and is the bias voltage. The key feature of the nonlinearity is that at higher voltages, the effective value of increases; depending on the chosen boundary conditions, this drives the circuit deeper into the topologically trivial or nontrivial regime.

Experimental results—The implemented nonlinear transmission line, shown Fig. 1(b), contains a total of 40 sites, or 20 unit cells. The linear circuit elements have and , so that . By fitting Eq. (5) to manufacturer data for the varactors at low bias voltages SM (), we obtain and (thus, in the linear limit, ). The fitted capacitance-voltage relation is shown in Fig. 1(c).

We supply a continuous-wave sinusoidal input voltage signal, with tunable frequency and amplitude , to either of the points labelled and in Fig. 1(a). This allows us to study the cases corresponding to Fig. 1(d) and (e), which we refer to as the “nontrivial” and “trivial” lattices respectively. In either case, the input site is denoted as .

Figure 2: (a)–(b) Magnitude of the first- and third-harmonic voltage signals measured at different lattice sites, for (a) the nontrivial lattice, which has an SSH-like edge state in the linear limit, and (b) the trivial lattice, which has no edge state in the linear limit. The sinusoidal input signal, applied at the lattice edge (site 0), has frequency and amplitude . (c) Measured spectrum at site 3 for the nontrivial lattice corresponding to (a). (d)–(e) Plot of the inverse participation ratio (IPR) versus input frequency and input voltage amplitude , calculated from experimental measurements of the first-harmonic signal in the (d) nontrivial and (e) trivial lattices. Here, is measured in steps of , and in steps of .

A typical set of measurement results is shown in Fig. 2(a)–(c), for and . On each site , the spectrum of the voltage signal is shown in Fig. 2(c), with prominent peaks at odd harmonics (, , , etc.); even harmonics are suppressed due to the symmetry of the capacitance-voltage relation Kozyrev and van der Weide (2005). Focusing on the first and third harmonics, we define the respective peak values as and , and use these to plot Fig. 2(a)–(b). In the Supplemental Material, we show that results from the SPICE circuit simulator agree well with the experimental data SM ().

From Fig. 2(a)–(b), we see that the nontrivial and trivial lattices exhibit very different behaviors for both the first- and third-harmonic signals. First, consider the first-harmonic signal. In both lattices, there is an exponential decay away from the edge, but the decay is sharper in the nontrivial lattice, which may be attributed to the enhanced intensity arising from the coupling of the input signal to the topological edge state. As a quantitative measure of the localization of the first-harmonic signal, Fig. 2(d)–(e) shows the inverse participation ratio (IPR) ; a larger IPR corresponds to a more localized profile Thouless (1974). We see that the IPR is substantially larger in the nontrivial lattice than in the trivial lattice, over a broad range of and . The strong difference in localization is a key signature of nonlinearity: in the linear regime, a driving voltage on the edges of the nontrivial and trivial lattices would produce different overall amplitudes, but the same exponential decay profile SM (). It is interesting to note that the region of enhanced IPR, shown in Fig. 2(d), closely resembles the nontrivial bandgap in Fig. 1(d).

Figure 3: (a)–(b) Normalized third-harmonic signal intensity versus input frequency and input voltage , for the (a) nontrivial and (b) trivial lattice. Here, is derived from experimental data using the definition (6). (c) Ratio of the trivial and nontrivial intensities, , within the region indicated by boxes in (a) and (b). (d) Measured third-harmonic intensities (normalized to the input signal) at different sites, for the three sets of input parameters indicated by stars in (c): , , and , with fixed .

We can also see from Fig. 2(a) and (c) that strong higher-harmonic signals are present in the nontrivial lattice. Moreover, Fig. 2(a) indicates that the third-harmonic signal is extended, not localized to the edge. To understand this in more detail, we define


which quantifies the intensity of the third-harmonic signal relative to the input intensity at the first harmonic. Here, denoting an average over the first ten lattice sites. Fig. 3(a)–(b) plots the variation of with and . In the nontrivial circuit, the maximum value of the normalized intensity is for and . The fact that peaks over a relatively narrow frequency range, as shown in Fig. 3(a), may be a finite-size effect: the high-frequency modes of the lattice form discrete sub-bands due to the finite lattice size [see Fig. 1(d)–(e)]. In computer simulations, we obtained a similar maximum value of for the nontrivial lattice, whereas a comparable left-handed NLTL of the usual design (containing only identical nonlinear capacitances) has maximum SM ().

The trivial lattice exhibits a much weaker third-harmonic signal. As indicated in Fig. 3(c), for certain choices of and , the value of in the nontrivial lattice is times that in the trivial lattice. Fig. 3(d) plots the normalized third-harmonic signal intensities versus the site index , showing that they do not decay exponentially away from the edge. In the nontrivial lattice, the normalized third-harmonic signal increases with (i.e., stronger nonlinearity).

Discussion—Our results point to a complex interplay between the topological edge state and higher-harmonic modes in the SSH-like NLTL. When a topological edge state exists in the linear lattice, it can be excited by an input signal at frequencies matching the linear lattice’s bandgap. The importance of the edge state is evident from the comparisons, in Fig. 2 and Fig. 3, between the topologically trivial and nontrivial lattices; note also that when the excitation frequency lies outside the linear bandgap, the two lattices behave similarly and the harmonic generation is relatively weak. In the topologically nontrivial lattice, the resonant excitation generates third- and higher-harmonic signals that penetrate deep into the lattice, unlike the first-harmonic mode which is localized to the edge. A few sites away from the edge, the higher-harmonic signals become stronger than the first harmonic, and hence dominate the effective value of the nonlinear parameter. In the linear lattice, is the parameter that “drives” the topological transition, and increasing leads to a larger bandgap and hence a more confined edge state. In the nonlinear regime, Fig. 2 shows an order-of-magnitude increase in the third-harmonic signal amplitude in the nontrivial lattice, relative to the trivial lattice; this implies an effective increase in , and indeed we see that the first-harmonic mode profile is more strongly localized. A more localized edge state, in turn, produces a stronger response to an input signal SM ().

The above interpretation is supported by a more detailed analysis of the coupled equations governing the different circuit mode harmonics, given in the Supplemental Material SM (). These equations involve an effective parameter whose approximate value, in the -th unit cell, is , where is the -th harmonic of the bias voltage on the nonlinear capacitor in the -th unit cell, and We further show that propagating waves can be self-consistently realized for higher () harmonics in the presence of non-linearity, even if the fundamental () mode only has decaying solutions. The first-harmonic mode is localized to the edge, with localization length decreasing with in a manner similar to the linear (SSH-like) lattice.

The generation of the higher-harmonic signals occurs mainly near the edge of the lattice, where the first-harmonic mode is largest. The nonlinearity-induced harmonic generation is aided by the well-known fact that the SSH edge state changes sign in each unit cell (corresponding to the fact that the gap closing in the SSH model takes place at the corner of the Brillouin zone). This feature increases the bias voltages across the nonlinear capacitors, which can thus exceed the values of the voltages at individual sites.

The role of higher-harmonic signals distinguishes our system from previous studies of nonlinear topological edge states, which were based on nonlinear self-modulation at a single harmonic. For instance, in a nonlinear SSH lattice where the coupling depends on the local intensity of a single mode, soliton-like edge states with anomalous mode profiles were predicted Hadad et al. (2016), and subsequently verified using a NLTL-like circuit with narrow frequency bands Hadad et al. (2018) (which, unlike our present circuit, did not support propagating higher-harmonic modes). Topological solitons based on nonlinear self-modulation are also predicted to exist in higher-dimensional lattices Lumer et al. (2013); Ablowitz et al. (2014); Leykam and Chong (2016); Lumer et al. (2016); Kartashov and Skryabin (2016); Gulevich et al. (2017). In our case, the effective value of away from the edge is dominated by the higher-harmonic signals; from the point of view of the first-harmonic mode, these act as a nonlocal nonlinearity, driving the entire lattice deeper into the topologically nontrivial regime, not just the sites with large first-harmonic intensity.

Our work opens the door to the application of topological edge states for enhancing harmonic generation, not just in transmission line circuits, but also a variety of other interesting systems. These include two-dimensional electronic lattices, where topological edge states have already been observed in the linear regime Jia et al. (2015), and the unidirectional nature of the edge states may be even more beneficial for frequency-mixing Peano et al. (2016). Higher dimensional circuit lattices can also possess different thresholds for bulk propagation in different directions, with an extreme generalization being that of a corner mode circuit constructed in Ref. Lee et al. ().

Methods—The nonlinear transmission line was implemented on a PCB (Seeed Tech. Co), with each nonlinear capacitor consisting of a pair of back-to-back varactors (Skyworks Solutions, SMV1253-004LF). A function generator (Tektronix AFG3102C) supplies the continuous-wave sinusoidal input voltage, and the voltages on successive lattice sites, , are measured by an oscilloscope (Rohde & Schwarz RTE1024) in high-impedance mode. Numerical results were obtained using the SPICE circuit simulator.

Acknowledgments—We are grateful to H. Wang and D. Leykam for helpful discussions. YW, LJL, BZ, and CYD were supported by the Singapore MOE Academic Research Fund Tier 2 Grant MOE2015-T2-2-008, and the Singapore MOE Academic Research Fund Tier 3 Grant MOE2016-T3-1-006.


  • Cullen (1958) A. L. Cullen, “A travelling-wave parametric amplifier,” Nature 181, 332–332 (1958).
  • Tien (1958) P. K. Tien, “Parametric amplification and frequency mixing in propagating circuits,” Journal of Applied Physics 29, 1347–1357 (1958).
  • Landauer (1960a) R. Landauer, “Shock waves in nonlinear transmission lines and their effect on parametric amplification,” IBM Journal of Research and Development 4, 391–401 (1960a).
  • Landauer (1960b) Rolf Landauer, “Parametric amplification along nonlinear transmission lines,” Journal of Applied Physics 31, 479–484 (1960b).
  • Lai et al. (2004) A. Lai, T. Itoh,  and C. Caloz, “Composite right/left-handed transmission line metamaterials,” IEEE Microw. Mag. 5 (2004).
  • Kozyrev and van der Weide (2005) A.B. Kozyrev and D.W. van der Weide, “Nonlinear wave propagation phenomena in left-handed transmission-line media,” IEEE Transactions on Microwave Theory and Techniques 53, 238–245 (2005).
  • Kozyrev and van der Weide (2008) A. B. Kozyrev and D. W. van der Weide, “Nonlinear left-handed transmission line metamaterials,” J. Phys. D: Appl. Phys. 41, 173001 (2008).
  • Powell et al. (2009) D. A. Powell, I. V. Shadrivov,  and Y. S. Kivshar, ‘‘Asymmetric parametric amplification in nonlinear left-handed transmission lines,” Appl. Phys. Lett. 94, 084105 (2009).
  • Vesalago (1968) V. G. Vesalago, “The electrodynamics of substances with simultaneously negative values of and ,” Sov. Phys. Usp. , 509 (1968).
  • Su et al. (1979) W. P. Su, J. R. Schrieffer,  and A. J. Heeger, “Solitons in polyacetylene,” Phys. Rev. Lett. 42 (1979).
  • Hadad et al. (2016) Y. Hadad, A. B. Khanikaev,  and A. Alù, “Self-induced topological transitions and edge states supported by nonlinear staggered potentials,” Phys. Rev. B 93, 155112 (2016).
  • Goren et al. (2018) T. Goren, K. Plekhanov, F. Appas,  and K. Le Hur, “Topological zak phase in strongly coupled lc circuits,” Phys. Rev. B 97, 041106(R) (2018).
  • Hadad et al. (2018) Y. Hadad, J. C. Soric, A. B. Khanikaev,  and A. Alù, “Self-induced topological protection in nonlinear circuit arrays,” Nat. Elect. 1, 178 (2018).
  • Bernevig and Hughes (2013) B.A. Bernevig and T.L. Hughes, Topological Insulators and Topological Superconductors (Princeton, 2013).
  • Jia et al. (2015) N. Jia, C. Owens, A. Sommer, D. Schuster,  and J. Simon, “Time- and site-resolved dynamics in a topological circuit,” Phys. Rev. X 5, 021031 (2015).
  • Albert et al. (2015) V. V. Albert, L. I. Glazman,  and L. Jiang, “Topological properties of linear circuit lattices,” Phys. Rev. Lett. 114, 173902 (2015).
  • (17) C. H. Lee, S. Imhof, C. Berger, F. Bayer, J. Brehm, L. W. Molenkamp, T. Kiessling,  and R. Thomale, ‘‘Topolectrical circuits,” arXiv:1705.01077 .
  • (18) S. Imhof, C. Berger, F. Bayer, J. Brehm, L. Molenkamp, T. Kiessling, F. Schindler, C. H. Lee, M. Greiter, T. Neupert,  and R. Thomale, “Topolectrical circuit realization of topological corner modes,” arXiv:1708.03647 .
  • Zhu et al. (2018) W. Zhu, S. Hou, Y. Long, H. Chen,  and J. Ren, “Simulating quantum spin hall effect in the topological lieb lattice of a linear circuit network,” Phys. Rev. B 97, 075310 (2018).
  • (20) M. Serra-Garcia, R. Süsstrunk,  and S. D. Huber, “Observation of quadrupole transitions and edge mode topology in an LC network,” arXiv:1806.07367 .
  • Wang et al. (2009) Z. Wang, Y. D. Chong, J. D. Joannopoulos,  and M. Soljačić, “Observation of unidirectional backscattering-immune topological electromagnetic states,” Nature 461, 772 (2009).
  • Lu et al. (2014) L. Lu, J. D. Joannopoulos,  and M. Soljačić, “Topological photonics,” Nat. Photon. 8, 821 (2014).
  • Khanikaev and Shvets (2017) A. B. Khanikaev and G. Shvets, “Two-dimensional topological photonics,” Nat. Photon. 11, 763 (2017).
  • (24) T. Ozawa, H. M. Price, A. Amo, N. Goldman, M. Hafezi, L. Lu, M. C. Rechtsman, D. Schuster, J. Simon, O. Zilberberg,  and I. Carusotto, “Topological photonics,” arXiv:1802.04173 .
  • Yang et al. (2015) Z. Yang, F. Gao, X. Shi, X. Lin, Z. Gao, Y. D. Chong,  and B. Zhang, “Topological acoustics,” Phys. Rev. Lett. 114, 114301 (2015).
  • Fleury et al. (2016) R. Fleury, A. B. Khanikaev,  and A. Alù, “Floquet topological insulators for sound,” Nat. Comm. 7, 11744 (2016).
  • He et al. (2016) C. He, X. Ni, H. Ge, X.-C. Sun, Y.-B. Chen, M.-H. Lu, X.-P. Liu,  and Y.-F. Chen, “Acoustic topological insulator and robust one-way sound transport,” Nat. Phys. 12, 1124 (2016).
  • Huber (2016) S. D. Huber, “Topological mechanics,” Nat. Phys. 12, 621 (2016).
  • Lee et al. (2018) C. H. Lee, G. Li, G. Jin, Y. Liu,  and X. Zhang, “Topological dynamics of gyroscopic and floquet lattices from newton’s laws,” Phys. Rev. B 97, 085110 (2018).
  • Lumer et al. (2013) Y. Lumer, Y. Plotnik, M. C. Rechtsman,  and M. Segev, “Self-localized states in photonic topological insulators,” Phys. Rev. Lett. 111, 243905 (2013).
  • Ablowitz et al. (2014) M. J. Ablowitz, C. W. Curtis,  and Y.-P. Ma, “Linear and nonlinear traveling edge waves in optical honeycomb lattices,” Phys. Rev. A 90, 023813 (2014).
  • Leykam and Chong (2016) D. Leykam and Y. D. Chong, “Edge solitons in nonlinear-photonic topological insulators,” Phys. Rev. Lett. 117, 143901 (2016).
  • Lumer et al. (2016) Y. Lumer, M. C. Rechtsman, Y. Plotnik,  and Segev M., “Instability of bosonic topological edge states in the presence of interactions,” Phys. Rev. A 94, 021801(R) (2016).
  • Kartashov and Skryabin (2016) Y. V. Kartashov and D. V. Skryabin, “Modulational instability and solitary waves in polariton topological insulators,” Optica 3, 1228 (2016).
  • Gulevich et al. (2017) D. R. Gulevich, D. Yudin, D. V. Skryabin, I. V. Iorsh,  and I. A. Shelykh, ‘‘Exploring nonlinear topological states of matter with exciton-polaritons: Edge solitons in kagome lattice,” Sci. Rep. 7, 1780 (2017).
  • Peano et al. (2016) V. Peano, M. Houde, F. Marquardt,  and A. A. Clerk, “Topological quantum fluctuations and traveling wave amplifiers,” Phys. Rev. X 6, 041026 (2016).
  • Zhou et al. (2017) X. Zhou, Y. Wang, D. Leykam,  and Y. D. Chong, ‘‘Optical isolation with nonlinear topological photonics,” New J. Phys. 19, 095002 (2017).
  • Rosenthal et al. (2018) E. I. Rosenthal, N. K. Ehrlich, M. S. Rudner, A. P. Higginbotham,  and K. W. Lehnert, “Topological phase transition measured in a dissipative metamaterial,” Phys. Rev. B 97, 220301(R) (2018).
  • (39) A. Chacón, W. Zhu, S. P. Kelly, A. Dauphin, E. Pisanty, A. Picón, C. Ticknor, M. F. Ciappina, A. Saxena,  and M. Lewenstein, “Observing topological phase transitions with high harmonic generation,” arXiv:1807.01616 .
  • Kruk et al. (2017) S. Kruk, A. Slobozhanyuk1, A. Shorokhov, I. Kravchenko, A. Poddubny, A. Poddubny,  and Y. Kivshar, ‘‘Third-harmonic generation from photonic topological states in zigzag arrays of silicon nanodisks,” Proc. Fontiers in Optics , FTh4B.2 (2017).
  • Kruk et al. (2018) S. Kruk, L. Wang, A. Slobozhanyuk1, A. Shorokhov, D. Smirnova, I. Kravchenko, A. Poddubny,  and Y. Kivshar, “Topological phase transition measured in a dissipative metamaterial,” Proc. CLEO/QELS , FM2E.4 (2018).
  • (42) See Supplemental Material.
  • Thouless (1974) D. J. Thouless, “Electrons in disordered systems and the theory of localization,” Physics Reports 13, 93 (1974).

Supplemental Material for:

Topologically Enhanced Harmonic Generation in a Nonlinear Transmission Line Metamaterial

Appendix S1 Linear circuit

Consider the circuit shown in Fig. S1 below. This is the same as Fig. 1(a) of the main text, but with the complex voltage variables relabeled for convenience. On unit cell , the voltages on the two sites are (to the right of the capacitor) and (to the right of the capacitor). Also, we let () denote the current through the inductor to the right of the () capacitor.

Figure S1: Schematic of the left-handed transmission line circuit.

Let , , and be constants, and take a harmonic mode with angular frequency . For , we apply Kirchhoff’s laws to the inductors and capacitors, with the phasor convention, and obtain


Combining these to eliminate and yields the following pair of coupled equations:


where , and .

Suppose we close the circuit by grounding the leftmost and rightmost sites. Consider the left edge (the right edge is handled similarly). There, the Kirchhoff equations simplify to


resulting in the boundary equation


Hence, we arrive at the modified SSH problem discussed in the main text:


This is the configuration referred to in the main text as the “nontrivial lattice”. The “trivial lattice” can be described by removing the first row and column of the matrix. In either case, the edges of the band gap are


and the angular frequency of the edge state is

Figure S2: Response of the linear circuit to an external voltage, calculated numerically from the circuit equations with the parameters and . (a)–(b) Resonant response of (a) the amplitude of the voltage on the first site (to the right of the leftmost capacitor), and (b) the amplitude of the stored charge on the first (leftmost) capacitor, versus source frequency , for . The resonance frequency predicted by Eq. (18) is indicated by the horizontal dashed line. (c) Spatial distribution of the voltage ampltidues for and . (d) Intensity of the nontrivial lattice’s resonant mode, as given by the value of (where is the voltage amplitude at site ) at resonance, versus the parameter. For each value of , the input voltage has amplitude and frequency given by Eq. (18).

Next, we consider the response of the circuit to a harmonic voltage source. Instead of grounding the left edge, we apply an input voltage of amplitude and frequency . Then Eq. (15) is replaced by


and the eigenvalue equation (16) is replaced by an inhomogenous equation.

Numerical solutions for this problem are shown in Fig. S2. Fig. S2(a)–(b) shows that the voltage amplitude in the nontrivial lattice is resonantly enhanced when matches . The trivial lattice, on the other hand, does not exhibit a resonant enhancement. However, when we plot the voltage distributions, they have the same decay constant, as shown in Fig. S2(c). This is due to the fact that they have the same bulk Hamiltonians.

Fig. S2(d) shows the modal intensity (the sum of squared voltage amplitudes over the lattice) at resonance, versus the parameter, for a nontrivial lattice. In this plot, the input frequency for each value of is adjusted to the edge state frequency given by Eq. (18). With increasing , the bandgap becomes larger (i.e., the lattice moves deeper into the topologically nontrivial phase); accordingly, the edge state is more strongly confined, and responds more strongly to the resonant excitation. The modal intensity scales exponentially with the bandgap size.

Appendix S2 Nonlinear circuit

s2.1 Circuit equations

We seek a set of time-domain equations for the circuit’s nonlinear regime, where the capacitors are nonlinear. Let and denote the charges stored in capacitors and , respectively, on site . These obey


Here, is the value of the nonlinear capacitance in unit cell . Using Kirchhoff’s laws, we can derive several additional equations. The time-dependent voltage-current relations on the inductors are


The current-charge relations on the capacitors, with the assumptions of current conservation and zero net charge, give


By combining Eqs. (20)–(25), we can eliminate the and variables, resulting in the following pair of time-domain circuit equations expressed in terms of the variables:




is the nonlinear capacitance ratio at site .

It is convenient to re-cast Eqs. (26)–(27) in terms of the variables




s2.2 Harmonic decomposition and nonlinearity model

If a harmonic signal is injected into the nonlinear circuit, higher harmonics are generated. Due to the symmetric C-V curve of the nonlinear capacitors, even-order harmonics are suppressed.

Let denote the frequency of the first harmonic. We will decompose the voltage variables in the following way:


On the right hand sides, the integer superscripts denote the harmonic index. The factor of is for later convenience; we expect the first harmonic mode to behave like an SSH edge state, which is characterized by alternating signs on adjacent unit cells (another way of saying this is that the band gap of the bulk SSH model is narrowest at the Brillouin zone boundary, , where is the lattice constant), and this factor ensures that the and variables act as smooth envelopes with the sign alternation taken out.

We now have to substitute the ansatz (33)–(34) into Eqs. (31)–(32). First, consider Eq. (31), which is easy to deal with since it is linear. Matching the individual harmonics, we obtain




Next, consider the nonlinear equation (32). The main complication here is the term involving


The time variation of gives rise to two classes of effects: (i) self-phase modulation and cross-phase modulation, which alter the effective value of “seen” by each given harmonic, and (ii) frequency-mixing processes, which couple the dynamical equations for the different harmonics. For now, let us try to pick out the contributions to category (i), neglecting (ii).

In our experiment, each nonlinear capacitor consists of a pair of back-to-back varactors. The nonlinear capacitance ratio was defined in Eq. (28). Let us make the assumption that


Here, is the capacitance ratio in the linear limit and is a Kerr-like parameter determining the strength of the lowest-order nonvanishing (cubic) nonlinearity. To obtain values for and , we use the manufacturer-supplied capacitance-voltage curve for the individual varactors to calculate and the bias voltage for a pair of back-to-back varactors. We then perform a linear least-squares fit of versus , using the subset of data points with voltage biases . The fitted parameters are and , and the fit is shown in Fig. 1(c) of the main text.

In (37), we can take the approximation of replacing with its time-independent part,


For each harmonic , this would then give rise to a term


with now playing the role of an ”effective” parameter.

This approximation does not capture all possible self-phase and cross-phase modulation terms. This can be seen in the low-intensity limit, where we can Taylor expand in the variables; in this expansion, there will be non-constant terms like , which couples to the harmonic term from to yield a term proportional to , and hence contributing to the self-phase or cross-phase modulation. We will not undertake a rigorous analysis of these terms, since the Taylor expansion is invalid anyway when the intensities are not small. Instead, our take-home message is as follows:

  1. Each harmonic contributes to the effective value of in direct proportion to its local intensity, like in Eq. (39).

  2. However, the precise numerical factor need not be exactly the same as in Eq. (39).

Based on this approximation, we can now deal with the nonlinear equation (32), which simplifies to


s2.3 Localized and traveling-wave solutions

Let us consider the case where is approximately constant in space, and look for solutions of the form


These are traveling-wave solutions if is real, and exponentially localized solutions if is complex. Substituting this into Eqs. (35) and (41) gives


Solving the characteristic equation gives


We can then easily show that, for , the domains over which the right-hand side has magnitude smaller than unity (i.e., is real) are:


For , this corresponds exactly to the bands shown in Fig. 1(d)–(e) of the main text. In particular, within the band gap between and , the right-hand side is larger than unity and hence is imaginary, in complete agreement with the linear analysis.

For the higher-harmonic modes, (45) is satisfied easily. For example, for the third harmonic, we require


For operating frequencies below the linear-regime band gap, , this is satisfied for


which is well within the regime considered in this experiment. This analysis thus confirms that the nonlinear circuit is capable of supporting traveling-wave higher-harmonic solutions.

s2.4 SPICE simulation results: voltage profiles

We performed simulations of the nonlinear circuit using the circuit simulation software SPICE. The simulations reproduce the basic features of the experimental results, though the results are not in exact agreement, probably due to imperfections in the circuit components.

Figure S3: On-site voltage amplitudes for the first-harmonic signal (orange circles) and third-harmonic signal (blue squares). The frequency of the input is . (a)–(b) Nontrivial and trivial lattices with . (c)–(d) Nontrivial and trivial lattices with .

Fig. S3 shows the on-site voltage amplitudes for the first- and third-harmonic signals. These are extracted from the simulation results in a manner similar to the experiment: after the simulation reaches steady-state, we take a time-dependent sample, Fourier transform, and extract the peak heights. To obtain simulation results matching the experimental results shown in Fig. 2(a)–(b) of the main text, we find that it is necessary to apply a higher input voltage amplitude than in the experiment, . The results are shown in Fig. S3(c)–(d). Similar to the experiment, the first-harmonic mode in the nontrivial lattice decays away from the edge, reaching values much lower than in the trivial lattice.

Figure S4: Squared bias voltage amplitudes on the nonlinear capacitors, for the first-harmonic signal (orange circles) and higher-harmonic signals (purple squares). The frequency of the input is . The higher-harmonic data is obtained by summing over the third, fifth, and seventh harmonic data (further harmonics are negligible). (a)–(b) Nontrivial and trivial lattices with . (c)–(d) Nontrivial and trivial lattices with .

Fig. S4 shows the bias voltage amplitudes on the nonlinear capacitors, which were not measured in the experiment. As discussed in Section S2.2, the bias voltages determine the effective value of the nonlinear parameter. To obtain this data from the simulations, we extract the time-dependent bias voltage samples (i.e., the time-dependent voltages between the ports of the nonlinear capacitors, denoted by in Section S2.1), Fourier transform, and extract the peak heights; this yields the components denoted by in Section S2.2. According to Eq. (39), the contribution of each harmonic to the effective local is proportional to . Fig. S4 shows a comparison between the first-harmonic contribution (orange circles) and the higher-harmonic contributions (purple squares). In particular, in Fig. S4(c), which correspond to the voltage plot of Fig. S3(c), the higher-harmonic signals are found to increasingly dominate the nonlinearity as we go deeper into the lattice.

s2.5 SPICE simulation results: comparison with conventional NLTL

Finally, we used SPICE simulations to compare the third harmonic intensity in this circuit to a conventional left-handed NLTL. In the conventional left-handed NLTL, all the linear capacitors are replaced with nonlinear capacitors (i.e., the lattice is no longer dimerized). In the linear limit, the conventional left-handed NLTL with is a high-pass filter with a Bragg cutoff frequency of .

Fig. S5 plots the simulation results for the normalized third-harmonic intensity (defined in the same way as in the main text), versus the input parameters and input voltage . The simulation results for the SSH-like lattice, shown in Fig. S5(a), are similar to the experimental results shown in Fig. 3(a) of the main text; in particular, the maximum value of is , comparable to the experimentally-obtained maximum value . By contrast, Fig. S5(b) shows that the conventional NLTL exhibits no comparable enhancement of the third-harmonic signal intensity, with throughout the entire parameter regime we investigated. Hence, the introduction of the topological edge mode has contributed to a five-fold increase in the intensity of the generated third-harmonic signal.

Figure S5: Plot of the normalized third-harmonic intensity versus input frequency and input voltage , for (a) the SSH-like lattice in its topologically nontrivial configuration, and (b) a conventional left-handed nonlinear transmission line (NLTL) with identical nonlinear capacitors. The figure of merit is defined in the same way as in the main text, as the mean squared third-harmonic amplitude on the first 10 sites relative to the squared input voltage amplitude. In the SSH case, we obtain , which is comparable to the experimental result of reported in the main text. In the conventional NLTL, we observe throughout the parameter range.
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description