Left-handed materials: Transfer matrix method studies Peter Markos and C. M. Soukoulis
Outline of Talk What are Metamaterials? An Example: Left-handed Materials Results of the transfer matrix method Negative n and FDTD results New left-handed structures Applications/Closing Remarks
What is an Electromagnetic Metamaterial? A composite or structured material that exhibits properties not found in naturally occurring materials or compounds. Left-handed materials have electromagnetic properties that are distinct from any known material, and hence are examples of metamaterials.
Electromagnetic Metamaterials Example: Metamaterials based on repeated cells
Material Properties All EM phenomena can be explained by Maxwell s four Equations. When materials are present, we average over the detailed properties to arrive at two parameters: m and e. Quantities e and m are typically positive. When e and m are both negative, there are dramatically new properties.
Why create metamaterials? The electromagnetic response of naturally occurring materials is limited. Metamaterials can extend the material properties, while simultaneously providing other advantageous properties (e.g., strength, thermal, conformal) Possibility of materials with ideal electromagnetic response over broad frequency ranges, that can be customized, made active, or made modulable or tunable.
Electromagnetic waves interact with materials Lenses Gratings Modulators Switches Antennas Multiplexers Sources
Related Terms Artificial Dielectrics Effective Medium Theory: Connects the macroscopic electromagnetic properties of a composite to its constituent components, whether atoms, molecules, or larger scale sub-composites. Results in e and m. Photonic Band Gap Periodic structures, periodicity is on the order of the wavelength. Metamaterials
What s New? Advances in Computational Ability Advances in Fabrication Ability Market Needs
An Example: Left-handed Materials Predict structure and properties Simulate and optimize structures Fabricate structures Measure and confirm properties New material properties!
Veselago We are interested in how waves propagate through various media, so we consider solutions to the wave equation. e,m space 2 E = 2 E em 2 t n = em (-,+) (+,+) k = w em (-,-) (+,-) Sov. Phys. Usp. 10, 509 (1968)
Left-Handed Media: Introduction In 1964, V. G. Veselago contemplated the consequences of a medium having simultaneously negative permeability and permittivity. Out of the three possibilities: 1. Such a condition is precluded by other physical laws. 2. Such a condition is possible, but has no effect on wave propagation, since em is positive. 3. If such a condition occurs, it makes a considerable difference in wave propagation! Out of the three possibilities, it is the third that is predicted.
Fundamental constraints on e and m While we may observe a certain range of values for the material parameters, the only fundamental limit seems to be set by causality, which imposes the ( conditions: we ) r > 1 e ( r w Æ ) = e/ e0 Æ + 1 w ( wm ) w r > 1 m w Æ m/ m0 1 r ( ) = Æ + There appears to be no restriction on negative material constants, other than that frequency dispersion is required.
Left-Handed Media: Introduction In the propagation of electromagnetic waves, the direction of energy flow is given by a right-hand rule, involving E, H, and S: S = E H The propagation, or phase velocity, is usually also determined by a right hand rule: =- B E t k~ E B = me H Thus, when e<0 and m<0, the medium is Left-Handed!
Left-Handed Media: Consequences Group and phase velocities reversed. A Left-Handed medium can be interpreted as having a negative refractive index, n. sinq I = ±n sinq T n = e m Snell s law still holds, but the n must be interpreted accordingly: Unusual and non-intuitive geometrical optics result!
Negative Refractive Index e eff w ( w) = 1 - w 2 p 2 m eff 2 Fw 0 ( w) = 1-2 2 w - w -iwg 0 n( w ) = 1 w ( ) 2 2 2 2 ( w - wb ) w - wp 2 ( w - w 2 0 )
Reversal of Snell s Law Right Handed Media Left Handed Media k T =nk I
Focusing in a Left-Handed Medium RH RH RH RH LH RH n=1 n=1.3 n=1 n=1 n=-1 n=1
Metamaterials Extend Properties e w ( ) = 1 - w w 2 p 2 J. B. Pendry ( ) = 1 - m w w 2 w 2 p - w 2 0
First Left-Handed Test Structure UCSD, PRL 84, 4184 (2000)
Transmission Measurements Transmitted Power (dbm) m>0 e<0 m<0 e<0 Wires alone Split rings alone e<0 m>0 e<0 Wires alone 4.5 5.0 5.5 6.0 6.5 7.0 Frequency (GHz) UCSD, PRL 84, 4184 (2000)
A 2-D Isotropic Structure UCSD, APL 78, 489 (2001)
Example of Utility of Metamaterial The transmission coefficient is an example of a quantity that can be determined simply and analytically, if the bulk material parameters are known. ikd t = exp( - ) s 1 Ê 1ˆ cos( nkd ) - z + 2 Ë z z = m e ( ) = 1 - e w w 2 sin( nkd ) n = 2 w ep 2 - w + e0 ige me 0 ( ) = 1 - m w w 2 2 w mp 2 - w + m0 igm 0 UCSD and ISU, PRB, 65, 195103 (2002)
Measurement of Refractive Index UCSD, Science 292, 77 2001
Measurement of Refractive Index UCSD, Science 292, 77 2001
Measurement of Refractive Index UCSD, Science 292, 77 2001
Transfer matrix is able to find: Transmission (p--->p, p--->s, ) p polarization Reflection (p--->p, p--->s, ) s polarization Both amplitude and phase Absorption Some technical details: Discretization: unit cell N x x N y x N z : up to 24 x 24 x 24 Length of the sample: up to 300 unit cells Periodic boundaries in the transverse direction Can treat 2d and 3d systems Can treat oblique angles Weak point: Technique requires uniform discretization
Structure of the unit cell EM wave propagates in the z -direction Periodic boundary conditions are used in transverse directions Polarization: p wave: E parallel to y s wave: E parallel to x For the p wave, the resonance frequency interval exists, where with Re m eff <0, Re e eff <0 and Re n p <0. For the s wave, the refraction index n s = 1. Typical size of the unit cell: 3.3 x 3.67 x 3.67 mm Typical permittivity of the metallic components: e metal = (-3+5.88 i) x 10 5
Structure of the unit cell: r w SRR EM waves propagate in the z-direction. Periodic boundary conditions are used in the xy-plane d c g LHM
Left-handed material: array of SRRs and wires 10 Resonance frequency as a function of metallic permittivity 0 [GHz] 9 8 7 6 complex e m 5 1000 10000 m Real em
Dependence of LHM peak on metallic permittivity Transmission 10 0 10 2 10 4 10 6 10 8 588000 300000+588000 i 1+588000 i 1 +168000 i 1 + 68000 i 38000 + 38000 i 1+38000 i 1+18000 i 10 10 10 12 6 7 8 9 10 Frequency [GHz] The length of the system is 10 unit cells
Dependence of LHM peak on metallic permittivity
Dependence of LHM peak on L and Im e m 10 1 10 0 10 1 Transmission peak 10 2 10 3 10 4 10 5 10 6 10 7 10 8 588000 300000+588000 i 1+588000 i 1+168000 i 1+68000 i 1+38000 i 1+18000 i 1+8000 i 0 1 2 3 4 5 6 7 8 9 10 11 Length of the system
7.5 8 8.5 9 9.5 10 10.5 Transmission depends on the orientation of SRR Transmission 10 0 10 2 10 4 10 6 10 8 10 10 10 12 # of unit cells N=5 N=10 N=15 N=20 Re(n)<0 Df=0.8GHz Frequency [GHz] Transmission 10 0 # of unit cells: 10 2 N=5 N=10 N=15 N=20 SRR: 10 4 3x3x3 mm 10 6 Re(n)<0 Df=1.2GHz 10 8 8 9 10 11 12 13 Frequency [GHz] Transmission properties depend on the orientation of the SRR: Lower transmission Narrower resonance interval Lower resonance frequency Higher transmission Broader resonance interval Higher resonance frequency
Determination of effective refraction index from transmission studies Refraction index 3 2 1 0 1 2 Im (n) = 10 2 e metal =( 3+5.88 i)x10 5 3 Re n Im n 4 Re(n) < 0 5 8 8.5 9 9.5 10 Frequency [GHz] Refraction index 4 3 2 1 0 1 2 3 4 Im (n) = 10 2 Re(n) <0 e metal =( 3+5.88 i)x10 5 Size of the unit cell: 3.3 x 3.67 x 3.67 mm Re n Im n 9 9.5 10 10.5 11 11.5 12 Frequency [GHz] From transmission and reflection data, the index of refraction n was calculated. Frequency interval with Re n<0 and very small Im n was found.
How do the Left-handed EM waves look like? Real part of the transmission 0.2 f=9.8 GHz n= 3.26 + 0.0155 i 0.1 0.0 0.1 0.2 200 225 250 275 300 System length [#unit cells] Real part of the transmission 1.0 0.5 0.0 0.5 f=11 GHz n= 0.378 + 0.008 i 1.0 0 100 200 300 System length [#unit cells] Real part of the transmission 1.0 f=10.5 GHz n= 1.31 + 0.005 i 0.5 0.0 0.5 1.0 200 225 250 275 300 System length [#unit cells] Propagation through 300 unit cells was simulated. The system length is 1.1 m. Transmission in the LH material is very good.
Length dependence of the transmission Transmission 10 0 10 2 10 4 8.00 8.50 8.70 10 6 9.00 9.30 10 8 9.50 10.00 10 10 10 12 0 5 10 15 20 System length [# of unit cells] Outside the resonance interval, transmission t pp decreases exponentially as the length of the system increases. However, it never decreases below a certain limit. This is due to the non-zero transmission of the p wave into the s wave: t pp =t pp (0) +t ps t ss t sp + The first term, t pp (0) decreases exponentially, the second term does not depend on the system length.
Anisotropy is important outside the LH interval Real (t) 1 0.5 0 0.5 1 1 0.5 0 0.5 1 1 0.5 0 0.5 1 f=9.7 GHz p wave p wave x 100 p wave has the frequency of s wave! s wave 0 10 20 30 40 50 System length [# unit cells] Transmission t pp decreases exponentially only for short system length. For longer systems, t pp =const. There is nonzero transmission t ps of the p wave into the s wave and back. The s wave propagates through the structure without decay. Re n s =1, Im n s =0 t pp =t pp (0) +t ps t ss t sp Therefore, only data for short systems are relevant for the estimation of refraction index of the p wave.
Dependence on the incident angle 10 0 10 0 Incident angle [p] Incident angle [p] 10 2 LHM, 0 LHM, 0.05 10 2 LHM, 0 LHM, 0.05 Transmission 10 4 10 6 10 8 LHM, 0.10 SRR, 0 SRR, 0.05 Transmission 10 4 10 6 LHM, 0.15 10 10 10 8 10 12 8 8.5 9 9.5 10 10.5 Frequency [GHz] 10 10 9 10 11 12 13 Frequency [GHz] Transmission peak does not depend on the angle of incidence! Transition peak strongly depends on the angle of incidence. This structure has an additional xz - plane of symmetry
Dependence of the LHM T peak on the Im e Board Transmission peak 10 0 10 2 10 4 10 6 10 0 10 4 10 8 10 12 Transmission 6 7 8 9 Frequency [GHz] 10 8 10 3 10 2 10 1 10 0 Im e Board Re e Board =3.4 In our simulations, we have: Periodic boundary condition, therefore no losses due to scattering into another direction. Very high Im e metal therefore very small losses in the metallic components. Losses in the dielectric board are crucial for the transmission properties of the LH structures.
Unit cell: 5x3.3x5 mm SRR LHM e metal =( 3+5.88)x10 5 SRR LHM Another 1D left-handed structure: 10 0 10 2 10 0 10 4 Transmission 10 6 10 8 10 10 10 12 Transmission 10 2 10 4 10 6 10 14 6 7 8 9 10 Frequency [GHz] 10 8 10 10 8 8.5 9 9.5 10 Frequency [GHz] Both SRR and wires are located on the same side of the dielectric board. Transmission depends on the orientation of SRR.
Photonic Crystals with negative refraction.
Photonic Crystals with negative refraction. Photonic Crystal vacuum FDTD simulations were used to study the time evolution of an EM wave as it hits the interface vacuum/photonic crystal. Photonic crystal consists of an hexagonal lattice of dielectric rods with e=12.96. The radius of rods is r=0.35a. a is the lattice constant.
Photonic Crystals with negative refraction. t 0 =1.5T T=l/c
Photonic Crystals with negative refraction.
Photonic Crystals with negative refraction.
Photonic Crystals with negative refraction.
Photonic Crystals: negative refraction The EM wave is trapped temporarily at the interface and after a long time, the wave front moves eventually in the negative direction. Negative refraction was observed for wavelength of the EM wave l= 1.64 1.75 a (a is the lattice constant of PC)
Conclusions Simulated various structures of SRRs & LHMs Predicted how resonance transmission frequency depends on the parameters of the system Calculated transmission, reflection and absorption Calculated m eff and e eff and refraction index (with UCSD) Analyzed transmission properties of 2d and 3d structures Found negative refraction in photonic crystals DOE, DARPA, NSF Publications: P. Markos and C. M. Soukoulis, Phys. Rev. B 65, 033401 (2002) P. Markos and C. M. Soukoulis, Phys. Rev. E 65, 036622 (2002) D. R. Smith, S. Schultz, P. Markos and C.M.Soukoulis, Phys. Rev. B 65, 195104 (2002) M. Bayindir, K. Aydin, E. Ozbay, P. Markos and C. M. Soukoulis, Appl. Phys. Lett. (2002) P. Markos, I. Rousochatzakis and C. M. Soukoulis, submitted (2002) S. Foteinopoulou, E. N. Economou and C. M. Soukoulis, submitted (2002)
Metamaterials: New Material Responses for Applications Compact cell phone antennas Improvements to base station antenna Multiple degree-of-freedom antennas for MIMO Active materials for beam steering/smart antennas (SDMA) Phased array antennas Military applications