Abstract

We analyze cloaking of transverse electric (TE) fields through homogenization of radially symmetric metallic structures. The two-dimensional circular cloak consists of concentric layers cut into a large number of small infinitely conducting sectors which is equivalent to a highly anisotropic permittivity. We find that a wave radiated by a magnetic line current source located a couple of wavelengths away from the cloak is almost unperturbed in magnitude but not in phase. Our structured cloak is shown to work for different wavelengths provided they are ten times larger than the outermost sectors.

© 2008 Optical Society of America

1. Introduction

Over a year ago, Pendry, Schurig and Smith theorized that a finite size object surrounded by a coating consisting of a meta-material might become invisible for electromagnetic waves [1]. The cornerstone of this work is the geometric transformation

r=R1+r(R2R1)R2,0rR2,
θ=θ,0<θ2π,
x3=x3,x3IR,

where r′, θ′, x3 are radially contracted cylindrical coordinates r, θ, x 3. The former transformation maps the disk DR2 with radius R 2 onto an annulus with outer radius R 2 and inner radius R 1. In other words, if a source located outside the disk DR2 radiates in vacuum, the electromagnetic field cannot reach the disk DR1 .

In terms of electromagnetic parameters, one should replace the material in the annulus by an equivalent one that is inhomogeneous and anisotropic [2]. The diagonalized form of the permittivity tensor (ε=) and the permeability tensor (μ=) was given by Pendry, Schurig and Smith [1]

εr=μr=rR1r,εθ=μθ=rrR1,ε3=μ3=(R2R2R1)2rR1r.

Note that there is no change in the impedance of the media since the permittivity and permeability undergo the same transformation.

An international team involving these authors subsequently implemented this idea using a meta-material consisting of concentric layers of SRR [3], which makes a copper cylinder invisible to an incident plane wave at a specific microwave frequency (8.5 GHz). The smooth behaviour of the electromagnetic field in the far field limit could be expected in view of the numerical evidence provided in [1] using a geometrical optics based software. This is also in agreement with the work of Leonhardt who independently studied conformal invisibility devices using the stationary Schrödinger equation [4]. To date, the only evidence that invisibility is preserved in the intense near field limit is purely numerical [5].

A very different route to invisibility is proposed by McPhedran, Nicorovici and Milton who studied a countable set of line sources using anomalous resonance when it lies in the close neighborhood of a cylindrical coating filled with negative permittivity material, which is nothing but a cylindrical version of the celebrated perfect lens of Sir John Pendry [6, 7, 8, 9, 10]. The former researchers attribute this cloaking phenomenon to anomalous localized resonances [11, 12].

Last, but not least, a team led by Shalaev [13] has recently shown the possibility to make an object nearly invisible in TE polarization using radially symmetric locally resonant inclusions at optical frequencies. The key point in their proposal is that when the magnetic field is polarized along the x 3-axis only two entries of ε¯¯ must satisfy the requirements in equation (2) (εr and εθ). Moreover µ 3 is the only entry of μ¯¯ which is involved in equation (2) and it can therefore be normalized to 1. The reduced set of material parameters is obtained by multiplying εr and εθ by µ 3 in (2) [13]:

εr=(R2R2R1)2(rR1r)2,εθ=(R2R2R1)2,μ3=1.

The electromagnetic parameters of the cloak (2) and (3) provide the same wave trajectory, but the latter one leads to impedance mismatch on the outer boundary r=R 2 of the cloak. Importantly, we note that the value of εθ increases with decreasing values of r in the annulus until it eventually becomes infinite on the inner boundary r=R 1 of the cloak in (2). The cloak proposed by Cai et al. does not exhibit a very large azimuthal anisotropy as it is clear that all parameters in (3) remain bounded (unlike the cloak proposed by Pendry et al.). To meet the criterion of (3), Cai et al. propose to use a locally resonant structure consisting of concentric layers of thin elongated ellipsoidal wires elongated along the radial direction. A derivation based on effective medium theory shows that such a cloak displays the prerequisite electromagnetic parameters around a given frequency (in the visible spectrum) at which the local resonators are excited by the TE wave.

In the current paper, we propose an alternative route to invisibility taking advantage of designs studied in [3] and [13]: we model a cloak consisting of concentric layers of increasing thickness, each of them being cut in a large number of small sectors. These inclusions are assumed to be infinitely conducting, which is an accurate model for waves ranging from Gigahertz to Terahertz frequencies.

2. Homogenization of the micro-structured cloak at fixed frequency

A rigorous and detailed derivation of the homogenized permittivity tensor is beyond the scope of this paper and involves subtle mathematics from the limit analysis domain. Thus, we will give here only the main homogenization results and the proof will be published elsewhere.

Homogenization theory predicts that this kind of ‘approximate invisibility’ is applicable for any TE incident wave (possibly out-of-plane) whose wavelength λ is large compared with the characteristic size d of sectors i.e. when η=λ/d≪1, unlike for the proposals in [3] (involving C-shaped Split ring Resonators) and [13] (involving elongated ellipsoidal wires). In the numerical experiments, we will see that the shadow region of the ‘F’-shaped object is already noticeably reduced when the parameter η~1/3. We notice that our micro-structured cloak which works through field averaging does not suffer from any blow-up of the electromagnetic field, unlike its locally resonant counterpart [14].

Throughout this paper, we work in cylindrical coordinates (r, θ, x 3), thereby taking into account the axi-symmetric geometry of the structure. From a two-scale model of the problem, assuming that the cloak lies in vacuum, we found that the homogenized permittivity tensor of the proposed structure is given by:

εhom¯¯=1area(Y*)(area(Y*)000area(Y*)000)(ϕrrϕ0ϕθrϕθθ0000),

where Y* denotes the elementary area around each scatterer B and ϕij represent corrective terms defined by:

i,j{r,θ},ϕij=<Vjyi>Y*=<Viyj>Y*=<Vi·Vj>Y*,

the brackets denoting averaging over Y*. Hence, thanks to the symmetry of the right matrix above (ϕij=ϕji), the homogenized permittivity is given by the knowledge of three terms ϕij. One can note that Vj is the unique Y-periodic solution with null mean of the following system, where derivatives are taken in the usual sense:

𝒦j:ΔVj=0,inY*=YB,andVjn=nj,onB,

where ∂B denotes the boundary of B, and nj, j ∈ {r,θ}, denotes the projection on the axis e j of a unit outward normal to ∂B. Although these results might at first glance look similar to those reported in [15], we emphasize that they differ since the earlier paper considered merely dielectric structures.

3. Numerical analysis of electromagnetic cloaking

In this section, we give some numerical illustrative examples of cloaking of an electric field coming from a source located closeby the obstacle through the infinitely conducting cloak. We replace the structured cloak by a cylinder surrounded by an effective coating whose homogenized permittivity is deduced from the numerical solution of the annex problems given by Eq. (6). This provides us with a qualitative picture of the cloaking mechanism (see Fig. 1). We then compare this asymptotic theory against the numerical solutions of the same scattering problems when we model the complete structured cloak.

3.1. The homogenized matrix

Using finite elements, we numerically solved the annex problems 𝓚r and 𝓚θ, for the geometry shown in Fig. 1 which provided us with two electrostatic potentials Vr and Vθ. Using (5) we found that ϕ and ϕθr vanish. From ϕθθ and ϕrr we deduced the following homogenized permittivity εhom¯¯=diag(1.7,8.2,) which displays a strong azimuthal anisotropy and moreover an infinite longitudinal anisotropy.

The corresponding scattering problem for a magnetic line current source is reported on Fig. 1 right. We notice that although very few back reflection takes place on the cloak, the shadow region behind the F-shaped obstacle has been nearly removed (compare with Fig. 3). This clearly demonstrates that the field follows an optical path surrounding the central region, thus making the object located inside invisible.

 

Fig. 1. Left: Geometry of the structured infinitely conducting cloak. The structure consists of 256 angular sectors. Right: Diffraction of the field radiated by a magnetic wire source by an infinitely conducting F-shaped obstacle surrounded by a homogeneous dielectric anisotropic coating with effective permittivity εhom¯¯=diag(1.7,8.2,) and effective transmission conditions n(εhom¯¯1(+iγ)Hhom)|r=n(+iγ)Hhom|r+ on its inner and outer boundaries (r=R 1=0.144 and r=R 2=0.4). The frequency ν is equal to 3.5.

Download Full Size | PPT Slide | PDF

3.2. Structured infinitely conducting cloak

On Figure 2, we display the real part ℜe(H 3) of the longitudinal component H 3 of the magnetic field along the x 2-axis. The black curve represents ℜe(H 3) in free space; the red curve represents ℜe(H 3) with a F-shaped infinitely conducting obstacle; the blue curve represents ℜe(H 3) with a F-shaped infinitely conducting obstacle surrounded by the cloak. We notice that the black and blue curves have nearly the same amplitude outside the cloak, but experience a phase shift. This phase shift can be attributed to the change in the optical path followed by the waves around the cloak: the path is longer than that of rays going straightforwardly in free space.

 

Fig. 2. Real part ℜe(H 3) of the longitudinal component H 3 of the magnetic field along the x 2-axis. The origin of the x 2 axis is taken at point (4,4) of Fig. 3 and it ends at point (0,0). A magnetic line source of wavelength λ=c/ν=c/3.5 is located at point (2.3,2.3). The vertical thick bold lines represent the outer boundary of the cloak.

Download Full Size | PPT Slide | PDF

Figure 3 shows maps of the real part ℜe(H 3) of the longitudinal component H 3 of the magnetic field when the F-shaped scatterer is coated (right) and uncoated (left) and for frequencies ν ranging from 3.5 to 5.5. Note that the modelling takes into account the actual structure of the infinitely conducting cloak. The diffraction is clearly reduced by the structure especially for the lowest frequencies. Moreover, the numerical results show that although not perfect the cloak allows us to reduce the scattering for frequencies in a relatively large domain thanks to the fact that the physical principle is not based on any resonant property of the structure. Of course, when the frequency increases the typical size of the sectors becomes significant when compared to the wavelength and thus the homogenized model is not longer valid.

Importantly, we observe a strong penetration of the field within the annular shell, similarly to recently reported simulations for annular cloaks, but in contrast with numerical simulations reported in [16] for a ferro-magnetic cloaking device whereby the electric field was mostly concentrated within the first concentric outer-ring.

These diagrams numerically show that the propagation of the longitudinal magnetic field H 3 (and hence the transverse electric field) is controlled by the cloak.

4. Conclusion

Invisibility was numerically studied using modelling of both the complete metallic structured device displaying around two hundred sub-wavelength sectors and its strongly anisotropic homogenized counterpart.

We numerically checked that a structured device with a large number of thin metallic sectors is already a reasonable cloak. Almost no back-scattering is observed, the field is well reconstructed behind the cloak with a much reduced shadow zone.

All our results are in dimensionless physical units, but for a micro-wave design, it is enough to take the size of the sector d to be of the order of magnitude of a millimeter, in which case results should hold for TE micro-waves propagating out-of-plane provided their propagation constant γ is such that γd ≪1. This suggests potential applications lie in improved telecommunication lines, for instance to transfer secure information. Interestingly, our cloak works well for the near field, and it is broad band to certain extent. Our design might prove useful for instance in certain problems of electromagnetic noise insulation for antennas applications. The influence of the absorption remains to be evaluated if the visible spectrum is aimed at.

 

Fig. 3. 2D plot of the real part ℜe(H 3) of the longitudinal component H 3 of the magnetic field radiated by a harmonic line current source of wavelength λ=c/ν. Left panel: Diffraction by an infinitely conducting F-shaped obstacle; Right panel: Diffraction by an infinitely conducting F-shaped obstacle surrounded by the structured cloak. When the frequency ν increases, the diffraction worsens and cloaking becomes less effective.

Download Full Size | PPT Slide | PDF

Acknowledgments

ABM and SG acknowledge funding from EPSRC under grant EP/F027125/1.

References and links

1. J.B. Pendry, D. Shurig, and D.R. Smith, “Controlling electromagnetic fields,” Science 3121780–1782 (2006). [CrossRef]   [PubMed]  

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

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

4. U. Leonhardt, “Optical conformal mapping,” Science 3121777–1780 (2006). [CrossRef]   [PubMed]  

5. F. Zolla, S. Guenneau, A. Nicolet, and J.B. Pendry, “Electromagnetic analysis of cylindrical invisibility cloaks and the mirage effect,” Opt. Lett. 32, 1069–1071 (2007). [CrossRef]   [PubMed]  

6. R.C. McPhedran, N.A. Nicorovici, and G.W. Milton,“Optical and dielectric properties of partially resonant composites,” Phys. Rev. B 49, 8479–8482 (1994). [CrossRef]  

7. J.B. Pendry, “Negative refraction makes a perfect lens,” Phys. Rev. Lett. 85, 3966–3969 (2000). [CrossRef]   [PubMed]  

8. J.B. Pendry and S.A. Ramakrishna, “Focussing light using negative refraction,” J. Phys. Cond. Matter 15, 6345–6364 (2003). [CrossRef]  

9. D. Maystre and S. Enoch, “Perfect lenses with left-handed material: Alice’s mirror?,” J. Opt. Soc. Am. A 21, 122 (2004). [CrossRef]  

10. S.A. Ramakrishna, “Physics of negative refractive index materials,” Rep. Prog. Phys. 68, 449–521 (2005). [CrossRef]  

11. G.W. Milton and N.A. Nicorovici, “On the cloaking effects associated with anomalous localised resonance,” Proc. Roy. Lond. A 462, 3027–3059 (2006). [CrossRef]  

12. N.A.P. Nicorovici, G.W. Milton, R.C. McPhedran, and L.C. Botten, “Quasistatic cloaking of two-dimensional polarizable discrete systems by anomalous resonance,” Opt. Express 15, 6314–6323 (2007). [CrossRef]   [PubMed]  

13. W. Cai, U.K. Chettiar, A.V. Kildiev, and V.M. Shalaev, “Optical Cloaking with metamaterials,” Nature 1, 224–227 (2007).

14. A. Greenleaf, Y. Kurylev, M. Lassas, and G. Uhlmann, “Improvement of cylindrical cloaking with the SHS lining,” Opt. Express 15, 12717–12734 (2007). [CrossRef]   [PubMed]  

15. S. Guenneau and F. Zolla, “Homogenization of three-dimensional finite photonic crystals,” JEWA14, 529–530 (2000) & Progress In Electromagnetics Research 27, 91–127 (2000).

16. D.P. Gaillot, C. Croenne, and D. Lippens, “An all-dielectric route for terahertz cloaking,” Opt. Express 16, 3986–3992 (2008). [CrossRef]   [PubMed]  

References

  • View by:
  • |
  • |

  1. J. B. Pendry, D. Shurig, D. R. Smith, "Controlling electromagnetic fields," Science 312, 1780-1782 (2006).
    [CrossRef] [PubMed]
  2. U. Leonhardt and T. G. Philbin, "General relativity in electrical engineering," New J. Phys.  8, 247 (2006).
    [CrossRef]
  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. U. Leonhardt, "Optical conformal mapping," Science 312, 1777-1780 (2006).
    [CrossRef] [PubMed]
  5. F. Zolla, S. Guenneau, A. Nicolet, and J. B. Pendry, "Electromagnetic analysis of cylindrical invisibility cloaks and the mirage effect," Opt. Lett. 32, 1069-1071 (2007).
    [CrossRef] [PubMed]
  6. R. C. McPhedran, N. A. Nicorovici, and G. W. Milton,"Optical and dielectric properties of partially resonant composites," Phys. Rev. B 49, 8479-8482 (1994).
    [CrossRef]
  7. J. B. Pendry, "Negative refraction makes a perfect lens," Phys. Rev. Lett. 85, 3966-3969 (2000).
    [CrossRef] [PubMed]
  8. J. B. Pendry and S. A. Ramakrishna, "Focussing light using negative refraction," J. Phys. Cond. Matter 15, 6345-6364 (2003).
    [CrossRef]
  9. D. Maystre and S. Enoch, "Perfect lenses with left-handed material: Alice’s mirror?," J. Opt. Soc. Am. A 21, 122 (2004).
    [CrossRef]
  10. S. A. Ramakrishna, "Physics of negative refractive index materials," Rep. Prog. Phys. 68, 449-521 (2005).
    [CrossRef]
  11. G. W. Milton and N. A. Nicorovici, "On the cloaking effects associated with anomalous localised resonance," Proc. Roy. Lond. A 462, 3027-3059 (2006).
    [CrossRef]
  12. N. A. P. Nicorovici, G. W. Milton, R. C. McPhedran, and L. C. Botten, "Quasistatic cloaking of two-dimensional polarizable discrete systems by anomalous resonance," Opt. Express 15, 6314-6323 (2007).
    [CrossRef] [PubMed]
  13. W. Cai, U. K. Chettiar, A. V. Kildiev and V. M. Shalaev, "Optical Cloaking with metamaterials," Nature 1, 224-227 (2007).
  14. A. Greenleaf, Y. Kurylev, M. Lassas, and G. Uhlmann, "Improvement of cylindrical cloaking with the SHS lining," Opt. Express 15, 12717-12734 (2007).
    [CrossRef] [PubMed]
  15. S. Guenneau and F. Zolla, "Homogenization of three-dimensional finite photonic crystals," JEWA 14, 529-530 (2000) &Progress In Electromagnetics Research 27, 91-127 (2000).
  16. D. P. Gaillot, C. Croenne and D. Lippens, "An all-dielectric route for terahertz cloaking," Opt. Express 16, 3986-3992 (2008).
    [CrossRef] [PubMed]

2008 (1)

2007 (4)

2006 (5)

J. B. Pendry, D. Shurig, D. R. Smith, "Controlling electromagnetic fields," Science 312, 1780-1782 (2006).
[CrossRef] [PubMed]

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

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]

U. Leonhardt, "Optical conformal mapping," Science 312, 1777-1780 (2006).
[CrossRef] [PubMed]

G. W. Milton and N. A. Nicorovici, "On the cloaking effects associated with anomalous localised resonance," Proc. Roy. Lond. A 462, 3027-3059 (2006).
[CrossRef]

2005 (1)

S. A. Ramakrishna, "Physics of negative refractive index materials," Rep. Prog. Phys. 68, 449-521 (2005).
[CrossRef]

2004 (1)

2003 (1)

J. B. Pendry and S. A. Ramakrishna, "Focussing light using negative refraction," J. Phys. Cond. Matter 15, 6345-6364 (2003).
[CrossRef]

2000 (2)

J. B. Pendry, "Negative refraction makes a perfect lens," Phys. Rev. Lett. 85, 3966-3969 (2000).
[CrossRef] [PubMed]

S. Guenneau and F. Zolla, "Homogenization of three-dimensional finite photonic crystals," JEWA 14, 529-530 (2000) &Progress In Electromagnetics Research 27, 91-127 (2000).

1994 (1)

R. C. McPhedran, N. A. Nicorovici, and G. W. Milton,"Optical and dielectric properties of partially resonant composites," Phys. Rev. B 49, 8479-8482 (1994).
[CrossRef]

J. Opt. Soc. Am. A (1)

J. Phys. Cond. Matter (1)

J. B. Pendry and S. A. Ramakrishna, "Focussing light using negative refraction," J. Phys. Cond. Matter 15, 6345-6364 (2003).
[CrossRef]

Nature (1)

W. Cai, U. K. Chettiar, A. V. Kildiev and V. M. Shalaev, "Optical Cloaking with metamaterials," Nature 1, 224-227 (2007).

New J. Phys. (1)

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

Opt. Express (3)

Opt. Lett. (1)

Phys. Rev. B (1)

R. C. McPhedran, N. A. Nicorovici, and G. W. Milton,"Optical and dielectric properties of partially resonant composites," Phys. Rev. B 49, 8479-8482 (1994).
[CrossRef]

Phys. Rev. Lett. (1)

J. B. Pendry, "Negative refraction makes a perfect lens," Phys. Rev. Lett. 85, 3966-3969 (2000).
[CrossRef] [PubMed]

Proc. Roy. Lond. A (1)

G. W. Milton and N. A. Nicorovici, "On the cloaking effects associated with anomalous localised resonance," Proc. Roy. Lond. A 462, 3027-3059 (2006).
[CrossRef]

Progress In Electromagnetics Research (1)

S. Guenneau and F. Zolla, "Homogenization of three-dimensional finite photonic crystals," JEWA 14, 529-530 (2000) &Progress In Electromagnetics Research 27, 91-127 (2000).

Rep. Prog. Phys. (1)

S. A. Ramakrishna, "Physics of negative refractive index materials," Rep. Prog. Phys. 68, 449-521 (2005).
[CrossRef]

Science (3)

J. B. Pendry, D. Shurig, D. R. Smith, "Controlling electromagnetic fields," Science 312, 1780-1782 (2006).
[CrossRef] [PubMed]

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]

U. Leonhardt, "Optical conformal mapping," Science 312, 1777-1780 (2006).
[CrossRef] [PubMed]

Cited By

OSA participates in CrossRef's Cited-By Linking service. Citing articles from OSA journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (3)

Fig. 1.
Fig. 1.

Left: Geometry of the structured infinitely conducting cloak. The structure consists of 256 angular sectors. Right: Diffraction of the field radiated by a magnetic wire source by an infinitely conducting F-shaped obstacle surrounded by a homogeneous dielectric anisotropic coating with effective permittivity ε h o m ¯ ¯ = diag ( 1.7 , 8.2 , ) and effective transmission conditions n ( ε h o m ¯ ¯ 1 ( + i γ ) H hom ) | r = n ( + i γ ) H hom | r + on its inner and outer boundaries (r=R 1=0.144 and r=R 2=0.4). The frequency ν is equal to 3.5.

Fig. 2.
Fig. 2.

Real part ℜe(H 3) of the longitudinal component H 3 of the magnetic field along the x 2-axis. The origin of the x 2 axis is taken at point (4,4) of Fig. 3 and it ends at point (0,0). A magnetic line source of wavelength λ=c/ν=c/3.5 is located at point (2.3,2.3). The vertical thick bold lines represent the outer boundary of the cloak.

Fig. 3.
Fig. 3.

2D plot of the real part ℜe(H 3) of the longitudinal component H 3 of the magnetic field radiated by a harmonic line current source of wavelength λ=c/ν. Left panel: Diffraction by an infinitely conducting F-shaped obstacle; Right panel: Diffraction by an infinitely conducting F-shaped obstacle surrounded by the structured cloak. When the frequency ν increases, the diffraction worsens and cloaking becomes less effective.

Equations (8)

Equations on this page are rendered with MathJax. Learn more.

r = R 1 + r ( R 2 R 1 ) R 2 , 0 r R 2 ,
θ = θ , 0 < θ 2 π ,
x 3 = x 3 , x 3 I R ,
ε r = μ r = r R 1 r , ε θ = μ θ = r r R 1 , ε 3 = μ 3 = ( R 2 R 2 R 1 ) 2 r R 1 r .
ε r = ( R 2 R 2 R 1 ) 2 ( r R 1 r ) 2 , ε θ = ( R 2 R 2 R 1 ) 2 , μ 3 = 1 .
ε h o m ¯ ¯ = 1 area ( Y * ) ( area ( Y * ) 0 0 0 area ( Y * ) 0 0 0 ) ( ϕ rr ϕ 0 ϕ θr ϕ θθ 0 0 0 0 ) ,
i , j { r , θ } , ϕ ij = < V j y i > Y * = < V i y j > Y * = < V i · V j > Y * ,
𝒦 j : Δ V j = 0 , in Y * = Y B , and V j n = n j , on B ,

Metrics