APPARATUS FOR NUCLEAR MAGNETIC RESONANCE UTILIZING METAMATERIALS OR DIELECTRIC MATERIALS

Abstract

An apparatus for increasing efficiency in the transmission phase and sensitivity in the reception phase, in specific regions of space, of magnetic resonance imaging technique by using at least one metamaterial or dielectric material is provided. Placing the metamaterial or dielectric material in a suitable geometry, in the space delimited by an RF coil and a sample, allows using the surface plasmonic resonances or equivalent dielectric resonances, induced in the metamaterial or dielectric material by the RF coil, to amplify the intensity of the magnetic field in the spatial region of the sample, improving the intensity of the signal transmission and/or the sensitivity of detection. The metamaterial or dielectric material is positioned outside the RF coil to maximize the amplification effect.

Claims

1. An apparatus for nuclear magnetic resonance analysis of nuclear or electronic spin of a sample containing at least one nucleus and/or one electronic spin of interest, comprising means for producing a static magnetic field and the following elements positioned along an axis (z): induction means at a predefined position along the axis (z) and having a maximum transverse dimension ρ.sub.0>0 perpendicularly to said axis (z), said induction means being tuned around a Larmor frequency defined on the basis of said static magnetic field and of the at least one nucleus and/or electronic spin of interest; at least one sample housing; and at least one metamaterial or dielectric material, having a dimension l.sub.m>0 along said axis (z) between a first plane and a second plane perpendicular to the axis (z), the first plane being further from said at least one sample housing and the second plane being closer to said at least one sample housing, along the axis (z); wherein: the at least one sample housing is bounded by a plane, perpendicular to the axis (z), said plane being the closest to the at least one metamaterial or dielectric material; a real value quantity d.sub.m is defined which represents a difference between a position along the axis (z) of the induction means and the position along the axis (z) of the first plane of the at least one metamaterial or dielectric material, and a real value quantity d.sub.s which represents the difference between the position along the axis (z) of the at least one sample housing and the position along the axis (z) of the induction means; and wherein: the at least one metamaterial is configured to develop a magnetic surface plasmonic regime; the at least one metamaterial has a relative magnetic permeability with negative real part; the at least one dielectric material is configured to develop a dielectric resonances regime; the at least one dielectric material has a relative dielectric permittivity ε.sub.d with positive real and imaginary parts; said induction means face said first or said second plane; and a condition d.sub.s+d.sub.m≥0 applies.

2. The apparatus of claim 1, wherein the real value quantity d.sub.m ranges from 0 to the maximum traverse dimension ρ.sub.0 of said induction means.

3. The apparatus of claim 1, wherein the real value quantity d.sub.m is comprised between 0 and 1/10 of the maximum traverse dimension ρ.sub.0 of said induction means.

4. The apparatus of claim 1, wherein the real value quantity d.sub.s ranges from 0 to the maximum traverse dimension ρ.sub.0 of said induction means.

5. The apparatus of claim 4, wherein the real value quantity d.sub.s is comprised between 0 and 1 cm.

6. The apparatus of claim 1, wherein said at least one metamaterial or dielectric material consists of a flat slab with two opposite faces, lying on said first and second planes, wherein a minimum radius of curvature of at least one of the two opposite faces of the flat slab is greater than the maximum transverse dimension ρ.sub.0 of the induction means.

7. The apparatus of claim 6, wherein said at least one metamaterial has a relative magnetic permeability μ.sub.m such that Re(μ.sub.m) is in a range of about −1, said range having a width equal to 2.Math.Im(μ.sub.m).

8. The apparatus of claim 6, wherein: the at least one metamaterial has a thickness l.sub.m between the two opposite faces such that l.sub.m> 1/10 of the maximum transverse dimension ρ.sub.0 of the induction means; and the maximum transverse dimension ρ.sub.0 of the induction means is ρ.sub.0<2/[l.sub.m.sup.−1 log(2/Im())].

9. The apparatus of claim 1, wherein said at least one metamaterial or dielectric material has a spherical or spheroidal shape.

10. The apparatus of claim 9, wherein said at least one metamaterial has a relative magnetic permeability μ.sub.m such that Re(μ.sub.m) is negative and approximate to the first order by ${\mu }_m=-\frac{1+L}{L}$ wherein L is a positive integer which identifies a magnetic plasmonic regime.

11. The apparatus of claim 9, wherein said at least one metamaterial has a relative magnetic permeability μ.sub.m such that Im(μ.sub.m) is less than 0.3.

12. The apparatus of claim 11, wherein said at least one metamaterial has a relative magnetic permeability μ.sub.m such that Im(μ.sub.m) is less than 0.1.

13. The apparatus of claim 9, wherein the relative dielectric permittivity ε.sub.d of said at least one dielectric material with relative magnetic permeability μ.sub.m=1 satisfies the equation:
φ.sub.L.sup.(1)(k.sub.mρ.sub.m)−μ.sub.mφ.sub.L.sup.(1)(k.sub.dρ.sub.m)=0. wherein φ.sub.L.sup.(1)(ξ)=(d[ξj.sub.L(ξ)]/dξ)/j.sub.L(ξ), j.sub.L are spherical Bessel functions, ξ=k.sub.mρ.sub.m, ρ.sub.m is the radius of the sphere, k.sub.m=√{square root over (ϵ.sub.mμ.sub.m)}k.sub.0 and k.sub.d=√{square root over (ϵ.sub.d)}k.sub.0 with k.sub.0=ω/c, where ω=2πv and v is a working frequency of the induction means, and wherein μ.sub.m is the relative magnetic permeability of a metamaterial sphere.

14. The apparatus of claim 10, wherein the real value quantity d.sub.m is chosen as a function of a geometry of the at least one metamaterial or dielectric material and said magnetic plasmonic regime or dielectric resonances regime, respectively.

15. An apparatus for nuclear magnetic resonance analysis of a sample containing at least one nucleus and/or one electronic spin of interest, comprising means for producing a static magnetic field and the following elements: induction means tuned around a Larmor frequency defined on the basis of said static magnetic field and of the at least one nucleus and/or electronic spin of interest; at least one sample housing; and at least one metamaterial or dielectric material; wherein the induction means, the at least one sample housing, and the at least one metamaterial or dielectric material have a development along respective concentric arcs of circumference; said induction means are located between said at least one metamaterial or dielectric material and the at least one sample housing; the at least one metamaterial is configured to develop a magnetic plasmonic regime; the at least one metamaterial has a relative magnetic permeability μ.sub.m with negative real part; the at least one dielectric material is configured to develop a dielectric resonances regime; and the at least one dielectric material has a relative dielectric permittivity ε.sub.d with positive real and imaginary parts.

16. The apparatus of claim 15, wherein said respective concentric arcs of circumference are 360° arcs.

17. The apparatus of claim 15, wherein said at least one metamaterial and said induction means consist respectively of a plurality of metamaterials or dielectric materials and a plurality of induction means, positioned in consecutive and separate portions of their respective concentric arcs of circumference.

18. The apparatus of claim 1, wherein said at least one metamaterial has a relative magnetic permeability μ.sub.m such that Im(μ.sub.m) is smaller than 10.sup.−1.

19. The apparatus of claim 18, wherein Im(μ.sub.m) is smaller than 10.sup.−2 or 10.sup.−3.

20. The apparatus of claim 1, wherein said at least one metamaterial or dielectric material displays at least two poles tuned to two different Larmor frequencies of at least two corresponding nuclei of interest.

21. The apparatus of claim 1, wherein at least one induction coil is inserted between said at least one metamaterial or dielectric material and said at least one sample housing.

22. The apparatus of claim 1, wherein the real value quantities d.sub.s and d.sub.m are both positive.

23. The apparatus of claim 15, wherein said at least one metamaterial has a relative magnetic permeability μ.sub.m such that Im(μ.sub.m) is smaller than 10.sup.−2.

24. The apparatus of claim 23, wherein Im(μ.sub.m) is smaller than 10.sup.−2 or 10.sup.−3.

25. The apparatus of claim 15, wherein said at least one metamaterial or dielectric material displays at least two poles tuned to two different Larmor frequencies of at least two corresponding nuclei of interest.

26. The apparatus of claim 15, wherein at least one induction coil is inserted between said at least one metamaterial or dielectric material and said at least one sample housing.

Description

DETAILED DESCRIPTION OF EXAMPLES OF PREFERRED EMBODIMENTS OF THE INVENTION

List of Figures

[0024] The invention will now be described by way of example, with particular reference to the drawings of the accompanying figures, in which:

[0025] FIG. 1 shows in (a) the geometry A of a configuration according to the invention which comprises: a conventional RF coil (of circular shape, reference “C”) positioned in the center of the reference system (x,y,z) and whose principal axis (z) is perpendicular to the applied static magnetic field {right arrow over (B.sub.0)}; the sample (reference “S”, thickness l.sub.s, permittivity ϵ.sub.s, conductivity σ.sub.s, permeability μ.sub.s, assumed to be of transverse dimension greater than the dimension of the coil along axes x and y, positioned at distance (or “quantity”) d.sub.s from the RF coil); a metamaterial slab (referenced by “MM”, thickness l.sub.m, permittivity ϵ.sub.m, permeability μ.sub.m=[Re (μ.sub.m)+Im (μ.sub.m)]), supposedly of dimensions larger than the RF coil along x- and y-axes, positioned at distance (or “quantity”) d.sub.m from the RF coil. (b) Construction detail of the conventional RF coil with radius ρ.sub.0 and radial width w; the spatial coordinates (ρ, ϕ, z) are used to identify points of interest for calculating (or measuring) the detection (counter-rotating) RF magnetic field per current unit of the RF coil |B.sub.1.sup.(−)|/μ.sub.0 (in A/m), the excitation (co-rotating) RF magnetic field per unit of RF oil current |B.sub.1.sup.(+)|/μ.sub.0 (in A/m), the RF electric field per unit of RF coil current [E], in absolute (V/m) or normalized units |E.sup.(n)|=|E|/[E.sub.0], where [E.sub.0] is the maximum value of the electric field in the sample calculated in the configuration without the MM, the specific absorption rate (SAR) per current unit in the RF coil in (W/kg) and the normalized (SNR.sup.(n)=SNR.sup.(m)/SNR.sup.(V)) signal-to-noise ratio (SNR), where SNR.sup.(m), SNR.sup.(V) are calculated in presence and absence of the metamaterial, respectively.

[0026] FIG. 2 shows the graph of |B.sub.1.sup.(−)|/μ.sub.0 (curves) and SNR.sup.(n) (curves with square symbols) as a function of: (a) Re(μ.sub.m) (assuming Im(μ.sub.m)=0.01); (b) Im(μ.sub.m)(assuming Re(μ.sub.m)=−1). In both cases, |B.sub.1.sup.(−)|/μ.sub.0 and SNR.sup.(n) are calculated at

[00001]$\rho =2\mathrm{cm};\varphi =\frac{\pi }{2};z=1\mathrm{mm};$

having assumed d.sub.m=d.sub.s=0 mm.

[0027] FIG. 3 shows a graph as in FIG. 2 where |B.sub.1.sup.(−)|/μ.sub.0 and SNR.sup.(n) are calculated in the point: ρ=2 cm; ϕ=π/2; z=3 mm; for d.sub.m=0 mm; d.sub.s=2 mm.

[0028] FIG. 4 shows a graph as in FIG. 2 where |B.sub.1.sup.(−)|/μ.sub.0 and SNR.sup.(n) are calculated in the point:

[00002]$\rho =2\mathrm{cm};\varphi =\frac{\pi }{2};z=3\mathrm{mm};$

for d.sub.m=2 mm; d.sub.s=2 mm.

[0029] FIG. 5 shows a two-dimensional map in the plane (ρ,z) of SNR.sup.(n) in the presence of the metamaterial slab assuming

[00003]$\mathrm{Re}\left({\mu }_m\right)=-1,d_m=d_s=0\mathrm{mm},\varphi =\frac{\pi }{2},$$l_m=5.7\mathrm{cm},l_s=20\mathrm{cm}$$\mathrm{with}:\left(a\right)\mathrm{Im}\left({\mu }_m\right)=0.1;\left(b\right)\mathrm{Im}\left({\mu }_m\right)=0.01;\left(c\right)\mathrm{Im}\left({\mu }_m\right)=0.001;$

panel (d) same parameters as (c) but with d.sub.s=4 mm; the white dotted curves correspond to level SNR.sup.(n)=1.

[0030] FIG. 6 shows a two-dimensional map in the plane (ρ, z) of log.sub.10|B.sub.1.sup.(+,n)|, with |B.sub.1.sup.(+,n)|=|B.sub.1.sup.(+,MM)|/|B.sub.1.sup.(+,y)|, where |B.sub.1.sup.(+,MM)| is the excitation RF magnetic field calculated in presence of the metamaterial, |B.sub.1.sup.(+,V)| is the maximum of the absolute value of the same quantity, within the sample, in the configuration without the metamaterial; in the figure, it is assumed

[00004]$d_m=d_s=0\mathrm{mm},\varphi =\frac{\pi }{2},l_m=5.7\mathrm{cm},l_s=20\mathrm{cm}$

and with: (a) μ.sub.m=1, (i.e. vacuum instead of metamaterial); (b) μ.sub.m=−1+i0.1; (c) μ.sub.m=1+i0.01; (d) μ.sub.m=−1+i 0.001. The white dotted curves correspond to level log.sub.10|B.sub.1.sup.(+,n)|=0.

[0031] FIG. 7 shows a two-dimensional map in the plane (ρ,z) of SNR.sup.(n) and log.sub.10|B.sub.1.sup.(+,n)| assuming

[00005]$\varphi =\frac{\pi }{2},{\mu }_m=-1+i0.01,l_m=5.7\mathrm{cm},l_s=20\mathrm{cm}$

and with: (a) d.sub.m=d.sub.s=0 mm; (b) d.sub.m=0 mm, d.sub.s=2 mm; (c) d.sub.m=d.sub.s=2 mm. The white dotted curves correspond to level SNR.sup.(n)=1.

[0032] FIG. 8 shows the level curves of the normalized electric field |E.sup.(n)|=|E|/|E.sub.0| in the plane (ρ,z), where |E.sub.0| is the maximum value of the electric field in the sample calculated in the configuration without metamaterial, l.sub.m=5.7 cm, l.sub.s=20 cm and with: (a) d.sub.m=d.sub.s=0 mm; (b) d.sub.m=0 mm, d.sub.s=2 mm; (c) d.sub.m=d.sub.s=2 mm.

[0033] FIG. 9 shows the dependency along axis z of: (a) |B.sub.1.sup.(+)|/μ.sub.0 and (b) both calculated for

[00006]${\mu }_m=-1+i0,1,\varphi =\frac{\pi }{2},\rho =0\mathrm{cm},d_m=d_s=0\mathrm{mm}$

metamaterial thickness values l.sub.m between 0 cm and 11 cm and l.sub.s=20 cm; in (c) the maximum value of |B.sub.1.sup.(+)(z)/μ.sub.0 and in (d) the maximum value of SNR.sup.(n)(z) are shown, both calculated in the corresponding coordinate z inferred from panels (a) and (b), calculated as a function of l.sub.m (between 0 cm and 25 cm) for

[00007]$\varphi =\frac{\pi }{2},\rho =0\mathrm{cm}.$

[0034] FIG. 10 shows a graph as in FIG. 9 where d.sub.m=0 mm, d.sub.s=2 mm.

[0035] FIG. 11 shows a graph as in FIG. 9 where d.sub.m=2 mm, d.sub.s=2 mm.

[0036] FIG. 12 shows a graph as in FIG. 9 where ρ=2 cm.

[0037] FIG. 13 shows a graph as in FIG. 9 where ρ2 cm, d.sub.m=0 mm, d.sub.s=2 mm.

[0038] FIG. 14 shows a graph as in FIG. 9 where ρ=2 cm, d.sub.m=d.sub.s=2 mm.

[0039] FIG. 15 shows a graph as in FIG. 9 where ρ=3 cm.

[0040] FIG. 16 shows a graph as in FIG. 9 where ρ=3 cm, d.sub.m=0 mm, d.sub.s=2 mm.

[0041] FIG. 17 shows a graph as in FIG. 9 where ρ=3 cm, d.sub.s=d.sub.m=2 mm.

[0042] FIG. 18 in (a) shows a layout similar to the one in FIG. 1(a) where the three constituent elements (MM, C, S) are deformed according to a given radius of curvature; in (b) shows a layout similar to the one in FIG. 1 (a) where two constituent parts are deformed to a given radius of curvature and the sample has a circular (or nearly circular) cross-section.

[0043] FIG. 19 in (a) shows a layout similar to the one in FIG. 1(a) where the metamaterial is deformed to a given radius of curvature, the sample has a circular (or nearly circular) cross-section and there are at least two RF coils which can operate in parallel mode in transmission and/or reception; in (b) shows a layout similar to the one in FIG. 1(a) where the metamaterial is deformed to a given radius of curvature and separated into two independent sections, the sample has a circular (or nearly circular) cross-section and there are at least two radio-frequency coils which can operate in parallel mode in transmission and/or reception.

[0044] FIG. 20 in (a) shows a layout similar to the one in FIG. 1(a) where the metamaterial completely surrounds a sample of circular (or nearly circular) cross-section and there are at least two RF coils which can operate in parallel mode in transmission and/or reception; in (b) shows a layout similar to the one in FIG. 1(a) where the metamaterial completely surrounds a sample of circular (or nearly circular) cross-section and at least one RF volume coil is present (e.g. of the birdcage, multiple transmission line type) which can operate in parallel mode in transmission and/or reception; both configurations in (a) and (b) have similar advantages even if the sample does not have a circular cross-section.

[0045] FIG. 21 in (a) shows a layout similar to the one in FIG. 20 (b) where the metamaterial partially surrounds the RF volume coil and the circular (or nearly circular) cross-section sample; in (b) it shows a layout similar to the one in FIG. 18 (a) where there are at least two layers of metamaterial facing the RF coil and the sample; the same configuration has similar advantages even if the sample does not have a circular cross-section; In (c) it shows a layout similar to the one in FIG. 18 (a) where there are at least two layers of metamaterial facing the sample, with the RF coil comprised between the two layers of metamaterial, the same configuration has similar advantages even if the sample does not have a circular cross-section.

[0046] FIG. 22 shows a layout similar to the one in FIG. 1(a) with the MRI configuration considered in geometry A which comprises: a RF coil of standard surface (reference “C”, with radius ρ.sub.0 and radial width w), positioned on the plane z=0 and placed between the magnetic metamaterial (MM) sphere (with radius ρ.sub.m and permeability μ.sub.m=[Re (μ.sub.m)+Im (μ.sub.m)] and the cylindrical sample (reference “S”, with radius ρs, thickness l.sub.s, relative permittivity ε.sub.s, conductivity σ.sub.s, permeability μ.sub.s, positioned at a distance d, from the RF coil. B.sub.0 is a homogeneous static magnetic field applied along the x-axis and d.sub.m is the minimum distance between the magnetic metamaterial and the plane of the RF coil.

[0047] FIG. 23 shows the magnetic field graph |B.sub.1.sup.(−)|/μ.sub.0 (solid line) and the amplitude of SNR.sup.(n) (solid line with star symbols), evaluated within the sample (for ρ=0 mm, z=6 mm, d.sub.m=0 mm, d.sub.s=2 mm), as a function of the real part of the permeability μ.sub.m, when Im (μ.sub.m)=0.01. The star markers highlight the first five local peaks of SNR.sup.(n). As a reference, the dark dashed lines show the permeability values which allow the existence of some MLSP, defined by the equation (20), having considered a magnetic metamaterial isolated in the vacuum in the static limit approximation, with the magnetic mode index L which varies from 3 to 7 (from left to right).

[0048] FIG. 24 shows the profiles of (a) |B.sub.1.sup.(−)|/μ.sub.0 and (b) SNR.sup.(n), for the geometry described in FIG. 22, evaluated along the z-axis of the sample (ρ=0 mm, z≥2 mm) for the configuration without metamaterial and for five exemplified configurations, with the metamaterial having Im(μ.sub.m)=0.01 and Re(μ.sub.m) corresponding to the values marked with a star in FIG. 23.

[0049] FIG. 2S shows, for the geometry described in FIG. 22: in (a) the SNR map within the sample (z≥2 mm, ϕ=π/2) without the magnetic metamaterial sphere; in (b-f) the map of SNR.sup.(n), with the magnetic metamaterial sphere, evaluated for Im(μ.sub.m)=0.01 and Re(μ.sub.m) corresponding to the values marked with a star in FIG. 23; the dashed black lines are curve lines for SNR.sup.(n)=1.

[0050] FIG. 26 shows the transmission field maps |B.sub.1.sup.(+)|/μ.sub.0 with geometry and parameters as in FIG. 25.

[0051] FIG. 27 shows the electric transmission field maps |E| with geometry and parameters as in FIG. 25.

[0052] FIG. 28 shows an example layout of the MRI configuration considered with geometry A according to an embodiment of the invention. A standard surface RF coil (“C”, with radius ρ.sub.0 and radial width w), positioned on the plane z=0, is placed between a positive dielectric constant sphere (“uHDC”, with radius ρ.sub.m and permittivity ϵ.sub.d=[Re (ε.sub.d)+Im (ϵ.sub.d)] and the cylindrical sample (“S”, with radius ρ.sub.s, thickness I.sub.s and relative permittivity ε.sub.s). B.sub.0 is a homogeneous static magnetic field applied along the x-axis and d.sub.m (d) is the minimum distance between the metamaterial (the sample) and the plane of the RF coil.

[0053] FIG. 29 shows, for the geometric configuration of FIG. 28 with Re (ϵ.sub.d)=1200 and radii ρ.sub.m=3.38; 6.21; 8.82 cm which support the resonances L=1; 3; 5, respectively, at the Larmor frequency of 127.74 MHz (B.sub.0=3 T): (a) the reception field |B.sub.1.sup.(−)|/μ.sub.0 and (b) the SNR.sup.(n) in point a=2 mm, ρ=0 mm as a function of the loss tangent tan δ. The lines show the evaluated values without the sphere and with several uHOC spheres, the black dashed vertical lines show the case of tan δ=0.04. The gray area at the top of the panel (b) corresponds to a tan δ range which produces a useful SNR.sup.(n)>1.

[0054] FIG. 30 shows the map in the plane (ρ,z) of the specific absorption rate SAR for the unit current at 127.74 MHz (B.sub.0=3 T) with the same parameters as FIG. 29 without (a) and with (b) the uHDC sphere and for tan δ=0.04 and ρ=8.82 cm.

[0055] FIG. 31 shows the map in the plane (ρ,z) of |B.sub.1,eff.sup.(+)|(ϕ=π/2) with the same parameters as in FIG. 29 without (a) and with (b) the uHDC sphere for tan δ=0.04 and ρ.sub.m=8.82 cm.

[0056] FIG. 32 shows the same geometric configuration as FIG. 28 per Re (ϵ.sub.d)=3300, radii ρ.sub.m=4.08, 7.49, 10.64 cm supporting resonances L=1, 3, 5, respectively, at Larmor frequency of 63.87 MHz (80=1.5 T). (a) |B.sub.1.sup.(−)|/μ.sub.0 and (b) SNR.sup.(n) at point z=2 mm, ρ=0 mm as a function of the tangent of tan losses d. The horizontal lines without symbols are evaluated without the uHDC sphere and vertical dashed black lines highlight the case of tan δ=0.04. The gray area in (b) corresponds to a tan δ range which produces a useful SNR.sup.(n)>1.

[0057] FIG. 3 shows the comparison between the geometries in FIGS. 22 and 28, respectively, with the metamaterial magnetic sphere (μ.sub.m=−1.20+i0.01) and the uHDC sphere (ϵ.sub.d=1324+i 1.65), with reference to resonance L=5 in both cases. The profiles of |B.sub.1.sup.(−)|/μ.sub.0 (a) and SNR.sup.(n) (b) are as a function of ρ within the cylindrical sample (z≥2 mm, ϕ=π/2) for different z-values in the presence of the magnetic sphere of metamaterial (continuous lines) or dielectric (dashed lines). The maps of SNR.sup.(n) for the magnetic sphere of metamaterial (c) or dielectric (d) refer to the plane (ρ, z) within the sample (z>=2 mm, ϕ=π/2). Black dashed lines are level lines for SNR.sup.(n)=1, white dashed lines for SNR.sup.(n)=3, and SNR.sup.(n)=1.5.

[0058] FIG. 34 shows an example layout of an MRI configuration with geometry B according to an embodiment of the invention, with the magnetic sphere MM, or uHDC positioned between a standard surface RF coil and the cylindrical sample.

[0059] FIG. 35 shows the maps, in the plane (ρ,z), of |B.sub.1.sup.(−)|/μ.sub.0 inside the cylindrical sample (d.sub.m=0 mm; z≥2 mm) in the presence of: magnetic sphere MM (μ.sub.m=−1.2+i0.01) with (a) geometry A (d.sub.s=2 mm) or (b) geometry B (d.sub.s=ρ.sub.m+2 mm); the uHDC sphere (ε.sub.d=1324+i 1.65) with (c) geometry A (d.sub.s=2 mm) or (d) geometry B (d.sub.s=ρ.sub.m+2 mm). Profile comparison of |B.sub.1.sup.(−)|/μ.sub.0 on axis ρ=0 mm for geometries A and in the presence of the MM sphere (a) or uHDC sphere (t) defined above.

[0060] FIG. 36 shows the graphs as in FIG. 35 for the field |E|.

[0061] FIG. 37 shows the graphs as in FIG. 35 for the S.

[0062] It is worth noting that hereinafter elements of different embodiments may be combined together to provide further embodiments without restrictions respecting the technical concept of the invention, as a person skilled in the art will effortlessly understand from the description.

[0063] The present description also makes reference to the prior art for its implementation, with regard to the detail features which not described, such as, for example, elements of minor importance usually used in the prior art in solutions of the same type.

[0064] When an element is introduced it is always understood that there may be “at least one” or “one or more”.

[0065] When a list of elements or features is given in this description it is understood that the invention according to the invention “comprises” or alternatively “consists of” such elements.

[0066] In the description of the embodiments, reference will generally be made to a sample to be subjected to magnetic resonance imaging (NMR/MRI/EPR/EPRI) and containing at least one electronic or nuclear spin of interest.

[0067] Furthermore, reference will be made to an “induction coil” or “RF coil” or even just “coil” meaning a coil that generates a non-static electric and/or magnetic field at radio frequencies or even microwaves or other useful frequencies. The term “RF coil” is also used in literature for frequencies other than radio-frequencies to distinguish this coil from other coils present in magnetic resonance equipment, such as coils for static magnetic fields, coils for magnetic field gradients necessary for spatial localization of the resonance signal.

[0068] Furthermore, the coil can have any cross-section shape (plane x,y in the figures) and thus in general we will speak of maximum transverse dimension instead of diameter in the circular case.

[0069] In general, the coil is tuned (e.g. In a bandwidth) about the Larmor frequency defined based on the static magnetic field and at least the electronic or nuclear spin of interest.

[0070] In this context, according to the invention, that illustrated as a technical effect for the case of magnetic plasmons also applies to electric plasmons. Indeed, the metamaterial can develop a surface plasmonic regime with electrical resonances (see [1, 3] and references cited therein), by appropriately selecting a negative dielectric permittivity value (ε.sub.m).

[0071] In the case of a slab of infinite transverse dimension (dimension x, y of FIG. 1 (a)) and finite thickness 4, this value is equal to Re(ε.sub.m)=−1. An additional geometric configuration of metamaterial capable of supporting the electrical plasmonic regime is that of a sphere, the negative dielectric permittivity value of which must satisfy, with a given degree of approximation, the following condition Re(ε.sub.m)=−[(1+L)/L)].

[0072] The excitation means of the magnetic and/or electrical plasmonic resonance must be appropriately chosen from the possible configurations which can be divided between methods with an internal or external metamaterial source. For example, a method may be used with a small circular RF coil (or other shapes) which has its axis oriented at a given angle variable between 0° and 90° with respect to the surface of the stab of MM (i.e. relative to an axis lying in the x-y plane in FIG. 1(a)). Alternatively, an RF coil may be used which has at least one linear current element in the plane of the coil itself (eight-shaped coil, or double-O coil). In the prior art, resonant transmission lines (microstrip transmission lines) are also used, which have at least one linear conductive element terminating on a capacitor, the axis of which must be appropriately oriented relative to the z-axis of FIG. 1(a).

[0073] Therefore, the coil or the excitation means (or more generally “induction means”) can also or only perform the function of excitation of electrical surface plasmons. Furthermore, as reported in the prior art of antenna theory, such excitation can occur by means of the use of a linear dipole induction coil, the main axis of which must be appropriately aligned with the electrical modes that the metamaterial can support.

[0074] In this respect, the excitation procedures of the magnetic metamaterial (finished slab, cylinder, sphere, spheroid, cube, parallelepiped, etc.) also apply to the resonance excitation method of the dielectric material, and the choice of method depends on the shape of the dielectric itself and the chosen resonance mode. The implementation details in individual cases can be obtained analytically, as in the examples below, or numerically, following methodologies well known in the literature [2] and verified by the inventors.

Embodiments

[0075] With the present invention, a step forward is made in the use of metamaterials for magnetic resonance imaging, by suggesting the use of excited surface plasmons on at least one surface of a magnetic type metamaterial (e.g. for a slab that has Re (μ.sub.m)=−1) and of electrical type (e.g. for a slab having Re (ε.sub.m)=−1), as far as it is possible to approximate these conditions in reality.

[0076] For the first time, to the knowledge of the inventors, it is shown that the resonant nature of magnetic surface plasmons can be appropriately exploited to improve the efficiency of magnetic resonance imaging. Here, by way of example, a metamaterial slab will be considered characterized by Re (μ.sub.m)=−1 and incorporated in a magnetic resonance configuration as shown in FIG. 1. In this configuration, we will show that the metamaterial supports magnetic surface plasmons and their excitations can increase the magnetic field useful to excite the sample (in general, containing at least one active nuclear spin and/or an electronic spin of interest) and/or increase the magnetic resonance signal-to-noise ratio (SNR) relative to the current settings.

[0077] In an attempt to exploit the high local fields associated with surface plasmons by keeping the RF coil (or in general a coil or induction means which can also generate microwaves or other frequencies) on the surface as close as possible to the sample, there is suggested the configuration shown in FIG. 1a), wherein the coil C is located between the metamaterial slab MM and the sample S. The considered configuration geometry has the added advantage of not introducing limitations to the relative position between the coil C and the sample S by placing the metamaterial slab MM in a region usually free in many magnetic resonance configurations. In the situation in which the distances (or real value quantities in general, because they can be negative) d.sub.m and d.sub.s are small (compared to the dimensions of the RF col, d.sub.m≅d.sub.s≅0 mm) and the thickness of the metamaterial is large (compared to the dimensions of the RF coil), we expect, based on known theories [1,2],that the metamaterial slab with Re(μ.sub.m)=−1 supports the magnetic surface plasmons, which are located in a reduced thickness on the surface of the metamaterial MM facing the coil C and on the surface of the opposite metamaterial MM, away from the coil C. The magnetic surface plasmons provide, following the excitation of such resonance, a considerable increase in the electromagnetic field within the sample. By performing appropriate full-wave numerical simulations of the electromagnetic field, the configuration shown in FIG. 1 a) was analyzed and the spatial distribution of the non-static magnetic field, as well as the spatial distribution of the SNR, was assessed. In the numerical examples, the chosen frequency v.sub.0=63.866 MHz (where v.sub.0 is the Larmor frequency of the hydrogen nucleus spin corresponding to astatic magnetic field |.sub.0|=1.5 T), l.sub.m=5.7 cm (thickness of the metamaterial slab MM or in general dimension along said axis z between a first plane and a second plane perpendicular to the axis Z which define the ends of the material along the same axis, the first plane being farther from said at least one housing of the sample S and the second plane being closer to said at least one housing of the sample S, along the axis z), l.sub.s=20 cm (thickness of the sample slab), the relative permittivity of the sample ε.sub.s=90 and a conductivity equal to σ.sub.x=0.69 S/m (the latter two values corresponding to the average of the known values for human tissues at the considered frequency). The coil C is modeled with a negligible thickness along the z-axis and a surface current density which has only one azimuthal component, i.e. J.sub.φ=Kδ(z), where K.sub.ϕ(ρ)=b.sub.0.Math.ρ.Math.exp[−(ρ−ρ.sub.0).sup.2/w.sup.2] being δ(⋅) the Dirac delta function, ρ.sub.0=2 cm, w=2 mm, boa constant whose value allows a unit current to be defined on the coil C. Further tests were done with l.sub.m between 1 cm and 5.7 cm still achieving an increase in signal-to-noise ratio. For values smaller than 1 cm, we noted that the improvement introduced by one side of the slab was canceled by the contribution of the other. In general, it can be said that l.sub.m> 1/10 of the transverse dimension (relative to the z or the cod axis) of the maximum induction coil C, however this is a preferred value and the minimum quantity depends on the whole system configuration: It can be calculated each time with analytical and/or numerical methods or by experimentally verifying the existence of plasmonic regimes and the effect of the electromagnetic field produced in the sample in each position of interest (it could affect only a very narrow area of the sample and consequently only some configurations of the magnetic plasmonic regime or dielectric regime resonances).

[0078] Although the distance d.sub.s should ideally be close to or equal to 0 mm, for safety reasons it is still set to a few mm, in any case preferably less than 1 cm. More in general, the maximum distanced between said at least one induction coil C and said at least one sample S housing (relative to a plane tangent to its end along the axis z closer to the metamaterial or dielectric material) is comprised in the range from 0 to the maximum transverse dimension of the induction coil.

[0079] In general, the distance d.sub.m is defined between at least one metamaterial MM and at least one induction coil C, or also as the difference between the position along the axis z of the induction means C and the position along the z axis of the first plane of the metamaterial MM or the dielectric material uHDC. d.sub.s can instead be defined as a real values quantity which represents the difference between the position along the z axis of the sample S housing and the position along the z axis of the induction means C. Both can be comprised in the range from 0 to the maximum transverse dimension of the induction coil, preferably between 0 and 1/10 of the maximum transverse dimension of the induction coil. d.sub.s can be comprised between 0 and 1 cm. It is worth noting that in all embodiments the various components S, MM, uHDC, C are positioned one after the other (in the order given each time or claimed) along the z-axis, but this does not mean that they must have symmetry relative to this axis or cannot be offset in directions perpendicular to such axis.

[0080] It is possible to obtain the aforesaid static magnetic field through a permanent magnet, or an electromagnet, or a superconducting magnet, or in general by means of a static magnetic field.

[0081] In a generic configuration for magnetic resonance, the signal from the sample detected by an RF coil is given by S∞|B.sub.1.sup.(−)(ρ, ϕ, z)|. Considering the geometry shown in FIG. 1 a), where .sub.0 is along the x-axis, it holds B.sub.1.sup.(−=(Bsin(ϕ)+iB)/2, where the RF magnetic field in cylindrical coordinates is given by .sub.1 (ρ, ϕ, z)=Re[(B{circumflex over (ρ)}+B{circumflex over (z)})], with ω=2πv and v the Larmor frequency of the spin of interest. On the other hand, the noise received by the RF coil is proportional to the square root of the power P dissipated in the system, so that the SNR of the receiving RF coil is ∝|B.sub.1.sup.(−)/√{square root over (P)}. After the RF coil losses, the power dissipation is expressed as P=P.sub.s+P.sub.m, where P.sub.s and P.sub.m are the power dissipated in the sample and the metamaterial, respectively. To highlight the advantages related to the presence of the metamaterial slab, hereinafter we will consider the normalized signal-to-noise ratio SNR.sup.(n) defined above.

[0082] The advantages related to the presence of the metamaterial slab are apparent in FIGS. 2, 3, 4 which show, as provided by the known theory, a resonant behavior of the RF magnetic field (continuous curve) and, more importantly, the SNR.sup.(n) (curve with square symbols) becomes greater than one in the plasmonic resonance condition (e.g. in FIG. 2 b) SNR.sup.(n)≅5 for μ.sub.m=−1+i10.sup.3).

[0083] To physically understand the role of surface plasmons and the results shown in FIG. 24, we have analytically solved Maxwell's equations in the configuration in FIG. 1 a) where, for simplicity of calculation, we considered that the sample is a semi-infinite slab in the direction of the z-axis (i.e. l.sub.s.fwdarw.∞). Taking advantage of the rotational symmetry of the system about the z-axis, the complex amplitudes of the electric and magnetic fields can be written as E.sub.ϕ=iωA.sub.ϕ, B.sub.1=∇×A.sub.ϕ, where A.sub.ϕ is the azimuthal component of the electromagnetic potential vector. We then look for solutions to Maxwell's equations using the Hankel transform to express A.sub.ϕ(ρ,z)=∫.sub.0.sup.+∞dk.sub.ρk.sub.ρJ.sub.1(k.sub.ρρ)Ã.sub.ϕ(k.sub.ρ, z). Within the static limit, two relevant regimes can be highlighted:

(1) the Pendry regime, in which the metamaterial slab with Re (μ.sub.m)=−1 can behave like a Pendry lens, if the spatial spectrum of the {circumflex over (K)}.sub.ϕ(k.sub.ρ) {Hankel transform of density current K.sub.ϕ(φ)} is different from aero in the region of

k.sub.ρ<<k.sub.l,  (1)

(ii) the plasmonic regime, in which a metamaterial slab with Re(μ.sub.m)=−1 supports surface plasmonic excitations, in the situation where {tilde over (K)}ϕ[k.sub.ρ] is not null in the region of

k.sub.ρ>>k.sub.l,  (2)

being the parameter k.sub.l=l.sub.m.sup.−1 log [2/Im(μ.sub.m)] defined by the geometry of the metamaterial slab and its losses, identified by the imaginary part of μ.sub.m. In general, we may also have a less important plasmonic regime for k.sub.p>k.sub.1.

[0084] In regime (i), the overall amplitude of the potential vector in the region occupied by the sample is

[00008]$\begin{array}{cc}{A}_{\varphi }\left(\rho ,z>0\right)\simeq -\frac{{\mu }_0{\mu }_m}{2}{\int }_0^{+\infty }{\mathrm{dk}}_{\rho }{J}_1\left(k_p\rho \right){\tilde{K}}_{\varphi }{e}^{-k_{\rho }z}& \left(3\right)\end{array}$

while, in the plasmonic regime (ii),we obtain

[00009]$\begin{array}{cc}{A}_{\varphi }\left(\rho ,z>0\right)\simeq \frac{{\mu }_0{\mu }_m}{1+{\mu }_m}{\int }_0^{+\infty }{\mathrm{dk}}_{\rho }{J}_1\left(k_{\rho }\rho \right){\tilde{K}}_{\varphi }{e}^{-k_{\rho }z}.& \left(4\right)\end{array}$

[0085] The Pendry mechanism applies to plane waves whose transverse wave number satisfies condition (i) k.sub.p<k.sub.l and this corresponds to a minimum resolution relative to the image Δ=2πl.sub.ω/log [2/Im(μ.sub.m)] [2]. Considering the definition of K.sub.ϕ(ρ), it can be understood that this regime is achieved when the coil size (ρ) is very large. On the other hand, if the RF coil is small enough, a significant portion of the spatial spectrum of K.sub.ϕ(ρ) can be found in the region k.sub.p>k.sub.l, where surface plasmons can be excited. It is worth noting that equation (3) with μ.sub.m=−1 coincides with the expression of the field potential if the metamaterial is absent. Consequently, in the regime (i) and configuration of FIG. 1 considered here, the metamaterial does not influence the spatial distribution of the electromagnetic field within the sample. From Eq. (4), it is apparent that a very large increase in the amplitude of the field A.sub.ϕ can be obtained in the condition Re(μ.sub.m)=−1 and Im(μ.sub.m)<1. This condition corresponds to the existence of electromagnetic modes located on the surface of the metamaterial at z=0 (with a skin depth which depends on the losses of the metamaterial). If we suppose to use the RF coil shown in FIG. 1 for the reception of the magnetic resonance signal, the intensity of the received signal, as indicated above, depends on the RF field B.sub.1.sup.(−) and condition (ii) may lead to increase it resulting in increased performance of the magnetic resonance system.

[0086] From the theoretical analysis (in the approximation of static regime), the amplitude of the magnetic resonance signal is proportional to the function

[00010]$\begin{array}{cc}f\left(\mathrm{Re}\left({\mu }_m\right).Math.\rho ,\varphi ,z\right)=.Math.\frac{\mathrm{Re}\left({\mu }_m\right)+i\mathrm{Im}\left({\mu }_m\right)}{1+\mathrm{Re}\left({\mu }_m\right)+i\mathrm{Im}\left({\mu }_m\right)}.Math..Math.B_1^{(-)}\left(\rho ,\varphi ,z\right).Math.& \left(5\right)\end{array}$

[0087] which has a maximum, once the spatial position has been fixed, for Re(μ.sub.m)=−1. The full width at half maximum (FWHM) of the function depends on the imaginary part of the relative magnetic permeability and is about 2Im(μ.sub.m). From this, it follows that in an optimized configuration the relative magnetic permeability μ.sub.m of the metamaterial is such that Re(μ.sub.m) is in a range around the value −1, said range being equal to 2-Im(μ.sub.m)

[0088] The physical mechanism considered here is very different from that suggested by Pendry. The Pendry mechanism is due to the fact that the evanescent waves show an exponential, non-intuitive growth within the metamaterial so that the wave modes emitted by the source, which satisfy condition (i), can be transmitted without diffraction for an adequate lens thickness. Instead, the surface plasmons located near the surface of the metamaterial exist in the opposite regime (li) in which the wave modes satisfy the condition given by Eq. (2). In this regime, the metamaterial with Re (μ)=−1 does not behave like a lens and can produce a hyperfocusing of the electromagnetic field near the surface of the metamaterial.

[0089] The spatial visualization of the mechanism is given by FIGS. 5,6, 7 which show the spatial distribution of SNR.sup.(n) and of B.sub.1.sup.(+,n). For Im(μ.sub.m)=10.sup.−1, we obtain a significant spatial modulation of SNR.sup.(n) and the losses, within the metamaterial, are responsible for an overall reduced performance of the receiving system (SNR.sup.(n)<1 throughout the explored region). Considering Im(μ.sub.m)=10.sup.−3 (Im(μ.sub.m)=10.sup.−2), near the surface of the slab, the plasmonic excitation results in a more intense RF electromagnetic field and a strong improvement in SNR.sup.(n), i.e. SNR.sup.(n)≅7.5 (SNR.sup.(n)≅2.5). In FIG. 5-7, dashed lines indicate isolines with SNR.sup.(n)=1. From FIG. 5 it is apparent that, for the considered geometry, with Im(μ.sub.m)<10, the value of SNR which can be achieved in a magnetic resonance experiment, in the presence of metamaterial, may be increased, by a high factor, in the region 0≤ρ<7 cm and 0≤z<2.5 cm.

[0090] To evaluate the impact of surface plasmon excitations on the signal transmitted by an RF coil, we will consider the spatial distribution of the excitation field (transmission) B.sub.1.sup.(+)=(B.sub.1,ρsin ϕ−iB.sub.1,)/2 normalized to the current in the RF coil (C)(FIGS. 6,7). As predicted by theoretical analysis, as Im(μ.sub.m) decreases, the value of k.sub.l Increases (see Eq. [2)], the excited wave modes are more closely confined near the metamaterial-vacuum interface, and the amplification factor increases (Eq. (4)). Consequently, a possible application of the setup according to the invention is related to the transmission phase of the magnetic resonance signal. The high field increase |B.sub.1.sup.(+,n)| can increase the transmission performance of the system by allowing much shorter RF pulses and/or the use of less powerful RF amplifiers (with cost savings and system management), the flip angle of the macroscopic magnetization of the sample in presence of the static magnetic field being equal. Such an effect may be possibly beneficial, also when multiple RF transmission coils are available, to implement parallel transmission magnetic resonance imaging techniques.

[0091] FIG. 8, for the sake of completeness, shows the spatial trend of the normalized electric field which is observed at three geometric configurations in which the mutual distance between coil and/or metamaterial and/or sample varies by a few millimeters.

[0092] Diagrams for the use of metamaterials are given below, but the conclusions are also valid when using a dielectric material appropriately shaped and with a permittivity value (uHOC) selected ad hoc so that it reproduces an equivalent electromagnetic field (with good approximation) relative to that produced by a magnetic and/or electrical metamaterial (see the example of the sphere of magnetic MM below). To better quantify the advantages which can be obtained with the present invention. FIGS. 9-17 show quantities |B.sub.1.sup.(+)| and SNR.sup.(n) as a function of the depth z in the sample, for different geometric configurations and their dependence on the thickness of the metamaterial slab.

[0093] Now with reference to FIGS. 18-21, we will illustrate some examples of possible geometries of the inventive layout of FIG. 1(a).

[0094] In the first of these embodiments, the three basic constituent elements are deformed according to a given radius of curvature, in particular in the variant of FIG. 18(a) the RF induction coil C, sample S and metamaterial MM all have substantially circular ring sections, although the lengths of the sections may vary, e.g. the length of coil C is less than the length of the other two elements. In the variant of FIG. 18 (b, the sample section is circular. In the case of the two opposite faces of the slab, the minimum radius of curvature of at least one of the two opposite faces of the slab is greater than the maximum transverse dimension of the induction coil.

[0095] Furthermore, in FIG. 18 as in each of the other figures and embodiments and variants, the sample must be understood as a volume of interest in the matter placed in a housing (not shown) of the magnetic resonance apparatus according to the invention. So, for example, the sample in FIG. 18 (b) may be an ROI within a body with a cubic outer shape, without loss of generality.

[0096] The embodiment in FIG. 19 comprises the elements as in FIG. 18 (b), in which the sample has a circular cross-section but the RF induction coil C is double, with shorter sections and the metamaterial MM is an arc of circumference in a single piece or is double, in this case consisting of identical or different MM elements based on the local properties of the sample adjacent to each one. The same applies to a subdivision of the RF coil elements into several parts, beyond the two shown, which can operate in parallel mode in transmission and/or reception.

[0097] The embodiment in FIG. 20 comprises in (a) a whole or almost whole-ring metamaterial and a circular section sample, while a plurality of curved RF coils C#1, C#2, C#3, etc. is present and in (b) instead there is a ring-shaped RF coil, which is also a whole or an almost whole ring. For these cases, the specifications on the plasmonic regime—and therefore on the choice of the metamaterial—can be calculated by numerical simulation.

[0098] Again, in the embodiment in FIG. 21(a) only the metamaterial is not circular in section, but extends for an arc of circumference, while in the embodiment in FIG. 21(b) there are two metamaterials (MM #1 and MM #2) along two concentric arcs of circumference and now the sample also extends along an arcofcircumference.Ingeneral,theremaybemorethantwometamaterialseven nonconcentric, and also flat. The combination of surface plasmonic resonances determined by the geometry and magnetic permeability of the two metamaterials (MM #1 and MM #2), coupled through at least one RF coil, will produce an RF magnetic field distribution inside the sample which can be modulated appropriately with beneficial effects for the magnetic resonance experiment.

[0099] Following the principle described in FIG. 1(a), the case of a slab of thickness l.sub.m and finite transverse dimension can also be considered. In this case, the geometric figure of the slab is transformed into that of a cylinder with one of the bases facing towards the circular coil. Such geometry may be modified, without losing the effectiveness of the metamaterial, by rotating the cylinder by an angle between 0° and 90°, i.e. by orienting its axis of symmetry in a direction which goes from parallel to perpendicular to the z-axis. Here, too, the specifications on the plasmonic regime can be calculated numerically.

Examples of Study of Operation

I. Details on the Numerical Simulation

[0100] The full-wave numerical results shown in FIGS. 2, 3, 4 are obtained by means of the commercial software package COMSOL Multiphysics. Taking advantage of the invariance by rotation around the axis of the RF induction coil (i.e. the axis z), we have performed 2D simulations the validity of which has been confirmed, in some specific geometries, by 3d full-wave simulations performed with Ansys Electromagnetic Desktop software. In the simulations, we will consider a finite dimension spatial domain in which, along the axis z, we considered two (not shown) vacuum regions of thickness l.sub.v=8.5 cm, the first at the metamaterial surface far from the RF coil and the second beyond the sample. Furthermore, perfect electrical conductor (PEC) boundary conditions haw been imposed on the spatial domain frontier. We used appropriate non-homogeneous spatial domain discretization with a maximum grid dimension of 1.5 mm (about 8×10.sup.5 degrees of freedom).

II. RF Coil Signal Calculation in an NMR/MRI Apparatus

[0101] In the configuration considered in FIG. 1, the RF coil can be used to transmit an RF pulse or receive the induction signal caused by the spin of the sample. Bearing in mind that the static magnetic field .sub.0, in FIG. 1 a), is along the axis x, the RF magnetic .sub.1=Re [B.sub.1e.sup.−1ωt] (ω is the angular frequency of the radiation) can be broken down into two contributions

[00011]$\begin{array}{cc}{B}_1^{(+)}=\frac{\left(\text{?}\right)}{2},B_1^{(-)}=\frac{\left(\text{?}\right)}{2},& \left(6\right)\end{array}$$\text{?}\text{indicates text missing or illegible when filed}$

where B.sub.l=B.sub.1x{tilde over (e)}.sub.x+B.sub.1yė.sub.y+B.sub.1z{tilde over (e)}.sub.z, is the alternating magnetic field per unit current flowing in the RF coil. Here we will use the symbols B.sub.1.sup.() to distinguish the two circular polarizations which rotate in opposite directions: B.sub.1.sup.(+) is the polarized field rotating in the same direction as the spin precession (transmission), B.sub.1.sup.(−) is the counter-rotating component (reception). Considering the cylindrical coordinates (ρ, ϕ, z), as defined in FIG. 1, the previous equations become

[00012]$\begin{array}{cc}{B}_1^{(+)}=\frac{(\text{?}\sin \varphi +\text{?}\cos \varphi -i\text{?}}{2},B_1^{(-)}=\frac{\text{?}\sin \varphi +\text{?}\cos \varphi +i\text{?}}{2}.& \left(7\right)\end{array}$$\text{?}\text{indicates text missing or illegible when filed}$

[0102] In our simulations, the surface current density has only one azimuthal component and the system has rotational symmetry, so we can write

[00013]$\begin{array}{cc}{B}_1^{(+)}\left(\rho ,\varphi ,z\right)=\frac{{\text{?}\sin \varphi -i\text{?})}^{*}}{2},B_1^{(-)}\left(\rho ,\varphi ,z\right)=\frac{\text{?}\sin \varphi +i\text{?}}{2}.& \left(8\right)\end{array}$$\text{?}\text{indicates text missing or illegible when filed}$

[0103] The co-rotating component B.sub.1.sup.(+) is the relevant component for the transmission of RF signals which causes the sample spin transitions. On the other hand, considering the principle of reciprocity, the received RF signals am proportional to B.sub.1.sup.(−)* (i.e. the complex conjugate of the counter-rotating RF magnetic field component per current unit), so the signal of the receiving RF coil is simply given by

S∝|B.sub.1.sup.(−)(ρ,ϕ,z)|.  (9)

III. Analytical Expression of the Electromagnetic Vector Potential

[0104] Here, from Maxwell's equations, we can obtain the analytical expression of the electromagnetic vector potential generated by the current flowing in the RF coil in the configuration described in FIG. 1. We will take into consideration the case in which the sample is a semi-infinite slab (l.sub.s.fwdarw.∞), d.sub.m=d.sub.s=0 cm and assume a negligible thickness for the RF coil. Furthermore, we will consider the dimensions of the metamaterial and the sample, along the directions orthogonal to the axis of symmetry z, much larger than the diameter of the RF coil.

[0105] Maxwell's equations admit a monochromatic solution of the shape =Re[A.sub.ϕ(ρ, z)e.sup.(−iωt)], where A.sub.ϕ=A.sub.ϕ{circumflex over (ϕ)} is the azimuthal component of the electromagnetic vector potential. Considering the Lorenz gauge (i.e. the electric and magnetic field are given by E.sub.ϕ=iωA.sub.ϕ, B.sub.1=∇×A, respectively), the spatial dynamics of the potential vector A.sub.ϕ(ρ,z) is ruled by the equation

∇.sup.2A.sub.ϕ+μ.sup.−1∇μ×(∇×A.sub.ϕ)+εμk.sub.0.sup.2A.sub.ϕ=−μ.sub.0μJ.sub.ϕ  (10)

[0106] where k.sub.0=ωc, J.sub.ϕ(ρ, z) is the current density of the RF coil, ε, μ represent the complex dielectric permittivity and the complex magnetic relative permeability, respectively, of the materials considered (c is the vacuum light speed, μ.sub.0 is the vacuum magnetic permeability). Considering the configuration shown in FIG. 1, in which the metamaterial and the sample are assumed to be homogeneous, permittivity and permeability depend only on the coordinate z. The current density distribution of the RF coil is given by J.sub.ϕ=K.sub.ϕ(ρ)δ(z){circumflex over (ϕ)}. By using the Hankel transform we can write A.sub.ϕ(ρ, z)=∫.sub.0.sup.+∞dk.sub.ρk.sub.ρJ.sub.1(k.sub.ρρ) Ā.sub.ϕ(k.sub.ρ, z), and the potential vector equation becomes:

[00014]$\begin{array}{cc}\frac{d}{\mathrm{dz}}\left({\mu }^{-1}\frac{d{\tilde{A}}_{\varphi }}{\mathrm{dz}}\right)+ϵ\text{?}{\tilde{A}}_{\varphi }=-{\mu }_0\delta \left(z\right){\tilde{A}}_{\varphi }.& \left(11\right)\end{array}$$\text{?}\text{indicates text missing or illegible when filed}$

wherein k.sub.z.sup.2 k.sub.0.sup.2 k.sub.ρ.sup.2 and K.sub.ϕ is the Hankel transform of the RF coil surface current.

[0107] Solving the previous equation, we obtain, for the regions occupied by the vacuum (v), the metamaterial (m) and the sample (s):

[00015]$\begin{array}{cc}{\tilde{A}}_{\varphi }=\{\begin{array}{cc}\text{?}& \mathrm{if}z<-l_m,\\ F_m\text{?}+C_m\text{?}& \mathrm{if}{l}_m\le z\le 0,\\ \text{?}& \mathrm{if}z>0.\end{array}& \left(12\right)\end{array}$$\text{?}\text{indicates text missing or illegible when filed}$

where k.sub.0.sup.()=√{square root over (k.sub.0.sup.2−k.sub.ρ.sup.2)}, k.sub.z.sup.(m)=√{square root over (k.sub.0.sup.2e.sub.mμ.sub.m−k.sub.ρ.sup.2)}, k.sub.z.sup.(s)=√{square root over (k.sub.0.sup.2ε.sub.s−k.sub.ρ.sup.2)}, k.sub.±.sup.(v)=k.sub.±.sup.(m)±μ.sub.mk.sub.z.sup.(v) and k.sub.±.sup.(s)=k.sub.2.sup.(m)±μ.sub.mk.sub.s.sup.(s), C.sub.v, C.sub.m, F.sub.m and F.sub.s are given by

[00016]$\begin{array}{cc}\text{?}=\frac{i2{\mu }_0{\mu }_m\text{?}}{k_{+}^{\left(v\right)}{k}_{+}^{\left(s\right)}\text{?}-k_{-}^{\left(v\right)}{k}_{-}^{\left(s\right)}\text{?}},C_m=\frac{i{\mu }_0{\mu }_m{k}_{+}^{\left(v\right)}{\tilde{K}}_{\varphi }\text{?}}{k_{+}^{\left(v\right)}{k}_{+}^{\left(s\right)}\text{?}-k_{-}^{\left(v\right)}{k}_{-}^{\left(s\right)}\text{?}},F_m=\frac{i{\mu }_0{\mu }_m{k}_{-}^{\left(v\right)}{\tilde{K}}_{\varphi }\text{?}}{k_{+}^{\left(v\right)}{k}_{+}^{\left(s\right)}\text{?}-k_{-}^{\left(v\right)}{k}_{-}^{\left(s\right)}\text{?}},F_s=\frac{i{\mu }_0{\mu }_m{\tilde{K}}_{\varphi }\left[k_{-}^{\left(v\right)}\text{?}+k_{+}^{\left(v\right)}\text{?}\right]}{k_{+}^{\left(v\right)}{k}_{+}^{\left(s\right)}\text{?}-k_{-}^{\left(v\right)}{k}_{-}^{\left(s\right)}\text{?}},& \left(13\right)\end{array}$$\text{?}\text{indicates text missing or illegible when filed}$

[0108] In the example given here, we are interested in the solution in the static limit (i.e., k.sub.p>>|ϵ.sub.mμ.sub.m|k.sub.0, k.sub.ρ>>k.sub.0 and k.sub.ρ>>|ϵ.sub.s|k.sub.0), to discuss the excitation of magnetic surface plasmons. In this limit, considering a metamaterial with Re(μ)=−1 and low electromagnetic losses (i.e., μ.sub.m≅1+iIm(μ.sub.m) and Im(μ.sub.m)<<1), the preceding relationships are reduced to:

[00017]$\begin{array}{cc}\text{?}\simeq \frac{2{\mu }_0{\mu }_m{\tilde{K}}_{\varphi }\text{?}}{k_{\rho }\left[{\left(1+p_m\right)}^2\text{?}-4\text{?}\right]},C_m\simeq \frac{{\mu }_0{\mu }_m\left(1+{\mu }_m\right){\tilde{K}}_{\varphi }\text{?}}{k_{\rho }\left[{\left(1+p_m\right)}^2\text{?}-4\text{?}\right]},F_m\simeq \frac{2{\mu }_0{\mu }_m{\tilde{K}}_{\varphi }\text{?}}{k_{\rho }[{\left(1+p_m\right)}^2\text{?}-4\text{?}},F_s\simeq \frac{{\mu }_0{\mu }_m{\tilde{K}}_{\varphi }\left[2\text{?}+\left(1+{\mu }_m\right)\text{?}\right]}{k_p\left[{\left(1+{\mu }_m\right)}^2\text{?}-4\text{?}\right]},& \left(14\right)\end{array}$$\text{?}\text{indicates text missing or illegible when filed}$

[0109] The expressions obtained highlight two relevant regimes, namely the Pendry regime for k.sub.ρ<<k.sub.1 and the plasmonic regime for k.sub.ρ>>k.sub.1, being k.sub.1=Im.sup.−1 log [2/Im(μ.sub.m)]. In the Pendry regime, when the support of {circumflex over (K)}.sub.ϕ is in the region k.sub.ρ<<k.sub.1, the potential vector within the sample (for z>0) is given by expression (3).

[0110] On the contrary, in the plasmonic regime, when the support of K.sub.ϕ is in the region k.sub.ρ>>k.sub.1, the potential vector, within the sample, is given by expression (4).

[0111] From the comparison of Eq. (3) and the Eq. (4), the resonant nature of the solution in the plasmonic regime is apparent: |1+μ.sub.m=Im(μ.sub.m) and, as Im(μ.sub.m) decreases, (4) shows a divergent trend.

[0112] The data above are provided as examples. It is worth noting that in general, in addition to the spatial arrangement of sample, coil, and metamaterial, it is sufficient to obtain the improvement effect of the invention that the metamaterial is chosen so that it is adapted to develop a surface plasmonic regime, the rest of the values of the parameters being related to optimized configurations of the basic concept of the invention.

[0113] Although the examples given refer to magnetic surface plasmons, the technical concept of the invention is also applicable to electric surface plasmons, as described above.

IV. Further Embodiment

[0114] According to the invention, an apparatus for the nuclear magnetic resonance analysis of a sample containing at least one nucleus of interest, comprising means of producing a static magnetic field, at least one induction coil C with a maximum transverse dimension ρ.sub.0 and tuned in a pass-band around the Larmor frequency defined on the basis of said static magnetic field and at least one nucleus of interest, at least one metamaterial MM, and at least one sample S housing. In the apparatus: [0115] said at least one induction coil C is inserted between said at least one metamaterial MM and said at least one sample S housing; [0116] the distance d.sub.m between said at least one metamaterial MM and said at least one. Induction coil C is in the range from 0 to the maximum transverse dimension of the induction coil; and [0117] the metamaterial (MM) is chosen so that it is capable of developing a magnetic or electric surface plasmonic regime;

[0118] According to an aspect of the invention, the distance d.sub.m is between 0 and 1/10 of the maximum transverse dimension of the induction coil.

[0119] According to a different aspect of the invention, the distance d.sub.s between said at least one induction coil C and said at least one sample S housing is in the range from 0 to the maximum transverse dimension of the induction coil. The distance d.sub.s can be comprised between 0 and 1 cm.

[0120] According to an aspect of the invention, said at least one metamaterial MM is a slab with two opposite faces (e.g. lying substantially on said first and second plane), wherein the minimum radius of curvature of at least one of the two opposite sides of the slab is greater than the maximum transverse dimension of the. Induction coil C.

[0121] According to a different aspect of the invention, said at least one metamaterial MM is characterized by a relative magnetic permeability pi such that Re(μ.sub.m) is in a range about the value −1, said range having a width equal to 2.Math.Im(μ.sub.m). Preferably: the at least one metamaterial MM has a thickness l.sub.m between the two opposite faces such that l.sub.m> 1/10 of the maximum transverse dimension of the induction coil C; the metamaterial MM has a relative magnetic permeability μ.sub.m; the maximum transverse dimension of the coil ρ.sub.0<2π/[l.sub.m.sup.−1 log(2/Im(μ.sub.m))]; and the condition that Re(μ.sub.m) is in an amplitude range of 2-Im(μ.sub.m) about the value −1 holds.

[0122] Said at least one metamaterial MM and said at least one induction coil C can have a development substantially along their respective concentric arcs of circumference. Preferably, the respective concentric arcs of circumference are arcs of 360°. According to another aspect of the invention, said at least one metamaterial MM and/or said at least one induction means C respectively consist of a plurality of metamaterials MM #1, MM #2, MM #3 and/or dielectric materials and induction means C#1, C#2, C#3, positioned in consecutive and separate portions of the respective arcs of circumference.

[0123] According to the invention, said at least one metamaterial MM can be characterized by a relative magnetic permeability μ.sub.m such that Im(μ.sub.m) is less than 10.sup.−1, preferably (μ.sub.m) is less than 10.sup.−2 or 10.sup.−3.

[0124] According to an aspect of the invention, said at least one metamaterial MM is chosen so as to present at least two poles tuned to two different Larmor frequencies of at least two corresponding nuclei of interest.

Embodiment with Sphere

I. MLSPS Excitations and Improved Signal-to-Noise Ratio

[0125] In the sections above (see also Ref. [3]), the inventors suggested excited magnetic surface plasmons on the surface of a negative permeability MM slab to increase the SNR values of the magnetic resonance. It is worth considering that surface plasmon polaritons (SPP) and magnetic and/or electrical surface plasmons may exist in geometries other than the slab (e.g, particles with dimensions below the wavelength or empty cavities with different topologies) and can be applied in the magnetic resonance according to the invention. Here, for example, we will discuss the existence of magnetic localized surface plasmons (MLSPs), hosted by a sphere (of radius ρ.sub.m), which in reference to the previous embodiments can be identified as l.sub.m/2; or a spheroid with two semi-axes) of MM with negative permeability. Exploiting both the spherical symmetry of the MM device considered and the rotational invariance relative to the axis z of the apparatus shown in FIG. 22, we will focus our attention on monochromatic solutions of the form Ā=Re[A.sub.ϕ(r,θ)e.sup.−iω] with angular frequency a and where A, A.sub.ϕ{circumflex over (ϕ)} is the azimuthal component of the electromagnetic vector potential, r=√{square root over (ρ.sup.2+z.sup.2)} and θ=arccos (z/r)(see FIG. 22). From Maxwell's equations, we can obtain

∇.sup.2A.sub.ϕ+μ.sup.−1∇μ×(∇×A.sub.ϕ)+εμk.sub.0.sup.2A.sub.ϕ=0  (15)

where ε and μ are dielectric permittivity and magnetic permeability, respectively, and k.sub.0=ω/c (c is the speed of radiation in vacuum). We will assume a homogeneous magnetic MM sphere (with radius ρ.sub.m) with relative permeability and permittivity μ=μ.sub.m, ε.sub.m=1 within the sphere and μ=1, ε=1 otherwise. Following Mie's approach, considering the expansion in spherical waves and imposing the connection conditions on the surface of the sphere, it results:

[00018]$\begin{array}{cc}\begin{array}{cc}{A}_{\varphi }=A_L\text{?}{j}_L\left(k_{m}r\right)P_L^{\left(1\right)}\left(\cos \theta \right)& \mathrm{for}r\le {\rho }_m,\\ A_{\varphi }=A_L{h}_L^{(+)}\left(k_{0}r\right)P_L^{\left(1\right)}\left(\cos \theta \right)& \mathrm{for}r>{\rho }_m,\end{array}& \left(16\right)\end{array}$$\text{?}\text{indicates text missing or illegible when filed}$

where A.sub.L is a constant, k.sub.m=√{square root over (ϵ.sub.mμ.sub.m)}k.sub.0, L a positive integer (L=1, 2, 3, . . . ), P.sub.L.sup.(1) is the Legendre polynomial P.sub.L.sup.(m) with m=1, j.sub.L the spherical Bessel functions and .sub.L.sup.(+) the output spherical Hankel functions. These solutions represent localized magnetic waves characterized by the dispersion relation

φ.sub.L.sup.(1)(k.sub.mρ.sub.m)−μφ.sub.L.sup.(+)(k.sub.0ρ.sub.m)=0,  (17)

wherein:

φ.sub.L.sup.(+)(ξ)=(d[ξh.sub.L.sup.(+)(ξ)]/dξ)/h.sub.L.sup.(+)(ξ), φ.sub.L.sup.(1)(ξ)=(d[ξj.sub.L(ξ)]/dξ)/j.sub.L(ξ)   (18)

with ξ=k.sub.mρ.sub.m.

[0126] To physically grasp the main features of these solutions, we will consider the static limit k.sub.0.fwdarw.0 where

[00019]$\begin{array}{cc}\begin{array}{cc}{A}_{\varphi }\simeq {A_L\left(\frac{r}{{\rho }_m}\right)}^L{P}_L^{\left(1\right)}\left(\cos \theta \right)& \mathrm{for}r\le {\rho }_m,\\ A_{\varphi }\simeq {A_L\left(\frac{{\rho }_m}{r}\right)}^{L+1}{P}_L^{\left(1\right)}\left(\cos \theta \right)& \mathrm{for}r>{\rho }_m,\end{array}& \left(19\right)\end{array}$

and the dispersion relation Eq. (17) becomes

[00020]$\begin{array}{cc}{\rho }_m=\frac{1+L}{L}.& \left(20\right)\end{array}$

[0127] It is worth noting that fora specific L and, therefore, a specific value of μ.sub.m, the second equation of the Eq. (19) coincides with a term of the standard multipole expansion. Equation (20) is the magnetic counterpart of the condition of the existence of electric localized surface plasmons [1] and makes these resonances exist only for discrete magnetic permeability values. It is worth noting that the excitation of an electromagnetic surface mode generally shows a resonant behavior [1], so an adequate MLSP excitation can produce a significant improvement in the RF electromagnetic field.

[0128] The improving effect obtained by using a sphere of MM applies to any value of the sphere radius ρ.sub.m once the μ.sub.m of the sphere is chosen according to one of the values determined by the equation (20) which is valid in the case of an isolated sphere, or by means of numerical simulations if the presence of the sample S and the RF coil C and/or in the case of the spheroid are to be taken into account.

[0129] Preferably, the metamaterial MM with spherical shape has a relative magnetic permeability μ.sub.m such that Im(μ.sub.m) is less than 0.2, even more preferably less than 0.1.

[0130] To numerically test the improvement of the electromagnetic field due to the excitation of these surface plasmons located in a magnetic resonance configuration, we will consider the case in which a surface RF coil is located between the MM sphere with negative permeability and the sample, as shown in FIG. 22. Exploiting the rotational invariance of the setup along the z-axis, we evaluate the electromagnetic field and SNR using full-wave 2D simulations in cylindrical coordinates (i.e. In the plane (ρ, z)). In the numerical examples, we will set the frequency v=127.74 MHz (corresponding to a static magnetic field |B.sub.0|=3 T), ρ.sub.m=8.4 cm, d.sub.m=0 mm, l.sub.m=l.sub.s=3ρ.sub.m, d.sub.x=2 mm. The electromagnetic response of the sample is that of muscle tissue (ϵ.sub.s=63.5+i101.2). The RF surface coil has negligible thickness along the z-axis and is described by the azimuth current density J.sub.ϕ(ρ, z)=K.sub.ϕ(ρ)δ(z) (δ(⋅) is the Dirac delta function), where K(ρ) K.sub.0 for ρ.sub.0−w/2<ρ<ρ.sub.0+w/2 and K.sub.ϕ(ρ)=0 otherwise (ρ.sub.0=ρ.sub.m/2=4.2 cm, w=ρ.sub.m/10=ρ.sub.0/5=8.4 Im and K.sub.0 is a constant chosen to obtain a unit current in the coil).

[0131] In FIG. 23, we trace |B.sub.1()|/ρ.sub.0 (continuous simple line) and SNR.sup.(n) (continuous gray line with star symbol) at the spatial point ρ=0 mm and z=6 mm (on the z-axis of the RF col) depending on the real part of the permeability μ.sub.m of the metamaterial, assuming Im (μ.sub.m)=0.1. For comparison, the permeability values obtained from Eq. (20) are shown in the same figure as vertical black dotted lines. It is worth noting that both |B.sub.1.sup.(−)|/ρ.sub.0 and SNR.sup.(n) show several peaks (for the latter, highlighted by the star markers in FIG. 23) and each peak is due to the excitation of a specific MLSP. On the other hand, in FIG. 23 you can see that the peaks are shifted relative to the MLSP existence conditions provided by Eq. (20). This can be explained because Eq. (20) neglects the delay effects in Maxwell's equations and applies in the absence of both the sample and the RF surface coil. However, FIG. 23 clearly demonstrates that MLSPs are excited and support a significant improvement of the RF signal |B.sub.1.sup.(−)|/ρ.sub.0 and SNR.sup.(n).

[0132] Hereinafter, we focus our attention on the values Re(μ.sub.m) highlighted by the star markers in close correspondence with cases in which Re(μ.sub.m)=−1.39; −1.26; −1.20; −1.16; −1.13. From the comparison of the spatial distribution of the analytical solutions of Eq. (19) with the numerical results, it is apparent that the values Re(μ.sub.m) highlighted by the star markers in FIG. 23 correspond to the MLSPs of the Eq. (19) with L=3, 4, 5, 6, 7. Both the resonant behavior of MLSP excitations and their multipolar structure can be exploited to improve the SNR of the magnetic resonance. The first can be improved by reducing MM losses as previously demonstrated for planar configuration (above and [8]). In the case of the MM sphere, by choosing a specific permeability value μ.sub.m, the desired L mode can be excited and then, by increasing the value by L, narrower spatial confinement of the magnetic field and its intensity are obtained compared to the standard case of a surface RF coil in which the field has a dipolar distribution. In FIG. 24, for the sphere with the above geometry and Im (μ)=0.01, we report |B.sub.1.sup.(−)|/ρ.sub.0 (panel a) and SNR.sup.(n) (panel b) along the axis z and within the sample (i.e. ρ=0 cm and z>0.2 cm) for the permeability values corresponding to the resonance modes marked by the star symbols in FIG. 23. In FIG. 24 it can be seen that the greater Re (μ.sub.m) the greater the values of |B.sub.1.sup.(−)|/ρ.sub.0 and SNR.sup.(n); for Re (μ.sub.m)=−1.13, the SNR.sup.(n) near the sample interface (z=0.2 cm) is ≅10, FIG. 25, we report the spatial distribution of SNR.sup.(n) without (panel a) and with the sphere MM having the magnetic permeability values as in FIG. 24 (panels b-f), where the dashed isolevel lines correspond to SNR.sup.(n)=1. For the considered system, an increase of SNR in the sample (e.g. SNR.sup.(n)>1) is obtained within a region with a longitudinal dimension (z-axis) of about 5 cm and a transverse dimension (radius ρ) of about 3.5 cm.

[0133] FIG. 25 shows SNR maps.sup.(n) in the presence of the MM sphere compared to the SNR map obtained in the standard configuration. It is worth noting that for the various modes, obtained with negative values of μ.sub.m, the value of SNR.sup.(n) increases up to about 10 times with a clear application advantage in the receiving phase of the magnetic resonance experiment. In FIG. 26 a similar advantage is observed for the RF excitation field, which implies advantages in RF pulse duration and/or maximum RF amplification power.

[0134] Furthermore, FIG. 27(a) shows a maximum centered around the position of the RF coil (ρ.sub.0=42 cm). In FIG. 27 (b)-(f) it is observed that the presence of the sphere of MM (Re(μ.sub.m)=−1.39, −1.26, −1.20, −1.16, −1.13) introduces an asymmetric distribution of the electric field with respect to the plane z=0, concentrating it more inside the sphere MM (values of z<0), near its surface, and shifting the maximum electric field towards smaller radial positions (FIG. 27(f), ρ=2 cm).

II. MLSPS Mimicking by a Dielectric

[0135] In the previous section, we studied and characterized MLSP hosted by a spherical MM with negative permeability, suitably inserted in a magnetic resonance configuration. As a matter of fact, the desired magnetic behavior (i.e. a resonant magnetic response at Larmor frequency and a negligible magnetic response at the static limit [3]) is not available in nature. However, a specific magnetic response can be achieved by means of an appropriate composite structure. For example, Freire et al. [6] made a slab of MM having ρ=−1 In the RF field with a periodic ring resonator structure [6]. The use of such repeated structures makes the manufacture of such devices complex. Furthermore, their theoretical description, based on effective medium theories, has imitations due to the intrinsic uneven response of such materials on scales comparable with those of their composite structure.

[0136] According to the invention, these limits can be exceeded by demonstrating that the electromagnetic field generated by MLSPs outside the sphere can be mimicked using dielectrics with preferably high relative dielectric constant (typical values of 100-4000 at the frequencies of the previous example are provided in the literature) already available in nature [5]. It is worth noting that several research teams have studied the inclusion of high-ε dielectric materials in a standard magnetic resonance scanner to manipulate the local RF field distribution [5]. Such materials support intense displacement currents capable of modifying the RF field distribution and this effect was taken into account for the impedance adaption, shimming, and focusing the RF field distribution to different static field values (3, 4, 7, 9, 4 T).

[0137] In this invention, on the contrary, we will show in detail, by way of example, the equivalence (mimicking), with good approximation, between the external scattering field of a homogeneous dielectric sphere with high permittivity and the electromagnetic field of a specific MLSP produced by a MM sphere of the same radius outside it. We will demonstrate that this dielectric sphere. In turn, produces the significant magnetic resonance SNR enhancement we have already shown for the MM sphere.

[0138] For this purpose, referring to FIG. 28, to test the electromagnetic equivalence between the isolated sphere of MM and a dielectric sphere with the same radius ρ.sub.m, we compare the localized waves for both configurations. Considering a dielectric sphere in vacuum the electromagnetic response of which is described by the permittivity ε.sub.d and permeability μ.sub.d=1, the condition of existence and the distribution of the electromagnetic vector potential are given by Eq. (17) and Eq. (16) (replacing μ.sub.m con μ.sub.d=1, k.sub.m con k.sub.d=√{square root over (ϵ.sub.d)}k.sub.0), respectively. From Eq. (16),It is apparent that the resonant surface mode of order L (one hosted by the dielectric sphere and the other by the MM sphere with negative permeability) have the same electromagnetic field profile outside the sphere while they differ inside. Equivalence is guaranteed by the condition of existence for both modes on the surface, i.e.

φ.sub.L.sup.(1)(k.sub.mρ.sub.m)−μ.sub.mφ.sub.L.sup.(+)(k.sub.0ρ.sub.m)=0,

φ.sub.L.sup.(1)(k.sub.dρ.sub.m)−φ.sub.L.sup.(+)(k.sub.0ρ.sub.m)=0.  (21)

Consequently, the equivalent dielectric permittivity ϵ.sub.d may be evaluated by solving the complex transcendent equation:

φ.sub.L.sup.(1)(k.sub.mρ.sub.m)−μ.sub.mφ.sub.L.sup.(1)(k.sub.dρ.sub.m)=0.  (22)

A similar approach was initially considered by Devilez et al. to mimic surface electric plasmons hosted by a spherical metal particle by means of a spherical dielectric particle [4]. Eq. (21) can only be met exactly for real permeability and permittivity values, i.e. for loss-free materials. For a low-loss magnetic material (the permeability value of which satisfies by first approximation the first of the Eq. (21)), we can still determine an equivalent complex permittivity which satisfies Eq. (22) and the accuracy of the electromagnetic equivalence can be verified by means of a numerical simulation taking into account all the implementation parameters.

[0139] The value of ϵ.sub.d determined by equation (22), with the parameter k.sub.0 implicitly contained in the equation by means of the definition of k.sub.m and k.sub.d, depends on the chosen working frequency.

[0140] Therefore, the resulting value of ϵ.sub.d will depend on the selected value of μ.sub.m of the MM sphere whose electromagnetic field one wishes to mimic, the radius of the sphere (or the radii for the spheroid) and the working frequency.

[0141] In the presence of magnetic metamaterial losses and/or in the presence of sample S and (RF) coil C the solution of equation (22) no longer guarantees the exact correspondence between the electromagnetic field generated outside the magnetic MM sphere and outside the uHDC sphere. In this case, the verification of the accuracy of the approximation between the electromagnetic fields must be performed by numerical methods, as shown in the example of FIG. 33(a).

[0142] If the accuracy of the solution found by means of the equation (22) is deemed not satisfactory, it can be improved by numerical methods by determining complex values of ϵ.sub.d that minimize the differences between the electromagnetic fields of the MM sphere and the uHDC sphere within the sample S.

[0143] For the frequencies of interest, we consider all those of use in MRI/NMR/EPR ranging from 1 kHz to 300 GHz. In the MRI scope, we expect the range of values of the radius of the sphere of MM or the equivalent sphere of uHDC that have a practical utility to be comprised between 0 and 20 cm. In the MRI scope for frequencies close to 400 MHz, the preferred values of ϵ.sub.d of the uHDC sphere (spheroid) would be Re(ϵ.sub.d) about 8000, Im(ϵ.sub.d) less than 300 (i.e. tan δ less then 0.038).

[0144] By way of example, setting L=5 and numerically searching for solutions of Eq. (22) with ρ.sub.m=8.4 cm, v=127.74 MHz (3 T) and μ.sub.m=−1.20+i 0.01, it is obtained that Eq. (22) is satisfied for ϵ.sub.d=1324+i 1.65. This result suggests that the MLSP considered with μ.sub.m=−1.2+i 0.01 (the permeability value ensures an SNR improvement as shown in FIG. 33 and almost satisfies the first Eq. (20)) is reproduced by the electromagnetic field outside the dielectric sphere with the same radius and ϵ.sub.d=1324+i 1.65. Of practical relevance is what happens when the sample and the RF coil are in the immediate vicinity of the sphere.

[0145] FIG. 30 shows the mapping in the (ρ, z) plane of the SAR at 127.74 MHz (B.sub.0=3 T) with the same parameters as FIG. 29 without (a) and with (b) the uHOC sphere for tan δ=0.04 and ρ.sub.m=8.82 cm. Similarly, FIG. 31 shows the mapping in the plane (ρ,z) of the effective transmission field |B.sub.1,eff.sup.(+)|(ϕ=π/2) at 127.74 MHz (B.sub.0=3 T) with the same parameters as FIG. 29 without (a) and with (b) the uHDC sphere for tan δ=0.04 and ρ.sub.m=8.82 cm.

[0146] FIG. 32 shows results for a uHDC sphere with Re (ϵ.sub.d)=3300 and radii ρ.sub.m=4.08, 7.49, 10.64 cm supporting resonances L=1, 3, S, respectively, at the Larmor frequency of 63.87 MHz (B.sub.0=1.5 T). Again in this case, a gain of SNR.sup.(n) (estimated at z=2 mm, ρ=0 mm) is also observed in this case, but at a narrower range of values of the loss tangent tan δ (gray area of FIG. 32 b).

[0147] In FIG. 33 we compare the configuration of FIG. 22 (magnetic MM sphere with μ.sub.m=−1.20+i 0.01) and that of FIG. 28 (uHDC sphere with ϵ.sub.d=1324+i 1.65) both corresponding to the resonant mode L=5. In FIGS. 33 (a) and 33 (b), within the sample (z >0.2 cm, φ=π/2), we compare |B.sub.1.sup.(−)|/μ.sub.0 and SNR.sup.(n) along the a-axis to different values of the radial coordinate; while, in FIGS. 33 (c) and 33 (d), we will compare the maps of SNR.sup.(n) in the plane (ρ, z)(for φ=π/2). From FIG. 33, it is apparent that, in a realistic MRI configuration, the equivalence is valid to a good extent because the deviation is closely located on the z axis and near the interface between the sample and the air. Indeed, we note that the maximum deviation for |B.sub.1.sup.(−)/μ.sub.0 is about 24% at z=2 mm and ρ=0 mm, while the correspondence is more and more accurate in the other regions.

[0148] It is worth noting that the mimicking is more accurate when losses are lower. In the limit of the absence of losses, the shape of the field is dominated by the divergent displacement (magnetization) currents within the dielectric (magnetic) sphere, which become very large compared to those of the coil and of the currents in the sample, and consequently, we approach the condition in which the spheres are isolated and Eq. (20) can be fully satisfied.

III. MIE Resonances with Very High Permittivity Ceramics

[0149] To verify the feasibility of the suggested configuration in which the magnetic sphere MM is replaced by the dielectric sphere, we will study the Mie resonances and their effects on magnetic resonance applications by assuming the properties of dielectric materials already used in the context of nuclear magnetic resonance. We will not discuss the quality of the mimicking approach, related to material losses, sample, and RF coil presence hereinafter. Here, we will focus our analysis on the improvement in MRI performance achievable with the inclusion of a dielectric sphere when its radius is chosen to satisfy the second of Eq. (20) at the desired Larmor frequency.

[0150] A large class of ferroelectric materials has low losses and has a very high real part of dielectric permittivity, with values which can be customized using different physical-chemical factors [5] (e.g applied static electric field, temperature, chemical composition, doping and mixing with other dielectrics). However, the desired dielectric permittivity value may not be easily obtained at the operating frequency. On the other hand, it is worth noting that MLSP resonances are also highly dependent on material geometry. Indeed, by assigning a specific value of the dielectric constant, it is possible to satisfy the condition of existence by finely adjusting the radius of the sphere ρ.sub.m. As mentioned above, dielectric ceramics have been used to improve the different aspects of the magnetic resonance. High dielectric permittivity values ϵ.sub.d were made from high-concentration aqueous ceramic mixtures (Re (ϵ.sub.d)=475 at 7 T) or sintered ceramic beads (Re (ϵ.sub.d)=515 at 3 T). Rupprecht et al. [5] demonstrated improved RF coil sensitivity using materials with an ultra-high dielectric constant (uHOC) at 1.5 T and 3 T. In particular, they experimentally studied lead zirconate titanate-based ceramics (PZT) where Re (ϵ.sub.d)=1200 or Re (co)=3300 at 3 T and 1.5 T, respectively. Recently, to increase the SNR of the magnetic resonance, the use of ceramic materials was suggested, based on BaTiO.sub.3 with ZrO.sub.2 and CeO.sub.2 as additives, leading to uHOC with Re (ϵ.sub.d)=4500 at 1.5 T. Here we will study the performance in magnetic resonance considering the two permittivity values reported by Rupprecht et al. [5].

[0151] In the first example, we will fix the real part of the sphere permittivity to the value of Re (ϵ.sub.d)=1200 for the working frequency 127.74 MHz (|B.sub.0|=3 T) and, to study the effect of dielectric losses on magnetic resonance performance, we will vary the imaginary part of the permittivity. For this purpose, full-wave numerical simulations were performed, using axial symmetry again, choosing the same coil and the example parameters in FIG. 33 (i.e., d.sub.m=0 mm, ρ.sub.s=l.sub.s=25.2 cm, d.sub.s=2 mm, ρ.sub.0=4.2 cm, w=8.4 mm e ε.sub.s=63.5+i 101.2), except for the radius of the sphere ρ.sub.m adjusted to select three different resonant modes (i.e. L=1, 3, 5).

[0152] An additional example is shown in the results of FIG. 29, wherein we have a sphere of ultra-high dielectric constant (uHOC) with Re (ϵ.sub.d)=1200 and radii ρ.sub.m=3.38, 6.21, 8.82 cm supporting resonances L=1, 3, 5, respectively, at the Larmor frequency of 127.74 MHz (B.sub.0=3 T) and in the presence of RF coil and sample. Also in this case, a gain of SNR.sup.(n) is observed (estimated at point z=2 mm, ρ=0 mm) at a wide range of values of the loss tangent tan δ of the uHOC sphere (gray area of FIG. 29 b).

[0153] In FIG. 29, we report |B.sub.1.sup.(−)/ρ.sub.0 and SNR.sup.(n) at the spatial point z=2 mm, ρ=0 mm, as a function of the dielectric losses parameterized by tan δ=Im (ε.sub.d)/Re (ε.sub.d), for ρ.sub.m=3.38 cm, ρ.sub.m=6.21 cm, ρ.sub.m=8.82 cm (dielectric sphere radii supporting MLSP with L=1, 3, 5, respectively). The range of losses considered (5-10.sup.−3<tan δ<0.17) has been selected in accordance with the literature [5] for the frequencies corresponding to magnetic resonance imaging at fields of 3 T or less. From FIG. 29 (a), the RF signal enhancement produced by the uHOC sphere is apparent. Furthermore, in FIG. 29(b), we can observe an increase in SNR (SNR.sup.(n)>1) in the wide range 0.005<tan δ<0.167 (gray region). In case tan δ=0.04, i.e. the value of the MRI dielectric pod tested in [5], we observe SNR.sup.(n)=1.6, 1.5, 1.3 for ρ.sub.m=3.38, 6.21, 8.82 cm, respectively. From the results shown in FIG. 29 (b), using a material with tan δ=−5.Math.10.sup.−3, the SNR.sup.(n) would be 2.7 (ρ.sub.m=3.38 cm, L=1), 3.3 (ρ.sub.m=6.21 cm, L=3) and 3.1 (ρ.sub.m=8.82 cm, L=S), respectively.

[0154] In a second series of full-wave simulations, we assume the real part of the permittivity of the uHDC sphere Re ( )=3300 and the working frequencyv=63.87 MHz (CDIN=1.5 T). We will consider the same coil and geometric parameters as in FIG. 29, choosing ε.sub.s=72.3+i 193.7 corresponding to the dielectric constant of muscle tissue at 1.5 T. The MLSP with L=1, 3, 5 correspond to ρm=4.08, 7.49, 10.64 cm, respectively. FIG. 32 shows |B.sub.1.sup.(−)|/μ.sub.0 (a) and SNR .sup.(N) (b) at the spatial point z=2 mm, ρ=0 mm as a function of tan δ for the chosen modes. FIG. 32(a) shows an enhancement of the receiving RF signal |B.sub.1.sup.(−)|/μ.sub.0 throughout the range 5.Math.10.sup.−3<tan δ<0.09. In FIG. 32(b), we note that SNR.sup.(n)>1 (region in gray) for tan δ<0022. For tan δ=0.005, we get an SNR.sup.(n) of about 1.4, 1.6, 1.5 for ρ.sub.m=a 4.08, 7.49, 10.64 cm, respectively. However, for tan δ=0.05, as in the material wed previously, we have SNR.sup.(n)=0.74, 0.69, 0.62 for ρ.sub.m=4.08, 7.49, 10.64 cm, respectively. As in many photonic sub-wavelength devices, tan δ is a crucial parameter because high losses can drastically reduce or even eliminate the electromagnetic resonance of the dielectric.

[0155] For the sake of completeness, we will evaluate the SAR and transmission efficiency within the sample. The local specific absorption rate is given by SAR σ|E|.sup.2/(2ρ.sub.v),where E is the complex electric field amplitude, σ=ωIm() and ρ.sub.v are the electrical conductivity and mass density of the sample, respectively [5]. In FIG. 30, we compare the SAR for the unit current without (a) and with (b) the 3 T uHDC sphere using the same parameters as FIG. 29 with tan δ=0.04 and ρ.sub.v=3490 kg/m.sup.3 [5]. Clearly, the maximum SAR (e.g. SAR.sub.max=max(SAR)) is a critical parameter because it limits the maximum power to be applied to the drive RF coil. Here, it is very interesting to note that the SAR.sub.max and the SAR averaged over the whole volume (e.g. SAR.sub.a=∫.sub.sampledSAR/V, where V is the whole sample volume) are both reduced in presence of the dielectric sphere. More precisely, SAR.sub.max (SAR) decreases from 21.5 W/kg (3.5.Math.10.sup.−2 W/kg) to 18.6 W/kg (2.1.Math.10.sup.−2 W/kg) without and with the uHDC sphere, respectively. The reduction of SAR.sub.max by approximately 14% (SAR.sub.a reduced by 40%) in the presence of the uHDC sphere is an important advantage for 3 T magnetic resonance and could be useful for higher static field applications (7; 9.4 T).

[0156] In FIG. 31, we compare the maps of |B.sub.1/.sup.(+)|, defined as the ratio of the absolute value of the field B.sub.1.sup.(+) and the square root of SAR.sub.max, without (a) and with (b) the uHDC sphere under the same conditions as in FIG. 30. Despite the fact that the geometrical parameters of the considered configuration have not been completely optimized, from FIG. 31 both an improvement in RF efficiency and a significant focus of the magnetic field in the region near the axis of the RF coil, i.e. near the central volume of the sample under study, are apparent. Finally, we can observe, in the presence of the uHDC sphere, that the SAR is concentrated at about ρ=4 cm, i.e. close to the RF coil. As a result, by placing a relatively small sample close to the coil axis, our configuration makes it possible to improve the magnetic resonance performance by reducing SAR in the region of interest.

[0157] FIG. 34 shows the layout of an MRI configuration with geometry B according to an aspect of the invention, with the sphere (of magnetic MM or UHDC) positioned between a standard surface RF coil and the cylindrical sample.

[0158] FIG. 35 shows the maps |B.sub.1.sup.(−)/μ.sub.0 within the cylindrical sample (d.sub.m=0 mm; z≥2 mm) in the presence of: magnetic MM sphere (μ.sub.m=−1.2+i 0.02) with (a) geometry A (d.sub.s=2 mm) or (b) geometry B (d.sub.s=ρ.sub.m+2 mm); uHDC sphere (ε.sub.d=1324+i 1.65) with (c) geometry A (d.sub.s=2 mm) or (d) geometry B (d.sub.s=ρ.sub.m+2 mm), FIG. 35 shows, for a more immediate comparison, the profile (ρ=0 mm) of |B.sub.1.sup.(−)/μ.sub.0 for geometry A and B in the presence of the MM sphere (e) or of the uHDC sphere (f).

[0159] FIG. 36 shows the graphs as in FIG. 35 for the |E| field (per current unit).

[0160] FIG. 37 shows the graphs as in FIG. 35 for the SAR (for the unit current).

[0161] FIG. 38, assuming a working frequency of 127.74 MHz (B.sub.0=3 T), shows: in (a) the profile of |B.sub.1.sup.(−)|/μ.sub.0 along the z-axis (d.sub.m=0 mm) for the uHDC sphere (a=1200+i 48) with geometry A (d.sub.x=2 mm, solid line) and geometry B (d.sub.s=μ.sub.m+2 mm, dashed line); in (b) the specific SAR absorption rate at the plane (ρ, z) without uHDC; in (c) the SAR in the presence of the uHDC sphere (ε.sub.d=1200+i 48) with geometry A (d.sub.s=2 mm); in (d) the SAR in the presence of the uHDC sphere (ε.sub.d=1200+i 48) with geometry B (d.sub.s=ρ.sub.m+2 mm).

[0162] FIG. 35 shows the maps of |B.sub.1.sup.(−)/μ.sub.0 in the plane (ρ, z) for the MM sphere and the uHDC sphere. A reasonable similarity of spatial distribution between geometric configurations A (FIG. 1 a) and B (FIG. 34) is observed, with a slight decrease in amplitude in the case of the uHDC sphere in geometry B. Similar results are observed for the field maps |E| shown in FIG. 36 and SAR shown in FIG. 37. The results shown in FIG. 35-37 allow us to conclude that the uHDC sphere is a valid alternative to the use of the MM sphere, with a considerable practical simplification.

[0163] Although the case in which the induction coil is adjacent to or away from one end of the metamaterial or dielectric along dimension z was always treated above, it is also possible for the coil to surround at least part of the metamaterial or dielectric. In other words, two parallel planes can be defined between which the metamaterial or dielectric extends, the planes being parallel and perpendicular to the z-direction, in such a case, the distance of the coil from either plane can be both positive and negative. In case of negative distance, the plane of the coil crosses somewhere through the metamaterial or dielectric, and obviously, the coil must be wide enough to surround it on the xy plane, so that there is no interpenetration between the two elements.

[0164] It is also possible to express this configuration by saying that the distance module d.sub.m is comprised in the range specified below. The possibility of using a positive or negative distance (coil between the two planes above or outside the metamaterial or dielectric) depends on the geometry of the metamaterial or dielectric as well as on the plasmonic or dielectric resonance regime to be excited. A positive distance is, however, generally preferred.

[0165] Even more in general, the various configurations of the apparatus according to the invention, in terms of the aforesaid distances, can be included in the relation d.sub.s+d.sub.m≥0. A particular case of the invention is when both d.sub.s and d.sub.m are positive.

Advantages of the Invention

[0166] According to the invention, a new use of a magnetic metamaterial slab is provided to increase the performance of an RF coil in a magnetic resonance device useful for both spectroscopy (NMR) and imaging (MRI) applications. The approach of the invention is based on magnetic plasmonic resonances present on at least one surface of a metamaterial slab with Re (μ.sub.m)=−1 which are responsible for a strong increase of the RF magnetic field within a sample suited for magnetic resonance imaging. A further advantage of the suggested configuration is the positioning of the metamaterial slab, i.e., outside the RF coil and sample assembly, in a region in which free space is usually available.

[0167] In this respect, the present invention has the potential to be applied in most current situations of use with minimal additional requirements compared to available configurations. The results are based on an approximate description of the current density in the RF coil and do not assume losses in the RF coil itself. Furthermore, the described mode can be implemented also if one desires to detect the signal coming from two or more NMR or MRI active nuclear species present in the sample, i.e. In multi-nuclear mode, using a metamaterial able to support at least two distinct plasmonic resonances the resonance frequency of which coincides, or is close to, the one corresponding to the known Larmor frequencies (metamaterial chosen to present at least two poles tuned to two different Larmor frequencies of at least two corresponding nuclei of interest). A two-dimensional metamaterial configuration has been described in the literature which can be used to improve the detection of the proton .sup.1H and phosphorus .sup.31P nuclear signal. Such metamaterial supports Fabry-Perot resonances by means of a given number of metal strips appropriately separated from each other and arranged on a plane. Such device behaves like a set of electric dipoles, suited for the low frequencies corresponding to the signal of .sup.31P and a second set of magnetic dipoles necessary for the detection of the signal .sup.1H.

[0168] Finally, we can note that the invention could also be extended to the context of electronic paramagnetic resonance (EPR).

[0169] Furthermore, to use the prior art with dielectrics according to the invention,thevaluesoftherealandImaginarypartoftheelectricalpermittivityof the dielectric material should be appropriately selectable to satisfy the conditions of electromagnetic equivalence relative to the magnetic metamaterial of identical or similar geometry.

[0170] It Is interesting to observe the ability of the invention to replace a magnetic and/or electrical metamaterial with an equivalent dielectric material, because the practical making of the metamaterial may present limitations due to the physical dimensions of the constituent unitary cells (usually small inductive/capacitive resonant circuits of a circular shape, see [2]), which makes it difficult to achieve the spatial homogeneity condition.

[0171] More generally, the following beneficial effects of the invention are listed in a non-exhaustive manner [0172] 1. The metamaterial slab can support surface plasmonic resonances at the frequency of use of magnetic resonance (Larmor frequency) on at least one of its component surfaces. Such plasmonic resonances can be appropriately excited by an RF coil, tuned to the Larmor frequency of the magnetic resonance apparatus. Plasmonic resonances, characterized by the presence of intense concentrated currents near at least one of the surfaces of the metamaterial slab, have the effect of amplifying the intensity of the RF magnetic field in a specific region of the sample under examination, which is placed at a given distance from the surface of the metamaterial slab. [0173] 2. Plasmonic resonances useful for the purposes of the present invention can be located on the surface of structures other than the slab, such as a spherical shape [1], a semi-spherical shape, a cylindrical shape, an ellipsoidal shape, a toroidal shape, and even structures with an irregular surface [1, 2]. [0174] 3. The RF coil can also be used to detect the signal of the sample under examination which, in a similar manner as described in the preceding point, is amplified by the plasmonic resonances of the metamaterial. [0175] 4. The circular RF coil used in the resonance apparatus is described in FIG. 1, may be replaced by a square, rectangular, or triangular coil, or any other shape capable of exciting plasmonic resonances on at least one surface of the metamaterial. [0176] 5. The geometry and composition of the metamaterial can be appropriately chosen to generate a given spatial distribution of the RF field amplitude in the inner volume of the sample under examination. [0177] 6. The metamaterial is preferably, but not necessarily, positioned outside the RF excitation/detection coil facing the sample itself, to maximize the amplification effect. [0178] 7. The properties of the metamaterial (used for making the slab or other useful structures) can be adjusted to assume the desired value at the working frequency (Larmor frequency) for the specific application of magnetic resonance, e.g. the frequency of about 64 MHz could be chosen to detect the hydrogen signal (.sup.1H) present in the tissues when these are in the presence of a static magnetic field of 1.5 T. [0179] 8. The functionality of the metamaterial can only be used during the excitation operating phase, or only during both the excitation and signal detection phases. [0180] 9. The electrical and/or magnetic parameters of the metamaterial can be modified, even in dynamic mode, within a given range by means of an appropriate electrical and/or mechanical control to modulate the effects on the signal in a specific spatial position. [0181] 10. The geometric arrangement of the metamaterial relative to the RF coil and the sample can be modified within a given range of values by means of a mechanical control to modulate the effects on the signal also in dynamic mode. [0182] 11. The metamaterial can support more than one mode of surface plasmonic resonance (multi-nuclear mode), each corresponding to a distinct frequency able to excite and/or detect, either simultaneously or consecutively, the signal of at least two nuclear species useful for magnetic resonance, and by way of example we could consider hydrogen (.sup.1H) and sodium (.sup.23Na) of biological tissues exposed to the same static magnetic field. [0183] 12. The element comprising the metamaterial and its excitation/detection RF coil can be structured in a volume configuration (e.g. of the birdcage, or saddle, or TEM type), which surrounds and encloses all or part of the test sample. [0184] 13. The element which comprises the metamaterial and the respective excitation/detection coil can be replicated a given number of times (N), and be arranged near the sample to ensure multi-channel operation, with sequential or parallel acquisition both for single nucleus (e.g. .sup.1H) and multi-nuclear (e.g. .sup.1H and .sup.23Na). [0185] 14. The properties of the MM can be adjusted to allow paramagnetic electronic resonance (ESR, EPR) applications in a frequency range from radio frequencies to microwaves. [0186] 15. In the case of the magnetic MM sphere, there is an infinite number of resonance modes which can be excited and each of which corresponds to its own spatial trend of the transmission and/or reception electromagnetic field and SNR, which can be useful for specific applications, so the expert user can select them according to needs. [0187] 16. To use a specific magnetic resonance mode to the desired Larmor frequency, the geometry (sphere radius) and the value of pw of the sphere (negative) must be adapted. For this purpose, analytical and/or numerical electromagnetic simulation methods may be used to optimize such parameters. [0188] 17. The efficiency maps of the transmission RF magnetic field with the magnetic MM sphere show that there is an improvement in RF efficiency and also a significant focus of the magnetic field in the region near the axis of the RF coil, i.e. near the central volume of the sample under study. [0189] 18. Having demonstrated the electromagnetic equivalence between the magnetic MM sphere and a dielectric sphere (uHDC) of the same radius and with selected permittivity value, it follows that the preceding advantages in terms of excitation and/or detection field and/or SAR apply to the case of the dielectric sphere, this advantage being particularly important for 3 T magnetic resonance and useful applications in higher static fields can be expected (7; 9.4 T). [0190] 19. The use of the uHDC sphere simplifies the practical implementation of the detection system, as it is not necessary to build unit cells with conductive and insulated elements, with considerable cost savings. [0191] 20. A further advantage of the uHDC sphere is the absence of static magnetic field disturbance, which allows the acquisition of MRI data without the introduction of artifacts. [0192] 21. The suggested uHDC device makes it possible to avoid complex manufacturing procedures and the inhomogeneous response of the electromagnetic field present in a magnetic composite MM when the size of the constituent inclusions of the MM becomes comparable to the radius of the sphere or with the size of the plasmonic resonance modes because the intrinsic inhomogeneity of the MM can dramatically modify or even eliminate the presence of such modes, the existence of which is based on the effective medium theory. With the use of uHDC, this fundamental limit is completely overcome because the homogeneous macroscopic dielectrics do not present spatial inhomogeneity.

LITERATURE

[0193] [1] A. V. Zayats, I. I. Smolyaninov, A. A. Maradudin, Nano-optics of surface plasmon polaritons, Physics Reports, 408, 131-134 (2005). [0194] [2] N. Engheta and R. W. Ziolkowski, Engineering, physics, and applications of Metamaterials. John Wiley & Sons& IEEE Press (2006). [0195] [3] J. B. Pendry, L Martin-Moreno, F. J. Garcia-Vidal, “Mimicking Surface Plasmons with Structured Surfaces”, Science. 305 (5685): 847-848 (2004). [0196] [4] W. X. Tang, H. C. Zhang, H. F. Ma, W. X. Jiang, T. J. Cui. Concept, Theory, Design, and Applications of Spoof Surface Plasmon Polaritons at Microwave Frequencies, Adv. Optical Mater. 1800421, 2018. [0197] [5] S. Rupprecht, C. T. Sica, W. Chen, M. T. Uanagan, and Q. X. Yang. Improvements of transmit efficiency and receive sensitivity with ultrahigh dielectric constant (uHOC) ceramics at 1.5 T and 3 T, Magn. Reson. Med. 79, 2842(2018). [0198] [6] M. J. Freire, L. Jelinek, R. Marques, M. Lapine, On the applications of μ.sub.t=−1 metamaterial lenses for magnetic resonance imaging. J. Magn. Reson. 203, 81(2010). [0199] [7] A. Devilez, X. Zambrana-Puyalto, B. Stout, and N. Bonod, Mimicking localized surface plasmons with dielectric particles, Phys. Rev. B 92, 241412(R)(2015). [0200] [8] C. Rizza, M. Fantasia, E. Palange, M. Alecci, and A. Galante, Harnessing Surface Plasmons for Magnetic Resonance imaging Applications, Phys. Rev. Appl. 12, 044023 (2019).

[0201] Hereto, we have described the preferred embodiments and suggested some variants of the present invention, but it is understood that a person skilled in the art can make modifications and changes without departing from the respective scope of protection, as defined by the appended claims.