Electromagnetic cloaking method

Abstract

A method of constructing a concealing volume comprises constructing a plurality of concealing volume elements around a concealable volume. Each concealing volume element has a material parameter arranged to direct a propagating wave around the concealable volume. The material parameter can be refractive index, electrical permittivity, and magnetic permittivity. The concealing volume can be a metamaterial. The concealing volume diverts incoming propagating waves such that outgoing propagating waves appear to be unperturbed to an observer.

Claims

1. A method, comprising: deflecting an electromagnetic wave that would pass through a selected region to avoid the selected region; and restoring the deflected electromagnetic wave so that it appears to have passed through the selected region; wherein the electromagnetic wave is a radio-frequency or microwave-frequency electromagnetic wave.

2. The method of claim 1, wherein the deflecting and the restoring include deflecting and restoring by refracting the electromagnetic wave in a region that encloses the selected region.

3. The method of claim 2, wherein the refracting includes artificially magnetically responding to the electromagnetic field.

4. The method of claim 3, wherein the artificially magnetically responding includes inducing currents at a discrete plurality of locations within the enclosing region, the locations having a spacing much larger than an atomic spacing but much smaller than a wavelength of the electromagnetic wave.

5. The method of claim 2, wherein the refracting includes artificially electrically responding to the electromagnetic field.

6. The method of claim 5, wherein the artificially electrically responding includes inducing plasmons at a discrete plurality of locations within the enclosing region, the locations having a spacing much larger than an atomic spacing but much smaller than a wavelength of the electromagnetic wave.

7. The method of claim 2, wherein the refracting is equivalent to a distorting of a geometry of space as seen by the electromagnetic wave.

8. The method of claim 7, wherein the distorting of the geometry of space corresponds to a coordinate transformation that maps the selected region into the enclosing region.

9. The method of claim 8, wherein the enclosing region is a spherical or cylindrical annulus.

10. An apparatus, comprising: an artificially structured material configured to deflect an electromagnetic wave that would pass through a selected region to avoid the selected region and to restore the deflected electromagnetic wave so that it appears to have passed through the selected region; wherein the electromagnetic wave is a radio-frequency or microwave-frequency electromagnetic wave.

11. The apparatus of claim 10, wherein the artificially structured material is a refractive material in a region that encloses the selected region.

12. The apparatus of claim 11, wherein the artificially structured material includes a plurality of structured objects having a size and spacing much smaller than a wavelength of the electromagnetic wave.

13. The apparatus of claim 12, wherein the plurality of structured objects includes an array of non-magnetic wire loops producing an effective magnetic response.

14. The apparatus of claim 13, wherein the non-magnetic wire loops are split ring resonators.

15. The apparatus of claim 12, wherein the plurality of structured objects includes an array of thin wires providing an effective electric response.

16. The apparatus of claim 12, wherein the plurality of structured objects provides an effective permittivity and an effective permeability for the artificially structured material.

17. The apparatus of claim 16, wherein the effective permittivity and the effective permeability correspond to a coordinate transformation that maps the selected region into the enclosing region.

18. The apparatus of claim 17, wherein the enclosing region is a spherical or cylindrical annulus.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) FIG. 1 is a cross-sectional ray trajectory diagram for a sphere of radius r<R.sub.1 cloaked within an annulus of radius R.sub.1<r<R.sub.2, assuming that R.sub.2>>;

(2) FIG. 2a shows a field line in free space plotted on a Cartesian mesh;

(3) FIG. 2b shows the field line and background coordinates of FIG. 2a both distorted by the same transformation;

(4) FIG. 2c shows a small element resembling a parallelpiped;

(5) FIG. 2d shows an integration path for solving Maxwell's equations;

(6) FIG. 3 is a cross-sectional view of the electrostatic displacement field lines for a sphere of radius r<R.sub.1 cloaked within an annulus of radius R.sub.1<r<R.sub.2, located close to a point charge, assuming that R.sub.2<<;

(7) FIG. 4a is a three-dimensional depiction of the diversion of electromagnetic rays around a cloaked volume, according to the present invention;

(8) FIG. 4b is a three-dimensional depiction of the diversion of electromagnetic rays around a cloaked surface, according to the present invention;

(9) FIG. 5 shows a front perspective view of a Swiss-roll metamaterial;

(10) FIG. 6 shows a top view of the Swiss-roll metamaterial of FIG. 5;

(11) FIG. 7 shows a single split ring metamaterial;

(12) FIG. 8 shows a plan view of a split ring structure in a square array, lattice spacing ;

(13) FIG. 9 depicts the building of three-dimensional split ring resonator (SRR) at unit cell level;

(14) FIG. 10 is a prior-art negative index metamaterial comprising SRR's and wires;

(15) FIG. 11 is a cross-sectional ray trajectory diagram of a concentration device of radius R.sub.3 in which a spherical volume of space radius R.sub.2 is compressed to a smaller volume of space radius R.sub.1; and

(16) FIG. 12 is a cross-sectional ray trajectory diagram in which rays from a point source incident on a sphere radius R.sub.3 are concentrated into an inner sphere radius R.sub.1.

(17) FIG. 13 is a flow chart depicting the steps of method 1300: identifying a coordinate transform (1302); computing the material parameters (1304); selecting diverting volume elements (1306); and constructing the diverting volume (1308).

DETAILED DESCRIPTION

(18) In overview, according to a first embodiment the invention provides a method of concealing an object or volume of space by means of redirecting the electromagnetic field lines around it. Any field lines which normally would flow through the concealed volume are redirected through a surrounding concealing volume. Outside of the cloaking structure, the field lines remain unperturbed. The cloaking structure is formed from a metamaterial, that is a material designed to exhibit the specific material parameters required for cloaking the object or volume in question. In particular, the electric permittivity and magnetic permeability and hence the refractive index n of the cloaking metamaterial can be varied effectively continuously spatially in order to manipulate the electromagnetic field lines and hence effectively conceal the enclosed volume. As a result of the invention, the concealed volume is undetectable to an outside observer, and any electromagnetic radiation produced by an object within the concealed volume cannot escape from it.

(19) In a further embodiment the invention provides a method of concentrating the electromagnetic fields incident upon the outer surface of a concentrator device into an inner core region. The concentrator structure comprises a concentrating volume and a collecting volume. Any field lines which would normally flow through the concentrating volume are redirected into an inner collecting volume. Outside of the concentrator annulus surrounding the concentrating volume, the field lines remain unperturbed. The concentrator structure is formed from a metamaterial in a similar manner to that described above for the cloaking structure above.

(20) Whilst, as is outlined in more detail below, the transformation method described herein can be applied to conceal any desired volume using any desired cloaking structure, for simplicity it will first be described in relation to the simple spherical case as shown in FIG. 1. According to FIG. 1, the volume to be hidden is a sphere 100 of radius R.sub.1 and the cloaking material is contained within the annulus 102 of radius R.sub.1<r<R.sub.2. In order to create a secure, radiation-free volume within the sphere 100 of radius R.sub.1, all those electromagnetic field lines 104 which normally would flow through the sphere 106 of radius r<R.sub.2 must therefore be compressed into the cloaking annulus 102 R.sub.1<r<R.sub.2. That is, in order to remove the sphere 100 of radius R.sub.1 from view, the volume within that sphere 100 plus the volume within the annulus 102 must all be mapped onto the annulus 102 of radius R.sub.1<r<R.sub.2. In the medium outside of the sphere 106 r<R.sub.2, the electromagnetic field lines 104 must remain unperturbed.

(21) The theory of propagation of electromagnetic waves is well known and so only the relevant aspects are summarised here in the context of the invention. In free space, the magnetic field strength H is related to the magnetic field intensity B and the magnetic permeability of free space, .sub.0, via the equation B=.sub.0H. The electric field E in free space is related to the electric displacement field D and electric permittivity of free space .sub.0 via the equation D=.sub.0E. In a dielectric the value of and are modified by the relative permittivity and permeability value .sub.r, .sub.r and the refractive index of the invention is given by n.sub.r=.sup.{square root over ()}. As will be well known to the skilled person in the art, the three vectors E, H, and S are mutually orthogonal, with the Poynting vector S being calculable via the vector cross product S=EH. The Poynting vector S gives the direction and magnitude of energy flow per unit area per unit time of an electromagnetic wave, and in a ray trajectory diagram such as that shown in FIG. 1, the direction of the rays of light correspond to the electromagnetic field lines and follow the direction of the Poynting vector. In order to implement this invention therefore, it is necessary to focus or redirect any of the electric displacement field D, the magnetic field intensity B, and the Poynting vector S so that they avoid the sphere 100 of radius R.sub.1, and return undisturbed to their original trajectories.

(22) As described in more detail below the approach further requires solution of Maxwell's equations in relation to the redirected fields. For ease of depiction, redirection of the Poynting vector in particular will be addressed in the following description.

(23) In order to calculate the electromagnetic properties which a cloaking structure must possess in order to effectively implement the desired cloaking effects, the first step is to calculate the distortion which must be applied to the Poynting vector S so that the rays do not pass through the sphere 100 of radius R.sub.1. As shown in FIGS. 2a and 2b, this is achieved by recording the initial configuration of each ray, which would occur in the absence of a cloaked volume and cloaking material, on a known geometric mesh, such as a spherical or Cartesian mesh. For example in FIG. 2a, a single ray 200 is plotted on a Cartesian mesh 202: Once the appropriate rays 200 have been plotted on it, the Cartesian mesh 202 is stretched or distorted in such a manner that the rays plotted thereon become distorted to follow their new desired path, avoiding the cloaked volume. The distortions which are imposed on the Cartesian mesh 202 in order to achieve this can then be recorded as a coordinate transformation between the original Cartesian mesh 202 and a second transformed mesh 204, wherein the second mesh is defined by:
q.sub.1(x,y,z),q.sub.2(x,y,z),q.sub.3(x,y,z)(1)

(24) In this second mesh 204, lines of constant q.sub.2, q.sub.3 define the generalized q.sub.1 axis, and so on. It will be seen that conversely if the second mesh 204 comprised a set of points defined by equal increments along the q.sub.1, q.sub.2, q.sub.3 axes, it would appear distorted in the x, y, z Cartesian frame, as a comparison of FIGS. 2a and 2b clearly demonstrates.

(25) The benefit of performing this transformation is that, as described below, it can be shown that it is possible to solve Maxwell's equations in the same form in the transformed geometry as in the original geometry, with only a corresponding transform of the spatially distributed values of and being required. This can be understood intuitively from the recognition that variation of the permittivity E and permeability values, and hence the refractive index n in a space, is equivalent to distorting the geometry of the space as seen by a light ray. It will be noted that the approach specifically addresses modification of the value of and , individually rather than of the refractive index n=+{square root over ()}. In particular this allows reflection between interfaces of different media to be removed by ensuring impedance (Z) matching where

(26) $Z=\sqrt{\frac{}{\mathrm{.Math.}}}.$

(27) As will be known to the skilled person, Maxwell's equations in a system of Cartesian coordinates take the form:
E=.sub.0H/t,
E=.sub.0E/t(2)
where both electric permittivity and magnetic permeability depend on position. When a coordinate transformation is applied to these equations, in the new co-ordinate system they take the form:
.sub.q=.sub.0{circumflex over ()}/t,
.sub.q=.sub.0{circumflex over ()}/t(3)
where and are in general tensors, and , and are renormalized electric and magnetic fields. All four quantities are simply related to originals. In other words the form of Maxwell's equation is preserved by a co-ordinate transformation: the co-ordinate transformation does not change the fact that we are still solving Maxwell's equations, it simply changes the definition of , .

(28) In particular, the effect of the transformation of and is to scale them both by a common factor. A generalized solution to this is as follows.

(29) We define three units vectors, u.sub.1, u.sub.2, u.sub.3, to point along the generalized q.sub.1, q.sub.2, q.sub.3 axes. The length of a line element is given by,
ds.sup.2=dx.sup.2+dy.sup.2+dz.sup.2=Q.sub.11dq.sub.1.sup.2+Q.sub.22M+Q.sub.33dq.sub.3.sup.2+2Q.sub.12dq.sub.1dq.sub.2+2Q.sub.13dq.sub.1dq.sub.3+2Q.sub.23dq.sub.2dq.sub.3(G1)
where,

(30) $\begin{array}{cc}{Q}_{\mathrm{ij}}=\frac{x}{q_i}\frac{x}{q_j}+\frac{y}{q_i}\frac{y}{q_j}+\frac{z}{q_i}\frac{z}{q_j}& \left(\mathrm{G2}\right)\end{array}$
In particular we shall need the length of a line element directed along one of the three axes,
ds.sub.i=Q.sub.idq.sub.i(G3)
where for shorthand,
Q.sub.i.sup.2=Q.sub.ii(G4)
To calculate E consider a small element, small enough that it resembles a parallelepiped (FIG. 2c). In this we assume that the transformation has no singularities such as points or lines where the co-ordinate system suddenly heads off in a different direction.
First we calculate the projection of E onto the normal to the u.sub.1-u.sub.2 plane by taking a line integral around the u.sub.1-u.sub.2 parallelogram and applying Stokes' theorem (FIG. 3). We define,
E.sub.1=E.Math.u.sub.1,E.sub.2=E.Math.u.sub.2,E.sub.3=E.Math.u.sub.3(G5)
so that,

(31) $\begin{array}{cc}\left(E\right).Math.\left(u_1{u}_2\right)d{q}_1{Q}_{1}d{q}_2{Q}_2=d{q}_1\frac{}{q_1}\left(E_2{\mathrm{dq}}_2{Q}_2\right)-{\mathrm{dq}}_2\frac{}{q_2}\left(E_{1}d{q}_1{Q}_1\right)& \left(\mathrm{G6}\right)\end{array}$
or,

(32) $\begin{array}{cc}\left(E\right).Math.\left(u_1{u}_2\right)Q_1{Q}_2=\frac{{\hat{E}}_2}{q_1}-\frac{{\hat{E}}_1}{q_2}={\left({}_q\hat{E}\right)}^3& \left(\mathrm{G7}\right)\end{array}$
where we use the conventional superscript notation for contravariant components of a vector. We have defined,
.sub.1=Q.sub.1E.sub.1,.sub.2=Q.sub.2E.sub.2,.sub.3=Q.sub.3E.sub.3(G8)
Note that the right-hand side of equation (10) is simple component 3 of curl evaluated in the new co-ordinate system. Now applying Maxwell,

(33) $\begin{array}{cc}\left(E\right).Math.\left(u_1{u}_2\right)Q_1{Q}_2=-{}_0\frac{H}{t}.Math.\left(u_1{u}_2\right)Q_1{Q}_2& \left(\mathrm{G9}\right)\end{array}$
We write H in terms of the contravariant components,
H=H.sup.1u.sub.1+H.sup.2u.sub.2+H.sup.3u.sub.3(G10)
which in turn can be expressed in terms of the covariant components,

(34) $\begin{array}{cc}{g}^{-1}\left[\begin{array}{c}{H}^1\\ H^2\\ H^3\end{array}\right]=\left[\begin{array}{ccc}{u}_1.Math.u_1& u_1.Math.u_2& u_1.Math.u_3\\ u_2.Math.u_1& u_2.Math.u_2& u_2.Math.u_3\\ u_3.Math.u_1& u_3.Math.u_2& u_3.Math.u_3\end{array}\right]\left[\begin{array}{c}{H}^1\\ H^2\\ H^3\end{array}\right]=\left[\begin{array}{c}{H}_1\\ H_2\\ H_3\end{array}\right]& \left(\mathrm{G11}\right)\end{array}$
where the first part defines g, and
H.sub.1=H.Math.u.sub.1,H.sub.2=H.Math.u.sub.2,H.sub.3=H.Math.u.sub.3(G12)
Inverting g gives
H.sup.i=.sub.j=1.sup.3g.sup.ijH.sub.j.(G13)
Substituting equation (13) into (12) gives

(35) $\begin{array}{cc}\left(E\right).Math.\left(u_1{u}_2\right)Q_1{Q}_2=-{}_0\frac{H}{t}.Math.\left(u_1{u}_2\right)Q_1{Q}_2=-{}_0\underset{j=1}{\overset{3}{.Math.}}{g}^{3j}\frac{H_j}{t}{u}_3.Math.\left(u_1{u}_2\right)Q_1{Q}_2& \left(\mathrm{G14}\right)\end{array}$
Define
{circumflex over ()}.sup.ij=g.sup.ij|u.sub.1.Math.(u.sub.2u.sub.3)|Q.sub.1Q.sub.2Q.sub.3(Q.sub.iQ.sub.j).sup.1(G15)
And
.sub.j=Q.sub.jH.sub.j(G16)
so that,

(36) $\begin{array}{cc}\left(E\right).Math.\left(u_1{u}_2\right)Q_1{Q}_2=-{}_0\underset{j=1}{\overset{3}{.Math.}}{\hat{}}^{3j}\frac{{\hat{H}}_j}{t}& \left(\mathrm{G17}\right)\end{array}$
Hence on substituting from equation (10),

(37) $\begin{array}{cc}{\left({}_q\hat{E}\right)}^i=-{}_0\underset{j=1}{\overset{3}{.Math.}}{\hat{}}^{\mathrm{ij}}\frac{{\hat{H}}_j}{t}& \left(\mathrm{G18}\right)\end{array}$
and by symmetry between E and H fields,

(38) 0$\begin{array}{cc}{\left({}_q\hat{H}\right)}^i=+{\mathrm{.Math.}}_0\underset{j=1}{\overset{3}{.Math.}}{\widehat{\mathrm{.Math.}}}^{\mathrm{ij}}\frac{{\hat{E}}_j}{t}& \left(\mathrm{G19}\right)\end{array}$
where
{circumflex over ()}.sup.ij=g.sup.ij|u.sub.1.Math.(u.sub.2u.sub.3)|Q.sub.1Q.sub.2Q.sub.3(Q.sub.iQ.sub.j).sup.1(G20)
Note that these expressions simplify considerable if the new co-ordinate system is orthogonal, e.g., cylindrical or spherical, when
g.sup.ij|u.sub.1.Math.(u.sub.2u.sub.3)|=.sub.ij(G21)
In the case of an orthogonal geometry (Cartesian, spherical or cylindrical) but the renormalized values of the permittivity and permeability are:

(39) $\begin{array}{cc}{\mathrm{.Math.}}_u^{}={\mathrm{.Math.}}_u\frac{Q_u{Q}_v{Q}_w}{Q_u^2},{}_u^{}={}_u\frac{Q_u{Q}_v{Q}_w}{Q_u^2},\mathrm{etc}.E_u^{}=Q_u{E}_u,H_u^{}=Q_u{H}_u,\mathrm{etc}.& \left(4\right)\\ \mathrm{where},& \left(5\right)\\ Q_u^2={\left(\frac{x}{u}\right)}^2+{\left(\frac{y}{u}\right)}^2+{\left(\frac{z}{u}\right)}^2{Q}_v^2={\left(\frac{x}{v}\right)}^2+{\left(\frac{y}{v}\right)}^2+{\left(\frac{z}{v}\right)}^2{Q}_w^2={\left(\frac{x}{w}\right)}^2+{\left(\frac{y}{w}\right)}^2+{\left(\frac{z}{w}\right)}^2& \left(6\right)\end{array}$
And:
B=.sub.0H and D=.sub.0E(7)

(40) Referring back to the spherical case as shown in FIG. 1, the natural choice of mesh on which to plot the original ray configuration is a spherical one, which uses the co-ordinates r, and . The co-ordinate transformation operates by taking all the fields in the sphere 106 of radius r<R.sub.2 and compressing them into the annulus 102 of radius R.sub.1<r<R.sub.2. This is achieved by any transformation of the form:

(41) $\begin{array}{cc}{r}^{}=R_1+\frac{f\left(r\right)-f\left(0\right)}{f\left(R_2\right)-f\left(0\right)}\left(R_2-R_1\right),{}^{}=,{}^{}=& \left(8\right)\end{array}$
where f(r) is a monotonically increasing function of r, this will map the entire interior of a sphere 106 of radius R.sub.2 into an annulus bounded by two spheres 100, 106 radii R.sub.1 and R.sub.2 respectively. There is symmetry for and because the field lines are being compressed radially only, so their angular positioning is not distorted.

(42) The simplest solution to equation (8) arises when f(r)=r, which simplifies to give:
r=R.sub.1+r(R.sub.2R.sub.1)/R.sub.2,
=,
=(9)

(43) In order to ascertain the electromagnetic properties which will be required of the cloaking annulus 102 material so that it achieves the desired result of cloaking any object in the sphere 100 R.sub.1, the transformation which has been imposed on r, and must be applied in an analogous manner to and , as shown for the general case above. At each point (x, y, z) in space, and will have component (.sub.x, .sub.y, .sub.z), (.sub.x, .sub.y, .sub.z) which must be transformed to spherical coordinate valve (.sub.r, .sub., .sub.), (.sub.r, .sub., .sub.) in the transformed geometry.

(44) We do this in two steps:

(45) First transform to spherical coordinates without any compression,
x=r cos sin ,y=r sin sin ,Z=r cos (A1)
so that, as shown in equation (4) above, in the new frame,

(46) $\begin{array}{cc}{\widetilde{\mathrm{.Math.}}}_i={\mathrm{.Math.}}_i\frac{Q_1{Q}_2{Q}_3}{Q_i^2},{\stackrel{~}{}}_i={}_i\frac{Q_1{Q}_2{Q}_3}{Q_i^2}& \left(A2\right)\end{array}$
We calculate,

(47) $\begin{array}{cc}{Q}_r^2={\left(\frac{x}{r}=\cos \sin \right)}^2+{\left(\frac{y}{r}=\sin \sin \right)}^2+{\left(\frac{z}{r}=\cos \right)}^2=1Q_{}^2={\left(\frac{x}{}=r\cos \cos \right)}^2+{\left(\frac{y}{}=r\sin \cos \right)}^2+{\left(\frac{z}{}=-r\sin \right)}^2=r^2{Q}_{}^2={\left(\frac{x}{}=-r\sin \sin \right)}^2+{\left(\frac{y}{Z}=+r\cos \sin \right)}^2+{\left(\frac{z}{Z}=0\right)}^2=r^2{\sin }^2& \left(\mathrm{A3}\right)\end{array}$
and hence,

(48) $\begin{array}{cc}{\widetilde{\mathrm{.Math.}}}_r={\stackrel{~}{}}_r=\frac{Q_r{Q}_{}{Q}_{}}{Q_r^2}{\mathrm{.Math.}}_r=\frac{r^2\sin }{1}{\mathrm{.Math.}}_r=r^2\sin {\mathrm{.Math.}}_r{\widetilde{\mathrm{.Math.}}}_{}={\stackrel{~}{}}_{}=\frac{Q_r{Q}_{}{Q}_{}}{Q_{}^2}{\mathrm{.Math.}}_{}=\frac{r^2\sin }{r^2}{\mathrm{.Math.}}_{}=\sin {\mathrm{.Math.}}_{}{\widetilde{\mathrm{.Math.}}}_{}={\stackrel{~}{}}_{}=\frac{Q_r{Q}_{}{Q}_{}}{Q_{}^2}{\mathrm{.Math.}}_{}=\frac{r^2\sin }{r^2{\sin }^2}{\mathrm{.Math.}}_{}=\frac{1}{\sin }{\mathrm{.Math.}}_{}& \left(\mathrm{A4}\right)\end{array}$
where .sub.r, .sub., .sub. the three components of the permittivity tensor in the original Cartesian frame, and we assume that the permittivity and permeability are equal. Note that we can easily extract the Cartesian view of .sub.r, .sub., .sub. from (A4) by removing the appropriate factors of r.

(49) In order to create a protected space inside the sphere radius R.sub.1 we make a further transformation to a new cylindrical system in which the radius is compressed taking any ray trajectories with it. It is noted that computer programs work in Cartesian coordinates and therefore that it is possible to re-express the compressed radial coordinate system in terms of a compressed Cartesian system, xyz, by removing factors of r.

(50) Consider the transformation from equation (9)
r=R.sub.1+r(R.sub.2R.sub.1)/R.sub.2, =, =(A5)
which can be re-expressed as,

(51) $\begin{array}{cc}r=\frac{\left(r^{}-R_1\right)R_2}{R_2-R_1},{}^{}=,{}^{}=& \left(\mathrm{A6}\right)\end{array}$
We calculate,

(52) $\begin{array}{cc}{Q}_{r^{}}^2={\left(\frac{r}{r^{}}=\frac{R_2}{R_2-R_1}\right)}^2+{\left(\frac{}{r^{}}=0\right)}^2+{\left(\frac{}{r^{}}=0\right)}^2={\left[\frac{R_2}{R_2-R_1}\right]}^2{Q}_{{}^{}}^2={\left(\frac{r}{{}^{}}=0\right)}^2+{\left(\frac{}{{}^{}}=1\right)}^2+{\left(\frac{}{{}^{}}=0\right)}^2=1Q_{{}^{}}^2={\left(\frac{r}{{}^{}}=0\right)}^2+{\left(\frac{}{{}^{}}=0\right)}^2+{\left(\frac{}{{}^{}}=1\right)}^2=1& \left(\mathrm{A7}\right)\end{array}$
so that in the new frame,
for R.sub.1<r<R.sub.2

(53) $\begin{array}{cc}\begin{array}{c}{\widetilde{\mathrm{.Math.}}}_{r^{}}={\stackrel{~}{}}_{r^{}}\\ ={\widetilde{\mathrm{.Math.}}}_r\frac{Q_{r^{}}{Q}_{{}^{}}{Q}_{{}^{}}}{Q_{r^{}}^2}\\ =\frac{\left(R_2-R_1\right)}{R_2}{r}^2\sin {\mathrm{.Math.}}_r\\ ={\frac{\left(R_2-R_1\right)}{R_2}\left[\frac{\left(r^{}-R_1\right)R_2}{R_2-R_1}\right]}^2\sin {\mathrm{.Math.}}_r\\ =\frac{{R_2\left(r^{}-R_1\right)}^2}{R_2-R_1}\sin {\mathrm{.Math.}}_r\end{array}{\widetilde{\mathrm{.Math.}}}_{{}^{}}={\stackrel{~}{}}_{{}^{}}={\widetilde{\mathrm{.Math.}}}_{}\frac{Q_{r^{}}{Q}_{{}^{}}{Q}_{{}^{}}}{Q_{{}^{}}^2}=\frac{R_2}{R_2-R_1}\sin {\mathrm{.Math.}}_{}{\widetilde{\mathrm{.Math.}}}_{{}^{}}={\stackrel{~}{}}_{{}^{}}={\widetilde{\mathrm{.Math.}}}_{}\frac{Q_{r^{}}{Q}_{{}^{}}{Q}_{{}^{}}}{Q_{{}^{}}^2}=\frac{R_2}{R_2-R_1}\frac{1}{\sin }{\mathrm{.Math.}}_{}& \left(\mathrm{A8a}\right)\end{array}$

(54) Alternatively we can reinterpret these values in a Cartesian coordinate frame, x y z, which is easily done by comparison with (A4),

(55) $\begin{array}{cc}{\mathrm{.Math.}}_{r^{}}={}_{r^{}}=\frac{R_2}{R_2-R_1}\frac{{\left(r^{}-R_1\right)}^2}{r^2}{\mathrm{.Math.}}_r{\mathrm{.Math.}}_{{}^{}}={}_{{}^{}}=\frac{R_2}{R_2-R_1}{\mathrm{.Math.}}_{}{\widetilde{\mathrm{.Math.}}}_{{}^{}}={\stackrel{~}{}}_{{}^{}}=\frac{R_2}{R_2-R_1}{\mathrm{.Math.}}_{}& \left(10\right)\end{array}$
where we have assumed that the starting material is vacuum.
for R.sub.2<r

(56) 0$\begin{array}{cc}{\widetilde{\mathrm{.Math.}}}_{r^{}}={\stackrel{~}{}}_{r^{}}=r^2\sin {}^{}{\mathrm{.Math.}}_{r^{}}{\widetilde{\mathrm{.Math.}}}_{{}^{}}={\stackrel{~}{}}_{{}^{}}=\sin {}^{}{\mathrm{.Math.}}_{{}^{}}{\widetilde{\mathrm{.Math.}}}_{Z^{}}={\stackrel{~}{}}_{Z^{}}=\frac{1}{\sin {}^{}}{\mathrm{.Math.}}_{{}^{}}& \left(\mathrm{A7b}\right)\end{array}$
In free space,
.sub.r=.sub.r=.sub.=.sub.=.sub.=.sub.=1(11)
for r<R.sub.1 we may choose any value we please for the dielectric functions since radiation never penetrates into this region. In this region , are free to take any value without restriction whilst still remaining invisible, and not contributing to electromagnetic scattering. This will be the case for all cloaked volumes, regardless of their shape or size.

(57) If a suitable cloaking material can be implemented for the annulus 102 of radius R.sub.1<r<R.sub.2, this will exclude all fields from entering the central region 100 of radius R.sub.1 and, conversely, will prevent all fields from escaping from this region. The cloaking structure itself and any object placed in the concealed sphere 100 R.sub.1 will therefore be undetectable in either reflection or transmission. It will be noted that at all points in the medium the impedance

(58) $z=\sqrt{\frac{}{\mathrm{.Math.}}}=1$
such that no unwanted reflection in the medium, caused by impedance mismatch, takes place.

(59) It is noted that if the impedance z={square root over (/)}=1 is achieved at all points in the metamaterial, the cloaking structure and concealed volume will appear to contain only free space, and so will be invisible if embedded in free space. If the cloaking structure is embedded in another medium the effect will differ, for example if embedded in water it would have the appearance of a bubble. In order to make the cloaking structure invisible in another medium it is necessary to match the ratio of and in the metamaterial to that of the surrounding medium at the interface between the two. As the skilled person will appreciate, whilst this changes the numerical values of and in the cloaking structure metamaterial, the general transformation theory outlined above remains the same. That is, because the present invention enables and to be controlled in a metamaterial independently of one another, it is possible to achieve any desired value of z, to match the impedance of the metamaterial to that of its surroundings.

(60) Whilst in FIG. 1 the cloaking annulus 102 is depicted as being relatively thick, this does not have to be the case. The solution is not dependent on shell thickness or any other length scales, such that it does not suffer from any scaling issues. Additionally, the solution provided holds in both in the near and far fields, because the solution does not depend on the wavelength . This is demonstrated in FIGS. 1 and 3. For purposes of illustration, suppose that R.sub.2>> where is the wavelength so that we can use the ray approximation to plot the Poynting vector. If our system is then exposed to a source of radiation at infinity we can perform the ray tracing exercise shown in FIG. 1. Rays 104 in this figure result from numerical integration of a set of Hamilton's equations obtained by taking the geometric limit of Maxwell's equations with anisotropic, inhomogeneous media. This integration provides an independent confirmation that the configuration specified by (9) and (10) excludes rays from the interior region 100. Alternatively if R.sub.2>> and we locate a point charge 306 nearby, the electrostatic (or magnetostatic) approximation applies. A plot of the local electrostatic displacement field which results is shown in FIG. 3.

(61) The theoretical solution outlined above can be demonstrated using a ray-tracing program to calculate the ray trajectories. According to such a program, a specified material is assumed to exist around a central volume that is to be cloaked, or rendered invisible to the user. Using the specified parameters as input, the ray-tracing code confirms the theoretical prediction. Ray tracing programs of this type will be well known to those skilled in the art, although some existing programs will not be able to deal with the complexity of the problem for this invention. In FIGS. 4a and 4b, a ray tracing program has been used which was developed from a custom ray tracing code based on a Hamiltonian formulation to solve Maxwell's equations in the geometric limit (i.e. zero wavelength limit). The results are shown for a spherical structure 400 and a cylindrical structure 402 respectively, wherein the transformed parameters E and in the cylindrical case are calculated in a similar manner to those for the sphere, as shown below:

(62) First transform to cylindrical coordinates without any compression,
x=r cos ,y=r sin ,Z=z(A8)
so that in the new frame,

(63) $\begin{array}{cc}{\widetilde{\mathrm{.Math.}}}_i={\mathrm{.Math.}}_i\frac{Q_1{Q}_2{Q}_3}{Q_i^2},{\stackrel{~}{}}_i={}_i\frac{Q_1{Q}_2{Q}_3}{Q_i^2}& \left(\mathrm{A9}\right)\end{array}$
We calculate,

(64) $\begin{array}{cc}{Q}_r^2={\left(\frac{x}{r}=\cos \right)}^2+{\left(\frac{y}{r}=\sin \right)}^2+{\left(\frac{z}{r}=0\right)}^2=1Q_{}^2={\left(\frac{x}{}=-r\sin \right)}^2+{\left(\frac{y}{}=+r\cos \right)}^2+{\left(\frac{z}{r}=0\right)}^2=r^2{Q}_Z^2={\left(\frac{x}{Z}=0\right)}^2+{\left(\frac{y}{Z}=0\right)}^2+{\left(\frac{z}{Z}=1\right)}^2=1& \left(\mathrm{A10}\right)\end{array}$
and hence,

(65) $\begin{array}{cc}{\widetilde{\mathrm{.Math.}}}_r={\stackrel{~}{}}_r=\frac{Q_r{Q}_{}{Q}_Z}{Q_r^2}{\mathrm{.Math.}}_r=r{\mathrm{.Math.}}_r{\widetilde{\mathrm{.Math.}}}_{}={\stackrel{~}{}}_{}=\frac{Q_r{Q}_{}{Q}_Z}{Q_{}^2}{\mathrm{.Math.}}_{}=r^{-1}{\mathrm{.Math.}}_{}{\widetilde{\mathrm{.Math.}}}_Z={\stackrel{~}{}}_Z=\frac{Q_r{Q}_{}{Q}_Z}{Q_Z^2}{\mathrm{.Math.}}_Z=r{\mathrm{.Math.}}_Z& \left(\mathrm{A11}\right)\end{array}$
where .sub.r, .sub., .sub.Z the three components of the permittivity tensor in the original Cartesian frame, and we assume that the permittivity and permeability are equal. Note that we can easily extract the Cartesian view of .sub.r, .sub., .sub.Z from (A13) by removing the appropriate factors of r.

(66) In order to create a protected space inside a cylinder radius R.sub.1 we make a further transformation to a new cylindrical system in which the radius is compressed taking any ray trajectories with it. It is noted that computer programs work in Cartesian coordinates and therefore that it is possible to re-express the compressed radial coordinate system in terms of a compressed Cartesian system, xy z, by removing factors of r.

(67) Consider the transformation,
r=R.sub.1+r(R.sub.2R.sub.1)/R.sub.2,=,Z=Z(A12)
or,

(68) $\begin{array}{cc}r=\frac{\left(r^{}-R_1\right)R_2}{R_2-R_1},{}^{}=,Z^{}=Z& \left(\mathrm{A13}\right)\end{array}$
We calculate,

(69) $\begin{array}{cc}{Q}_{r^{}}^2={\left(\frac{r}{r^{}}=\frac{R_2}{R_2-R_1}\right)}^2+{\left(\frac{}{r^{}}=0\right)}^2+{\left(\frac{Z}{r^{}}=0\right)}^2={\left[\frac{R_2}{R_2-R_1}\right]}^2{Q}_{{}^{}}^2={\left(\frac{r}{{}^{}}=0\right)}^2+{\left(\frac{}{{}^{}}=1\right)}^2+{\left(\frac{Z}{{}^{}}=0\right)}^2=1Q_{Z^{}}^2={\left(\frac{r}{Z^{}}=0\right)}^2+{\left(\frac{}{Z^{}}=0\right)}^2+{\left(\frac{Z}{Z^{}}=1\right)}^2=1& \left(\mathrm{A14}\right)\end{array}$
so that in the new frame,
for R.sub.1<r<R.sub.2

(70) $\begin{array}{cc}{\widetilde{\mathrm{.Math.}}}_{r^{}}={\stackrel{~}{}}_{r^{}}={\widetilde{\mathrm{.Math.}}}_r\frac{Q_{r^{}}{Q}_{{}^{}}{Q}_{Z^{}}}{Q_{r^{}}^2}=\frac{\left(R_2-R_1\right)}{R_2}r{\mathrm{.Math.}}_r=\frac{\left(R_2-R_1\right)}{R_2}\frac{\left(r^{}-R_1\right)R_2}{R_2-R_1}{\mathrm{.Math.}}_r=\left(r^{}-R_1\right){\mathrm{.Math.}}_r{\widetilde{\mathrm{.Math.}}}_{{}^{}}={\stackrel{~}{}}_{{}^{}}={\widetilde{\mathrm{.Math.}}}_{}\frac{Q_{r^{}}{Q}_{{}^{}}{Q}_{Z^{}}}{Q_{{}^{}}^2}=\frac{R_2}{R_2-R_1}{r}^{-1}{\mathrm{.Math.}}_{}=\frac{R_2}{R_2-R_1}\frac{R_2-R_1}{\left(r^{}-R_1\right)R_2}{\mathrm{.Math.}}_{}=\frac{1}{\left(r^{}-R_1\right)}{\mathrm{.Math.}}_{}& \left(\mathrm{A15a}\right)\\ {\widetilde{\mathrm{.Math.}}}_{Z^{}}={\stackrel{~}{}}_{Z^{}}={\widetilde{\mathrm{.Math.}}}_Z\frac{Q_{r^{}}{Q}_{{}^{}}{Q}_{Z^{}}}{Q_{Z^{}}^2}=\frac{R_2}{R_2-R_1}r{\mathrm{.Math.}}_Z=\frac{R_2}{R_2-R_1}\frac{\left(r^{}-R_1\right)R_2}{R_2-R_1}{\mathrm{.Math.}}_Z={\left[\frac{R_2}{R_2-R_1}\right]}^2\left(r^{}-R_1\right){\mathrm{.Math.}}_Z& \left(\mathrm{A15a}\right)\end{array}$
for R.sub.2<r
{tilde over ()}.sub.r={tilde over ()}.sub.r=r{tilde over ()}.sub.r
{tilde over ()}.sub.={tilde over ()}.sub.=r.sup.1.sub.
{tilde over ()}.sub.Z={tilde over ()}.sub.Z=r.sub.Z(A15b)
for r<R.sub.1 we may choose any value we please for the dielectric functions since radiation never penetrates into this region.
Alternatively we can reinterpret these values in a Cartesian coordinate frame, xyz, which is easily done by comparison with (A4),

(71) $\begin{array}{cc}{\mathrm{.Math.}}_{r^{}}={}_{r^{}}=\frac{r^{}-R_1}{r^{}}{\mathrm{.Math.}}_{{}^{}}={}_{{}^{}}=\frac{r^{}}{\left(r^{}-R_1\right)}{\mathrm{.Math.}}_{Z^{}}={}_{Z^{}}={\left[\frac{R_2}{R_2-R_1}\right]}^2\frac{r^{}-R_1}{r^{}}& \left(\mathrm{A16}\right)\end{array}$
where we have assumed that the starting material is vacuum.

(72) FIGS. 4a and 4b show that the incident rays 404, 406 are diverted around a central region 400, 402 and emerge from the cloaking region 408, 410 apparently unperturbed, as predicted by the theoretical solution.

(73) Whilst the theoretical solution has been confirmed by ray-tracing methods, in order to put the invention into effect a suitable material to form the cloaking structure has been developed.

(74) The material parameters of the cloaking structure require that the components of equal those of , that and vary throughout space, and that the tensor components of and vary independently. The anticipated material is thus inhomogeneous and anisotropic. Moreover, the material parameters and must assume values less than unity and approach zero at the interface between the concealed region and the cloaking shell. Furthermore, it is noted that equation (6) is singular at r=R.sub.1. This is unavoidable as can be seen by considering a ray headed directly towards the center of the sphere 100 (FIG. 1). This ray does not know whether to be deviated up or down, left or right. Neighboring rays 104 are bent around tighter and tighter arcs closer to the critical ray they are. This in turn implies very steep gradients in and . Parameters and are necessarily anisotropic in the cloaking material because space has been compressed anisotropically. This set of constraints is not attainable by conventional materials for example because finding matching values for and in the same waveband does not occur. However, given the rapid progress that has occurred over the past several years in artificial materials, a metamaterial can be practically designed that satisfies the specifications.

(75) Metamaterials are artificially constructed materials which can exhibit electromagnetic characteristics that are difficult or impossible to achieve with conventional materials. From an electromagnetic point of view, the wavelength, , of a wave passing through a material determines whether a collection of atoms or other objects can be considered a material and properties at the atomic level determine and . However, the electromagnetic parameters and need not arise strictly from the response of atoms or molecules: Any collection of objects whose size and spacing are much smaller than can be described by an and . In that case, the values of and are determined by the scattering properties of the structured objects. Although such an inhomogeneous collection may not satisfy an intuitive definition of a material, an electromagnetic wave passing through the structure cannot tell the difference, and, from the electromagnetic point of view, we have created an artificial material, or metamaterial. These properties of metamaterials are explained further in Metamaterials and Refractive Index, D R Smith, J B Pendry, MCK. Wiltshire, VOL 305, SCIENCE 061081 (2004), which is incorporated herein by reference.

(76) One useful property exhibited by certain metamaterials is that of artificial magnetism. It has been found that structures consisting of non-magnetic arrays of wire loops in which an external field can induce a current consequently produce an effective magnetic response. One such design is the so-called Swiss-roll structure, as shown in FIG. 5. The roll 500 is manufactured by rolling an insulated metallic sheet around a cylinder 504. A design with about 11 turns on a 1-cm-diameter cylinder gives a resonant response at 21 MHz. The metamaterial is formed by stacking together many of these cylinders. In this structure, the coiled copper sheets have a self-capacitance and self-inductance that create a resonance. The currents that flow when this resonance is activated couple strongly to an applied magnetic field, yielding an effective permeability that can reach quite high values. No current can flow around the coil except by self capacitance. As is shown in FIG. 6, when a magnetic field parallel to the cylinder 600 is switched on it induces currents (j) in the coiled sheets 602, which are spaced a distanced apart from one another. Capacitance between the first 604 and last 606 turns of the coil enables current to flow. This is described further in Magnetism from Conductors and Enhanced Non-Linear Phenomena, J B Pendry, A J Holden, D J Robbins and W J Stewart IEEE Trans. Micr. Theory and Techniques, 47, 2075 (1999), which is incorporated herein by reference.

(77) One problem with the Swiss-roll structures is that the continuous electrical path provided by metal cylinders can cause the metamaterial to respond like an effective metal when the electric field applied is not parallel to the cylinders, hence restricting its usefulness in certain applications. An adaptation described in the same paper which avoids this undesirable effect is the split ring resonator (SRR). The SRR is built up from a series of split rings, as shown in FIG. 7. Each split ring 700 comprises at least two concentric thin metal rings 702, 704 of width c, spaced apart by a distance d, each one with a gap, 706, 708. By eliminating the continuous conducting path which the cylinders provide, an SRR eliminates most of the electrical activity along this direction.

(78) As shown in FIG. 8, it is possible to form a planar split ring structure 800 from a series of split rings 802, spaced a distance a apart. The two dimensional square array 800 shown in FIG. 8 can be made by printing with metallic inks. If each printed sheet is then fixed to a solid block of inert material, thickness a, the blocks can be stacked to give columns of rings. This would establish magnetic activity along the direction of stacking, z-axis.

(79) The unit cell 900 of the stacked SRR structure is shown in FIG. 9 on the left. It is possible to form a symmetrical three dimensional structure starting from such a structure comprising successive layers of rings stacked along the z-axis. This is achievable by cutting up the structure into a series of slabs thickness a, making incisions in the y-z plane and being careful to avoid slicing through any of the rings. Each of the new slabs contains a layer of rings but now each ring is perpendicular to the plane of the slab and is embedded within. The next stage is to print onto the surface of each slab another layer of rings and stack the slabs back together again. The unit cell 902 of this second structure is shown in the middle of FIG. 9. In the final step, a third set of slabs is produced by cutting in the x-z plane, printing on the surface of the slabs, and reassembling. This results in a structure with cubic symmetry whose unit cell 904 is shown on the right of FIG. 9.

(80) Similarly c can be governed by an array of thin wires as described in Extremely Low Frequency Plasmas in Metallic Mesostructures, J B Pendry, A J Holden, W J Stewart, I Youngs, Phys Rev Lett, 76, 4773 (1996) which is incorporated herein by reference.

(81) Hence known techniques can be applied to construct the metamaterial required to put the current invention as described above into practice. The materials comprise many unit cells, each of which is designed to have a predetermined value for and , such that they correspond to the values determined by equation (4) above (or for the non-orthogonal case, equation G23 above. As will be apparent to the skilled reader, this may involve each unit cell having slightly different values of E and to those of its radially neighboring cells. Each unit cell in the metamaterial acts as a concealing volume element. The concealing volume elements are assembled in the correct respective spatial positions to form the cloaking structure, as calculated by the methods outlined above. The cloaking structure will, as a result, possess the spatially distributed material parameter values required for electromagnetic concealment of the enclosed volume.

(82) For practical implementation of the theory outlined above, a three dimensional cloaking structure will be necessary in most cases. A solid composite metamaterial structure can be formed using fabrication techniques similar to those used to manufacture a gradient index lens. One such known structure 1000 is shown in FIG. 10, in which SRR's and wires are deposited on opposite sides lithographically on a standard circuit board to provide a desired material response. This is described in more detail in Reversing Light with Negative Refraction, JB Pendry, DR Smith, Physics Today, June 2004, which is incorporated herein by reference.

(83) The unit cell structure of the metamaterial used ensures that anisotropy and continuous variation of the parameters, as required by the theoretical solution, can be achieved. Unrelated to cloaking, another application of the transformation theory described above is a device for concentrating fields. FIG. 11 shows a concentrator that concentrates electromagnetic fields incident upon its outer surface (110) onto an inner core region (114). The transformation that characterizes the concentrator, for the spherical geometry, is shown in FIG. 11. The outer radius of the device (110) is R.sub.3. A concentrating volume of space (112) with radius R.sub.2<R.sub.3 is compressed to a collecting volume (114) of radius R.sub.1. The shell volume (116) that lies between R.sub.2 and R.sub.3 is expanded to fill the region (118) between R.sub.1 and R.sub.3, as FIG. 11 depicts. A simple example of this transformation that compresses the core (112) uniformly and stretches the shell region (116) with a linear radial function is given by

(84) $\begin{array}{cc}r=\{\begin{array}{cc}\frac{R_2}{R_1}{r}^{}& 0<r^{}<R_1\\ \frac{R_2-R_1}{R_3-R_1}{R}_3+\frac{R_3-R_2}{R_3-R_1}{r}^{}& R_1<r^{}<R_3\\ r^{}& R_3<r^{}\end{array}& \left(12\right)\\ ={}^{}={}^{}& \left(13\right)\end{array}$
The resulting material properties of permitting and permeability in the concentration device are then

(85) 0$\begin{array}{cc}\left[{\mathrm{.Math.}}_{r^{}}^{},{\mathrm{.Math.}}_{{}^{}}^{},{\mathrm{.Math.}}_{{}^{}}^{}\right]=\left[{}_{r^{}}^{},{}_{{}^{}}^{},{}_{{}^{}}^{}\right]=\{\begin{array}{cc}\frac{R_2}{R_1}\left[1,1,1\right]& 0<r^{}<R_1\\ \frac{R_3-R_2}{R_3-R_1}\left[{\left(1+\frac{R_2-R_1}{R_3-R_2}\frac{R_3}{r^{}}\right)}^2,1,1\right]& R_1<r^{}<R_3\\ \left[1,1,1\right]& R_3<r^{}\end{array}& \left(15\right)\end{array}$
where the material properties are given relative to the external environment. As with the cloaking structure, the external environment may be free space or any other medium.

(86) In order to concentrate fields into an arbitrarily small volume, we let R.sub.1.fwdarw.0 and obtain

(87) $\begin{array}{cc}\left[{\mathrm{.Math.}}_{r^{}}^{},{\mathrm{.Math.}}_{{}^{}}^{},{\mathrm{.Math.}}_{{}^{}}^{}\right]=\left[{}_{r^{}}^{},{}_{{}^{}}^{},{}_{{}^{}}^{}\right]=\{\begin{array}{cc}\frac{R_3-R_2}{R_3}\left[{\left(1+\frac{R_3}{R_3-R_2}\frac{R_2}{r^{}}\right)}^3,1,1\right]& 0<r^{}<R_3\\ \left[1,1,1\right]& R_3<r^{}\end{array}& \left(16\right)\end{array}$
If we wish the concentration cross section of the sphere to equal its outer radius, then we let R.sub.2.fwdarw.R.sub.3.

(88) $\begin{array}{cc}{\lim }_{R_2.fwdarw.R_3}\left[{\mathrm{.Math.}}_{r^{}}^{},{\mathrm{.Math.}}_{{}^{}}^{},{\mathrm{.Math.}}_{{}^{}}^{}\right]{\lim }_{R_2.fwdarw.R_3}\left[{}_{r^{}}^{},{}_{{}^{}}^{},{}_{{}^{}}^{}\right]=\{\begin{array}{cc}\left[+,0,0\right]& 0<r^{}<R_3\\ \left[1,1,1\right]& R_3<r^{}\end{array}& \left(17\right)\end{array}$

(89) This is the perfect field concentrator. The material parameters which are required in the concentrator can be achieved through the use of metamaterials as described in relation to the cloaking structure. The permittivity and permeability can be controlled independently of one another in each of three dimensions within each unit cell of the metamaterial.

(90) The concentrator device is unique, and different from focusing lenses and mirrors, in that it operates on fields incident from any direction. For electromagnetic sensing or energy harvesting applications, a detector or collector would be placed at the core (114) to interact with the intensified fields at that location. Part of the utility of this device is concentration of fields when the direction of their source is initially unknown or constantly changing e.g. solar energy collection or microwave energy beaming on mobile platforms.

(91) A detector placed at the concentrator's core (114) could easily be configured to provide accurate directional information as well as collect the energy. This could be accomplished by, for example, using a spherical shaped detector patterned into eight quadrants. By comparing the energy collected from each quadrant, accurate directional information will be obtained. The detector need not be spherical. A transformation could be used that results in a core region that is a thin circular disk. This would make the concentrator compatible with common planar detectors.

(92) The device can also provide impedance matching to the detector. This is achieved by applying the inverse transformation to the desired detector, expanding it from the core radius, R.sub.1, to the radius R.sub.2. It is then possible to design an impedance matching layer, (between R.sub.1 and R.sub.3), to match the transformed detector to its environment. This impedance matching layer is then transformed with the forward transformation expanding it to fill the region (118) between R.sub.1 and R.sub.3. The resulting shell layer between R.sub.1 and R.sub.3 provides both concentration and impedance matching functions.

(93) Like the cloak structures the concentrator (110) can be perfectly matched to its environment, so that if the fields are not disturbed in the core (114) by a detector, the device (110) is undetectable by external electromagnetic probing. This device (110) can operate in both near and far field regimes and at any size scale and wavelength. The basic formalism of this transformation can be applied to other basic geometries, (e.g. cylindrical) or applied conformally to an arbitrarily shaped object. FIG. 12 demonstrates the concentration of rays or field lines from a point source (120) incident upon a sphere (122) radius R.sub.3 into an inner core region (124) radius R.sub.1.

(94) The co-ordinate transfer mapping theory outlined above has been described with respect to diverting the propagation of electromagnetic waves. However it is noted that there is a significant analogy between electromagnetic and other propagating waves such as acoustic waves. In particular, this theory may be applied to obtain specifications for the mechanical properties of a material such as density and stiffness that yield acoustic cloaking. Such a cloak may be used for example to conceal bodies from sonar detection.

(95) Whilst the description of the cloaking theory above is directed mainly to concealing objects, it will be appreciated that the diversion of electromagnetic disturbances from the cloaked volume into a cloaking shell creates an electromagnetic quiet zone inside. Equation (9) dictates that any electromagnetic radiation produced by an object inside the quiet zone will not be able to emerge from it. This may be put to use in electromagnetic interference shielding applications, for example inside Mill machines. Analogously, acoustic quiet zones may also be achieved.

(96) The application of the present invention is not limited to the cloaking and concentrator structures as shown in the Figures. As demonstrated clearly in above, the mathematical solution provided herein applies to cloak a structure of any shape. Furthermore, if a volume is cloaked, an object of any shape or material placed therein will not be detectable by an outside observer, either in reflection mode or transmission mode.

(97) The methods described herein can be used to apply any appropriate coordinate transformation between any appropriate co-ordinate systems and using any appropriate function of the co-ordinates in the original frame. The cloaking structure used to conceal the material within it can be of any shape or thickness. In particular, the cloaking structure does not have to follow the same shape as the object within that is to be concealed. It will be apparent to the skilled reader that the methods described herein can be used to transform all frequencies of electromagnetic radiation, and for a given application can be applied to radiation in any specific frequency range of the electromagnetic spectrum, including infrared, radio, visible light, x-rays, microwaves and radar.

(98) In the examples shown impedance matching has been achieved by the ratio

(99) $\frac{}{\mathrm{.Math.}}=1$
As will be apparent to the skilled reader, impedance matching of the sort described herein can be achieved for

(100) $\frac{}{\mathrm{.Math.}}=m$
where m is any real constant value. In particular, the desired ratio of / will depend on the medium in which the cloaking structure is to be embedded.

(101) The reference J. B. Pendry et al., Magnetism from Conductors and Enhanced Nonlinear Phenomena, IEEE Trans. Micro. Theory Tech. 47 (1999), 2074-2084, above incorporated by reference, includes the following text:

(102) C Swiss Roll Capacitor

(103) In this instance, we find for the effective permeability

(104) $\begin{array}{cc}{}_{\mathrm{eff}}=1-\frac{F}{1+\frac{2i}{r{}_0\left(N-1\right)}-\frac{1}{2{}^2{r}^3{{}_0\left(N-1\right)}^2{}^{2}c}}=1-\frac{\frac{r^2}{a^2}}{1+\frac{2i}{r{}_0\left(N-1\right)}-\frac{{\mathrm{dc}}_0^2}{2{}^2{r}^3{}_0\left(N-1\right){}^2}}& \left(29\right)\end{array}$
where F is as before the fraction of the structure not internal to a cylinder, and the capacitance per unit area between the first and the last of the coils is

(105) $\begin{array}{cc}C=\frac{{\mathrm{.Math.}}_0}{d\left(N-1\right)}=\frac{1}{{}_0{\mathrm{dc}}_0^2\left(N-1\right)}& \left(30\right)\end{array}$

(106) IV. An Isotropic Magnetic Material

(107) . . . . We propose an adaptation of the split ring structure, in which the cylinder is replaced by a series of flat disks each of which retains the split ring configuration, but in slightly modified form . . . .

(108) The effective magnetic permeability we calculate, on the assumption that the rings are sufficiently close together and that the magnetic lines of force are due to currents in the stacked rings, are essentially the same as those in a continuous cylinder. This can only be true if the radius of the rings is of the same order as the unit cell side. We arrive at

(109) ${}_{\mathrm{eff}}=1-\frac{\frac{r^2}{a^2}}{1+\frac{2l{}_1}{r{}_0}i-\frac{3l}{{}^2{}_0{}^2{C}_1{r}^3}}=1-\frac{\frac{r^2}{a^2}}{1+\frac{2l{}_1}{r{}_0}i-\frac{3{\mathrm{lc}}_0^2}{{}^2\ln \frac{2c}{d}{r}^3}}$
where .sub.1 is the resistance of unit length of the sheets measured around the circumference.

(110) The reference J. B. Pendry et al., Extremely Low Frequency Plasmons in Metallic Mesostructures, Phys. Rev. Lett 76 (1996), 4773-4776, above incorporated by reference, includes the following text:

(111) The plasmons have a profound impact on properties of metals, not least on their interaction with electromagnetic radiation where the plasmon produces a dielectric function of the form

(112) $\begin{array}{cc}\mathrm{.Math.}\left(\right)=1-\frac{{}_p^2}{\left(+i\right)}& \left(2\right)\end{array}$
which is approximately independent of wave vector, and the parameter is a damping term representing dissipation of the plasmon's energy into the system . . . .

(113) In this Letter we show how to manufacture an artificial material in which the effective plasma frequency is depressed by up to 6 orders of magnitude. The building blocks of our new material are very thin metallic wires . . . .

(114) Having both the effective density, n.sub.eff, and the effective mass, m.sub.eff, on hand we can substitute into (1),

(115) $\begin{array}{cc}{}_p^2=\frac{n_{\mathrm{eff}}{e}^2}{{\mathrm{.Math.}}_0{m}_{\mathrm{eff}}}=\frac{2c_0^2}{a^2\ln \left(\frac{a}{r}\right)}{\left(8.2\mathrm{GHz}\right)}^2& \left(14\right)\end{array}$
Here is the reduction in the plasma frequency promised.

(116) Note in passing that, although the new reduced plasma frequency can be expressed in terms of electron effective mass and charge, these microscopy quantities cancel, leaving a formula containing only macroscopic parameters of the system: wire radius and lattice spacing.

(117) While this invention has been particularly shown and described with references to preferred embodiments thereof, it will be understood by those skilled in the art that various changes in form and details may be made therein without departing from the scope of the invention encompassed by the appended claims