Editing Halbach array

{{Short description|Special arrangement of permanent magnets}}
{{Use dmy dates|date=December 2019}}
[[File:HallbachArrayField.jpg|thumb|220x220px|The flux diagram of a Halbach array]]
[[File:Halbach2.svg|thumb|220x220px|A Halbach array, showing the orientation of each piece's magnetic field. This array would give a strong field underneath, while the field above would cancel.]]A '''Halbach array''' ({{IPA|de|ˈhalbax|lang}}) is a special arrangement of permanent [[magnet]]s that augments the magnetic field on one side of the array while cancelling the field to near zero on the other side.<ref name="Halbach1980"/><ref name="Halbach1985"/> This is achieved by having a spatially rotating pattern of magnetisation.

The rotating pattern of permanent magnets (on the front face; on the left, up, right, down) can be continued indefinitely and have the same effect. The effect of this arrangement is roughly similar to many horseshoe magnets placed adjacent to each other, with similar poles touching.

This magnetic orientation process replicates that applied by a magnetic recording tape head to the magnetic tape coating during the recording process.  The principle was further described by James (Jim) M. Winey of [[Magnepan]] in 1970, for the ideal case of continuously rotating magnetization, induced by a one-sided stripe-shaped coil.<ref>{{Cite web|url=https://worldwide.espacenet.com/patent/search?q=pn%3DUS3674946A|title=Electromagnetic Transducer, James Winey, Figure 29; United States Patent 3,674,946, filed Dec. 23, 1970|website=www.espacenet.com}}</ref>

The effect was also discovered by [[John C. Mallinson]] in 1973, and these "one-sided flux" structures were initially described by him as a "curiosity", although at the time he recognized from this discovery the potential for significant improvements in [[magnetic tape]] technology.<ref name=":0">{{cite journal | author = Mallinson J.C. | year = 1973| title = One-Sided Fluxes — A Magnetic Curiosity? | journal = IEEE Transactions on Magnetics | volume = 9 | issue = 4| pages = 678–682 | doi = 10.1109/TMAG.1973.1067714 | bibcode = 1973ITM.....9..678M}}</ref>

Physicist [[Klaus Halbach]], while at the [[Lawrence Berkeley National Laboratory]] during the 1980s, independently invented the Halbach array to focus particle accelerator beams.<ref>{{Cite web|url=https://www.eurekalert.org/features/doe/2004-11/ddoe-mlt111104.php|title=Magnetically levitated train takes flight &#124; US Department of Energy Science News &#124; EurekAlert! Science News|website=www.eurekalert.org|access-date=11 February 2019|archive-date=12 February 2019|archive-url=https://web.archive.org/web/20190212070507/https://www.eurekalert.org/features/doe/2004-11/ddoe-mlt111104.php|url-status=dead}}</ref>

==Linear arrays==
{{Multi image
| image1            = Linear Halbach Array (Strong side up).png
| image2            = Linear Halbach Array (Weak side up).png
| direction         = vertical
| total_width       = 220
| footer            = Orientation of strong and weak sides in a linear Halbach array
}}

===Magnetization===
[[File:HalbachArray1.png|thumb|Cancellation of magnetic components resulting in a one-sided flux|220x220px]]

The magnetic flux distribution of a linear Halbach array may seem somewhat counter-intuitive to those familiar with simple magnets or [[solenoids]]. The reason for this flux distribution can be visualised using Mallinson's original diagram (note that it uses the negative ''y'' component, unlike the diagram in Mallinson's article).<ref name=":0" /> The diagram shows the field from a strip of [[ferromagnetic material]] with alternating magnetization in the ''y'' direction (top left) and in the ''x'' direction (top right). Note that the field above the plane is in the ''same'' direction for both structures, but the field below the plane is in ''opposite'' directions. The effect of superimposing both of these structures is shown in the figure.

The crucial point is that the flux will ''cancel below the plane and reinforce itself above the plane''. In fact, any magnetization pattern where the components of magnetization are <math>\pi/2</math> out of phase with each other will result in a one-sided flux. The mathematical transform that shifts the phase of all components of some function by <math>\pi/2</math> is called a [[Hilbert transform]]; the components of the magnetization vector can therefore be any Hilbert-transform pair (the simplest of which is simply <math>\sin(x) \cos(y)</math>, as shown in the diagram above).

[[File:InfiniteHalbachArray.JPG|thumb|The magnetic field around an infinite Halbach array of cube magnets. The field does not cancel perfectly due to the discrete magnets used.|220x220px]]
The field on the non-cancelling side of the ideal, continuously varying, infinite array is of the form<ref>{{cite web |url=http://www.uta.edu/physics/main/resources/ug_seminars/papers/HalbachArrays.doc |title=Concerning the Physics of Halbach Arrays |first=James R. |last=Creel |date=2006 |access-date=31 August 2008  |archive-url=https://web.archive.org/web/20110604160523/http://www.uta.edu/physics/main/resources/ug_seminars/papers/HalbachArrays.doc |archive-date=4 June 2011  }}</ref>

: <math>F(x, y) = F_0 e^{ikx} e^{-ky},</math>
where
: <math>F(x, y)</math> is the field in the form <math>F_x + i F_y</math>,
: <math>F_0</math> is the magnitude of the field at the surface of the array,
: <math>k</math> is the [[wavenumber]] (i.e., the spatial frequency) <math>2\pi / \lambda.</math>

===Applications===
The advantages of one-sided flux distributions are twofold:
* The field is twice as large on the side on which the flux is confined (in the idealized case).
* There is no [[stray field]] produced (in the ideal case) on the opposite side. This helps with field confinement, usually a problem in the design of magnetic structures.

Thus they have a number of applications, ranging from flat [[refrigerator magnet|refrigerator magnets]] through industrial applications such as the brushless [[DC motor]], [[voice coil]]s,<ref>{{Cite web|url=https://patents.google.com/patent/US7368838B2/en|title=High efficiency voice coil motor}}</ref> magnetic drug targeting<ref name="Sarwar2012"/> to high-tech applications such as [[Wiggler (synchrotron)|wiggler]] magnets used in [[particle accelerator]]s and [[free-electron laser]]s.

The [[Inductrack]] [[maglev train]]<ref>{{cite web |author1=Richard F. Post |date=10 October 2005 |title=Toward More Efficient Transport: The Inductrack Maglev System |url=https://gcep.stanford.edu/pdfs/ChEHeXOTnf3dHH5qjYRXMA/09_Post_10_11_trans.pdf |access-date=1 December 2017 |publisher=Lawrence Livermore National Laboratory |archive-date=4 April 2023 |archive-url=https://web.archive.org/web/20230404162449/https://gcep.stanford.edu/pdfs/ChEHeXOTnf3dHH5qjYRXMA/09_Post_10_11_trans.pdf |url-status=dead }}</ref> and Inductrack rocket-launch system<ref name="Tung2001"/> utilize the Halbach array to lift the train by repelling loops of wire in the track.
[[File:Magnetic_viewing_film.jpg|thumb|[[Magnetic viewing film]] showing a flat refrigerator magnet's magnetization|220x220px]]
Flat flexible (not [[Ferrite (magnet)#Hard ferrites|hard ceramic ferrite]]) refrigerator magnets are created with a Halbach magnetization pattern for a stronger holding force when attached to a flat [[ferromagnetic]] surface (e.g. a fridge door) than the holding force from a uniform magnetization. They're made from powdered ferrite mixed in a flexible binder (e.g. plastic or rubber) that is exposed to a Halbach magnetization field pattern as it is [[Extrusion|extruded]], permanently giving the ferrite particles in the magnetic compound this one-sided flux distribution (which can be viewed with [[magnetic viewing film]]).

[[File:Fridge Magnet Halbach.svg|thumb|Flux distribution for a flat refrigerator magnet|220x220px]][[File:HalbachArrayFEL2.png|thumb|Schematic diagram of a free-electron laser|220x220px]]
Scaling up this design and adding a top sheet gives a [[wiggler magnet]], used in [[synchrotron]]s and [[free-electron laser]]s. Wiggler magnets wiggle, or oscillate, an electron beam perpendicular to the magnetic field. As the electrons are undergoing acceleration, they radiate electromagnetic energy in their flight direction, and as they interact with the light already emitted, photons along its line are emitted in phase, resulting in a "laser-like" monochromatic and coherent beam.

The design shown above is usually known as a Halbach wiggler. The magnetization vectors in the magnetized sheets rotate in the opposite senses to each other; above, the top sheet's magnetization vector rotates clockwise, and the bottom sheet's magnetization vector rotates counter-clockwise. This design is chosen so that the ''x'' components of the magnetic fields from the sheets cancel, and the ''y'' components reinforce, so that the field is given by

: <math>H_y \approx \cos(kx),</math>

where ''k'' is the [[wavenumber]] of the magnetic sheet given by the spacing between magnetic blocks with the same magnetization vector.

{{Clear}}

===Variable linear arrays===

[[File:HalbachArraySchematic1.png|thumb|Schematic of a Halbach array consisting of a series magnetized rods|220x220px]]
[[File:HalbachArraySchematic2.png|thumb|Equal-gearing arrangement for a variable Halbach array|220x220px]]
A series of magnetic rods, magnetized perpendicular to their axes, can be arranged into a Halbach array. If each rod is then rotated alternately through 90°, the resultant field moves from one side of plane of the rods to the other, as shown schematically in the figure.

This arrangement allows the field to effectively be turned on and off above or below the plane of the rods, depending on the rotation of the rods. Such a device makes an efficient mechanical magnetic latch requiring no power. A detailed study of this arrangement has shown that each rod becomes a subject to a strong torque from its neighboring rods when rotated.<ref name="Hilton2012"/> However, a simple and efficient solution, providing both stabilization and the ability to rotate each rod alternately, is simply to provide an equal-gearing arrangement on each rod, as shown in the figure.
{{Clear}}

==Cylinder==
[[File:Halbach cylinder.png|thumb|A ferromagnetic cylinder showing various magnetization patterns and magnetic field|220x220px]]
[[File:Halbach cylinder2.png|thumb|right|Cylinder magnetization]]
A '''Halbach cylinder''' is a magnetized cylinder composed of [[ferromagnetic]] material producing (in the idealized case) an intense magnetic field confined entirely within the cylinder, with zero field outside. The cylinders can also be magnetized such that the magnetic field is entirely outside the cylinder, with zero field inside. Several magnetization distributions are shown in the figures.

The direction of magnetization within the ferromagnetic material, in plane perpendicular to the axis of the cylinder, is given by

: <math>M = M_r\left[\cos\left((k - 1)\left(\varphi - \frac{\pi}{2}\right)\right) \widehat{\rho} + \sin\left((k - 1)\left(\varphi - \frac{\pi}{2}\right)\right) \widehat{\varphi}\right],</math>

where ''M<sub>r</sub>'' is the ferromagnetic [[remanence]] (A/m). A positive value of ''k''&nbsp;−&nbsp;1 gives an internal magnetic field, and a negative one gives an external magnetic field.

Ideally, these structures would be created from an infinite-length cylinder of magnetic material with the direction of magnetization continuously varying. The magnetic flux produced by this ideal design would be perfectly uniform and be entirely confined to either the bore of the cylinder or the outside of the cylinder. Of course, the ideal case of infinite length is not realizable, and in practice the finite length of the cylinders produces ''end effects'', which introduce non-uniformities in the field.<ref name="Mhiochain1999"/><ref name="Bjoerk2011"/> The difficulty of manufacturing a cylinder with a continuously varying magnetization also usually leads to the design being broken into segments.

===Applications===
These cylindrical structures are used in devices such as brushless AC motors, magnetic couplings and high-field cylinders. Both brushless motors and coupling devices use multipole field arrangements:
* Brushless motors or alternators typically use cylindrical designs in which all the flux is confined to the centre of the bore (such as ''k''&nbsp;=&nbsp;4 above, a 6-pole rotor) with the AC coils also contained within the bore. Such self-shielding motor or alternator designs are more efficient and produce higher torque or output than conventional motor or alternator designs.
* Magnetic-coupling devices transmit torque through magnetically transparent barriers (that is, the barrier is non-magnetic or is magnetic but not affected by an applied magnetic field), for instance, between sealed containers or pressurised vessels. The optimal torque couplings consists of a pair of coaxially nested cylinders with opposite +''k'' and −''k'' flux magnetization patterns, as this configuration is the only system for infinitely long cylinders that produces a torque.<ref name="Bjoerk2010"/> In the lowest-energy state, the outer flux of the inner cylinder exactly matches the internal flux of the outer cylinder. Rotating one cylinder relative to the other from this state results in a restoring torque.
*Cylindrical Halbach Arrays are used in Portable [[MRI]] Scanners.<ref>{{cite web |url=https://www.stanfordmagnets.com/halbach-magnets-history-types-and-uses.html |title=Halbach Magnets: History, Types, and Uses |last=Marchio |first=Cathy |date=May 23, 2024 |website=Stanford Magnets |access-date=July 31, 2024}}</ref> They offer the potential for a relatively lightweight, low to mid-field system with no [[cryogenics]], a small fringe field, and no electrical power requirements or heat dissipation needs.<ref>{{cite journal |last1=Cooley |first1=C. Z. |last2=Haskell |first2=M. W. |year=2018 |title=Design of Sparse Halbach Magnet Arrays for Portable MRI Using a Genetic Algorithm |journal=IEEE Transactions on Magnetics |volume=54 |issue=1 |pages=1-12 |doi=10.1109/TMAG.2017.2751001|pmc=5937527 }}</ref> The reduced stray fields also enhance safety and minimize interference with surrounding electronic devices.<ref>{{cite journal |last1=Sarwar |first1=A. |last2=Nemirovski |first2=A. |year=2012 |title=Optimal Halbach permanent magnet designs for maximally pulling and pushing nanoparticles |journal=Journal of Magnetism and Magnetic Materials |volume=234 |issue=5 |pages=742-754 |doi=10.1016/j.jmmm.2011.09.008|pmc=3547684 }}</ref>

===Uniform fields===
[[File:Halbach array by Zureks.png|thumb|Uniform field inside Halbach cylinder|220x220px]]
For the special case of ''k'' = 2, the field inside the bore is uniform and is given by

: <math>H = M_r \ln\left(\frac{R_\text{o}}{R_\text{i}}\right) \widehat{y},</math>

where the inner and outer cylinder radii are ''R''<sub>i</sub> and ''R''<sub>o</sub> respectively. ''H'' is in the ''y'' direction. This is the simplest form of the Halbach cylinder, and it can be seen that if the ratio of outer to inner radii is greater than [[e (mathematical constant)|''e'']], the flux inside the bore actually exceeds the [[remanence]] of the magnetic material used to create the cylinder. However, care has to be taken not to produce a field that exceeds the coercivity of the permanent magnets used, as this can result in demagnetization of the cylinder and the production of a much lower field than intended.<ref name="Bjoerk2015"/><ref name="Insinga2016"/>

[[File:Halbach cylinder3.png|thumb|Three designs (A) (B) (C) producing uniform magnetic fields within their central air gap|220x220px]]
This cylindrical design is only one class of designs that produce a uniform field inside a cavity within an array of permanent magnets. Other classes of design include wedge designs, proposed by Abele and Jensen, in which wedges of magnetized material are arranged to provide uniform field within cavities inside the design as shown.

The direction of magnetization of the wedges in (A) can be calculated using a set of rules given by Abele and allows for great freedom in the shape of the cavity. Another class of design is the magnetic mangle (B), proposed by Coey and Cugat,<ref name="Coey2003"/><ref name="Cugat1998"/> in which uniformly magnetized rods are arranged such that their magnetization matches that of a Halbach cylinder, as shown for a 6-rod design. This design greatly increases access to the region of uniform field, at the expense of the volume of uniform field being smaller than in the cylindrical designs (although this area can be made larger by increasing the number of component rods). Rotating the rods relative to each other results in many possibilities, including a dynamically variable field and various dipolar configurations. It can be seen that the designs shown in (A) and (B) are closely related to the ''k''&nbsp;=&nbsp;2 Halbach cylinder. Other very simple designs for a uniform field include separated magnets with soft iron return paths, as shown in figure (C).

In recent years, these Halbach dipoles have been used to conduct low-field [[NMR]] experiments.<ref>{{cite journal |author=Raich, H. |author2=Blümler, P. |title=Design and construction of a dipolar Halbach array with a homogeneous field from identical bar magnets: NMR Mandhalas |journal=Concepts in Magnetic Resonance Part B: Magnetic Resonance Engineering |date=21 October 2004 |volume=23B |pages=16–25 |doi=10.1002/cmr.b.20018|s2cid=58309210 }}</ref> Compared to commercially available ([[Bruker]] Minispec) standard plate geometries (C) of permanent magnets, they, as explained above, offer a huge bore diameter, while still having a reasonably homogeneous field.

====Derivation in the ideal case====
The method used to find the field created by the cylinder is mathematically very similar to that used to investigate a uniformly magnetised sphere.<ref>{{cite web
 | title = Uniformly Magnetized Sphere
 | url = https://farside.ph.utexas.edu/teaching/jk1/Electromagnetism/node61.html
 | last1 = Fitzpatrick
 | first1 = Richard
 | date = 2014-06-27
 | access-date =  2022-03-28
 | publisher = [[The University of Texas at Austin]]
}}</ref>

Because of the symmetry of the arrangement along the cylinder's axis, the problem can be treated as two-dimensional. Work in [[Polar coordinate system|plane-polar coordinates]] <math>(r, \theta)</math> with associated unit vectors <math>\hat\mathbf{r}</math> and <math>\hat\boldsymbol{\theta}</math>, and let the cylinder have radial extent <math>r_\mathrm{i} < r < r_\mathrm{o}</math>. Then the [[Magnetization|magnetisation]] in the cylinder walls, which has magnitude <math>M_0</math>, rotates smoothly as
:<math>\mathbf{M}_\mathrm{cyl} = M_0 (\cos\theta\, \hat\mathbf{r} + \sin\theta\, \hat\boldsymbol{\theta}),</math>
while the magnetisation vanishes outside the walls, that is for the bore <math>r < r_\mathrm{i}</math> and surroundings <math>r > r_\mathrm{o}</math>.

By definition, the auxiliary [[Magnetic field|magnetic field strength]] <math>\mathbf{H}</math> is related to the magnetisation and [[magnetic flux density]] <math>\mathbf{B}</math> by <math>\mathbf{B} = \mu_0 (\mathbf{H} + \mathbf{M})</math>. Using [[Gauss's law for magnetism|Gauss' law]] <math>\nabla\cdot\mathbf{B} = 0</math>, this is equivalently
{{NumBlk
 | :
 | <math>\nabla\cdot\mathbf{H} = -\nabla\cdot\mathbf{M}.</math>
 | {{EquationRef|1}}
}}
Since the problem is static there are no free currents and all time derivatives vanish, so [[Ampère's circuital law|Ampère's law]] additionally requires <math>\nabla\times\mathbf{H}=0 \implies \mathbf{H} = \nabla\varphi</math>, where <math>\varphi</math> is the [[magnetic scalar potential]] (up to a sign under some definitions). Substituting this back into the previous Equation {{EquationNote|1}} governing <math>\mathbf{H}</math> and <math>\mathbf{M}</math>, we find that we need to solve
{{NumBlk
 | :
 | <math>\nabla^2\varphi = -\nabla\cdot\mathbf{M},</math>
 | {{EquationRef|2}}
}}
which has the form of [[Poisson's equation]].

Consider now the boundary conditions at the cylinder-air interfaces <math>r = r_\mathrm{i}</math> and <math>r = r_\mathrm{o}</math>. Integrating <math>\nabla \times \mathbf{H} = 0</math> over a small loop straddling the boundary and applying [[Stokes' theorem]] requires that the parallel component of <math>\mathbf{H}</math> is continuous. This in turn requires that <math>\varphi</math> is continuous across the boundary. (More properly this implies that <math>\varphi</math> must differ by a ''constant'' across the boundary, but since the physical quantities we are interested in depend on gradients of this potential, we can arbitrarily set the constant to zero for convenience.) To obtain a second set of conditions, integrate Equation {{EquationNote|1}} across a small volume straddling the boundary and apply the [[divergence theorem]] to find
:<math>\left[\frac{\partial \varphi}{\partial r}\right] = \pm \mathbf{M} \cdot \hat{\mathbf{r}},</math>
where the notation <math>[f]</math> denotes a jump in the quantity <math>f</math> across the boundary, and in our case the sign is negative at <math>r = r_\mathrm{i}</math> and positive at <math>r = r_\mathrm{o}</math>. The sign difference is due to the relative orientation of the magnetisation and the surface normal to the part of the integration volume inside the cylinder walls being opposite at the inner and outer boundaries.

In plane-polar coordinates, the [[divergence]] of a [[vector field]] <math>\mathbf{F} = F_r\hat\mathbf{r} + F_\theta \hat\boldsymbol{\theta}</math> is given by
{{NumBlk
 | :
 | <math>\nabla\cdot\mathbf{F} = \frac{1}{r}\frac{\partial}{\partial r}(r F_r) + \frac{1}{r}\frac{\partial F_\theta}{\partial\theta}.</math>
 | {{EquationRef|3}}
}}
Similarly, the [[gradient]] of a [[scalar field]] <math>f</math> is given by
{{NumBlk
 | :
 | <math>\nabla f = \frac{\partial f}{\partial r}\hat\mathbf{r} + \frac{1}{r}\frac{\partial f}{\partial\theta}\hat\boldsymbol{\theta}.</math>
 | {{EquationRef|4}}
}}
Combining these two relations, the [[Laplace operator|Laplacian]] <math>\nabla^2 f = \nabla\cdot(\nabla f)</math> becomes
{{NumBlk
 | :
 | <math>\nabla^2 f = \frac{1}{r}\frac{\partial}{\partial r}\left(r \frac{\partial f}{\partial r}\right) + \frac{1}{r^2}\frac{\partial^2 f}{\partial\theta^2}.</math>
 | {{EquationRef|5}}
}}
Using Equation {{EquationNote|3}}, the magnetisation divergence in the cylinder walls is
:<math>
\begin{align}
\nabla \cdot \mathbf{M}_\mathrm{cyl} &= \frac{1}{r}\frac{\partial}{\partial r}(r M_0 \cos\theta) + \frac{1}{r} \frac{\partial}{\partial\theta}(M_0 \sin\theta)\\[5pt]
&= \frac{M_0 \cos\theta}{r} + \frac{M_0 \cos\theta}{r}\\[5pt]
&= \frac{2M_0 \cos\theta}{r}.
\end{align}
</math>
Hence Equation {{EquationNote|2}}, which is that we want to solve, becomes by using Equation {{EquationNote|5}}
{{NumBlk
 | :
 | <math>
\frac{1}{r}\frac{\partial}{\partial r}\left(r \frac{\partial \varphi}{\partial r}\right) + \frac{1}{r^2} \frac{\partial^2 \varphi}{\partial\theta^2} =
\begin{cases}
0 & r < r_\mathrm{i}, \\
-\frac{2M_0 \cos\theta}{r} & r_\mathrm{i} < r < r_\mathrm{o}, \\
0 & r > r_\mathrm{0}.
\end{cases}
</math>
 | {{EquationRef|6}}
}}
Look for a [[particular solution]] of this equation in the cylinder walls. With the benefit of hindsight, consider <math>\varphi_\mathrm{p} = r \ln r \cos\theta</math>, because then we have
:<math>
\begin{align}
\frac{1}{r}\frac{\partial}{\partial r}\left(r \frac{\partial\varphi_\mathrm{p}}{\partial r}\right) &= \frac{1}{r}\frac{\partial}{\partial r}\left[r \frac{\partial}{\partial r}(r \ln r \cos\theta)\right]\\[5pt]
&= \frac{\cos\theta}{r}\frac{\partial}{\partial r}\left[r(\ln r + 1)\right]\\[5pt]
&= \frac{\cos\theta}{r}(\ln r + 1 + 1)\\[5pt]
&= \frac{\ln r\cos\theta}{r} + \frac{2 \cos\theta}{r}
\end{align}
</math>
and also
:<math>
\begin{align}
\frac{1}{r^2}\frac{\partial^2 \varphi_\mathrm{p}}{\partial\theta^2} &=  \frac{1}{r^2}\frac{\partial^2}{\partial\theta^2}(r \ln r \cos\theta)\\[5pt]
&= -\frac{\ln r \cos\theta}{r}.
\end{align}
</math>
Hence <math>\nabla^2 \varphi_\mathrm{p} = \frac{2 \cos\theta}{r}</math>, and comparison with Equation {{EquationNote|6}} shows that <math>-M_0 \varphi_\mathrm{p} = - M_0 r \ln r \cos\theta</math> is the appropriate particular solution.

Now consider the homogeneous equation for Equation {{EquationNote|6}}, namely <math>\nabla^2 \varphi_\mathrm{h} = 0</math>. This has the form of [[Laplace's equation]]. Through the method of [[separation of variables]], it can be shown that the general homogeneous solution whose gradient is periodic in <math>\theta</math> (such that all the physical quantities are single-valued) is given by
:<math>\varphi_\mathrm{h} = A_0 + B_0 \theta + C_0 \ln r + \sum_{n=1}^\infty \left(A_n r^n + B_n r^{-n}\right)\cos n\theta + \sum_{n=1}^\infty \left(C_n r^n + D_n r^{-n}\right)\sin n\theta,</math>
where the <math>A_i,\ B_i,\ C_i,\ D_i</math> are arbitrary constants. The desired solution will be the sum of the particular and homogeneous solutions that satisfies the boundary conditions. Again with the benefit of hindsight, let us set most of the constants to zero immediately and assert that the solution is
:<math>
\varphi =
\begin{cases}
\alpha r \cos\theta & r < r_\mathrm{i}, \\
\beta r \cos \theta - M_0 r \ln r \cos\theta & r_\mathrm{i} < r < r_\mathrm{o}, \\
0 & r > r_\mathrm{o},
\end{cases}
</math>
where now <math>\alpha,\ \beta</math> are constants to be determined. If we can choose the constants such that the boundary conditions are satisfied, then by the [[uniqueness theorem for Poisson's equation]], we must have found the solution.

The continuity conditions give
{{NumBlk
 | :
 | <math>\alpha = \beta - M_0 \ln r_\mathrm{i}</math>
 | {{EquationRef|7}}
}}
at the inner boundary and
{{NumBlk
 | :
 | <math>\beta - M_0 \ln r_\mathrm{o} = 0</math>
 | {{EquationRef|8}}
}}
at the outer boundary. The potential gradient has non-vanishing radial component <math>\beta\cos\theta - M_0 \ln r \cos\theta - M_0 \cos\theta</math> in the cylinder walls and <math>\alpha\cos\theta</math> in the bore, and so the conditions on the potential derivative become
:<math>(\beta\cos\theta - M_0 \ln r_\mathrm{i} \cos\theta - M_0\cos\theta) - \alpha \cos\theta = -M_0\cos\theta \implies \beta - M_0 \ln r_\mathrm{i} - \alpha = 0 \implies \alpha = \beta - M_0 \ln r_\mathrm{i}</math>
at the inner boundary and
:<math>0 - (\beta\cos\theta - M_0 \ln r_\mathrm{0} \cos\theta - M_0\cos\theta) = M_0\cos\theta \implies \beta - M_0 \ln r_\mathrm{o} = 0</math>
at the outer boundary. Note that these are identical to Equations {{EquationNote|7}} and {{EquationNote|8}}, so indeed the guess was consistent. Hence we have <math>\beta = M_0 \ln r_\mathrm{o}</math> and <math>\alpha = M_0 \ln\left(\frac{r_\mathrm{o}}{r_\mathrm{i}}\right)</math>, giving the solution
:<math>
\varphi =
\begin{cases}
M_0 \ln\left(\frac{r_\mathrm{o}}{r_\mathrm{i}}\right) r \cos\theta & r < r_\mathrm{i},\\
M_0 \ln\left(\frac{r_\mathrm{o}}{r}\right) r \cos\theta & r_\mathrm{i} < r < r_\mathrm{o},\\
0 & r > r_\mathrm{o}.
\end{cases}
</math>
Consequently, the magnetic field is given by
:<math>
\mathbf{H} = \nabla\varphi =
\begin{cases}
M_0 \ln\left(\frac{r_\mathrm{o}}{r_\mathrm{i}}\right) \cos\theta\,\hat\mathbf{r} - M_0 \ln\left(\frac{r_\mathrm{o}}{r_\mathrm{i}}\right) \sin\theta\,\hat\boldsymbol{\theta} & r < r_\mathrm{i},\\
M_0 \cos\theta \left[\ln\left(\frac{r_\mathrm{o}}{r}\right)-1\right]\,\hat\mathbf{r} - M_0 \ln\left(\frac{r_\mathrm{o}}{r}\right) \sin\theta\,\hat\boldsymbol{\theta} & r_\mathrm{i} < r < r_\mathrm{o},\\
0 & r > r_\mathrm{o},
\end{cases}
</math>
while the magnetic flux density can then be found everywhere using the previous definition <math>\mathbf{B} = \mu_0 (\mathbf{H} + \mathbf{M})</math>. In the bore, where the magnetisation vanishes, this reduces to <math>\mathbf{B}_\mathrm{bore} = \frac{1}{\mu_0}\mathbf{H}_\mathrm{bore}</math>. Hence the magnitude of the flux density there is
:<math>B_\mathrm{bore} = \left|\mathbf{B}_\mathrm{bore}\right| = \sqrt{\cos^2\theta + \sin^2\theta}\frac{M_0}{\mu_0} \ln\left(\frac{r_\mathrm{o}}{r_\mathrm{i}}\right) = \frac{M_0}{\mu_0} \ln\left(\frac{r_\mathrm{o}}{r_\mathrm{i}}\right),</math>
which is independent of position. Similarly, outside the cylinder the magnetisation also vanishes, and since the magnetic field vanishes there, the flux density does too. So indeed the field is uniform inside and zero outside the ideal Halbach cylinder, with a magnitude depending on its physical dimensions.

===Varying the field===
Halbach cylinders give a static field. However, cylinders can be nested, and by rotating one cylinder relative to the other, cancellation of the field and adjustment of the direction can be achieved.<ref>{{Cite web |url=http://www.tipmagazine.com/tip/INPHFA/vol-4/iss-3/p34.pdf |archive-url=https://web.archive.org/web/20060328054612/http://www.tipmagazine.com/tip/INPHFA/vol-4/iss-3/p34.pdf  |title=Tip Magazine: Magnets, Markets, and Magic Cylinders The Industrial Physicist by Michael Coey and Denis Weaire |archive-date=28 March 2006}}</ref> As the outside field of a cylinder is quite low, the relative rotation does not require strong forces. In the ideal case of infinitely long cylinders, no force would be required to rotate a cylinder with respect to the other.

[[File:Halbach_planar_array.svg|thumb|Magnetic levitation using a planar Halbach array and concentric structure windings|220x220px]]

==Sphere==

If the two-dimensional magnetic distribution patterns of the Halbach cylinder are extended to three dimensions, the result is the Halbach sphere. These designs have an extremely uniform field within the interior of the design, as they are not affected by the "end effects" prevalent in the finite-length cylinder design. The magnitude of the uniform field for a sphere also increases to 4/3 the amount for the ideal cylindrical design with the same inner and outer radii. However, for a spherical struction, access to the region of uniform field is usually restricted to a narrow hole at the top and bottom of the design.

The equation for the field in a Halbach sphere is<ref>[http://www.magnetocaloric.com/Code/pub01.htm Permanent magnet based sources of magnetic field] {{Webarchive|url=https://web.archive.org/web/20120424225940/http://www.magnetocaloric.com/Code/pub01.htm |date=24 April 2012 }}.</ref>

: <math>B = \frac{4}{3} B_0 \ln\left(\frac{R_\text{o}}{R_\text{i}}\right).</math>

Higher fields are possible by optimising the spherical design to take account of the fact that it is composed of point dipoles (and not line dipoles). This results in the stretching of the sphere to an elliptical shape and having a non-uniform distribution of magnetization over the component parts. Using this method, as well as soft pole pieces within the design, 4.5&nbsp;[[Tesla (unit)|T]] in a working volume of 20&nbsp;mm<sup>3</sup> was achieved,<ref name="Bloch1998"/> and this was increased further to 5&nbsp;T in 2002,<ref>{{Cite web | url=http://cerncourier.com/cws/article/cern/28598 | title=Record-breaking magnet has five-tesla field |website=CERN Courier| date=22 March 2002 }}</ref> although over a smaller working volume of 0.05&nbsp;mm<sup>3</sup>. As hard materials are temperature-dependent, refrigeration of the entire magnet array can increase the field within the working area further. This group also reported development of a 5.16&nbsp;T Halbach dipole cylinder in 2003.<ref name="Kumada2004"/>

== See also ==
*[[Neodymium magnet]]
*[[Permanent magnet]]
*[[Strong focusing]]
*[[Inductrack]] uses Halbach arrays to generate strong fields for maglev trains
*[[Helmholtz coil]] can give very even magnetic fields

== References ==
{{Reflist|refs=

<ref name="Bloch1998">{{cite journal
  | author = Bloch, F. and Cugat, O. and Meunier, G. and Toussaint, J.C.  
  | title = Innovating approaches to the generation of intense magnetic fields: design and optimization of a 4 Tesla permanent magnet flux source
  | journal = IEEE Transactions on Magnetics
  | volume = 34
  | issue = 5
  | pages = 2465–2468
  | year = 1998
  | doi = 10.1109/20.717567
|bibcode = 1998ITM....34.2465B }}</ref>

<ref name="Kumada2004">{{cite journal
  | author = Kumada, M. and Antokhin, E.I. and Iwashita, Y. and Aoki, M. and Sugiyama, E.
  | title = Super Strong Permanent Magnet Quadrupole for a Linear Collider
  | journal = IEEE Transactions on Applied Superconductivity
  | volume = 14
  | issue = 2
  | pages = 1287–1289
  | year = 2004
  | doi = 10.1109/TASC.2004.830555
  | url = http://www.slac.stanford.edu/pubs/slacpubs/10000/slac-pub-10248.pdf
| bibcode = 2004ITAS...14.1287K
  | s2cid = 23698444
 }}</ref>

<ref name="Halbach1980">{{cite journal
  | author = Klaus Halbach
  | title = Design of permanent multipole magnets with oriented rare earth cobalt material
  | journal = Nuclear Instruments and Methods
  | volume = 169
  | number = 1
  | pages = 1–10
  | year = 1980
  | issn = 0029-554X
  | doi = 10.1016/0029-554X(80)90094-4
  | url = http://www.askmar.com/Magnets/Permanent%20Multipole%20Magnets%20Design.pdf
  | bibcode = 1980NucIM.169....1H
  | s2cid = 53486791
  | access-date = 12 January 2016
  | archive-date = 31 December 2018
  | archive-url = https://web.archive.org/web/20181231103639/http://www.askmar.com/Magnets/Permanent%20Multipole%20Magnets%20Design.pdf
  | url-status = dead
  }}</ref>

<ref name="Halbach1985">{{cite journal
  | author = Klaus Halbach
  | title = Applications of Permanent Magnets in Accelerators and Electron Storage Rings
  | journal = Journal of Applied Physics
  | volume = 57
  | number = 1
  | pages = 3605–3608
  | year = 1985
  | issn = 0029-554X
  | doi = 10.1063/1.335021
  | url = http://www.askmar.com/Magnets/Permanent%20Magnet%20Applications.pdf
  | bibcode = 1985JAP....57.3605H
  | access-date = 13 June 2017
  | archive-date = 6 April 2012
  | archive-url = https://web.archive.org/web/20120406184154/http://www.askmar.com/Magnets/Permanent%20Magnet%20Applications.pdf
  | url-status = dead
  }}</ref>

<ref name="Mhiochain1999">{{cite journal
  | author1 = T. R. Ni Mhiochain
  | author2 = D. Weaire 
  | author3 = S. M. McMurry
  | author4 = J. M. D. Coey
  | title = Analysis of torque in nested magnetic cylinders
  | journal = Journal of Applied Physics
  | volume = 86
  | issue = 11 
  | pages = 6412–6424
  | year = 1999
  | doi=10.1063/1.371705
| bibcode = 1999JAP....86.6412N
  }}</ref>

<ref name="Bjoerk2011">{{cite journal
  | author1 = R. Bjørk
  | title = The ideal dimensions of a Halbach cylinder of finite length
  | journal = Journal of Applied Physics
  | volume = 109
  | issue = 1
  | pages = 013915–013915–6
  | year = 2011
  | doi = 10.1063/1.3525646
  | arxiv = 1410.0496| bibcode = 2011JAP...109a3915B
  | s2cid = 119168717
 }}</ref>

<ref name="Bjoerk2010">{{cite journal
  | author1 = R. Bjørk
  | author2 = A. Smith
  | author3 = C. R. H. Bahl
  | title = Analysis of the magnetic field, force, and torque for two-dimensional Halbach cylinders
  | journal = Journal of Magnetism and Magnetic Materials
  | volume = 322
  | issue = 1
  | pages = 133–141
  | year = 2010
  | doi = 10.1016/j.jmmm.2009.08.044
  | url = http://orbit.dtu.dk/files/100219185/Analysis_of_the_magnetic_field_force_and_torque_for_two_dimensional_Halbach_cylinders.pdf
| arxiv = 1409.1712
  | bibcode = 2010JMMM..322..133B
  | s2cid = 56325133
 }}</ref>

<ref name="Bjoerk2015">{{cite journal
  | author1 = R. Bjørk
  | author2 = A. Smith
  | author3 = C. R. H. Bahl
  | title = The efficiency and the demagnetization field of a general Halbach cylinder
  | journal = Journal of Magnetism and Magnetic Materials
  | volume = 384
  | pages = 128–132
  | year = 2015
  | doi = 10.1016/j.jmmm.2015.02.034
  | url = http://orbit.dtu.dk/files/106679364/The_efficiency_and_the_demagnetization_postprint.pdf
| arxiv = 1502.06700
  | bibcode = 2015JMMM..384..128B
  | s2cid = 54826296
 }}</ref>

<ref name="Insinga2016">{{cite journal
  | author1 = A. R. Insinga
  | author2 = C. R. H. Bahl
  | author3 = R. Bjørk
  | author4 = A. Smith
  | title = Performance of Halbach magnet arrays with finite coercivity
  | journal = Journal of Magnetism and Magnetic Materials
  | volume = 407
  | pages = 369–376
  | year = 2016
  | doi = 10.1016/j.jmmm.2016.01.076
| bibcode = 2016JMMM..407..369I
  | s2cid = 124300587
 }}</ref>

<ref name="Coey2003">{{cite book
  |author1=J. M. D. Coey |author2=T.R. Ní Mhíocháin | chapter = Permanent Magnets
  | title = High Magnetic Fields: Science and Technology
  |editor1=F. Herlach |editor2=N. Miura | publisher = World Scientific Publishing
  | volume = 1
  | pages = 25–47
  | year = 2003
  | isbn = 978-981-02-4964-9
}}</ref>

<ref name="Cugat1998">{{cite journal
  | author1 = O. Cugat
  | author2 = F. Bloch 
  | author3 = J.C. Toussaint
  | title = 4-Tesla Permanent Magnetic Flux Source
  | journal = Proc. 15th International Workshop on Rare Earth Magnets and Their Applications
  | page = 807
  | year = 1998
}}</ref>

<ref name="Hilton2012">{{cite journal
  | author1 = J. E. Hilton
  | author2 = S. M. McMurry
  | title = An adjustable linear Halbach array
  | journal = Journal of Magnetism and Magnetic Materials
  | volume = 324
  | issue = 13
  | pages = 2051–2056
  | year = 2012
  | doi = 10.1016/j.jmmm.2012.02.014
  | bibcode = 2012JMMM..324.2051H
| hdl = 2262/63909
  | url = http://www.tara.tcd.ie/bitstream/2262/63909/1/An%20adjustable%20linear%20Halbach%20array.pdf
  | hdl-access = free
  }}</ref>

<ref name="Sarwar2012">{{cite journal
  | author1 = A. Sarwar
  | author2 = A. Nemirovski
  | author3 = B. Shapiro
  | title = Optimal Halbach permanent magnet designs for maximally pulling and pushing nanoparticles
  | journal = Journal of Magnetism and Magnetic Materials
  | volume = 324
  | number = 5
  | pages = 742–754
  | year = 2012
  | doi = 10.1016/j.jmmm.2011.09.008
  | pmid = 23335834
  | bibcode = 2012JMMM..324..742S
  | url = http://www.controlofmems.umd.edu/publications/journal-papers/2011/SarwarEtAl_OptimalHalbachArrays_JMMM_Mar11.pdf
| pmc = 3547684
  }}</ref>

<ref name="Tung2001">{{cite journal
  | author1 = L. S. Tung
  | author2 = R. F. Post
  | author3 = J. Martinez-Frias
  | title = Final progress report for the NASA Inductrack model rocket launcher at the Lawrence Livermore National Laboratory
  | date = 27 June 2001
  | url = http://www.askmar.com/Inductrack/2001-6-27%20Inductrack%20NASA.pdf
  | id = UCRL-ID-144455
  | access-date = 12 January 2016
  | archive-url = https://web.archive.org/web/20160305122350/http://www.askmar.com/Inductrack/2001-6-27%20Inductrack%20NASA.pdf
  | archive-date = 5 March 2016
  
  }}</ref>

}}

== External links ==
*{{Commons category-inline}}
* [https://web.archive.org/web/20040405022851/http://dmtwww.epfl.ch/~jsandtne/GPMB/project1.htm Passive levitation of a shaft]
* [https://web.archive.org/web/20080214163302/http://goldeneye.ethz.ch/motoren/electric/inrunner/index_EN/ Electric model aircraft motor]
* [https://web.archive.org/web/20040324213841/http://reeserail.com/ magnet levitated vehicle switching system]

{{DEFAULTSORT:Halbach Array}}
[[Category:Magnetic levitation]]
[[Category:Magnetic devices]]
[[Category:Types of magnets]]