OSA's Digital Library

Optics Express

Optics Express

  • Editor: C. Martijn de Sterke
  • Vol. 16, Iss. 10 — May. 12, 2008
  • pp: 6809–6814
« Show journal navigation

A cone-shaped concentrator with varying performances of concentrating

Lan Lin, Wei Wang, Chunlei Du, and Xiangang Luo  »View Author Affiliations


Optics Express, Vol. 16, Issue 10, pp. 6809-6814 (2008)
http://dx.doi.org/10.1364/OE.16.006809


View Full Text Article

Acrobat PDF (267 KB)





Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

A three-dimensional cone-shaped concentrator was designed and analyzed through an approach of coordinate transformation theory. The device can provide varying performances for concentrating along the symmetric axis. The physical picture regarding concentrating ability of this structure was revealed and quantitative analyses were performed for the purpose of investigating the dependence of the concentrating properties on the structural parameters. Moreover, reduced material parameters were theoretically derived and the corresponding mismatched impedance at boundaries was analyzed. Finite element method-based numerical simulations results of the device were further presented to verify our theoretical design.

© 2008 Optical Society of America

1. Introduction

Since control of electromagnetic waves using the coordinate transformation was proposed by Ward, et al., [1

1. A. J. Ward and J. B. Pendry, “Refraction and geometry in Maxwell’s equation,” J. Mod. Opt. 43, 773–793 (1996). [CrossRef]

], various types of devices, recently termed transformation media, have been suggested for the purpose of realization of some specific electromagnetic functions. The cylindrical and sphere cloaks were put forth firstly based on this theory [2

2. J. B. Pendry, D. Schurig, and D. R. Smith, “Controlling Electromagnetic Fields,” Science 312, 1780–1782 (2006). [CrossRef] [PubMed]

, 3

3. D. Schurig, J. J. Mock, B. J. Justice, S. A. Cummer, J. B. Pendry, A. F. Starr, and D. R. Smith, “Metamaterial electromagnetic cloak at microwave frequencies,” Science 314, 977–980 (2006). [CrossRef] [PubMed]

]. Other optical transformation media such as magnifying perfect lens and superlens were designed by applying this method [4

4. M. Tsang and D. Psaltis, “Magnifying perfect lens and superlens design by coordinate transformation,” http://www.arxiv.org:physics/0708.0262v1, (2007).

]. This approach was confirmed by interpreting the coordinate transformation in terms of Jacobean matrix called form-invariant transformation of Maxwell’s equation by D.Schurig later [5

5. D. Schurig, J. B. Pendry, and D. R. Smith, “Calculation of material properties and ray tracing in transformation media,” Opt. Express 14, 9794–9804 (2006) [CrossRef] [PubMed]

]. Consequently, an electromagnetic fields rotator, a square-shaped cloak were designed using this method [6

6. H. Chen and C. T. Chan, “Transformation media that rotate electromagnetic fields,” Appl. Phys. Lett , 90, 241105 (2007). [CrossRef]

, 7

7. M. Rahm, D. Schurig, D. A. Roberts, S. A. Cummer, and D. R. Smith, “Design of electromagnetic cloaks and concentrators using form-invariant coordinate transformations of Maxwell’s Equations,” http://www.arxiv.org:physics/0706.2452v1, (2007).

]. Most recently, the appearance of metamaterials provides a way to realize these anisotropic and inhomogeneous materials, and a cloak in microwave range has been experimentally demonstrated [3

3. D. Schurig, J. J. Mock, B. J. Justice, S. A. Cummer, J. B. Pendry, A. F. Starr, and D. R. Smith, “Metamaterial electromagnetic cloak at microwave frequencies,” Science 314, 977–980 (2006). [CrossRef] [PubMed]

, 8

8. W. Cai, U. K. Chettiar, A. Kildishev, and V. M. Shalaev, “Optical cloaking with metamaterials,” Nat. Photonics 1, 224–227 (2007). [CrossRef]

]. Since the material parameters are complex, a way to derive an imperfect but simpler version of the materials parameters was offered to make the fabrication easier [3

3. D. Schurig, J. J. Mock, B. J. Justice, S. A. Cummer, J. B. Pendry, A. F. Starr, and D. R. Smith, “Metamaterial electromagnetic cloak at microwave frequencies,” Science 314, 977–980 (2006). [CrossRef] [PubMed]

, 9

9. S. A. Cummer, B. I. Popa, D. Schurig, and D. R. Smith, “Full-wave simulations of electromagnetic cloaking structures,” Phy. Rev. E 74, 036621 (2006). [CrossRef]

].

A cylindrical electromagnetic concentrator that makes the incident wave power flow concentrate within the inner region was proposed [7

7. M. Rahm, D. Schurig, D. A. Roberts, S. A. Cummer, and D. R. Smith, “Design of electromagnetic cloaks and concentrators using form-invariant coordinate transformations of Maxwell’s Equations,” http://www.arxiv.org:physics/0706.2452v1, (2007).

]. It is expected to be of potential use in the applications where high field intensities are required. However, the constant enhancement of the energy along the symmetric axis in the cylindrical concentrator can not satisfy some requirements where various enhancement of energy are needed. In this work, a cone-shaped concentrator was proposed to provide varying energy distribution along the symmetry axis. The design method and approach to derive the material parameters based on basic theory of Ward’s coordinate transformation are given in the following section. Then the numerical simulation are presented and analyzed. Finally, several sets of reduced material parameters are discussed.

2. Principle

Concentrator is a device to collect the energy into a small region and the schematic diagram of the cone-shaped concentrator is represented in Fig. 1(a). When the electromagnetic wave propagates into the region enclosed by the dashed circle, the energy is concentrated into the yellow region. For this purpose, anisotropic and inhomogeneous material is needed which can be designed by the approach of coordinate transformation according to the change of the space.

Fig. 1. (Color online) (a). A schematic diagram of our proposed cone-shaped concentrator. The yellow region is the concentrating area. The incident wave is normal to the symmetric axis. (b). A schematic diagram of the space transformation in a cross-section where region I (the concentrating area in yellow color) and region II (blue) are marked respectively.

To redistribute the field, a two-steps space deformation is necessary as suggested in Ref. [7

7. M. Rahm, D. Schurig, D. A. Roberts, S. A. Cummer, and D. R. Smith, “Design of electromagnetic cloaks and concentrators using form-invariant coordinate transformations of Maxwell’s Equations,” http://www.arxiv.org:physics/0706.2452v1, (2007).

]: the first one is compressing the space within (0≤r<R 2) into a circle with a radius of R 1, named as region I hereinafter; the second one is expanding the space between R2 and R3 to a region between R1 and R3 so as to make the space continuous, named region II hereinafter. Figure 1(b) is the schematic diagram of the space transformation. Note that Ri(i=1, 2, 3) is a function of z in this device. The process can be mathematically described as:

r={R1R2r(0rR2)R3R1R3R2rR2R1R3R2R3(R2rR3)
θ=θ
z=z
(1)

According to the theory of coordinate transformation in Ref. [1

1. A. J. Ward and J. B. Pendry, “Refraction and geometry in Maxwell’s equation,” J. Mod. Opt. 43, 773–793 (1996). [CrossRef]

], we firstly calculate the transformation coefficients with respect to different regions in the cylindrical system:

In region I :

Qr=Qθ=R2R1,Qz=1
(2)

In region II :

Qr=R3R2R3R1
Qθ=[1+R3(R2R1)r(R3R2)]R3R2R3R1
Qz=1
(3)

So the material parameters in the cylinder coordinate are

ηij=ηQ1Q2Q3(QiQj)1η=ε,μ
(4)
εij=μij={(10001000(R2R1)2)(0rR1)(R2R1R3R2R3r+1000(R2R1R3R2R3r+1)1000(R3R2R3R1)2×R2R1R3R2R3r+1)(R1rR3)
(5)
k=R2R1R3R2R3r+1
(6)

Here, we designed a structure with parameters of R1=R10-zt, R2=R20-zt, R3=R30-zt, where z is the symmetric axis, Ri0(i=1, 2, 3) is the radius at z=0, and t is the slope of the cone. The ratio R2/R1 that affects the concentrating performance greatly can be expressed as

R2R1=R20z×tR10z×t=1+R20R10R10z×t
(7)

It can be seen from Eq. (5) that the parameters in region II are anisotropic and inhomogeneous which is difficult to fabricate using current technologies, so reduced parameters are required. For a complete EM mode concentrator, we have 6 material parameters taken into account, including both the permittivity and permeability in the radial (r), azimuthally (φ) and z directions which are hard to satisfy at the same time. Here a TE mode design is considered and the TM mode design follows the same principle by making εµ and µε substitutions. In the TE mode design case, only εz, µφ and µr enter into the Maxwell’s equations. Therefore, we just consider the εz, µφ and µr, and assume that the other three parameters, µz, εφ and εr, are units which means the TM mode energy will remain the same as in free space. The reduced parameters that make the same solutions for the wave equations are obtained as long as the product of µrεz and µφεz remain the same value [3

3. D. Schurig, J. J. Mock, B. J. Justice, S. A. Cummer, J. B. Pendry, A. F. Starr, and D. R. Smith, “Metamaterial electromagnetic cloak at microwave frequencies,” Science 314, 977–980 (2006). [CrossRef] [PubMed]

, 9

9. S. A. Cummer, B. I. Popa, D. Schurig, and D. R. Smith, “Full-wave simulations of electromagnetic cloaking structures,” Phy. Rev. E 74, 036621 (2006). [CrossRef]

]. When εz and µφ are constants, only µr is dependent on the coordinate; and there are numerous sets of the parameters, theoretically. Considering the convenience of fabrication and the matching condition at the boundary, three sets of the reduced material parameters are derived. For the first set, εz=1, which is easy to realize; for the second set, the parameters in φ and z directions are the same, and the impedance matches at the boundary of r=R3; and the third set has an impedance matched at the boundary of r=R 1. They can be written in the following equations respectively.

Set1:εz=1,μφ=(R3R2R3R1)2,μr=(R3R2R3R1)2k2
(8)
Set2:εz=μφ=R3R2R3R1,μr=k2R3R2R3R1
(9)
Set3:εz=R3R2R3R1×R2R1,μφ=R3R2R3R1×R1R2,μr=k2×R3R2R3R1×R1R2
(10)

3. Ideal material parameters of the cone-shaped concentrator

Numerical simulation results for the ideal concentrator using finite element method are given in Fig. 2(a). The simulation area is 30 cm×30 cm×15 cm. Origin of the coordinate is at the center of the cone, and z axis is the symmetric axis. The incident wave is TE polarized with the frequency of 4GHz, which propagates along the y axis with parameters of R10=1.5 cm, R20=3.5 cm, R30=5.5 cm and t=1/15. Figure 2(a) shows the distribution of Ez in x-y plane and y-z plane respectively. And the distribution of Ez in a cross-section of the cone parallel to the x-y plane at z=6 cm with R1=1.1 cm, R2=3.1 cm and R3=5.1 cm is presented in Fig. 2(b). The electric field is concentrated to the region with a radius of R1. The wavelength in free space is about R2/R1 times of that in region I.

Fig. 2. (Color online) (a). Distribution of Ez in the calculated region. The x-y plane is at z=0 and the y-z plane is at x=0. (b). Distribution of Ez at x-y plane at site of z=6 cm.

It can be seen from Fig. 2 that the phase velocity in region II which is caused by the distribution of the inhomogeneous refraction index defined as

nφ=μrεz=R3R2R3R1k
(11)
nr=μφεz=R3R2R3R1
(12)

Refraction index in the region II decrease with the increase of the radius, so a smaller r will lead to a larger wave number and a smaller velocity.

Considering the space deformation, the enhancement of energy in region I is proportion to the square of R2/R1, because we have compressed the space with radius of R2 into a region with radius of R1 in the first step of the transformation. Figure 3(a) shows the energy distribution along the z axis at x=y=0. It can be seen that the energy varies with different position of z, which means that the ability of concentrating is not constant along the symmetrical axis. Figure 3(b) represents the relationship between R2/R1 and the concentrating ability. The energy distribution in free space is marked by the red line, and the average value along the line is 4.5×10-12 J/m2. The green curve represent energy distribution along the y axis of the concentrator at site of x=0 and z=7.5 cm with R1=1 cm, R2=3 cm and R3=5 cm, so R2/R1=3. Obviously, the energy in the center, whose average value is 4.25×10-11 J/m2, is about 9 times larger than that of the value in free space. A similar result can be obtained from the yellow curve which represents energy distribution along y axis at x=0 and z=0 with R1=1.5 cm, R2=3.5 cm and R3=5.5 cm. In this case R2/R1=2.3, the averaged energy value in region I is 2.4×10e-11 J/m2, so the enhancement obtained from Fig. 3(b) is ~5.3. The averaged energy in region I of the blue line at x=0 and z=-7.5 cm with R1=2 cm, R2=4 cm and R3=6 cm is 1.68×10-11 J/m2 which is almost 4 times larger than the case in free space. The coarseness of the lines is caused by the finite number from the finite elements.

Fig. 3. (Color online) (a). Energy distribution along the z axis at the site of x=0 and y=0. (b). Energy distribution in free space is shown by the red line. The other three lines are obtained through the concentrator along the y axis at x=0 with different positions of z.

4. Reduced parameters of the cone-shaped concentrator

The designed parameters are all complicated and position-relevant, and difficult to realize in practice. Therefore, sets of reduced parameters are proposed here based on the Maxwell’s equations and wave equations. The physical insight of this change is that the reduced parameters have the same refraction index which controls the trace of the light as the ideal ones with the penalty of causing the scattering due to the impedance mismatch at the interfaces. In order to keep the product of µrεz and µφεz to be the same, three sets of reduced parameters have been derived in section2 for their specificity. Figure 4 shows the impedance of four conditions at the site of x=z=0. For the ideal case, the impedance is Z=(R3-R1)/[(R3-R2)k], which is matched at both boundaries with r=R3 and r=R1. The impedance of the reduced parameters is Z=(R3-R2)/(R3-R1) for set 1 (green line), and it does not match the impedance at both boundaries of r=R1 and r=R3 in this case; it can be matched at the inner boundary at r=R1 only in the case of R3=R2+R1. The impedance is Z=1 for set 2 which is matched at r=R3 and Z=R1/R2 for set 3 matched at r=R1. Hence, we can choose the proper set of the reduced material parameters as desired.

Fig. 4. (Color online) Impedance of four conditions at the site of x=z=0. The blue line represents the ideal case. The other three lines are derived by three sets of the reduced parameters obtained in section 2.

5. Conclusion

We have presented a cone-shaped concentrator with various abilities for concentrating along the symmetric axis based on the coordinate transformation in this paper. Both cases of ideal and reduced material parameters are given respectively. Electromagnetic behavior of the device was simulated and analyzed by means of a finite element method. Our numerical simulation results demonstrated that the proposed configuration is an appropriate structure for obtaining varying energy enhancements. This structure can be extended to design a concentrator with various functioning effects along the symmetric axis as we needed by specifying the generatrix of the structure.

Acknowledgments

This work was supported by 973 Program of China (No.2006CB302900) and 863 Program of China (2006AA04Z310).

References and links

1.

A. J. Ward and J. B. Pendry, “Refraction and geometry in Maxwell’s equation,” J. Mod. Opt. 43, 773–793 (1996). [CrossRef]

2.

J. B. Pendry, D. Schurig, and D. R. Smith, “Controlling Electromagnetic Fields,” Science 312, 1780–1782 (2006). [CrossRef] [PubMed]

3.

D. Schurig, J. J. Mock, B. J. Justice, S. A. Cummer, J. B. Pendry, A. F. Starr, and D. R. Smith, “Metamaterial electromagnetic cloak at microwave frequencies,” Science 314, 977–980 (2006). [CrossRef] [PubMed]

4.

M. Tsang and D. Psaltis, “Magnifying perfect lens and superlens design by coordinate transformation,” http://www.arxiv.org:physics/0708.0262v1, (2007).

5.

D. Schurig, J. B. Pendry, and D. R. Smith, “Calculation of material properties and ray tracing in transformation media,” Opt. Express 14, 9794–9804 (2006) [CrossRef] [PubMed]

6.

H. Chen and C. T. Chan, “Transformation media that rotate electromagnetic fields,” Appl. Phys. Lett , 90, 241105 (2007). [CrossRef]

7.

M. Rahm, D. Schurig, D. A. Roberts, S. A. Cummer, and D. R. Smith, “Design of electromagnetic cloaks and concentrators using form-invariant coordinate transformations of Maxwell’s Equations,” http://www.arxiv.org:physics/0706.2452v1, (2007).

8.

W. Cai, U. K. Chettiar, A. Kildishev, and V. M. Shalaev, “Optical cloaking with metamaterials,” Nat. Photonics 1, 224–227 (2007). [CrossRef]

9.

S. A. Cummer, B. I. Popa, D. Schurig, and D. R. Smith, “Full-wave simulations of electromagnetic cloaking structures,” Phy. Rev. E 74, 036621 (2006). [CrossRef]

OCIS Codes
(160.1190) Materials : Anisotropic optical materials
(260.2110) Physical optics : Electromagnetic optics
(260.2710) Physical optics : Inhomogeneous optical media

ToC Category:
Physical Optics

History
Original Manuscript: February 4, 2008
Revised Manuscript: March 22, 2008
Manuscript Accepted: April 19, 2008
Published: April 28, 2008

Citation
Lan Lin, Wei Wang, Chunlei Du, and Xiangang Luo, "A cone-shaped concentrator with varying performances of concentrating," Opt. Express 16, 6809-6814 (2008)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-16-10-6809


Sort:  Author  |  Year  |  Journal  |  Reset  

References

  1. A. J. Ward and J. B. Pendry, "Refraction and geometry in Maxwell's equation," J. Mod. Opt. 43, 773-793 (1996). [CrossRef]
  2. J. B. Pendry, D. Schurig, and D. R. Smith, "Controlling Electromagnetic Fields," Science 312, 1780-1782 (2006). [CrossRef] [PubMed]
  3. D. Schurig, J. J. Mock, B. J. Justice, S. A. Cummer, J. B. Pendry, A. F. Starr, and D. R. Smith, "Metamaterial electromagnetic cloak at microwave frequencies," Science 314, 977-980 (2006). [CrossRef] [PubMed]
  4. M. Tsang and D. Psaltis, "Magnifying perfect lens and superlens design by coordinate transformation," http://www.arxiv.org:physics/0708.0262v1 (2007).
  5. D. Schurig, J. B. Pendry, and D. R. Smith, "Calculation of material properties and ray tracing in transformation media," Opt. Express 14, 9794-9804 (2006) [CrossRef] [PubMed]
  6. H. Chen and C. T. Chan, "Transformation media that rotate electromagnetic fields," Appl. Phys. Lett.  90, 241105 (2007). [CrossRef]
  7. M. Rahm, D. Schurig, D. A. Roberts, S. A. Cummer, and D. R. Smith, " Design of electromagnetic cloaks and concentrators using form-invariant coordinate transformations of Maxwell's Equations," http://www.arxiv.org:physics/0706.2452v1. (2007).
  8. W. Cai, U. K. Chettiar, A. Kildishev, and V. M. Shalaev, "Optical cloaking with metamaterials," Nat. Photonics 1, 224-227 (2007). [CrossRef]
  9. S. A. Cummer, B. I. Popa, D. Schurig, and D. R. Smith, "Full-wave simulations of electromagnetic cloaking structures," Phy. Rev. E 74, 036621 (2006). [CrossRef]

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

Figures

Fig. 1. Fig. 2. Fig. 3.
 
Fig. 4.
 

« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited