Abstract:
Methods of designing an active shield for substantially zeroing an electro-magnetic field on one side of a predetermined boundary in open magnetic resonance imaging systems and in open electric systems are provided. The methods include defining a finite length geometry for a primary structure which, in the case of zeroing a magnetic field, carries a first current distribution on its surface. A finite length geometry is also defined for a secondary structure which, in the case of zeroing a magnetic field, carries a second current distribution on its surface. In the case of zeroing an electric field, the current distributions are replaced with charge distributions. A total magnetic or electric field resulting from a combination of the first and second current or charge distributions respectively is constrained such that normal components (in the magnetic field case) or tangential components (in the electric field case) thereof substantially vanish at the surface of one of the primary and secondary structures. The first current or charge distribution is constrained to the surface of the primary structure, and the second current or charge distribution is constrained to the surface of the secondary structure. The first and second current or charge distributions are then calculated concurrently allowing both the first and second current or charge distributions to vary while observing the constraints such that a predetermined magnetic or electric field is achieved in a first region and a magnetic or electric field on one side of a predetermined boundary is substantially zeroed.

Description:
BACKGROUND OF THE INVENTION 
     The present invention relates to the art of active electro-magnetic shielding. It finds particular application in conjunction with self-shielded magnetic gradient coils used in magnetic resonance imaging (MRI), and will be described with particular reference thereto. However, it is to be appreciated that the present invention is also amenable to other like applications where high quality efficient active shielding of magnetic and/or electric fields is desired. 
     In MRI systems, gradient coil assemblies are commonly pulsed with electrical current pulses to produce magnetic gradients across a main magnetic field in the vicinity of an imaging region. Commonly, the magnetic field gradients interact with external metallic structures such as the magnet cold shields, the magnet dewar, and the like. This interaction tends to generate eddy currents in the affected structures. These eddy currents, in turn, generate eddy magnetic fields which have a deleterious effect on the temporal and spatial quality of the magnetic field in the vicinity of the imaging region and, hence, in the resultant image quality. 
     The eddy current problem is often addressed by placing an active shielding coil between the primary gradient coil and the affected structure. The shielding coils are designed to substantially zero or cancel the residual magnetic field external to itself thereby preventing the formation of eddy currents in potentially vulnerable structures. 
     Previously methods for production for magnetic gradients in MRI systems consisted of winding discrete coils in a bunched or distributed fashion on an electrically insulating hollow cylindrical former and driving the coils with a current source of limited voltage. Conventional bunched coil designs include the Maxwell and the Modified Maxwell Pair for z-gradient production (i.e. axial gradient production), and the Golay or Modified Golay (multi-arc) Saddle Coils for x and/or y-gradient production (i.e. transverse gradient production). Typically, these methods consisted of iteratively placing coil loops or arcs on the cylindrical formers until the desired gradient strength, gradient uniformity, and inductance (related to stored energy) were achieved. These previous designs were generally developed in a “forward approach” whereby a set of initial coil positions were defined (i.e., the initial coil distribution), the fields and the inductance/energy calculated, and if not within particular design parameters, the coil positions would be shifted (statistically or otherwise) and results re-evaluated. This iterative procedure continued until a suitable design was obtained. 
     More recent methods of generating magnetic fields in magnetic resonance imaging systems utilize an “inverse approach.” In the “inverse approach” method, the gradient magnetic field is forced to match predetermined values at specified spatial locations inside the imaging volume and a continuous current density is calculated which is capable of producing such a field. The “inverse approach” method assumes that the primary gradient coil has finite dimensions while those of the secondary or shield coil are left unrestricted (infinite). After the generation of continuous current distributions for both the primary and the shield or secondary coils, an apodization algorithm is performed on the continuous current density of the shield coil in order to restrain it to desirable dimensions. Following the modification of the shielding coil&#39;s continuous current, the Stream Function technique is employed in order to obtain discrete current patterns for both coils. Application of the Biot-Savart law to the discrete current pattern ensures that the discretization procedure was proper. This approach creates generally more energy efficient gradient coil assemblies with higher gradient strengths and faster slew rates as compared to the “forward approach” method. 
     One particular prior art approach is discussed in U.S. Pat. No. 4,794,338 to Roemer, et al. The approach of designing a shielded gradient coil assembly introduced by Roemer, et al. is based on the “forward approach” method. The outcome is a shielded gradient coil assembly with a moderate to low efficiency rate in terms of gradient strength and slew rate. Further, there is no precondition to the method for controlling eddy current effects inside the imaging region. 
     Another particular prior art gradient coil assembly is described in U.S. Pat. No. 5,296,810 to Morich. Morich describes a cylindrically shaped shielded gradient coil assembly for MRI applications. Morich uses the “inverse approach” method of designing gradient coil assemblies where the primary coil has a finite length while the length of the shielding coil is considered infinite. This configuration generates coils with high gradient strengths and slew rates while at the same time reduces the eddy current effects when the length of the shield coil is substantially longer (20% or more) than the length of the primary coil. In order to restrain the current of the shielding coil within desired dimensional boundaries, apodization techniques (e.g., guassian apodization) are employed. In this manner, the overall length of the shielding coil is approximately 20% longer than the total length of the primary gradient coil. Since the shielding coil was modeled as having an infinite length, in order to confine the length of the shielding coil to finite dimensions, current truncation is employed. The current apodization process towards the ends of the shielding coil cause disturbances in the shielding field and ultimately result in unwanted eddy current effects inside the imaging region. These effects become more deleterious as the length of the shielding coil approaches the length of the primary coil, or as both coil dimensions approach the imaging region. 
     In another particular prior art shielded gradient coil assembly described in U.S. Pat. No. 5,132,618 to Sugimoto, the design is based on the “inverse approach” method. In this design, both the primary and the shielded coil lengths were assumed infinite and continuous current densities for both the primary and the secondary coil are modeled based on this assumption. In order to restrain the current densities on both the primary and secondary coils, truncation is again employed. Although the outcome of this method is similar to that of the Morich patent discussed earlier, the additional truncation of the primary coil&#39;s current in this case introduces increased levels of eddy current effects inside the imaging region. 
     For interventional procedures and like applications where patient access is desirable, it is advantageous to design the gradient shielding coil such that its dimensions do not exceed those of the primary gradient coil. In this manner, patient access is maximized and the feeling of openness reduces patient claustrophobia. However, in general, the previous methods and prior art discussed above suffer the drawback that as the shielding coil length approaches that of the primary coil, increased levels of eddy current effects within the imaging region deteriorate image quality. Conversely, when sufficient shielding is achieved, the dimensions of the shielding coil are sufficiently larger than those of the primary coil that the level of patient access is encumbered and an increased level of patient claustrophobia is experienced. 
     Additionally, prior methods commonly determine the current density for the primary coil first, and then determine the appropriate current density for the shielding or secondary coil which will produce the zeroing of the residual magnetic field outside itself. These approaches, in varying degrees, often fail to appreciate the cooperative and interactive nature of the two-coil system. Consequently, they achieve generally less favorable results. 
     The present invention contemplates a new and improved technique for actively shield magnetic and/or electric fields which overcomes the above-referenced problems and others. 
     SUMMARY OF THE INVENTION 
     In accordance with one aspect of the present invention, a method of designing an active shield for substantially zeroing a magnetic field on one side of a predetermined boundary is provided. It includes defining a geometry for a primary coil which carries a first current distribution on its surface. A geometry for a secondary coil is also defined which carries a second current distribution on its surface. The total magnetic field resulting from a combination of the first and second current distributions is constrained such that normal components thereof vanish at the surface of one of the primary and secondary coils. The first current distribution is constrained to the surface of the primary coil of arbitrary length, and the second current distribution is constrained to the surface of the secondary coil of arbitrary length. The first and second current distributions are calculated concurrently while observing the constraints such that a predetermined magnetic field is achieved in a first region and a magnetic field on one side of a predetermined boundary is substantially zeroed. 
     In accordance with a more limited aspect of the present invention, the geometries of the primary and secondary coils define finite length coaxial cylinders with the primary coil being inside the secondary coil. 
     In accordance with a more limited aspect of the present invention, the lengths of the primary coil and the secondary coil are substantially the same. 
     In accordance with a more limited aspect of the present invention, the total magnetic field resulting from the combination of the first and second current distributions is constrained such that normal components thereof vanish at the surface of the secondary coil. The substantially zeroed magnetic field is outside the predetermined boundary which is defined by a cylinder of infinite length that is coincident with the secondary coil. 
     In accordance with a more limited aspect of the present invention, the total magnetic field resulting from the combination of the first and second current distributions is constrained such that normal components thereof vanish at the surface of the primary coil. The substantially zeroed magnetic field is inside the predetermined boundary which is defined by a cylinder of infinite length that is coincident with the primary coil. 
     In accordance with a more limited aspect of the present invention, an amount of energy associated with the first and second current distributions is minimized. 
     In accordance with another aspect of the present invention, a gradient coil assembly is included in a magnetic resonance scanner. The magnetic resonance scanner includes a main magnet for generating a main magnetic field through an examination region, said main magnet being arranged such that its geometry defines said examination region; a couch which supports a subject to be examined within said examination region; a radio frequency coil disposed adjacent said examination region for transmitting radio frequency signals into said examination region and selectively exciting dipoles disposed therein; a radio frequency transmitter for driving said radio frequency coil; a receiver which receives magnetic resonance signals from resonating dipoles within said examination region; and, an image processor which reconstructs an image representation from received magnetic resonance signals for display on a human-readable display. The gradient coil assembly generates a substantially linear predetermined magnetic gradients across said main magnetic field in said examination region. The gradient coil assembly has at least one primary coil disposed about said examination region including a first array of current carrying coil loops arranged on a first former to define a first current density. The gradient coil assembly also has at least one secondary coil disposed about said primary coil between said primary coil and said main magnet including a second array of current carrying coil loops arranged on a second former to define a second current density. The first current density is a discretized version of a first continuous current density, and said second current density is a discretized version of a second continuous current density. A residual magnetic field outside a boundary coincident with said second former is substantially zeroed by concurrently calculating said first and second continuous current densities while: (i) constraining a magnetic field generated by a combination of said first and second continuous current densities such that normal components thereof vanish at a surface defined by said second former, (ii) constraining said first continuous current density to a surface defined by said first former, and (iii) constraining said second continuous current density to said surface defined by said second former. 
     In accordance with a more limited aspect of the present invention, the first and second formers are finite length coaxial cylinders. 
     In accordance with a more limited aspect of the present invention, first and second formers have substantially the same lengths. 
     In accordance with a more limited aspect of the present invention, the boundary coincident with said second former defines an infinitely long cylinder. 
     In accordance with a more limited aspect of the present invention, an amount of energy used by said primary and secondary coils is minimized. 
     In accordance with a more limited aspect of the present invention, the gradient coil assembly includes three pairs of said primary and secondary coils arranged such that each pair generates substantially linear predetermined magnetic gradients along one of three mutually orthogonal axes. Said axes arranged such that one of the axes is coincident with a central axis of the primary coil. 
     In accordance with another aspect of the present invention, a method of designing an active shield for substantially zeroing an electric field on one side of a predetermined boundary is provided. It includes defining a geometry for a primary structure which carries a first charge distribution on its surface. A geometry for a second structure is also defined which carries a second charge distribution on its surface. The total electric field resulting from a combination of the first and second charge distributions is constrained such that tangential components thereof vanish at the surface of one of the primary and secondary structures. The first charge distribution is constrained to the surface of the primary structure, and the second charge distribution is constrained to the surface of the secondary structure. The first and second charge wed distributions are calculated concurrently while observing the constraints such that a predetermined electric field is achieved in a first region and an electric field on one side of a predetermined boundary is substantially zeroed. 
     In accordance with a more limited aspect of the present invention, the geometries of the primary and secondary structures define finite length coaxial cylinders with the primary strucuture being inside the secondary structure. 
     In accordance with a more limited aspect of the present invention, the lengths of the primary structure and the secondary structure are substantially the same. 
     In accordance with a more limited aspect of the present invention, the total electric field resulting from the combination of the first and second charge distributions is constrained such that tangential components thereof vanish at the surface of the secondary structure. The substantially zeroed electric field is outside the predetermined boundary which is defined by a cylinder of infinite length that is coincident with the secondary structure. 
     In accordance with a more limited aspect of the present invention, the total electric field resulting from the combination of the first and second charge distributions is constrained such that tangential components thereof vanish at the surface of the primary structure. The substantially zeroed electric field is inside the predetermined boundary which is defined by a cylinder of infinite length that is coincident with the primary structure. 
     One advantage of the present invention is that high quality electro-magnetic shielding is achieved. 
     Another advantage of the present invention is that it provides an analytical solution to a wide range of shielding applications. 
     Yet another advantage of the present invention is the improved design of a gradient coil assembly for use in MRI scanners. 
     Another advantage of the present invention is that high quality magnetic gradients are achieved in the examination region of an MRI scanner while reducing eddy current effects caused by residual magnetic fields outside the gradient coil assembly. 
     Still further advantages and benefits of the present invention will become apparent to those of ordinary skill in the art upon reading and understanding the following detailed description of the preferred embodiments. 
    
    
     BRIEF DESCRIPTION OF THE DRAWINGS 
     The invention may take form in various components and arrangements of components, and in various steps and arrangements of steps. The drawings are only for purposes of illustrating preferred embodiments and are not to be construed as limiting the invention. 
     FIG. 1A is a diagrammatic illustration of a magnetic resonance imaging apparatus including a shielded gradient coil assembly designed in accordance with the present invention; 
     FIGS. 1B and 1C are diagrammatic illustrations of a shielded gradient coil assembly designed in accordance with the present invention; 
     FIG. 2 is a graph showing J φ  components of a current density on primary and secondary coils in accordance with aspects of the present invention; 
     FIG. 3 is a graph showing in the z direction the behavior of the B z  component of the magnetic field from the primary coil in comparison with the total magnetic field at a radial distance ρ=0.0 m in accordance with aspects of the present invention; 
     FIG. 4A is a graph showing the magnetic field from the primary coil as a function of z at ρ=0.48 m in the region 0.0 m&lt;z&lt;0.55 m in accordance with aspects of the present invention; 
     FIG. 4B is a graph showing the total magnetic field as a function of z at ρ=0.48 m in the region 0.0 &lt;z&lt;0.55 m in accordance with aspects of the present invention; 
     FIG. 5A is a graph showing the magnetic field from the primary coil as a function of z at ρ=0.48 m in the region 0.55 m&lt;z &lt;1.0 m in accordance with aspects of the present invention; 
     FIG. 5B is a graph showing the total magnetic field as a function of z at ρ=0.48 m in the region 0.55 m&lt;z&lt;1.0 m in accordance with aspects of the present invention; 
     FIG. 6A is a graph showing the magnetic field from the primary coil as a function of ρ at z=0.1 m in accordance with aspects of the present invention; 
     FIG. 6B is a graph showing the total magnetic field as a function of ρ at z=0.1 m in accordance with aspects of the present invention; 
     FIG. 7 is a graph showing J φ  components of a current density on primary and secondary coils in accordance with aspects of the present invention when the length of the shield coil is longer than the length of the primary coil; 
     FIGS. 8A and 8B are graphs showing the z dependence of the J φ -components of the current density on the primary and secondary coils, respectively, in accordance with aspects of the present invention; 
     FIGS. 9A and 9B are graphs showing the z dependence of the magnetic field outside the coil system at a radial distance ρ=0.55 m within the length of the coil system (0.0&lt;z&lt;0.4 m) in accordance with aspects of the present invention; 
     FIG. 10 is a graph showing the behavior of the B z  components of the total magnetic field as a function of z inside the primary coil at ρ=0.0 m in accordance with aspects of the present invention; 
     FIGS. 11A and 11B are graphs illustrating the linearity and uniformity, respectively, of the B z  component of the total magnetic field in accordance with aspects of the present invention; 
     FIG. 12 is a graph showing the z dependence of a charge density on both primary and secondary surface structures in accordance with aspects of the present invention; 
     FIG. 13A is a graph showing the z dependence of two non-zero components of the total electric field inside a primary surface structure in accordance with aspects of the present invention; 
     FIG. 13B is a graph showing the z dependence of two components of the electric field between a primary and secondary surface structure at a radial distance ρ=0.4 m in accordance with aspects of the present invention; 
     FIGS. 14A and 14B are graphs showing the contribution of the primary structure to the electric field and the total electric field at ρ=0.55 m, respectively, in accordance with aspects of the present invention; 
     FIGS. 15A through 15C are graphs illustrating the radial ρ dependence of the electric field within the surface structure system in accordance with aspects of the present invention; 
     FIG. 16 is a graph showing the ρ dependence of both components of the electric field at the position z=0.6 m which is beyond the length of the surface structure system in accordance with aspects of the present invention; and, 
     FIG. 17 is a graph showing a continuous charge density for a surface structure system wherein the secondary surface structure is substantially twice as long as the primary surface structure in accordance with aspects of the present invention. 
    
    
     DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
     With reference to FIG. 1A, a main magnetic field control  10  controls superconducting or resistive magnets  12  such that a substantially uniform, temporally constant main magnetic field is created along a z axis through an examination region  14 . A couch (not illustrated) suspends a subject to be examined within the examination region  14 . A magnetic resonance echo means applies a series of radio frequencies (RF) and magnetic field gradient pulses to invert or excite magnetic spins, induce magnetic resonance, refocus magnetic resonance, manipulate magnetic resonance, spatially and otherwise encode the magnetic resonance, to saturate spins, and the like to generate magnetic resonance imaging and spectrography sequences. More specifically, gradient pulse amplifiers  20  apply current pulses to selected ones or pairs of gradient coil assemblies  22  to create magnetic field gradients along mutually orthogonal x, y, and z axes of the examination region  14 . A digital radio frequency transmitter  24  transmits radio frequency pulses or pulse packets to a whole-body RF coil  26  to transmit RF pulses into the examination region. A typical radio frequency pulse is composed of a packet of immediately contiguous pulse segments of short duration which taken together with each other and any applied gradients achieve a selected magnetic resonance manipulation. For whole-body applications, the resonance signals are commonly picked up by the whole-body RF coil  26 . 
     For generating images of local regions of the subject, specialized radio frequency coils are placed continuous to the selected region. For example, an insertable RF coil may be inserted surrounding a selected region at the isocenter of the bore. The insertable RF coil is used to excite magnetic resonance and receive magnetic resonance signals emitting from the patient in the region being examined. Alternatively, the insertable RF coil can be used to only receive resonance signals introduced by whole-body coil RF transmissions. In any event, the resultant radio frequency signals are picked up by the whole-body RF coil  26 , the insertable RF coil, or other specialized RF coils and demodulated by a receiver  30 , preferably a digital receiver. 
     A sequence control circuit  40  controls the gradient pulse amplifiers  20  and the transmitter  24  to generate any of a plurality of multiple echo sequences such as echo planar imaging, echo volume imaging, gradient and spin echo imaging, fast spin echo imaging, and the like. For the selected sequence, the receiver  30  receives a plurality of data lines in rapid succession following each RF excitation pulse. Ultimately, the radio frequency signals received are demodulated and reconstructed into an image representation by a reconstruction processor  50  which applies a two dimensional Fourier transform or other appropriate reconstruction algorithm. The image may represent a planar slice through the patient, an array of parallel planar slices, a three dimensional volume, or the like. The image is then stored in an image memory  52  where it may be accessed by a display, such as a video monitor  54  which provides a human readable display of the resultant image. 
     With reference to FIGS. 1B and 1C, the geometrical configuration of the actively shielded gradient coil assembly  22  is shown. In a preferred embodiment, the gradient coil assembly  22  includes three coaxial primary coil  60  and secondary coil  62  pairs (in the interest of simplification only one pair is shown) arranged such that each pair generates substantially linear predetermined magnetic gradients along one of the three mutually orthogonal axes, x, y, and z. The z axis preferably coincides with the central axis of the coils. For this configuration, the total length of the primary coil  60  is denoted as 2L ρ , and the length of the secondary coil  62  is denoted as 2L s . The radius of the primary coil  60  is denoted as R ρ  while the radius of the secondary or shielding coil  62  is denoted as R s . 
     Heretofore, the conventional design for shielded gradient coil assemblies has generally been based upon the assumption that the total length of the shielding coil is unrestricted and thus its current density is assumed as infinitely extended. Consequently, as “infinite” shielding coils are not available for practical applications, the current density was subsequently apodized or truncated to fit on appropriate finite length elements. 
     On the other hand, the technique employed herein involves the design of a shielded gradient coil configuration whereby the current for the shielding coil is initially restricted within the boundaries of a shielding coil having a finite length. Moreover, the technique employed herein generates the current densities on both the primary and secondary coils concurrently to thereby account for the cooperative and interactive nature of the two coils. 
     In the present invention, for open systems such as two coaxial cylinders of arbitrary radii (R ρ  for the inner cylinder or primary coil and R s  for the outer cylinder or secondary coil) and arbitrary lengths (2L ρ  for the primary coil and 2L s  for the secondary coil), current distributions (in the case of magnetic field shielding) or alternately charge distributions (in the case of electric field shielding) on both cylinders are simultaneously found such that a quasi-static magnetic field or quasi-static electric field vanishes in the entire region outside an infinitely long cylindrical space having a radius coincident with that of the secondary coil (defined and referred to hereinafter as an exterior super shielding effect). Similarly, current or charge distributions (different from those for exterior super shielding) are also found such that a quasi-static magnetic field or quasi-static electric field vanishes in the entire region inside an infinitely long cylindrical space having a radius coincident with that of the primary coil (defined and referred to hereinafter as an interior super shielding effect). In both circumstances, for the exterior and/or interior super shielding, there is sufficient freedom such that the desired field behavior inside the primary coil, or outside the secondary coil, respectively, is achieved. 
     The following details the method of the present invention for determining the current or charge distributions of primary and secondary coils to achieve the aforementioned super shielding effects. More specifically, in an exemplary gradient coil assembly  22  having the geometric configuration illustrated in FIGS. 1B and 1C, the method is applied to achieve the exterior super shielding effect for a magnetic field. 
     SUPER SHIELDING OF THE MAGNETIC FIELD 
     I. Current Distribution 
     First, consider a current density distribution {overscore (J)}({overscore (r)}) on a cylinder having a radius R and a total length 2L. For the cylindrical components on the surface of the cylinder given by 
     
       
           {overscore (J)} (ρ,φ, z )={ ê   φ   J   φ (φ, z )+ ê   z   J   z (φ, z )}δ(ρ− R )  (1), 
       
     
     where δ(ρ−R) is the restriction that confines the electrical current on a cylinder surface with radius R, the current density continuity equation {overscore (∇)}·{overscore (J)}({overscore (r)})=0 has the form:                    1   R                       ∂     J   φ         ∂   φ         +       ∂     J   z         ∂   z         =   0.           (   2   )                                
     A Fourier series representation follows from the periodicity in φ. Current conservation eliminates the J z -component of the zeroth harmonic, and reduces the current density to the form:                    J   →          (     τ   →     )       =       (           J   →     0          (   z   )       +       ∑     i   =   1     2            ∑     n   =   1     ∞              J   →     n     (   i   )            (     φ   ,   z     )             )          δ        (     ρ   -   R     )           ;           (   3   )                                
     where 
     
       
           {overscore (J)}   0 ( z )= ê   φ   f   0 ( z )  (4), 
       
     
     
       
           {overscore (J)}   η   (1) (φ, z )= ê   φ   f   η ( z )cos(ηφ)+ ê   x   q   η ( z )sin(ηφ)  (5), 
       
     
     
       
           {overscore (J)}   η   (2) (φ, z )= ê   φ   g   η ( z )sin(ηφ)+ ê   z ρ η ( z )cos(ηφ)  (6), 
       
     
     and where f 0 (z) is an arbitrary function and f η (z), q η (z), g η (z), and ρ η (z) satisfy the following relations for η≧1:                      n   R            f   n          (   z   )         -       q   n   ′          (   z   )         =   0     ,           n   R            g   n          (   z   )         +       p   n   ′          (   z   )         =   0     ,           (   7   )                                
     where q η   1 (z) and p η   1 (z) are the first derivatives of q η (z) and ρ η (z) respectively, and f η (z), q η   1 (z), g η (z), and ρ η   1 (z) are adjusted to achieve desired field qualities. The end-condition is that there can be no abrupt change in the J φ -component of the current density, yielding 
     
       
           q   η (| z|≧L )=ρ η (| z|≧L )=0  (8). 
       
     
     However, there can be an abrupt change in the J φ -component of the current density, yielding 
     
       
           f   η (| z|&lt;L )=0  (9). 
       
     
     For an azimuthally directed current density, only the term {overscore (J)} 0 (z) in equation (4) is employed. For other applications, a combination of the terms of equations (4)-(6) is employed. 
     II. Magnetic Field 
     Next, consider the vector potential {overscore (A)}({overscore (r)}) corresponding to the current distribution from equation (3), which takes the form:                  A   →          (     r   →     )       =         μ   o       4                 π            ∫       G        (       r   →     ,       r   →     ′       )              J   →          (         r   →                  ′     )                   r   →     ′                     (   10   )                                
     where G({overscore (r)},{overscore (r)} 1 ) is the Green function of the Maxwell equations which in cylindrical coordinates has the following form:                G        (       r   →     ,       r   →     ′       )       =       2   π            ∑     m   =     -   ∞       ∞                 im        (     φ   -     φ   ′       )                ∫   0   ∞            cos        (     k        (     z   -     z   ′       )       )              T   m          (     k   ,   ρ   ,     ρ   ′       )               k                       (   11   )                                
     where 
     
       
           T   m ( k, ρ,ρ   1 )=θ(ρ−ρ 1 ) I   m ( kρ   1 ) K   m ( kρ )+θ(ρ 1 −ρ) I   m ( kρ ) K   m (ρ 1 )  (12) 
       
     
     and where I m (kρ) and K m (kρ) are the modified Bessel functions of the first and second kind, and θ(x) is the Heaviside step function. The vector potential then has the following representation:                  A   →          (     r   →     )       =           A   →     0          (     r   →     )       +       ∑     i   =   1     2            ∑     n   =   1     ∞              A   →     n     (   i   )            (     r   →     )                     (   13   )                                
     where {overscore (A)} 0 ({overscore (r)}) corresponds to the current density {overscore (J)} 0 (z) and the terms with the indices i=1, 2 correspond to the current densities {overscore (J)} η   (1)  (φ,z) and {overscore (J)} η   (2)  (φ,z) respectively. 
     Starting with the magnetic field corresponding to the term {overscore (J)} 0 (z) in equation (3), there are only two components of the magnetic field in this case, B 0,ρ ({overscore (r)}) and B 0,z ({overscore (r)}), which take the following forms:                    B     0   ,   ρ            (     r   →     )       =       -         μ   o        R     π              ∫   0   ∞              kT     1   ,   0            (     k   ,   ρ   ,   R     )              ∫     -   L     L              f   0          (     z   ′     )            sin        (     k        (     z   -     z   ′       )       )                 z   ′               k                                 and   ,             (   14   )                     B     0   ,   z            (     r   →     )       =       -         μ   o        R     π              ∫   0   ∝              kT     2   ,   0            (     k   ,   ρ   ,   R     )              ∫     -   L     L              f   0          (     z   ′     )            cos        (     k        (     z   -     z   ′       )       )                 z   ′               k                               where           (   15   )                     T     1   ,   n            (     k   ,   ρ   ,   R     )       =         θ        (     ρ   -   R     )              I   n   ′          (   kR   )              K   n   ′          (     k                 ρ     )         +       θ        (     R   -   ρ     )              I   n   ′          (     k                 ρ     )              K   n   ′          (   kR   )                             and   ,             (   16   )                   T     2   ,   n            (     k   ,   ρ   ,   R     )       =         θ        (     ρ   -   R     )              I   n   ′          (   kR   )              K   n          (     k                 ρ     )         +       θ        (     R   -   ρ     )              I   n          (     k                 ρ     )                K   n   ′          (   kR   )       .                 (   17   )                                
     The magnetic field expressions that correspond to the current densities {overscore (J)} η   (1)  (φ,z), for η≧1 are                    B     n   ,   ρ       (   1   )            (     r   →     )       =       -         μ   o          R   2       π                         cos        (     n                 φ     )       n            ∫   0   ∞            k   2            T     1   ,   n            (     k   ,   ρ   ,   R     )              ∫     -   L     L              q   n          (     z   ′     )            cos        (     k        (     z   -     z   ′       )       )                 z   ′               k                 ,           (   18   )                     B     n   ,   φ       (   1   )            (     r   →     )       =       -         μ   o          R   2       πρ                       sin        (   nφ   )              ∫   0   ∞              kT     2   ,   n            (     k   ,   ρ   ,   R     )              ∫     -   L     L              q   n          (     z   ′     )            cos        (     k        (     z   -     z   ′       )       )                 z   ′               k                 ,              and           (   19   )                   B     n   ,   z       (   1   )            (     r   →     )       =           μ   o          R   2       π                       cos        (     n                 φ     )       n            ∫   0   ∞            k   2            T     2   ,   n            (     k   ,   ρ   ,   R     )              ∫     -   L     L              q   n          (     z   ′     )            sin        (     k        (     z   -     z   ′       )       )                 z   ′                 k     .                       (   20   )                                
     The expressions for the magnetic field that correspond to the current density {overscore (J)} η   (2) (φ,z), for η≧1 are found from equations (18)-(20) by the following substitutions: 
     
       
         cos(ηφ)→−sin(ηφ),sin(ρφ)→cos(ηφ), q   η ( z )→ρ θ ( z )  (21). 
       
     
     To shield the coil system, consider two coaxial cylinders of arbitrary radii, R ρ  for the inner cylinder or primary coil and R x  for outer cylinder or secondary coil, and arbitrary lengths of, 2L ρ  and 2L z , respectively. On each cylinder, the current distribution has the form given by equation (3). The total magnetic field then is a vector sum of the magnetic fields that correspond to the primary and secondary coils. In what follows all quantities that correspond to the primary coil are labeled with an index &#39;ρ&#39;, while all quantities that correspond to the secondary coil are labeled with an index &#39;s&#39;. 
     III. Exterior Super Shielding of the Magnetic Field 
     For excellent shielding of the magnetic field outside the secondary coil, a constraint is introduced which simulates the boundary condition expected for a very good conductor (or even a &#39;superconductor&#39;) by making the normal component of the total magnetic field vanish at the surface of the secondary coil, expressed analytically as 
     
       
           B   η,ρ   tot(1) ( {overscore (r)} )| ρ=R     s   =0  (22) 
       
     
     which holds for any η≧0. Sufficient conditions satisfying equation (22) are                    F     0   ,     L   s       ±          (   k   )       +         R   p       R   s                           I   0   ′          (     kR   p     )           I   0   ′          (     kR   s     )                F     0   ,     L   p       ±          (   k   )           =   0           (   23   )                                
     for axial coils (such as z-gradient coils or main magnet coils), and                      Q     n   ,     L   s       ±          (   k   )       +         R   p   2       R   s   2                I   n   ′          (     kR   p     )           I   n   ′          (     kR   s     )                Q     n   ,     L   p       ±          (   k   )           =   0     ,     n   ≥   1             (   24   )                                
     for transverse coils (such as x-gradient coils or y-gradient coils); where,                    F     0   ,   L     ±          (   k   )       =       ∫     -   L     L              f   0          (   z   )            (           cos        (   kz   )                 sin        (   kz   )             )             z                         and           (   25   )                     Q     n   ,   L     ±          (   k   )       =       ∫     -   L     L              q   n          (   z   )            (           cos        (   kz   )                 sin        (   kz   )             )             z           ,     n   ≥   1             (   26   )                                
     and the sign &#39;+&#39; refers to an even function of z and the sign &#39;−&#39; refers to an odd function of z. The full implications of these conditions determine those current densities on both coils that satisfy the boundary condition of equation (22). 
     An inverse Fourier transform of equations (23) and (24) together with the definitions from equations (25) and (26) lead to two types of additional consistency conditions                    1   π            ∫   0   ∞              F     0   ,     L   p       ±          (   k   )            (           cos        (   kz   )                 sin        (   kz   )             )             k           =       θ        (       L   p     -        z          )              f   0     ±     ,     (   p   )                (   z   )                         and           (   27   )                   1   π            ∫   0   ∞              F     0   ,     L   s       ±          (   k   )            (           cos        (   kz   )                 sin        (   kz   )             )             k           =       θ        (       L   s     -        z          )              f   0     ±     ,     (   s   )                (   z   )                 (   28   )                                
     for axial coils, and                    1   π            ∫   0   ∞              Q     n   ,     L   p       ±          (   k   )            (           cos        (   kz   )                 sin        (   kz   )             )             k           =       θ        (       L   p     -        z          )              q   n     ±     ,     (   p   )                (   z   )           ,       n   ≥   1                   and             (   29   )                     1   π            ∫   0   ∞              Q     n   ,     L   s       ±          (   k   )            (           cos        (   kz   )                 sin        (   kz   )             )             k           =       θ        (       L   s     -        z          )              q   n     ±     ,     (   s   )                (   z   )           ,     n   ≥   1             (   30   )                                
     for transverse coils, where 
     
       
           f   0   ±,(ρ,s) ( z )=½( f   0   (ρ,s) ( z )± f   0   ρ,s) (− z )) q   η   ±,(ρ,s) ( z )=½( q   η   (ρ,s) ( z )±i q η   (ρ,s) (− z )  (31). 
       
     
     When |z|&gt;L ρ  the equations (27) and (29) represent constraints for the current density on the primary coil. When |z|≧L s  the equations (28) and (30) represent constraints for the current density on the secondary coil. These consistency conditions represent the fact that the current does not flow off the ends of the coils. That is to say, the current densities are confined to finite cylinders. 
     For exterior shielding of the axial coils, the sufficiency condition (23) plus the consistency conditions (27)-(28) constitute the supershielding set of equations. By an inverse Fourier transform of (23), subject to (27)-(28), and an application of the Fourier convolution theorem, the supershielding set of equations is found to be equivalent to                    θ        (       L   s     -        z          )              f   0     ±     ,     (   s   )                (   z   )         +       1   π            ∫     -     L   p         L   p                r   a          (     z   -   y     )              f   0     ±     ,     (   p   )                (   y   )                          y             =   0           (23a)             where                             r   a          (     z   -   y     )       =       ∫   0   ∞              R   a          (   k   )          cos                   k        (     z   -   y     )                          k                 (     23                 b     )                                
     which is the Fourier transform of the coefficient in (23)                  R   a          (   k   )       =         R   p       R   s                  I   0   ′          (     kR   p     )           I   0   ′          (     kR   s     )         .               (23c)                                
     For exterior shielding of the transverse coils, the sufficiency condition (24) plus the consistency conditions (29)-(30) constitute the supershielding set of equations. By an inverse Fourier transform of (24), subject to (29)-(30), and an application of the Fourier convolution theorem, the supershielding set of equations is found to be equivalent to                    θ        (       L   s     -        z          )              q   n     ±     ,     (   s   )                (   z   )         +       1   π            ∫     -     L   p         L   p                r   t          (     z   -   y     )              q   n     ±     ,     (   p   )                (   y   )                          y             =   0           (24a)             where                             r   t          (     z   -   y     )       =       ∫   0   ∞              R   t          (   k   )          cos                   k        (     z   -   y     )                          k                 (     24      b     )                                
     which is the Fourier transform of the coefficient in (24)                  R   t          (   k   )       =         R   p   2       R   s   2                  I   n   ′          (     kR   p     )           I   n   ′          (     kR   s     )         .               (24c)                                
     The shielding conditions from equations (23) and (24) are now used to write the expressions for all components of the total magnetic field. Inside an infinite cylinder with a radius given by the primary coil radius (i.e. ρ&lt;R ρ ) this is expressed as follows:                    B     0   ,   ρ     tot          (     r   →     )       =       -         μ   o          R   p       π              ∫   0   ∞              kI   0   ′          (     k                 ρ     )              G   0          (     k   ,     R   p     ,     R   s       )              ∫     -     L   p         L   p                f   0     (   p   )            (     z   ′     )            sin        (     k        (     z   -     z   ′       )       )                 z   ′               k                 ,           (   32   )                     B     0   ,   z     tot          (     r   →     )       =       -         μ   o          R   p       π              ∫   0   ∞              kI   0          (     k                 ρ     )              G   0          (     k   ,     R   p     ,     R   s       )              ∫     -     L   p         L   p                f   0     (   p   )            (     z   ′     )            cos        (     k        (     z   -     z   ′       )       )                 z   ′               k                 ,           (   33   )                     B     n   ,   ρ       tot        (   1   )              (     r   →     )       =       -         μ   o          R   p   2       π                         cos        (     n                 φ     )       n     ×       ∫   0   ∞            k   2                       I   n   ′     (                k                 ρ     )                       G   n     (                k   ,                R   p     ,                R   s       )                       ∫     -     L   p         L   p                q   n     (   p   )            (     z   ′     )          cos                   (     k        (     z   -     z   ′       )       )               z   ′               k                 ,                n   ≥              1     ,           (   34   )                     B     n   ,   φ       tot        (   1   )              (     r   →     )       =           μ   o          R   p   2       πρ                     sin        (     n                 φ     )       ×       ∫   0   ∞          k                     I   n     (                k                 ρ     )                       G   n     (                k   ,                R   p     ,                R   s       )                       ∫     -     L   p         L   p                q   n     (   p   )            (     z   ′     )          cos                   (     k        (     z   -     z   ′       )       )               z   ′               k                 ,                n   ≥              1     ,           (   35   )             and                               B     n   ,   z       tot        (   1   )              (     r   →     )       =           μ   o          R   p   2       π                       cos        (     n                 φ     )       n     ×       ∫   0   ∞          k                     I   n     (                k                 ρ     )                       G   n     (                k   ,                R   p     ,                R   s       )                       ∫     -     L   p         L   p                q   n     (   p   )            (     z   ′     )          sin                   (     k        (     z   -     z   ′       )       )               z   ′               k                 ,                n   ≥              1             (   36   )             where                             G   n          (     k   ,     R   p     ,     R   s       )       =         K   n   ′          (     kR   p     )       -           I   n   ′          (     kR   p     )           I   n   ′          (     kR   s     )                  K   n   ′          (     kR   s     )       .                 (   37   )                                
     In the region between two cylinders of infinite length whose radii are the same as those of the primary and secondary coils (i.e. R ρ &lt;ρ&lt;R s ), the expressions for the total magnetic field are as follows:                    B     0   ,   ρ     tot          (     r   →     )       =       -         μ   o          R   p       π                         ∫   0   ∞          k            I   0   ′          (     kR   p     )           I   0   ′          (     kR   s     )         ×     {           I   0   ′          (     k                 ρ     )              K   0   ′          (     kR   s     )         -         I   0   ′          (     kR   s     )              K   0   ′          (     k                 ρ     )           }                       ∫     -     L   p         L   p                f   0     (   p   )            (     z   ′     )                       sin        (     k        (     z   -     z   ′       )       )                 z   ′               k                 ,           (   38   )                     B     0   ,   z     tot          (     r   →     )       =           μ   o          R   p       π                       ∫   0   ∞          k            I   0   ′          (     kR   p     )           I   0   ′          (     kR   s     )         ×     {           I   0          (     k                 ρ     )              K   0   ′          (     kR   s     )         -         I   0   ′          (     kR   s     )              K   0          (     k                 ρ     )           }                       ∫     -     L   p         L   p                f   0     (   p   )            (     z   ′     )                       cos        (     k        (     z   -     z   ′       )       )                 z   ′               k                 ,           (   39   )                     B     n   ,   ρ       tot        (   1   )              (     r   →     )       =       -         μ   o          R   p   2       π              cos        (     n                 φ     )       n            ∫   0   ∞            k   2              I   n   ′          (     kR   p     )           I   n   ′          (     kR   s     )         ×     {         -       I   n   ′          (     k                 ρ     )                K   n   ′          (     kR   s     )         +         I   n   ′          (     kR   s     )              K   n   ′          (     k                 ρ     )           }            ∫     -     L   p         L   p                q   n     (   p   )            (     z   ′     )            cos        (     k        (     z   -     z   ′       )       )                 z   ′               k                 ,     n   ≥   1     ,           (   40   )                     B     n   ,   φ       tot        (   1   )              (     r   →     )       =           μ   o          R   p   2       πρ          sin        (     n                 φ     )              ∫   0   ∞          k            I   n   ′          (     kR   p     )           I   n   ′          (     kR   s     )         ×     {         -       I   n          (     k                 ρ     )                K   n   ′          (     kR   s     )         +         I   n   ′          (     kR   s     )              K   n          (     k                 ρ     )           }                       ∫     -     L   p         L   p                q   n     (   p   )            (     z   ′     )                       cos        (     k        (     z   -     z   ′       )       )                 z   ′               k                 ,     n   ≥   1     ,           (   41   )             and                               B     n   ,   z       tot        (   1   )              (     r   →     )       =           μ   o          R   p   2       π            cos        (     n                 φ     )       n            ∫   0   ∞          k            I   n   ′          (     kR   p     )           I   n   ′          (     kR   s     )         ×     {         -       I   n          (     k                 ρ     )                K   n   ′          (     kR   s     )         +         I   n   ′          (     kR   s     )              K   n          (     k                 ρ     )           }                       ∫     -     L   p         L   p                q   n     (   p   )            (     z   ′     )                       sin        (     k        (     z   -     z   ′       )       )                 z   ′               k                 ,     n   ≥   1.             (   42   )                                
     In the region outside an infinite cylinder whose radius is the same as that of the secondary coil (i.e. ρ&gt;R s ), the total magnetic field vanishes, which is expressed as follows: 
     
       
           B   η,ρ   tot(1) ( {overscore (r)} )= B   η,φ   tot(1) ( {overscore (r)} )= B   η,z   tot(1) ( {overscore (r)} )=0,η≧0  (43). 
       
     
     Equation (43) represents the principal result of this section. The shielding conditions given by equations (23)-(24), together with the consistency conditions from equations (27)-(30), lead to the complete cancellation of the magnetic field in the space outside the infinite cylinder coincident with the secondary coil. Moreover, the field components inside this infinite cylinder all vanish reasonably rapidly in the z-direction as the distance increases from the ends of the coils. 
     IV. Interior super shielding of the Magnetic Field 
     For excellent shielding of the magnetic field inside the primary coil, a constraint is introduced which simulates the boundary condition expected for a very good conductor (or even a &#39;superconductor&#39;) by making the normal component of the total magnetic field vanish at the surface of the primary coil, expressed analytically as 
     
       
           B   η,ρ   tot(1) ( {overscore (r)} )| η=R     ρ   =0  (44) 
       
     
     which holds for any mode η≧0. Sufficient conditions satisfying equation (44) are                      F     0   ,     L   p       ±          (   k   )       +         R   s       R   p                K   0   ′          (     kR   s     )           K   0   ′          (     kR   p     )                F     0   ,     L   s       ±          (   k   )           =   0                     and             (   45   )                       Q     n   ,     L   p       ±          (   k   )       +         R   s   2       R   p   2                K   n   ′          (     kR   s     )           K   n   ′          (     kR   p     )                Q     n   ,     L   s       ±          (   k   )           =   0     ,                n              ≥   1.             (   46   )                                
     Here the same notations have been used for the Fourier transforms for the current densities as in the previous section. The consistency conditions for the system given by equations (45) and (46) have the same form as in equations (27)-(30). 
     For interior shielding of the axial coils, the sufficiency condition (45) plus the consistency conditions having the same form as (27)-(28) constitute the supershielding set of equations. By an inverse Fourier transform of (45), subject to the consistency conditions, and an application of the Fourier convolution theorem, the supershielding set of equations is found to be equivalent to                    θ        (       L   p     -        z          )              f   0     ±     ,     (   p   )                (   z   )         +       1   π            ∫     -     L   p         L   p                r   A          (     z   -   y     )              f   0     ±     ,     (   s   )                (   y   )                          y             =   0           (45a)             where                             r   A          (     z   -   y     )       =       ∫   0   ∞              R   A          (   k   )          cos                   k        (     z   -   y     )                          k                 (     45      b     )                                
     which is the Fourier transform of the coefficient in (45)                  R   A          (   k   )       =         R   s       R   p                  K   0   ′          (     kR   s     )           K   0   ′          (     kR   p     )         .               (45c)                                
     For interior shielding of the transverse coils, the sufficiency condition (46) plus the consistency conditions having the same form as (27)-(28) constitute the supershielding set of equations. By an inverse Fourier transform of (46), subject to the consistency conditions, and an application of the Fourier convolution theorem, the supershielding set of equations is found to be equivalent to                      Θ        (       L   p     -        z          )              q   n     ±     ,     (   p   )                (   z   )         +       1   π            ∫     -     L   s         L   s                r   T          (     z   -   y     )              q   n     ±     ,     (   s   )                (   y   )               y             =   0                   where           (46a)                   r   T          (     z   -   y     )       =       ∫   0   ∞              R   T          (   k   )          cos                   k        (     z   -   y     )               k                 (46b)                                
     which is the Fourier transform of the coefficient in (46)                  R   T          (   k   )       =         R   s   2       R   p   2                  K   n   ′          (     kR   s     )           K   n   ′          (     kR   p     )         .               (46c)                                
     The interior shielding conditions are now used to write the expressions for all components of the total magnetic field. Inside an infinite cylinder with a radius given by the primary coil radius (i.e. ρ&lt;R ρ ), the total magnetic field exactly vanishes, which is expressed analytically as 
     
       
           B   η,ρ   tot(1) ({overscore (r)})= B   η,φ   tot(1) ( {overscore (r)} )= B   η,z   tot(1) ( {overscore (r)} )=0,η≧0  (47). 
       
     
     In the region between two cylinders of infinite length whose radii are the same as those of the primary and secondary coils (i.e. R ρ &lt;ρ&lt;R s ), the expressions for all components of the total magnetic field are                    B     0   ,   p     tot          (     r   ⇀     )       =         -         μ   0          R   s       π              ∫   0   ∞          k            K   0   ′          (     kR   s     )           K   0   ′          (     kR   p     )         ×     {            I   0   ′          (     k                 ρ     )                         K   0   ′     (                kR   p     )           -         I   0   ′          (     kR   p     )              K   0   ′          (     k                 ρ     )            }                       ∫     -     L   s         L   s                f   0     (   s   )            (     z   ′     )                       sin        (     k        (     z   -     z   ′       )       )                 z   ′               k               ,           (   48   )                     B     0   ,   z     tot          (     r   ⇀     )       =         -         μ   0          R   s       π              ∫   0   ∞          k            K   0   ′          (     kR   s     )           K   0   ′          (     kR   p     )         ×     {            I   0   ′          (     k                 ρ     )                         K   0   ′     (                kR   p     )           -         I   0   ′          (     kR   p     )              K   0          (     k                 ρ     )            }                       ∫     -     L   s         L   s                f   0     (   s   )            (     z   ′     )                       cos        (     k        (     z   -     z   ′       )       )                 z   ′               k               ,           (   49   )                     B     n   ,   p       tot        (   1   )              (     r   ⇀     )       =             μ   0          R   s   2       π            cos        (     n                 φ     )       n            ∫   0   ∞            k   2              K   n   ′          (     kR   s     )           K   n   ′          (     kR   p     )         ×     {           -                    I   n   ′          (     k                 ρ     )                           K   n   ′     (                  kR     p   )       +         I   n   ′          (     kR   p     )              K   n   ′          (     k                 ρ     )             }                       ∫     -     L   s         L   s                q   n     (   s   )            (     z   ′     )                       cos        (     k        (     z   -     z   ′       )       )                 z   ′               k               ,                n   ≥   1     ,           (   50   )                     B     n   ,   φ       tot        (   1   )              (     r   ⇀     )       =             μ   0          R   s   2       πρ          sin        (     n                 φ     )              ∫   0   ∞          k            K   n   ′          (     kR   s     )           K   n   ′          (     kR   p     )         ×     {           -                    I   n          (     k                 ρ     )                         K   n   ′     (                kR   p     )       +         I   n   ′          (     kR   p     )              K   n          (     k                 ρ     )            }                       ∫     -     L   s         L   s                q   n     (   s   )            (     z   ′     )                       cos        (     k        (     z   -     z   ′       )       )                 z   ′               k               ,                n   ≥   1     ,                and                                                                                                                    (   51   )                     B     n   ,   z       tot        (   1   )              (     r   ⇀     )       =             μ   0          R   s   2       π            cos        (     n                 φ     )       n            ∫   0   ∞          k            K   n   ′          (     kR   s     )           K   n   ′          (     kR   p     )         ×     {           -                    I   n          (     k                 ρ     )                           K   n   ′     (                  kR     p   )       +         I   n   ′          (     kR   p     )              K   n          (     k                 ρ     )             }                       ∫     -     L   s         L   s                q   n     (   s   )            (     z   ′     )                       sin        (     k        (     z   -     z   ′       )       )                 z   ′               k               ,                n   ≥   1.             (   52   )                                
     In the region outside an infinite cylinder whose radius is the same as that of the secondary coil (i.e. ρ&gt;R s ), all components of the total magnetic field have the following form:                    B     0   ,   ρ     tot          (     r   ⇀     )       =       -         μ   0          R   s       π              ∫   0   ∞              kK   0   ′          (     k                 ρ     )              D   0          (     k   ,                R   p     ,                R   s       )              ∫     -     L   s         L   s                f   0     (   s   )            (     z   ′     )            sin        (     k        (     z   -     z   ′       )       )                 z   ′               k                 ,           (   53   )                     B     0   ,   z     tot          (     r   ⇀     )       =       -         μ   0          R   s       π              ∫   0   ∞              kK   0          (     k                 ρ     )              D   0          (     k   ,                R   p     ,                R   s       )              ∫     -     L   s         L   s                f   0     (   s   )            (     z   ′     )            cos        (     k        (     z   -     z   ′       )       )                 z   ′               k                 ,           (   54   )                     B     0   ,   ρ       tot        (   1   )              (     r   ⇀     )       =       -         μ   0          R   s   2       π              cos        (     n                 φ     )       n     ×       ∫   0   ∞            k   2                       K   n   ′          (     k                 ρ     )                         D   n          (     k   ,                R   p     ,                R   s       )                         ∫     -     L   s         L   s                q   n     (   s   )            (     z   ′     )            cos        (     k        (     z   -     z   ′       )       )                 z   ′               k                 ,     n   ≥              1     ,           (   55   )                     B     n   ,   φ       tot        (   1   )              (     r   ⇀     )       =           μ   0          R   p   2       πρ          sin        (     n                 φ     )       ×       ∫   0   ∞          k                     K   n          (     k                 ρ     )                         D   n          (     k   ,                R   p     ,                R   s       )                         ∫     -     L   s         L   s                q   n     (   s   )            (     z   ′     )            cos        (     k        (     z   -     z   ′       )       )                 z   ′               k                 ,     n   ≥              1     ,                  and                                                                                                                  (   56   )                     B     0   ,   ρ       tot        (   1   )              (     r   ⇀     )       =       -         μ   0          R   s   2       π              cos        (     n                 φ     )       n     ×       ∫   0   ∞            k   2                       K   n   ′          (     k                 ρ     )                         D   n          (     k   ,                R   p     ,                R   s       )                         ∫     -     L   s         L   s                q   n     (   s   )            (     z   ′     )            cos        (     k        (     z   -     z   ′       )       )                 z   ′               k                 ,     n   ≥              1     ,           (   57   )               where                                          D   n          (     k   ,                R   p     ,                R   s       )       =         I   n   ′          (     kR   s     )       -           K   n   ′          (     kR   s     )           K   n   ′          (     kR   s     )                  I   n   ′          (     kR   p     )       .                 (   58   )                                
     Equation (47) is the principal result of this section. The interior shielding conditions from equations (45) and (46) together with their consistency conditions lead to the complete cancellation of the magnetic field in the space inside the infinite cylinder coincident with the primary coil. Again, the field components inside this infinite cylinder all vanish reasonably rapidly in the z-direction as the distance increases from the ends of the coils. 
     Functional Analysis 
     A definite solution for the current distributions on the primary and secondary coils is obtained by posting an energy-like functional and incorporating the shielding constraints, along with other design constraints, in an extreme calculation. 
     I. Exterior Super Shielding of the Magnetic Field 
     The functional analysis begins with the introduction of an energy-like functional having the form:                W   g     =       ∑     n   =   0     ∞                     (       W     n   ,   g       (   -   )       +     W     n   ,   g       (   +   )         )               (   59   )             where                             W     0   ,   g       (   ±   )       =         μ   0     2            ∫   0   ∞            g        (   k   )              W   0          (     k   ,                R   p     ,                R   s       )              F     0   ,     L   p           (   ±   )        2            (   k   )                          k             ,           (   60   )                     W   0   ′                     (     k   ,                R   p     ,                R   s       )       =       -     R   p   2              I   0   ′          (     kR   p     )              Q   0          (     k   ,                R   p     ,                R   s       )           ,           (   61   )                   W     n   ,   g       (   ±   )       =         μ   0     2            ∫   0   ∞            g        (   k   )              W   n          (     k   ,                R   p     ,                R   s       )              Q     n   ,     L   p           (   ±   )        2            (   k   )                          k             ,           (   62   )               and                                 W   n   ′                     (     k   ,                R   p     ,                R   s       )       =       -         R   p   4          k   2         n   2                I   n   ′          (     kR   p     )              G   n          (     k   ,                R   p     ,                R   s       )           ,     n   ≥   1     ,           (   63   )                                
     and where g(k) is a regularization function for controlling the energy and the current behavior and will be introduced later. This function grows sufficiently fast with k in order to eliminate singular currents at the edges of the coils. For g=1, on the other hand, the functional W g  reduces to a Parseval&#39;s representation of the energy functional when the shielding conditions have been implemented. This takes the following form: 
     
       
           W   g=1 =½∫ {right arrow over (J)} ( {right arrow over (r)} )· {right arrow over (A)} ( {right arrow over (r)} ) d{right arrow over (r)}   (64). 
       
     
     The diagonal form (in the mode index η as well as in the symmetry index ±) of the functional W g  represents the fact that, when the shielding conditions are imposed the mutual inductance between two coils with different mode indices and symmetry indices ±, is zero. Next, in order to minimize the energy, the minimum of the functional W g  subject to the constraints imposed is found through the use of Lagrange multipliers. The constraints used in this variational principle are the shielding conditions from equations (23)-(24), its consistency conditions from equations (27)-(30), and desired values of the magnetic field at some predetermined points. The consistency conditions represented by equations (27)-(30) are replaced with their discretized versions as follows:                    1   π            ∫   0   ∞              F     0   ,     L   p       ±          (   k   )            (           cos        (     k   a     (   p   )       )                 sin        (     k   a     (   p   )       )             )             k           =   0     ,     a   =   1     ,     .           .           .                ,     N   p     ,               a     (   p   )            &gt;     L   p       ,           (   65   )                     1   π            ∫   0   ∞              F     0   ,     L   s       ±          (   k   )            (           cos        (     k   b     (   s   )       )                 sin        (     k   b     (   s   )       )             )             k           =   0     ,     b   =   1     ,     .           .           .                ,     N   s     ,               b     (   s   )            &gt;     L   p       ,           (   66   )                     1   π            ∫   0   ∞              Q     n   ,     L   p       ±          (   k   )            (           cos        (     k   a     (   p   )       )                 sin        (     k   a     (   p   )       )             )             k           =   0     ,     b   =   1     ,     .           .           .                ,     N   s     ,               b     (   p   )            ≥     L   p       ,              and           (   67   )                     1   π            ∫   0   ∞              Q     n   ,     L   s       ±          (   k   )            (           cos        (     k   b     (   s   )       )                 sin        (     k   b     (   s   )       )             )             k           =   0     ,     b   =   1     ,     .           .           .                ,     N   s     ,               b     (   s   )            ≥     L   s               (   68   )                                
     where N ρ , N s  are the number of points where the current densities vanish outside the primary and secondary coils (in the z direction). When |Z a   (ρ) |&lt;L ρ  and |Z b   (s) |&lt;L s  the expressions for the currents density can be obtained from equations (27)-(30) as follows:                    f   0     ±     ,     (   p   )                (   z   )       =       1   π            ∫   0   ∞              F     0   ,     L   p       ±          (   k   )            (           cos        (   kz   )                 sin        (   kz   )             )             k             ,           (   69   )                     f   0     ±     ,     (   s   )                (   z   )       =       1   π            ∫   0   ∞              F     0   ,     L   s       ±          (   k   )            (           cos        (   kz   )                 sin        (   kz   )             )             k             ,           (   70   )                     q   n     ±     ,     (   p   )                (   z   )       =       1   π            ∫   0   ∞              Q     n   ,     L   p       ±          (   k   )            (           cos        (   kz   )                 sin        (   kz   )             )             k             ,              and           (   71   )                   q   n     ±     ,     (   s   )                (   z   )       =       1   π            ∫   0   ∞              Q     n   ,     L   s       ±          (   k   )            (           cos        (   kz   )                 sin        (   kz   )             )               k     .                   (   72   )                                
     The corresponding constrained functionals that correspond to a particular mode is given by:                  E   0     (   ±   )       =       W     0   ,   g       (   ±   )       -       μ   o            ∫   0   ∞                λ   0   ±          (   k   )            [         F     0   ,     L   s         (   ±   )            (   k   )       +         R   p       R   s                I   0   ′          (     kR   p     )           I   0   ′          (     kR   s     )                F     0   ,     L   p         (   ±   )            (   k   )           ]               k           -       ∑     j   =   1     N                         λ   _       0      j     ±     ·     [           B   _     0   ±          (       ρ   j     ,     φ   j     ,     z   j       )            |       p   j     &lt;     R   p              -       B   _       0   ,   j     target         ]         -       μ   o            ∑     a   =   1       N   p                         λ     0   ,   a       ±     (   p   )                ∫   0   ∞              F     0   ,     L   p         (   ±   )            (   k   )            (           cos        (     k   a     (   s   )       )                 sin        (     k   a     (   s   )       )             )             k               -       μ   o            ∑     b   =   1       N   s                         λ     0   ,   b       ±     (   s   )                ∫   0   ∞              F     0   ,     L   s         (   ±   )            (   k   )            (             cos   (     k   b     (   s   )         }               sin        (     k   b     (   s   )       )             )             k                                            and           (   73   )                 E   n     (   ±   )       =       W     n   ,   g       (   ±   )       -       μ   o            ∫   0   ∞          k                       λ   n   ±          (   k   )            [         Q     n   ,     L   s         (   ±   )            (   k   )       +         R   p   2       R   s   2                I   n   ′          (     kR   p     )           I   n   ′          (     kR   s     )                Q     n   ,     L   p         (   ±   )            (   k   )           ]               k           -       ∑     j   =   1     N                         λ   _     nj   ±     ·     [           B   _     n   ±          (       ρ   j     ,     φ   j     ,     z   j       )            |       p   j     &lt;     R   p              -       B   _       n   ,   j     target         ]         -       μ   o            ∑     a   =   1       N   p                         λ     n   ,   a       ±     (   p   )                ∫   0   ∞              Q     n   ,     L   p         (   ±   )            (   k   )            (           cos        (     k   a     (   s   )       )                 sin        (     k   a     (   s   )       )             )             k               -       μ   o            ∑     b   =   1       N   s                         λ     n   ,   b       ±     (   s   )                ∫   0   ∞              Q     n   ,     L   s         (   ±   )            (   k   )            (             cos   (     k   b     (   s   )         }               sin        (     k   b     (   s   )       )             )             k                         (   74   )                                
     where λ η   ± (k) are Lagrange multipliers for the shielding conditions from equations (23), (24), λ η,a   i(ρ)  and λ η,b   i(s)  are the Lagrange multipliers for the points in the consistency conditions represented by equations (27) and (28), and equations (29) and (30) respectively, {overscore (B)} η,j   target  are the field constraints, and {overscore (λ)} η,j   ±  are their Lagrange multipliers. The variational principle for each mode of the functional E η   i  is given by:                    δ                   E   0   ±         δ                     F     0   ,     L   p       ±          (   k   )           =   0     ,         δ                   E   0   ±         δ                     F     0   ,     L   s       ±          (   k   )           =   0     ,         δ                   E   0   ±         δ                     λ   0   ±          (   k   )           =   0     ,           (   75   )                     ∂     E   0   ±         ∂       λ   _       0   ,   j     ±         =   0     ,         ∂     E   0   ±         ∂     λ     0   ,   a       ±     (   p   )             =   0     ,         ∂     E   0   ±         ∂     λ     0   ,   a       ±     (   s   )             =   0     ,           (   76   )                     δ                   E   n   ±         δ                     Q     n   ,     L   p       ±          (   k   )           =   0     ,         δ                   E   n   ±         δ                     Q     n   ,     L   s       ±          (   k   )           =   0     ,         δ                   E   n   ±         δ                     λ   n   ±          (   k   )           =   0     ,   ,     n   ≧   1     ,              and           (   77   )                     ∂     E   n   ±         ∂       λ   _       n   ,   j     ±         =   0     ,         ∂     E   n   ±         ∂     λ     n   ,   a       ±     (   p   )             =   0     ,         ∂     E   n   ±         ∂     λ     n   ,   a       ±     (   s   )             =   0     ,     n   ≧   1     ,           (   78   )                                
     which leads to the expressions for the Fourier transforms of the current densities in terms of the Lagrange multipliers {overscore (λ)} η,j   ± , λ η,a   ±(ρ) , λ η,b   ±(s)  as follows:                    F     0   ,     L   p       ±          (   k   )       =         ∑     j   =   1     N                         (     μ   o     )       -   1              (             δ            B   _     0   ±          (       ρ   j     ,     φ   j     ,     z   j       )              |       p   j     &lt;     R   p             δ                     F     0   ,     L   p       ±          (   k   )           /     g        (   k   )                W   0   ±          (     k   ,     R   p     ,     R   s       )         )     ·       λ   _       0   ,   j     ±           +       ∑     a   =   1       N   p                           (           cos        (     k   a     (   p   )       )                 sin        (     k   a     (   p   )       )             )         g        (   k   )              W   0   ±          (     k   ,     R   p     ,     R   s       )                λ     0   ,   a       ±     (   p   )             -         R   p       R   s                I   0   ′          (     kR   p     )           I   0   ′          (     kR   s     )                ∑     b   =   1       N   s                (           cos        (     k   b     (   s   )       )                 sin        (     k   b     (   s   )       )             )         g        (   k   )              W   0   ±          (     k   ,     R   p     ,     R   s       )                λ     0   ,   b       ±     (   s   )                                 and           (   79   )                     F     0   ,     L   s       ±          (   k   )       =       -       R   p       R   s                  I   0   ′          (     kR   p     )           I   0   ′          (     kR   s     )                F     0   ,     L   p       ±          (   k   )                             for                 axial                 coils     ,              and             (   80   )                     Q     n   ,     L   p       ±          (   k   )       =         ∑     j   =   1     N                         (     μ   o     )       -   1              (             δ            B   _     n   ±          (       ρ   j     ,     φ   j     ,     z   j       )              |       p   j     &lt;     R   p             δ                     F     n   ,     L   p       ±          (   k   )           /     g        (   k   )                W   n   ±          (     k   ,     R   p     ,     R   s       )         )     ·       λ   _       n   ,   j     ±           +       ∑     a   =   1       N   p                           (           cos        (     k   a     (   p   )       )                 sin        (     k   a     (   p   )       )             )         g        (   k   )              W   n   ±          (     k   ,     R   p     ,     R   s       )                λ     n   ,   a       ±     (   p   )             -         R   p   2       R   s   2                I   n   ′          (     kR   p     )           I   n   ′          (     kR   s     )                ∑     b   =   1       N   s                (           cos        (     k   b     (   s   )       )                 sin        (     k   b     (   s   )       )             )         g        (   k   )              W   n   ±          (     k   ,     R   p     ,     R   s       )                λ     n   ,   b       ±     (   s   )                                 and           (   81   )                     Q     n   ,     L   s       ±          (   k   )       =       -       R   p   2       R   s   2                  I   n   ′          (     kR   p     )           I   n   ′          (     kR   s     )                Q     n   ,     L   p       ±          (   k   )           ,     n   ≧   1             (   82   )                                
     for transverse coils. 
     Using the properties of the Bessel function, in order for the functional W g  to be finite, the function g(k) grows as k→∞. One possible choice for the function g(k) is: 
     
       
           g ( k )=exp( kγ ( R   s   −R   ρ ))  (83) 
       
     
     where γ is a positive constant parameter. 
     Substituting the expressions from equations (80)-(82) into equations (76) and (78) the following system of algebraic equations for the Lagrange multipliers is obtained:                      ∑     l   =   1     N            Z     pl   ,   n     ±          Λ     l   ,   n     ±         =     Q     p   ,   n     ±       ,        =     N   +     N   p     +     N   s                         where           (   84   )                     Λ     l   ,   n     ±     =     (                   λ   →       n   ,   j     ±               λ     n   ,   a       ±     (   p   )                         λ     n   ,   b       ±     (   s   )               )       ,       Q     p   ,   n     ±     =     (                   μ   o     -   1              B   →       n   ,   j     target               0                 0         )                       and           (   85   )                   -     Z     pl   ,   n     ±       =     (           C     ij   ,   n     ±           A     ia   ,   n     ±           B     ib   ,   n     ±               D         a   1        j     ,   n     ±           E         a   1        a     ,   n     ±           F         a   1        b     ,   n     ±               G         b   1        j     ,   n     ±           N         b   1        a     ,   n     ±           M         b   1        b     ,   n     ±           )       ,     n   ≥   0             (   86   )                                
     and where the matrices involved here are expressed through the integrals over k. Assuming that the matrix Z ρl,η   1  is not degenerate, the Lagrange multipliers are found in the following form:                Λ     l   ,   n     ±     =       ∑     p   =   1                   (     Z   n     -   1       )     lp   ±            Q     p   ,   n     ±     .                 (   87   )                                
     Once the Lagrange multipliers Λ l,η   ±  are found, the expressions for the Fourier transforms F 0,L     ρ     (L     s     )   ± (k) and Q n,L     ρ     (L     s     )   ± (k) are obtained using the equations (80)-(82). Then the equations (69)-(72) define those current distributions on both coils that lead to the exterior super shielding effect represented by equation (43). 
     II. Interior Super Shielding of the Magnetic Field 
     Next, the extremum calculation for the interior super shielding of the magnetic field is addressed. The energy-like functional considered in this case is                  W   g     =       ∑     n   =   0     ∞          (       W     n   ,   g       (   -   )       +     W     n   ,   g       (   +   )         )                       where           (   88   )                   W     0   ,   g       (   ±   )       =         μ   0     2            ∫   0   ∞            g        (   k   )              W   0          (     k   ,     R   p     ,     R   s       )              F     0   ,     L   s           (   ±   )        2            (   k   )                          k             ,           (   89   )                     W   0          (     k   ,     R   p     ,     R   s       )       =       -     R   s   2              K   0   ′          (     k                   R   s       )              D   0          (     k   ,     R   p     ,     R   s       )           ,           (   90   )                   W     n   ,   g       (   ±   )       =         μ   0     2            ∫   0   ∞            g        (   k   )              W   n          (     k   ,     R   p     ,     R   s       )              Q     n   ,     L   s           (   ±   )        2            (   k   )                          k             ,              and           (   91   )                     W   n          (     k   ,     R   p     ,     R   s       )       =       -         R   s   4          k   2         n   2                K   n   ′          (     k                   R   s       )              D   n          (     k   ,     R   p     ,     R   s       )           ,     n   ≥   1     ,           (   92   )                                
     and where g(k) is the same function that was introduced in the previous section. 
     The constrained functionals E η   (i)  have the same form as in equations (73) and (74) where exterior super shielding conditions from equations (23) and (24) are replaced by the interior super shielding conditions from equations (45) and (46), respectively. 
     SUPER SHIELDING OF THE ELECTRIC FIELD 
     I. Electric Field 
     First, consider a charge distribution σ({overscore (r)}) on a cylinder having a radius R and total length 2L. A Fourier series representation follows from the periodicity in φ and therefore                σ        (     r   →     )       =       ∑     n   =   0     ∞          (           σ   n          (   z   )            cos        (     n                 φ     )         +         κ   n          (   z   )            sin        (     n                 φ     )           )               (   93   )                                
     where σ η (z) and k η (z) are arbitrary functions of z. The end-condition is that there can be an abrupt change in the charge density, yielding 
     
       
         σ η (| z|&gt;L )= k   η (| z|&gt;L )=0  (94). 
       
     
     The scalar potential of the electric field is defined as                Φ        (     r   →     )       =       1     4      π                   ε   o              ∫       G        (       r   →     ,       r   →     ′       )            σ        (       r   →     ′     )                   r   →     ′                     (   95   )                                
     where the Green function G({overscore (r)},{overscore (r)} 1 ) is given by equation (11). Then the scalar potential can be written as                  Φ        (     r   →     )       =       ∑     i   =   1     2            ∑     n   =   0     ∞            Φ   n     (   i   )            (     r   →     )                           where           (   96   )                     Φ   n     (   1   )            (     r   →     )       =       R     π                   ε   o              cos        (     n                 φ     )              ∫   0   ∞              T   n          (     k   ,   ρ   ,   R     )              ∫     -   L     L              σ   n          (     z   ′     )            cos        (     k        (     z   -     z   ′       )       )                            z   ′                          k                               and           (   97   )                   Φ   n     (   2   )            (     r   →     )       =       R     π                   ε   o              sin        (     n                 φ     )              ∫   0   ∞              T   n          (     k   ,   ρ   ,   R     )              ∫     -   L     L              κ   n          (     z   ′     )            cos        (     k        (     z   -     z   ′       )       )                            z   ′                          k                       (   98   )                                
     and where the function T η ( k,ρ,R) is defined in equation (12). Then the electric field is given by                    E   →          (     r   →     )       =       ∑     i   =   1     2            ∑     n   =   0     ∞              E   →     n     (   i   )            (     r   →     )                           where           (   99   )                       E   →     n     (   i   )            (     r   →     )       =     -       ∇   →            Φ   n     (   i   )            (     r   →     )             ,           (   100   )                                
     in terms of the cylindrical coordinates can be presented in the following form (for each mode η):                    E     n   ,   ρ       (   1   )            (     r   →     )       =       -     R     π                   ε   o                cos        (     n                 φ     )              ∫   0   ∞          k                     T   n   ′          (     k   ,   ρ   ,   R     )              ∫     -   L     L              σ   n          (     z   ′     )            cos        (     k        (     z   -     z   ′       )       )                            z   ′               k                 ,           (   101   )                     E     n   ,   φ       (   1   )            (     r   →     )       =       R     π                   ε   o        ρ          n                   sin        (     n                 φ     )              ∫   0   ∞          k                     T   n          (     k   ,   ρ   ,   R     )              ∫     -   L     L              σ   n          (     z   ′     )            cos        (     k        (     z   -     z   ′       )       )                            z   ′               k                 ,              and           (   102   )                     E     n   ,   z       (   1   )            (     r   →     )       =       R     π                   ε   o              cos        (     n                 φ     )              ∫   0   ∞          k                     T   n          (     k   ,   ρ   ,   R     )              ∫     -   L     L              σ   n          (     z   ′     )            sin        (     k        (     z   -     z   ′       )       )                            z   ′               k                               where           (   103   )                   T   n   ′          (     k   ,   ρ   ,   R     )       =         θ        (     ρ   -   R     )              I   n          (     k                 R     )              K   n   ′          (     k                 ρ     )         +         θ   ′          (     R   -   ρ     )              I   n   ′          (     k                 ρ     )                K   n          (     k                 R     )       .                 (   104   )                                
     The expressions for {overscore (E)} η   (2) ({overscore (r)}) is obtained from equations (101)-(103) by the following substitutions: 
     
       
         cos(ηφ)→sin(ηφ), sin(ηφ)→−cos(ηφ), σ η ( z )→k η ( z )  (105). 
       
     
     To shield the system, consider two coaxial cylinders of arbitrary radii, R ρ  for the inner cylinder or primary structure and R s  for the outer cylinder or secondary structure, and arbitrary lengths of these cylinders, 2L ρ  and 2L s  respectively. On each cylinder, the charge distribution has the form given by equation (93). The total electric field then is a vector sum of the electric fields that correspond to that of the primary and secondary structures. In what follows all quantities that correspond to the primary structure will be labeled by an index &#39;ρ&#39;, while all quantities that correspond to the secondary structure will be labeled by an index &#39;s&#39;. 
     II. Exterior Super Shielding of the Electric Field 
     For excellent shielding of the electric field outside the secondary structure, the boundary conditions for the total electric field are that the tangential components of the total electric field vanish on the surface of the secondary structure, which is analytically represented as follows: 
     
       
           E   η,φ   tot(1) ( {overscore (r)} )| ρ=R     s     =E   η,z   tot,(1) ( {overscore (r)} )| ρ=R     s   =0,ρ≧0  (106). 
       
     
     A sufficient condition satisfying this equation is [for simplicity herein, in what follows only the electric field labeled by the index i=1 is considered; the corresponding expressions for the electric field labeled by the index i=2 are obtained using the substitutions of equation (105)]:                      ∑     n   ,     L   s       ±          (   k   )       +         R   p       R   s                I   n          (     k                   R   p       )           I   n          (     k                   R   s       )                ∑     n   ,     L   p       ±          (   k   )           =   0                   where           (   107   )                   ∑     n   ,   L     ±          (   k   )       =       ∫     -   L     L              σ   n   ±          (   z   )            (           cos        (   kz   )                 sin        (   kz   )             )                        z                 (   108   )                                
     and the sign &#39;+&#39; refers to an even function of z and the sign &#39;−&#39; refers to an odd function of z. The full implications of these conditions determine those charge densities on both coils that satisfy the boundary condition from equation (106). 
     An inverse Fourier transform of equation (107) together with the definitions from equation (108) lead to two types of additional consistency conditions.                    1   π            ∫   0   ∞            ∑     n   ,     L   p       ±            (   k   )          (           cos        (   kz   )                 sin        (   kz   )             )                        k             =       θ        (       L   p     -        z          )              σ   n     ±     ,     (   p   )                (   z   )                         and           (   109   )                     1   π            ∫   0   ∞            ∑     n   ,     L   s       ±            (   k   )          (           cos        (   kz   )                 sin        (   kz   )             )                        k             =       θ        (       L   s     -        z          )              σ   n     ±     ,     (   s   )                (   z   )                         where           (   110   )                   σ   n     ±     ,     (     p   ,   s     )                (   z   )       =       1   2            (         σ   n     (     p   ,   s     )            (   z   )       ±       σ   n     (     p   ,   s     )            (     -   z     )         )     .               (   111   )                                
     When |z|&gt;L ρ  the equation (109) represents some constraint for the charge density on the primary structure. When |z|&gt;L s  the equation (110) represents some constraint for the current density on the secondary structure. These consistency conditions represent the fact that the charge does not exist off the ends of the structures. That is to say, the charge densities are confined to finite cylinders. 
     For exterior shielding, the sufficiency condition (107) plus the consistency conditions (109)-(110) constitute the supershielding set of equations. By an inverse Fourier transform of (107), subject to (109)-(110), and an application of the Fourier convolution theorem, the supershielding set of equations is found to be equivalent to                    Θ        (       L   s     -        z          )              σ   n     ±     ,     (   s   )                (   z   )         +       1   π            ∫     -     L   p         L   p                r   e          (     z   -   y     )              σ   n     ±     ,     (   p   )                (   y   )               y             =     0                 where             (107a)                   r   e          (     z   -   y     )       =       ∫   0   ∞              R   e          (   k   )          cos                   k        (     z   -   y     )               k                 (107b)                                
     which is the Fourier transform of the coefficient in (107)                  R   e          (   k   )       =         R   p       R   s                  I   n          (     kR   p     )           I   n          (     kR   s     )         .               (107c)                                
     The shielding conditions from equation (107) are now used to write the expressions for all components of the total electric field. Inside an infinite cylinder with a radius given by that of the primary structure radius (i.e. ρ&lt;R ρ ), the following is obtained:                  E     n   ,   ρ       tot   ,     (   1   )         =       -       R   p       π                   ε   o                cos        (     n                 φ     )              ∫   0   ∞          k                     I   n   ′          (     k                 ρ     )              M   n          (     k   ,     R   p     ,     R   s       )              ∫     -     L   p         L   p                σ     (   p   )            (     z   ′     )            cos        (     k        (     z   -     z   ′       )       )                            z   ′                          k                 ,           (   112   )                   E     n   ,   φ       tot   ,     (   1   )         =         R   p       π                   ε   o        ρ            sin        (     n                 φ     )              ∫   0   ∞              I   n          (     k                 ρ     )              M   n          (     k   ,     R   p     ,     R   s       )              ∫     -     L   p         L   p                σ     (   p   )            (     z   ′     )            cos        (     k        (     z   -     z   ′       )       )                            z   ′                          k                 ,              and           (   113   )                   E     n   ,   z       tot   ,     (   1   )         =         R   p       π                   ε   o              cos        (     n                 φ     )              ∫   0   ∞          k                     I   n          (     k                 ρ     )              M   n          (     k   ,     R   p     ,     R   s       )              ∫     -     L   p         L   p                σ     (   p   )            (     z   ′     )            sin        (     k        (     z   -     z   ′       )       )                            z   ′                          k                               where           (   114   )                   M   n          (     k   ,     R   p     ,     R   s       )       =         K   n          (     k                   R   p       )       -           I   n          (     k                   R   p       )           I   n          (     k                   R   s       )                  K   n          (     k                   R   s       )       .                 (   115   )                                
     In the region between two cylinders having infinite lengths and whose radii are the same as those of the primary and secondary structures (i.e. R ρ &lt;ρ&lt;R s ), the expressions for all components of the total electric field are                  E     n   ,   ρ       tot   ,     (   1   )         =       -       R   p       π                   ε   o                cos        (     n                 φ     )              ∫   0   ∞          k            I   n          (     k                   R   p       )           I   n          (     k                   R   s       )         ×     {           K   n   ′          (     k                 ρ     )              I   n          (     k                   R   s       )         -         I   n   ′          (     k                 ρ     )              K   n          (     k                   R   s       )           }            ∫     -     L   p         L   p                σ     (   p   )            (     z   ′     )            cos        (     k        (     z   -     z   ′       )       )                 z   ′               k                 ,           (   116   )                   E     n   ,   φ       tot   ,     (   1   )         =         R   p       π                   ε   o        ρ            cos        (     n                 φ     )              ∫   0   ∞                I   n          (     k                   R   p       )           I   n          (     k                   R   s       )         ×     {           K   n          (     k                 ρ     )              I   n          (     k                   R   s       )         -         I   n          (     k                 ρ     )              K   n          (     k                   R   s       )           }            ∫     -     L   p         L   p                σ     (   p   )            (     z   ′     )            sin        (     k        (     z   -     z   ′       )       )                 z   ′               k                 ,              and           (   117   )                 E     n   ,   z       tot   ,     (   1   )         =       -       R   p       π                   ε   o                sin        (     n                 φ     )              ∫   0   ∞          k            I   n          (     k                   R   p       )           I   n          (     k                   R   s       )         ×     {           K   n          (     k                 ρ     )              I   n          (     k                   R   s       )         -         I   n          (     k                 ρ     )              K   n          (     k                   R   s       )           }            ∫     -     L   p         L   p                σ     (   p   )            (     z   ′     )            cos        (     k        (     z   -     z   ′       )       )                 z   ′                 k     .                       (   118   )                                
     In the region outside an infinite cylinder whose radius is the same as that of the secondary structure (i.e. ρ&gt;R s ), the total electric field exactly vanishes, which is expressed as follows: 
     
       
           E   η,ρ   tot,(1) ( {overscore (r)} )= E   η,φ   tot,(1) ( {overscore (r)} )= E   η,z   tot,(1) ( {overscore (r)} )=0  (119). 
       
     
     Equation (119) is the principal result of this section. The shielding condition from equation (107) leads to the complete cancellation of the electric field in the space outside the infinite cylinder coincident with the secondary structure. Moreover, the field components inside this infinite cylinder all vanish reasonably rapidly in the z-direction as the distance increases from the ends of the structures. 
     III. Interior Super Shielding of the Electric Field 
     For excellent shielding of the electric field inside the primary structure, the following boundary conditions are used for the total electric field: 
     
       
           E   η,φ   tot ( {overscore (r)} )| ρ=R     ρ     =E   η,z   tot ( {overscore (r)} )| ρ=R     ρ   0,η≧0  (120) 
       
     
     which constrains the total electric field such that tangential components thereof vanish at the surface of the primary structure. A sufficient condition satisfying this equation is [again for simplicity herein, in what follows, only the electric field labeled by the index i=1 is considered; the corresponding expressions for the electric field labeled by the index i=2 are obtained using the substitutions from equation (105)]:                    ∑     n   ,     L   p       ±          (   k   )       +         R   s       R   p                K   n          (     k                   R   s       )           K   n          (     k                   R   p       )                ∑     n   ,     L   s       ±          (   k   )           =   0.           (   121   )                                
     The consistency conditions of this equation have exactly the same form as equations (109) and (110). 
     For interior shielding, the sufficiency condition (121) plus the consistency conditions of the same form as (109)-(110) constitute the supershielding set of equations. By an inverse Fourier transform of (121), subject to the consistency conditions, and an application of the Fourier convolution theorem, the supershielding set of equations is found to be equivalent to                    Θ        (       L   p     -        z          )              σ   n     ±     ,     (   p   )                (   z   )         +       1   π            ∫     -     L   s         L   s                r   i          (     z   -   y     )              σ   n     ±     ,     (   s   )                (   y   )               y             =     0                 where             (121a)                   r   i          (     z   -   y     )       =       ∫   0   ∞              R   i          (   k   )          cos                   k        (     z   -   y     )               k                 (121b)                                
     which is the Fourier transform of the coefficient in (121)                  R   i          (   k   )       =         R   s       R   p                  K   n          (     kR   s     )           K   n          (     kR   p     )         .               (121c)                                
     The interior shielding condition from equation (121) is now used to write the expressions for all components of the total electric field. Inside the infinite cylinder with a radius given by the primary structure radius (i.e. ρ&lt;R ρ ), the total electric field exactly vanishes, which is represented by the following: 
     
       
           E   η,ρ   tot,(1) ( {overscore (r)} )= E   η,φ   tot,(1) ( {overscore (r)} )= E   η,z   tot,(1) ( {overscore (r)} )=0  (122). 
       
     
     In an infinite region between infinitely long cylinders whose radii are the same as those of the primary and secondary structures (i.e. R ρ &lt;ρ&lt;R s ), the expressions for all components of the total electric field are                  E     n   ,   ρ       tot   ,     (   1   )         =       -       R   s       π                   ε   o                cos        (     n                 φ     )              ∫   0   ∞          k            K   n          (     k                   R   s       )           K   n          (     k                   R   p       )         ×     {           I   n   ′          (     k                 ρ     )              K   n          (     k                   R   p       )         -         K   n   ′          (     k                 ρ     )              I   n          (     k                   R   p       )           }            ∫     L   s       L   s                σ     (   s   )            (     z   ′     )            cos        (     k        (     z   -     z   ′       )       )                 z   ′               k                 ,           (   123   )                   E     n   ,   φ       tot   ,     (   1   )         =         R   s       π                   ε   o        ρ            sin        (     n                 φ     )              ∫   0   ∞                K   n          (     k                   R   s       )           K   n          (     k                   R   p       )         ×     {           I   n          (     k                 ρ     )              K   n          (     k                   R   p       )         -         K   n          (     k                 ρ     )              I   n          (     k                   R   p       )           }            ∫     L   s       L   s                σ     (   s   )            (     z   ′     )            cos        (     k        (     z   -     z   ′       )       )                 z   ′               k                 ,              and           (   124   )                 E     n   ,   z       tot   ,     (   1   )         =       -       R   s       π                   ε   o                cos        (     n                 φ     )              ∫   0   ∞                K   n          (     k                   R   s       )           K   n          (     k                   R   p       )         ×     {           I   n          (     k                 ρ     )              K   n          (     k                   R   p       )         -         K   n          (     k                 ρ     )              I   n          (     k                   R   p       )           }            ∫     L   s       L   s                σ     (   s   )            (     z   ′     )            sin        (     k        (     z   -     z   ′       )       )                 z   ′                 k     .                       (   125   )                                
     In the region outside the infinite cylinder whose radius is the same as that of the secondary structure (i.e. ρ&gt;R s ), all components of the total electric field are found to be                  E     n   ,   ρ       tot   ,     (   1   )         =       -       R   s       π                   ε   o                cos        (     n                 φ     )              ∫   0   ∞          k                     K   n   ′          (     k                 ρ     )              E   n          (     k   ,     R   p     ,     R   s       )              ∫     L   s       L   s                σ     (   s   )            (     z   ′     )            cos        (     k        (     z   -     z   ′       )       )                 z   ′               k                 ,           (   126   )                   E     n   ,   φ       tot   ,     (   1   )         =         R   p       π                   ε   o   -        ρ            sin        (     n                 φ     )              ∫   0   ∞                         K   n          (     k                 ρ     )              E   n          (     k   ,     R   p     ,     R   s       )              ∫     L   s       L   s                σ     (   s   )            (     z   ′     )            cos        (     k        (     z   -     z   ′       )       )                 z   ′               k                 ,              and           (   127   )                   E     n   ,   z       tot   ,     (   1   )         =       -       R   s       π                   ε   o                cos        (     n                 φ     )              ∫   0   ∞              K   n          (     k                 ρ     )              E   n          (     k   ,     R   p     ,     R   s       )              ∫     L   s       L   s                σ     (   s   )            (     z   ′     )            sin        (     k        (     z   -     z   ′       )       )                 z   ′               k                               where           (   128   )                   E   n          (     k   ,     R   p     ,     R   s       )       =         I   n          (     kR   s     )       -           K   n          (     kR   s     )           K   n          (     kR   p     )                  I   n          (     kR   p     )       .                 (   129   )                                
     Equation (122) is the principal result of this section. The shielding condition from equation (121) and its consistency condition have led to the complete cancellation of the electric field in the space inside the infinite cylinder coincident with the primary structure. Again, the field components inside this infinite cylinder all vanish reasonably rapidly in the z-direction as the distance increases from the ends of the structures. 
     Functional Analysis 
     A definite solution for the charge distributions on the primary and secondary structures are obtained by posting an energy-like functional and incorporating the shielding constrains, along with other design constraints, in an extremum calculation. 
     I. Exterior Super Shielding of the Electric Field 
     In this case, the energy-like functional introduced take the form:                  W   g     =       ∑     n   =   0     ∞          (       W     n   ,   g       (   -   )       +     W     n   ,   g       (   +   )         )                       where           (   130   )                   W     n   ,   g       (   ±   )       =       1     2        ε   o                ∫   0   ∞            g        (   k   )              W   n          (     k   ,     R   p     ,     R   s       )              ∑     n   ,     L   p           (   ±   )        2              (   k   )             k                               and           (   131   )                   W   n          (     k   ,     R   p     ,     R   s       )       =       R   p   2            I   n          (     kR   p     )              M   n          (     k   ,     R   p     ,     R   s       )                 (   132   )                                
     and where g(k) is a regularization function for controlling the energy and the current behavior. For g=1, on the other hand, the functional W g  reduces to a Parseval&#39;s representation of the energy functional when the shielding conditions have been implemented. This is represented as follows: 
     
       
           W   g=1   =½∫σ(   {overscore (r)} )φ( {overscore (r)} ) d{overscore (r)}   (133). 
       
     
     To find those current distributions that lead to the exterior super shielding of the electric field, the same constraint functional technique that was used for exterior super shielding of the magnetic field is used. The constraints are the exterior shielding constraints from equation (107), its consistency conditions, and the desired values of the electric field at some particular points. 
     II. Interior super Shielding of the Electric Field 
     In this case the introduced energy-like functional takes the form:                  W   g     =       ∑     n   =   0     ∞          (       W     n   ,   g       (   -   )       +     W     n   ,   g       (   +   )         )                       where           (   134   )                   W     n   ,   g       (   ±   )       =       1     2        ε   o                ∫   0   ∞            g        (   k   )              W   n          (     k   ,     R   p     ,     R   s       )              ∑     n   ,     L   s           (   ±   )        2              (   k   )             k                               and           (   135   )                   W   n          (     k   ,     R   p     ,     R   s       )       =       R   s   2            K   n          (     kR   p     )              E   n          (     k   ,     R   p     ,     R   s       )                 (   136   )                                
     and where g(k) is a regularization function for controlling the energy and the current behavior. For g=1, on the other hand, the functional W g  reduces to a Parseval&#39;s representation of the energy functional when the shielding conditions have been implemented. Again, this takes the form: 
     
       
           W   g=1   =½∫σ(   {right arrow over (r)} )Φ( {right arrow over (r)} ) d{right arrow over (r)}   (137). 
       
     
     To find those charge distributions that lead to the interior super shielding of the electric field, the same constraint functional technique that was used for interior super shielding of the magnetic field is used. The constraints are the interior shielding constraints from equation (121), its consistency conditions, and the desired values of the electric field at some particular points. 
     EXAMPLES OF THE SUPER SHIELDING FOR THE MAGNETIC FIELD 
     I. B 0  Coil Example 
     In order to illustrate the method described above, consider the super shielding for the magnetic field produced by a azimuthally symmetric coil system. In this case, as it was mentioned previously, only the term {overscore (J)} 0 (z) with an even function f 0 (z) is employed. The geometrical dimensions of the primary coil  60  are R ρ =0.3528 m for the radius and 2L ρ =1.1 m for the total length, while for the secondary coil  62  has a radius R s =0.417151 m and 2L s =1.1 m for the total length. Only one constraint, set out in table 1 below, is imposed for the B z  component of the magnetic field. 
     
       
         
               
               
               
               
               
             
           
               
                   
                 TABLE 1 
               
               
                   
                   
               
               
                   
                 ρ (in m) 
                 φ (in rad) 
                 z (in m) 
                 B z  (in mT) 
               
               
                   
                   
               
             
             
               
                   
                 0.000 
                 0.000 
                 0.00 
                 1.87 
               
               
                   
                   
               
             
          
         
       
     
     Six constraint points for the J φ -components of the current densities outside the primary and secondary coils, set out in tables 2 and 3 respectively, are used in order to satisfy the consistency conditions. 
     
       
         
               
               
               
             
           
               
                 TABLE 2 
               
               
                   
               
               
                 ρ = R ρ  (in m) 
                 z (in m) 
                 J φ   ρ  (in A/m) 
               
               
                   
               
             
             
               
                 0.3528 
                 0.5709 
                 0.0 
               
               
                 0.3528 
                 0.5918 
                 0.0 
               
               
                 0.3528 
                 0.6127 
                 0.0 
               
               
                 0.3528 
                 0.6336 
                 0.0 
               
               
                 0.3528 
                 0.6545 
                 0.0 
               
               
                 0.3528 
                 0.6754 
                 0.0 
               
               
                   
               
             
          
         
       
     
     
       
         
               
               
               
             
           
               
                 TABLE 3 
               
               
                   
               
               
                 ρ = R s  (in m) 
                 z (in m) 
                 J φ   s  (in A/m) 
               
               
                   
               
             
             
               
                 0.417151 
                 0.5709 
                 0.0 
               
               
                 0.417151 
                 0.5918 
                 0.0 
               
               
                 0.417151 
                 0.6127 
                 0.0 
               
               
                 0.417151 
                 0.6336 
                 0.0 
               
               
                 0.417151 
                 0.6545 
                 0.0 
               
               
                 0.417151 
                 0.6754 
                 0.0 
               
               
                   
               
             
          
         
       
     
     The discrete current density for both coils, with a resolution of 1 mm, is found by numerical evaluation of the analytical expressions. Excellent shielding level occurs when the parameter γ=3.8 and only six constraint points outside the primary and secondary coils  60  and  62 , respectively, that correspond to the consistency conditions, are equally spaced with a step of 0.0209 m. With reference to FIG. 2, the results for the z-dependence of the J φ -components of the current density on both coils are shown. From the Biot-Savart law, the calculated current densities are used to find all components of the magnetic field. FIG. 3 illustrates the z-dependence of the B z -component of the magnetic field that corresponds to the contribution of the primary coil  60  only and the total magnetic field inside the primary coil  60  at the radial distance ρ=0.0 m. FIGS. 4A and 4B show the z dependence of the B ρ - and B z -components of the magnetic field within the coil system (0.0&lt;z&lt;0.55 m) at the radial distance ρ=0.48 mm (outside the secondary coil). FIG. 4A shows the contribution of the primary current to all components of the magnetic field of order 1.5 mT, and FIG. 4B shows the total shielded magnetic field of order 0.06 μT. FIGS. 5A and 5B illustrate the z dependence of the magnetic field outside the coil system (0.55&lt;z&lt;1.0 m) at the radial distance ρ=0.48 m. FIG. 5A shows the contribution of the primary current to all components of the magnetic field and FIG. 5B shows the total magnetic field. The shielding defined as follows:          shielding  factor =             max                   B   pr                     (     z   ,     ρ   =     0.48                 m         )         max                   B   total                     (     z   ,     ρ   =     0.48                 m         )                                
     in the case under consideration is equal to 2.3×10 4 . 
     FIGS. 6A and 6B present the ρ dependence of the magnetic field at z=0.1 m. FIG. 6A shows the contribution of the primary current to all component of the magnetic field. It illustrates that the B ρ -component of the magnetic field is continuous at ρ=R ρ , while the B z -component is discontinuous at the same point. FIG. 6B shows the ρ dependence of the total magnetic field. It illustrates that the B ρ -component of the magnetic field is continuous at ρ=R ρ  and ρ=R s , while the B z -component is discontinuous at the same points. All these results are consistent with the boundary conditions for the magnetic field behavior. 
     With reference to FIG. 7, An additional comparison is made when the length of the secondary coil L s  is twice as long as the length L ρ  of the primary coil. In this case, the field constraint is per Table 1, the primary current constraints are per Table 2, and the shield current constraints are per Table 3 but shifted further out along z by 55 cm. In this case, the current density on the primary coil is positive everywhere, while the current density on the secondary coil is negative everywhere. Comparing FIG.  2  and FIG. 7, in the case when the length of the primary coil is comparable with the length of the secondary coil, the current densities on both the primary and secondary coils have negative and positive lobes. In the case when the length of the secondary coil is considerably longer than the length of the primary coil, the current densities for both the primary and secondary coils do not go through a sign change within the length of the coils. 
     II. Z-Gradient Coil Example 
     Next, consider the example of an extra short z-gradient coil with an active shield. The gradient coil system consists of two coaxial cylindrical formers whose radii are: R ρ =0.3900 m and R s =0.4725 m, and half-lengths are L ρ =0.35 m and L s =0.4 m, respectively. Note that for this coil system, the ratios R ρ /L ρ  and R s /L s are both greater than one R ρ /L ρ =1.1143 and R s /L s =1.18125. In designing the z-gradient coil, four field constraints were chosen. These constraints are set forth in Table 4 below. 
     
       
         
               
               
               
               
             
           
               
                 TABLE 4 
               
               
                   
               
               
                 ρ (in m) 
                 φ (in rad) 
                 z (in m) 
                 B z  (in mT) 
               
               
                   
               
             
             
               
                 0.000 
                 0.000 
                 0.001 
                 1.78e-5 
               
               
                 0.185 
                 0.000 
                 0.001 
                 1.99329e-5 
               
               
                 0.225 
                 0.000 
                 0.001 
                 1.99329e-5 
               
               
                 0.000 
                 0.000 
                 0.220 
                 0.00325 
               
               
                   
               
             
          
         
       
     
     Nine constraint points for the J φ -components of the current densities outside the primary and secondary coils are used in order to satisfy the consistency conditions. The nine constraint points outside the primary and secondary coils are set forth in following Tables 5 and 6, respectively. 
     
       
         
               
               
               
             
           
               
                 TABLE 5 
               
               
                   
               
               
                 ρ = R ρ  (in m) 
                 z (in m) 
                 J φ   ρ  (in A/m) 
               
               
                   
               
             
             
               
                 0.3900 
                 0.354 
                 0.0 
               
               
                 0.3900 
                 0.374 
                 0.0 
               
               
                 0.3900 
                 0.397 
                 0.0 
               
               
                 0.3900 
                 0.419 
                 0.0 
               
               
                 0.3900 
                 0.442 
                 0.0 
               
               
                 0.3900 
                 0.464 
                 0.0 
               
               
                 0.3900 
                 0.487 
                 0.0 
               
               
                 0.3900 
                 0.509 
                 0.0 
               
               
                 0.3900 
                 0.531 
                 0.0 
               
               
                   
               
             
          
         
       
     
     
       
         
               
               
               
             
           
               
                 TABLE 6 
               
               
                   
               
               
                 ρ = R s  (in m) 
                 z (in m) 
                 J φ   s  (in A/m) 
               
               
                   
               
             
             
               
                 0.4725 
                 0.407 
                 0.0 
               
               
                 0.4725 
                 0.422 
                 0.0 
               
               
                 0.4725 
                 0.443 
                 0.0 
               
               
                 0.4725 
                 0.465 
                 0.0 
               
               
                 0.4725 
                 0.487 
                 0.0 
               
               
                 0.4725 
                 0.508 
                 0.0 
               
               
                 0.4725 
                 0.530 
                 0.0 
               
               
                 0.4725 
                 0.551 
                 0.0 
               
               
                 0.4725 
                 0.572 
                 0.0 
               
               
                   
               
             
          
         
       
     
     The discrete current density for both coils, with resolution of 1 mm, is found by numerical evaluation of the analytical expressions. Excellent shielding levels occur when the parameter γ=0.30. In FIGS. 8A and 8B, the results for the z-dependence of the J φ -components of the current density on both coils are shown. From the Biot-Savart law, the calculated current densities are used to find all components of the magnetic field. 
     The gradient strength at the isocenter of the coil system is G z =17.73 mT/m. FIGS. 9A and 9B illustrate the z dependence of the magnetic field outside the coil system at the radial distance ρ=0.55 m. Both graphs illustrate this behavior within the length of the coil system (0.0&lt;z&lt;0.4 m). In particular, FIG. 9A shows the total B z  component of the magnetic field. In this region, the shielding factor is equal to 199.94. FIG. 10 represents the behavior of the B z  components of the total magnetic field as a function of z inside the primary coil at ρ=0.0 m. FIG. 11A illustrates the linearity, and FIG. 11B, the uniformity of the B z  component of the total magnetic field. The linearity is defined as          L   in     =             G   z          (     z   ,                ρ   =   0.0       )       -     G   z         G   z       ×   100      %                            
     where G z (z,ρ=0.0) is the global gradient strength of the magnetic field. The uniformity is defined as        Unif   =             B   z          (       z   =   0.001     ,              ρ     )       -       B   z          (       z   =   0.001     ,                ρ   =   0.0       )             B   z          (       z   =   0.001     ,                ρ   =   0.0       )         ×   100      %                            
     This uniformity illustrates how the gradient strength of the B z  component of the magnetic field changes as ρ increases. 
     III. Electric Field Shielding Example 
     The surface structure system that is considered consists of two concentric cylinders: an inner cylinder whose radius is R ρ =0.39 m and total length is 2L ρ =1.1 m, and an outer cylinder whose radius is R s =0.43 m and total length is 2L s =1.1 m. This representative case uses charge distributions that produce non-zero longitudinal electric fields in the isocenter of the surface structure system. This means that the charge distributions are antisymmetric functions of z. One field constraint for the electric field at the isocenter is used. The cosntraint is shown in the following Table 7. 
     
       
         
               
               
               
               
               
             
           
               
                   
                 TABLE 7 
               
               
                   
                   
               
               
                   
                 ρ (in m) 
                 φ (in rad) 
                 z (in m) 
                 E z  (in V/m) 
               
               
                   
                   
               
             
             
               
                   
                 0.0 
                 0.000 
                 0.00 
                 840.0 
               
               
                   
                   
               
             
          
         
       
     
     Nine contraint points for the charge distributions O (ρ) (z) and O (s) (z) outside the primary and secondary surface structures are used in order to satisfy the consistency conditions employed for the method of supershielding. These are given in Tables 8 and 9, respectively. 
     
       
         
               
               
               
             
           
               
                 TABLE 8 
               
               
                   
               
               
                 ρ = R ρ  (in m) 
                 z (in m) 
                 o (ρ)  (in A/m) 
               
               
                   
               
             
             
               
                 0.39 
                 0.5600 
                 0.0 
               
               
                 0.39 
                 0.6120 
                 0.0 
               
               
                 0.39 
                 0.7016 
                 0.0 
               
               
                 0.39 
                 0.8058 
                 0.0 
               
               
                 0.39 
                 0.9207 
                 0.0 
               
               
                 0.39 
                 1.0444 
                 0.0 
               
               
                 0.39 
                 1.1755 
                 0.0 
               
               
                 0.39 
                 1.3131 
                 0.0 
               
               
                 0.39 
                 1.4565 
                 0.0 
               
               
                   
               
             
          
         
       
     
     
       
         
               
               
               
             
           
               
                 TABLE 9 
               
               
                   
               
               
                 ρ = R s  (in m) 
                 z (in m) 
                 o (s)  (in A/m) 
               
               
                   
               
             
             
               
                 0.43 
                 0.5600 
                 0.0 
               
               
                 0.43 
                 0.6120 
                 0.0 
               
               
                 0.43 
                 0.7016 
                 0.0 
               
               
                 0.43 
                 0.8058 
                 0.0 
               
               
                 0.43 
                 0.9207 
                 0.0 
               
               
                 0.43 
                 1.0444 
                 0.0 
               
               
                 0.43 
                 1.1755 
                 0.0 
               
               
                 0.43 
                 1.3131 
                 0.0 
               
               
                 0.43 
                 1.4565 
                 0.0 
               
               
                   
               
             
          
         
       
     
     Discrete charge densities for both surface structures, with a resolution of 1 mm, are found by numerical evaluation of the analytical expressions. With the experience gained in the previous work, a significant degree of shielding is quickly obtained even prior to any systematic tweaking of the numerical method parameters. For example, it is seen below that shielding corresponding to field reductions by a factor of 3500 is obtained in the case where the cutoff parameter is chosen to be γ=15.0. 
     FIGS. 12-17 illustrate the numerical results. In FIG. 12, the results for the z-dependence of the charge density on both surface structures are shown. It is seen that the charge densities on both surface structures have different polarities as z changes from z=0 to the edges of the surface structures. The calculated charge densities are used to find all components of the electric field. FIG. 13A shows the z dependence of the two non-zero components of the total electric field inside the primary surface structure. In particular, it illustrates the behavior of the E z  component at the radial distance ρ=0.0 m and the E ρ  component at the radial distance of ρ=0.1 m. FIG. 13B presents the z dependence of both E z  and E ρ  components of the electric field between the primary and secondary surface structures at the radial distance ρ=0.4 m. It is seen from these figures that the electric field is stronger between the surface structures rather than inside the primary surface structure. This result is consistent with principles of shielding. In order to illustrate the shielding effect, the contributions of the charge density on the primary surface structure to all components of the electric field outside the surface structure system (at the radial distance ρ=0.55 m) are isolated and compared with the corresponding behavior of the total electric field at the same radial distance. FIG. 14A shows this behavior for the primary electric field and FIG. 14B shows similar behavior for the total shielded electric field. The &#39;shielding factor&#39; is defined as          shielding  factor =             max                   E     (   pr   )                       (     z   ,     ρ   =   0.55       )         max                   E     (   total   )                       (     z   ,     ρ   =   0.55       )                                
     In the case under consideration, the shielding factor=6.7×10 3 . 
     FIGS. 15A-15C and  16  show the radial dependence of all components of the total electric field at various z-positions. FIGS. 15A and 15C specifically illustrate the radial ρ dependence of the electric field within the surface structure system. They show that the E z  component of the electric field is continuous as ρ crosses the surface structure former, while the E ρ  component of the electric field is discontinuous in accordance with the presence of a surface charge distribution (Gauss&#39;s law). FIG. 16 shows the ρ dependence of both components of the electric field at the position z=0.6 m which is beyond the length of the surface structure system. Here both components of the electric fields are continuous as the radial distance ρ crosses the two cylinders of the surface structure system. Note that beyond ρ=0.43 m which is the radius of the secondary surface structure, both components of the electric field are well shielded. 
     An additional comparison is made when the length of the secondary surface structure L s  is almost twice as long as the length L ρ  of the primary surface structure. In this case the charge density on both surface structures still change sign as z changes from z=0.0 m to the edges of the surface structures as it is shown in FIG.  17 . Recall that in the case of the magnetic field, the currents&#39; densities on both surface structures do not change their signs in the case when the secondary surface structure is much longer than the primary surface structure. As it is seen from FIGS. 13A and 14B, the electric field is concentrated mostly within the surface structure system and, outside the surface structures, falls off rapidly as z increases. 
     The invention has been described with reference to the preferred embodiments. Obviously, modifications and alterations will occur to others upon reading and understanding the preceding detailed description. It is intended that the invention be construed as including all such modifications and alterations insofar as they come within the scope of the appended claims or the equivalents thereof.