OSA's Digital Library

Optics Express

Optics Express

  • Editor: C. Martijn de Sterke
  • Vol. 16, Iss. 17 — Aug. 18, 2008
  • pp: 13414–13420
« Show journal navigation

Two-dimensional electromagnetic cloaks with arbitrary geometries

Chao Li and Fang Li  »View Author Affiliations


Optics Express, Vol. 16, Issue 17, pp. 13414-13420 (2008)
http://dx.doi.org/10.1364/OE.16.013414


View Full Text Article

Acrobat PDF (1063 KB)





Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

Transformation optics opens an exciting gateway to design electromagnetic “invisibility” cloaks with anisotropic and inhomogeneous medium. In this paper, we establish a generalized transformation procedure to highly improve the flexibilities for the design of two-dimensional (2D) cloaks. The general expressions for the complex medium parameters are developed, which can be readily applied to design 2D cloaks with arbitrary geometries. An example of 2D cloak with irregular cross section is designed and studied by full-wave simulations. The Huygens’ Principle is applied to quantitatively evaluate its unusual electromagnetic behaviors. All the theoretical and numerical results verify the effectiveness of the proposed approach. The generalization in this Paper makes a great step forward for the flexible design of electromagnetic cloaks with arbitrary shapes.

© 2008 Optical Society of America

1. Introduction

The methods of reducing the electromagnetic (EM) scattering of objects have always been a subject of keen interest in microwave and optical society. Recently, an intriguing idea was proposed by Pendry et. al. for the design of invisible cloaking devices based on a forminvariant transformations of Maxwell’s equations [1–4

J. B. Pendry, D. Schurig, and D. R. Smith, “Controlling electromagnetic fields,” Science 312, 1780 (2006). [CrossRef] [PubMed]

]. The basic principle is to squeeze a volume in a virtual space into a shell with complex medium parameters in the physical space, which excluding EM waves in the concealment volume. This has been demonstrated by experiment in microwave regime using a composite metamaterial made of split ring resonators (SRRs) [2

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 (2006). [CrossRef] [PubMed]

]. At optical wavelengths, a nonmagnetic cloak design has also been proposed that incorporates metallic nanowires [5

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

]. Lately, extensive researches were performed about the theoretical and practical possibility of this ‘invisibility’ [6–9

Z. Ruan, M. Yan, C. W. Neff, and M. Qiu, “Ideal Cylindrical Cloak : Perfect but Sensitive to Tiny Perturbations,” Phys. Rev. Lett. 99, 113903 (2007). [CrossRef] [PubMed]

]. For example, the sensitivity of the field to the tiny perturbations of the cloak’s parameters was studied in Ref. 6

Z. Ruan, M. Yan, C. W. Neff, and M. Qiu, “Ideal Cylindrical Cloak : Perfect but Sensitive to Tiny Perturbations,” Phys. Rev. Lett. 99, 113903 (2007). [CrossRef] [PubMed]

and 7

H. S. Chen, B.-I. Wu, B. L. Zhang, and J. A. Kong, “Electromagnetic wave Interactions with a Metamaterial Cloak,” Phys. Rev. Lett. 99, 063903 (2007). [CrossRef] [PubMed]

, and the method to extending the bandwidth of the cloaks was studied in Ref. 8

H. Y. Chen, Z. X. Liang, P. J. Yao, X. Y. Jiang, H. Ma, and C. T. Chan, “Extending the bandwidth of electromagnetic cloaks,” Phys. Rev. B 76, 241104 (2007). [CrossRef]

. The transformation approach has also been extended to other electromagnetic applications, such as the concentrators [10

M. Rahm, D. Schurig, D. A. Roberts, S. A. Cummer, D. R. Smith, and J. B. Pendry, “Design of Electromagnetic Cloaks and Concentrators Using Form-Invariant Coordinate Transfromations of Maxwell’s Equations,” Photon. Nanostruct.: Fundam. Applic. 6, 87 (2008). [CrossRef]

] and rotating fields [11

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

]. Up to date, all of the previous cloak designs have symmetrical or regular shapes. The mostly often studied structures are the spherical and circular cylindrical cloaks [1–9

J. B. Pendry, D. Schurig, and D. R. Smith, “Controlling electromagnetic fields,” Science 312, 1780 (2006). [CrossRef] [PubMed]

] with rotational symmetry. The elliptic and eccentric elliptic cloaks [12

H. Ma, S. Qu, Z. Xu, J. Q. Zhang, B. W. Chen, and J. F. Wang, “Material parameter equation for elliptical cylindrical cloaks,” Phys. Rev. A 77, 013825 (2007). [CrossRef]

,13

D. H. Kwon and D. H. Werner, “Two-dimensional eccentric elliptic electromagnetic cloaks,” Appl. Phys. Lett. 92, 013505 (2008). [CrossRef]

] were latterly investigated to reduce the symmetry and increase the flexibility of the cloak shapes. The square cloak with flat surfaces and sharp corners was also developed to exemplify the strength of the transformation method [10

M. Rahm, D. Schurig, D. A. Roberts, S. A. Cummer, D. R. Smith, and J. B. Pendry, “Design of Electromagnetic Cloaks and Concentrators Using Form-Invariant Coordinate Transfromations of Maxwell’s Equations,” Photon. Nanostruct.: Fundam. Applic. 6, 87 (2008). [CrossRef]

].

Toward the practical and flexible realizations of EM cloaks, this Paper presents a general transformation method for designing cylindrical cloaks with arbitrarily regular or irregular 2D cross sections. The fundamental strategy is to introduce a new coordinate system conformal to the surface of the scatterer to be concealed. The general expressions for the complex medium parameters are deduced. A 2D cloak device with irregular geometry is designed. Full-wave simulation results are provided for verification. The cloaking performance is then quantitatively evaluated based on the Huygens’ Principle, in which the scatter width is calculated from the simulated near field. The generalization introduced here can be readily specialized to most of the previously designed cloaks, and represents an important step toward the realization of arbitrarily shaped cloaks.

2. Medium transformation for 2D electromagnetic cloaks with arbitrary geometries

In the original space, a point can be described by (x, y, z) in Cartesian coordinate or (r, θ, z) in Cylindrical coordinate with relationship,

x=rcosθ, y=r sinθ.
(1)

For a cylinder with an arbitrary cross section enclosed by a contour r=R0(θ), we first introduce a new coordinate variable defined as

ρ= r R0 (θ)= x2+ y2 R0 (θ)
(2)

to normalize the length of the position vector r⃑ in terms of R0(θ). Then we obtain

x=ρ R0 (θ)cosθ, y=ρ R0 (θ)sinθ.
(3)

Lines of constant-ρ represent a family of contours with similar shapes as r=R0(θ). Hence, (ρ, θ, z) forms a conformal coordinate system with the contour r=R0(θ), as shown in Fig.1(a). To create a cloak, we define a spatial transformation that compresses the cylindrical volume with 0<ρ<1 in the (ρ, θ, z) system into an annular volume with τ<ρ′<1 (τ<1) in the (ρ′, θ′, z′) system via

ρ = τ+ ( 1τ)ρ, θ =θ,z =z.
(4)

Outside this domain the identity transformation is adopted. Finally,

x =r cosθ =ρ R0 (θ )cosθ , y =r sinθ =ρ R0 (θ )sinθ
(5)
Fig. 1. Generalized coordinate transformation. (a) The original coordinate system. (b) The transformed coordinate system. The region with 0<ρ<1 (shaded) in (a) is transformed to the region with τ<ρ′<1 (shaded) in (b).

is applied to complete the whole transformation. The transformed coordinate system is shown in Fig.1(b), where the cylindrical volume with ρ′<τ is completely excluded. Based on above transformation, we can derive the relationship between the transformed Cartesian coordinate (x′, y′, z′) and the original Cartesian coordinate (x, y, z) as

{ x = [τ R0 ( tan 1 yx )+ ( 1τ) x2+ y2] x x2+ y2 y =[τ R0 ( tan 1 yx )+ ( 1τ) x2+ y2] y x2+ y2 z =z
(6)

According to the procedure stated by Pendry et al., the associated permittivity and permeability tensors of the transformation media become [3

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

]

ε i j= det ( Λi i) 1 Λi i Λj j εij, μ i j= det ( Λi i) 1 Λi i Λj i μij,
(7)

where Λi i= q i qi is the Jacobian transformation matrix between the transformed and the original coordinates. After a nontrivial derivation, the tensor components of the transformation media can be written as

εxx= [ r τ R0 (θ )]2 cos2θ + τ2 [ d R0 (θ ) dθ ]2 cos2 θ 2τr d R0 (θ ) dθ sinθ cosθ + r 2 sin2θ r [ r τR0 (θ )]
(8a)
εxy= εyx= τ R0 (θ ) [ 2r τ R0 (θ )]sinθ cosθ + τ2sinθ cosθ [ d R0 (θ ) dθ ]2+τr d R0 (θ ) dθ ( cos2θ sin2θ ) r [ r τ R0 (θ )]
(8b)
εyy= [ r τ R0 (θ )]2 sin2θ + τ2 [ d R0 (θ ) dθ ]2 sin2 θ 2τr d R0 (θ ) dθ sinθ cosθ + r 2 cos2θ r [ r τR0 (θ )]
(8c)
εzz= ( 1 1τ)2 r τR0 (θ ) r
(8d)

and εxz=εyz=εzx=εzy=0. The permeability tensor μ is equal to ε . Here d R0 (θ ) dθ represents the first order derivative of R0(θ′) over θ′, and τ represents the linear compressing ratio. Eq.(8) gives the general expressions of the medium parameters for 2D cloaks with outer boundary defined by r′=R0(θ′) and inner boundary defined by r′=τ R0(θ′). For the special case R0(θ′)=b and τ=a/b, the tensors in Eq.(8) are simplified to the medium parameters of the cylindrical-symmetry cloak [1–4

J. B. Pendry, D. Schurig, and D. R. Smith, “Controlling electromagnetic fields,” Science 312, 1780 (2006). [CrossRef] [PubMed]

]. For elliptical cylindrical cloaks with semi-major axis b and semi-minor axis a, we can insert the contour equation

R0 (θ )= ab a2 cos2θ + b2 sin2θ
(9)

and a desired compressing ratio τ into Eq.(8) to get the media parameters. In fact, R0(θ′) can be chosen as arbitrary continuous functions with period 2π to represent closed contours with arbitrary shapes. It can be generally expressed by a Fourier series as

R0 (θ )= n=0 Ancos (nθ )+ n=1 Bnsin (nθ )
(10)

If d R0 (θ ) dθ is continuous, such as the cylindrical-symmetry cloak and the elliptical cylindrical cloak, the medium parameters will be continuously varying in the cloak region. If d R0 (θ ) dθ is discontinuous in certain θd , which means the cloaks have sharp corners, the medium parameters will also be discontinuous at the corresponding positions. Fortunately, it has been verified by the square cloak in Ref.10

M. Rahm, D. Schurig, D. A. Roberts, S. A. Cummer, D. R. Smith, and J. B. Pendry, “Design of Electromagnetic Cloaks and Concentrators Using Form-Invariant Coordinate Transfromations of Maxwell’s Equations,” Photon. Nanostruct.: Fundam. Applic. 6, 87 (2008). [CrossRef]

that, such discontinuity does not break any fundamental cloaking properties. For a square cloak with a side length 2s 1 for the inner square and a side length 2s 2 for the outer square, we can substitute the contour equation

R0 (θ )= s2 cos ( θ θ0), where θ0= { 0 0<θ < π4, 7π4<θ <2π π2 π4<θ < 3π4 π 3π4<θ < 5π4 3π2 5π4<θ < 7π4
(11)

and τ=s 1/s 2 into Eqs.(8) to obtain the medium parameters. The above discussion means the generalization in this Paper can be specialized to all of the formerly designed cloaks with conformal inner and outer boundaries.

3. Electromagntic properties of a 2D cloak with irregular geometry

To show the flexibility of the approach to design 2D irregular cloaks without any symmetry,

R0 (θ)= 0.7+0.1sin (θ)+0.3sin (3θ)+0.2cos (5θ), τ= 13
(12)

is chosen as an example. To verify the properties of the cloak, full wave simulations based on finite-element method (FEM) are performed. The exciting TM plane wave has an electric filed polarized in z direction with unit amplitude, and incident upon the cloak along the +x direction. The frequency of the harmonic wave is set to be 1 GHz. In simulation, we fill the inner region of the cloak with perfect electric conductor (PEC) and see whether it can be “seen” from outside. Perfect matched layers (PML) are applied to terminate the computational domain in ±x and ±y directions. Fig.2 shows the numerical results for the electric field distribution near the cloaked structure, which is computed with 98456 elements and 394480 unknowns. As can be seen, the wave is smoothly bent around the cloaked area and the phase fronts are perfectly restored when the wave exits the cloak.

Since the proposed cloak has no symmetry in any directions, it’s necessary to study its interaction with EM waves from different orientations. An effective way is to investigate its property under the illumination of a cylindrical wave, which can be decomposed to different planar wave components. In the simulation, a line source is set at the position x=-1m, y=-1m to generate the cylindrical wave. The results are given in Fig.3. It’s seen that the cylindrical wave is perfectly guided around the cloaked object without any obvious scattering.

Fig. 2. Electric field distribution in the vicinity of the cloaked PEC cylinder with irregular shape (excited by a TM plane wave).
Fig. 3. Electric field distribution in the vicinity of the cloaked PEC cylinder with irregular shape (excited by a cylindrical wave).

To quantitatively evaluate the cloaking performance, the scatter width σ(the 2D equivalent of a radar cross section) is calculated based on the Huygens’ Principle. To determine σ, the scattered electric field in far field region is calculated by the integration of the simulated near field along the outer boundary of the scattering object or any other contours which enclosing the scattering object. The integration expression for σ in terms of the near field is

σ= k0 r̂0×C[ ( n̂× Ec) η0 r̂0× ( n̂× Hc)]exp ( ik r· r̂0)dl2 4 Ei2
(13)

where E⃑c and H→s is the EM fields on the integration contour C, 0 is the unit vector of the scattering direction, r→′ is the position vector on the contour C, and η 0 is the free space wave impedance. Considering the non-symmetry of the cloak structure introduced above, the scatter widths for four different incident directions are calculated. The cases with and without cloak are both investigated. The scattering patterns are plotted in Fig.4. Table1 lists some parameters to describe and compare the scattering properties, including the averaged and the maximum scatter widths, and the ratios between the cases with and without cloaks. It’s seen that the cloak greatly reduces the scatter width in different scattering angles. The total scatter power (equivalent to the averaged scatter width) of the irregular PEC cylinder is reduced more than 20 times and the maximum scatter width is reduced more than 90 times. No doubt the ratios could be pushed even larger with finer meshes in simulation. A more interesting phenomenon is that the scattering power of the cloaked structure is almost isotropy over all the angles, which is very different from conventional scattering from objects with irregular shapes.

Fig.4. The scatter width of the cloaked and uncloaked PEC cylinder for different incident directions. The dotted lines are for PEC cylinder without cloak. The solid lines are for PEC cylinder with cloak. The blue, red, black, and green lines are for the incident direction (angle) +x(θ i=0°), -x(θ i=180°), +y(θ i=90°), -y(θ i=270°), respectively.
Table 1.  Some parameters to describe and compare the scatter properties
Incident angle (direction) θ i=0°
(+x)
θ i=180°
(+x)
θ i=90°
(+y)
θ i=270°
(-y)
Averaged Scatter width (m)no cloak (σavg1)1.53561.51511.56041.5610
with cloak (σavg2)0.07270.07240.07330.0732
σavg1 /σavg2 21.12 20.93 21.29 21.32
Maximum Scatter width (m)no cloak (σmax1)9.4919.58710.359.669
with cloak (σmax2)0.10190.10180.09950.0973
σmax1 /σmax2 93.14 94.17 104.1 99.37

In this section, the invisibility of the cloaks to TM waves are numerically verified. Since the relative permeability tensor μ is identical to the relative permittivity tensor ε , the responses of the cloaks to TE waves are the dualtiy of their responses to TM waves. Hence the proposed cloak is also invisible to TE-porlarized waves. Here, the numerical results for TE cases are not included for brevity.

4. Conclusion

A general transformation procedure for designing 2D cloaks with arbitrarily geometries is demonstrated. The general expressions for the complex medium parameters are deduced, which can be readily specialized to most of the previously designed cloaks. A peculiar cloak device with irregular shape is designed as an example and studied by FEM simulations. The Huygens’ Principle is applied to quantitatively evaluate its cloaking properties. All the results verify the flexibility and effectiveness of the proposed method. Although we limit ourselves to 2D cases, the method can be readily extended to construct three dimensional (3D) cloaks with arbitrary shapes. The generalization proposed in this Paper represents an important step and provides a powerful tool toward the flexible design of EM cloaks with arbitrary shapes.

Acknowledgements

This work was supported by the National Natural Science Foundation of China (60501018), the National Basic Research Program of China under Grant (2004CB719800), and the Knowledge Innovation Program of Chinese Academy of Sciences. The authors also acknowledge Master student Kan Yao for his contribution in the numerical simulations.

References and links

1.

J. B. Pendry, D. Schurig, and D. R. Smith, “Controlling electromagnetic fields,” Science 312, 1780 (2006). [CrossRef] [PubMed]

2.

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 (2006). [CrossRef] [PubMed]

3.

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

4.

U. Leonhardt and T. G. Philbin, “General relativity in electrical engineering,” New J. Phys 8, 247 (2006). [CrossRef]

5.

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

6.

Z. Ruan, M. Yan, C. W. Neff, and M. Qiu, “Ideal Cylindrical Cloak : Perfect but Sensitive to Tiny Perturbations,” Phys. Rev. Lett. 99, 113903 (2007). [CrossRef] [PubMed]

7.

H. S. Chen, B.-I. Wu, B. L. Zhang, and J. A. Kong, “Electromagnetic wave Interactions with a Metamaterial Cloak,” Phys. Rev. Lett. 99, 063903 (2007). [CrossRef] [PubMed]

8.

H. Y. Chen, Z. X. Liang, P. J. Yao, X. Y. Jiang, H. Ma, and C. T. Chan, “Extending the bandwidth of electromagnetic cloaks,” Phys. Rev. B 76, 241104 (2007). [CrossRef]

9.

F. Zolla, S. Guenneau, A. Nicolet, and J. B. Pendry, “Cylindrical invisibility cloaks and the mirage effect,” Opt. Lett. 32, 1069 (2007). [CrossRef] [PubMed]

10.

M. Rahm, D. Schurig, D. A. Roberts, S. A. Cummer, D. R. Smith, and J. B. Pendry, “Design of Electromagnetic Cloaks and Concentrators Using Form-Invariant Coordinate Transfromations of Maxwell’s Equations,” Photon. Nanostruct.: Fundam. Applic. 6, 87 (2008). [CrossRef]

11.

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

12.

H. Ma, S. Qu, Z. Xu, J. Q. Zhang, B. W. Chen, and J. F. Wang, “Material parameter equation for elliptical cylindrical cloaks,” Phys. Rev. A 77, 013825 (2007). [CrossRef]

13.

D. H. Kwon and D. H. Werner, “Two-dimensional eccentric elliptic electromagnetic cloaks,” Appl. Phys. Lett. 92, 013505 (2008). [CrossRef]

OCIS Codes
(160.1190) Materials : Anisotropic optical materials
(230.0230) Optical devices : Optical devices
(260.2110) Physical optics : Electromagnetic optics
(160.2710) Materials : Inhomogeneous optical media
(230.3205) Optical devices : Invisibility cloaks

ToC Category:
Physical Optics

History
Original Manuscript: June 4, 2008
Revised Manuscript: August 5, 2008
Manuscript Accepted: August 5, 2008
Published: August 15, 2008

Citation
Chao Li and Fang Li, "Two-dimensional electromagnetic cloaks with arbitrary geometries," Opt. Express 16, 13414-13420 (2008)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-16-17-13414


Sort:  Author  |  Year  |  Journal  |  Reset  

References

  1. J. B. Pendry, D. Schurig, D. R. Smith, "Controlling electromagnetic fields," Science 312, 1780 (2006). [CrossRef] [PubMed]
  2. 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 (2006). [CrossRef] [PubMed]
  3. D. Schurig, J. B. Pendry, and D. R. Smith, "Calculation of material properties and ray tracing in transformation media," Opt. Express 14, 9794 (2006). [CrossRef] [PubMed]
  4. U. Leonhardt, T. G. Philbin, "General relativity in electrical engineering," New J. Phys 8, 247 (2006). [CrossRef]
  5. W. Cai, U. K. Chettiar, A. V. Kildishev, and V. M. Shalaev, "Optical cloaking with metamaterials," Nat. Photonics 1, 224 (2007). [CrossRef]
  6. Z. Ruan, M. Yan. C. W. Neff, and M. Qiu, "Ideal Cylindrical Cloak:Perfect but Sensitive to Tiny Perturbations," Phys. Rev. Lett. 99, 113903 (2007). [CrossRef] [PubMed]
  7. H. S. Chen, B.-I. Wu, B. L. Zhang, and J. A. Kong, "Electromagnetic wave Interactions with a Metamaterial Cloak," Phys. Rev. Lett. 99, 063903 (2007). [CrossRef] [PubMed]
  8. H. Y. Chen, Z. X. Liang, P. J. Yao, X. Y. Jiang, H. Ma, and C. T. Chan, "Extending the bandwidth of electromagnetic cloaks," Phys. Rev. B 76, 241104 (2007). [CrossRef]
  9. F. Zolla, S. Guenneau, A. Nicolet, J. B. Pendry, "Cylindrical invisibility cloaks and the mirage effect," Opt. Lett. 32, 1069 (2007). [CrossRef] [PubMed]
  10. M. Rahm, D. Schurig, D. A. Roberts, S. A. Cummer, D. R. Smith, and J. B. Pendry, "Design of Electromagnetic Cloaks and Concentrators Using Form-Invariant Coordinate Transfromations of Maxwell??s Equations," Photon. Nanostruct.: Fundam. Applic. 6, 87 (2008). [CrossRef]
  11. H. Chen and C. T. Chan, "Transformation media that rotate electromagnetic fields," Appl. Phys. Lett. 90, 241105 (2007). [CrossRef]
  12. H. Ma, S. Qu, Z. Xu, J. Q. Zhang, B. W. Chen, and J. F. Wang, "Material parameter equation for elliptical cylindrical cloaks," Phys. Rev. A 77, 013825 (2007). [CrossRef]
  13. D. H. Kwon, D. H. Werner, "Two-dimensional eccentric elliptic electromagnetic cloaks," Appl. Phys. Lett. 92, 013505 (2008). [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