Skip to main content
  • Research article
  • Open access
  • Published:

Fast and accurate simulations of transmission-line metamaterials using transmission-matrix method


Recently, two-dimensional (2D) periodically L and C loaded transmission-line (TL) networks have been applied to represent metamaterials. The commercial Agilent's Advanced Design System (ADS) is a commonly-used tool to simulate the TL metamaterials. However, it takes a lot of time to set up the TL network and perform numerical simulations using ADS, making the metamaterial analysis inefficient, especially for large-scale TL networks. In this paper, we propose transmission-matrix method (TMM) to simulate and analyze the TL-network metamaterials efficiently. Compared to the ADS commercial software, TMM provides nearly the same simulation results for the same networks. However, the model-process and simulation time has been greatly reduced. The proposed TMM can serve as an efficient tool to study the TL-network metamaterials.

PACS Codes:41.20.Jb


In 1968, Veselago first proposed the concept of left-handed material (LHM) [1]. The new medium exhibits many particular properties because of its negative permittivity and negative permeability. However, LHM does not exist in nature, which has to be realized using periodic structures with the unit cell of split-ring resonators (SRR) and conducting wires [25] or electric LC (ELC) resonators [6]. The two-dimensional (2D) L and C loaded transmission-line (TL) networks is another method to realize LHM [717]. Both isotropic [714] and anisotropic TL metamaterials [1517] have been studied. The particular 2D L and C loaded TL networks can be equivalent to relevant media. In order to obtain the voltage and current distributions on the periodically TL networks, which are equivalent to the field distributions in metamaterial, the Agilent's Advanced Design System (ADS) has been used [817]. If we want to observe the node voltages and currents of the LC-loaded TL networks when using ADS, the observation nodes should be labeled in ADS, which is extremely time and memory consuming for large-scale circuit structures. The numerical simulations were time consuming too for large-scale structures using ADS.

Are there some methods which can be used to set up the model quickly and to simulate the TL strucutures efficiently? The transmission-matrix method (TMM) is a good choice [18]. TMM can construct the LC-loaded networks automatically and quickly. The most important thing is that one does not need to label the observation nodes. Hence TMM can save much time and vigor in the construction of TL metamaterials and is more efficient than ADS. Caloz and Itoh have studied such a method in Ref. [19] to analyze TL metamaterials. However, they only considered the periodic structures composed of lumped inductance and capacitance, the composed right/left-handed circuits, in which the transmission lines are not involved [19]. Also, the voltage and current distributions obtained from this method cannot be equivalent to those from the ADS simulations, because the load impedances of observers are added to observation nodes.

In this paper, we generalize the TMM approach to study the 2D LC-loaded TL periodic structures, which have been widely used to represent metamaterials [717]. Both the voltage and current distributions on the LC-loaded TL structures can be computed accurately and efficiently using the proposed TMM, which have excellent agreements to the ADS simulation results. However, the time used in the model setup and numerical simulations have been greatly reduced. The proposed TMM approach can also solve large-scale TL metamaterials, which can hardly be realized using ADS.

Theoretical analysis of TMM

The arbitrary 2D TL network under consideration is divided into m × n unit cells, as shown in Fig. 1(a). The corresponding general T structure of each unit cell is illustrated in Fig. 1(b), in which Z0, β and l are the characteristic impedance, the propagation constant, and the length of TL section, respectively; and Z s and Y p are loaded inductance and capacitance. Using TMM, we divide the structure into columns of unit cells, one of which is shown as the gray line in Fig. 1(a). The m × 1 unit cells can be represented as a 2m × 2m matrix formed by m input ports and m output ports with the node voltages and port currents [19]. The equivalent column cell with terminated impedance R j is shown in Fig. 2. The input and output ports are indexed as 1 to m and m + 1 to 2m, respectively, and the corresponding 2m × 2m matrix is expressed as

Figure 1
figure 1

(a) Arbitrary 2D L and C loaded TL network. (b) General T network of each unit cell.

Figure 2
figure 2

The equivalent-circuit T network for a column cell and the corresponding decomposition. (a) T network for a column of cell. (b) The decomposition of three parts.

[ [ V i n ] [ I i n ] ] = [ T j ] [ [ V o u t ] [ I o u t ] ] , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaWaamWaaeaafaqabeGabaaabaGaei4waSLaemOvay1aaSbaaSqaaiabdMgaPjabd6gaUbqabaGccqGGDbqxaeaacqGGBbWwcqWGjbqsdaWgaaWcbaGaemyAaKMaemOBa4gabeaakiabc2faDbaaaiaawUfacaGLDbaacqGH9aqpcqGGBbWwcqWGubavdaWgaaWcbaGaemOAaOgabeaakiabc2faDjabgwSixpaadmaabaqbaeqabiqaaaqaaiabcUfaBjabdAfawnaaBaaaleaacqWGVbWBcqWG1bqDcqWG0baDaeqaaOGaeiyxa0fabaGaei4waSLaemysaK0aaSbaaSqaaiabd+gaVjabdwha1jabdsha0bqabaGccqGGDbqxaaaacaGLBbGaayzxaaGaeiilaWcaaa@5723@

which can be referred to Eq. (4.18) in Ref. [19].

By decomposing the column cell into three sub-column cells, as shown in Fig. 2, we can easily obtain the 2m × 2m transmission matrix [T j ] as the product of three matrices [Th1], [Tv] and [Th2]:

[ T j ] = [ [ A ] [ B ] [ C ] [ D ] ] = [ T h 1 ] [ T v ] [ T h 2 ] , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaei4waSLaemivaq1aaSbaaSqaaiabdQgaQbqabaGccqGGDbqxcqGH9aqpdaWadaqaauaabeqaciaaaeaacqGGBbWwcqWGbbqqcqGGDbqxaeaacqGGBbWwcqWGcbGqcqGGDbqxaeaacqGGBbWwcqWGdbWqcqGGDbqxaeaacqGGBbWwcqWGebarcqGGDbqxaaaacaGLBbGaayzxaaGaeyypa0Jaei4waSLaemivaqLaemiAaGMaeGymaeJaeiyxa0LaeyyXICTaei4waSLaemivaq1aaSbaaSqaaiabdAha2bqabaGccqGGDbqxcqGHflY1cqGGBbWwcqWGubavcqWGObaAcqaIYaGmcqGGDbqxcqGGSaalaaa@5A73@

where the computation of matrices [Th1], [Tv] and [Th2] is similar to that in Ref. [19].

Once [T j ] has been obtained, the n-column transmission matrix as shown in Fig. 1(a) will be determined as

[ T ] = j = 1 n [ T j ] . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaei4waSLaemivaqLaeiyxa0Laeyypa0ZaaebCaeaacqGGBbWwcqWGubavdaWgaaWcbaGaemOAaOgabeaakiabc2faDbWcbaGaemOAaOMaeyypa0JaeGymaedabaGaemOBa4ganiabg+GivdGccqGGUaGlaaa@3DF5@

Voltage and current distributions

We consider the case of an m × n network with a given source which is arbitrarily loacted at the kth row and between the ith and (i + 1)th columns of the network. We determine to obtain the output voltage and current at any point of the network. The whole network is demonstrated in Fig. 3, where R 1 , R 2 , , and  R m MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafmOuaiLbauaadaWgaaWcbaGaeGymaedabeaakiabcYcaSiqbdkfaszaafaWaaSbaaSqaaiabikdaYaqabaGccqGGSaalcqWIVlctcqGGSaalcaaMc8UaaGPaVlabbggaHjabb6gaUjabbsgaKjabbccaGiqbdkfaszaafaWaaSbaaSqaaiabd2gaTbqabaaaaa@3FC3@ are horizontally terminated impedances.

Figure 3
figure 3

The whole network under consideration, including m inputs and m outputs, which is divided into the input network T in , the center network T c , and the output network T out . A source is arbitrarily located at the kth row and between the ith and (i + 1)th columns of the network.


[ T i n ] = [ T i n ] [ T R ] = [ [ A i n ] [ B i n ] [ C i n ] [ D i n ] ] , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaei4waSLafmivaqLbauaadaWgaaWcbaGaemyAaKMaemOBa4gabeaakiabc2faDjabg2da9iabcUfaBjabdsfaunaaBaaaleaacqWGPbqAcqWGUbGBaeqaaOGaeiyxa0LaeyyXICTaei4waSLaemivaq1aaSbaaSqaaiabdkfasbqabaGccqGGDbqxcqGH9aqpdaWadaqaauaabeqaciaaaeaacqGGBbWwcuWGbbqqgaqbamaaBaaaleaacqWGPbqAcqWGUbGBaeqaaOGaeiyxa0fabaGaei4waSLafmOqaiKbauaadaWgaaWcbaGaemyAaKMaemOBa4gabeaakiabc2faDbqaaiabcUfaBjqbdoeadzaafaWaaSbaaSqaaiabdMgaPjabd6gaUbqabaGccqGGDbqxaeaacqGGBbWwcuWGebargaqbamaaBaaaleaacqWGPbqAcqWGUbGBaeqaaOGaeiyxa0faaaGaay5waiaaw2faaiabcYcaSaaa@6009@
[ T o u t ] = [ T c ] [ T o u t ] [ T R ] = [ [ A o u t ] [ B o u t ] [ C o u t ] [ D o u t ] ] , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaei4waSLafmivaqLbauaadaWgaaWcbaGaem4Ba8MaemyDauNaemiDaqhabeaakiabc2faDjabg2da9iabcUfaBjabdsfaunaaBaaaleaacqWGJbWyaeqaaOGaeiyxa0LaeyyXICTaei4waSLaemivaq1aaSbaaSqaaiabd+gaVjabdwha1jabdsha0bqabaGccqGGDbqxcqGHflY1cqGGBbWwcqWGubavdaWgaaWcbaGaemOuaifabeaakiabc2faDjabg2da9maadmaabaqbaeqabiGaaaqaaiabcUfaBjqbdgeabzaafaWaaSbaaSqaaiabd+gaVjabdwha1jabdsha0bqabaGccqGGDbqxaeaacqGGBbWwcuWGcbGqgaqbamaaBaaaleaacqWGVbWBcqWG1bqDcqWG0baDaeqaaOGaeiyxa0fabaGaei4waSLafm4qamKbauaadaWgaaWcbaGaem4Ba8MaemyDauNaemiDaqhabeaakiabc2faDbqaaiabcUfaBjqbdseaezaafaWaaSbaaSqaaiabd+gaVjabdwha1jabdsha0bqabaGccqGGDbqxaaaacaGLBbGaayzxaaGaeiilaWcaaa@70CB@

in which

[ T R ] = [ [ I ] [ X ] [ O ] [ I ] ] , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaei4waSLaemivaq1aaSbaaSqaaiabdkfasbqabaGccqGGDbqxcqGH9aqpdaWadaqaauaabeqaciaaaeaacqGGBbWwcqWGjbqscqGGDbqxaeaacqGGBbWwcqWGybawcqGGDbqxaeaacqGGBbWwcqWGpbWtcqGGDbqxaeaacqGGBbWwcqWGjbqscqGGDbqxaaaacaGLBbGaayzxaaGaeiilaWcaaa@43B5@


[ X ] = [ R 1 0 0 0 R 2 0 0 0 0 R m ] . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaei4waSLaemiwaGLaeiyxa0Laeyypa0ZaamWaaeaafaqabeabeaaaaaqaaiqbdkfaszaafaWaaSbaaSqaaiabigdaXaqabaaakeaacqaIWaamaeaacqWIVlctaeaacqaIWaamaeaacqaIWaamaeaacuWGsbGugaqbamaaBaaaleaacqaIYaGmaeqaaaGcbaGaeS47IWeabaGaeGimaadabaGaeSO7I0eabaGaeSO7I0eabaGaeSy8I8eabaGaeSO7I0eabaGaeGimaadabaGaeGimaadabaGaeGimaadabaGaemOuai1aaSbaaSqaaiabd2gaTbqabaaaaaGccaGLBbGaayzxaaGaeiOla4caaa@4D79@

In Eqs. (4) and (5), [ A i n ] , [ B i n ] , [ C i n ] , [ D i n ] , [ A o u t ] , [ B o u t ] , [ C o u t ] , [ D o u t ] MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaei4waSLafmyqaeKbauaadaWgaaWcbaGaemyAaKMaemOBa4gabeaakiabc2faDjabcYcaSiabcUfaBjqbdkeaczaafaWaaSbaaSqaaiabdMgaPjabd6gaUbqabaGccqGGDbqxcqGGSaalcqGGBbWwcuWGdbWqgaqbamaaBaaaleaacqWGPbqAcqWGUbGBaeqaaOGaeiyxa0LaeiilaWIaei4waSLafmiraqKbauaadaWgaaWcbaGaemyAaKMaemOBa4gabeaakiabc2faDjabcYcaSiabcUfaBjqbdgeabzaafaWaaSbaaSqaaiabd+gaVjabdwha1jabdsha0bqabaGccqGGDbqxcqGGSaalcqGGBbWwcuWGcbGqgaqbamaaBaaaleaacqWGVbWBcqWG1bqDcqWG0baDaeqaaOGaeiyxa0LaeiilaWIaei4waSLafm4qamKbauaadaWgaaWcbaGaem4Ba8MaemyDauNaemiDaqhabeaakiabc2faDjabcYcaSiabcUfaBjqbdseaezaafaWaaSbaaSqaaiabd+gaVjabdwha1jabdsha0bqabaGccqGGDbqxaaa@6CA1@ and [X] are all m × m matrices.

Using the series transmission-matrix theory, we can easily obtain the following equations

[ [ V i n ] [ I 11 ] ] = [ T i n ] [ [ O ] [ I ] ] , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaWaamWaaeaafaqabeGabaaabaGaei4waSLaemOvay1aaSbaaSqaaiabdMgaPjabd6gaUbqabaGccqGGDbqxaeaacqGGBbWwcqWGjbqsdaWgaaWcbaGaeGymaeJaeGymaedabeaakiabc2faDbaaaiaawUfacaGLDbaacqGH9aqpcqGGBbWwcuWGubavgaqbamaaBaaaleaacqWGPbqAcqWGUbGBaeqaaOGaeiyxa0LaeyyXIC9aamWaaeaafaqabeGabaaabaGaei4waSLaem4ta8Kaeiyxa0fabaGaei4waSLaemysaKKaeiyxa0faaaGaay5waiaaw2faaiabcYcaSaaa@4EA2@
[ [ V i n ] [ I 12 ] ] = [ T o u t ] [ [ O ] [ I ] ] , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaWaamWaaeaafaqabeGabaaabaGaei4waSLaemOvay1aaSbaaSqaaiabdMgaPjabd6gaUbqabaGccqGGDbqxaeaacqGGBbWwcqWGjbqsdaWgaaWcbaGaeGymaeJaeGOmaidabeaakiabc2faDbaaaiaawUfacaGLDbaacqGH9aqpcqGGBbWwcuWGubavgaqbamaaBaaaleaacqWGVbWBcqWG1bqDcqWG0baDaeqaaOGaeiyxa0LaeyyXIC9aamWaaeaafaqabeGabaaabaGaei4waSLaem4ta8Kaeiyxa0fabaGaei4waSLafmysaKKbauaacqGGDbqxaaaacaGLBbGaayzxaaGaeiilaWcaaa@503B@

here, [O] is m × 1 zero matrix, [I] and [I'] are the current as illustrates in Fig. 3. From Eqs. (8) and (9), we have

[ I 11 ] = [ D i n ] [ B i n ] 1 [ V i n ] = [ Y i n ] [ V i n ] , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaei4waSLaemysaK0aaSbaaSqaaiabigdaXiabigdaXaqabaGccqGGDbqxcqGH9aqpcqGGBbWwcuWGebargaqbamaaBaaaleaacqWGPbqAcqWGUbGBaeqaaOGaeiyxa0LaeyyXICTaei4waSLafmOqaiKbauaadaWgaaWcbaGaemyAaKMaemOBa4gabeaakiabc2faDnaaCaaaleqabaGaeyOeI0IaeGymaedaaOGaeyyXICTaei4waSLaemOvay1aaSbaaSqaaiabdMgaPjabd6gaUbqabaGccqGGDbqxcqGH9aqpcqGGBbWwcuWGzbqwgaqbamaaBaaaleaacqWGPbqAcqWGUbGBaeqaaOGaeiyxa0LaeyyXICTaei4waSLaemOvay1aaSbaaSqaaiabdMgaPjabd6gaUbqabaGccqGGDbqxcqGGSaalaaa@5EE7@
[ I 12 ] = [ D o u t ] [ B o u t ] 1 [ V i n ] = [ Y o u t ] [ V i n ] , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaei4waSLaemysaK0aaSbaaSqaaiabigdaXiabikdaYaqabaGccqGGDbqxcqGH9aqpcqGGBbWwcuWGebargaqbamaaBaaaleaacqWGVbWBcqWG1bqDcqWG0baDaeqaaOGaeiyxa0LaeyyXICTaei4waSLafmOqaiKbauaadaWgaaWcbaGaem4Ba8MaemyDauNaemiDaqhabeaakiabc2faDnaaCaaaleqabaGaeyOeI0IaeGymaedaaOGaeyyXICTaei4waSLaemOvay1aaSbaaSqaaiabdMgaPjabd6gaUbqabaGccqGGDbqxcqGH9aqpcqGGBbWwcuWGzbqwgaqbamaaBaaaleaacqWGVbWBcqWG1bqDcqWG0baDaeqaaOGaeiyxa0LaeyyXICTaei4waSLaemOvay1aaSbaaSqaaiabdMgaPjabd6gaUbqabaGccqGGDbqxcqGGSaalaaa@638A@


[I in ] = [I11] + [I12]

Furthermore, we obtain

[ V i n ] = ( [ Y i n ] + [ Y o u t ] ) 1 [ I i n ] = [ Y ] [ I i n ] . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaei4waSLaemOvay1aaSbaaSqaaiabdMgaPjabd6gaUbqabaGccqGGDbqxcqGH9aqpcqGGOaakcqGGBbWwcuWGzbqwgaqbamaaBaaaleaacqWGPbqAcqWGUbGBaeqaaOGaeiyxa0Laey4kaSIaei4waSLafmywaKLbauaadaWgaaWcbaGaem4Ba8MaemyDauNaemiDaqhabeaakiabc2faDjabcMcaPmaaCaaaleqabaGaeyOeI0IaeGymaedaaOGaeyyXICTaei4waSLaemysaK0aaSbaaSqaaiabdMgaPjabd6gaUbqabaGccqGGDbqxcqGH9aqpcqGGBbWwcqWGzbqwcqGGDbqxcqGHflY1cqGGBbWwcqWGjbqsdaWgaaWcbaGaemyAaKMaemOBa4gabeaakiabc2faDjabc6caUaaa@5EDC@

According to Eq. (4.32) in Ref. [19], the input current [I in ]m× 1 will be zero except the kth-row current Iin, k. The kth-row input voltage Vin, k= V s - R s ·Iin, k. Then we have

I i n , k = V s R s + y k k , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemysaK0aaSbaaSqaaiabdMgaPjabd6gaUjabcYcaSiabdUgaRbqabaGccqGH9aqpjuaGdaWcaaqaaiabdAfawnaaBaaabaGaem4CamhabeaaaeaacqWGsbGudaWgaaqaaiabdohaZbqabaGaey4kaSIaemyEaK3aaSbaaeaacqWGRbWAcqWGRbWAaeqaaaaakiabcYcaSaaa@3FBF@

in which y kk is the kth-row and kth-column element of matrix [Y].

We remark that the kth element of the matrix [I in ] is Iin, k, and all other elements are zero. Hence, once we obtain Iin, k, the matrix [I in ] will be known, and then matrices [V in ], [I11] and [I12] will also be determined from Eqs. (13), (10) and (11), respectively.

Once we know [V in ], [I11] and [I12], the output voltage and current at any point of the network which is located on the right side of the source will be determined easily from the following equation

[ [ V o u t ] [ I o u t ] ] = [ T C ] 1 [ [ V i n ] [ I 12 ] ] . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaWaamWaaeaafaqabeGabaaabaGaei4waSLaemOvay1aaSbaaSqaaiabd+gaVjabdwha1jabdsha0bqabaGccqGGDbqxaeaacqGGBbWwcqWGjbqsdaWgaaWcbaGaem4Ba8MaemyDauNaemiDaqhabeaakiabc2faDbaaaiaawUfacaGLDbaacqGH9aqpcqGGBbWwcqWGubavdaWgaaWcbaGaem4qameabeaakiabc2faDnaaCaaaleqabaGaeyOeI0IaeGymaedaaOGaeyyXIC9aamWaaeaafaqabeGabaaabaGaei4waSLaemOvay1aaSbaaSqaaiabdMgaPjabd6gaUbqabaGccqGGDbqxaeaacqGGBbWwcqWGjbqsdaWgaaWcbaGaeGymaeJaeGOmaidabeaakiabc2faDbaaaiaawUfacaGLDbaacqGGUaGlaaa@580F@

Similarly, we can obtain the output voltages and currents at points located on the left side of the source. Using the same method, we can also deal with the cases of two given sources which are located at the same row.

The above method can be used to analyze the electromagnetic properties of 2D TL metamaterials such as super lens [711], EM localizations [12, 13], and super waveguides [14] with high efficiencies and high accuracies.

Simulation results

In order to verify the correctness of the proposed TMM method, we re-compute the sub-wavelength focusing structure possessing the same TL network and parameters as those in Ref. [9] using TMM. The corresponding TMM-simulation results are shown in Figs. 4 and 5. Comparing with the ADS simulations [9], we notice that two results have excellent agreements in both amplitude and phase distributions. The results show that it is valid to replace ADS by TMM in numerical simulations of TL metamaterial structures. In this example, TMM does not exhibit its advantage obviously in saving time with the small-scale network.

Figure 4
figure 4

TMM-simulation results for TL super lens [9] at 1 GHz (voltage amplitude.

Figure 5
figure 5

TMM-simulation results for TL super lens [9] at 1 GHz (voltage phase).

To verify the accuracy of TMM further, we give a full comparison of TMM-simulation results to the commercial ADS-simulation results. We consider an example of super waveguide, which is a planar waveguide filled with air and left-handed material (LHM) [14, 20]. It has been shown that extremely high-power densities with opposite propagation directions can be generated and transmitted along the waveguide if the air and LHM have equal thickness and the permittivity and permeability of LHM satisfy ε = ε0(1 + δ) and μ = -μ0/(1 + δ), in which δ is a small parameter [20]. The above theoretical prediction has been realized using 2D LC-loaded TL metamaterials [14]. Now we re-analysis the TL-metamaterial super waveguide using the TMM technique.

Consider a TL-metamaterial super waveguide which is composed of 179 × 18 unit cells [14], where the right-handed TL (RHTL) and the left-handed TL (LHTL) parts extend 9 cells in the x direction and 179 cells in the y direction, respectively. A 1-V (0 dB) voltage source is connected to the node of the cell numbered as (4, 90) in the RHTL region. We choose the size of unit cell as d = 1 cm, the small parameter as δ = -0.053, the characteristic impedance in TL as Z0 = 533.1459 Ω, and the propagation constant in TL as β = 14.8096. The simulation results of voltage and current distributions computed by TMM are illustrated in Fig. 6. As a comparison, the corresponding voltage distributions simulated by ADS are shown in Fig. 7. From these two figures, excellent agreements have been observed for the voltage distributions. The current distributions have also excellent agreements, which are not repeated here.

Figure 6
figure 6

The voltage and current distributions on the planar TL-metamaterial super waveguide with 18 × 179 unit cells computed by TMM. (a) The amplitude distribution of voltage. (b) The phase distribution of voltage. (c) The amplitude distribution of current. (d) The phase distribution of current.

Figure 7
figure 7

The voltage and current distributions on the planar TL-metamaterial super waveguide with 18 × 179 unit cells simulated by the ADS commercial software. (a) The amplitude distribution of voltage. (b) The phase distribution of voltage.

We remark that much less time has been used in TMM than that in ADS, since it takes long time to label the voltage and current probes in ADS for the 179 × 18 network. At least several days were taken to set up and simulate such a large-scale TL metamaterial using ADS, while it only takes several minutes to solve the same problem using TMM.

To further verify the accuracy of TMM, the voltage and current distributions on the planar TL-metamaterial super waveguide along the x direction with y = 90 are computed using TMM and ADS, as shown in Fig. 8. Clearly, the two results have excellent agreements.

Figure 8
figure 8

The comparison of voltage and current distributions on the planar TL-metamaterial super waveguide along the line y = 90, where red lines indicate the ADS simulation results, and the dashed blue lines are TMM computation results. (a) The amplitude of voltage. (b) The phase of voltage. (c) The amplitude of current. (d) The phase of current.

Then we consider a longer TL-metamaterial super waveguide, which is composed of 599 × 18 unit cells. Similar to the previous case, the RHTL and LHTL parts extend 9 cells in the x direction but 599 cells in the y direction. A 1-V (0-dB) voltage source is located at the node numbered as (4, 300). The TMM-simulation results for the voltage and current distributions are shown in Fig. 9. We remark that ADS cannot be used to handle the problem in a 2GB-memory personal computer due to the large requirement of memory.

Figure 9
figure 9

The voltage and current distributions of a longer TL-metamaterial super waveguide with 18 × 599 unit cells computed by TMM. For the same structure, ADS cannot be used due to the memory limit of personal computer. (a) The amplitude distribution of voltage. (b) The phase distribution of voltage. (c) The amplitude distribution of current. (d) The phase distribution of current.

In another example, we consider an LHTL superlens and use it to localize electromagnetic waves and energies [12, 13]. To build up the TL-metamaterial superlens [12], 33 × 21 unit cells have been used, where the left RHTL part, the middle LHTL part, and the right RHTL part extend to 7 cells in the x direction and 33 unit cells in the y direction, respectively. Two 1-V (0 dB) voltage sources V1 and V2 with opposite phases are connected to the nodes numbered as (17, 4) and (17, 17) in the RHTL regions, respectively. Here, we chose d = 1 mm, δ = 10-4, the characteristic impedance Z0 = 533.1459 Ω, the propagation constant β = 14.8096, the series capacitance C = 10.0773 pF, and the shunt inductance L = 1432.2120 nH. The TMM-simulation results of the voltage and current distributions are demonstrated in Fig. 10, and the corresponding ADS-simulation results are shown in Fig. 11. Comparing Figs. 10 and 11, such two results have excellent agreements.

Figure 10
figure 10

The voltage and current distributions on the TL-metamaterial superlens for two sources computed by TMM, where 21 × 33 unit cells have been used. (a) The amplitude distribution of voltage. (b) The phase distribution of voltage. (c) The amplitude distribution of current. (d) The phase distribution of current.

Figure 11
figure 11

The voltage and current distributions on the TL-metamaterial superlens for two sources simulated by ADS, where 21 × 33 unit cells have been used. (a) The amplitude distribution of voltage. (b) The phase distribution of voltage. (c) The amplitude distribution of current. (d) The phase distribution of current.

Similar to the earlier example, we have observed the voltage and current distributions computed by the two methods along a line in the x direction with y = 17, as illustrated in Fig. 12. The results demonstrate that TMM and ADS provide nearly the same accuracy.

Figure 12
figure 12

The amplitudes and phases of voltages and currents along the line y = 17 in the superlens for two sources, where the red lines are ADS simulation results and the dashed blue lines are TMM computation results. (a) The amplitudes of voltages. (b) The phases of voltages. (c) The amplitudes of currents. (d) The phases of currents.

We also consider a large-scale TL-metamaterial superlens for two sources with 599 × 21 unit cells, which cannot be solved using ADS. In this example, two 1-V (0 dB) voltage sources with opposite phases are connect to the nodes numbered as (4, 300) and (17, 300) in the RHTL regions, respectively. The TMM-simulation results of the voltage and current distributions are illustrated in Fig. 13.

Figure 13
figure 13

The voltage and current distributions on the long TL-metamaterial superlens for two sources computed by TMM, where 21 × 599 unit cells have been used. (a) The amplitude distribution of voltage. (b) The phase distribution of voltage. (c) The amplitude distribution of current. (d) The phase distribution of current.


In this paper, a TMM approach has been proposed to determine the voltage and current distributions on 2D LC-loaded TL metamaterials. The proposed TMM has the same accuracy as ADS for numerical simulations of finite-size metamaterial structures. However, it is much easier to set up LC-loaded TL models and save much more time using TMM than ADS. For large-scale metamaterial structures, ADS can even not be used due to the large memory requirements. As an example, for the large-scale TL-metamaterial super waveguide [14] with 179 × 18 unit cells, whose results are shown in Figs. 6 and 7, the total computation time is only 36.543 seconds using TMM in a 2G-memory personal computer. The method is realized by MATLAB with computing 18 358 × 358 matrices. We also realized the same super waveguide using ADS. It took us nearly two weeks to set up the circuit model and label the observation nodes in the same 2G-memory personal computer. The boring work is labeling the observation nodes, which must be very careful and to avoid mistakes. When the circuit structure is constructed, the total simulation time using ADS is about fifteen minutes. Compared to ADS, TMM is much more efficient when treating particular circuits with large amount of same unit cells such as the TL metamaterials.


  1. Veselago VG: The electrodynamics of substances with simultaneously negative values of ε and μ. Sov Phys-Usp. 1968, 10: 509-514. 10.1070/PU1968v010n04ABEH003699.

    Article  ADS  Google Scholar 

  2. Pendry JB, Holden AJ, Stewart WJ, Youngs I: Extremely low frequency plasmons in metallic mesostructures. Phys Rev Lett. 1996, 76: 4773-4776. 10.1103/PhysRevLett.76.4773.

    Article  ADS  Google Scholar 

  3. Pendry JB, Holden AJ, Robbins DJ, Stewart WJ: Magnetism from conductors and enhanced nonlinear phenomena. IEEE Trans Mircowave Theory Tech. 1999, 47: 2075-10.1109/22.798002.

    Article  ADS  Google Scholar 

  4. Smith DR, Padilla WJ, Vier DC, Nemat-Nasser SC, Schultz S: Composite medium with simultaneously negative permeability and permittivity. Phys Rev Lett. 2000, 84: 4184-4187. 10.1103/PhysRevLett.84.4184.

    Article  ADS  Google Scholar 

  5. Shelby RA, Smith DR, Schultz S: Experimental verification of a negative index of refraction. Science. 2001, 292: 77-79. 10.1126/science.1058847.

    Article  ADS  Google Scholar 

  6. Liu RP, Degiron A, Mock JJ, Simth DR: Negative index material composed of electric and magnetic resonators. App Phys Lett. 2007, 90: 263504-10.1063/1.2752120.

    Article  ADS  Google Scholar 

  7. Grbic A, Eleftheriades GV: Overcoming the diffraction limit with a planar left-handed transmission-line lens. Phys Rev Lett. 2004, 92: 117403-10.1103/PhysRevLett.92.117403.

    Article  ADS  Google Scholar 

  8. Grbic A, Eleftheriades GV: Growing evanescent waves in negative-refractive-index transmission-line media. Appl Phys Lett. 2003, 82: 1815-1817. 10.1063/1.1561167.

    Article  ADS  Google Scholar 

  9. Grbic A, Eleftheriades GV: Subwavelength focusing using a negative-refractive-index transmission line lens. IEEE Antennas Wireless Propagat Lett. 2003, 2: 186-189. 10.1109/LAWP.2003.819042.

    ADS  Google Scholar 

  10. Grbic A, Eleftheriades GV: Periodic analysis of a 2-D negative refractive index transmission line structure. IEEE Trans Antennas Propagat. 2003, 72: 2064-2072. []

    Google Scholar 

  11. Eleftheriades G, Iyer AK, Kremer PC: Planar negative refractive index media using periodically L-C loaded transmission lines. IEEE Trans Microwave Theory Tech. 2002, 50: 2702-2712. 10.1109/TMTT.2002.805197.

    Article  ADS  Google Scholar 

  12. Cui TJ, Cheng Q, Huang ZZ, Feng YJ: Electromagnetic wave localization using a left-handed transmission-line superlens. Phys Rev B. 2005, 72: 035112-10.1103/PhysRevB.72.035112.

    Article  ADS  Google Scholar 

  13. Qin Y, Cui TJ: A general representation of left-handed materials using LC-loaded transmission lines. Microwave Opt Technol Lett. 2006, 48: 2167-2171. 10.1002/mop.21883.

    Article  Google Scholar 

  14. Yao YH, Cui TJ, Cheng Q, Liu RP, Huang D, Smith DR: The realiazation of super waveguide using left-handed transmission-line circuits. Phys Rev E. 2007, 76: 036602-10.1103/PhysRevE.76.036602.

    Article  ADS  Google Scholar 

  15. Caloz C, Itoh T: Positive/negative refractive index anisotropic 2-D metamaterials. IEEE Microwave Wireless Compon Lett. 2003, 13: 547-10.1109/LMWC.2003.820641.

    Article  Google Scholar 

  16. Eleftheriades GV, Siddiqui OF: Negative Refraction and Focusing in Hyperbolic Transmission-Line Periodic Grids. IEEE Trans Microwave Theory Tech. 2005, 53: 396-10.1109/TMTT.2004.839944.

    Article  ADS  Google Scholar 

  17. Wong JKH, Balmain KG, Eleftheriades GV: Fields in planar anisotropic transmission-line metamaterials. IEEE Trans Antenna Propag. 2006, 54: 2742-10.1109/TAP.2006.882158.

    Article  ADS  Google Scholar 

  18. Kim JH, Swaminathan M: Modeling of irregular shaped power distribution planes using transmission matrix method. IEEE Trans Advanced Packag. 2001, 24: 334-346. 10.1109/6040.938301.

    Article  Google Scholar 

  19. Caloz C, Itoh T: Electromagnetic Metamaterials: Transmission Line Theory and Microwave Applications. 2006, New York: Wiley, []

    Google Scholar 

  20. Cheng Q, Cui TJ: High-power generation and transmission through a left-handed material. Phys Rev B. 2005, 2: 113112-10.1103/PhysRevB.72.113112.

    Article  ADS  Google Scholar 

Download references


This work was supported in part by the National Basic Research Program (973) of China under Grant No. 2004CB719802, in part by the National Science Foundation of China under Grant Nos. 60671015, 60496317, 60225001, and 60621002, and in part by the National Doctoral Foundation of China under Grant No. 20040286010.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Tie Jun Cui.

Authors’ original submitted files for images

Rights and permissions

Open Access This is an open access article distributed under the terms of the Creative Commons Attribution Noncommercial License ( ), which permits any noncommercial use, distribution, and reproduction in any medium, provided the original author(s) and source are credited.

Reprints and permissions

About this article

Cite this article

Ma, H.F., Cui, T.J., Chin, J.Y. et al. Fast and accurate simulations of transmission-line metamaterials using transmission-matrix method. PMC Phys B 1, 10 (2008).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: