Analysis of Uniformity of Magnetic Field Generated by the Two-Pair Coil System

In this paper we use a simple analysis based on properties of the axial field generated by symmetrical multipoles to reveal all possible distributions of two coaxial pairs of circular windings, which result in systems featuring zero octupole and 32 pole magnetic moments (six-order systems). Homogeneity of magnetic field of selected systems is analyzed. It has been found that one of the derived systems generates homogenous magnetic field whose volume is comparable to that yielded by the eight-order system. The influence of the current distribution and the windings placement on the field homogeneity is considered. The table, graphs and equations given in the paper facilitate the choice of the most appropriate design for a given problem. The systems presented may find applications in low field electron paramagnetic resonance imaging, some functional f-MRI (nuclear magnetic resonance imaging) and bioelectromagnetic experiments requiring the access to the working space from all directions.

In many applications it is required that the field be highly homogenous over some specified volume. This is of particular importance in magnetic resonance imaging experiments. The systems used for in vivo medical diagnostic studies most often employ solenoidal superconducting electromagnets that are expensive and in certain applications pose disadvantages associated with the limited access to the region of uniform field. In electron paramagnetic resonance imaging (EPRI) [5] and in some functional nuclear magnetic resonance imaging (MRI) experiments [6], electromagnets generating low fields and/or allowing access to the working space from all directions and not just axial are desirable. A classic example of systems satisfying these conditions are air-core assemblies comprising a number of circular or square windings placed co-axially and distributed so that the leading perturbation terms in the field series expansion are eliminated.
In this paper, we consider the system consisting of two coaxial pairs of circular loops with the same radius. The use of properly distributed windings of the same radius makes the radial access to the uniform field possible and does not impose restriction on the axial access, which may sometime occur in systems based on spherical harmonic expansion [4,10] with outer pair of windings of smaller radius or having a single loop in the mid plane.
We use a simple analysis based on properties of the axial magnetic field with the aim to reveal the possible distributions of windings, that result in systems featuring zero octupole and 32 magnetic moments, i.e., generating the central magnetic field in which the sixth-order term is the first one non-vanishing in the field expansion. The table, formulae, and graphs given in the paper facilitate the choice of design, which is the most suitable for the problem at hand.
The system presented generates extended volume of uniform magnetic field, which can be accessed from all directions. It may be suitable for very-low field MRI and EPRI as well as bioelectromagnetic experiments [7]. The high-field system can be easily shielded by confinement in other with larger radius, which cancels the total dipole moment and results in reduction of the stray field at the expense of a slight decrease of strength of the very homogenous central field.

System Configuration Analysis
Consider a system of two circular current loops of the same radius R, with current I flowing in the same sense in each loop. Let the loops encircle z-axis, be located symmetrically with respect to the origin, and separated by the distance 2d. Using the Biot-Savart law we calculate that the axial field B z is: Let us now focus the attention on how the axial field behaves outside the coil system and in its central region. To show the behavior of the distant axial B z field, we may expand in the Taylor series the function in square brackets of Eq. (1) in inverse powers of z. This yields: It is seen the distant axial field may be viewed as the superposition of magnetic multipoles at the system origin. Due to symmetry of current distribution all coefficients c kÀ1 for an even k will vanish. This means that in the field expansion there will be no 2 k poles with k even, just dipole (k ¼ 1), octupole (k ¼ 3), 32 pole (k ¼ 5) and so on.
The first coefficients are found to be: The first term in the field expansion c 0 z À3 is evidently the axial field due to the dipole at origin (k ¼ 1). The second term c 2 z À5 is the octupole term (k ¼ 3), the third (k ¼ 5) 32 pole term and so on. But we see from Eqs. (2) and (3b) that the octupole term becomes zero if 2d ¼ R (the Helmholtz condition). Now as Purcell [9] pointed out, the field of any symmetrical multipole cannot be zero along the symmetry axis unless it is zero off the axis everywhere as well. Consequently, the entire octupole field, not just B z on the axis, vanishes, if 2d ¼ R. This means that every Helmholtz pair has the zero octupole moment.
By expanding in the Taylor series the function in square brackets of Eq. (1) around the center, we can readily show that the central field generated by the system is: where: Owing to the symmetry of current distribution the axial magnetic field in the central region of the system is an even function of z. But for the Helmholtz condition (2d ¼ R) the c 2 coefficient vanishes, which maximizes the volume and uniformity of the field in the central region of the system. It is possible to make this field still more uniform and to extend its useful volume by using additional pairs of coils, which help to nullify the next coefficient in the field expansion [1].
To eliminate both octupole and 32 pole terms, we consider a system of two coaxial pairs of circular current loops of the same radius R, i.e., a four-coil system with number of Amp-turns ðNIÞ 1 and ðNIÞ 2 , respectively, and current flowing in the same sense in each loop. Let the loops be located symmetrically with respect to the origin of z-axis at distances 2d 1 and 2d 2 , respectively. From the Biot-Savart law the axial field B z generated by the system is: Expanding in the Taylor series the function in square brackets around the center we have: where c ð1Þ kÀ1 and c ð2Þ kÀ1 are the coefficients in the multipole field expansion for respective coil pairs. By analogy to Eqs. (5a, 5b, 5c) we may write: It is seen that both octupole (k ¼ 3) and 32 pole (k ¼ 5) terms in the field expansion will be zero if: The set of Eqs. (9a, 9b) contains three variables d 1 ; d 2 ; ðNIÞ 1 =ðNIÞ 2 È É , one of which can be treated as a parameter. In our calculation the chosen parameter has been the Amp-turn ratio.
The numerical solutions of Eqs. (9a, 9b) may be divided into two families. The first one, which contracts the length of the coils system, and the second one, which extracts the length. Table 1 Coil spacings and design specifications for contracting systems Uniformity of Magnetic Field Generated by the Two-Pair Coil System 609 Table 1 continued  Table 1 lists the Amp-turn ratios and the corresponding coil spacing obtained for the family of solutions contracting the coil length. Using the LS approximation we have found that these solutions may be approximated by:

Numerical Results
¼ 0:25153 þ 0:06065e Àt À 0:00173t À 0:00001t 2 ; ð10aÞ d 2 R ¼ 0:96173 À 0:06781e Àt À 0:00466t þ 0:22810te Àt þ 0:00003t 2 ; ð10bÞ ðNIÞ 2 =ðNIÞ 1 ¼ 2:12 þ 0:02t ð10cÞ for 1 t 82. The obtained solutions are shown graphically in Fig. 1. For the family extracting the length, the solutions may be approximated by: for 1 t 394. The latter solutions are shown graphically in Fig. 2a, b. The derived solutions allow easy construction of the four-coil system featuring zero octupole and 32 pole terms. For the family of solutions contracting the coils length one set of the derived variables corresponds well to that obtained previously by means of the Bessel functions formalism by Lee-Whiting [8]. As it is seen from  Hereafter, we consider in more detail the NI ð Þ 2 = NI ð Þ 1 ¼ 9=4 system, which has the potential of being power supply in series. For this Amp-turn ratio the windings are positioned at d 1 =R ¼ 0:24483 and d 2 =R ¼ 0:94485. To compare performance of the system with that of Lee-Whiting [8], we have analyzed the spatial distribution of the magnetic field generated by both systems. The evaluation of the field homogeneity involved the use of the Biot-Savart relation applied to small segments of windings. Results of analysis are presented in Fig. 3a, b, which show contours of constant magnetic field relative to the field at center (field error contours) defined as (B given point -B center )/B center plotted at ±1, ±5, and ±10 ppm intervals for the ideal eight-order Lee-Whiting design and our 9/4 system, respectively. It is seen that that the field error contours of both systems compare excellently and the volume of the homogenous magnetic field is similar. In Ref. [11] the Amp-turn ratio of 2.2604 given by Lee-Whiting is approximated by 9/4 (2.25) and the windings placements are left unchanged. Such an approximation expressing the numerator and denominator of the Amp-turn ratio by integers allows the design of systems having coils connected in series and fed by one common source, which is most preferable in practical applications. Unfortunately, as we show in Fig. 4 this small change of only one parameter has a significant impact on the system performance, leading to a very substantial reduction of the generated field homogeneity.
In Table 1 we specify the coil's position with an accuracy of 5 significant digits. With the same accuracy the location of windings is given in the Lee-Whiting paper. Such a precision, however, would be difficult to achieve when constructing the real system. To evaluate the influence of the location of coils on the field homogeneity, we have calculated the field error contours of both systems assuming the coils placement specified with the accuracy of 3 significant digits, i.e., we have assumed d 1 =R ¼ 0:243 and d 2 =R ¼ 0:941 for the Lee-Whiting system and d 1 =R ¼ 0:245 and d 2 =R ¼ 0:945 for the system presented in this paper. The results of calculation are plotted in Fig. 5a, b for the Lee-Whiting arrangement and our 9/4 system, respectively. When comparing pair of Figs. 3a, 5a with the pair 3b and 5b, we see that the 9/4 restricted system generates the greater volume of homogenous magnetic field and is less sensitive to the precision of the location of coils. The real assemblies with coils connected in series consist of many turns of wire fed with equal current and their Amp-turn ratio (NI) 2 /(NI) 1 depends solely on the number of windings of the outer and inner coils. As an example of such assemblies, we have analyzed the arrangements comprising nine circular current loops in the outer coils and four loops in the inner coils. The loops are deposited side by side around z coordinates corresponding to the Lee-Whiting and our 9/4 contracted designs. The expected distributions of the magnetic field generated by the two designs are shown in Fig. 6a, b, respectively. It is seen that the setup suitable for the serial power supply based on our 9/4 system produces much greater volume of the homogeneous magnetic field than that using the Lee-Whiting coordinates.
The setup based on the 9/4 contracted system having the unit radius (R = 1 m) and fed with the current 1 A generates magnetic field whose value in the center is 89.5 9 10 -7 T (89.5 mG). The field homogeneous to 50 ppm, which meets, e.g., EPRI requirements, extends axially to 0.22R and radially to 0.3R from the center. In the case of setup that uses the Lee-Whiting coordinates, the field value is 89.9 9 10 -7 T and the field homogeneous to 50 ppm extends axially to 0.11R and radially to 0.17R from the center.
The magnetic field homogeneity has been analyzed for all systems listed in Table 1. The maximum volume of the homogeneous field is generated by the NI ð Þ 2 = NI ð Þ 1 ¼ 2:26 system. When the Amp-turn ratio changes the volume of the uniform magnetic field slowly decreases. For the last system listed in Table 1 the field homogeneous to 1 ppm extends axially to 0.11R and radially to 0.13R. The value of field is 142 9 10 -7 T.

Shielding of the System
The distant (stray) field of the two-pair coil system may be easily shielded. To show the behavior of the distant axial B z ; we may expand Eq. (6) in inverse powers of z. This yields: Uniformity of Magnetic Field Generated by the Two-Pair Coil System 615 The first three turns in the expansion are readily found to be those for k ¼ 1, k ¼ 3 and k ¼ 5 with coefficients: So, we see that the external (stray) field of the system is dominated by the dipolar term and decreases as z À3 . However, in some applications, it is required that the central field be very homogenous and at the same time the stray field be falling off rapidly. Could we somehow suppress the dipolar term of the four-coil system? Yes, and it is rather easy to do. We can surround the system with another one of larger radius with opposed currents arranged so dipolar moments of the systems given by Eqs. (13a, 13b, 13c) cancel.
For example, if the outer system has twice the radius of the inner one, we need the current flowing in it to be four times smaller than that in the inner system to achieve cancellation. Adding the outer system cancels the total dipole moment not effecting the homogeneity of the central field, but slightly reducing the field value. For the case considered it amounts to 12.5 %.The behavior of the stray field is shown in Fig. 7, it does fall off faster to zero when the whole nested system is energized compared with the inner four-coil system alone.
It is clear that the confinement provides a simple system with the reduced stray field, which in some circumstances may prove to be useful.

Conclusions
The possible distributions of two coaxial pairs of circular windings resulting in systems featuring zero octupole and 32 pole moments are given. Analysis given in the paper shows that one of the derived system, which has the Amp-turn ratio of 9/4, generates homogeneous magnetic field whose volume compares excellently with that of the eight-order Lee-Whiting design, but is less sensitive to the impact of current distribution and the precision of windings location. Moreover, the system may be more easily adopted for the serial power supply. The ideal 9/4 system of the unit radius fed with 1 A current generates the magnetic field of 89.5 mG, which is  about 10 times stronger than that yielded by the Helmholtz pair. The field is proportional to the number of Amp-turns and scales with R -1 , hence systems of smaller radius and/or carrying higher currents will generate higher magnetic fields. Yet, to avoid heating problems they will require thermal and current stabilization. This refers particularly to long-time experiments. Based on the design presented, we plan to build 90 G magnet for micro-EPRI experiments, which should shed some light on the problems associated with practical implementation of the system. Finally, we would like to add that the high-field version of the system can be easily shielded by confinement of the system inside another one with bigger radius. This exerts no significant influence on the center field, but reduces greatly the stray field outside.