Rank-2 Magnetic Polarizability Tensor

• The rank-2 magnetic polarizability tensor is a complex symmetric 3×3 matrix that describes the dipole response of small conducting permeable objects.
• It encodes key features such as shape, conductivity, permeability, size, and frequency dependence via asymptotic perturbation analysis.
• Its compact formulation underpins target characterization, inverse problem-solving, and metal detection, while also extending to higher-order tensor models.

Searching arXiv for the cited papers on magnetic polarizability tensors and related formulations. {"query":"all:magnetic polarizability tensor rank-2 Ledger Lionheart Ammari deuteron twisted electron", "max_results": 10} Searching arXiv for relevant papers. The rank-2 magnetic polarizability tensor is, in the standard eddy-current theory of small conducting permeable objects, a complex symmetric second-order tensor that appears as the leading-order coefficient in the asymptotic perturbation of the magnetic field. In Cartesian representation it is a $3\times 3$ tensor, typically denoted $\mathcal M$, with at most six independent complex coefficients, and it encodes shape, size, conductivity, permeability contrast, frequency dependence, and orientation through its tensor transformation law rather than through explicit position dependence (Ledger et al., 2017). Within this literature, “rank-2” refers to a spatial tensor on $\mathbb R^3$, not to a relativistic spacetime tensor, and the object serves as a compact descriptor for target characterization, inverse problems, and metal detection (Wilson et al., 2021).

1. Definition and conceptual scope

In the metal-detection and eddy-current asymptotics literature, the magnetic polarizability tensor is the tensor $\mathcal M$ in the leading perturbation formula

$$(\boldsymbol H_\alpha-\boldsymbol H_0)(\boldsymbol x)_i = (\boldsymbol D_x^2 G(\boldsymbol x,\boldsymbol z))_{ij} (\mathcal M)_{jk} (\boldsymbol H_0(\boldsymbol z))_k +O(\alpha^4),$$

for a small object $B_\alpha=\alpha B+\boldsymbol z$ (Wilson et al., 2020). This identifies $\mathcal M$ as the object-dependent factor in a dipole-type response: it maps the local background magnetic field at the object location into the leading perturbation measured away from the object (Ledger et al., 2019).

The tensor is explicitly treated as complex symmetric and rank 2, equivalently a $3\times 3$ matrix in an orthonormal basis, with six independent complex coefficients in the generic case (Wilson et al., 2021). Its entries scale like $\alpha^3$, consistent with dipolar volume scaling, and its dependence on frequency is governed by the parameter

$\nu=\omega\mu_0\sigma_*\alpha^2,$

or the corresponding piecewise definition for inhomogeneous objects (Ledger et al., 2015). In this formulation, the tensor depends on shape, conductivity, permeability, size, and excitation frequency, but it is independent of object position; translation enters separately through the Green function and background field evaluation (Ledger et al., 2018).

A central point in the later literature is that the rank-2 MPT is the lowest-order member of a larger hierarchy of generalised magnetic polarizability tensors. In the complete asymptotic expansion, the classical rank-2 tensor is recovered by taking the lowest-order indices $\mathcal M$0 and $\mathcal M$1, so the familiar MPT is the leading-order truncation of a higher-order object-characterization framework (Ledger et al., 2022).

2. Eddy-current asymptotics and the leading-order field perturbation

The governing regime is the time-harmonic eddy-current approximation of Maxwell’s equations for a small conducting permeable inclusion. The fields satisfy

$\mathcal M$2

with piecewise-constant conductivity and permeability in the object and its exterior (Ledger et al., 2017). The asymptotic regime assumes $\mathcal M$3 with $\mathcal M$4, so the quasi-static eddy-current scaling is retained while the object size shrinks (Ledger et al., 2015).

The leading-order perturbation formula has the standard dipole structure

$\mathcal M$5

with

$\mathcal M$6

This makes the tensor the coefficient of the leading-order magnetic dipole response (Ledger et al., 2019). In the notation of the GMPT literature, the same term is obtained from the complete expansion by retaining only the lowest-order contribution in derivatives of both the Green tensor and the background field (Ledger et al., 2022).

Historically, the engineering literature often used a rank-2 tensor directly for object characterization, whereas earlier rigorous asymptotics had produced a rank-4 object. The reduction to a complex symmetric rank-2 tensor in orthonormal coordinates is one of the major conceptual clarifications in this line of work (Ledger et al., 2015). This reduction is not merely notational; it establishes that the object signature relevant to leading-order hidden-target characterization can be represented by a second-order tensor with a standard transformation law under rotations (Ledger et al., 2017).

3. Coefficient formulas, transmission problems, and spectral structure

The rank-2 tensor is commonly written as

$\mathcal M$7

with $\mathcal M$8 the conductivity-driven part and $\mathcal M$9 the permeability-contrast part (Ledger et al., 2018). For a homogeneous object, the coefficients are

$\mathbb R^3$0

$\mathbb R^3$1

where the auxiliary fields $\mathbb R^3$2 solve a vector transmission problem on the reference domain and its exterior (Ledger et al., 2019).

Those auxiliary fields satisfy, for $\mathbb R^3$3,

$\mathbb R^3$4

together with continuity of tangential fields, the appropriate jump condition for $\mathbb R^3$5, a divergence constraint, and decay at infinity (Ledger et al., 2017). The tensor is therefore not postulated phenomenologically; it is computed from a well-posed transmission problem tied directly to the eddy-current model.

A further structural decomposition isolates a magnetostatic part and two conductive frequency-dependent pieces: $\mathbb R^3$6 where $\mathbb R^3$7, $\mathbb R^3$8, and $\mathbb R^3$9 are each real symmetric rank-2 tensors (Ledger et al., 2019). This decomposition is especially important because $\mathcal M$0 is the low-frequency magnetostatic contribution, while $\mathcal M$1 and $\mathcal M$2 encode the dispersive conductive response. In the homogeneous case, $\mathcal M$3, and for homogeneous $\mathcal M$4 one has $\mathcal M$5, identifying the low-frequency MPT with the Pólya–Szegő tensor (Ledger et al., 2015).

The spectral analysis of the tensor coefficients is one of the most developed parts of the subject. For homogeneous conductivity, the fields admit an eigenfunction expansion in terms of an auxiliary spectral problem, and the conductive parts satisfy modal formulas of the form

$\mathcal M$6

$\mathcal M$7

with

$\mathcal M$8

These formulas explain why the real part tends to show bounded monotone behavior and the imaginary part tends to exhibit a single peak in the homogeneous case (Ledger et al., 2019).

Because $\mathcal M$9, $$(\boldsymbol H_\alpha-\boldsymbol H_0)(\boldsymbol x)_i = (\boldsymbol D_x^2 G(\boldsymbol x,\boldsymbol z))_{ij} (\mathcal M)_{jk} (\boldsymbol H_0(\boldsymbol z))_k +O(\alpha^4),$$0, and $$(\boldsymbol H_\alpha-\boldsymbol H_0)(\boldsymbol x)_i = (\boldsymbol D_x^2 G(\boldsymbol x,\boldsymbol z))_{ij} (\mathcal M)_{jk} (\boldsymbol H_0(\boldsymbol z))_k +O(\alpha^4),$$1 are real symmetric, each is orthogonally diagonalizable and has real eigenvalues. This makes eigenvalues and principal invariants natural orientation-independent descriptors, whereas the raw tensor entries depend on the lab-frame orientation (Wilson et al., 2021).

4. Generalisations, multiple objects, and computational use

The rank-2 MPT formalism extends beyond a single homogeneous target. For $$(\boldsymbol H_\alpha-\boldsymbol H_0)(\boldsymbol x)_i = (\boldsymbol D_x^2 G(\boldsymbol x,\boldsymbol z))_{ij} (\mathcal M)_{jk} (\boldsymbol H_0(\boldsymbol z))_k +O(\alpha^4),$$2 sufficiently well-separated homogeneous objects $$(\boldsymbol H_\alpha-\boldsymbol H_0)(\boldsymbol x)_i = (\boldsymbol D_x^2 G(\boldsymbol x,\boldsymbol z))_{ij} (\mathcal M)_{jk} (\boldsymbol H_0(\boldsymbol z))_k +O(\alpha^4),$$3, the leading perturbation becomes a sum of $$(\boldsymbol H_\alpha-\boldsymbol H_0)(\boldsymbol x)_i = (\boldsymbol D_x^2 G(\boldsymbol x,\boldsymbol z))_{ij} (\mathcal M)_{jk} (\boldsymbol H_0(\boldsymbol z))_k +O(\alpha^4),$$4 dipole-type terms,

$$(\boldsymbol H_\alpha-\boldsymbol H_0)(\boldsymbol x)_i = (\boldsymbol D_x^2 G(\boldsymbol x,\boldsymbol z))_{ij} (\mathcal M)_{jk} (\boldsymbol H_0(\boldsymbol z))_k +O(\alpha^4),$$5

so each object contributes its own position-independent complex symmetric rank-2 tensor (Ledger et al., 2018). When objects are closely spaced or materially inhomogeneous but scale and translate as a single composite, the cluster is instead characterized by one effective rank-2 tensor computed from a coupled transmission problem over the composite geometry (Ledger et al., 2018).

The complete asymptotic expansion later developed for generalised MPTs shows precisely when the classical rank-2 description ceases to be sufficient. Higher-order tensors enter when the background field varies appreciably over the object or when greater discrimination power is needed beyond the uniform-field dipole approximation (Ledger et al., 2017). This suggests that the rank-2 tensor is the correct first descriptor, but not necessarily the complete one in near-field or strongly nonuniform excitation settings.

Computationally, the tensor is especially attractive because it compresses the forward electromagnetic response into a small number of coefficients. Proper-orthogonal-decomposition reduced-order models were developed to accelerate the evaluation of the spectral signature $$(\boldsymbol H_\alpha-\boldsymbol H_0)(\boldsymbol x)_i = (\boldsymbol D_x^2 G(\boldsymbol x,\boldsymbol z))_{ij} (\mathcal M)_{jk} (\boldsymbol H_0(\boldsymbol z))_k +O(\alpha^4),$$6, with full-order finite element solutions used offline and reduced projections used online (Wilson et al., 2020). That work also established scaling identities,

$$(\boldsymbol H_\alpha-\boldsymbol H_0)(\boldsymbol x)_i = (\boldsymbol D_x^2 G(\boldsymbol x,\boldsymbol z))_{ij} (\mathcal M)_{jk} (\boldsymbol H_0(\boldsymbol z))_k +O(\alpha^4),$$7

and

$$(\boldsymbol H_\alpha-\boldsymbol H_0)(\boldsymbol x)_i = (\boldsymbol D_x^2 G(\boldsymbol x,\boldsymbol z))_{ij} (\mathcal M)_{jk} (\boldsymbol H_0(\boldsymbol z))_k +O(\alpha^4),$$8

which make size and conductivity sweeps inexpensive once one reference signature is known (Wilson et al., 2020).

For classification, the spectral MPT is treated as a function of frequency,

$$(\boldsymbol H_\alpha-\boldsymbol H_0)(\boldsymbol x)_i = (\boldsymbol D_x^2 G(\boldsymbol x,\boldsymbol z))_{ij} (\mathcal M)_{jk} (\boldsymbol H_0(\boldsymbol z))_k +O(\alpha^4),$$9

and machine-learning pipelines typically use rotation-invariant features derived from $B_\alpha=\alpha B+\boldsymbol z$0 and $B_\alpha=\alpha B+\boldsymbol z$1, rather than raw coefficients (Wilson et al., 2021). The principal invariants

$B_\alpha=\alpha B+\boldsymbol z$2

are preferred because they are orientation invariant, avoid eigenvalue-ordering ambiguity, and are numerically smoother than root-finding-based features (Wilson et al., 2021).

5. Distinct meanings of “rank-2 magnetic polarizability tensor”

The term is not used uniformly across all branches of electromagnetic and particle theory. Three distinct meanings appear in the cited literature.

Context	Rank-2 object	Role
Eddy-current metal detection	$B_\alpha=\alpha B+\boldsymbol z$3	Complex symmetric $B_\alpha=\alpha B+\boldsymbol z$4 descriptor of a small conducting permeable object
Relativistic polarization theory	$B_\alpha=\alpha B+\boldsymbol z$5	Polarization-magnetization tensor, not a standalone magnetic polarizability tensor
Spin- or OAM-dependent particle response	$B_\alpha=\alpha B+\boldsymbol z$6 via quadratic operators	Tensor magnetic polarizability in spin or orbital degrees of freedom

In relativistic field theory for structural microparticles, the central second-rank object is the polarization tensor

$B_\alpha=\alpha B+\boldsymbol z$7

defined through

$B_\alpha=\alpha B+\boldsymbol z$8

and reconstructed as

$B_\alpha=\alpha B+\boldsymbol z$9

In that paper, magnetic polarizability is not a rank-2 tensor; it is a scalar coefficient $\mathcal M$0 entering the constitutive relation

$\mathcal M$1

and the interaction Lagrangian

$\mathcal M$2

The paper is explicit that there is no anisotropic tensor $\mathcal M$3 there (Maksimenko et al., 2021).

A different meaning appears in the deuteron literature, where tensor magnetic polarizability is a genuine rank-2 spin-dependent response of a spin-1 system. There the interaction is

$\mathcal M$4

and $\mathcal M$5 couples to the quadratic spin operator $\mathcal M$6 in the frozen-spin storage-ring geometry (Baryshevsky et al., 2010). This is tensorial in spin space rather than a $\mathcal M$7 spatial object characterizing a hidden conductor.

A third variant arises for twisted electrons, where an orbital, rather than spin, tensor magnetic polarizability appears through

$\mathcal M$8

Here the rank-2 structure is carried by bilinears in intrinsic orbital angular momentum, not by the eddy-current object descriptor $\mathcal M$9 (Silenko et al., 2019).

These distinctions are essential because the same phrase can refer either to a spatial object-characterization tensor, a relativistic polarization tensor, or a quadratic spin/OAM interaction coefficient.

6. Applications, invariants, and limitations

In metal detection and hidden-object identification, the rank-2 MPT is used as the sole descriptor layer between the electromagnetic forward problem and the classifier. The tensor’s spectral signature can be obtained either from simulated eddy-current solutions or, in practice, from induced-voltage measurements over a frequency band followed by inversion to the tensor coefficients (Wilson et al., 2021). Because the tensor is position independent and transforms covariantly under rotation, it is well suited to dictionary-based identification once rotationally invariant features are extracted.

The practical value of the tensor follows from several structural properties. First, it is compact: six complex coefficients in the generic case, fewer under symmetry (Wilson et al., 2021). Second, it is interpretable: low-frequency limits connect to the Pólya–Szegő tensor, while high-conductivity limits approach $3\times 3$0 for simply connected objects and exteriors (Ledger et al., 2015). Third, it is computationally accessible: explicit coefficient formulas, finite element transmission problems, reduced-order surrogates, and scaling laws all exist in the literature (Wilson et al., 2020).

The limitations are equally clear. The rank-2 model assumes a small-object regime and, at leading order, effectively a background field that is approximately uniform over the object (Ledger et al., 2017). When the background field varies significantly across the target, the classical tensor may be insufficient and higher-order GMPTs become relevant (Ledger et al., 2022). Measurement noise, nonuniform background fields, capacitive coupling, soil and background effects, parasitic voltages, and filtering also perturb recovered coefficients; reported measurement errors for MPT coefficients are typically about $3\times 3$1 to $3\times 3$2 in the classification study (Wilson et al., 2021).

A common misconception is that the rank-2 MPT is a universal notion independent of modeling regime. The literature does not support that view. In the dominant eddy-current usage, it is a complex symmetric $3\times 3$3 tensor associated with small conducting permeable objects. In relativistic polarization theory, the analogous second-rank object is not itself a magnetic polarizability tensor. In nuclear and vortex-beam physics, “tensor magnetic polarizability” usually denotes a quadratic spin or OAM response coefficient rather than a spatial object-characterization tensor (Maksimenko et al., 2021).

The mature interpretation, therefore, is domain-specific. In low-frequency induction sensing, the rank-2 magnetic polarizability tensor is the leading-order, orientation-covariant, frequency-dependent object signature of a small conducting permeable target. It forms the classical entry point to a broader hierarchy of generalised tensors, and it remains the foundational descriptor for asymptotic modeling, reduced-order computation, and spectral classification of hidden metallic objects (Ledger et al., 2022).