Numerical determination of frequency behavior in cloaking structures based on LC distributed networks with TLM method

The increasing interest in metamaterials with negative refractive index has been prompted by a variety of promising optical and microwave applications. Often, the resulting electromagnetic problems to be solve are not analytically derivable; therefore, numerical modeling must be employed and the Transmission Line Modeling (TLM) method constitutes a possible choice. After having greatly simplified the existing TLM techniques for the modeling of metamaterials, we propose in this paper to carry out a frequency study of cloaking structure. ©2008 Optical Society of America OCIS codes: (230.3205) Invisibility cloaks; (050.1755) Computational electromagnetic methods; (160.3918) Metamaterials. References and links 1. P. B. Johns and R. L. Beurle, “Numerical solution of 2-dimensional scattering problems using a transmission-line matrix,” Proc. Inst. Elec. Eng. 118, 1203-1208 (1971). 2. C. Christopoulos, The Transmission-Line Modeling Method, The Institute of Electrical and Electronic Engineers (New York and Oxford University Press, Oxford, 1995). 3. C. Blanchard, J. A. Portí, J. A. Morente, A. Salinas, and E. A. Navarro, “Determination of the effective permittivity of dielectric mixtures with the transmission line matrix method,” J. Appl. Phys. 102, 064101 (2007). 4. P. P. M. So, H. Du, and W. J. R. Hoefer, “Modeling of metamaterials with negative refractive index using 2-D shunt and 3-D SCN TLM networks,” IEEE Trans. Microwave Theory Tech. 53, 1496-1505 (2005). 5. C. Blanchard, J. Portí, B-I Wu, J. A. Morente, A. Salinas, and J. A. Kong, “Time domain simulation of electromagnetic cloaking structure with TLM method,” Opt. Express 16, 6461-6470 (2008). 6. J. B. Pendry, D. Schurig, and D. R. Smith, “Controlling electromagnetic fields,” Science 312, 1780-1782 (2006). 7. U. Leonhardt, “Optical conformal mapping,” Science 312, 1777-1780 (2006). 8. G. V. Eleftheriades, A. K. Iyer, and P. C. Kremer, “Planar negative refractive index media using periodically L-C loaded transmission lines,” IEEE Trans. Microwave Theory Tech. 50, 2702-2712 (2002). 9. Y. Huang, Y. Feng, and T. Jiang, “Electromagnetic cloaking by layered structure of homogeneous isotropic materials,” Opt. Express 15, 11133 (2007). 10. J. A. Portí, J. A. Morente, A. Salinas, M. Rodríguez-Sola, and C. Blanchard, “On the circuit description of TLM nodes,” Int. J. Electron. 93, 479-491 (2006). 11. J. A. Portí, J. A. Morente, and M. C. Carrión, “Simple derivation of scattering matrix for TLM nodes,” Electron. Lett. 34, 1763-1764 (1998). 12. P. B. Johns, “A symmetrical condensed node for the TLM method,” IEEE Trans. Microwave Theory Tech. 35, 370-377 (1987). 13. Z. Liang, P. Yao, X. Sun, and X. Jiang, “The physical picture and the essential elements of the dynamical process for dispersive cloaking structures,” Appl. Phys. Lett. 92, 131118 (2008).


Introduction
TLM is a well-known time-domain numerical method which can be regarded as a pure modeling procedure rather than a Maxwell's equations solver.Indeed, TLM consists of modeling a certain medium and its electromagnetic (EM) properties by filling the field space with a network of transmission lines formed by interconnecting unitary circuits or cells, termed nodes [1].Any variation of the EM parameters, such as permittivity, permeability, or conductivity is achieved by adding extra stub transmission lines to the nodes [2].This characteristic confers to TLM an interesting versatility, which allows simple and elegant modeling of non-homogeneous media [3].
The increasing interest in the study of metamaterials (MM) has made necessary to adapt TLM to them.So et al. have presented a novel technique for modeling two-(2D) and threedimensional (3D) MM with TLM [4].However, as it has been pointed out by the authors themselves, the new node presented in this innovative paper requires a significant modification of the impulse scattering process when compared to the classical simple version.Considering that these artificial media are on the cutting-edge of current research, it is of first interest to render more comfortable the modeling involving MM.In this sense, it has been shown recently that such a drastic modification is not necessary for the 2D case; a simple 2D condensed node with only three new stubs is perfectly able to model materials with negative parameters [5].This new technique offers the possibility of easily modeling complex structures involving MM whilst maintaining the potentiality of a time domain method.After a complete presentation of the new nodes (TE and TM modes), we propose in this article to model cloaking structures in order to determine their behavior in terms of the frequency [6,7].It is worth noting that the TLM technique presented here is nothing more than the numerical incarnation of the well-known dual transmission line approach of MM [8]; this renders the frequency study particularly interesting since the TLM mesh scheme is exactly the same as the real dual network that could be employed to construct an effective cloaking structure.
First, we propose in this paper a solid description of the series node for MM associated with a TM mode, the scattering matrix is in particular entirely given.Second, the description, already given in [5], of the parallel node associated with a TE mode is completed.Third, by using the approximation presented by Huang et al. [9], a cloaked Perfectly Electric Conducting (PEC) cylinder is modeled, while the time-domain nature of TLM is exploited so as to study, with a single simulation run, its frequency behavior over a wide spectral range.Although the far field pattern for a simple PEC cylinder is different depending on the TE or TM mode nature of the illuminating EM wave, the results are shown to be the same for the two polarizations once the cylinder is concealed in the cloaking shell.

Series node for TM modes
For propagation on the x-y plane, the only non-zero field components for a TM mode are H z , E x , and E y .This polarization may be simulated by a TLM series node formed with seven transmission lines [10].The link lines 1 to 4 are of identical impedance Z 0 =1/Y 0 .Line 5, with impedance Z z Z 0 , adds extra inductance (shown to be L z = Z z Z 0 Δt/2) to the node, allowing an independent control of the relative permeability μ z ; while lines 6 and 7, with admittance Y x Y 0 and Y y Y 0 respectively, are stubs which add extra capacitance (shown to be C x = Y x Y 0 Δt/2 and C y = Y y Y 0 Δt/2 respectively) to the node, allowing an independent control of the relative permittivity ε x and ε y respectively, Δt being the TLM time-step.
As shown in Fig. 1, this node can be split into three parts: a series circuit made up of five transmission lines describing the H z component, and two parallel circuits, each one constituted by three transmission lines, for the E x and E y components.Each sub-circuit is described by an equation providing Y x , Y y , and Z z : Usually, the impedance Z 0 is chosen to make Z z of Eq. ( 1) equal to zero in free space (with impedance η 0 ) in the case Δx = Δy = Δz.This yields Z 0 = η 0 /√2.
The scattering matrix S is obtained by following the procedure described in [11]: Concerning the calculation of the EM field at the center of a stubbed SCN, its three components may be obtained by using Thevenin equivalents of the circuits shown in Fig. 1 [2],   ( ) A particular field component at a node is excited by injecting voltage into the appropriate lines.The hereunder set of voltages excites the desired EM field: 0.

Metamaterials modeling with a series node
The above development allows the modeling of usual material.Adapting it to MM is very natural.Instead of equipping the series sub-circuit of Fig. 1(a) with an inductive stub, we provide it with a capacitive stub.As a result, the introduced capacitance into the node, In the same way, the capacitive stubs of the parallel sub-circuits of Fig. 1 are substituted by inductive stubs, similar expressions applying for these other two lines.Finally, Y x , Y y , and Z z of Eq. ( 1) are now given by: This set of equations allows modeling, for a certain frequency, relative ε and μ below unity and even the zero value, which becomes a natural value with this approach.It is worth noting that the impedance and admittances given by Eq. ( 1) diverge only by the factor -Δt 2 ω 2 /4 respectively to those of Eq. ( 6).The scattering matrix, as well as Eqs.(3)(4)(5), is absolutely unaltered with respect to the classical node, which renders the technique very comfortable to use.For propagation on the x-y plane, the only non-zero field components for a TE mode are E z , H x , and H y .The parallel node is required for such a polarization, and it is obtained as it has been done hereinbefore for the series node.The three circuits describing each component of the EM field are depicted in Fig. 2. Similarly to Eq. ( 1), the corresponding equations giving Y z , Z x , and Z y are:

Parallel node for TE modes
In Eq. ( 7), Z 0 is usually chosen to be Z 0 = η 0 √2, so that there is no stub if free space is modeled with Δx = Δy = Δz.
The scattering matrix, as well as its elements, may in this case be written as The three components of the EM fields at the center of each node are while the following set of incident pulses allows exciting the desired EM source: 0.5 Finally, the procedure adopted in Section 2.2 may be employed for the modeling of MM with the parallel node.This lets Eqs.(8)(9)(10)(11) absolutely unchanged, while Eq. ( 7) becomes:

Modeling of metamaterials with the 3D Symmetrical Condensed Node
Finally, we would like to make a comment on the 3D modeling of MM.The well-known Symmetrical Condensed Node (SCN) is usually employed in TLM for 3D simulation because it offers very high performance [12].In its original version, the SCN has six extra lines (three capacitive open circuits for the E-field, and three inductive short circuits for the H-field), in addition to the regular 12 link lines.So et al. have claimed that the corresponding scattering matrix for modeling metamaterials with the 3D SCN would be of size 27×27, which is unnecessarily too big to be easily implemented, and pushes the authors to introduce an intercell approach to get around the problem [4].Actually, in the light of the development presented hereinbefore, it appears that the original SCN described by the 18×18 scattering matrix is perfectly well suited for the implementation of MM; there is in particular no need to add neither extra lines nor changes in the scattering matrix or inter-cell positions.The only requirement is to substitute the inductive stubs by capacitive stubs, and vice versa.

Numerical results
Employing a time domain method like TLM for the simulation of a cloaked object presents a certain advantages.For instance, the actual dynamic process can be modeled [13], or the frequency behavior can be reached from a simple Fourier Transform.A cloaked PEC infinite cylinder is modeled with the series node for TM mode, following the same approach presented by Huang et al. [9].This approach consists basically in substituting the anisotropic medium of the shell by a concentric layered structure of alternating homogeneous isotropic materials.Furthermore, the permittivity and permeability of the layers are calculated very close to their inner boundary (but not rigorously in the vicinity of the boundary in order to avoid divergence of the EM parameters in the first layer) since this procedure has been shown to give the best results [5].The cloaking entity is illuminated by a Gaussian pulse, which theoretically contains all the frequencies.Nevertheless, 2 GHz is chosen to be the working frequency in order to fix the parameters of Eq. ( 6), which are frequency dependent.The inner and the outer radius of the cloaking shell are R 1 = 0.1 m and R 2 = 0.2 m, respectively, and it is made up of 20 layers.The far field pattern for 2 GHz is depicted in Fig. 3; moreover, by using a simple Fourier Transform, the far field in terms of the frequency for five different angles (from 0º to 180º using 45º steps) is obtained and shown in the same Fig. 3.As expected, the cloaking shell is manifestly efficient for the 2 GHz functional frequency, but it is worth noting that a frequency band appears around this frequency for all directions with the noticeable exception of 0º.This last direction is characterized by a very narrow low radiation region, which shows that it is the most conflictive direction.On the other hand, simulating the TE mode requires the use of the parallel node presented in section 2.3.It turns out that the results for this polarization are exactly the same as for the TM mode.This means that, although the cloaking is not perfect, the PEC cylinder is well concealed by the cloaking shell.Indeed, it is well known, and moreover verified by our TLM simulations, that the Scattering Width of a simple circular conducting cylinder depends strongly on the polarization.In order to get a decisive confirmation of this observation, the same cloaking structure, but without the cylinder inside, is modeled, and once again, the same results are obtained.To explain this, let us note that the first layer of the cloaking shell in our simulation has a relative permittivity ε = 2.10 5 , and relative permeability μ = 4.10 -5 .It is reasonable to think that these extreme values shield the central region.Nevertheless, the other layers are not perfectly able to steer the radiation around themselves due to both numerical discretization and to the approximation consisting on substituting the theoretical anisotropic material by isotropic layers, which produces the observable forward scattering.
Finally, it is worth noting that the technique presented above corresponds to an actual situation for the whole frequency range.Indeed, it is well known that a distributed system with series capacitors and shunt inductors is a possible way to build materials with permittivity and permeability below unity.Therefore, the TLM nodes presented in this paper for modeling MM, are completely equivalent to such a dual L-C network, and are thus expected to perfectly describe it, not only at the design frequency, but also for the whole frequency band.

Conclusion
TLM nodes for MM modeling have been reported in this article; one for the 2D TM mode, another one for the TE mode.We have proposed a detailed description of each one; the corresponding scattering matrices have been given.These new nodes have the merit to leave the original ones almost unaltered.
As an application, we have built a cloaking structure for a TM polarization.The advantages offered by a time domain method such as TLM have been applied to propose a behavior study in terms of the frequency in the cloaked PEC cylinder.In this sense, it is worth noting that modeling a cloaking structure with TLM is more than simply solving Maxwell's equations, it is indeed equivalent to substituting the cloak shell by a dual distributed L-C network, which is well known as an actual technique for implementing the exotic material parameters associated with MM.

Fig. 1 .
Fig. 1.Splitting of the TLM series node into 3 sub-circuits: (a) series sub-circuit for H z , (b) parallel sub-circuits for E y , and (c) parallel sub-circuits for E x .

Fig. 2 .
Fig. 2. Splitting of the TLM parallel node into 3 sub-circuits: (a) parallel sub-circuit for E z , (b) series sub-circuits for H y , and (c) series sub-circuits for H x .