{"id":"3b1758b6-fe50-443a-9ac1-3c5e126142f9","arxiv_id":"2507.07478","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"The paper derives and experimentally checks stability boundaries for a magnet levitating above a rotating magnet with copper-plate damping, showing that stable levitation requires a bounded window of spin speed and damping.","lead":"This paper studies why a small magnet can float above a spinning magnet when a copper plate sits between them, and calculates the rotation speeds and friction levels that keep the floater stable. The results could help engineers design contact-free magnetic particle handlers and levitation devices.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The linearized rotational block omits the gyroscopic coupling from the conserved spin pψ=I3ωr; the quantitative stability diagrams are therefore not justified.","rationale":"The reader correctly identifies that the empirical drag model is a weak point, and also notes that the linearized matrix appears to contain errors. The most decisive issue, however, is internal rather than empirical: the Routh reduction of a spinning symmetric top necessarily produces gyroscopic velocity coupling in the two tilt-angle equations. The paper includes the translational Coriolis terms but sets the rotational block of C to zero, and the corresponding K entries are also zero. Because the central quantitative claim is a bounded stable region in (α,β,ωr) computed from eigenvalues of this linearized system, an omitted gyroscopic term can move the stability boundaries. This does not automatically invalidate the qualitative conclusion that damping is essential, nor the experimental observations, which are useful and independently calibrated in part. But the phase diagrams and the claimed agreement with frequency data depend on a matrix that appears incomplete. A direct re-derivation and eigenvalue recomputation would settle the issue; until then, the verdict should remain conditional rather than accepting the quantitative diagrams as correct.","tokens_in":8849,"tokens_out":29853,"duration_ms":314936,"concrete_test":"Re-derive the rotational Routhian from Eq. (1) and Eq. (A1), keeping terms first order in φx, φy and their velocities with pψ=I3ωr fixed, and compute the resulting 2×2 gyroscopic block G (expected off-diagonal entries ±I3ωr/I12, or ±ωr if I3≈2I12). Add this block to the C matrix in Sec. III, recompute the eigenvalues, and re-plot the stability diagrams for the parameters of Fig. 4(a). If the stable β-interval at fixed (α,ωr) changes by more than about 10%, or if the qualitative shape of the stable region changes, the published quantitative phase diagrams are unreliable and the paper needs a corrected derivation plus re-analysis of the experiments.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The stability analysis in Sec. III reduces to the eigenvalues of M p¨ + C p˙ + K p = 0, so every entry of C and K is load-bearing. After the Routh reduction for the cyclic spin coordinate ψ, the Routhian must contain a velocity coupling between φx and φy. Using the body angular velocity in Eq. (A1) and expanding to first order gives, up to sign convention, ω_body ≈ (φ̇x − ωr φy, φ̇y, ωr + ψ̇). The rotational kinetic energy then contains −I12 ωr φy φ̇x, so the linearized equations acquire off-diagonal gyroscopic terms ±I12 ωr in the rotational block of C. The printed C matrix has zeros in the entire 2×2 rotational block, and the printed K has zeros at positions (4,5) and (5,4). For a symmetric top, this gyroscopic coupling is the standard mechanism that stabilizes nutational motion; omitting it changes which eigenvalues have Re(λ)<0 and therefore changes the stability boundaries in the (α,β,ωr) diagrams of Fig. 4. This is an internal omission in the central derivation, distinct from the empirical question of whether the copper-board damping is well modeled by scalar coefficients α(h), β(h).","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":9090,"tokens_out":15240,"duration_ms":157290,"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":[{"comment":"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.","section":"Sec. III (C matrix); Appendix A (Eq. A1)"},{"comment":"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.","section":"Sec. III (K matrix)"},{"comment":"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.","section":"Sec. IV, Fig. 6"},{"comment":"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.","section":"Appendix B"}],"minor_comments":[{"comment":"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.","section":"Sec. IIA"},{"comment":"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.","section":"Sec. IIB, Fig. 4"},{"comment":"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.","section":"Sec. IVB, Fig. 3"},{"comment":"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.","section":"Eq. (6)"},{"comment":"Reference [3] duplicates reference [2]; also, Ref. [18] has a formatting error in the volume/page string.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The paper addresses a phenomenon of considerable pedagogical and practical interest, and the experimental data appear to be taken in good faith. However, the missing gyroscopic coupling in the linearized rotational block is a fundamental error in the central derivation; all stability diagrams must be recomputed. I also recommend that the editor ask the authors to provide the explicit linearization and a direct theory-experiment comparison. If these are addressed, the paper could become a useful contribution."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This paper makes a real addition to the rotating-magnet levitation line: it keeps six degrees of freedom, computes nonzero horizontal equilibrium positions X,Y, includes eddy-current damping via measured coefficients, and predicts a bounded stable region in (α,β,ωr). The qualitative observation that both too little and too much damping kill levitation is a genuine result, and the authors earn credit for measuring α(h) and β(h) independently rather than fitting the stability boundary. That part is worth reading.\n\nBut the central linearization has a load-bearing hole. After the Routh reduction for the conserved spin coordinate, the rotational kinetic energy contains a cross term −I12 ωr φy φ̇x (up to sign). That term should produce off-diagonal gyroscopic entries ±I12 ωr in the 2×2 rotational block of C. The printed C has zeros there. This is exactly the coupling that stabilizes nutation of a symmetric top, and omitting it changes which eigenvalues have negative real parts, so the boundaries in Fig. 4 are not justified as printed. The stress-test note holds up on reading the paper.\n\nSmaller problems point the same way: the K matrix has α (a damping coefficient) sitting in stiffness entries, and the matrix is printed non-symmetric with no derivation. Appendix B's pure linear translational drag model, with rotational damping set to zero and field distortion neglected, might be acceptable if the derivation were clean, but it adds uncertainty.\n\nThe experimental section is honest but thin: no error bars, a rough boundary line, and the FFT comparison to theoretical natural frequencies is qualitative. It supports the claim that damping is essential and the stable region is bounded, but it cannot validate the specific diagrams.\n\nNet: this is a serious, readable paper with a fixable central error. The authors need to redo the Routh reduction and show the corrected rotational block. If the gyroscopic terms preserve a bounded stable region, the paper is genuinely useful for practical design rules. I would send it to peer review, but the referee must insist on the corrected derivation before the quantitative results are accepted.","headline":"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.","tokens_in":9602,"tokens_out":5468,"would_cite":false,"duration_ms":57758,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Damping is essential for stable rotating magnetic levitation, and the stable region in rotation speed and damping is bounded.","keywords":["magnetic levitation","rotating magnet","eddy-current damping","stability analysis","six-degree-of-freedom model","phase diagram","precessional motion","copper plate"],"falsifier":"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.","tokens_in":8663,"feed_emoji":"🧲","tokens_out":10378,"duration_ms":90535,"temperature":0.7,"pith_summary":"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.","feed_headline":"Stable levitation needs damping—but only within a bounded range","feed_subtitle":"A six-degree-of-freedom model and experiments map the window of rotation speed and damping that keeps the magnet aloft.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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$."],"forward_implications":["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."],"supporting_citations":[{"why":"First reported the orthogonal-configuration levitation phenomenon that this paper analyzes.","marker":"[13]"},{"why":"Supplies the systematic experimental baseline whose stability behavior the model must reproduce.","marker":"[17]"},{"why":"Provides the reduced-degree-of-freedom scaling laws that this paper extends by including horizontal equilibrium and damping.","marker":"[18]"},{"why":"Supplies the eigenvalue criterion used to decide whether the linearized system is stable.","marker":"[19]"}],"fun_headline_variants":["Damping sweet spot: too little or too much kills levitation","Stable levitation requires damping within a bounded window","Rotation speed and damping delimit magnetic levitation stability","No damping, no levitation; too much, same result","Bounded damping range discovered for stable levitation"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Damping sweet spot: too little or too much kills levitation","Stable levitation requires damping within a bounded window","Rotation speed and damping delimit magnetic levitation stability","No damping, no levitation; too much, same result","Bounded damping range discovered for stable levitation"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00063,"raw_usage":{"total_tokens":2875,"prompt_tokens":873,"completion_tokens":2002,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":489,"completion_tokens_details":{"reasoning_tokens":1922}},"tokens_in":489,"tokens_out":2002,"duration_ms":14954,"temperature":1.0,"reasoning_tokens":1922,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T18:40:06.572394+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"First reported the orthogonal-configuration levitation phenomenon that this paper analyzes."},{"cited_title":"Committee, Problems for the 37th iypt 2024, https:// www.iypt.org/problems/problems-iypt-2024/ (2024), accessed: 27 May 2025","cited_arxiv_id":null,"evidence_quote":"Supplies the systematic experimental baseline whose stability behavior the model must reproduce."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the reduced-degree-of-freedom scaling laws that this paper extends by including horizontal equilibrium and damping."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the eigenvalue criterion used to decide whether the linearized system is stable."}],"review_version":1}