{"id":"9a95d2a1-9679-4426-a951-c8176857d644","arxiv_id":"1908.05064","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"For a spherical elastic inclusion, the paper gives explicit finite-frequency eigenvalues and eigenfunctions of the Neumann-Poincaré operator and constructs core-shell metamaterials whose anomalous resonances cloak sources inside a critical radius.","lead":"This paper derives the complete spectrum of the elastic Neumann-Poincaré operator on a sphere at finite frequency and uses it to construct elastic core-shell metamaterials that cloak by anomalous localized resonance, without the usual quasi-static approximation. It matters because it gives a rigorous finite-frequency route to invisibility cloaking in elasticity and a reusable spectral tool.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Theorem 5.1's restricted-source CALR proof relies on unverified determinant asymptotic estimates (5.17)-(5.20); without them the energy lower bound and boundedness estimates do not follow.","rationale":"The reader's weakest assumption was the restriction of the source to the expansion (5.7) and the unproved Remark 5.2 generalization. That is a real limitation, but the theorem as stated already restricts to that class. My stress-test identifies a more internal and, in my view, more load-bearing gap: the asymptotic estimates for the determinant d_{n,m} in (5.19)-(5.20) and the energy bound (5.22) are asserted without proof, and they are needed even for the restricted source class. The paper repeatedly relies on 'straightforward though tedious calculations' (Theorem 3.1, Proposition 3.6, Theorem 5.1) and ships no code or symbolic verification, so these omitted derivations cannot be checked by the reader. There is no evidence the claims are false; the structure is plausible and the quasi-static reductions match known results. But a conditional acceptance should require either a published derivation of the determinant asymptotics or a reproducible symbolic/numeric verification of (5.17)-(5.22). I therefore keep the reader's CONDITIONAL verdict, with the condition sharpened to explicitly cover the omitted determinant estimates. My agreement is partial because I emphasize a different concern than the reader's stated weakest assumption, though both point to gaps in Theorem 5.1's proof.","tokens_in":23790,"tokens_out":9563,"duration_ms":93839,"concrete_test":"Use a computer algebra system to compute d_{n,m} exactly from (5.9) for the representative parameters of Remark 5.3 (n0=50, omega=5, ri=0.8, re=1, mu=mu_breve=1, and p2 chosen so that (5.15) holds), and verify: (i) |d_{n0,m}| is within, say, 10% of rho^{2n0}; (ii) |d_{n,m}| >= rho^{2n0} + rho^{2n} for n=N,...,n0+50; and (iii) the energy inequality (5.22) with the same parameters. If any of these fail, the proof of Theorem 5.1 collapses; if they pass, the central estimates are at least numerically corroborated in the claimed regime.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The CALR conclusion of Theorem 5.1 rests on four asymptotic estimates stated without derivation: (5.17)-(5.20). In particular, (5.19)-(5.20) assert that the determinant d_{n,m} of the 4x4 transmission system (5.9) satisfies |d_{n0,m}| ≈ ρ^{2n0} and |d_{n,m}| ≥ ρ^{2n0} + ρ^{2n} for all n ≥ N, conditional only on the parameter choice (5.14)-(5.15). Everything downstream depends on these: the energy lower bound (5.22) and the boundedness estimate (5.25) both invoke the size of d_{n,m}. But d_{n,m} is a multiline expression in products of Bessel and Hankel functions and their derivatives, and the paper gives no derivation or error control for the claimed two-sided bounds; the proof of Theorem 5.1 merely says 'direct calculations show'. The same pattern appears earlier: Theorem 3.1 and Proposition 3.6 defer key identities to 'tedious calculations'. If the determinant asymptotics are off by even a polynomial factor in n or an extra ρ^n factor, the energy lower bound (5.22) and the critical-radius threshold r* would not follow. Thus even the restricted-source CALR claim, the paper's headline result, is not yet established. The source restriction (5.7) is an additional acknowledged limitation, but the omitted determinant estimates are more load-bearing because they are needed even for that restricted class.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies polariton resonances and cloaking due to anomalous localized resonance (CALR) for the time-harmonic elastic system at finite frequency, dropping the quasi-static assumption ω diam(Ω) ≪ 1. For spherical geometries it derives (Theorem 3.2) a claimed complete spectral system for the elastic Neumann–Poincaré operator (K^ω_SR)^* acting on vector spherical harmonics, with eigenvalues λ_{1,n}, λ_{2,n}, λ_{3,n} expressed through Bessel/Hankel coefficients. It then uses these spectra to construct material configurations: Theorem 4.1 claims polariton resonance for a homogeneous inclusion with no core when the shell Lamé parameter μ̂ is chosen with suitably large imaginary part, and Theorem 5.1 claims CALR for a core-shell-matrix structure with ˘μ = μ and μ̂ = −μ + iρ^{n0} + p_{2,n0} when the source's Newtonian potential has the restricted form (5.7) and is supported inside the critical radius r* = sqrt(r_e^3/r_i), with no resonance for sources outside r*.","tokens_in":24056,"tokens_out":7117,"duration_ms":66675,"significance":"The novelty is in moving the elastic ALR/CALR construction beyond the quasi-static limit using spectral information of the finite-frequency N-P operator; the claimed construction only requires violation of one of the two strong-convexity conditions, and taking the static limit is said to recover earlier results. The spectral system of Theorem 3.2, if correct, is a substantial independent contribution. The paper includes numerical illustrations (Figures 1–3) supporting the parameter conditions. However, the paper does not supply machine-checked derivations, and the two central theorems depend on asymptotic estimates that are quoted as 'direct calculations' rather than proved, which currently prevents verification of the headline claims.","major_comments":[{"comment":"The proof of Theorem 5.1 is built on four asymptotic estimates for the numerators and the determinant d_{n,m} of the 4×4 transmission system (5.9) that are stated without derivation immediately before (5.17). In particular, (5.19)–(5.20) assert the two-sided bounds |d_{n0,m}| ≈ ρ^{2n0} and |d_{n,m}| ≥ ρ^{2n0}+ρ^{2n} for all n ≥ N under the parameter choice (5.14)–(5.15). These bounds are load-bearing: the energy lower bound (5.22) and the boundedness estimate (5.25) both invoke the size of d_{n,m}. Since d_{n,m} is a multilinear combination of Bessel and Hankel functions and derivatives, an undetected polynomial-in-n factor or an extra ρ^n factor would invalidate the construction. The authors should either supply complete error-controlled derivations or state the determinant estimates as explicit lemmas with proofs.","section":"§5, (5.17)–(5.20)"},{"comment":"Equation (4.18) writes ~ψ_{1,n0,m} = C(μ̂ + μ + q_{1,n0}) with q_{1,n0} = O(1/n0), but the constant C is dropped in the subsequent estimate (4.21). The claimed choice μ̂ = −μ + i/M + p_{1,n0} satisfying (4.10) only makes the parenthesis of order O(1/M); it does not imply ℑ(μ̂)/|~ψ_{1,n0,m}|² ≥ M unless C is bounded away from zero and absorbed into the O-term. This gap affects the proof that condition (4.8) is achievable and should be repaired.","section":"§4, (4.18)–(4.21)"},{"comment":"The theorem is proved only for Newtonian potentials consisting of T_n^m components with n ≥ N, as in (5.7). The statement of Theorem 5.1 and the surrounding discussion advertise CALR for general sources, but no argument is supplied for low-order T_n^m modes or for the I_n^m and N_n^m components that appear in the full expansion (4.7). Remark 5.2 asserts that this is 'just a technical issue', but the proof of the boundedness condition (1.11) for a general source is not given. The claimed general-source CALR and the critical-radius dichotomy are therefore not established by the present proof.","section":"Theorem 5.1 and Remark 5.2"},{"comment":"The energy lower bound (5.22) is asserted without derivation from the representation (5.21), and the displayed limsup in (5.23) does not transparently follow from (5.7): for a source at distance r_s the coefficient growth in (5.7) should be governed by k r_s, whereas (5.23) involves the reciprocal of k r* and contains a square root. The subsequent combination of (5.22) and (5.23) to prove E(u) > M is therefore not justified as written; this step needs a detailed derivation.","section":"§5, (5.22)–(5.23)"}],"minor_comments":[{"comment":"The manuscript contains numerous typographical errors, including 'ovelrine' for 'overline', 'indentities' for 'identities', 'formual' for 'formula', 'hod' for 'hold', and repeated 'the the'; a careful proofreading is needed.","section":"Throughout"},{"comment":"The notation P_{λ̂/μ̂,1}(u,u) is used without being defined; please specify the Lamé parameters and the meaning of the subscript 1.","section":"§4, (4.15)"},{"comment":"The shorthand jn0i, jn1i, etc. is easy to misread because the subscripts mix the mode number n with the radius labels i and e; a table or clearer notation (e.g., j_n(k_s r_i)) would improve readability.","section":"§5, (5.2)"},{"comment":"In the 'furthermore' part of Theorem 4.1, it should be stated explicitly that n0 is assumed large enough for the asymptotic (2.8) and that the choice (4.9) depends on both n0 and M; the current phrasing 'p_{1,n0} should satisfy' is imprecise.","section":"Theorem 4.1"},{"comment":"The caption 'The absolute value of the LHS quantity in (4.10)' refers to a condition, not a quantity; it should specify what is plotted, presumably p_{1,n0} + q_{1,n0} as a function of p_{1,n0}.","section":"Figure 2"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing to know: this is a serious extension of the elastic Neumann–Poincaré spectral program, and the finite-frequency spectral system in Theorem 3.2 is a genuinely useful result. The CALR construction is a real step beyond the quasi-static analyses in [16] and [24], and the authors are honest about the source restriction. But the proof of the headline cloaking theorem is not fully written: the determinant estimates in (5.17)–(5.20) are asserted without derivation, and there is a concrete slip in Theorem 4.1 where a multiplicative constant is dropped. Treat the cloaking result as conditional, not as a closed proof.\n\nWhat is new and good: the complete spectral system for the elastic N-P operator on a sphere at finite frequency. The derivation uses vector spherical harmonics and the jump relations, and the structure — the T_n^m family decoupling and the I/N families mixing through a 2x2 block — is plausible. Taking the quasi-static limit recovers the static result of [16], which is the right sanity check. On the application side, the construction needs only one of the two strong convexity conditions to fail, and the no-resonance statement outside the critical radius is the expected physical behavior. The paper also flags its own limitations in Remark 5.2 rather than hiding them. The self-citations are appropriate: the spectral system builds directly on the authors' verified scalar result [23], and the lineage to [16] and [24] is clear.\n\nSoft spots: (i) In Theorem 4.1, (4.18) writes \\tilde\\psi = C(\\hat\\mu+\\mu+q). Then (4.21) claims the choice \\hat\\mu = -\\mu + i/M + p with p+q = O(1/M) gives the resonance inequality. That only follows if |C| <= 1 in the relevant regime, which is not established and is not obviously true. (ii) In Theorem 5.1, the estimates (5.17)–(5.20) on the determinants and numerators are load-bearing. The energy lower bound (5.22) and the boundedness estimate (5.25) both depend on the size of d_{n,m}. The paper says 'direct calculations show', but no error control in n and \\rho is provided. If the polynomial-in-n factors are off, the critical radius r* would shift. (iii) The source is restricted to the expansion (5.7), so Remark 5.2's 'general source' statement is asserted rather than proved. All three are fixable, and none makes me suspect the approach is wrong.\n\nWho it is for: researchers in plasmonic/polaritonic cloaking and spectral theory of layer potentials. The spectral system alone is worth having and worth citing. No code is shipped, so the numerical panels are illustrative.\n\nRecommendation: send to peer review. The paper deserves referee time, but the referees should push for the omitted asymptotics and the constant in Theorem 4.1. If those are supplied, the result stands.","headline":"Finite-frequency elastic N-P spectral system is a real contribution, but the CALR proof as written is conditional on unprovided determinant estimates and a dropped constant in Theorem 4.1.","tokens_in":24663,"tokens_out":4188,"would_cite":true,"duration_ms":38780,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35R30","35B30","35Q60","47G40"],"pacs":[],"model":"deepseek-v4-flash","headline":"Theorem proves elastic invisibility cloak at finite frequencies.","keywords":["anomalous localized resonance","polariton resonance","Neumann-Poincaré operator","finite frequency","beyond quasi-static limit","core-shell structure","negative material","elastic cloaking"],"falsifier":"Place a point force inside the critical radius $r_*$ so that its multipole expansion necessarily contains low-order vector spherical harmonics, and set the shell parameters as in (5.14)-(5.15). Compute the dissipation energy $E(u)$ numerically as the imaginary part of $\\hat\\mu$ is sent to zero: if $E(u)$ does not exceed every prescribed $M$ while the exterior field stays bounded, then the claim in Remark 5.2 that the restriction (5.7) is only technical would be false.","tokens_in":1777,"feed_emoji":"🛡️","tokens_out":2731,"duration_ms":87749,"temperature":0.7,"pith_summary":"This paper proves that the cloaking effect known as anomalous localized resonance, previously established for elastic waves only in the static or quasi-static regime, also occurs at finite frequencies. Working with a spherical core-shell-matrix structure in which the shell is a metamaterial with a chosen negative Lamé parameter, the authors show that a source placed inside a critical radius $r_* = \\sqrt{r_e^3/r_i}$ triggers a resonance that dissipates unbounded energy while the displacement field stays bounded outside the structure, so the source and structure become invisible. If the source lies outside that radius, no resonance occurs. The argument rests on a new complete spectral decomposition of the elastic Neumann-Poincaré operator at finite frequency, which reduces the resonance condition to a single mode whose amplitude is controlled by a small determinant.","feed_headline":"Theorem proves elastic invisibility cloak at finite frequencies","feed_subtitle":"A shell with one negative constant hides sources inside a critical radius, even at wavelengths comparable to its size.","key_machinery":"The central object is the complete spectral system of the finite-frequency Neumann-Poincaré operator, the boundary integral operator whose eigenvalues govern polariton resonances, derived on a sphere in Theorem 3.2. It yields three infinite families of eigenfunctions: the tangential vector spherical harmonics $T_n^m$ with eigenvalue $\\lambda_{1,n}$, and two mixed families $U_n^m$ and $V_n^m$ built from the normal-type harmonics $I_{n-1}^m$ and $N_{n+1}^m$, with eigenvalues $\\lambda_{2,n}$ and $\\lambda_{3,n}$ expressed through combinations of spherical Bessel and Hankel functions. This spectral decomposition reduces the core-shell transmission problem to a $4\\times 4$ algebraic system per spherical mode, and the resonance mechanism is the near-vanishing of the determinant $d_{n,m}$ for a selected large mode $n_0$ once the shell shear modulus is tuned appropriately. The critical radius $r_*$ emerges from comparing the decay rate of the source's multipole coefficients with the growth of the Bessel functions inside the shell, so that inside $r_*$ the energy lower bound grows faster than any prescribed $M$.","core_discovery":"The central discovery is that polariton resonance and cloaking due to anomalous localized resonance in elasticity do not require the quasi-static approximation. There exist core-shell-matrix configurations, with the shell shear modulus tuned to a slightly imaginary negative value of the form $\\hat\\mu = -\\mu + i\\rho^{n_0} + p_{2,n_0}$ satisfying the smallness condition (5.15), such that any source whose Newtonian potential consists of high-order vector spherical harmonic modes and lies inside the critical radius $r_* = \\sqrt{r_e^3/r_i}$ produces unbounded energy dissipation while the exterior displacement field remains bounded. Outside $r_*$ the energy stays bounded, so no resonance occurs. The proof uses the complete finite-frequency spectral system of the Neumann-Poincaré operator: for each spherical mode the eigenvalues are explicit combinations of Bessel and Hankel functions, and resonance is driven by a single large mode $n_0$ whose algebraic denominator $d_{n_0,m}$ behaves like $\\rho^{2n_0}$ and can be made small by the parameter choice.","pith_inferences":["The proof restricts the source to high-order $T_n^m$ vector spherical harmonic modes and states in Remark 5.2 that this is a technical issue, but no argument covers low-order modes or the $I_n^m$/$N_n^m$ components; numerical tests with generic point sources inside $r_*$ would reveal whether the advertised general-source cloaking actually holds.","The critical radius formula $r_* = \\sqrt{r_e^3/r_i}$ is derived from Bessel asymptotics in the large-mode limit, so at moderate frequencies or for low-order sources the effective cloaking boundary may shift or develop frequency corrections.","Because only the shear modulus needs to be negative, shells made of anisotropic or fluid-like materials that violate one convexity condition might realize the same cloak, widening the class of physical realizations.","The mode-by-mode determinantal condition suggests a tunable resonance: by choosing which $n_0$ is amplified through the imaginary part of $\\hat\\mu$, one could in principle select which source frequencies are cloaked."],"forward_implications":["Elastic cloaking via anomalous localized resonance works at frequencies comparable to the structure size, not just in the static or small-wavelength regimes.","The metamaterial design only requires tuning a single Lamé parameter, the shear modulus, leaving the other parameter free and relaxing the conditions for constructing the shell.","The critical radius gives a precise spatial threshold: sources inside it are cloaked, sources outside it generate no resonance and remain observable.","The complete finite-frequency spectral system provides a tool for analyzing other resonance-based elastic phenomena, such as field enhancement or superlensing-like effects.","Prior static and quasi-static polariton resonance constructions are recovered as limits of this spectral framework, so the present analysis unifies them."],"supporting_citations":[{"why":"Supplies the spectral system of the scalar single-layer potential used as the building block for the elastic Neumann-Poincaré spectral system in Lemma 3.1.","marker":"[23]"},{"why":"Contains the static spectral system of the elastic Neumann-Poincaré operator that Theorem 3.2 extends to finite frequencies via the quasi-static limit.","marker":"[16]"},{"why":"Provides the prior construction of polariton resonance with finite frequencies, which the present work generalizes by requiring only one convexity condition to be violated.","marker":"[24]"},{"why":"Introduced cloaking due to anomalous localized resonance and supplies the phenomenon which the paper re-derives for the elastic finite-frequency setting.","marker":"[28]"},{"why":"Establishes the rigorous framework for anomalous localized resonance through Neumann-Poincaré spectral analysis, which the paper transfers to elasticity.","marker":"[3]"},{"why":"Gives the fundamental solution, radiation condition, and traction formulas for the elastic system that define the single-layer potential and the transmission problem.","marker":"[20]"},{"why":"Provides the large-order asymptotic expansions of spherical Bessel and Hankel functions used throughout the resonance estimates in Theorems 4.1 and 5.1.","marker":"[15]"}],"fun_headline_variants":["Elastic cloaking works beyond quasi-static limit","Finite-frequency resonance enables elastic invisibility","Polariton resonances cloak at finite frequencies","No quasi-static needed for elastic resonance cloaking","Theoretical proof: elastic invisibility at finite freq"],"cache_read_input_tokens":26624,"weakest_assumption_plain":"The CALR proof assumes the source's Newtonian potential consists only of high-order $T_n^m$ vector spherical harmonic modes with $n \\geq N$; the paper calls this restriction technical in Remark 5.2, but gives no argument covering sources with low-order or $I_n^m$/$N_n^m$ components, so the advertised general-source behavior is not established.","fun_headline_variants_meta":{"raw":{"variants":["Elastic cloaking works beyond quasi-static limit","Finite-frequency resonance enables elastic invisibility","Polariton resonances cloak at finite frequencies","No quasi-static needed for elastic resonance cloaking","Theoretical proof: elastic invisibility at finite freq"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000344,"raw_usage":{"total_tokens":1864,"prompt_tokens":891,"completion_tokens":973,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":507,"completion_tokens_details":{"reasoning_tokens":902}},"tokens_in":507,"tokens_out":973,"duration_ms":9387,"temperature":1.0,"reasoning_tokens":902,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:26:32.519496+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Place a point force inside the critical radius $r_*$ so that its multipole expansion necessarily contains low-order vector spherical harmonics, and set the shell parameters as in (5.14)-(5.15). Compute the dissipation energy $E(u)$ numerically as the imaginary part of $\\hat\\mu$ is sent to zero: if $E(u)$ does not exceed every prescribed $M$ while the exterior field stays bounded, then the claim in Remark 5.2 that the restriction (5.7) is only technical would be false.","supporting_citations":[{"cited_title":"Li and H","cited_arxiv_id":null,"evidence_quote":"Supplies the spectral system of the scalar single-layer potential used as the building block for the elastic Neumann-Poincaré spectral system in Lemma 3.1."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Contains the static spectral system of the elastic Neumann-Poincaré operator that Theorem 3.2 extends to finite frequencies via the quasi-static limit."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the prior construction of polariton resonance with finite frequencies, which the present work generalizes by requiring only one convexity condition to be violated."},{"cited_title":"Milton and N.-A.P","cited_arxiv_id":null,"evidence_quote":"Introduced cloaking due to anomalous localized resonance and supplies the phenomenon which the paper re-derives for the elastic finite-frequency setting."},{"cited_title":"Ammari, G","cited_arxiv_id":null,"evidence_quote":"Establishes the rigorous framework for anomalous localized resonance through Neumann-Poincaré spectral analysis, which the paper transfers to elasticity."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the fundamental solution, radiation condition, and traction formulas for the elastic system that define the single-layer potential and the transmission problem."},{"cited_title":"Colton and R","cited_arxiv_id":null,"evidence_quote":"Provides the large-order asymptotic expansions of spherical Bessel and Hankel functions used throughout the resonance estimates in Theorems 4.1 and 5.1."}],"review_version":1}