Pith. sign in

REVIEW 4 major objections 5 minor 26 references

Stability of rotating magnetic levitation

T0 review · 4 major / 5 minor · reviewed 2026-08-06 · deepseek-v4-flash

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

Pith's one-line read Damping is essential for stable rotating magnetic levitation, and the stable region in rotation speed and damping is bounded.

desk verdict A useful six-DOF extension of the rotating-levitation program, but the printed stability analysis drops the gyroscopic coupling that stabilizes a spinning top, so the phase diagrams need a corrected derivation. read the letter →

arxiv 2507.07478 v1 pith:Q3WPVJBG submitted 2025-07-10 physics.class-ph physics.app-ph

Mingjun Fan , Jinyu Chen , Yongquan Ji , Long Li , Chichuan Ma , Yu-Han Ma This is my paper · ORCID
classification physics.class-phphysics.app-ph
keywords magneticlevitationrotatingmagneteddy-currentdampingstabilityanalysissix-degree-of-freedommodelphasediagramprecessionalmotioncopperplate
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 argues that a magnet floating above a rotating magnet and copper plate is stable only inside a bounded region in rotation speed and damping coefficients, not merely at high speed or with strong damping. The authors build a six-degree-of-freedom model with eddy-current damping from the copper plate, find the equilibrium, linearize around it, and compute eigenvalues to map stability. The central result is that damping is essential—without the copper plate's drag the floater does not levitate stably—and that too much damping also destroys stability. The predicted shape of the stable region and the oscillation frequencies near the rotor speed match the authors' experiments qualitatively. This matters because it gives a stability criterion for a contactless levitation scheme that does not require diamagnets or superconductors.

What carries the argument

The carrying object is the linearized five-degree-of-freedom perturbation system $M\ddot{p}+C\dot{p}+Kp=0$ in the rotor's co-rotating frame, where $C$ carries the copper-plate damping ($\alpha$ horizontally, $\beta$ vertically, zero rotational damping) together with Coriolis terms, and $K$ collects magnetic, centrifugal and damping-force gradients. Stability is read off the eigenvalues of this system: all eigenvalues with negative real part define stable levitation. The mechanism that makes damping load-bearing is that the non-conservative drag enters $C$ off-diagonally with the rotation terms, reshaping the eigenvalue spectrum; a copper plate, or a viscous fluid, supplies this damping, while a wooden plate does not.

What would settle it

Measure the damping force on a floater magnet moving at constant velocity in several directions, including rotation about its own axis, at the same heights used in levitation, and compare with the linear $\alpha(h)v$ and $\beta(h)v$ laws; if the horizontal drag differs between sliding on a tilted plate and moving over the plate in levitation geometry, or if rotating the floater produces measurable drag torque, the eigenvalue stability diagram would not apply.

Watch

Extended reading notes

Core claim

The central claim is that stable magnetic levitation of a dipole floater above a horizontal rotating dipole rotor, with a copper plate between them, occurs only within clearly delimited ranges of the rotation speed $\omega_r$ and the horizontal and vertical damping coefficients $\alpha$ and $\beta$. In the co-rotating frame, the balance point is found by balancing magnetic forces with centrifugal and damping forces; then small perturbations around that point are studied through a linearized matrix equation $M\ddot{p}+C\dot{p}+Kp=0$. Stability is decided numerically from the eigenvalues: levitation is stable whenever all real parts are negative. The resulting phase diagrams show a bounded stable region in the $(\alpha,\beta)$ plane at fixed $\omega_r$, so damping is necessary but excessive damping is harmful, and the minimum $\omega_r$ rises as the floater's moment of inertia shrinks. The paper states that these predictions are in qualitative agreement with experiments, including the observed frequency peaks near $\omega_r$ during stable levitation.

Load-bearing premise

The copper plate is represented as a purely linear translational drag with two position-dependent coefficients $\alpha(h)$ and $\beta(h)$, measured in steady sliding or in simulation, while rotational damping is ignored and the plate's distortion of the magnetic field is neglected; if the real eddy-current force during levitation is nonlinear, direction-dependent, or exerts spin torque, the computed stability boundaries may not describe the actual system.

Editorial extensions

If this is right

  • For a given rotor, floater, and copper plate, there is a concrete lower bound on rotation speed and a finite window of damping coefficients; engineers can place the floater at a height $h$ that lands inside this window at the selected $\omega_r$.
  • Floaters with larger moments of inertia $I_{12}$ can levitate at lower rotation speeds and over a wider damping window, matching the experimental observation that smaller floater magnets need higher $\omega_r$.
  • Increasing the rotor's magnetic moment does not always help: it can shrink the stable region in the $(\alpha,\beta)$ plane, so an optimal rotor strength exists rather than 'stronger is better.'
  • Overdamping is as fatal as underdamping: pushing the floater too close to the plate raises $\alpha$ and $\beta$ beyond the stable window and destroys levitation even at high $\omega_r$.
  • The natural oscillation frequencies of the floater near the stability boundary cluster close to the rotor frequency $\omega_r$, and the FFT peaks of the measured $x(t)$ align with the theoretical eigenfrequencies.

Reading between the lines

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

  • If the linear drag model is replaced by a measured nonlinear eddy-current force, the stability window will likely shift or become asymmetric; an impulse-decay test of the floater at several heights would reveal whether the decay is exponential as predicted.
  • The same phase-diagram logic should apply to other rotating-dipole systems in viscous or conductive surroundings, meaning controlled fluid viscosity could act as a tunable stabilizer for contactless manipulation.
  • Because the paper sets rotational damping to zero, spin-damping effects on the nutation mode remain untested; spinning the floater in a viscous fluid would expose whether an additional stability boundary exists.
  • A practical design rule follows: measure $\alpha(h)$ and $\beta(h)$ for a given plate, then choose the plate–rotor distance so the equilibrium height falls inside the stable window at the intended $\omega_r$.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

4 major / 5 minor

Summary. The manuscript studies stable magnetic levitation of a permanent-magnet floater above a rotating dipole magnet with a stationary copper plate. The authors formulate a Lagrangian in the co-rotating frame, eliminate the cyclic spin coordinate via Routh's procedure, derive equilibrium positions, linearize the five-degree-of-freedom dynamics, and compute stability boundaries numerically in the parameter space (ωr, α, β). They measure horizontal and vertical damping coefficients of the copper plate and compare predicted stability regions and oscillation frequencies with high-speed-camera observations, reporting qualitative agreement. The central claim is that stable levitation occurs only in bounded ranges of rotation speed and damping coefficients, with damping essential for stability.

Significance. The topic is timely and the paper is ambitious: it attempts the first full linearized stability analysis of this recently discovered levitation configuration, including independently measured damping coefficients rather than fitting the stability boundary. The experimental effort, including measurement of α(h) and β(h) and FFT identification of natural frequencies, is a genuine strength. If the derivation were correct, the resulting phase diagrams would be a useful reference for designing levitation experiments and educational demonstrations. However, the central linearized equations contain a clear internal omission that affects the computed stability boundaries, so the main theoretical result cannot be accepted in its present form.

major comments (4)
  1. [Sec. III (C matrix); Appendix A (Eq. A1)] The linearized damping matrix C printed in Sec. III has an identically zero 2x2 rotational block. This is inconsistent with the kinematic expression in Eq. (A1): expanding the rotational kinetic energy 1/2 I12(ω1^2 + ω2^2) to second order in the small angles and their derivatives gives the cross terms -I12 ωr φy φ̇x + I12 ωr φx φ̇y. The Euler-Lagrange equations for φx and φy therefore contain gyroscopic terms -2 I12 ωr φ̇y and +2 I12 ωr φ̇x, i.e., C45 = -2ωr and C54 = +2ωr after dividing by I12. The translational block of the same matrix correctly includes the analogous Coriolis terms ±2ωr, so the omission is not a matter of convention. Because gyroscopic coupling is the standard stabilizing mechanism for nutational motion of a spinning body, omitting it changes the eigenvalues and thus the stability boundaries in Figs. 3 and 4. Please re-derive the linearized equations from the Lagrangian and recompute all stability diagrams with the corrected C matrix.
  2. [Sec. III (K matrix)] The stiffness matrix K is presented without a derivation, and several entries contain the damping coefficient α (for example K12 = -(αωr/m - 15μY/(mZ^6)) and K21 = +(αωr/m + 15μY/(mZ^6))). These entries are not obviously stiffnesses; they originate from the position-dependent terms in the damping law in Eq. (4), -α(ẋ - ωr y) and -α(ωr x + ẏ), when written in the co-rotating frame. Since the eigenvalue problem Mλ^2 + Cλ + K = 0 is the entire basis for the stability diagrams, please show the linearization steps explicitly so the reader can verify both C and K, including the correct gyroscopic terms.
  3. [Sec. IV, Fig. 6] The theoretical stability regions are computed in (α, β, ωr) space, while the experimental boundary in Fig. 6 is plotted as a curve in (h, ωr). The paper claims qualitative agreement but does not overlay the theoretical boundary on the experimental data. Because α(h) and β(h) are measured (Appendix B), the theoretical stability region can be mapped to (h, ωr) for the actual apparatus; without such a comparison, the agreement between theory and experiment is asserted rather than demonstrated.
  4. [Appendix B] The damping model assumes scalar, purely translational drag coefficients α(h) and β(h) with no rotational damping, and neglects distortion of the magnetic field by the copper plate. The horizontal coefficient is measured in steady sliding contact (with paper spacers) and the vertical coefficient is obtained from COMSOL, not from the levitating motion. Because the stability boundaries in Sec. III are sensitive to the functional form of the damping, please validate that these coefficients describe oscillatory motion at the levitation heights, for example by comparing measured decay rates of the floater oscillations with the eigenvalues of the linearized model, or by quantifying the uncertainty this introduces in the boundary curves.
minor comments (5)
  1. [Sec. IIA] The phrase "we retain all six degrees of freedom" is misleading because the cyclic spin ψ is eliminated via Routh reduction, leaving five second-order equations; please rephrase.
  2. [Sec. IIB, Fig. 4] There are typos in the text, e.g., "mathb f mf" should be "mf", and in Fig. 4 the two moments are both labeled mf; presumably the second is the floater moment mr.
  3. [Sec. IVB, Fig. 3] The FFT peak comparison in Fig. 3(a) would be more informative with error bars or the number of trials; as printed, single red points are hard to evaluate.
  4. [Eq. (6)] The scaling law in Eq. (6) is presented without a derivation; please state the approximations under which X is negative and the exponents are obtained.
  5. [References] Reference [3] duplicates reference [2]; also, Ref. [18] has a formatting error in the volume/page string.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the stability analysis is self-contained, with damping coefficients measured independently and no load-bearing self-citations.

full rationale

The paper derives the stability conditions from a stated Lagrangian and Routh-reduced equations of motion, then substitutes independently measured damping coefficients α(h) and β(h) obtained from separate friction–velocity experiments and COMSOL simulations (Appendix B). These damping parameters are not fitted to the target stability boundaries or to the observed stable/unstable regions; they are pre-existing inputs to the linearized eigenvalue analysis in Sec. III. The observed synchronous precession is used as an empirical justification for adopting the co-rotating frame, not as a derived prediction, and the central claim that damping is essential for stable levitation follows from the eigenvalue analysis rather than being imposed. The paper cites prior experimental and theoretical work (e.g., Refs. [17, 18]) but does not rely on self-citations, and no uniqueness theorem or imported ansatz from the authors' own prior work is load-bearing. The skeptic's concern about the omitted gyroscopic coupling in the 2×2 rotational block of C is an internal modeling/correctness issue, not a circularity: it does not make the result equivalent to its inputs. Overall, the derivation chain is self-contained and externally benchmarked against qualitative experimental observations, so no circular step is identified.

Assumptions & free parameters 2 free parameters · 6 assumptions · 0 invented entities

No parameter is fitted to the levitation stability result itself. The damping coefficients are calibrated from independent force-velocity measurements or COMSOL. The main caveats are modeling assumptions about point dipoles, synchronous precession, and scalar linear damping, none of which are directly verified under levitation conditions.

free parameters (2)
  • horizontal damping coefficient α(h) = function of height h, from linear fits to friction-velocity curves
    Measured by sliding the floater on a tilted copper plate (Fig. 8). It is an independently calibrated input, not tuned to the stability boundary, but it is fitted to data.
  • vertical damping coefficient β(h) = function of height h, from COMSOL simulation and linear fits
    Computed in COMSOL because direct vertical drag measurement was not available (Fig. 8). This is a simulation-calibrated input to the model.
assumptions (6)
  • domain assumption The floater and rotor are treated as point magnetic dipoles interacting through the dipole-dipole potential.
    Used in Eq. 2 for the magnetic potential; valid only when separation is much larger than magnet size, which is not explicitly verified for the experimental geometry.
  • domain assumption During stable levitation the floater precesses synchronously with the rotor with zero phase lag, so the rotor is static in the co-rotating frame.
    Stated as experimental observation in Sec. II, item ii. If phase lag is nonzero, the time-dependent magnetic field must be retained.
  • domain assumption Eddy-current damping from the copper plate is a linear translational drag with coefficients α and β, with rotational damping neglected and plate-induced field distortion below 5%.
    Appendix B describes this simplification. It is load-bearing because damping is essential for stability in the model.
  • domain assumption The rotor angular velocity is externally maintained and unaffected by the copper plate or the floater.
    Appendix B reports measured field variation below 5%, which is used to justify neglecting back-action on the rotor.
  • standard math Perturbations around equilibrium are small enough that higher-order terms in X, Y, φX, φY can be neglected and linearization is valid.
    Introduced in Sec. II B and used throughout Sec. III to justify the linearized equation.
  • standard math The cyclic coordinate ψ can be eliminated via the Routh transform, with conserved m3 set to I3ωr from stationary initial conditions.
    Sec. II A uses this to reduce the equations of motion before the stability analysis.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Stability of rotating magnetic levitation." pith.science (2026). https://pith.science/paper/Q3WPVJBG

@misc{pith2026250707478,
  author       = {Pith},
  title        = {Pith review of: Stability of rotating magnetic levitation},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/Q3WPVJBG}},
  note         = {Machine review of arXiv:2507.07478}
}
read the original abstract

Dynamical magnetic levitation has attracted broad interest in the realm of physics and engineering. The stability analysis of such system is of great significance for practical applications. In this work, we investigate the stable magnetic levitation of a floater magnet above a rotating magnet and copper board system. The conditions for stable levitation are analyzed through both theoretical modeling and experimental observation. This study focuses on the interplay between magnetic forces, damping effects from the copper board, and rotational dynamics. We derive the equilibrium conditions, perform stability analysis, and present phase diagrams in parametric spaces of rotation speed and damping coefficients. The theoretical predictions show qualitative agreement with experimental results, particularly in demonstrating how damping is essential for stable levitation and how the stability region depends on the geometric and magnetic parameters of the system.

Figures

Figures reproduced from arXiv: 2507.07478 by the authors.

Figure 1
Figure 1. FIG. 1. Sketch of the system [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. Stable horizontal position [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3. Imaginary and real part of eigenvalues of linearized equation [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figures from the paper (6 more)
Figure 6
Figure 6. Figure 6: FIG. 6. Experimental Stability Diagram. The dots represent col [PITH_FULL_IMAGE:figures/full_fig_p005_6.png]
Figure 4
Figure 4. Figure 4: FIG. 4. Stability diagram with parametric space [PITH_FULL_IMAGE:figures/full_fig_p005_4.png]
Figure 5
Figure 5. Figure 5: FIG. 5. (a) Side view of the system’s configuration. The parallel [PITH_FULL_IMAGE:figures/full_fig_p005_5.png]
Figure 7
Figure 7. Figure 7: FIG. 7. Trajectories of the center of mass from top views. [PITH_FULL_IMAGE:figures/full_fig_p006_7.png]
Figure 8
Figure 8. Figure 8: FIG. 8. (a) Sketch of measurement of horizontal damping coefficient (b) Sketch of measurement of vertical damping coefficient. (c) Horizontal [PITH_FULL_IMAGE:figures/full_fig_p007_8.png]
Figure 9
Figure 9. Figure 9: FIG. 9 [PITH_FULL_IMAGE:figures/full_fig_p007_9.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

26 extracted references · 25 canonical work pages

  1. [1]

    Because the floater slides in uniform velocity, the friction can be calculated by force balance

    Horizontal Damping Coefficient In order to measure the damping coefficient in horizontal direction, we slant the copper board, place the floater magnet on top of it, and let the floater magnet to naturally slide down. Because the floater slides in uniform velocity, the friction can be calculated by force balance. Adjusting the slant angle, the constant ve...

    work page

  2. [2]

    One method involves moving the floater vertically at constant speed and measuring the re- sulting damping force

    Vertical Damping Coefficient Unlike the horizontal case, measuring the vertical damping coefficient is more challenging. One method involves moving the floater vertically at constant speed and measuring the re- sulting damping force. However, this requires precise control and sensing equipment beyond our current capabilities. As an alternative, we impleme...

    work page

  3. [3]

    Turker and A

    E. Turker and A. Arslan-Yildiz, Recent advances in magnetic levitation: a biological approach from diagnostics to tissue en- gineering, ACS Biomaterials Science & Engineering 4, 787 (2018)

    work page 2018

  4. [4]

    S. Ge, A. Nemiroski, K. A. Mirica, C. R. Mace, J. W. Hen- nek, A. A. Kumar, and G. M. Whitesides, Magnetic levitation in chemistry, materials science, and biochemistry, Angewandte Chemie International Edition 59, 17810 (2020)

    work page 2020

  5. [5]

    S. Ge, A. Nemiroski, K. A. Mirica, C. R. Mace, J. W. Hennek, A. A. Kumar, and G. M. Whitesides, Magnetic levitation in chemistry, materials science, and biochemistry, Angewandte Chemie International Edition 59, 17810 (2020), https://onlinelibrary.wiley.com/doi/pdf/10.1002/anie.201903391

    doi:10.1002/anie.201903391 2020

  6. [6]

    M. D. Simon, L. O. Heflinger, and A. K. Geim, Diamagnetically stabilized magnet levitation, American Journal of Physics 69, 702 (2001). 8

    work page 2001

  7. [7]

    Supreeth, S

    D. Supreeth, S. I. Bekinal, S. R. Chandranna, and M. Dod- damani, A review of superconducting magnetic bearings and their application, IEEE Transactions on Applied Superconduc- tivity 32, 1 (2022)

    work page 2022

  8. [8]

    M. D. Simon, L. O. Heflinger, and S. Ridgway, Spin stabi- lized magnetic levitation, American Journal of Physics 65, 286 (1997)

    work page 1997

Show all 26 references
  1. [9]

    A. T. Pérez and P. García-Sánchez, Dynamics of a levitron un- der a periodic magnetic forcing, American Journal of Physics 83, 133 (2015)

    2015

  2. [10]

    K. A. Baldwin, J.-B. de Fouchier, P. Atkinson, R. Hill, M. Swift, and D. Fairhurst, Magnetic levitation stabilized by streaming fluid flows, Physical Review Letters 121, 064502 (2018)

    2018

  3. [11]

    Perdriat, C

    M. Perdriat, C. Pellet-Mary, T. Copie, and G. Hétet, Planar magnetic paul traps for ferromagnetic particles, Phys. Rev. Res. 5, L032045 (2023)

    2023

  4. [12]

    O. N. Kirillov, Nonconservative Stability Problems of Modern Physics (De Gruyter, Berlin, Boston, 2021)

    2021

  5. [13]

    R. W. Farquhar and A. A. Kamel, Quasi-periodic orbits about the translunar libration point, Celestial mechanics 7, 458 (1973)

    1973

  6. [14]

    O. N. Kirillov and M. Levi, Rotating saddle trap as foucault’s pendulum, American Journal of Physics 84, 26 (2016)

    2016

  7. [15]

    Ucar, Polarity free magnetic repulsion and magnetic bound state, Symmetry 13, 442 (2021)

    H. Ucar, Polarity free magnetic repulsion and magnetic bound state, Symmetry 13, 442 (2021)

    2021

  8. [16]

    Committee, Ipt 2023 problem list, https://2023.iptnet

    I. Committee, Ipt 2023 problem list, https://2023.iptnet. info/problems/ (2023), accessed: 27 May 2025

    2023

  9. [17]

    Committee, Problems for the 37th iypt 2024, https:// www.iypt.org/problems/problems-iypt-2024/ (2024), accessed: 27 May 2025

    I. Committee, Problems for the 37th iypt 2024, https:// www.iypt.org/problems/problems-iypt-2024/ (2024), accessed: 27 May 2025

    2024

  10. [18]

    Magnetic Games, Maglev train with powerful magnets (2021), accessed: 2025-06-27

    2021

  11. [19]

    J. M. Hermansen, F. L. Durhuus, C. Frandsen, M. Beleggia, C. R. Bahl, and R. Bjørk, Magnetic levitation by rotation, Phys- ical Review Applied 20, 044036 (2023)

    2023

  12. [20]

    Le Lay, S

    G. Le Lay, S. Layani, A. Daerr, M. Berhanu, R. Dolbeault, T. Person, H. Roussille, and N. Taberlet, Magnetic levitation in the field of a rotating dipole, Physical Review E110, 045003 (2024)

    2024

  13. [21]

    D. R. Merkin, Introduction to the Theory of Stability , V ol. 24 (Springer Science & Business Media, 2012)

    2012

  14. [22]

    Bardin and A

    B. Bardin and A. Chekin, Orbital stability of planar oscilla- tions of a satellite in a circular orbit, Cosmic Research 46, 273 (2008)

    2008

  15. [23]

    Wesson, Hydromagnetic stability of tokamaks, Nuclear Fu- sion 18, 87 (1978)

    J. Wesson, Hydromagnetic stability of tokamaks, Nuclear Fu- sion 18, 87 (1978)

    1978

  16. [24]

    S. A. Balbus and J. F. Hawley, A powerful local shear instability in weakly magnetized disks. iv. nonaxisymmetric perturbations, Astrophysical Journal v. 400, p. 610-621 400, 610 (1992)

    1992

  17. [25]

    Meerkov, Principle of vibrational control: theory and ap- plications, IEEE Transactions on Automatic Control 25, 755 (1980)

    S. Meerkov, Principle of vibrational control: theory and ap- plications, IEEE Transactions on Automatic Control 25, 755 (1980)

    1980

  18. [26]

    R. M. Murray, Z. Li, and S. S. Sastry, A mathematical introduc- tion to robotic manipulation (CRC press, 2017)

    2017

Pith tools

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