Generalised Magnetic Polarizability Tensors

We present a new complete asymptotic expansion for the low frequency time-harmonic magnetic field perturbation caused by the presence of a conducting (permeable) object as its size tends to zero for the eddy current regime of Maxwell's equations. The new asymptotic expansion allows the characterisation of the shape and material properties of such objects by a new class of generalised magnetic polarizability tensors and we provide an explicit formula for their calculation. Our result will have important implications for metal detectors since it will improve small object discrimination and, for situations where the background field varies over the inclusion, this information will be useable, and indeed useful, in identifying their shape and material properties. Thus, improving the ability of metal detectors to locate landmines and unexploded ordnance, sort metals in recycling processes, ensure food safety as well as enhancing security screening at airports and public events.


Introduction
The characterisation of highly conducting objects from low frequency magnetic field perturbations has important applications in metal detection where the goal is to locate and identify concealed inclusions in an otherwise low conducting background. Metal detectors are used in the search for artefacts of archaeological significance, the detection of landmines and unexploded ordnance, the recycling of metals, ensuring food safety as well as in security screening at airports and at public events. The ability to better characterise objects offers considerable advantages in reducing the number of false positives in metal detection and, in particular, to accelerate and improve the accuracy of object location and discrimination. For a range of electromagnetic and acoustic phenomena similar findings have been found where, in each case, an asymptotic expansion of the field perturbation caused by the presence of an inclusion as its size, α, tends to zero results in formula which permits the low-cost characterisation of an object. In particular, in electromagnetics, the expansion has been found to be of the form pU α´U 0 qpxq i " pD 2 x Gpx, zqq ij A jk pU 0 pzqq k`p Rq i , when orthonormal coordinates and Einstein summation convention is used. In this expression pU α´U 0 qpxq represents the perturbed field at location x, the object is assumed to have the form B α " αB`z, where B is its shape and z denotes its position, U 0 pxq is the background field, D 2 x Gpx, zq is the Hessian of an appropriate Green's function, R is a residual term and A is a symmetric rank 2 polarizability tensor. The polarizability tensor is independent of z and is a function of B, and hence its topology 1 , it is also a function of the object's material characteristics and, thus, provides a means for its characterisation. Importantly, (1) separates an object characteristics from the background field and, consequently, has applications in the low-cost identification of hidden targets in inverse problems [7,6]. Explicit formulae for the calculation of polarizability tensors have been found in a range of different electromagnetic applications. These include the leading order term in an expansion of the perturbed magnetic field pH α´H 0 qpxq as α Ñ 0 for a (multiply connected) permeable object with β 0 pBq " 1, β 1 pBq ě 0, β 2 pBq ě 0 in magnetostatics [14] and expansions of pH α´H 0 qpxq, and of the perturbed electric field pE α´E0 qpxq, as α Ñ 0 in electromagnetic scattering by simply connected dielectric, permeable or conducting objects with β 0 pBq " 1, β 1 pBq " β 2 pBq " 0 [8,9,13]. In the aforementioned cases, A " T pcq is a suitably parameterised rank 2 Póyla-Szegö tensor and its coefficients can be computed by solving 3 scalar transmission problems or through the solution of 3 integral equations [7] for given α, B and material contrast c. More recently, for the metal detection problem, Ammari, Chen, Chen, Garnier and Volkov [2] have obtained the leading order term in an expansion of pH α´H 0 qpxq as α Ñ 0 for a highly conducting (multiply connected) object placed in a low-frequency time harmonic background field, H 0 pxq. This expansion was obtained for the eddy current regime of the Maxwell system, which is the relevant case for metal detection, and they showed that the object is characterised by a rank 4 tensor and, therefore, is of a different form to (1). However, for the case of orthonormal coordinates, we have shown that their leading order term does reduce to the same form as in (1) and, in this case, A " } } M is a new complex symmetric rank 2 tensor [12]. In this notation, a single check indicates a reduction in a tensor's rank by one, which is achieved due to skew symmetry of the tensor's coefficients in two of its indices, a double check indicates a reduction in a tensor's rank by two. The coefficients of this tensor are computed by solving 3 vectorial transmission problems and are a function of B, α, the object's conductivity, σ˚, its permeability contrast, µ˚{µ 0 , as well as the frequency of excitation, ω. This result also provides a solid mathematical footing for denoting } } M as the rank 2 magnetic polarizability tensor (MPT), which the electrical engineering community predict for describing the characteristics of a conducting object in metal detection e.g. [17,10,16]. In [14], we have obtained further results which relate the coefficients of } } M to T pµ˚{µ 0 q in the limiting case of ω Ñ 0, and, independent of the value of β 1 pBq. For the limiting case of σ˚Ñ 8, we have also shown that the coefficients of } } M tend to those of T p0q if β 1 pBq " 0. These results allow permeable and non-permeable objects to be distinguished and some topology information to be extracted. Furthermore, we have computed } } M for a range of simply and multiply connected objects using a hp-finite element approach and explored how their coefficients vary over a range of frequencies, within the validity of the eddy current model [14].
Although the leading order term in the expansion of pH α´H 0 qpxq as α Ñ 0 and } } M provide useful information about an object, they impose limitations due to the assumption that H 0 pxq is uniform over the inclusion, with it only being evaluated at z, and that } } M has at most 6 independent complex coefficients. In practical magnetic induction metal detection and testing, however, the H 0 pxq generated by coil arrays is significantly non-uniform over the object unless the distance from the coils is very large. For instance, this is the case in walk-through metal detectors, when there is little space between the coil arrays and the person being tested for a security threat, and in subsurface metal detection, when a metallic object is buried close to the surface. In such situations, the leading order term in the expansion of pH α´H 0 qpxq as α Ñ 0 will not accurately describe the field perturbation and } } M will not provide an accurate object characterisation. Still further, there are difficulties in separating geometrical information from the material contrast c in T pcq [7], and hence the limiting cases of } } M. Finally, if an object has rotational or reflectional symmetries the number of non-zero independent coefficients in a symmetric rank 2 tensor are greatly reduced [12] making discrimination between objects difficult (eg. the independent nonzero coefficients of } } M for a cylinder and a cone are the same due to rotational and reflectional symmetries that are present in both objects, even though the cylinder has an additional mirror symmetry normal to a rotation axis that is not present in a cone). In order to describe the field perturbation more accurately, and better characterise a conducting permeable object, higher order terms in the asymptotic expansion are required. For the electrical impedance tomography (EIT) problem, where the perturbed electric field due to the presence of a small conducting inclusion can be described in terms of the gradient of a scalar potential, Ammari and Kang have obtained a complete asymptotic expansion as α Ñ 0 [4,7]. Here, the object is described by generalised polarizability tensors (GPTs) with the lowest order case agreeing with the rank 2 tensor T pcq. This class of GPTs satisfies symmetry conditions on the space of harmonic polynomials. Complete asymtopic expansions of the perturbed field for a small object have also been obtained for acoustics [5] and the elasticity problem [7] where the object is again characterised by GPTs. In this work, we provide a new complete asymptotic expansion of pH α´H 0 qpxq for a highly conducting (possibly permeable and multiply connected) object as α Ñ 0. Thus, extending the result in [2,12], which provided only the leading order term. We write our result in terms of a new class of (higher rank) generalised magnetic polarizability tensors (GMPTs), which characterise the object's shape and its material characteristics. The GMPTs we obtain are quite different to the GPTs previously presented by Ammari and Kang. The explicit expression for their coefficients are with respect to the standard orthonormal basis rather than the space of harmonic polynomials. They are functions of B, α, σ˚, µ˚{µ 0 , ω and can be computed by solving a generalised form of the vectorial transmission problem obtained in [2,12]. Moreover, the leading order term in our new expansion agrees with our previous result [12] and here the GMPT agrees with } } M. Our new complete expansion will overcome the aforementioned difficulties of just using the leading order term for pH α´H 0 qpxq as α Ñ 0 and describing the object using } } M when H 0 pxq is non-uniform over the object, such as in a walk-through metal detector for a security threat and in subsurface metal detection for a metallic object buried close to the surface. For such applications, it will improve the accuracy of pH α´H 0 qpxq, by including more terms in the expansion, and will improve the characterisation of an object's shape and its material properties, by the increased number of independent parameters in the GMPTs. Consequently, improving object identification and location. The paper is organised as follows: In Section 2, the problem configuration is briefly described and some notation is introduced for the presentation of our new results. Section 3 summaries the previously known results about } } M and the leading order term of pH α´H 0 qpxq as α Ñ 0 due to the presence of a highly conducting object [2,12]. Section 4 states our new main result and Section 5 contains the lemmas associated with the derivation of our asymptotic formula. Finally, Section 6 is concerned with the representation of the asymptotic formula in terms of a new class of higher order GMPTs.

Problem configuration
The problem configuration has already been described in [2,12] and is briefly recalled. We consider an electromagnetic inclusion in R 3 of the form B α " z`αB, where B Ă R 3 is a bounded, smooth domain. Let Γ and Γ α denote the boundary of B and B α , respectively, and µ 0 the permeability of free space. We continue to follow the previous notation and write where µ˚and σ˚denote the object's permeability and conductivity, respectively, which we assume to be constant. The time harmonic fields E α and H α that result from a compactly supported time varying current source, J 0 , located away from B α and satisfying ∇¨J 0 " 0 in R 3 , and their interaction with the object B α , satisfy the eddy current equations in a weak sense [1] where ω denotes the angular frequency and i :" ?´1 . Letting α " 0 in (3) we obtain the corresponding fields, E 0 and H 0 , that result from time varying current source in the absence of an object. As explained in [1], the eddy current model is completed by ∇¨E α " 0 in R 3 zB α and the uniqueness of E α in R 3 zB α is achieved by additionally specifying where n is the outward normal to Γ α . Furthermore, in practice, the decay of the fields is actually faster than the |x|´1 stated in the original eddy current model [1]. The task is to develop an asymptotic expansion for pH α´H 0 qpxq as α Ñ 0 for the case where ν :" kα 2 " Op1q, k :" ωµ 0 σ˚, which includes the case of fixed σ˚, µ˚, ω as α Ñ 0 (since in this case |ν| ď Cα 2 ď C). Notice that the condition on ν is required to ensure the eddy current model is not violated as the object size vanishes.
For what follows it is beneficial to introduce the following notation: Definition 2.1. We will use boldface for vector quantities (e.g. u) and denote by e j , j " 1, 2, 3 the units vectors associated with an orthonormal coordinate system. We denote the j-th component of a vector u in this coordinate system by puq j " u¨e j " u j .
Definition 2.2. We will use calligraphic symbols to denote rank 2 tensors e.g. N " N ij e i b e j and denote their coefficients by N ij . Definition 2.3. By symbols in the Fraktur font, e.g. A, we shall denote higher order tensors and, to describe their coefficients with respect to an orthonormal coordinate basis, it will be useful to introduce the p-tuple of positive integers Jppq :" rj 1 , j 2 ,¨¨¨, j p s, the m-tuple of positive integers Kpmq :" rk 1 , k 2 ,¨¨¨, k m s and to introduce the pp`1q-and pm`1q-tuple of positive integers Jpp`1q ": rj, Jppqs and Kpm`1q " rk, Kpmqs, respectively. Thus by N rj,Jppq,k,Kpmqs " N Jpp`1qKpm`1q 2 we denote the coefficients of the rank 2`p`m tensor 2 When no confusion arises, we will drop the square parentheses on the lists of indices For p " m " 0 this reduces to the rank 2 tensor N " N kj e k b e j " N kj e k b e j . Consider the rank 4`m`p tensor often we will write A rh,i,k,Kpmq,j,Jppqs " A rrh,i,k,Kpmqs,rj,Jppqss " A rrh,i,Kpm`1qs,Jpp`1qs to group indices and assist when considering products with other terms as well as when considering the skew symmetries of this tensor. However, by the introduction of such brackets, we do not imply skew systematization over these indices. Using skew symmetries, we will denote by a single check (i.e. q A) a reduction in rank by 1 and by a double check (i.e. q q A) a reduction in its rank by 2.
where ξ¨e j " ξ j is the jth spatial coordinate measured from an origin contained in B. Furthermore, when p " 0 (or m " 0 ) then J " H (respt. K " H) and, in this case, pΠpξqq H " 1.
Using this notation, we shall imply the Einstein summation convention for repeated sets of indices as appropriate.
Definition 2.5. We recall that for 0 ď ă 8, 0 ď p ă 8, where the derivatives are defined in a weak sense and 3 Leading order term of pH α´H 0 qpxq as α Ñ 0 In [2], Ammari et al. have obtained the leading order term in an asymptotic expansion of pH αH 0 qpxq as α Ñ 0 for ν " Op1q and x away from the location z of the inclusion. In [12] we have previously shown that their result can be conveniently expressed using Einstein summation notation and, in the case of orthonormal coordinates, that it reduces to with |Rpxq| ď Cα 4 }H 0 } W 2,8 pBαq . In the above, pD 2 x Gpx, zqq ik are the components of the rank 2 tensor D 2 x Gpx, zq " 1{p4π|r| 3 qp3r br´Iq " pD 2 x Gpx, zqq ik e i b e k . This is obtained from pD 2 x Gpx, zqq ik :" B x i B x k pGpx, zqq where Gpx, zq :" 1{p4π|x´z|q " 1{p4π|r|q, r " x´z, r " r{|r| and pIq ik " δ ik are the components of the identity tensor, which are equal to the Kronceker delta. Furthermore, we have shown that } } M kj :"´q C kj`Nkj are the coefficients of a complex symmetric magnetic polarizability (MPT) rank 2 tensor, which describes the shape, conductivity, and permeability (contrast) of the object, and is computed using These, in turn, rely on the vectoral solutions θ j , j " 1, 2, 3 to the transmission problem where r¨s Γ denotes the jump of the function over Γ. Note that θ j ‰ pθq j , the latter being the jth component of a vector. Instead, the subscript j should be interpreted as the jth solution of the transmission problem corresponding to the source terms in B and on Γ being constructed from the jth unit vector e j . In the above, we have dropped the subscript α on the position dependent µ as this problem is formulated for the object B rather than B α .
4 Complete asymptotic expansion of pH α´H 0 qpxq as α Ñ 0 Our main result is Theorem 4.1. The magnetic field perturbation in the presence of a small conducting object B α " αB`z for the eddy current model when ν is order one and x is away from the location z of the inclusion is completely described by the asymptotic formula Rpxqq i , (9) Jpp`1q :"rj, Jppqs " rj, j 1 , j 2 ,¨¨¨, j p s, Kpm`1q :"rk, Kpmqs " rk, k 1 , k 2 ,¨¨¨, k m s, In the above, Jppq and Kpmq are pand m-tuples of integers, respectively, with each index taking values 1, 2, 3. Also and the coefficients of a rank 2`p`m generalised magnetic polarizability tensor (GMPT) are defined by where ξˆ`pΠpξqq Kpmq pθ Jpp`1q`p Πpξqq Jppq e jˆξ q˘dξ, In the above, θ Jpp`1q satisfy the transmission problem Proof. The expansion follows from the asymptotic formula in Theorem 5.6, and the results in Lemmas 6.3 and 6.4 by noting that the coefficients of the rank 4`m`p tensor A can be expressed in terms of the coefficients of a rank 3`m`p and then a rank 2`m`p tensor by using the skew symmetry of their coefficients. Explicitly, we find that A rri, ,k,Kpmqs,Jpp`1qs "ε ikr C rrr, ,Kpmqs,Jpp`1qs " ε rs ε ikr q C rrs,Kpmqs,Jpp`1qs "pδ k δ si´δ i δ sk q q C rrs,Kpmqs,Jpp`1qs " δ k q C rri,Kpmqs,Jpp`1qs´δ i q C rrk,Kpmqs,Jpp`1qs where ε is as defined in (53) and we have used ε rs ε ikr "´ε r s ε rik " δ k δ si´δ i δ sk . Finally, by using x Gpx, zqq " 0, and by a term by term application of the above arguments, (9) is obtained.
Remark 4.2. Theorem 4.1 extends the asymptotic expansion obtained by [2], which provides the leading order term for pH α´H 0 qpxq as α Ñ 0. We have previously shown in [12] that this leading order term can be written in the alternative form presented in (6). In this case, B, α, σå nd µ r " µ˚{µ 0 are described by a complex symmetric rank 2 MPT } } M, which is also a function of ω. However, this description can only provide limited amounts of information about an object. Our new result reduces to this case when M " 1. For M " 2, our new result gives with |Rpxq| ď Cα 5 }H 0 } W 3,8 pBαq . In the above, | | M rk,rj,j 1 ss "´q C rk,rj,j 1 ss`Nrk,rj,j 1 ss where q C rk,rj,j 1 ss "´i να 4 6 e k¨ż B ξˆ`θ rj,j 1 s`p ξq j 1 e jˆξ˘d ξ, (12a) and these can be computed using the solution of (8) already found for the computation of } } M kj . For the computation of (12) the solution of is also required. In the case of M Ñ 8, Theorem 4.1 provides a complete description of the field perturbation pH α´H 0 qpxq caused by the presence of a permeable conducting object as α Ñ 0. The object's shape and material properties in our new result are described by | | M, which are GMPTs of increasing rank up to a maximum of M`1 and are again functions of B, α, σ˚, µ r and ω. By applying Theorem 4.1 with M ą 1, pH α´H 0 qpxq can be more accurately described by including more of the higher order terms. Complete asymptotic field expansions for small objects and GPTs and have previously been obtained for the EIT problem, acoustics and elasticity [4,5,7]. But, to the best of the authors' knowledge, this is the first time they have been obtained for a Maxwell problem. Like in [2], our analysis makes the assumption that B has a smooth boundary. The extension of the analysis to non-smooth boundaries will form part of our future work. However, numerical evidence from computing M for objects with edges indicates that our results are also likely to hold for such objects [3,12,14].

Remark 4.3.
To be able to characterise an unknown conducting permeable object from measurements of pH α´H 0 qpxq, using Theorem 4.1, a range of alternative approaches are possible, which include adapting the algorithms described by Ammari and Kang [6] for the EIT problem or using a statistical classifier [3]. In the latter case, we assume that we have a set of possible candidate objects and we follow Ammari and Kang [6, pg. 80] to put these in canonical form such that the description B α " αB`z, for each object, implies that the origin for ξ coincides with the object's centre of mass and that the determinant of the Póyla-Szegö tensor associated with B (i.e. T pµ r qrBs for µ r ‰ 1 and T p0qrBs for µ r " 1), it is equal to 1 [6, pg. 80]. In an off-line stage, the coefficients of | | M are then computed numerically for these objects for a range of frequencies ω by solving the transmission problem for θ Jpp`1q to form a dictionary. In an on-line stage, the unknown object's position z can be found by rotating the candidate objects in the dictionary (and hence their } } M) and determining the best statistical fit for z by using measurements of pH α´H 0 qpxq and (6). To find the unknown object's size, it may be necessary to ensure the dictionary includes } } M computed at very small ω, or frequencies at the limit of the eddy current model, such that } } M reduces to a suitably parameterised Póyla-Szegö tensor for the candidate objects [14], and then repeat the above process to find the best fit for α. To determine further geometrical and material parameter information, measurements of pH α´H 0 qpxq will be compared against Theorem 4.1 by using the known z and rotating the candidate objects in the dictionary (and hence their | | M) to find the best fit.
Remark 4.4. Currently, in practical magnetic induction metal detection, rather than solving (3) with α " 0 to obtain the true background magnetic field, engineers frequently approximate the field at a position x obtained from a small coil centred at y as that of a magnetic dipole where m e j is a constant vector that is a function of the coil's dimensions and the current flowing in it. However, in walk through metal detectors, where there is little space between the coil arrays and the person being tested for a security threat, this does not provide an accurate representation of the field as the coils dimensions are no longer small compared to |x´y| and the background field can vary considerably over the object. Engineers also assume that measurement coil, if treated as an emitter, will act as a dipole source. This means that for a single emittermeasurement coil arrangement the induced voltage in a measurement coil located at x would be of the form of the leading order term for pH α´H 0 qpxq [14] pm m q i pH α´H 0 qpxq i « pm m q i pD 2 Gpx, zqq ij } } M jk pD 2 Gpz, yqq k pm e q .
However, again the assumption of a dipole field for the measurement coil breaks down for walk through metal detectors and similar problems also exist in subsurface metal detection, when a metallic object is buried close to the surface. Theorem 4.1 can improve the characterisation of hidden objects in magnetic induction metal detection by the following: Instead of (14) the system (3) should be solved with α " 0 to obtain the true H 0 pxq; Rather than (15) pH α´H 0 qpxq should be integrated over an appropriate volume [15] to obtain the correct induced voltage; Instead of just using the leading order term for pH α´H 0 qpxq as α Ñ 0, and an object characterisation using } } M, more terms in (9) should be used, and an object characterised by | | M. Furthermore, object location and identification can then be improved by using the approach described in Remark 4.3.

Eliminating the current source
We will build on [2,12], but, in order to give a physical interpretation, it is first instructive to rewrite the original problem described in Section 2 with α " 0 and α ‰ 0 as transmission problems and to eliminate the current source. To do this, we note that in absence of an object, H 0 " µ´1 0 ∇ˆA 0 and A 0 solves Then, in the presence of the object, we can write H α " µ´1 α ∇ˆA α , E α " iωA α where, after appropriate gauging, A α solves A α " Op|x|´1q as |x| Ñ 8 .
Next, introducing ξ :" x´z α , A ∆ pξq :" αpA α´A0 q´x´z α¯a nd rescaling we see that A ∆ solves the following transmission problem where the current source no longer appears and, instead, is replaced by source terms in (18a) and (18c). Electrical engineers would call A ∆ pξq a scaled reduced vector potential. We now need to represent a polynomial vector field as the curl of another; we call this an uncurling formula.

Uncurling formula
Lemma 5.1. Given a smooth divergence free polynomial vector field in the form spxq " where pD p z psqpzqq Jpp`1q :"˜p satisfies s " ∇ xˆt .
Proof. We consider the pth term in (20) and apply the standard identity ∇ˆpuˆvq " u∇¨vv ∇¨u`pv¨∇qu´pu¨∇qv where the differentiation is with respect to x and puq j " 1 p!pp`2q pD p z pspzqqq Jpp`1q pΠpx´xqq Jppq , pvq j " px´zq j .
It is obvious that ∇¨v " 3, p∇vq ji " δ ji and we can deduce ∇¨u " 1 p!pp`2q`p D p z pspzqqq Jpp`1q`δj 1 j px j 2´z j 2 q¨¨¨px jp´zjp q`¨¨¨p x j 1´z j 1 qpx j 2´z j 2 q¨¨¨px j p´1´z j p´1 qδ jpj˘" 1 p!pp`2q`p D p z pspzqqq rj,j,j 2 ,j 3 ,¨¨¨,jps px j 2´z j 2 qpx j 3´z j 3 q¨¨¨px jp´zjp q`¨¨¨p D p z pspzqqq rj,j 1 ,j 2 ,j 3 ,¨¨¨,j p´1 ,js px j 1´z j 1 qpx j 2´z j 2 q¨¨¨px j p´1´z j p´1 q˘" 0, by interchanging the order of differentiation of s (e.g. pD p z pspzqqq rj,j 1 ,j 2 ,j 3 ,¨¨¨,j p´1 ,js " pD p z pspzqqq rj,j,j 1 ,j 2 ,j 3 ,¨¨¨,j p´1 s ) and noting the repeated j index, which imples pD z pspzqqq jj " trpD z pspzqqq " ∇ z¨s pzq " 0. Note also Thus, and, by a term by term application of the above arguments, (19) immediately follows.
2. An immediate consequence of Lemma 5.1 and the smoothness of the divergence free ∇ˆE 0 " iωµ 0 H 0 in B α is that we can introduce a vector field F pxq as where Jpp`1q " rj, Jppqs and whose curl is the polynomial vector field ∇ˆF " which is also the P th order Taylor series expansion of ∇ˆE 0 about z for |x´z| Ñ 0. Note that (21) and (22) generalise the expressions for F pxq and ∇ˆF stated in [2], which are for the case of P " 1.

Accuracy of Taylor series approximations
The smoothness of H 0 " 1 iωµ 0 ∇ˆE 0 in B α enables us to deduce that where, here and throughout the following, the constant C is independent of α. Note that in the case of P " 1 (25) is analogous to the bound (3.6) in [2]. It also follows from (25) that
In a departure from [2], we define w PX α pR 3 q as the weak solution to and, if we compare the above with their (3.7), we see it reduces to the latter for P " 1 and also find that their Lemma 3.2 generalises to: There exists a constant C such that }∇ˆpE α´E0´w q} L 2 pBαq ď C`|1´µ´1 r |`ν˘α where µ r :" µ˚{µ 0 and ν is as defined in (5).
Using this result and wpxq " αw 0`x´z α˘w e find that Theorem 3.1 in [2] immediately generalises to: Theorem 5.4. There exists a constant C such that Repeating their steps for the higher order terms we find that where Remark 5.5. The indices on the solution θ Jpp`1q pξq to the auxiliary problem (40) should be interpreted differently to the tensoral indices previously presented. They should be interpreted as the vector valued solution when the source terms in B and on Γ contain the product e j pΠpξqq Jppq . The transmission problem for θ Jpp`1q pξq is independent of the object's position and is independent of the background excitation. It depends only on the shape of the object, its size, material properties and the frequency of the excitation. It generalises the transmission problem stated in (8), obtained in [2], and reduces to this case when p " 0. We will examine the transformation of θ Jpp`1q pξq under rotations and/or reflections of the object in the proof of Lemma 6.1.
Alternatively, by using (23) and (24), and substituting in to the source terms in B and on Γ in (18), we see that it is possible to write which, by truncating at P terms and multiplying by iω, coincides with the weak solution w 0 pξq and a hence provides a physical interpretation for the latter.

Integral representation formula
Ammari et al. [2] have derived the following Stratton-Chu type formula for x exterior to B α , which relates the magnetic field perturbation outside the object to the magnetic field in its interior. By introducing the representation for B α , and using the eddy current equations (3), we have the alternative form x Gpx, yqH α pyqdy "I`II.

Asymptotic formula
One approach to approximating integrals in (42) is to transform the domain of integration from B α to B, to express E α pαξ`zq and H α pαξ`zq in terms of A ∆ pξq and A 0 pαξ`zq (and their curls) and then to substitute in truncated expansions of (23) and (41). However, as rigorous estimates for these approximations are not available, we would not be able to quantify the remainder and so, instead, we pursue the previously presented approach in [2], which uses weak solutions and has a rigorous theoretical framework. We extend this approach to the higher order case by using the bounds we have derived in Theorem 5.4 and the result is the following theorem, which generalises their Theorem 3.2.
Theorem 5.6. Let ν be order one and let α be small. For x away from the location z of the inclusion we have where |Rpxq| ď Cα 3`M }H 0 } W M`1,8 pBαq .
Proof. The result follows immediately from Lemmas 5.10 and 5.12 presented in Section 5.7.

Results for the proof of the asymptotic formula
It is useful to note for x away from B α and y in B α that Gpx, yq is smooth and analytic where we plan to use it and so the Taylor series expansions converge as |y´z| Ñ 0. Consequently, we have the estimates
Lemma 5.8. We can bound I 1 , I 2 and I 3 as Proof. Using Theorem 5.4 we have where the second inequality follows from ∇ˆE 0 " iωµ 0 H 0 and (5). The final result for I 1 follows by recalling µ r " Op1q and ν " kα 2 " Op1q. Next, using (36) we find Lemma 5.9. The term corresponding to m " 0 in I 4 vanishes so that pD p z pH 0 pzqqq Jpp`1q pθ Jpp`1q`p Πpξqq Jppq e jˆξ q¸dξ.

Approximation of II
In similar way to [2], we write II "´1´µμ 0¯p II 1`I I 2`I I 3`I I 4 q where pΠpy´zqq Jppq e j¸d y, and so for II 1 we find that For II 2 , we have from (50) and [2] that Finally, for |II 3 |, in a similar manner to (26) we have Lemma 5.12. The integral II can be expressed as Proof. We recall that II "´1´µμ 0¯p II 1`I I 2`I I 3`I I 4 q and, in light of Lemma 5.10, we need to choose P " M . Making the choice of S " M´1 and T " M we see from Lemma 5.11 that the terms associated with II 1 , II 2 and II 3 all form part of Rpxq with conservative estimates in the power of α in II 1 and II 3 , but all involving }H 0 } W M`1,8 pBαq . Also A rri, ,Kpm`1qs,Jpp`1qs rJ pBqs "J ii 1 J 1 J kk 1 J jj 1 J KpmqK 1 pmq J JppqJ 1 ppq A rri 1 , 1 ,K 1 pm`1qs,J 1 pp`1qs rBs, N Kpm`1qJpp`1q rJ pBqs "J kk 1 J jj 1 J KpmqK 1 pmq J JppqJ 1 ppq N K 1 pm`1qJ 1 pp`1q rBs, where the term inside the square parenthesis indicates the object for which the tensor is evaluated, J is an orthogonal transformation matrix and It follows that these arrays of functions are the coefficients of the rank 4`m`p and 2`m`p tensors, A and N, respectively.
Having verified that A and N are tensors we now investigate whether the former rank 4`m`p tensor can instead be represented by a rank 2`m`p tensor.

Reduction of A
Without loss of generality, we assume a positively orientated orthogonal frame. Any change in sign that follows from the reduction in rank for a different frame will cancel as we shall reduce the rank of A by 2.