{"id":"42944c10-07ec-4664-93f5-ea611cbe5f45","arxiv_id":"2412.06798","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"Alternating the sign of Willis coupling in a two-layer elastic unit cell creates wavenumber bandgaps with exactly known width and limits, while the overall band structure becomes reciprocal.","lead":"This paper shows that a rod with alternating sign of its Willis coupling, a bias that makes waves travel differently forward and backward, develops gaps in the allowed wavenumbers. This gives an analytic recipe for filtering waves by spatial scale rather than by frequency.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claimed wavenumber bandgap is a temporal instability: Eq. (9) forces cos Ω < −1 for q inside the gap, so Ω has positive imaginary part and the mode grows in time rather than being a passive stop band.","rationale":"The reader's verdict is CONDITIONAL and I agree that the paper should not be accepted unconditionally, but my central concern differs from the reader's stated weakest assumption. The reader focused on physical realizability of alternating Willis coupling; I focus on an internal physical interpretation issue. The dispersion formula and gap-width formula are correct, and the finite-element agreement is plausible though not reproducible from the text. The load-bearing issue is that the 'bandgap' is actually a region of linear temporal instability: for q in the gap, Eq. (9) yields complex Ω with both signs, and the positive-imaginary solution grows exponentially. This is not a passive stop band, and it undermines the claim that alternating Willis coupling 'opens wavenumber bandgaps' in the usual physical sense. Since the mathematical result itself survives and the paper could be revised to analyze stability and reframe the phenomenon, a conditional verdict remains appropriate. My concrete test would settle the interpretation by direct evaluation of the complex frequency or time-domain simulation; if it confirms growth, the paper must explicitly discuss the instability before the bandgap language can be accepted.","tokens_in":17,"tokens_out":24434,"duration_ms":743711,"concrete_test":"Evaluate the complex eigenfrequency at the gap center q = π using Eq. (9): compute cos Ω = −(1+ν²)/(1−ν²) for a representative ν = 0.5, giving Ω = π ± 1.0986 i. Then run a time-domain simulation of the periodic PDE with a Bloch initial condition u(x,0) = e^{iπx}φ(x) and verify whether the amplitude grows at the predicted rate. If the growing solution is present, the 'bandgap' is a temporal instability and the paper must reframe or provide a stability analysis.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The algebraic derivation of Eq. (9) is internally consistent: with equal-magnitude opposite-sign Willis couplings, the unit-cell transfer matrix has trace 2[ν² + (1−ν²)cos Ω], giving the stated reciprocal dispersion relation. The load-bearing problem is the physical meaning of the resulting 'wavenumber bandgap'. For any real q inside the gap, Eq. (9) implies (cos q − ν²)/(1 − ν²) < −1, so cos Ω < −1 and Ω must be complex. At the gap center q = π, cos Ω = −(1+ν²)/(1−ν²), whence Ω = π ± iα with α = arcosh((1+ν²)/(1−ν²)) > 0. Because the underlying PDE is real, the two signs form a conjugate pair, and the + sign is a linearly growing mode e^{ατ}. Thus every real Bloch wavenumber in the purported bandgap is temporally unstable for any ν ∈ (0,1), even though the uniform rod is stable for ν < 1. A true bandgap would instead suppress real-frequency propagation via spatial evanescence; here real-frequency excitation always yields real q because the argument of arccos in Eq. (9) lies in [−1,1]. The paper's statement that 'complex frequencies emerge in wavenumber bandgap' stops short of noting that one of these frequencies has positive imaginary part. Therefore the central claim 'wavenumber bandgaps open' is better described as an instability band, and the analogy to passive phononic-crystal frequency bandgaps in Table 1 is formal rather than physical.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies a one-dimensional Willis-type (axially moving rod) wave equation and shows that a unit cell made of two identical segments with equal-magnitude, opposite-sign Willis coupling yields the reciprocal dispersion relation q = ± arccos(ν² + (1−ν²) cos Ω), Eq. (9). The author derives explicit lower and upper limits of the resulting wavenumber gaps, Eqs. (10a)-(10b), and their width, Eq. (11), draws a formal analogy to frequency bandgaps of a special bi-layered phononic crystal, and reports finite-element verification. The algebraic derivation is transparent and parameter-free, with the expected limit ν=0 recovering the linear reciprocal dispersion.","tokens_in":6611,"tokens_out":13188,"duration_ms":131814,"significance":"If the result is correctly interpreted, it provides a clean, exactly solvable example of complex-frequency bands in a periodic Willis-type medium, which could be of interest for active metamaterial designs. The derivation is self-contained and does not rely on fitted parameters or opaque numerics. However, the central physical interpretation is currently misleading: the predicted 'wavenumber bandgaps' are in fact temporal instability bands, not passive stop bands. The paper explicitly notes that complex frequencies appear, but it does not state that one of the conjugate pair has a positive imaginary part, nor does it reconcile this with the claim of a bandgap. A revision that reframes the phenomenon as an instability/complex-frequency band and adds a stability analysis would make the contribution sound and publishable.","major_comments":[{"comment":"The claimed 'wavenumber bandgap' is a temporal instability, not a passive stop band. For any real q inside the gap, Eq. (9) implies cos Ω = (cos q − ν²)/(1−ν²) < −1, so Ω = π ± iα with α > 0. Since the underlying PDE is real, the two signs form a conjugate pair, and the solution with positive imaginary part grows exponentially in time as e^{ατ}. By contrast, for any real driving frequency Ω, the right-hand side of Eq. (9) lies in [2ν²−1, 1], so q is always real and no frequency is suppressed. The paper's statement that 'complex frequencies emerge in wavenumber bandgap' stops short of noting the positive imaginary part. The footnote on page 1 restricting ν to [0,1) 'to avoid possible dynamical instabilities' is therefore in direct tension with the phenomenon presented. The central claim should be revised to state that the configuration produces instability bands, with a stability analysis.","section":"Elastic rod with periodic Willis coupling, Eq. (9) and Fig. 3"},{"comment":"The label 'driven-wave dispersion relation' for Eq. (9) is misleading. When Ω is real (the usual driven situation), Eq. (9) gives real q for all Ω; the complex frequencies appear only if one instead fixes a real q and solves for Ω, which is an eigenvalue/initial-value problem. The manuscript does not distinguish these two settings. This distinction is load-bearing because it determines whether the gap acts as a filter for real-frequency excitation (it does not) or as an instability band for real-wavenumber perturbations (it does). Please clarify the problem statement and the meaning of 'driven' in this context.","section":"Elastic rod with periodic Willis coupling, Eqs. (6) and (9)"},{"comment":"The analogy in Table 1 is formal rather than physical. In a bi-layered phononic crystal, a frequency bandgap is a range of real frequencies with no real wavenumber, and the corresponding Floquet exponent is complex, giving spatial evanescence. Here, the 'wavenumber bandgap' is a range of real wavenumbers with no real frequency, and the corresponding temporal frequency is complex, giving temporal growth. These are dual in a purely algebraic sense, but their physical consequences are opposite: spatial decay versus temporal instability. The paper should explicitly discuss this difference before claiming equivalence of the roles of ν and β.","section":"Similarities between reversed-sign Willis couplings and bi-layered periodicity, Table 1"}],"minor_comments":[{"comment":"The finite-element verification is stated but not documented: no mesh parameters, element type, boundary conditions, or error metrics are given, and the procedure used to obtain complex-frequency branches is not described. Since the analytic derivation is central and convincing, this is a reproducibility issue rather than a scientific flaw, but the details should be provided or a reference given.","section":"Fig. 3 and 'finite element method' statement"},{"comment":"In Eq. (8), the notation ν+ and ν− is ambiguous: it is not immediately clear that these are magnitudes of the signed couplings rather than signed values. Please define them explicitly, e.g., as positive magnitudes, so that 'setting ν+ = ν− = ν' cannot be misread as both segments having the same sign.","section":"Eq. (8) and notation"},{"comment":"The 'Real Imag.' inset in Fig. 3 is difficult to read. It would help to label the branches with positive and negative imaginary parts explicitly and to state in the caption that the positive-imaginary branch corresponds to temporal growth.","section":"Fig. 3"},{"comment":"The nondimensionalization that leads from Eq. (1) to Eq. (2) is correct but terse; showing the cancellation of the ω0 and c factors would make the derivation easier to follow for readers not familiar with the cited references.","section":"Introduction, Eq. (2)"}],"recommendation":"major_revision","confidential_remarks":"The algebraic result is sound and likely correct, and the paper is short and readable. The main risk is that the term 'wavenumber bandgap' will be understood as a passive stop band, which it is not; the revision must reframe the result as a complex-frequency/instability band or demonstrate some practical regime where the instability is avoided. The lack of FEM reproducibility is secondary but should also be addressed. I would not reject the paper, since the central derivation is clean and the phenomenon is well-defined once correctly interpreted."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The one thing you should know: Eq. (9) is derived correctly, and the formulas (10)-(11) for the gap limits and width are new and exact. But the paper's central physical claim—that alternating Willis signs open wavenumber bandgaps—does not survive contact with the dispersion relation. For any real wavenumber inside the alleged gap, Eq. (9) forces cos Ω < -1, so Ω is complex. Because the underlying PDE is real, the two roots are a conjugate pair, and one of them has positive imaginary part. The mode grows in time. This is an instability band, not a passive phononic bandgap. The paper itself notes that 'complex frequencies emerge' but stops short of recognizing that one of those frequencies represents growth. The analogy to bilayered phononic crystals in Table 1 is formal: those structures have real frequency with complex wavenumber (spatial evanescence); here it's real wavenumber with complex frequency (temporal growth).\n\nWhat is genuinely new is the observation that two identical layers with opposite Willis coupling cancel the phase shift and give a reciprocal dispersion relation, and the exact closed-form description of where the q-gaps appear. That is a clean, checkable result, and it may be useful if reframed as a non-Hermitian amplification mechanism. But as written, the interpretation is wrong.\n\nOther soft spots: the finite-element verification is stated but no code, mesh, or error metrics are provided, so the numerical claim is unverifiable. The physical realization is deferred to feedback controllers, which is honest but leaves the model speculative. Also, the paper restricts ν ∈ [0,1) to avoid 'dynamical instabilities' in the uniform rod, yet the alternating structure itself is unstable for any ν ≠ 0—this point is never addressed.\n\nWho is this for? Researchers working on non-reciprocal wave propagation and non-Hermitian phononics might find the algebra useful, but the paper needs major revision before it can be taken as a valid bandgap result. I would send it to peer review—a good referee could get the author to reframe it as an instability and include a stability analysis—but I would not let it pass as is.\n\nRecommendation: engage with it, but treat the central claim as incorrect until the author reanalyzes the sign of the imaginary frequency.","headline":"The algebra is right, but the headline claim is wrong: the 'wavenumber bandgap' is a temporal instability, not a passive stop band.","tokens_in":7125,"tokens_out":6184,"would_cite":false,"duration_ms":57679,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["74J05","74Q10"],"pacs":[],"model":"deepseek-v4-flash","headline":"Periodic sign-flips of a wave-bias term create wavenumber bandgaps in an elastic rod.","keywords":["wavenumber bandgaps","Willis coupling","non-reciprocal waves","periodic structures","transfer matrix","dispersion diagram","phononic crystals","momentum bias"],"falsifier":"Measure the driven-wave dispersion of a rod whose unit cell alternates between $+\\nu$ and $-\\nu$ Willis coupling segments (for example, via feedback controllers or moving segments) and check whether a stop band in wavenumber appears around $q=\\pm\\pi$ at frequencies $\\Omega=(2p-1)\\pi$ for $|\\nu|>0$; absence of that gap would contradict Equation (9). A simpler numerical check is to integrate Equation (1) with piecewise-constant $\\nu(x)$ and compare the Bloch wavenumber to the closed-form relation.","tokens_in":6088,"feed_emoji":"🌊","tokens_out":12893,"duration_ms":104009,"temperature":0.7,"pith_summary":"An elastic rod whose unit cell is split into two identical halves with opposite signs of Willis coupling—a momentum-bias term in the wave equation—has a reciprocal dispersion relation even though each half is non-reciprocal on its own. The paper derives the exact dispersion relation $q = \\pm \\cos^{-1}\\left(\\nu^2 + (1-\\nu^2)\\cos\\Omega\\right)$ and shows that wavenumber bandgaps open for any nonzero normalized modulation speed $\\nu$, with gap width $\\Delta q = 2(\\pi - \\cos^{-1}(2\\nu^2-1))$. This gives a static, parameter-tunable route to wavenumber bandgaps that requires neither damping, temporal modulation, nor spatiotemporal stiffness variation. The same formulas are shown to mirror the frequency-bandgap expressions of bi-layered phononic crystals with zero frequency contrast and nonzero impedance contrast.","feed_headline":"Alternating wave-bias signs open wavenumber bandgaps","feed_subtitle":"Two identical non-reciprocal layers with opposite bias form a reciprocal rod with tunable, analytic gaps.","key_machinery":"The load-bearing tool is the transfer matrix of the unit cell, obtained by multiplying the transfer matrices of the two half-cells, each of normalized length $\\xi=1/2$. Because the Willis couplings of the two halves are equal in magnitude and opposite in sign, the exponential phase factor $e^{-\\frac{i}{2}(\\nu_+-\\nu_-)\\Omega}$ reduces to $1$, so the product becomes $T = Y(\\nu)Y(-\\nu)$, a $2\\times 2$ matrix of unit determinant. The dispersion relation follows from writing the eigenvalues as $\\lambda = e^{iq}$ and solving $\\det(T-\\lambda I)=0$, which yields the cosine relation for $q$. The bandgap limits are then read directly from the argument of the inverse cosine, and the width formula follows by subtraction. The paper uses the same transfer-matrix setting to derive the skew angle $\\phi=\\tan^{-1}\\nu$ of the single-layer non-reciprocal dispersion, showing that the alternation removes this skew entirely.","core_discovery":"The paper's central claim is that a unit cell made of two identical elastic layers with equal-magnitude, opposite-sign Willis couplings has the reciprocal driven-wave dispersion relation $q = \\pm \\cos^{-1}\\left(\\nu^2 + (1-\\nu^2)\\cos\\Omega\\right)$, so that wavenumber bandgaps open at odd multiples of $q=\\pm\\pi$ for every $\\nu\\neq 0$. The lower and upper limits of the $n$-th gap are $q_- = \\pm[2(n-1)\\pi + \\cos^{-1}(2\\nu^2-1)]$ and $q_+ = \\pm[2n\\pi - \\cos^{-1}(2\\nu^2-1)]$, giving the identical width $\\Delta q = 2(\\pi-\\cos^{-1}(2\\nu^2-1))$, which grows with $|\\nu|$ and reaches $2\\pi$ as $\\nu\\to 1$. The overall band structure is reciprocal even though each constitutive layer is non-reciprocal, and it is insensitive to the sign of $\\nu$. The paper also finds that complex frequencies appear inside the wavenumber gaps and that the formulas are numerically confirmed by finite-element simulations. It closes by drawing a formal analogy between this design and bi-layered phononic crystals with zero frequency contrast and nonzero impedance contrast, identifying $\\nu$ with the impedance contrast $\\beta$ and $q$ with $\\Omega$.","pith_inferences":["If the analogy to bi-layered phononic crystals is more than formal, the wavenumber gap should close continuously as $\\nu\\to 0$ with a square-root-like scaling, a behavior a transmission experiment on a finite stack could test.","The same sign-alternation idea might transfer to other momentum-biased wave systems (gyroscopic, rotating, or magneto-elastic), where a piecewise-constant bias of opposite signs could open wavenumber gaps without temporal modulation.","Since each half is non-reciprocal but the unit cell is reciprocal, inserting a defect that reverses the alternation pattern might create a localized mode with net non-reciprocal response; this is an extension the paper does not explore.","A lumped-parameter analog of Equations (9)-(11) could be built from a Willis monatomic lattice with alternating coupling signs, making the mechanism testable in a tabletop experiment."],"forward_implications":["Wavenumber bandgaps become achievable with a purely static spatial pattern, without damping, temporal modulation, or complex spatiotemporal stiffness functions.","Each half of the unit cell is non-reciprocal, yet the assembled cell is reciprocal; this offers a route to cancel non-reciprocity locally while keeping a gap.","The gap width is controlled by the single parameter $\\nu$ and saturates at the full Brillouin-zone width $2\\pi$ in the limit $\\nu\\to 1$.","The analytical limits and width match finite-element simulations, so the formulas can be used directly for design.","The formal mapping to bi-layered phononic crystals implies that the same mathematics governs frequency gaps from impedance contrast and wavenumber gaps from Willis-sign alternation."],"supporting_citations":[{"why":"Provides the axially-moving-rod wave equation and the transfer-matrix method used to build the unit-cell transfer matrix.","marker":"[1]"},{"why":"Establishes the concept of Willis materials and identifies the mixed-derivative term as Willis coupling arising from momentum bias.","marker":"[3]"},{"why":"Introduces the wavenumber phase shift $q_s$ for a Willis monatomic lattice, which Equation (6) borrows.","marker":"[5]"},{"why":"Provides the bi-layered phononic-crystal bandgap limits and width formulas that the paper compares with its own in Table 1.","marker":"[8]"},{"why":"Supports the statement that complex frequencies emerge inside wavenumber bandgaps.","marker":"[28]"}],"fun_headline_variants":["Flip Willis sign per layer to open wavenumber gaps","Opposite Willis signs make reciprocal gaps appear","Non-reciprocal layers pair into reciprocal bandgaps","Alternating Willis coupling unlocks analytic wavenumber gaps","Reversed Willis signs yield tunable wavenumber bandgaps"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The result holds only if a medium can actually be built in which the Willis coupling alternates in sign from one half of the unit cell to the next while density and stiffness stay constant; the paper defers such a realization to feedback-controlled designs, so the predicted gaps may not exist in any currently available passive material.","fun_headline_variants_meta":{"raw":{"variants":["Flip Willis sign per layer to open wavenumber gaps","Opposite Willis signs make reciprocal gaps appear","Non-reciprocal layers pair into reciprocal bandgaps","Alternating Willis coupling unlocks analytic wavenumber gaps","Reversed Willis signs yield tunable wavenumber bandgaps"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000192,"raw_usage":{"total_tokens":1348,"prompt_tokens":946,"completion_tokens":402,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":562,"completion_tokens_details":{"reasoning_tokens":323}},"tokens_in":562,"tokens_out":402,"duration_ms":4209,"temperature":1.0,"reasoning_tokens":323,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T14:24:26.419499+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the driven-wave dispersion of a rod whose unit cell alternates between $+\\nu$ and $-\\nu$ Willis coupling segments (for example, via feedback controllers or moving segments) and check whether a stop band in wavenumber appears around $q=\\pm\\pi$ at frequencies $\\Omega=(2p-1)\\pi$ for $|\\nu|>0$; absence of that gap would contradict Equation (9). A simpler numerical check is to integrate Equation (1) with piecewise-constant $\\nu(x)$ and compare the Bloch wavenumber to the closed-form relation.","supporting_citations":[],"review_version":1}