Magnetic Dipole in a Cuboidal Superconducting Trap

Technical Summary: Magnetic Dipole in a Cuboidal Superconducting Trap

Problem Statement
The paper addresses the theoretical challenge of determining the exact trapping potential and equilibrium orientation of a point magnetic dipole confined within a closed, rectangular superconducting trap. While the image method for a dipole between two infinite parallel superconducting plates is well-established, it leaves the dipole free to drift parallel to the plates. Real-world superconducting traps, however, are fully enclosed, confining the dipole in all three translational degrees of freedom and introducing a complex orientational degree of freedom determined by the trap's geometry. The authors seek to derive an exact analytical solution for this fully enclosed geometry that satisfies the Meissner boundary condition ($B \cdot \hat{n} = 0$) on all six walls simultaneously.

Methodology
The authors employ the method of images to construct the solution. Starting from the single-plane case, they generalize the construction to a cuboidal domain $\Omega = [-a, a] \times [-b, b] \times [-c, c]$.

1. Image Lattice Construction: By reflecting the source dipole and its subsequent images across all six orthogonal walls, the authors generate an infinite lattice of image dipoles. The positions $\mathbf{r}_n$ and moments $\boldsymbol{\mu}_n$ of these images are indexed by integer triples $\mathbf{n} = (n_x, n_y, n_z)$, with the moment components alternating signs based on the parity of the indices.
2. Boundary Condition Verification: The authors prove that this specific lattice construction satisfies the Neumann boundary condition on all six walls simultaneously. This relies on a "pair identity" showing that any dipole and its mirror image across a specific wall produce a vanishing normal field component on that wall. The proof utilizes the reflection symmetry of the lattice to demonstrate that the total field cancels exactly on the boundaries.
3. Convergence Analysis: The authors address the conditional convergence of the dipole field sum in the closed 3D trap. They demonstrate that because the image moments alternate in sign such that the net magnetization of the lattice is zero, the sum converges to a shape-independent limit when evaluated symmetrically about the source.
4. Energy Formulation: The trapping potential is derived as the self-energy of the dipole interacting with the induced field from the image lattice. For a dipole centered at the origin, the orientational energy is reduced to a diagonal quadratic form, $U(\theta, \phi) = U_x \sin^2\theta \cos^2\phi + U_y \sin^2\theta \sin^2\phi + U_z \cos^2\theta$. The coefficients $U_x, U_y, U_z$ are expressed as Epstein-zeta-type lattice sums.
5. Verification: The analytical results are verified against finite-element method (FEM) solutions of the same boundary-value problem using two independent solvers (FEniCSx and scikit-fem).

Key Results

• Exact Potential: The paper provides a closed-form expression for the trap potential (Eq. 10) valid for arbitrary dipole position and orientation within the cuboid.
• Orientational Flip: In the limit of an infinite rectangular tube ($c \to \infty$), the authors identify a non-monotonic dependence of the equilibrium orientation on the aspect ratio ($a/b$). Contrary to simple intuition, a centred dipole aligns with the short cross-sectional axis for a finite range of aspect ratios ($1 < a/b < 1.2249$). Beyond this threshold, the alignment flips to the long axis. This behavior arises from a competition between contributions from images on the coordinate axes (favoring the long axis) and off-axis images (favoring the short axis).
• Phase Diagram: For the fully closed trap, the equilibrium orientation is mapped as a function of the two aspect ratios ($a/b$ and $c/b$). The authors identify exactly four "triple points" where the orientational energy is degenerate ($U_x = U_y = U_z$):
  1. The cubic case ($a=b=c$).
  2. A symmetry orbit of three accidental points with aspect ratios $1 : \lambda_c : \lambda_c$, where $\lambda_c \approx 2.346182$. This value is a root of a difference of Epstein zeta functions and has no known closed algebraic form.
• Parallel-Plate Limit: The solution correctly reduces to the known result for infinite parallel plates when the side walls are sent to infinity, recovering the potential where the dipole aligns parallel to the plates.
• Numerical Agreement: The analytical image-lattice results agree with the finite-element solutions to better than 0.16% across the tested parameter space.

Significance and Claims
The paper claims to provide the first exact image-dipole potential for a dipole in a closed cuboidal superconducting trap, generalizing previous parallel-plate results to a geometry that confines all translational degrees of freedom. The primary significance lies in the rigorous proof that the constructed image lattice satisfies the Meissner boundary conditions on all six walls simultaneously, resolving potential concerns about boundary condition violations in multi-wall geometries.

The work reveals that the selection of the "easy axis" for a levitated magnet is not strictly monotonic with trap aspect ratio; instead, it exhibits a counter-intuitive preference for the short axis in specific geometric regimes. The authors provide a complete phase diagram for the orientational equilibrium, classifying the degeneracies and identifying the specific geometric ratios where transitions occur. Finally, the paper establishes that the closed-form solution supplies the full translational and orientational stiffness, enabling the calculation of all trap and librational frequencies without resorting to numerical approximations, while maintaining high precision verified against direct numerical solutions.