PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
 
Magn Reson Med. Author manuscript; available in PMC 2010 June 23.
Published in final edited form as:
PMCID: PMC2890274
NIHMSID: NIHMS209985

The Design and Test of a New Volume Coil for High Field Imaging

Abstract

A major problem in the development of high field (>100 MHz) large volume (>6000 cm3) MR coils is the interaction of the coil with the subject as well as the radiation loss to the environment. To reduce subject perturbation of the coil resonance modes, a volume coil that uses an array of freely rotating resonant elements radially mounted between two concentric cylinders was designed for operation at 170 MHz. Substantial electromagnetic energy is stored in the resonant elements outside the sample region without compromising the efficiency of the overall coil. This stored energy reduces the effect of the subject on the circuit and maintains a high Q, facilitating the tuning and matching of the coil. The unloaded O of the coil is 680; when loaded with a head, it was 129. The ratio of 5.3 of the unloaded to loaded Q supports the notion that the efficiency of the coil was maintained in comparison with previous designs. The power requirement and signal-to-noise performance are significantly improved. The coil is tuned by a mechanism that imparts the same degree of rotation on all of the elements simultaneously, varying their degree of mutual coupling and preserving the overall coil symmetry. A thin radiofrequency shield is an integral part of the coil to reduce the radiation effect, which is a significant loss mechanism at high fields. MR images were collected at 4T using this coil design with high sensitivity and B1 homogeneity.

Keywords: 4T, head, imaging, efficiency

INTRODUCTION

Magnetic resonance studies at high magnetic fields are fundamentally advantageous because of the inherent high sensitivity (1, 2), greater spectroscopic resolution (35), increased susceptibility effects for contrast with both external and internal agents (deoxyhemoglobin) (68), and “black blood” flow quantification schemes that improve with the extended T1 at high fields (9). To realize the signal to noise advantage of high magnetic fields, the front end of the system (including the radiofrequency coil and receiver) must be optimized. The challenging issue for high field (>100 MHz) large volume coils (>6000 cm3) is the variable coil-subject interaction that significantly changes the field distribution and resonance frequency of the coil. In addition, the radiation loss from the coil to the environment increases drastically with field strength (10, 11). Due to the short wavelength at high Larmor frequencies (20 cm for protons at 4T), far field effects cannot be ignored. A useful theoretical tool for coil optimization is the reciprocity relation between coil efficiency and signal-to-noise (SNR) ratio (1, 12). This relation must be generalized beyond quasistatic approximation and simple coil-sample geometry to apply to higher field strength, as discussed in the Appendix.

The conventional “birdcage” design that has proven to be successful at lower field strength (13, 14) encounters problems at higher fields. At 4 Tesla, the radiation loss of a birdcage head coil without a RF shield degrades the Q value to ~20. After an RF shield is added, the Q value is increased to ~100. However, when loaded with a human head, Q again drops to ~35. This is accompanied by significant perturbations of the resonance modes of the coil indicated by the frequency shift (about 1 MHz for our conventional birdcage coil with a shield). This is probably due to the fact that the human head is placed at the region of the coil where most of the field energy is stored. These specific problems of high field large volume coils call for a coil design that reduces the subject perturbation effects and the radiation loss while maintaining coil efficiency.

METHODS

To overcome the loss factors, subject perturbation effects, and loaded Q degradation, a quadrature-driven volume head coil was designed and constructed with the structure of an array of individually tuned, inductively coupled resonant elements in a relatively enclosed cage to form a low-loss resonator (Fig. 1). The electromagnetic energy stored in the elements outside the loading region of the coil reduces the perturbation of the subject on the resonance modes of the coil while maintaining the efficiency. This design is named “the free element design” (FE design).

FIG. 1
A picture of the prototype head coil.

The design enables convenient coil tuning with respect to each patient to achieve optimal performance, without disturbing its symmetry. The initial tuning of the coil is greatly simplified to the tuning of individual elements. The design also enables asymmetric adjustments to compensate for the specific part of the body being imaged.

The coil designed and tested has two coaxial cylinders (Figs. 1, ,2).2). Both cylinders are made of 0.635-cm-thick fiberglass. The outer diameter of the exterior cylinder is 40.64 cm, the clear bore of the inner cylinder (the imaging region of the coil) is 27.94 cm, and the length of the coil is 36.83 cm. The inner surface of the outer cylinder is lined with a contiguous piece of thin conductor (in our case, it is 0.0025-cm-thick copper sheet) as the RF shield. Terminal end rings are placed in the interspace between the two concentric cylinders at both ends of the coil.

FIG. 2
A schematic of a section of the coil and the resonant element. For symmetric tuning, all elements are rotated simultaneously by the angle.

Each resonant element is a rectangle of 0.635-cm-diameter copper tube spaced with series capacitors. One of the capacitors is a variable capacitor that tunes the resonant element (Fig. 2). The dimension of the copper tube rectangle is 4.44 cm by 20 cm. Each element is suspended in the space between the two concentric cylinders by attaching its ends to axles mounted on the terminal end rings filling the interspace. The axles are 1.27-cm-diameter lexan rods. One of the end rings contains a mechanism that permits uniform rotation of all of the elements around their axes (Fig. 2). This rotation uniformly changes the mutual inductances between the resonant elements and thereby changes the resonance frequency of the coil.

Alternatively, the coil can be tuned asymmetrically by turning each element individually around its axis or by adjusting the variable capacitor on the element. This asymmetric adjustment may be used to generally compensate for the asymmetric perturbation introduced by the part of the body being imaged.

Denoting the self inductance and capacitance of each resonant element as L and C and the mutual inductance between two elements n units apart as Mn, the resonance modes of the coil are solved from the following equation,

iωLIk+IkiωC+jkiωMkjIj=0,
[1]

where Ik is the current in the kth element. The resonance modes of interest are the two degenerate modes with the current distribution

Ik=I0cos(2πkN+φ),
[2]

where N is the total number of elements (in our case, it is 16) and the phase shift [var phi] is 0 and π/2 for the two modes, respectively. These two modes have a uniform magnetic field distribution in the loading region. The resonance frequency of these two modes is solved from Eq. [1] as

ω={CL[1+j0MjLcos(2πj/N)]}1/2=ω0[1+j0MjLcos(2πj/N)]1/2,
[3]

where ω0 is the resonance frequency of each individual element.

To determine the initial resonance frequency of each individual element such that the desired frequency ωL falls into the tuning range of the coil, the following procedure is used:

Two resonant elements are tuned to an identical frequency ωt not far away from the desired frequency ωL. The two elements are placed in the coil at n units apart. Due to the mutual inductance of the two elements, two resonance modes exist, with the following frequencies:

ω±=ωt[1±Mn/L]1/2.
[4]

These two frequencies are measured for each n, and the ratio Mn/L is then calculated from Eq. [4]. Knowing the ratios Mn/L and assuming they do not change dramatically with frequency, they are used in Eq. [3] to predict the resonance frequency ω0. This method has been found to be accurate within 2%.

The coil is driven in quadrature at two elements 90° apart. The actual driving circuit is shown in Fig. 3. For each driven element, a coaxial cable is introduced into the coil with its outer conductor grounded to the shield and its inner conductor connected to a tap point on the element via a flexible copper band. This enables the driven element to pivot freely just like all of the other elements. A series capacitor is used for matching the high impedance of the tap point to 50 ohms.

FIG. 3
A sketch of the driving circuit to one driven element.

RESULTS

The prototype head coil operates at 170.74 MHz (4T protons). The tuning range is ±6 MHz for symmetric tuning with the uniform change of mutual inductances between elements, as mentioned above. The Q values in both loaded and unloaded conditions were measured with two 2-cm-diameter shielded loops. The transmission coefficient between the two loops was measured with an HP4195A network analyzer (Fig. 4). The Q values from Fig. 4 are 680 without loading and 129 when loaded with a human head. Also indicated by Fig. 4 is that the resonance frequency increased by ~120 kHz when loaded, indicating inductive loading of the coil. The high loaded Q signifies a significant reduction in the subject perturbation of the whole coil, whereas the ratio of 5.2 between unloaded Q and loaded Q shows that the efficiency of the coil is well maintained (15).

FIG. 4
Transmission coefficient between two shielded probe loops in the coil. The Q of unloaded and loaded conditions were calculated from this measurement.

The coil can be conveniently tuned and matched with the symmetric tuning mechanism and the matching capacitors. Figure 5 shows the reflection coefficient from one of the two driving ports in the three states: tuned and matched to empty load, loaded with a head, and retuned and matched with the head. The 3dB width of the empty load condition was 345 KHz. When retuned and matched to the head, the 3dB width broadened to 1.44 MHz. The relatively high value of loaded Q made the resonance width narrow enough to accurately tune the coil. When tuned and matched with a head, the isolation between the two quadrature channels was between 25 and 30 dB.

FIG. 5
Reflection coefficient at one drive port in loaded and unloaded conditions.

As a test of the B1 homogeneity of the coil, a spin echo image was acquired on a phantom made of a matrix of test tubes containing saline water (Fig. 6). An array of test tubes was used to avoid any long current paths, which give rise to sample resonance effects at 170 MHz (4). This segmented sample more realistically reflected the field distribution of the coil when compared to a solid sample. The signal intensity of the test tube array across the diameter and along the axis of the coil is plotted in Fig. 7. In the transaxial plane, the variation over 20 cm distance is ~10%. Along the axial direction, the signal stayed within 10% over 15 cm, which is 3/4 of the resonant element’s length.

FIG. 6
A spin echo transaxial phantom image taken as a measure of the B1 homogeneity of the coil. TE = 21 ms; TR = 1 s, voxel = 1 × 1 × 4-mm.
FIG. 7
Signal intensity in the transaxial and axial directions of the coil; (a) is of the transaxial plane, (b) is of the axial direction.

To evaluate the B1 homogeneity and SNR of the coil in a human subject, a single acquisition proton density spin echo image was collected on a normal volunteer (Fig. 8). The voxel size was 1×1×4 mm, TE = 26 ms, TR = 3 s. The SNR averaged 80:1 in a single white matter pixel. The improved efficiency of the coil was demonstrated by the modest 229 watts peak power required for a 90° 3-ms sinc pulse (sincc = 2) used in this study. As a comparison, a conventional birdcage coil with an RF shield similar to the FE coil was used for a spin echo image with the same imaging parameters. The RF shield, not usually used on birdcage coils at low field strengths, was necessary at 170 MHz to bring the loaded Q to a reasonable level for tuning purpose. The unloaded Q of the coil was 96; it was 35 when loaded with a head. The peak power required for a 90° 3-ms sinc pulse (sincc = 2) was 620 watts.

FIG. 8
Spin echo proton density image of a head at 4T. NEX = 1, TR = 3 s, TE = 26 ms, voxel = 1 × 1 × 4 mm.

The homogeneity of the coil was perturbed slightly by the subject, requiring a small amount of asymmetric tuning of the individual elements located near the back of the head. The actual asymmetric tuning was done empirically by adjusting the resonance frequencies of the four elements at the back of the head. Their resonance frequencies were tuned down by ~0.6% to compensate for the decreased B1 there.

DISCUSSION

The cylindrical arrangement of resonant elements is an effective solution to the problems encountered with large volume NMR coils at high fields. This arrangement permitted accurate planning of the coil’s tuning and matching characteristics, reduced the net perturbation of the subject on the coil, enabled precise tuning and matching with each patient, as well as reduction in dielectric and radiation losses of the coil. These factors resulted in a very efficient MR volume coil, which produced improved MR images at 4T. A simple mathematical model was shown to be adequate for planning and construction of the coil within reasonable range of the target values. This is an important aspect because it reduces the empirical iterations usually required. For high field coil optimization problems, a generalized form of the reciprocity relation between coil efficiency and sensitivity was derived. It is applicable to all frequencies and coil-subject geometries.

The net perturbation of the subject on the coil was reduced significantly, as illustrated by the high loaded Q value. However, the efficiency and sensitivity of the coil is well maintained. The B1 distribution of the coil without sample dielectric resonance was demonstrated to be quite uniform with a segmented phantom. The B1 homogeneity with a head was slightly perturbed and could be easily corrected with a small amount of asymmetric tuning of a few resonant elements located near the back of the head.

The radiation loss of the coil is significantly reduced by the RF shield as shown in the high unloaded Q. The self-shielding nature of each resonant element may also contribute to the decrease in radiation loss (16).

Because the design consists of individual tuned elements, multiple tuning can be easily realized by making each element multiple tuned. Because each resonant element is relatively small, many multiple tuning schemes can be applied to it without losing much of its efficiency. This, in turn, also guarantees the efficiency of the whole coil in the multiple tuned configuration.

Whereas the prototype head coil has cylindrical symmetry, the general design is not confined to this form. In the case of thoracic or abdominal imaging, for example, a coil of elliptical cross section with the elements asymmetrically arranged is conceivably preferential.

APPENDIX

The General Reciprocity of an NMR Transceiver

Relations between transceiver sensitivity and efficiency have been derived for simple coil geometries, often under quasistatic approximations (1, 12). At 4T, the Larmor frequency of protons is 170 MHz, high enough so that the wavelength of 20 cm in water approaches the scale of the subject in human imaging studies. Theoretical description beyond quasistatic approximation and simple coil-sample geometry is needed. It is shown here theoretically that a general reciprocal relation exists between the sensitivity and the efficiency of a transceiver (CGS units):

SNR=18sin(α)eτnVoxN1/22ω3B01P1G1/2(kT)3/2(Δf)1/2,
[A1]

where SNR is the signal to noise level of a voxel in the image, α is the flip angle of the experiment, e−τ represents the spin relaxation factor of the nuclear magnetization, n is the spin density, Vox is the volume of the voxel, N is the number of scans in the image acquisition, P1G is the power needed to generate a circularly polarized RF magnetic field of 1 Gauss magnitude at the voxel, Δf is the frequency bandwidth of each voxel, ω is the Larmor frequency, B0 is the main static magnetic field, T is the temperature, and [variant Planck's over 2pi] is the Planck’s constant.

In the following discussion, all electromagnetic fields and moments are oscillating at frequency ω, and the complex notation a(t) = aeiωt is used, a being a constant complex vector. The factor eiωt will be omitted unless specified with “(t)”, and the following notation is also used:

1TTTa(t)·b(t)dt=a·b+b·a2=ab.
[A2]

Apparently

ab=ab.
[A3]

Consider the simple case of two general electromagnetic “current” elements (electric dipole, magnetic dipole, etc.) I1 and I2 at coordinates r1 and r2 in free space, the conjugate field of I2 at r2 generated by I1 is represented by the retarded free space. Green function operator G21,

v21(t,r2)=G21(t,r2,t,r1)I1(t)dt,
[A4]

vice versa,

v12(t,r1)=G12(t,r1,t,r2)I2(t)dt.
[A5]

From the explicit form of the free space retarded Green functions, it is easy to see that

I1G12I2=I2G21I1.
[A6]

Now consider the general case of two current elements Is and Id at rs and rd. The two current elements are emersed in a linear dissipative environment, which can be represented by many discreet point “elements,” each having a linear response to the local electromagnetic field:

I1(t)=εiv1(t,ri),
[A7]

where Ii is the general current of the element i in response to its conjugate field Vi at its locale ri, εi, is the response operator. If the linear environment is at or near thermal equilibrium, Onsager’s relation means that for any oscillating fields V1 (ri) and V2 (ri), the following is satisfied by the response operator,

εiv1(ri)v2(ri)=v1(ri)εiv2(ri).
[A8]

Taking into account the effect of the medium, the conjugate field of Id generated by Is at rd is

vds=vds+ivdi=GdsIs+iGdiIi,
[A9]

where lower case vds is the field at rd generated by Is in the absence of the medium, vdi is the field contribution of the current Ii of medium “element” i. Because

Ii=εvi=ε(vis+jivij)=ε(vis+jiGijIj),
[A10]

substituting Eq. [A10] into Eq. [A9] iteratively gives

Vds=vds=ivdi=GdsIs+iGdiIi=GdsIs+iGdiεiGisIs+iGdiεijiGijIj=GdsIs+iGdiεiGisIs+iGdiεiijGijεjGjsIs+iGdiεiijGijεjjkGjkεkGksIs+=FdsIs,
[A11]

where Fds denotes the Green function in the presence of the medium. With the Eq. [A3], [A6], and [A8],

IaGabεbVb=GbaIa(εbVb)=(GbaIa)εbVb=εbGbaIaVb.
[A12]

Thus

IdGdiεiGijεjGmsIs=εjGidId(GijεjGmsIs)=(εiGidId)GijεjGmsIs=εjGjiεiGidId(GmsIs)=εjGjiεiGidId(GmsIs)=(εjGjiεiGidId)(GmsIs)=GsmεjGjiεiGidIdIs=IsGsmεmεjGjiεiGidId.
[A13]

Equations [A11] and [A13] together lead to the general reciprocal relation of the Green function in a linear environment at or near thermal equilibrium:

IdFdsIs=IsFadId.
[A14]

To apply the general reciprocity to the discussion of sensitivity versus efficiency, the two specific current elements are chosen to represent the nuclear magnetic moment at rd in the subject, and the RF power input/signal output port of the transceiver. The element at rd is thus the circularly polarized nuclear magnetic dipole moment

Id=m=m0(ex+iey),
[A15]

where m0 is the amplitude of the magnetic moment, assuming the main static magnetic field is in the z direction. The other element is simply a current source Is at the I/O port of the transceiver. The conjugate field of Id generated by Is is

Vds=×E1(t)=dB1(t)dt=iωB1
[A16]

The conjugate field of Is generated by Id is simply the electric voltage across the current source. If Vs is the voltage across Is from the nuclear magnetic moment m* = m0 (exiey), then Eq. [A14] results in

IsVs=iωB1m.
[A17]

Note ex+iey as the clockwise polarization, Is the magnitude of Is, Vs,cc the amplitude of the voltage generated by the counter-clockwise polarized m*, and B1,c the amplitude of the clockwise component of B1. It is easy to see from Eq. [A17] that

IsVs,cc=2ωB1,ccm0.
[A18]

Because B1,c/Is reflects the efficiency of the transceiver and Vs,cc/m0 reflects the sensitivity of the transceiver, Eq. [A18] is the basis for the eventual relation [A1].

Notice that Eq. [A18] relates the coil’s efficiency to generate a circularly polarized magnetic field to its sensitivity toward the nuclear magnetic moment of the other circular polarization. This change of polarization makes it necessary to distinguish linear transceivers from quadrature transceivers, as discussed below.

A linear transceiver generates linearly polarized B1. Because a linear field is the sum of the two circular polarizations of equal magnitude, Eq. [A18] can be written as

IsVs,cc=2ωB1,ccm0,
[A19]

where B1,cc is the counter-clockwise component of B1. If the impedance into the output port of the transceiver (where the element Is lies) is Zm, its resistive part Rm, then the RF power needed to generate B1,cc is

Pt=Is2Rm2.
[A20]

With Eq. [A20], Eq. [A19] becomes

Vs,cc=2ωB1,ccm02Pt/Rm=2ωRmm0P1G,
[A21]

where P1g is the power needed to generate a circularly polarized B1 of 1 Gauss magnitude at the voxel. If the preamplifier’s input impedance is Zp the actual voltage on the preamplifier is

Vp(t)=Vs,cc(t)ZpZp+Zm.
[A22]

The noise voltage Vn from the transceiver is also divided between Zp and Zm,

Vnp(t)=Vn(t)ZpZp+Zm.
[A23]

The SNR is thus

SNR2=Vp2(t)Vnp2(t)=Vs,cc2(t)Vn2(t)=Vs,cc22Vn2(t).
[A24]

If the frequency bandwidth of each voxel in the data acquisition is Δf, then

Vn2(t)=4ΔfRmkT,
[A25]

where k is the Boltzman’s constant. Substitution of Eq. [A21] and Eq. [A25] into Eq. [A24] gives

SNR2=ω2n024ΔfkTP1G.
[A26]

For a voxel of volume Vox and spin density n, the magnitude of the transverse magnetic moment can be expressed as

m0=eτsin(α)M0nVox,
[A27]

where M0 is the average magnetization of each spin, and the other parameters are defined under Eq. [A1]. It has been shown (17] that for magnetic dipoles,

M0=2ω24B0kT.
[A28]

Substituting Eq. [A28] and Eq. [A27] into Eq. [A26], taking into account the factor N1/2 if N scans are performed to acquire the image, the general relation Eq. [A1] is obtained.

A quadrature transceiver is essentially two linear coils combined at a quadrature hybrid (see, e.g., ref. 15). The hybrid produces 90° phase difference between the two linear coils and a 180° phase flip on one of the linear coils between the transmission and the reception. Thus, the transceiver transmits and receives in the same circular polarization, and Eq. [A1] still holds. The difference is that a quadrature transceiver generates only the correct circular polarization when transmitting, whereas the linear transceiver generates an equal amount of both circular polarizations. Thus, everything being equal, a quadrature coil needs only half of the power to produce the same magnitude of B1,cc, which is reflected in the value of P1G in Eq. [A1].

An important noise source that is not included in Eq. [A1] is the coherent noise in the environment. Depending on the shielding of the magnetic room and the general geometry of the coil and the subject, etc., the coherent noise can be picked up by the coil to various degrees. The general reciprocity relation [A1] therefore predicts an SNR that is usually higher than the measured value.

References

1. Hoult DI, Richards RE. The signal-to-noise ratio of the nuclear magnetic resonance experiment. J Magn Resort. 1976;24:71–85. [PubMed]
2. Hoult DI, Lauterbur PC. The sensitivity of the zeugmatographic experiment involving human samples. J Magn Reson. 1979;34:425–433.
3. Barfuss H, Fischer H, Hentschel D, Ladebeck R, Vetter J. Whole-body MR imaging and spectroscopy with 4-T system. Radiology. 1988;169:811–816. [PubMed]
4. Barfuss H, Fischer H, Hentschel D, Ladebeck R, Oppelt A, Wittig R, Duerr W, Oppelt R. In vivo magnetic resonance imaging and spectroscopy of humans with a 4T whole-body magnet. NMR Biomed. 1990;3:31–45. [PubMed]
5. Hardy CJ, Bottomley PA, Roemer PB, Redington RW. Rapid 31P spectroscopy on a 4-T whole-body system. Magn Reson Med. 1988;8:104–109. [PubMed]
6. Kwong KK, Belliveau JW, Chesler DA, Goldberg IE, Weisskoff RM, Poncelet BP, Keddedy DN, Hoppel BE, Cohen MS, Turner R, Cheng HM, Brady TJ, Rosen BR. Dynamic magnetic resonance imaging of human brain activity during primary sensory stimulation. Proc Natl Acad Sci (USA) 1992;89:5675. [PubMed]
7. Turner R, Bihan DL, Moonen CTW, Despres D, Frank J. Echo-Plana Time course MRI of cat brain oxygenation changes. Magn Reson Med. 1991;22:159–166. [PubMed]
8. Turner R, Jezzard P, Wen H, Kwong KK, Bihan DL, Zeffiro T, Balaban RS. Functional mapping of the human visual cortex at 4 and 1.5 tesla using deoxygenation contrast EPI. Magn Reson Med. 1993;29:277–279. [PubMed]
9. Detre JA, Leigh JS, Williams DS, Koretsky AP. Perfusion imaging. Magn Reson Med. 1992;23:37–45. [PubMed]
10. Harpen MD. Radiative losses of a birdcage resonator. Magn Reson Med. 1993;29:713–716. [PubMed]
11. Jackson JD. Classical Electrodynamics. John Wiley & Sons; New York: 1975.
12. Tropp J. An extended reciprocity principle for NMR reception. Proc., SMRM, 12th Annual Meeting; New York. 1993. p. 1320.
13. Hayes CE, Edelstein WA, Schenck JF, Mueller OM, Eash M. An efficient, highly homogeneous radiofrequency coil for whole-body NMR imaging at 1.5 T. J Magn Reson. 1985;63:622–628.
14. Hinshaw WS, Gauss RC. European Patent Application. P1303522.7
15. Chen C, Hoult DI. Biomedical Magnetic Resonance Technology. Adam Hilger; Bristol and New York: 1989.
16. Mansfield P, Chapman B. Active magnetic screening of coils for static and time dependent magnetic field generation in NMR imaging. J Phys E: Sci Instrum. 1986;19:540–545.
17. Abragam A. The Principles of Nuclear Magnetism. Oxford University Press; London: 1961.