Magnetic Dipole Trapping Potential between Infinite Superconducting Plates

Francis J. Headley

REVIEW 2 major objections 4 minor 1 cited by

Magnetic Dipole Trapping Potential between Infinite Superconducting Plates

T0 review · 2 major / 4 minor · reviewed 2026-08-16 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read The paper derives a closed-form magnetic dipole trapping potential between two infinite superconducting plates, built from an infinite series of image dipoles and expressible through the Riemann zeta function.

desk verdict The two-plate potential in Eq. (16) is a correct and useful new closed form for the idealized point-dipole/ideal-Meissner geometry, but the paper's claim that it models real experiments is not quantified. read the letter →

arxiv 2504.18852 v4 pith:QK6E7XY6 submitted 2025-04-26 cond-mat.supr-con physics.app-ph

Francis J. Headley This is my paper · ORCID

classification cond-mat.supr-conphysics.app-ph

keywords magneticdipolesuperconductingplatesimagemethodMeissnereffecttrappingpotentiallevitatedmagnetRiemannzetafunctionresonancefrequencies

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper derives the exact analytic form of the potential felt by a point magnetic dipole sitting between two infinite, parallel, ideal superconducting plates. The derivation extends the image method to magnetism: an infinite ladder of mirror dipoles enforces the Meissner boundary condition on both plates, and summing their fields gives a closed-form potential that involves the Riemann zeta function. The resulting formula yields explicit resonance frequencies for the vertical and tilt modes of a levitated magnet and reduces to the known single-plate result as one plate recedes. It is intended as a benchmark for numerical simulations and as a building block for models of levitated magnets in precision experiments.

What carries the argument

The load-bearing object is the infinite image-dipole ladder: a one-dimensional stack of equal-strength dipoles with alternating orientation, placed at the mirror positions $z_{n+}=2n(b-a)+z_0$ and $z_{n-}=2n(b-a)+2b-z_0$ for all integers $n$. It is generated iteratively because the first mirror on one plate leaves a nonzero normal field at the other plate, and each cancelling image creates a new residual that must itself be cancelled. Summing the scalar potentials of all images and evaluating the self-energy $U=-\tfrac12\boldsymbol{\mu}\cdot\mathbf{B}_I$ produces Eq. (16), with the Riemann zeta function $\zeta(3)$ (the infinite sum $\sum_{n=1}^\infty n^{-3}$) emerging from the $1/|b-a|^3$ spacing of the ladder.

What would settle it

Measure the vertical and tilt resonance frequencies of a small spherical magnet levitated between two planar superconducting surfaces whose lateral size is much larger than their separation, with the superconductor well below its critical temperature; if the frequencies do not approach Eq. (19) as the plates are made wider and the penetration depth smaller, the image-dipole potential is missing a physical contribution.

Watch

Extended reading notes

Core claim

The central claim is Eq. (16), which states that the trapping potential of a dipole $\boldsymbol{\mu}$ between plates at $z=a$ and $z=b$ is $U = \frac{\mu_0\mu^2}{8\pi}\left( \frac{\cos^2\beta_0 - 2\sin^2\beta_0}{4|b-a|^3}\zeta(3) + \sum_{n\in\mathbb{Z}} \frac{\sin^2\beta_0 + 1}{|z_0 - z_{n-}|^3} \right)$, where the image positions are $z_{n+}=2n(b-a)+z_0$ and $z_{n-}=2n(b-a)+2b-z_0$. The paper argues this is the unique solution of the ideal Meissner boundary-value problem because each added image cancels the residual normal field left by the previous ones, and the series converges everywhere except at the plates. For symmetric plates this gives the resonance frequencies (19), and the rotational potential (20) is a double well whose minima place the dipole parallel to the plates.

Load-bearing premise

The derivation is exact only for point dipoles and perfect superconductors with zero London penetration depth; real magnets have finite size and real superconductors have a finite penetration depth.

Editorial extensions

If this is right

Any numerical solver for magnetostatic trapping in geometries that approach two parallel plates can be checked against a closed-form benchmark instead of relying only on convergence studies.
The resonance frequencies in Eq. (19) let an experimenter predict the vertical and tilt oscillation frequencies of a levitated magnet directly from its magnetic moment, mass, radius, and plate separation.
The double-well rotational potential in Eq. (20) implies the magnet has two stable orientations parallel to the plates, with an energy barrier between them that can be tuned by the plate separation.
In the limit of infinite plate separation, the two-plate potential reproduces the known single-plate result, providing a consistency check on the image construction.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

Inference: a finite London penetration depth should soften the image ladder, and the leading correction to Eq. (16) would plausibly scale as $\lambda/(b-a)$; a numerical study varying $\lambda$ could quantify when the ideal-Meissner formula stops being accurate.
Inference: the simulated tilt-mode discrepancy of about 5% is consistent with side-wall contributions in the finite cylinder used for the simulation; extrapolating finite-element results to infinitely wide plates would test whether the analytic tilt frequency is the true infinite-plate limit.
Inference: the same mirror-ladder method might produce closed-form potentials for a dipole near two perpendicular superconducting planes or inside a rectangular superconducting box, where finite sums over images would again condense into zeta-type series.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

The real new thing here is Eq. (16), the closed-form potential for a point dipole between two ideal superconducting plates, with the zeta(3) term and the resulting resonance frequencies. I checked the derivation against the known single-plate limit and the image construction; the image positions close, the orientations alternate correctly, and Eq. (16) reduces to Eq. (8) as b goes to infinity. The reader's re-derivation found no errors, and I agree. The paper does this cleanly, with a self-contained derivation and a plausible FEM sanity check.

The soft spots are not in the math but in the packaging. First, the derivation skips the summation that produces zeta(3) and the second derivative evaluations for Eq. (19). Those are routine but should be shown; a referee would want them. Second, and more substantively, the paper claims the result 'should serve as a good approximation of the trapping potential observed in real world experiments,' but it never quantifies that. The model is a point dipole with zero London penetration depth. The Numerics section uses r = 0.24 mm and b = 1 mm, so r/b is about 0.24; at that ratio higher multipoles of a real magnet are not obviously negligible. The FEM benchmark also models a point dipole, so it can only validate the image method mathematics, not the finite-size approximation. London penetration depth is also left out. These idealizations are acknowledged, but their quantitative impact is never estimated, so the experimental claim is more wish than result.

The comparison with FEM is honestly reported, including the 5% beta-mode discrepancy, and the explanation (finite cylinder walls) is reasonable. But the simulation geometry is a finite cylinder, not infinite plates, so the agreement is a useful sanity check, not a rigorous validation of the physical model.

Bottom line: the central formula is correct and genuinely new, though it is a direct extension of a standard method. It will serve as a benchmark for numerical methods in idealized geometries. The paper deserves peer review, but it needs revision: show the skipped derivations, add more detail or code for the FEM, and either quantify or sharply circumscribe the point-dipole/ideal-Meissner approximation before claiming experimental relevance.

Referee Report

2 major / 4 minor

Summary. The manuscript derives the magnetostatic potential of a point magnetic dipole between two infinite parallel superconducting plates in the ideal Meissner state, using an infinite series of image dipoles. The main result is Eq. (16), which is used to obtain the rotational potential (20) and the resonance frequencies (19). The analytic result is compared with FEniCS finite-element simulations in a cylindrical geometry.

Significance. If the derivation is correct, Eq. (16) provides a useful exact benchmark for numerical solvers and a starting point for levitated-magnet dynamics. The derivation is self-contained, uses no fitted parameters, and the image construction closes consistently. The FEM comparison is an independent check rather than an input, which strengthens the paper.

major comments (2)

[Introduction and §Numerics] The claim that Eq. (16) 'should serve as a good approximation of the trapping potential observed in real world experiments' is not supported by the evidence in the manuscript. The FEM benchmark in the Numerics section treats the levitated particle as a point dipole, so it cannot quantify the error from the finite size of a real magnet; with r = 2.4e-4 m and b = 1e-3 m, r/b ≈ 0.24, so higher multipole moments are not obviously negligible. The authors should either quantify finite-size corrections (for example by integrating the dipole field over the magnet volume) or explicitly restrict the stated applicability to the idealized point-dipole, ideal-Meissner model.
[§Numerics, Fig. 2 and Table I] The FEM simulation uses a hollow cylinder of radius 25 mm and height 2 mm, while the analytic model assumes infinite plates. The text asserts that increasing the cylinder radius systematically reduces the discrepancy, but no data are shown. Because the paper proposes the analytic result as a benchmark for numerical methods, the comparison should include a convergence study with respect to cylinder radius and mesh resolution to confirm that the 4.78% discrepancy in fβ is attributable to the finite geometry rather than numerical error.

minor comments (4)

[Between Eqs. (15) and (16)] The evaluation of the infinite sums leading to the ζ(3) coefficient is not shown; please include the sums or provide an appendix with the relevant identities.
[Eq. (19)] The second derivatives of the potential that give the resonance frequencies are not shown; please include the calculation or a reference.
[Fig. 2] The figure caption and labels are garbled in the manuscript (e.g., '□0.05' and 'Interpolated'), and the units and the meaning of the fitted curves should be stated clearly.
[References] The single-plate result Eq. (8) is cited to [9] but also appears in [14]; citing [14] explicitly at Eq. (8) would be helpful.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: Eq. (16) is derived from the image-dipole construction and the Meissner boundary condition, with no fitted parameter or self-citation serving as a load-bearing input.

full rationale

The derivation is self-contained given standard magnetostatics and the ideal-Meissner assumption. The single-plate case is re-derived by imposing n·B=0 at the plate, choosing the mirrored image position z1=2a-z0 and orientation β1=-β0, and invoking the Laplace uniqueness theorem from the external textbook Jackson [8]. The two-plate result is obtained by iterating the same reflection condition: the positions in Eq. (14) are the closure of reflections across the two plates, and the alternating orientations are chosen to cancel the residual normal fields shown in Eqs. (12)-(13). Potential (16) is then evaluated from U=-1/2 μ·B_I, with B_I=-μ0∇Φ_I; that energy formula is imported from the external image-method literature [9], not from the present authors. No parameter is fitted to produce Eq. (16): μ, |b-a|, and z0 are inputs, and the ζ(3) coefficient follows from the closed image sum. The resonance frequencies in Eq. (19) are second derivatives of the same closed-form potential, not independently fitted outputs. The FEniCS FEM simulation is an independent numerical check under the same point-dipole, ideal-Meissner assumptions; the paper explicitly notes in the Numerics section that 'it is not possible to simulate two infinite plates using FEniCS' and that the residual discrepancy is due to the finite cylinder geometry, so the benchmark is not an input to the derivation. The acknowledged idealizations (point dipole, zero London penetration depth, infinite plates) limit experimental applicability but do not make the derivation circular. Prior work [14] is cited only for the known one-plate limit and as application context; the load-bearing uniqueness and energy arguments are external or re-derived here. No circular step is present.

Assumptions & free parameters 0 free parameters · 4 assumptions · 0 invented entities

The derivation has no free parameters and no invented entities. It relies on standard magnetostatics, the ideal Meissner boundary condition, the uniqueness of the Laplace solution, and the energy formula U = -1/2 mu dot B_I taken from Ref. [9].

assumptions (4)

domain assumption The superconductor is in an ideal Meissner state with zero magnetic induction inside, so the normal component of B vanishes on each plate.
Invoked in the Introduction ('ideal Meissner state, a perfect diamagnet') to justify the boundary condition dz Phi = 0 used throughout the image construction.
standard math The region outside the plates is current-free, so the magnetic field is the gradient of a scalar potential satisfying Laplace's equation.
Used to write Phi0 and PhiI in Eqs. (1) through (6); a standard magnetostatics reduction.
standard math The image construction together with the uniqueness theorem for Laplace's equation yields the unique physical solution.
Stated after Eq. (6) and used to justify summing the infinite image series in Eq. (15).
domain assumption The energy of a permanent dipole in the induced field is U = -1/2 mu dot B_I.
Taken from Ref. [9] and used in Eq. (7) and Eq. (16); the factor 1/2 encodes the work done against the self-consistent induced field.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Magnetic Dipole Trapping Potential between Infinite Superconducting Plates." pith.science (2026). https://pith.science/paper/QK6E7XY6

@misc{pith2026250418852,
  author       = {Pith},
  title        = {Pith review of: Magnetic Dipole Trapping Potential between Infinite Superconducting Plates},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/QK6E7XY6}},
  note         = {Machine review of arXiv:2504.18852}
}

read the original abstract

We derive the exact analytic form of the potential experienced by a magnetic dipole trapped between two infinite parallel superconducting plates using the method of image dipoles, providing a benchmark for numerical methods and a foundation for studying the stability and dynamics of magnetically levitated systems in precision measurements and fundamental physics experiments.

Figures

Figures reproduced from arXiv: 2504.18852 by the authors.

**Figure 2.** FIG. 2. Finite element method (FEM) simulations of a magnetic dipole trapped in a hollow superconducting cylinder near [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

Enhancement of the effects due to the Schr\"odinger-Newton equation
quant-ph 2025-07 conditional novelty 6.0 of 10

Periodic trap frequency modulation can amplify Schrödinger-Newton deviations in the position variance by up to six orders of magnitude, enabling a possible experimental test.

Reference graph

Works this paper leans on

14 extracted references · 11 canonical work pages · cited by 1 Pith paper

[1]

M. S. Alnaes, J. Blechta, J. Hake, A. Johansson, B. Kehlet, A. Logg, C. N. Richardson, J. Ring, M. E. Rognes, and G. N. Wells , The FEniCS project version 1.5 , Archive of Numerical Software, 3 (2015)

work page 2015
[2]

D. W. P. Amaral, D. G. Uitenbroek, T. H. Oost- erkamp, and C. D. Tunnell , First Search for Ultra- light Dark Matter Using a Magnetically Levitated Parti- cle, (2024)

work page 2024
[3]

S. Bose, I. Fuentes, A. A. Geraci, S. M. Khan, S. Qvarfort, M. Rademacher, M. Rashid, H. Ul- bricht, and C. C. W anjura , Massive quantum sys- tems as interfaces of quantum mechanics and gravity , arXiv:2311.09218, (2023)

arXiv 2023
[4]

Carney, S

D. Carney, S. Ghosh, G. Krnjaic, and J. M. Tay- lor, Proposal for gravitational direct detection of dark matter, Phys. Rev. D, 102 (2020), p. 072003

work page 2020
[5]

T. M. Fuchs, D. Uitenbroek, J. Plugge, N. van Halteren, A. Vinante, H. Ulbricht, and T. H. Oosterkamp, Magnetic Zeppelin: Detection of gravita- tional drive in the Hz regime , arXiv:2303.03545, (2023)

arXiv 2023
[6]

Gonzalez-Ballestero, M

C. Gonzalez-Ballestero, M. Aspelmeyer, L. Novotny, R. Quidant, and O. Romero-Isart , Levitodynamics: Levitation and control of microscopic objects in vacuum , Science, 374 (2021), p. eabg3027

work page 2021
[7]

F. J. Headley, A. Belenchia, M. Paternostro, and D. Braun , Quantum Metrology of Newton ’s Con- stant with Levitated Mechanical Systems , Mar. 2025. arXiv:2503.16215 [quant-ph]

arXiv 2025
[8]

J. D. Jackson , Classical Electrodynamics, Wiley, New York, 3rd ed., 1998

work page 1998

Show all 14 references

[9]

Q.-G. Lin , Theoretical development of the image method for a general magnetic source in the presence of a su- perconducting sphere or a long superconducting cylinder , Physical Review B, 74 (2006), p. 024510

2006
[10]

Prat-Camps, C

J. Prat-Camps, C. Teo, C. Rusconi, W. Wieczorek, and O. Romero-Isart, Ultrasensitive Inertial and Force Sensors with Diamagnetically Levitated Magnets , Physi- cal Review Applied, 8 (2017), p. 034002

2017
[11]

Timberlake, G

C. Timberlake, G. Gasbarri, A. Vinante, A. Set- ter, and H. Ulbricht , Acceleration sensing with mag- netically levitated oscillators above a superconductor , Ap- plied Physics Letters, 115 (2019), p. 224101

2019
[12]

Tinkham , Introduction to Superconductivity , McGraw-Hill, New York, 2nd ed., 1996

M. Tinkham , Introduction to Superconductivity , McGraw-Hill, New York, 2nd ed., 1996

1996
[13]

Vinante, P

A. Vinante, P. F alferi, G. Gasbarri, A. Set- ter, C. Timberlake, and H. Ulbricht , Ultralow mechanical damping with Meissner-levitated ferromag- netic microparticles, Physical Review Applied, 13 (2020), p. 064027

2020
[14]

Vinante, C

A. Vinante, C. Timberlake, and H. Ulbricht , Lev- itated Micromagnets in Superconducting Traps: A New Platform for Tabletop Fundamental Physics Experiments , Entropy, 24 (2022), p. 1642

2022

Pith tools

Reviewed August 16, 2026 · model on record in the stance chip above.

[1] [1]

M. S. Alnaes, J. Blechta, J. Hake, A. Johansson, B. Kehlet, A. Logg, C. N. Richardson, J. Ring, M. E. Rognes, and G. N. Wells , The FEniCS project version 1.5 , Archive of Numerical Software, 3 (2015)

work page 2015

[2] [2]

D. W. P. Amaral, D. G. Uitenbroek, T. H. Oost- erkamp, and C. D. Tunnell , First Search for Ultra- light Dark Matter Using a Magnetically Levitated Parti- cle, (2024)

work page 2024

[3] [3]

S. Bose, I. Fuentes, A. A. Geraci, S. M. Khan, S. Qvarfort, M. Rademacher, M. Rashid, H. Ul- bricht, and C. C. W anjura , Massive quantum sys- tems as interfaces of quantum mechanics and gravity , arXiv:2311.09218, (2023)

arXiv 2023

[4] [4]

Carney, S

D. Carney, S. Ghosh, G. Krnjaic, and J. M. Tay- lor, Proposal for gravitational direct detection of dark matter, Phys. Rev. D, 102 (2020), p. 072003

work page 2020

[5] [5]

T. M. Fuchs, D. Uitenbroek, J. Plugge, N. van Halteren, A. Vinante, H. Ulbricht, and T. H. Oosterkamp, Magnetic Zeppelin: Detection of gravita- tional drive in the Hz regime , arXiv:2303.03545, (2023)

arXiv 2023

[6] [6]

Gonzalez-Ballestero, M

C. Gonzalez-Ballestero, M. Aspelmeyer, L. Novotny, R. Quidant, and O. Romero-Isart , Levitodynamics: Levitation and control of microscopic objects in vacuum , Science, 374 (2021), p. eabg3027

work page 2021

[7] [7]

F. J. Headley, A. Belenchia, M. Paternostro, and D. Braun , Quantum Metrology of Newton ’s Con- stant with Levitated Mechanical Systems , Mar. 2025. arXiv:2503.16215 [quant-ph]

arXiv 2025

[8] [8]

J. D. Jackson , Classical Electrodynamics, Wiley, New York, 3rd ed., 1998

work page 1998