{"id":"d7776a02-6b5f-46fe-af66-6b0b37955d9d","arxiv_id":"2608.07072","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"In periodic two-layer media the effective mass density depends on the shear modulus, transitioning from the isotropic elastic arithmetic average to the anisotropic acoustic values (arithmetic normal, harmonic parallel) once band-gap intervals are excluded.","lead":"This paper derives an analytical expression for the effective mass density of a periodic two-layer elastic material as a function of the shear modulus, showing how that density changes as the material softens from elastic to acoustic behavior. It also maps out band gaps and exceptional points in shear-modulus space where no well-defined effective density exists.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The acoustic-limit conclusion is a pass-band-filtered statement, not a genuine limit, and the claimed domain of validity of Eq. (4.21) is internally contradicted by the paper's own avoided-crossing exclusions.","rationale":"The reader's weakest assumption already identifies the pass-band-filtered nature of the acoustic-limit conclusion and the accumulation of excluded intervals as μ→0. My reading of the manuscript confirms this and adds a sharper internal inconsistency: the conclusions claim validity of Eq. (4.21) for all real K_p, but Section 4.3 explicitly excludes avoided-crossing neighbourhoods from the well-defined effective-density domain, even though K_p can be real there. This is not a disagreement with the reader; it strengthens the same concern with a concrete contradiction. The transfer-matrix/Bloch derivation itself appears standard, and the recovery of classical acoustic and elastic effective densities in Sections 4.1–4.2 is a good internal check. The self-consistent validation in Figure 17 supports Eq. (4.21) in well-behaved pass bands. However, the central transition claim depends on the excluded-set filtering, so the conditional verdict is appropriate: the derivation is defensible, but the stated domain and the unqualified 'tends to' conclusion need revision. My proposed test directly checks the domain claim at an avoided crossing, where the paper's own text predicts failure. I therefore recommend no change to the reader's conditional verdict.","tokens_in":28377,"tokens_out":4327,"duration_ms":40994,"concrete_test":"Using the paper's own BGIIb parameters, locate the avoided-crossing value μ_AC where the eigenvector components cross (as in Figure 14), verify that the compressional Bloch wavenumber K_p(μ_AC) is real, and compute ρ_x^eff from Eq. (4.21) for θ_i = π/6, π/4, and π/3. If the three values differ by more than 0.1%, the claim that Eq. (4.21) is valid for all real K_p is falsified; if they agree, that domain concern is refuted.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim is that Eq. (4.21) gives a μ-dependent effective density that, as μ→0, tends to the acoustic anisotropic value ρ_g = (φ/ρ1 + (1−φ)/ρ2)^{-1}, and the conclusions state that (4.21) is valid \"over the entire range of values of μ where the Bloch wavenumber for compressional waves is real.\" Both parts are load-bearing and both are overclaimed. In Section 4.4 the authors explicitly say the analytic μ→0 limit \"cannot be taken\" because μ-gaps, EP gaps, μ-P-gaps, and avoided-crossing neighbourhoods accumulate at μ=0; the plotted approach to ρ_g is obtained by sampling only pass bands and discarding all intervals where ρ_x^eff is complex or angle-dependent. This is a filtered statement, not a limit in the standard sense: every punctured neighbourhood of 0 contains excluded intervals, so the claim \"ρ_x^eff tends to ρ_g\" is not well-defined without specifying the filter. Moreover, the domain-of-validity claim is internally contradicted by Section 4.3, where avoided crossings (e.g., the BGIIb example of Figure 14) occur outside μ-gaps with K_p still real, yet ρ_x^eff is shown to depend on the incident angle and the authors state that \"effective properties also cannot be defined around these crossings.\" Thus Eq. (4.21) is not valid for all real K_p; acceptable μ-values must also exclude AC neighbourhoods. Since the paper's headline contribution is precisely the elasto-acoustic transition, this domain/limit gap weakens the central result as stated.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript derives a shear-modulus-dependent effective mass density for a periodic two-layer elastic medium by combining a transfer-matrix analysis with Bloch's theorem. It recovers the classical anisotropic acoustic effective density and the isotropic elastic effective density in the low-frequency regime, obtains Eq. (4.21) as an analytical expression for the effective density rho_x^eff, and studies the associated mu-gaps, exceptional points, and avoided crossings. A comparison with a 5000-cell self-consistent calculation is used to validate the expression in pass bands.","tokens_in":28646,"tokens_out":5111,"duration_ms":46609,"significance":"If the main claims are taken in the qualified sense advocated in Sections 4.3 and 4.4, the paper provides a useful analytical bridge between elastic and acoustic effective descriptions. The transfer-matrix/Bloch derivation is transparent, the low-frequency eigenvalue formulas in Appendix C are checkable, and the comparison with the dynamic self-consistent method in Figure 17 is a meaningful internal consistency check. The paper is also commendably explicit about the intervals in which an effective density cannot be defined. The two overclaims identified below concern the stated domain of validity of Eq. (4.21) and the status of the mu-to-0 limit, and they should be corrected before publication.","major_comments":[{"comment":"The Conclusion states that Eq. (4.21) is valid \"over the entire range of values of mu where the Bloch wavenumber for compressional waves is real.\" This is contradicted by Section 4.3, where avoided-crossing neighbourhoods are excluded from the admissible domain even though K_p is real there (see Figure 14 and the text following Eq. (4.24)). The domain of validity should be restated as real-K_p values outside EP gaps, mu-P-gaps, and AC neighbourhoods, or a precise characterisation of the excluded set should be given.","section":"Section 5 / Section 4.3"},{"comment":"The statement that rho_x^eff \"tends to\" the acoustic geometric average rho_g is not a limit in the standard sense. As the paper itself acknowledges, the analytic mu-to-0 limit \"cannot be taken\" because mu-gaps, EP gaps, mu-P-gaps, and avoided-crossing neighbourhoods accumulate at mu = 0; Figure 18b is obtained by sampling pass bands and discarding all excluded intervals. The claim should be reformulated as a pass-band-filtered statement, specifying the filter (for example, evaluation along pass-band midpoints) and ideally supplemented by a quantitative convergence statement.","section":"Section 4.4 / Figure 18"},{"comment":"Equation (4.21) is derived by equating the Bloch phase speed to the phase speed of a homogeneous anisotropic medium. For this procedure to define an intrinsic effective density, the resulting rho_x^eff must be independent of the incidence angle. The paper demonstrates this numerically in pass bands (Figure 16b), but the main-result statement in Section 5 does not mention this restriction. The authors should state explicitly that angle-independence is part of the definition or add a condition excluding incidence-angle-dependent regimes.","section":"Eq. (4.21)"},{"comment":"The low-frequency eigenvalue approximations in Eqs. (C.10)-(C.13) are stated without derivation. Since these formulas underpin the effective-density results in Eqs. (4.8), (4.13), (4.17), and (4.20), a derivation or a precise reference for the approximation should be included, or at least the leading-order error should be quantified.","section":"Appendix C"}],"minor_comments":[{"comment":"In Eq. (4.12), the definition rho_a = phi*rho_1 + (1-phi) is missing the factor rho_2; it should read phi*rho_1 + (1-phi)*rho_2.","section":"Eq. (4.12)"},{"comment":"In Eq. (4.19), the subscript p in phi_gamma_p appears to be a typo: the expression is derived from the shear wavenumber k_s and should be phi_gamma_s.","section":"Eq. (4.19)"},{"comment":"The notation gamma_1^* = 1/gamma_1 is used for reciprocal pairs, but the asterisk conventionally denotes complex conjugation; these two notions coincide only when |gamma| = 1. The notation should be clarified.","section":"Section 3.1"},{"comment":"The figure axis labels should specify clearly that the horizontal axis is mu/lambda, as in the caption, to avoid confusion with dimensional mu.","section":"Figure 16"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal and the underlying approach is sound, but the stated domain of validity and the acoustic-limit claim need to be qualified carefully. The two overclaims identified in the main report are fixable by rewording and by making the filtered nature of the limit explicit. The novelty claim \"for the first time\" in the Abstract should be checked against Ref. [21] and the related literature on exceptional points in laminates."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's my read. The paper extends the authors' earlier single-cell numerical study [21] to the infinite periodic case and, for the first time, gives an analytical shear-modulus-dependent effective density, Eq. (4.21). That is a real result, and the transfer-matrix/Bloch machinery is standard but competently executed. The low-frequency eigenvalue calculations in Appendix C check out, the acoustic and elastic limits reproduce the known anisotropic/isotropic densities, and the self-consistent validation with 5000 unit cells is a good consistency check. The classification of mu-gaps, EP gaps, and avoided crossings is careful and the figures support it.\n\nThe soft spots are about how the conclusions are phrased, not about the core derivation. The claim that (4.21) is valid over the entire range where the compressional Bloch wavenumber is real is contradicted by the paper's own Section 4.3: avoided crossings occur with K_p real, yet the effective density is angle-dependent there and the authors state effective properties cannot be defined. So the formula is valid only where K_p is real and one is away from AC/EP neighbourhoods. Second, the mu->0 acoustic limit is not a genuine limit: mu-gaps, EP gaps, and ACs accumulate at zero, and the paper itself says the analytic limit cannot be taken (Section 4.4, Figure 18). The plotted approach to rho_g is obtained by sampling pass bands and discarding the excluded intervals. That is a filtered statement, not a limit in the usual sense. The authors acknowledge this, but the abstract and conclusion still state the transition as if it were a clean limit.\n\nNeither issue destroys the paper. The fix is wording: qualify the domain of validity and describe the acoustic limit as a pass-band-filtered tendency rather than a limit. I would send this to a serious referee; the central analytical result is reproducible and the caveats are already latent in the text. Minor: no code or data, but the derivations are checkable and the validation is described in enough detail to reproduce.\n\nFor whom: people working on homogenization, metamaterials, and effective densities in layered media. It belongs in a good applied maths journal, after revision.","headline":"Solid extension to the periodic case with a genuinely new mu-dependent effective density, but the domain-of-validity and acoustic-limit claims are overbroad and need qualification.","tokens_in":29216,"tokens_out":2348,"would_cite":true,"duration_ms":21014,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["74J20","74Q15","35B27"],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper derives an analytical, shear-modulus-dependent effective mass density for a periodic two-layer elastic medium and shows that, in pass bands, it tends to the anisotropic acoustic effective density as shear vanishes, while band…","keywords":["effective mass density","Bloch wavenumber","transfer matrix","periodic layered media","elasto-acoustic transition","exceptional points","mu-gaps","anisotropic effective density"],"falsifier":"Evaluate Eq. (4.21) on a sequence of shear moduli tending to zero that lies entirely inside pass bands; the central claim fails if $\\rho_x^{\\rm eff}$ does not converge to $1/(\\phi/\\rho_1 + (1-\\phi)/\\rho_2)$. A complementary check is to compute transmission through a finite stack of several thousand unit cells in a narrow pass band and compare the density extracted from the scattering coefficients with Eq. (4.21) at the same parameters.","tokens_in":28144,"feed_emoji":"🌊","tokens_out":5851,"duration_ms":50407,"temperature":0.7,"pith_summary":"This paper sets out to resolve a mismatch: in a layered medium the effective mass density is isotropic in elasticity and anisotropic in acoustics, yet the standard elastic expression does not depend on the shear modulus, so taking the no-shear limit gives no information. It combines the transfer matrix method with Bloch's theorem to obtain the Bloch wavenumber of a periodic two-layer medium, and from it derives an explicit shear-modulus-dependent effective density in the direction parallel to the layers. The central conclusion is that, in the pass bands where the Bloch wavenumber is real, this effective density approaches the acoustic anisotropic values as the shear modulus tends to zero. The paper also shows that the approach to that limit is obstructed by shear-wave band gaps in shear-modulus space, exceptional points, and avoided crossings, which leave intervals where the effective density is complex, negative, or incidence-angle-dependent. A sympathetic reader should therefore read the acoustic limit as a pass-band-filtered statement, not as an unfiltered limit.","feed_headline":"New formula traces layered-media density from elastic to acoustic regime","feed_subtitle":"As shear modulus drops to zero, effective density becomes anisotropic, with singular gaps near band edges.","key_machinery":"The load-bearing object is the compressional Bloch wavenumber $K_p = \\ln(\\gamma_p)/(ih)$, computed from the eigenvalues $\\gamma$ of the $4\\times 4$ symplectic transfer matrix for the two-layer unit cell. The eigenvalues come in reciprocal pairs, so $\\eta_p$ and $\\eta_s = \\gamma + 1/\\gamma$ determine pass bands when $|\\eta|<2$, $\\mu$-gaps when $|\\eta|>2$, and exceptional points when the two pairs coalesce. Converting $K_p$ into an effective phase speed and an effective propagation angle via Snell's law, and comparing with the phase speed of a medium with anisotropic density, yields Eq. (4.21), the formula whose behaviour carries the whole paper.","core_discovery":"The central claim, stated in Eq. (4.21), is that for a periodic two-layer elastic medium with identical Lamé parameters in both layers, low-frequency oblique incidence, and perfect contact, the x-component of the effective mass density is a closed-form function of the shear modulus obtained by matching the anisotropic elastic phase-speed relation to the compressional Bloch wavenumber. The formula reproduces the isotropic elastic result $\\rho_x^{\\rm eff} = \\phi\\rho_1 + (1-\\phi)\\rho_2$ when $\\mu$ is comparable to $\\lambda$, and it approaches the acoustic harmonic average $1/\\rho_x^{\\rm eff} = \\phi/\\rho_1 + (1-\\phi)/\\rho_2$ along pass-band sequences as $\\mu \\to 0$, while $\\rho_z^{\\rm eff}$ stays at the arithmetic average for all $\\mu$. The paper explicitly notes that the analytic $\\mu \\to 0$ limit cannot be taken because $\\mu$-gaps, exceptional-point gaps, and avoided-crossing neighbourhoods accumulate as $\\mu$ tends to zero; the acoustic limit is therefore a statement about the dominant pass-band behaviour, not a limit in the ordinary sense.","pith_inferences":["Because the $\\mu$-gaps accumulate near $\\mu=0$, the phrase 'tends to' should be read as convergence along pass-band sequences; the paper's own analysis suggests an unfiltered limit does not exist.","The same transfer-matrix construction should extend to anisotropic or slightly viscoelastic layers; a testable prediction is that small material damping smooths the singularities near exceptional-point gaps and restores a well-defined effective density across the transition.","The coexistence of acoustic-like pass bands with evanescent $\\mu$-intervals suggests that a complete homogenised description of the elasto-acoustic transition should carry extra bookkeeping for forbidden intervals, not just a single effective density.","The angle-dependence of the exceptional points could be exploited experimentally: tuning the incidence angle moves the exceptional points in $\\mu$-space, which changes where the effective density breaks down."],"forward_implications":["For any incident angle inside a pass band, the effective density from Eq. (4.21) is independent of angle to within about 0.1%, and much less in the interior of pass bands.","The elastic effective density is not necessarily isotropic: at finite small $\\mu$ it becomes anisotropic before reaching the acoustic values.","Outside the $\\mu$-gaps, the analytical formula agrees with a dynamic self-consistent calculation for 5000 unit cells, and it explains why the self-consistent method fails to converge inside the gaps.","The widths of $\\mu$-gaps and adjacent pass bands both shrink as $\\mu \\to 0$, but pass bands are on average about 18 times wider, so pass-band behaviour dominates.","Inside $\\mu$-P-gaps and exceptional-point gaps the effective density can become complex or negative, so no effective medium description exists in those intervals."],"supporting_citations":[{"why":"Established the elastic-to-acoustic transition of effective mass density for a single unit cell, which this paper extends to periodic media.","marker":"[21]"},{"why":"Supplies the classical anisotropic acoustic effective density values that the $\\mu \\to 0$ pass-band limit is compared against.","marker":"[17]"},{"why":"Provides the standard result that the elastic effective density of a stratified medium is isotropic, the starting point this paper refines.","marker":"[20]"},{"why":"Dynamic generalized self-consistent model used as the numerical validation benchmark for Eq. (4.21).","marker":"[22]"},{"why":"Identified exceptional points in lossless elastic layered media, the framework used here for $\\mu$-gaps, EP gaps, and avoided crossings.","marker":"[47]"},{"why":"Supplies the transfer-matrix and symplectic eigenvalue machinery used to derive the Bloch wavenumber.","marker":"[26]"},{"why":"Early observation that acoustic effective density is dynamic and differs from the static density, motivating the whole question.","marker":"[3]"}],"fun_headline_variants":["Transition in effective density mapped in layered media","Layered media reveal density shift from elastic to acoustic","Effective mass density bridges elastic and acoustic regimes","Shear modulus zero reveals anisotropic density gaps"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The acoustic-limit conclusion depends on discarding every interval of shear modulus where the Bloch wavenumber is complex or the effective density is angle-dependent, and these intervals accumulate as the shear modulus tends to zero; without that filtering, the limit is not defined.","fun_headline_variants_meta":{"raw":{"variants":["Transition in effective density mapped in layered media","Layered media reveal density shift from elastic to acoustic","Effective mass density bridges elastic and acoustic regimes","Shear modulus zero reveals anisotropic density gaps"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000121,"raw_usage":{"total_tokens":1121,"prompt_tokens":999,"completion_tokens":122,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":615,"completion_tokens_details":{"reasoning_tokens":65}},"tokens_in":615,"tokens_out":122,"duration_ms":1963,"temperature":1.0,"reasoning_tokens":65,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T15:23:20.703154+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Evaluate Eq. (4.21) on a sequence of shear moduli tending to zero that lies entirely inside pass bands; the central claim fails if $\\rho_x^{\\rm eff}$ does not converge to $1/(\\phi/\\rho_1 + (1-\\phi)/\\rho_2)$. A complementary check is to compute transmission through a finite stack of several thousand unit cells in a narrow pass band and compare the density extracted from the scattering coefficients with Eq. (4.21) at the same parameters.","supporting_citations":[{"cited_title":"N´ u˜ nez, W","cited_arxiv_id":null,"evidence_quote":"Established the elastic-to-acoustic transition of effective mass density for a single unit cell, which this paper extends to periodic media."},{"cited_title":"Schoenberg and P","cited_arxiv_id":null,"evidence_quote":"Supplies the classical anisotropic acoustic effective density values that the $\\mu \\to 0$ pass-band limit is compared against."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the standard result that the elastic effective density of a stratified medium is isotropic, the starting point this paper refines."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Dynamic generalized self-consistent model used as the numerical validation benchmark for Eq. (4.21)."},{"cited_title":"Alizadeh and A","cited_arxiv_id":null,"evidence_quote":"Identified exceptional points in lossless elastic layered media, the framework used here for $\\mu$-gaps, EP gaps, and avoided crossings."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the transfer-matrix and symplectic eigenvalue machinery used to derive the Bloch wavenumber."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Early observation that acoustic effective density is dynamic and differs from the static density, motivating the whole question."}],"review_version":1}