Research Article | Open Access

# FDTD-SPICE for Characterizing Metamaterials Integrated with Electronic Circuits

**Academic Editor:**Alejandro Lucas Borja

#### Abstract

A powerful time-domain FDTD-SPICE simulator is implemented and applied to the broadband analysis of metamaterials integrated with active and tunable circuit elements. First, the FDTD-SPICE modeling theory is studied and details of interprocess communication and hybridization of the two techniques are discussed. To verify the model, some simple cases are simulated with results in both time domain and frequency domain. Then, simulation of a metamaterial structure constructed from periodic resonant loops integrated with lumped capacitor elements is studied, which demonstrates tuning resonance frequency of medium by changing the capacitance of the integrated elements. To increase the bandwidth of the metamaterial, non-Foster transistor configurations are integrated with the loops and FDTD-SPICE is applied to successfully bridge the physics of electromagnetic and circuit topologies and to model the whole composite structure. Our model is also applied to the design and simulation of a metasurface integrated with nonlinear varactors featuring tunable reflection phase characteristic.

#### 1. Introduction

In recent years, the subject of metamaterials integrated with electronic circuits, nonlinear elements, and non-Foster transistor configurations have been drawing significant attention in research communities. Metamaterials are typically periodic structures created by array of resonant elements which are small and stationary compared to the wavelength of electromagnetic wave. The combination of metamaterials with circuit elements can bring novel and tunable characteristics to these materials. For instance combining loop-based metamaterials with varactors provides a means to electronically tune the performance by applying a DC voltage or a high-power signal without modifying the structure [1–4]. Alternatively, a non-Foster transistor configuration can be introduced to enhance the metamaterial bandwidth [5]. Active components can be also integrated with transmission line metamaterials to compensate the loss [6], and finally metasurface can be integrated with circuit elements to tailor their reflection phase performance [7].

The complexity of the structure requires a powerful computational technique to bridge the physics of electromagnetic metamaterial with circuit configuration. In this paper, the Finite Difference Time Domain (FDTD) technique is combined with SPICE-like software to provide a capable numerical technique for time-domain full wave analysis of the composite structure. The technique is applied to the design and simulation of metamaterials integrated with circuit elements. Broadband performance is obtained with the focus on bandwidth and tunable behavior engineering.

Section 2 provides the details of FDTD-SPICE hybridization and shows the validation through some examples. In Section 3, the following three applications are explored: (a) loop-based metamaterial integrated with lumped capacitors to tune the resonance performance, (b) metamaterial integrated with non-Foster transistor configuration for bandwidth manipulation, and (c) a metasurface combined with varactors featuring tunable reflection phase.

#### 2. FDTD-SPICE Computational Model

##### 2.1. Hybridization of FDTD with SPICE

As a major computational electromagnetic modeling scheme, Finite Difference Time Domain technique can tackle problems by providing a full wave solution in time. It is a capable method for modeling a variety of scientific and engineering problems [8, 9]. On the other hand, the general-purpose analog electronic circuit simulator, SPICE [10], can solve complex circuit topologies in time domain. It makes great sense to combine FDTD with SPICE to feature a powerful numerical engine for solving metamaterials integrated with circuit configurations.

FDTD uses Maxwell’s equations (1) and (2) in time domain to recursively calculate the electromagnetic fields at each position in three-dimensional space: While both equations can provide the coupling between FDTD and SPICE according to reference [8], the first equation is used in this paper. After some simple manipulations, Ampere’s Law of (1) can be written as [11] Equation (3) describes a circuit model as shown in Figure 1, where is a capacitance that abstracted from the model (note that here is not the intrinsic lattice capacitance of FDTD cell), is the cell size, and is the area of one FDTD unit cell. is the current across the SPICE element, and is the voltage across the element, is the whole current of the current source. Thus, for the cell containing SPICE element, instead of using (1) to update the electric field, SPICE can carry out (3) and obtain the voltage across the element and then pass it to FDTD program. From the voltage, it is easy to obtain field afterward.

To determine the voltage across the corresponding element, SPICE software needs to run a transient simulation for a period of , which should be the same as FDTD time step. Then, in the next time step, regular FDTD method will update all the fields except the cell with SPICE element. For the SPICE element, SPICE program needs to run another , starting from where it stopped last time. Before each run of SPICE, FDTD sends data that contains circuit information to SPICE. After each time step of SPICE simulation, it passes results back to FDTD. In this way, the FDTD and SPICE engines run alternatively until the end of the whole simulation time. With the help of SPICE, FDTD is ready to simulate any complex integrated configuration as long as it is provided in SPICE’s library. Since conventional FDTD can only model very simple circuit elements, SPICE function adds a big advantage to FDTD.

The SPICE software used here is SPICE3f5. It has an interactive mode that can receive commands from and send results back to the interface. It also has “stop” and “resume” commands that can set a breakpoint after running some time steps in a transient simulation and then resume it. These features make it very suitable to be connected to FDTD.

To exchange data between FDTD and SPICE, interprocess communication should be developed. Figure 2 presents a complete picture of our FDTD-SPICE model. When a simulation begins, FDTD main program calls fork function which spawns two processes: the son process that will run SPICE, and the parent process which is the original FDTD process itself. Two sockets: Socket[0] and Socket[1], are created in the main program before it forks so that both of the two children have access to the sockets. Here these sockets act like two-way pipe, where FDTD process sends SPICE commands to Socket[0] and reads results from Socket[1], while SPICE process reads message from Socket[0] and sends results back to Socket[1]. FDTD process will be waiting for SPICE outputs while SPICE process is running and vice versa. These two processes then communicate through the two sockets at each time step until the end of the simulation.

While DC solution is needed in SPICE for many cases but not very trivial to be obtained with FDTD, literature [12, 13] addressed DC solution problem by separating the AC and DC components, which is helpful for the metasurface simulations of diode varactors with bias that presented in Section 3.

Our FDTD code has been written in a very powerful manner and it has many features such as modeling periodic structures and linear dispersive metamaterials [14]. By applying auxiliary differential equation (ADE) method [15], it can also model nonlinear dispersive materials as well as gain materials. The absorbing wall is based on CPML [16] where one can truncate the computational domain very close to the structure. The code has been applied successfully to a variety of metamaterials composites [9, 17–19].

##### 2.2. Test Results

Two cases are studied to validate the FDTD-SPICE model. All simulations in this paper use uniform gridding and are run on Sun Ultra45 workstation with 2 GHz CPU and 4 GB memory.

In the first case as depicted in Figure 3(a), one parallel RLC circuit is connected to a Gaussian pulse voltage source in order to obtain the impedance. There are about 8000 time steps in the calculation with each time step , cell size mm. Including CPML boundary, the whole computational domain is modeled by Yee’s cells. The FDTD process takes 456 seconds, while, since the circuit is simple, SPICE process takes only 10 seconds. Figure 3(b) shows both the FDTD-SPICE results and the analytical calculation for the real and imaginary parts of the impedance , which agree well with each other.

**(a)**

**(b)**

In the second case as shown in Figure 4(a), one BFG403W transistor from Discrete Semiconductor is utilized in a single frequency simulation with a sinusoidal voltage source. The BJT has a saturation current aA and the forward active current gain is . The calculation has 18000 time steps with each time step fs and cell size *μ*m. Because of the complex property of BJT circuit, the cell size and time step are much smaller than those used in previous simulation. The whole computational domain has a size of Yee’s cells. FDTD and SPICE processes take 1734 seconds and 41 seconds, respectively. Figure 4(b) shows the comparison for the waveform of voltage across resistor R2 in time domain, between FDTD-SPICE and complete SPICE simulation. An excellent agreement is demonstrated.

**(a)**

**(b)**

#### 3. Applications on Metamaterials

This section presents applications of FDTD-SPICE method to some recent metamaterial designs, namely, loop-based composite metamaterial with capacitors, metamaterial integrated with non-Foster circuits, and nonlinear elements loaded metasurface. In all of the cases, uniform gridding, perfectly matched layers, and periodic boundary conditions are used. Since they all deal with plan wave incidence, total field/scattered field method is also utilized [17].

##### 3.1. Metamaterial Loops Integrated with Passive Elements

A periodic structure made from metallic loops can artificially create a medium with permeability characteristics (positive and negative) [19]. The metamaterial structure to be simulated here is shown in Figure 5(a). It is a four-layer structure in direction where each layer has loops periodic in - directions and the loops are terminated to capacitor . The unit loop dimension is 4.8 mm and 1.8 mm in and directions, respectively. The excitation is a plane wave that has Ez/Hy components and propagates in direction with a Gaussian signature. A transmission line analogy can be obtained for plane wave propagation through the medium [19], which is embedded with a loop that has a self-inductance and is terminated to a lumped capacitor . By changing the circuit elements ( and/or ) in the equivalent transmission line model, the effective parameter of the medium can be tuned, as well as the property of the metamaterial.

**(a)**

**(b)**

FDTD-SPICE method is applied to model the structure and obtain the transmission coefficient as shown in Figure 5(b). The loops are modeled in FDTD and the capacitors in SPICE. There are 15700 time steps in the calculation with each time step ps, cell size mm. The whole computational domain has a size of Yee’s cells. FDTD and SPICE processes take 1669 seconds and 1319 seconds, respectively. Because of the multiple-port interprocess communication, SPICE process needs considerable simulation time, which may be improved by introducing batch data communication between processes. Figure 5(b) shows 4 simulation results. The red dots represent results from pure FDTD simulation with pF ( is also modeled in FDTD). Then the blue line comes from FDTD-SPICE simulation with the same parameter, and it matches well with red dots, which show the accuracy of our model. Changing the capacitor allows tuning the resonance and as well as the equivalent permeability of the metamaterial. Figure 5(b) also shows the transmission coefficients for the metamaterial loops terminated to pF and pF lumped capacitors, depicted in black and magenta, respectively. Doubling capacitance would half the resonance frequency. As can be seen from Figure 5(b), the resonance frequency is shifted from 1.1 GHz down to 0.78 GHz and 0.55 GHz, respectively, according the value of capacitors.

##### 3.2. Integration with Non-Foster Circuit Elements

Non-Foster elements have recently attracted significant interests due to their novel properties in offering negative impedance performance and successfully engineering the frequency dispersion characteristics of antennas, transmission lines, and metamaterials towards the goals of interests [20–23]. A non-Foster element can be realized using a Negative Impedance Converter (NIC) circuit by translating a given positive load into a negative one. NIC circuit is made from transistor elements and one would need a comprehensive modeling scheme to characterize them while they are integrated with metamaterials. It is also worth mentioning that the design of NIC circuit involves several challenging issues such as noise, biasing, and stability and special care must be taken to obtain a wideband performance for non-Foster circuit. The NIC stability (which is not the subject of this paper) has also been checked with the help of pole-zero analysis to ensure that the system has absolutely no pole in the Right-Half-Plane (RHP) of the complex S-plane.

Here we explore the integration of NIC circuit with metamaterial loops [5] and its FDTD-SPICE simulation to investigate the bandwidth enhancement of the structure. A two-layer loop-based composite metamaterial with capacitors as shown in Figure 5(a) is utilized. The loops are terminated to capacitors and inductors (4 nH). By integrating NIC designs offering −2 nH inductance, we reduce the effect of the positive inductance and enhance the bandwidth. The NIC circuit used in this paper is based on Linvill’s model [24], and it is an example of Voltage Inversion NIC (VINIC). As shown in Figure 6(a), the circuit consists of two N-type BJTs and a lumped passive inductor nH as the load. After adding the non-Foster circuit to the loop, the equivalent circuit model [19] for the transmission line is depicted in Figure 6(a), where the NIC configuration is included in a negative inductor .

**(a)**

**(b)**

The simulation has 11600 time steps with each time step ps, cell size mm. The whole computational domain has a size of Yee’s cells. FDTD and SPICE processes take 686 seconds and 133 seconds, respectively. Because of the complexity of the NIC circuit model, SPICE needs more time than the test examples from previous section. Here three cases are performed. Firstly, metamaterial with loop terminated to capacitor pF and lumped inductor nH is simulated. As depicted by the blue curve in Figure 6(b), it has a resonance around 2.4 GHz. Secondly, to broaden the bandwidth, a NIC circuit nH is added to inductor nH, making the whole inductance 2 nH. In order to make a fair comparison, the capacitor is made pF so that the two simulations have the same resonance. As can be observed from the red curve in Figure 6(b), the metamaterial integrated with non-Foster element has a wider bandwidth than the case without it (blue curve). Increasing the value of can enable a wider bandwidth performance (as has been comprehensively detailed in [5]). Lastly, to validate the results of the second case, NIC circuit nH together with inductor nH is replaced by one lumped inductor nH. As shown in Figure 6(b), the black curve follows the behavior of the second case.

At the end we want to mention that the reason we have considered a positive 4 nH inductor in series with the −2 nH non-Foster element is to make sure the time-domain SPICE modeling is stable (the equivalent total inductance is positive, nH).

##### 3.3. Metasurface Integrated with Varactor Diodes

Metasurfaces have great potential for making low profile antennas [25], and achieving an electronically tunable reflection phase is of remarkable interest in this regard. This section presents the simulation for a structure of metasurface integrated with varactors. The metasurface used in this paper is similar to the one in [7], which is designed to act as an artificial magnetic conductor (AMC) and uses active NIC circuits as negative elements to increase the bandwidth. Here varactor diodes are integrated in the metamaterial and make the metasurface tunable in resonance frequency.

Figure 7(a) shows a metasurface structure consisting of periodic square metal patches and PEC ground plane with dielectric slab in between. It also presents the unit cell of the metasurface loaded with varactor diodes. The structure is periodic in and directions and the incidence plane wave is direction. The FDTD-SPICE method is applied to simulate the structure. There are 7000 time steps in the calculation with each time step ps, cell size mm. The whole computational domain has a size of Yee’s cells. FDTD and SPICE processes take 386 seconds and 8 seconds, respectively. Figure 7(b) shows the reflection coefficient phase for the metasurface without varactor diodes, as well as those with diodes that have different biases. The green line, as a reference frequency response, represents the one without varactor diode, while the others are for results that have diode with different bias. The varactor diode model used here is Skywork’s SMV1232-079 Hyperabrupt Tuning Varactors [26]. When the diode is added to the structure, it increases the total capacitance therefore decreases the resonance frequency of the metasurface. As can been seen from Figure 8, typical capacitance values of the varactor diode decrease when the reverse bias voltage increases. The capacitance is a maximum while there is no reverse bias voltage. Thus for reverse bias with smaller voltage it should give larger shift in resonance frequency, which agrees well with the simulation results.

**(a)**

**(b)**

#### 4. Conclusion

This paper develops a hybridized FDTD-SPICE model and applies it to the simulations of some novel metamaterials integrated with electronic circuits. The Maxwell’s equations are linked with SPICE through the interchange between electromagnetic fields and voltages across the SPICE element nodes. Interprocess communication technique with socket is used as it is crucial for FDTD and SPICE to exchange data at proper running time. Diagrammatic flowchart is discussed that demonstrates the detail of the execution process as well as the data exchanging with the help of sockets. To show the validity of the model, two simple cases that have analytical solution or pure SPICE results are investigated. Then, three kinds of metamaterial are studied and modeled by the FDTD-SPICE engine. First, loop-based metamaterial with lumped capacitor is investigated where changing the capacitance can vary the resonance frequency. Then, integration of metamaterial with non-Foster circuit element, which provides negative impedance and increases the bandwidth, is studied and characterized. Lastly, metasurface simulation with integrated nonlinear varactor elements shows a tunable reflection phase property for this metamaterial. The technique is very powerful and its accuracy is validated well. We expect our scheme to be of great engine for comprehensive simulation of metamaterials integrated with active and nonlinear circuits and allow exploring new artificial materials combined with a lot of potential complex circuits.

#### Acknowledgments

This work was supported in part by the US Office of Naval Research (ONR) Grant no. N00014-10-1-0264. The authors would like to thank Allen Taflove from Northwestern University for his kind discussion on the FDTD-SPICE model theory and acknowledge Vincent Thomas from Los Alamos National Laboratory for his helpful information on SPICE and sockets programming.

#### References

- D. A. Powell, I. V. Shadrivov, Y. S. Kivshar, and M. V. Gorkunov, “Self-tuning mechanisms of nonlinear split-ring resonators,”
*Applied Physics Letters*, vol. 91, no. 14, Article ID 144107, 3 pages, 2007. View at: Publisher Site | Google Scholar - I. V. Shadrivov, A. B. Kozyrev, D. W. Van Der Weide, and Y. S. Kivshar, “Tunable transmission and harmonic generation in nonlinear metamaterials,”
*Applied Physics Letters*, vol. 93, no. 16, Article ID 161903, 3 pages, 2008. View at: Publisher Site | Google Scholar - D. Huang, E. Poutrina, and D. R. Smith, “Analysis of the power dependent tuning of a varactor-loaded metamaterial at microwave frequencies,”
*Applied Physics Letters*, vol. 96, no. 10, Article ID 104104, 3 pages, 2010. View at: Publisher Site | Google Scholar - I. V. Shadrivov, S. K. Morrison, and Y. S. Kivshar, “Tunable split-ring resonators for nonlinear negative-index metamaterials,”
*Optics Express*, vol. 14, no. 20, pp. 9344–9349, 2006. View at: Publisher Site | Google Scholar - S. Saadat, M. Adnan, H. Mosallaei, and E. Afshari, “Composite Metamaterial and Metasurface Integrated with Non-Foster Active Circuit Elements: A Bandwidth-Enhancement Investigation,”
*IEEE Transactions on Antennas and Propagation*. In press. View at: Google Scholar - T. Jiang, K. Chang, L.-M. Si, L. Ran, and H. Xin, “Active microwave negative-index metamaterial transmission line with gain,”
*Physical Review Letters*, vol. 107, no. 20, Article ID 205503, 5 pages, 2011. View at: Publisher Site | Google Scholar - D. J. Kern, D. H. Werner, and M. J. Wilhelm, “Active negative impedance loaded EBG structures for the realization of ultra-wideband artificial magnetic conductors,” in
*Proceedings of the IEEE International Antennas and Propagation Symposium and USNC/CNC/URSI North American Radio Science Meeting*, vol. 2, pp. 427–430, June 2003. View at: Google Scholar - A. Taflove and S. C. Hagness,
*Computational Electrodynamics: The Finite-Difference Time-Domain Method*, House, Norwood, Mass, USA, 3rd edition, 2005. - H. Mosallaei and Y. Rahmat-Samii, “Broadband characterization of complex periodic EBG structures: an FDTD/Prony technique based on the split-field approach,”
*Electromagnetics*, vol. 23, no. 2, pp. 135–151, 2003. View at: Publisher Site | Google Scholar - http://bwrc.eecs.berkeley.edu/classes/icbook/spice/.
- V. A. Thomas, M. E. Jones, M. Piket-May, A. Taflove, and E. Harrigan, “Use of SPICE lumped circuits as sub-grid models for FDTD analysis,”
*IEEE Microwave and Guided Wave Letters*, vol. 4, no. 5, pp. 141–143, 1994. View at: Publisher Site | Google Scholar - G. Kobidze, A. Nishizawa, and S. Tanabe, “Ground bouncing in PCB with integrated circuits,” in
*Proceedings of the IEEE International Symposium on Electromagneti Compatibility*, vol. 1, pp. 349–352, August 2000. View at: Google Scholar - N. Orhanovic, R. Raghuram, and N. Matsui, “Full wave analysis of planar interconnect structures using FDTD-SPICE,” in
*Proceedings of the 51st Electronic Components and Technology Conference*, pp. 489–494, June 2001. View at: Google Scholar - H. Mosallaei, “FDTD-PLRC technique for modeling of anisotropic-dispersive media and metamaterial devices,”
*IEEE Transactions on Electromagnetic Compatibility*, vol. 49, no. 3, pp. 649–660, 2007. View at: Publisher Site | Google Scholar - M. Fujii, M. Tahara, I. Sakagami, W. Freude, and P. Russer, “High-order FDTD and auxiliary differential equation formulation of optical pulse propagation in 2-D Kerr and Raman nonlinear dispersive media,”
*IEEE Journal of Quantum Electronics*, vol. 40, no. 2, pp. 175–182, 2004. View at: Publisher Site | Google Scholar - J. A. Roden and S. D. Gedney, “Convolutional PML (CPML): an efficient FDTD implementation of the CFS-PML for arbitrary media,”
*Microwave and Optical Technology Letters*, vol. 27, pp. 334–339, 2000. View at: Google Scholar - H. Mosallaei and Y. Rahmat-Samii, “Periodic bandgap and effective dielectric materials in electromagnetics: characterization and applications in nanocavities and waveguides,”
*IEEE Transactions on Antennas and Propagation*, vol. 51, no. 3, pp. 549–563, 2003. View at: Publisher Site | Google Scholar - H. Mosallaei and K. Sarabandi, “Magneto-dielectrics in electromagnetics: concept and applications,”
*IEEE Transactions on Antennas and Propagation*, vol. 52, no. 6, pp. 1558–1567, 2004. View at: Publisher Site | Google Scholar - H. Mosallaei and K. Sarabandi, “Design and modeling of patch antenna printed on magneto-dielectric embedded-circuit metasubstrate,”
*IEEE Transactions on Antennas and Propagation*, vol. 55, no. 1, pp. 45–52, 2007. View at: Publisher Site | Google Scholar - W. Jendernalik, “A low-voltage CMOS negative impedance converter for analogue filtering applications,”
*Bulletin of the Polish Academy of Sciences*, vol. 55, no. 4, pp. 419–423, 2007. View at: Google Scholar - R. M. Rudish and S. E. Sussman-Fort, “Progress in use of non-Foster impedances to. match electrically-small antennas and arrays,” in
*Proceedings of the Antenna Applications Symposium*, pp. 89–108, Allerton Park, Ill, USA, September 2005. View at: Google Scholar - P. Jin and R. W. Ziolkowski, “Broadband, efficient, electrically small metamaterial-inspired antennas facilitated by active near-field resonant parasitic elements,”
*IEEE Transactions on Antennas and Propagation*, vol. 58, no. 2, pp. 318–327, 2010. View at: Publisher Site | Google Scholar - S. Hrabar, I. Krois, and A. Kiricenko, “Towards active dispersionless ENZ metamaterial for cloaking applications,”
*Metamaterials*, vol. 4, no. 2-3, pp. 89–97, 2010. View at: Publisher Site | Google Scholar - J. G. Linvill, “Transistor negative-impedance converters,”
*Proceedings of the IRE*, vol. 41, pp. 725–729, 1953. View at: Google Scholar - H. Mosallaei and K. Sarabandi, “Antenna miniaturization and bandwidth enhancement using a reactive impedance substrate,”
*IEEE Transactions on Antennas and Propagation*, vol. 52, no. 9, pp. 2403–2414, 2004. View at: Publisher Site | Google Scholar - Skywork’s SMV1232-079 datasheet.

#### Copyright

Copyright © 2012 Zhengwei Hao et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.