Quantifying MRI frequency shifts due to structures with anisotropic magnetic susceptibility using pyrolytic graphite sheet

Magnetic susceptibility is an important source of contrast in magnetic resonance imaging (MRI), with spatial variations in the susceptibility of tissue affecting both the magnitude and phase of the measured signals. This contrast has generally been interpreted by assuming that tissues have isotropic magnetic susceptibility, but recent work has shown that the anisotropic magnetic susceptibility of ordered biological tissues, such as myelinated nerves and cardiac muscle fibers, gives rise to unexpected image contrast. This behavior occurs because the pattern of field variation generated by microstructural elements formed from material of anisotropic susceptibility can be very different from that predicted by modelling the effects in terms of isotropic susceptibility. In MR images of tissue, such elements are manifested at a sub-voxel length-scale, so the patterns of field variation that they generate cannot be directly visualized. Here, we used pyrolytic graphite sheet which has a large magnetic susceptibility anisotropy to form structures of known geometry with sizes large enough that the pattern of field variation could be mapped directly using MRI. This allowed direct validation of theoretical expressions describing the pattern of field variation from anisotropic structures with biologically relevant shapes (slabs, spherical shells and cylindrical shells).


Theory
PGS is strongly diamagnetic (with magnetic susceptibility χ ⊥ ) when a magnetic field is applied perpendicular to the sheet, but only weakly diamagnetic (with magnetic susceptibility χ || ), when the field is applied parallel to the sheet 24 . This susceptibility anisotropy can be described by using a cylindrically-symmetric tensor characterize the isotropic and anisotropic parts of the susceptibility, respectively. If a structure made from PGS is exposed to an applied magnetic field, H, the magnetization, M, is given by = χ M H and the field perturbation, H d , due to this magnetization can be described as the gradient of a scalar potential: where the volume integral is taken over the volume of the PGS structure and the surface integral extends over the internal and external surfaces of the structure. If the structure is embedded in a medium then the value of the isotropic susceptibility of the PGS that we use in solving Eq. (2) should be taken relative to that of the medium. We are interested in the NMR frequency perturbation generated when PGS structures are exposed to a spatially uniform magnetic field, Ĥ z 0 (where B 0 = μ 0 H 0 ) and since χ || , χ ⊥ ≪ 1 we can neglect the effect of terms of order χ 2 and above, by setting = χM H z 0 when calculating Φ(r), and also only need to consider the component of the field perturbation that is parallel to the applied field (i.e. along z). As a result of the susceptibility anisotropy, the magnitude and direction of the induced magnetization depend on the orientation of the PGS sheet with respect to the uniform field direction. Considering a small volume, ΔV, of PGS oriented such that the normal to the sheet lies in the x-z plane and makes an angle Θ to the z-axis, the induced magnetization is given by is generates a spatially varying perturbation of the z-component of the magnetic flux density outside the PGS which is described by 8 sin2 (sin2 cos ) 4 (1 3 cos2 ) (3 cos 1) , , with r, θ and φ representing standard spherical polar co-ordinates. Compared with a volume element of material of isotropic susceptibility, the effect of the susceptibility anisotropy is to produce an orientation-dependent modulation of the strength of the dipole field due to magnetization aligned with the applied field and to introduce an additional field perturbation due to a dipole moment that is perpendicular to the field. The amplitude of the perpendicular dipole moment is maximized when the normal to the PGS is oriented at 45° to the field.
The susceptibility anisotropy gives rise to interesting behavior when PGS is shaped to form structures in which the direction of the principal axis of the susceptibility tensor varies in space. In particular, a thin spherical shell Scientific REPORTs | (2018) 8:6259 | DOI:10.1038/s41598-018-24650-2 of PGS of radius, R, and thickness, t ≪ R, arranged such that the axis of strong diamagnetism is radial, produces a field perturbation of the form where r and θ in represent spherical polar co-ordinates (see Supplementary Material for derivation of Eqs (5 and 6)).
Similarly, a long cylindrical annulus of radius, R, and thickness, t ≪ R oriented at angle α to the applied B 0 -field produces a field perturbation that is described by where ρ, φ and z′ represent cylindrical polar co-ordinates defined relative to the cylindrical annulus. Figure 1 shows simulations of the patterns of field variation (in ppm) produced by thin spherical (Fig. 1A) and cylindrical ( Fig. 1B) shells of material that were generated using Eqs (5) and (6). The field perturbations due to the isotropic and anisotropic components of the susceptibility are shown separately and then combined for both geometries. For these simulations, we used values of t = 25 μm, R = 10 mm, χ I = −135 ppm, and χ A = −260 ppm, which are representative of the PGS structures used in the experiments reported below, and it is assumed that the cylindrical shell is oriented perpendicular to B 0 .
For both spherical and cylindrical shells, the anisotropic component of the magnetic susceptibility generates a uniform field offset inside the structure, whose magnitude depends on the ratio of the shell's thickness to its radius 7 . The isotropic magnetic susceptibility generates fields of the well-known form outside the shell (dipolar for the spherical shell, and varying as cos 2φ/ρ 2 for the cylindrical annulus). For the spherical shell, the anisotropic susceptibility does not contribute to the external field perturbation, whilst for the cylindrical shell the anisotropic part of the susceptibility makes a four times smaller contribution to the external field than the isotropic susceptibility. Inspection of Eq. (6) shows that simplistically representing the effect of susceptibility anisotropy by setting χ A to zero and making the magnitude of χ I vary with the orientation of the cylindrical shell with respect to the field as sin 2 α, will produce a spatially varying external field whose amplitude varies as sin 4 α, with zero field perturbation inside the shell.

Methods
MRI-based field mapping experiments were carried out at 3T to demonstrate the effects of the anisotropic magnetic susceptibility of shaped structures formed from PGS on the spatial pattern of NMR frequency offsets. Different structures were produced from 25 µm thick PGS (Panasonic EYGS121803). Three-dimensional field maps with 1.5 mm isotropic resolution were acquired using a dual-echo field mapping sequence. For the spherical and cylindrical shell phantoms, the imaging parameters were: TE 1 /TR = 3/40 ms, ΔTE = 1 ms, flip angle = 18°. For the PGS stacks, the imaging parameters were: TE 1 /TR = 3/7.4 ms, ΔTE = 2 ms, flip angle = 8°. All experiments were carried out at room temperature, which was in the range of 18-21 °C. Measuring the external field perturbation due to a small PGS stack. To evaluate the field perturbation produced by a small volume of material of anisotropic susceptibility, disks of nominally 25-μm-thick PGS with a diameter of 4.95 mm were stuck together using a small amount of cyanoacrylate glue to form a cylindrical stack of ~5 mm height, using 65.96 ± 0.03 mg of pre-cut disks with a mean mass of 0.984 ± 0.006 mg per disk. The mean thickness of the PGS disks was measured to be 29 ± 0.5 μm and the volume, ΔV, of PGS in the stack was calculated to be (3.74 ± 0.22) × 10 −8 m 3 . The stack was embedded in an agar-filled, perspex sphere which had a diameter of 180 mm. The agar phantom was constructed from 2% (by total weight) agar (Sigma Aldrich) in distilled water and 0.5% NaCl. The values of isotropic magnetic susceptibility which we determined here and in the other experiments described below (which were also made on PGS structures embedded in agar) are therefore measured relative to the susceptibility of agar, which is around −9 ppm. Field maps were generated with the stack oriented at 5 different angles (Θ) to the field. The field variation in a spherical shell centered on the stack with an internal radius of r i = 10 mm and external radius r o = 15 mm was analyzed in each of these maps to identify the amplitudes of the field contributions due to M x and M z . From Eq. (4), the z-component of the field perturbation (in Hz) may be expressed as , with 42 576 MHzT 1 γ = .
− . To find the values of A x and A z at each stack orientation, the voxels in the spherical shell were first binned based on their θ values (using 10 bins each of width π/10) and the field at each voxel was multiplied by the value of r 3 to eliminate the dependence on radial co-ordinate. The average value of the scaled field offset (〈ΔB z r 3 〉), in each bin was then plotted against (3 cos 2 θ − 1) and linear regression used to estimate the value of A z . Averaging the value of the scaled field multiplied by the value of cos φ for each voxel in a bin and scaling by the average value of cos 2 φ, and fitting the resulting values of B r cos / cos , against 3/2 sin 2θ yielded an estimate of A x . Fitting the variation of these amplitudes with stack orientation to 3/2 sin 2Θ or (1 − 3/2 sin 2 Θ) allowed evaluation of the magnitudes of χ I and χ A using Eqs (3) and (4).
Measuring the field perturbation due to spherical PGS shells. The field perturbation due to spherical shells with anisotropic magnetic susceptibility was assessed by using five thin-walled plastic spheres with radii, R, of 5, 10, 12.5, 19 and 25 mm, which were covered with a layer of nominally 25-µm-thick PGS, filled with water, and then set in a 180-mm-diameter spherical agar phantom. Each spherical shell was formed from 10 identical pieces of PGS, with small gaps left between the pieces to reduce the RF shielding of the inside of the sphere by the electrically conducting graphite. A field map of the phantom was generated and the field variation was measured inside each sphere and over a spherical shell outside each sphere with inner/outer radii = 10.5/18, 18/25.5, 22.5/30, 31.5/39 and 40.5/48 mm for the 5, 10, 12.5, 19 and 25 mm radius spheres, respectively. The difference of the mean field over the internal ROI relative to that in the external ROI was calculated for each sphere. A linear fit of these ΔB z values with respect to 1/R was then calculated, reflecting the field variation expected from Eq. (5). In addition, the external field perturbation was analyzed within a spherical shell outside each sphere with inner/outer radii = 7.5/16.5, 12/18, 15/19.5, 24/28.5 and 30/33 mm for the 5, 10, 12.5, 19 and 25 mm radius spheres, respectively. These radii were chosen to avoid artifacts due to small air bubbles trapped on the surface of the shell and to ensure adequate sampling of the field surrounding each sphere. The field map was additionally filtered using SHARP filtering 26 . This technique selectively removes background fields originating from sources outside of a region of interest, which may have a significant effect on these measurements due to the relatively lower amplitude and larger length-scale spatial variation of the external fields. The field at each point within each annulus was multiplied by the cube of the radius (r 3 ) and then its variation was fitted against 3 cos 2 θ − 1, reflecting the variation expected from Eq. (5). The amplitude of this variation was then used to estimate χ I .
Measuring the internal and external field perturbations due to cylindrical PGS shells. To measure the field perturbation due to cylindrical shells with anisotropic susceptibility, a phantom was constructed by covering glass tubes with three different external diameters (5, 10 and 15 mm) with a layer of nominally 25-µm-thick PGS, producing hollow cylinder structures. A 7.5 cm length of the tube was covered with PGS, but narrow slots running lengthwise down each cylinder were left uncovered to reduce the RF shielding of the inside of the tubes by the graphite coating. The tubes were filled with agar and embedded in a 180-mm-diameter spherical agar phantom. Field maps were generated with the tubes oriented at angles (α) of 1°, 17°, 33°, 49°, 59°, 73° and 89° with respect to the magnetic field. The field perturbations were measured inside each cylinder and over a cylindrical annulus surrounding each cylinder (length = 15 mm, inner/outer radii ≈7.5/15, 6.0/13.5, and 4.5/12 mm for the 15, 10 and 5 mm diameter cylinders, respectively) at each angle to the field. To minimize errors due to other sources of field variation, mean internal field shifts were calculated relative to the average field in the surrounding cylindrical annulus at each angle and the results were compared to the behavior predicted by Eq. (6). In addition, the external field perturbation was analyzed within a cylindrical annulus of inner radius 7.5 mm and outer radius 15 mm around each tube. The field map was additionally filtered using SHARP filtering 26 . The field at each point within the annulus was scaled by multiplication by the square of the radius (ρ 2 ) and then its variation was fitted against cos 2φ, reflecting the variation expected from Eq. (6). This value was then plotted as a function of sin 2 α, where α is the angle between the axis of the tube and the applied B 0 -field.

Results
Measuring the external field perturbation due to a small PGS stack. Figure 2A shows field maps measured in the x-z plane with the PGS stack at three different orientations to the B 0 -field with the normal to the stack also always lying in the x-z plane. Each image shows a 48 × 48 mm 2 region cutting through the center of the PGS stack. Simulated field maps are also shown for the cases where the full effect of the anisotropic susceptibility is accounted for (Fig. 2B) and where only the effect of the induced magnetization that is parallel to the magnetic field is considered (Fig. 2C) voxel. A field perturbation following a standard dipolar field pattern (3 cos 2 θ − 1)/r 3 , aligned with B 0 , is evident in the experimental measurements made when the normal to the stack is nearly parallel (2°) or perpendicular (90°) to the B 0 -field, with the amplitude being greater in the parallel case. At the intermediate angle (64°) the field contribution from M x , which varies as sin 2θ cosφ, is evident from the rotation of the experimentally measured dipolar field pattern. The rotated field pattern is seen in the simulations where the full effect of the anisotropic susceptibility is accounted for (Fig. 2B), but when only the effect of the magnetization that is parallel to the applied field (M z ) is considered, the field variation corresponds to a standard dipolar field with an orientation-dependent amplitude (Fig. 2C). The data shown in Fig. 2B was generated using Eq. (4), while the data shown in Fig. 2C were modelled using only the second of the two summed terms in Eq. (4), which represents the contribution of z-directed magnetization only. Figure 3 shows the variation of B r cos / cos measured in a spherical shell around the PGS stack with 3/2 sin 2θ ( Fig. 2A), along with the variation of the value of 〈ΔB z r 3 〉 with 3 cos 2 θ − 1 (Fig. 3B), for five different orientations of the stack with respect to the B 0 -field. The linearity of these plots indicates that the spatially-varying fields produced by the PGS follow the form predicted by Eq. (4). The calculated variation of A x scaled by γB 0 ΔV/4π = 0.38 Hzm 3 with 3/2 sin Θ cos Θ is shown in Fig. 4A. The slope of this plot gave a value of χ A = − 257 ± 21 ppm. Figure 4B shows the variation of A z scaled by γB 0 ΔV/4π with 1 − 3/2 sin 2 Θ. The slope of this plot gave a value of χ A = −234 ± 18 ppm, while the intercept gave a value of χ I = −121 ± 10 ppm.
Measuring the field perturbation due to spherical PGS shells. Figure 5A shows an axial cross-section through the field map spanning the centers of the five spheres. It shows a uniform negative field offset inside each sphere, with a magnitude that decreases with increasing radius, as predicted by Eq. (5). Figure 5B shows the linear variation of the mean internal field shift with 1/R. Linear regression yields a fit with R 2 > 0.99, indicating that inside the shells ΔB ∝ R −1 , as predicted by Eq. (5). The slope of the plot is γ χ , yielding a value of χ A = − 212 ± 16 ppm. Inspection of Fig. 5A shows that the field is relatively homogeneous outside the spheres, indicating that the external dipolar field perturbation in an axial plane cutting approximately through the center of a sphere is not as large as the internal field offset. However, inset Fig. 5C shows the variation of the field in a ρ-z plane cutting through the 10-mm-diameter sphere. The dipolar field variation due to the isotropic part of the PGS susceptibility is more evident in this field map, which uses a different scaling to Fig. 5A. From these data, we estimate a value of χ I = −111 ± 39 ppm.
Measuring the internal and external field perturbations due to cylindrical PGS shells. Figure 6A shows a 46.5 × 149 mm 2 section of the field map acquired with the three PGS covered tubes oriented perpendicular to B 0 . A negative frequency offset, whose magnitude increases with decreasing shell radius is evident in the map, as predicted by Eq. (6). The field perturbation outside the cylinders, which varies as φ ρ cos2 2 is also evident in this map. Figure 6B shows plots of the average field offset inside the three different cylinders as a function of sin 2 α, where α characterizes the orientation of the cylinder with respect to the applied field. These plots show that the frequency offset scales linearly with sin 2 α with a constant of proportionality which increases in magnitude as R  Figure 6C shows the gradient of the slopes shown in Fig. 6B as a function of the inverse radius of the cylinders. Linear regression gave a good fit to the data (R 2 > 0.99), confirming the linear relationship and yielded a value of χ A = −221 ± 11 ppm. Figure 7 shows variation of the amplitude of the fit of the external field variation to cos 2φ/ρ 2 with sin 2 α for the three PGS covered tubes. As expected from Eq. (6) the amplitude varies linearly with α sin 2 (R 2 > 0.99) and increases with the radius of the tube. An F-test showed no evidence of any significant improvement in the model with the addition of a sin 4 α term which would be expected to occur if an orientation-dependent value for χ I were used in Eq. (6) (P = 0.9018/0.9998/0.1352 for the 5/10/15 mm cylinders). Based on the slopes of the plots in Fig. 7 and the value of χ A calculated from the internal frequency offsets, we estimated a value of χ I = −131 ± 31 ppm.

Discussion
In this work, we exploited MRI-based field mapping in conjunction with the flexible nature and strong susceptibility anisotropy of pyrolytic graphite sheet to visualize directly the magnetic field perturbations generated by simply-shaped structures composed of anisotropic material when they are exposed to a uniform magnetic field. In all cases the measured field perturbations were in good agreement with analytic expressions that are valid in the regime where the magnitude of the volume susceptibility is much less than unity (which is the case for most paramagnetic and diamagnetic materials at ambient temperatures).
The field perturbation due to an elemental volume of anisotropic material arranged at different angles to the field was experimentally characterized and shown to be different from that produced by purely isotropic material. As shown in Eq. (3) this behavior results from the generation of a magnetization component that is perpendicular to the applied magnetic field, as well as from the orientation dependent modulation of the amplitude of the component of magnetization that is parallel to the field. The perpendicular component of the magnetization is largest in amplitude when the principal component of the cylindrically-symmetric susceptibility tensor is oriented at an angle of 45° or 135° with respect to the magnetic field.
The field perturbations due to spherical and cylindrical shells of anisotropic material in which the principal axis of the susceptibility tensor is oriented normal to the shell surface were also experimentally investigated and shown to display behavior that is very different from that produced by similar structures formed from magnetically isotropic material: the most significant difference being the generation of spatially uniform field offsets inside the shells. The experimental data is consistent with the theoretically predicted field patterns, shown in Fig. 1. In the case of the cylinder, the internal field offset is orientation dependent and largest when the cylinder is perpendicular to the applied field. Even at this orientation, the offset produced by the cylindrical shell is smaller by a factor of 3/4 than that produced by a spherical shell. For both types of structure, the internal frequency offset is proportional to the anisotropic susceptibility, χ A , and to the ratio of the shell thickness to its radius. Detailed analysis (Supplementary Material) shows that the internal frequency offset scales as the logarithm of the ratio of the outer and inner diameters of the shell when the shell's thickness is not much smaller than its radius. In the case of the spherical shell, the anisotropic susceptibility does not contribute to the external field perturbation, while for the cylindrical shell the anisotropic susceptibility contributes by a factor of 4 less than the isotropic susceptibility to the spatial varying external field. This interesting behavior also results from the generation of a magnetization component that is perpendicular to the magnetic field and whose amplitude varies coherently within the shell. In the case of the cylindrical shell, this produces a component of the magnetization distribution which follows the pattern of variation found in a Halbach cylinder 27 , (cos sin ), with Λ = H 0 χ A /2. Similarly, the magnetization distribution in the spherical shell due to the anisotropic part of the susceptibility tensor produces  Hz) showing the field variation due to five plastic spheres (radii, R = 5, 10, 12.5, 19 and 25 mm) coated in nominally 25-μm-thick PGS and embedded in a 180 mm diameter spherical agar phantom. A uniform field offset, which increases with decreasing sphere radius, can be seen inside each shell, with little field perturbation outside each sphere. (B) Plot showing the variation of the mean internal field shift due to the spherical shells 〈ΔB z 〉 as a function of 1/R. Error bars show the standard deviation of the field measured inside each sphere. (C) Field variation in field in a ρ-z plane cutting through the 10-mm-diameter sphere. Note that the frequency range shown here is a factor of 6 smaller than that shown in A (with values clipped to ±25 Hz to aid visualization of the external dipolar field pattern). with Λ = H 0 χ A , as is found in a Halbach sphere 27 .
The values of the anisotropic susceptibility, χ A , of the 25-μm-thick PGS measured from the two different forms of spatially varying field found around the slab (−257 ± 21 ppm and −234 ± 18 ppm) were consistent with one another and with those measured from the internal field offset in the spherical (−212 ± 16 ppm) and cylindrical (−221 ± 11 ppm) shells. The χ A -values measured from the sphere and cylinder would be expected to be slight underestimates because of the small gaps which were left between the PGS strips forming the shells to reduce RF screening effects. This meant that there was less PGS material present than was assumed in deriving expressions for the field perturbation based on continuous shells. Using the value of χ I = −121 ± 10 ppm estimated from the measurement of the dipolar field around the slab and the average of the above χ A -values, yields χ ⊥ = −450 ± 15 ppm and χ || = −5 ± 15 ppm. The magnetic susceptibility of commercially supplied PGS has not previously been reported as far as we are aware, but these values are consistent with previous literature on the susceptibility of pyrolytic graphite 25 . The values of χ I estimated from analysis of the external fields produced by the spherical (111 ± 39 ppm) and cylindrical shells (131 ± 31 ppm) are in general agreement with the value calculated from the measurements on the PGS stack.
The measurements reported here provide some supporting evidence for the models that have previously been used to explain the effects of anisotropic magnetic susceptibility in biological samples. The largest such effects are seen in the white matter of the brain in which myelinated nerves are the source of susceptibility anisotropy. The myelin sheaths of these nerves are formed from lipid bi-layers in which the long lipid chains are radially oriented, such that the principal axis of the susceptibility tensor in the myelin sheath is radially directed and the value of χ A is negative 6,7 . This is similar to the geometry of the cylindrical PGS shells investigated here, and as a consequence of the larger size and much greater anisotropy of the PGS structures, we were able to directly visualize the uniform negative field offset inside the cylinder (Fig. 6A) and experimentally to confirm the predicted dependence of the field perturbations on the shell geometry and orientation (Figs 6 and 7). This dependence is important since it underpins the interpretation of the time dependence of the gradient echo signal generated from white matter in vivo 21,28 . In contrast to the PGS-based model system used here, some NMR signal is generated by the cylindrical shell in myelinated nerves, but the low concentration and short T 2 * -relaxation time of the myelin water means that this signal does not contribute significantly to images acquired at standard echo times (>10 ms) in gradient echo imaging. It is consequently the negative intra-axonal frequency that generally makes the dominant contribution to the negative average NMR signal frequency (and related phase effects) seen in white matter regions where the nerve fibers are oriented perpendicular to the applied field 1,7 . This effect also underlies the increased range of signal frequencies and consequent reduction of T 2 * found in these white matter regions compared with those in which the nerve fibers are parallel to the field 29 . More generally, Figs 5A and 6A make it evident that the field variation and associated NMR signal variation produced by microstructural tissue features formed from components of anisotropic magnetic susceptibility will not usually be well characterized by use of a model that represents each image voxel as a homogeneous piece of material of isotropic or anisotropic susceptibility. Since this is the modelling approach that is used in QSM and STI 29 , there are potential problems with accurate characterization of the susceptibility of white matter tracts and adjacent tissues using these techniques 30 , so the results which they produce must be interpreted with care.
The measurements of the field perturbations generated by the PGS slabs make evident the importance of accounting for the effect of the component of the induced magnetization that is perpendicular to the applied field when considering the effect of anisotropic susceptibility -the maximum field perturbation that this component of the magnetization generates at a given radius for fixed χ A is more than half of the maximum field generated by the parallel magnetization, and the amplitude of field variation due to the perpendicular magnetization is larger than that due to the parallel magnetization for values of Θ between ~34 and 71° (and ~109 and 146°). Approaches that analyze the effect of the anisotropic susceptibility in terms of an apparent isotropic susceptibility whose amplitude varies with the orientation of the material with respect to the magnetic field 5,15 do not therefore properly characterize the field perturbations and may yield erroneous results. For example, applying this approach to a cylindrical shell would yield no internal field offset and a contribution to the cos 2φ/ρ 2 external field variation whose amplitude varies with sin 4 α 14 -no such variation is of course seen experimentally.
The maps which show the uniform field offsets inside spherical shells and the scaling of this offset with the inverse of the shell radius provide direct experimental visualization of the model previously used by Lounila et al. 10 to explain size-dependent spectral shifts in the 13 C NMR frequency measured from compounds confined within naturally occurring lipoproteins. These shifts result from the anisotropic magnetic susceptibility of the lipid chains that form the lipoprotein shell in such a way that the long axis of the chains, and thus the principal axis of the susceptibility tensor, is always radially oriented in the shell: a geometry that is replicated by the PGS-covered spherical shells. It is worth noting however that the anisotropic susceptibility of lipid membranes is estimated to be around −0.2 ppm 7,10 , which is more than 1000 times less than the value in PGS. The high magnetic anisotropy and relatively low density of PGS means that large field perturbations can be generated with small amounts of material. The work described here highlights new opportunities for improved local passive shimming in MRI 31 , where the additional degrees of freedom provided by the possibility to change the orientation of the induced magnetization could allow better cancellation of field inhomogeneity by appropriate placement of PGS material around the outside of the body. In addition, microscopic PGS shells could potentially be exploited in producing "multi-spectral" contrast agents for MRI. Recent work by Zabow et al. 32,33 has shown how nano-fabricated structures formed from pairs of nickel or iron disks can be used as contrast agents that encode distinguishable spectral signatures. This is possible because of the geometry-sensitive field shift produced between the paramagnetic disks. The capability to produce a geometry-sensitive field shift inside a spherical shell formed from PGS could offer advantages over planar metallic structures, although the RF screening effects of PGS would need to be addressed.