{"id":"9af90f47-c789-4a6c-aa6c-d5b0fdb6718f","arxiv_id":"2505.22844","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Negative stiffness plus ground stiffness tunes phonon dispersion to hit zero frequency at chosen wavenumbers, opening passive wavenumber band gaps demonstrated in 1D and 2D magnetic lattices.","lead":"This paper shows that negative stiffness springs can bend wave-frequency curves so that waves disappear at chosen speeds, creating wavenumber band gaps without time modulation. The authors demonstrate this with floating magnetic disks in one and two dimensions, claiming the first experimental two-dimensional wavenumber gap.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Bi-stable stabilization changes the unit-cell potential, so the measured wavenumber gaps may be an artifact of the stabilization rather than evidence for the predicted negative-stiffness dispersion.","rationale":"The theory of the dispersion anomalies is elementary and internally consistent: a linear chain with K_g > 0 and K_in < 0 satisfying the ratio condition has omega^2(k) < 0 on an interval of k, giving a wavenumber band gap. The numerical simulations are consistent with this linear model. The single most load-bearing assumption in the experimental claim is that the physical magnetic lattices with local bi-stability actually implement the same linear Hamiltonian when probed dynamically. The paper's own Methods introduce a modification precisely to avoid the instability that the negative-stiffness model predicts, but no measurement or simulation establishes that the modified unit cell has the claimed negative K_in at the operating point. The reader's weakest assumption already identifies this. My proposed static force-displacement measurement would settle whether the experiment's gap is intrinsic. Since the reader already made the verdict conditional on this kind of evidence, my stress-test does not shift the verdict.","tokens_in":10933,"tokens_out":9092,"duration_ms":103732,"concrete_test":"Measure the static force-displacement relation of a single disk in the actual bi-stable unit cells (R=-3.4, R=-0.13, 2D R=-0.5), linearize at the static equilibrium, and recompute the dispersion of the chain/lattice using the measured stiffness values. If the recomputed band structure has no imaginary-frequency interval at the reported gap wavenumbers (or the gap shifts), then the experimental gap is an artifact of the stabilization, not the negative-stiffness model.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The experimental wavenumber-gap cases (1D R=-3.4 and R=-0.13, 2D R=-0.5) are all realized by modifying the unit cell with local bi-stability to suppress global instability. The Methods say: 'To avoid the global instability due to the wavenumber band gaps ... we introduce local instability within each unit-cell,' via 'an almost constant potential energy at the middle of the unit-cell' or 'two potential wells with very small energy barrier in-between.' The analytical dispersion in Fig.1(b) is derived from a linear spring-mass Hamiltonian with constant stiffnesses K_in and K_g. A flat-bottomed or double-well potential has zero or positive linearized stiffness at the stable equilibrium where the disk sits, not the negative K_in assumed in the model. The experimental 2D-FFT therefore samples a different Hamiltonian from the one whose dispersion is plotted analytically. The paper does not report force-displacement curves of the modified unit cells, nor does it simulate the full nonlinear potential to show that the measured branches coincide with the negative-stiffness branches. Since the title and abstract claim an experimental observation of the predicted anomalies, the missing equivalence between the stabilized system and the linear model is the load-bearing weakness.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a passive mechanism for engineering phonon dispersion anomalies by combining negative physical stiffness with ground stiffness in spring-mass lattices. The authors derive analytical dispersion relations for monoatomic, diatomic, locally resonant, nonlocal, and two-dimensional square lattices, and identify stiffness-ratio conditions under which the dispersion branch touches zero frequency at nonzero wavenumbers or develops a wavenumber band gap (an interval of imaginary frequencies). They support the analytical results with numerical 2D-FFT simulations and with experiments on magnetic disk lattices in one and two dimensions, claiming the first experimental observation of wavenumber band gaps in higher dimensions.","tokens_in":11139,"tokens_out":4735,"duration_ms":55529,"significance":"If fully substantiated, the work is significant for dispersion engineering: it offers a purely passive route to wavenumber band gaps, without time modulation, electron-phonon coupling, or long-range interactions, and it provides explicit analytical design conditions. The stiffness ratios in the experiments are computed from geometry and independently characterized magnetic constants rather than fitted to the measured dispersion, which is a strength. The analytical dispersion relations are direct mass-spring results, and the numerical stable-mode-superposition procedure is transparent. However, the experimental validation of the headline wavenumber-gap claims relies on locally bi-stabilized unit cells whose linearized dynamics are not shown to match the constant-negative-stiffness model, and the two-dimensional evidence covers only one row of the Brillouin zone. These gaps currently prevent the experimental claims from being fully supported.","major_comments":[{"comment":"The experimental wavenumber-gap cases (1D R=-3.4 and R=-0.13, and 2D R=-0.5) are realized by modifying the unit cell with 'an almost constant potential energy at the middle of the unit-cell' or 'two potential wells with very small energy barrier in-between,' as stated in the Methods. These modifications change the unit-cell potential from the linear negative-stiffness potential used in the analytical dispersion (Fig.1(b)) and in the numerical simulations. The paper does not report force-displacement curves of the modified unit cells, nor does it simulate the full nonlinear potential to show that the measured branches coincide with the constant-K_in branches of the model. Without this equivalence, the experimental dispersion curves in Figs.2(f,g) and 6(c) cannot be attributed to the predicted negative-stiffness dispersion, and the central experimental claim of the paper remains unsupported.","section":"Methods, Unit-cell design"},{"comment":"The experimental validation of the two-dimensional wavenumber gap is based on the Γ-X row only: the Methods state that in 2D the authors 'excite a disk on one of the structure edges ... and analyze the motion of all disks in the considered row.' The theoretical analysis in 'Dispersion anomalies in 2D lattices' (Fig.5) shows that the conditions for anomalies differ among the Γ-X, X-M, and M-Γ segments and between branches. A measured gap along one row of a 9×11 lattice is not sufficient to establish a 'complete wavenumber band gap' or the 'first experimental observation of wavenumber band gaps in higher dimensions.' Data along at least the other high-symmetry directions, or a full 2D FFT over the entire lattice, are required to support that claim.","section":"Experimental validation in 2D; Fig.6"},{"comment":"Equation (2) defines initial conditions as superpositions of stable modes only, deliberately excluding the κ-gap (imaginary-frequency) modes. The numerical 2D-FFT overlays therefore validate the stable branches but cannot validate the existence or location of the wavenumber gap itself; the gap evidence rests entirely on the experimental stabilization discussed in the first major comment. The manuscript should state this limitation explicitly and explain what, if anything, in the numerical or experimental data directly confirms the imaginary-frequency interval rather than only the stable branches that border it.","section":"Methods, Numerical simulations"}],"minor_comments":[{"comment":"There are several typographical errors, including 'wavemnubers' for 'wavenumbers' and 'Wills lattice' for 'Willis lattice'; these should be corrected throughout.","section":"Introduction"},{"comment":"The phrase 'passive- or active- experimental observation' has unusual hyphens and should be rephrased, for example as 'first experimental observation, passive or active.'","section":"Abstract and Discussion"},{"comment":"Equation (1) uses f(d), e_i, and d without defining all symbols and the sign convention; define the inter-disk distance, the unit vectors, and the force function before using the formula.","section":"Methods, Eq. (1)"},{"comment":"The text refers to scenarios 1, 2, and 3, but the panels in Fig.1(b) are labeled I-VI without an explicit mapping; add labels or a sentence connecting each panel to its scenario.","section":"Fig.1"},{"comment":"The phrase 'complete wavenumber band gaps' is used ambiguously: clarify whether it means a gap that exists for all frequencies at a given wavevector interval, or a gap that spans all wavevectors in the Brillouin zone; the current wording invites the latter reading, which is not what is measured.","section":"Discussion"}],"recommendation":"major_revision","confidential_remarks":"The analytical and numerical portions are solid and the parameter-free stiffness-ratio characterization is a real strength. The main risk is that the experimental validation of the headline k-gap claims rests on bi-stabilized unit cells whose linearized dynamics are not characterized, and on a single row of the 2D Brillouin zone. If the authors can supply force-displacement measurements or full nonlinear simulations establishing the equivalence, and measure additional Brillouin-zone directions, the paper would be much stronger. There is also a scope question: the contribution is a mechanics/physics metamaterials result rather than an applied-physics device demonstration, but it is plausibly within the journal's broad scope."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The theory is straightforward and correct. For a monoatomic chain with inter-stiffness K_in and ground stiffness K_g, omega^2 = (K_g + 4K_in sin^2(ka/2))/m; setting K_in < 0 and -4 < K_g/K_in < 0 gives a band of imaginary frequencies on a tunable wavenumber interval. That is a genuine wavenumber gap with a simple passive mechanism, and the authors generalize it to diatomic, locally resonant, non-local, and 2D lattices. The numerical 2D-FFT overlays match the analytical curves, and the stiffness ratios used in experiment are derived from geometry and magnetic characterizations, not fit to the data. All of that is solid.\n\nThe experimental apparatus is impressive too: magnetic disks on an air table with tunable boundary and neighbor magnet arrangements. The stable-case measurements (R=-5.4, R=+8.5) look plausible.\n\nThe soft spot is the bi-stability workaround. In the Methods, the authors say that for the wavenumber-gap designs (R=-3.4, R=-0.13, and the 2D R=-0.5 case), they introduce \"local instability within each unit-cell\" via an almost constant potential energy or two potential wells with a small barrier. The linear theory assumes constant negative inter-stiffness. A disk sitting in a double-well or flat-bottomed potential has non-negative linearized stiffness at its stable equilibrium. So the measured dispersion of the stabilized structure may correspond to a different Hamiltonian than the one plotted analytically. The paper gives no force-displacement curves for the modified unit cells and no nonlinear simulation showing that the measured branches coincide with the negative-stiffness branches. That gap is load-bearing for the experimental claim.\n\nA second weakness: the 2D experiment only measures the Gamma-X row of a 9x11 lattice. That supports a directional gap, not a complete wavenumber gap in the first Brillouin zone.\n\nThe theoretical contribution alone is worth a referee's time, but the experimental claim needs another round of characterization. I would ask for force-displacement data and either a nonlinear model or a convincing argument that the bi-stability does not alter the effective linear stiffness at the operating point.\n\nVerdict: worth engaging, but the current version is not credible as evidence for the headline claim. Send to peer review and let a good referee push on the stabilization equivalence.","headline":"The theoretical framework is clean, but the experimental wavenumber-gap evidence likely rests on a stabilization scheme that changes the equilibrium stiffness, so the first 2D observation claim is not yet supported.","tokens_in":11682,"tokens_out":4470,"would_cite":false,"duration_ms":46673,"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":"This paper shows that pairing negative physical stiffness with ground stiffness can make phonon frequency reach zero at any chosen nonzero wavenumber, and reports the first 2D wavenumber band gaps in passive lattices.","keywords":["dispersion engineering","wavenumber band gap","negative stiffness","phonon dispersion","magnetic metamaterials","zero-frequency phonon","phononic crystals","locally resonant metamaterials"],"falsifier":"A reader could settle this by measuring the single-unit-cell force-displacement curve to confirm the inter-stiffness is negative at the equilibrium the model assumes, or by extracting the two-dimensional dispersion along the $X$–$M$ and $M$–$\\Gamma$ segments of the $R=-0.5$ lattice and checking whether the imaginary-frequency interval appears there.","tokens_in":10679,"feed_emoji":"🧲","tokens_out":13288,"duration_ms":125014,"temperature":0.7,"pith_summary":"By pairing negative physical stiffness with ground stiffness in spring-mass lattices, the paper shows that phonon dispersion can be sculpted so the frequency reaches zero at any chosen nonzero wavenumber, not just at the edges of the first momentum zone. When the ratio $R=K_g/K_{in}$ of ground to inter stiffness lies between $-4$ and $0$, the squared frequency becomes negative on a tunable interval of wavenumbers, meaning those modes are wavenumber-gapped: they grow instead of propagating. The authors realize these conditions in magnetic disk lattices and report the first experimental observation of such wavenumber band gaps in a two-dimensional lattice, a regime previously accessible only through time modulation or other active mechanisms. They also extend the design rules to diatomic phononic crystals, locally resonant metamaterials, and nonlocal chains. If correct, the work turns a historically exotic dispersion anomaly into a design variable for passive materials.","feed_headline":"First 2D wavenumber band gap opens in passive lattice","feed_subtitle":"Negative stiffness plus ground stiffness lets phonon frequency hit zero at chosen wavenumbers","key_machinery":"The load-bearing object is the dispersion relation of the grounded monoatomic chain, $\\omega^2=(K_g+4K_{in}\\sin^2(ka/2))/m$, together with the sign structure of the stiffness ratio $R=K_g/K_{in}$. The wavenumber band gap appears exactly when $\\omega^2<0$, which happens for a tunable interval of $k$ when $-4<R<0$; the landing wavenumber $\\kappa_-$ is the point where $\\sin^2(ka/2)=-K_g/(4K_{in})$. Negative physical stiffness — a coupling that pushes in the direction of its own displacement — is realized experimentally through magnetic repulsion between floating disks and fixed boundary magnets, and a local bi-stability inside each unit cell is introduced to suppress the global instability that the imaginary-frequency modes would otherwise cause.","core_discovery":"The central claim is that a passive lattice with one negative physical stiffness — either the inter-site coupling or the ground coupling — can support zero-frequency phonon anomalies at arbitrary nonzero wavenumbers and can open wavenumber band gaps, including fully in two dimensions. In the grounded monoatomic chain, the dispersion is $\\omega^2 = (K_g + 4K_{in}\\sin^2(ka/2))/m$; with $K_g/K_{in}$ between $-4$ and $0$, the right-hand side is negative for a band of wavenumbers, giving imaginary frequencies. The critical wavenumber where the curve lands on the wavenumber axis is determined by the stiffness ratio, so both the location and width of the gap are tunable. The paper claims to verify this analytically, numerically, and experimentally in one-dimensional magnetic lattices at ratios $R=-5.4,-3.4,-0.13,+8.5$, and in two-dimensional square magnetic lattices at $R=0,+4,-0.5$, the last being the first experimental wavenumber band gap in higher dimensions, passive or active.","pith_inferences":["Our inference: the same ratio condition should transfer to any physical mechanism that provides negative stiffness — electrostatic, geometric buckling, or pre-stressed elements — not only magnetic disks, so the design rule may define a general passive route to wavenumber gaps.","Our inference: because the wavenumber-gap modes grow exponentially in time, the local bi-stability used in the experiments likely changes the effective linearization; a direct measurement of the single-unit-cell force-displacement curve at the operating point would show whether the observed gap is the linear gap or a nonlinear artifact.","Our inference: the two-dimensional demonstration only analyzes the $\\Gamma$–$X$ row of the lattice; checking the $X$–$M$ and $M$–$\\Gamma$ segments would test whether the complete two-dimensional wavenumber gap is actually closed as the theory states.","Our inference: wavenumber-gapped modes do not propagate but grow in time, so a practical use could be spatial region selection or energy localization in finite structures, where the unstable growth is confined by the local bi-stability."],"forward_implications":["Zero-frequency phonon anomalies can be placed at arbitrary chosen wavenumbers in a passive lattice, without electron-phonon coupling or long-range interactions.","Wavenumber band gaps — intervals of wavenumber with imaginary frequency — can be opened in passive media, a regime previously requiring time modulation or other active mechanisms.","The design rules carry over to diatomic phononic crystals, locally resonant metamaterials, and nonlocal chains, each with explicit conditions on the stiffness ratio.","In two dimensions, a complete wavenumber band gap is claimed to be achievable and experimentally observed for the first time, based on the $\\Gamma$–$X$ row of a square lattice.","Transmission through a finite lattice can be engineered to begin at zero frequency even when ground stiffness is present, because the gap is in wavenumber rather than frequency."],"supporting_citations":[{"why":"Supplies the classic zero-frequency phonon softening at a finite wavenumber that the paper reproduces by design rather than through electron-phonon coupling.","marker":"[10]"},{"why":"The long-standing speculation that such anomalies require long-range interactions, which the paper's negative-stiffness mechanism avoids.","marker":"[16]"},{"why":"One active mechanism (time modulation) previously used to open wavenumber band gaps, the baseline the passive design must match.","marker":"[23]"},{"why":"The other active route to wavenumber band gaps, which the paper contrasts with its passive approach.","marker":"[25]"},{"why":"Shows nonlocal coupling that softens dispersion but does not reach zero frequency, the gap filled by the nonlocal generalization.","marker":"[26]"},{"why":"Establishes ground stiffness as a dispersion-engineering knob, the second ingredient of the design rule.","marker":"[40]"},{"why":"Prior all-flat-band phononic metamaterial framework that the zero-frequency anomaly design extends.","marker":"[42]"},{"why":"Supplies the duality mapping used to relate stiffness ratios $R$ and $-(R+4)$ and the self-dual point $R=-2$.","marker":"[50]"},{"why":"Provides the magnetic repulsion model used to compute inter and ground stiffnesses from magnet geometry.","marker":"[53]"}],"fun_headline_variants":["Magnets create tunable wavenumber band gaps in passive lattices","Negative stiffness from magnets shapes phonon dispersion","Zero-frequency phonons tuned to any wavenumber via magnetic coupling","First passive 2D wavenumber band gap from magnetic lattices","Magnetic negative stiffness opens band gaps at any wavenumber"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the bi-stable magnetic lattices, at the operating point used in the experiments, respond to small motions exactly like the negative-stiffness spring-mass model whose dispersion is plotted, so the observed wavenumber gaps come from the design rather than from the stabilization trick.","fun_headline_variants_meta":{"raw":{"variants":["Magnets create tunable wavenumber band gaps in passive lattices","Negative stiffness from magnets shapes phonon dispersion","Zero-frequency phonons tuned to any wavenumber via magnetic coupling","First passive 2D wavenumber band gap from magnetic lattices","Magnetic negative stiffness opens band gaps at any wavenumber"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000811,"raw_usage":{"total_tokens":3568,"prompt_tokens":967,"completion_tokens":2601,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":583,"completion_tokens_details":{"reasoning_tokens":2514}},"tokens_in":583,"tokens_out":2601,"duration_ms":21781,"temperature":1.0,"reasoning_tokens":2514,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T13:00:01.967822+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A reader could settle this by measuring the single-unit-cell force-displacement curve to confirm the inter-stiffness is negative at the equilibrium the model assumes, or by extracting the two-dimensional dispersion along the $X$–$M$ and $M$–$\\Gamma$ segments of the $R=-0.5$ lattice and checking whether the imaginary-frequency interval appears there.","supporting_citations":[{"cited_title":"Kohn, Image of the fermi surface in the vibration spectrum of a metal, Physical Review Letters2, 393 (1959)","cited_arxiv_id":null,"evidence_quote":"Supplies the classic zero-frequency phonon softening at a finite wavenumber that the paper reproduces by design rather than through electron-phonon coupling."},{"cited_title":"Brillouin, Wave propagation in periodic structures, McGraw- Hill Book Company google schola2, 2022 (1946)","cited_arxiv_id":null,"evidence_quote":"The long-standing speculation that such anomalies require long-range interactions, which the paper's negative-stiffness mechanism avoids."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"One active mechanism (time modulation) previously used to open wavenumber band gaps, the baseline the passive design must match."},{"cited_title":"Onset of wavenumber bandgaps via alternating Willis coupling signs","cited_arxiv_id":"2412.06798","evidence_quote":"The other active route to wavenumber band gaps, which the paper contrasts with its passive approach."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Shows nonlocal coupling that softens dispersion but does not reach zero frequency, the gap filled by the nonlocal generalization."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes ground stiffness as a dispersion-engineering knob, the second ingredient of the design rule."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Prior all-flat-band phononic metamaterial framework that the zero-frequency anomaly design extends."},{"cited_title":"Vitelli, M","cited_arxiv_id":null,"evidence_quote":"Supplies the duality mapping used to relate stiffness ratios $R$ and $-(R+4)$ and the self-dual point $R=-2$."},{"cited_title":"Jiao and S","cited_arxiv_id":null,"evidence_quote":"Provides the magnetic repulsion model used to compute inter and ground stiffnesses from magnet geometry."}],"review_version":1}