{"id":"131f8530-c5a0-4c52-a604-c8fe06e5c9fb","arxiv_id":"2506.01670","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Multicontinuum splitting schemes for the wave equation with high-contrast coefficients yield partially explicit time discretizations with contrast-independent stability bounds.","lead":"The paper constructs time-stepping schemes that split a high-contrast wave problem into fast and slow parts, treating the fast part implicitly and the slow part explicitly. The result is a method whose stability bound is claimed to be independent of the material contrast, which could make large multiscale wave simulations cheaper.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Contrast-independence of the stability condition is asserted, not proved; the wave-specific c22 term enters scheme 2, and the cited flow preprint [37] does not obviously cover it.","rationale":"The reader's weakest assumption identifies the same load-bearing point: the paper defers the contrast-independence of the relevant norms to the unpublished, overlapping-author preprint [37]. My reading confirms that this is the single most load-bearing concern. The time-stepping algebra in Theorem 1 is coherent, and the stability inequalities of Theorem 2 would be standard if the contrast-independence of m22, a22, and c22 and the behavior of gamma were established. But the paper does not prove these properties, and the numerical experiments do not test them by varying contrast. The concern is not that the schemes are unstable; it is that the advertised contrast-free time-step restriction is not supported within this paper. Since the issue is a gap in support rather than a demonstrated contradiction, the conditional verdict is appropriate. No change to the reader's verdict is needed.","tokens_in":17142,"tokens_out":15089,"duration_ms":178185,"concrete_test":"Reproduce the cell problems (2.7)-(2.8) for the layered geometry of Example 1 on a fixed coarse and fine mesh, with kappa_high/kappa_low = 10^2, 10^4, 10^6, and 10^8. Assemble M22, A22, C22 from (3.5) and compute lambda1 = smallest generalized eigenvalue of A22 x = lambda M22 x and lambda2 = smallest generalized eigenvalue of (A22 + C22) x = lambda M22 x on V2,H. If lambda1 or lambda2 changes by more than about 10% across this contrast range, the claimed contrast-independent stability condition (3.26)/(3.29) fails. Also evaluate gamma in (3.24) for the same range and check that 2(1-gamma^2) does not degrade with contrast.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central advertised benefit is that the stability conditions (3.26) and (3.29) are contrast-independent when fast modes are assigned to Vmc,1,H and slow modes to Vmc,2,H. The only support for this is the sentence after Theorem 3: the m22-, a22-, and c22-norms are 'all independent of the contrast ... thanks to the construction of the cell problems [37]'. No estimate of these norms as functions of kappa_max/kappa_min is given, and [37] is an overlapping-author preprint on flow problems; the paper gives no indication that it proves contrast-independence for c22, which is exactly the new term entering the scheme-2 condition (3.29). In addition, the stability conditions contain the coupling constant gamma from (3.24), which is not shown to be contrast-independent; a gamma approaching 1 as contrast grows would degrade (1-gamma^2) even if the Rayleigh ratios are constant. The numerical section does not close this gap: all examples use a single contrast ratio 10^3 and a single time step tau = 1e-3, and no maximum stable time step is reported as a function of contrast. Thus the load-bearing claim currently rests on an unverified deferred assertion.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes two three-level partially explicit time discretization schemes, (3.8) and (3.9), for the high-contrast wave equation, built on a multicontinuum homogenization space decomposition. Fast continua are updated implicitly and slow continua explicitly. A discrete energy identity is proved for scheme 1 (Theorem 1), and stability is derived as condition (3.26) (Theorem 2); an analogous stability condition (3.29) for scheme 2 is stated as Theorem 3 without proof. The paper claims that, with appropriately chosen continua, the stability conditions are independent of the coefficient contrast, citing the cell-problem construction in the overlapping-author preprint [37]. An optimized decomposition via a tensor Rayleigh quotient and a simplified generalized eigenvalue problem is discussed in Section 4. Numerical examples for layered and point fields with two and three continua show that the proposed schemes match the accuracy of the fully implicit scheme for contrast 10^3.","tokens_in":17402,"tokens_out":11547,"duration_ms":107729,"significance":"If the contrast-independence assertion were proved, the schemes would provide a practically useful way to remove the high-contrast CFL restriction while retaining accuracy comparable to implicit methods. The self-contained algebraic derivations of Theorem 1 and Theorem 2 are a genuine contribution, as is the formulation of the stability conditions in terms of the explicitly treated subspace. However, the headline property is not established within the manuscript: it is deferred to an unpublished preprint, and the numerical section does not vary the contrast. The paper is therefore conditionally acceptable after the central gap is closed.","major_comments":[{"comment":"The claim that the stability conditions (3.26) and (3.29) are contrast-independent is the central advertised result, but the only support is the sentence following Theorem 3 stating that the m22-, a22-, and c22-norms are 'all independent of the contrast ... thanks to the construction of the cell problems [37]'. No estimate of these norms as functions of κmax/κmin is given, and [37] is an overlapping-author preprint on flow problems that does not obviously cover the wave-specific c22 term entering (3.29). Moreover, the constant γ in (3.24) appears in both conditions through the factor (1−γ2), and no argument shows that γ is bounded away from 1 uniformly in the contrast. Please provide a proof, or at minimum a quantitative estimate, of contrast-independence for all quantities appearing in (3.26) and (3.29).","section":"Section 3, after Eq. (3.29)"},{"comment":"Theorem 3 is stated without proof: the text says only that 'similar energy conservation can be established' for scheme 2. Since scheme 2 is one of the two proposed time discretizations and its stability condition (3.29) is used in the numerical section, the energy identity and the stability argument should be written out in full or explicitly reduced to the proof of Theorem 2 with the c22 terms accounted for.","section":"Theorem 3"},{"comment":"The optimized decomposition section defers its central assertions to [37]: 'It can be shown that γ = 0' and 'The same stability results hold under certain assumption' are not proved, and the localization assumption that 'the decompositions for different coarse blocks are similar' is stated without justification. Because Section 4 is presented as a way to relax the stability conditions and decouple the schemes, these claims need to be either proved in the manuscript or clearly stated as assumptions with supporting numerical evidence.","section":"Section 4"},{"comment":"The numerical experiments do not test the contrast-independence claim. All examples use the single contrast ratio 10^3 and the single time step τ = 10^-3, and no maximum stable time step is reported as a function of κmax/κmin. To support the headline property, please add experiments with varying contrast (for example 10^2, 10^4, 10^5) and report the largest stable τ for each case, or otherwise verify empirically that the stability conditions are contrast-independent.","section":"Section 5"}],"minor_comments":[{"comment":"The sentence 'the variation in material properties ... necessities extremely fine spatial discretization' should read 'necessitates', and the sentence is grammatically incomplete.","section":"Section 1"},{"comment":"The first line, 'Throughout this work, We begin by summarizing...', contains a stray comma and an uppercase 'We' in the middle of the sentence.","section":"Section 3"},{"comment":"The definition of the Rayleigh quotient contains unmatched parentheses, and the contraction A : (v⊗w) is not defined explicitly.","section":"Equation (4.4)"},{"comment":"The sentence 'requiring less computationally effort' should be 'requiring less computational effort'.","section":"Section 6"},{"comment":"The tables report eigenvalues and eigenvectors but not the resulting stability restriction τmax; reporting this value would allow readers to check conditions (3.26) and (3.29) directly.","section":"Tables 1-4"}],"recommendation":"major_revision","confidential_remarks":"The principal unresolved issue is the reliance on the authors' own preprint [37] for the contrast-independence property that is the paper's headline. This is not a question of correctness of Theorem 1, but of the completeness of the contribution for a journal publication. The editor may also wish to check that the overlap with [37] and with the existing partially explicit wave papers [34,35] is sufficiently differentiated."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nThe honest headline: this is a real extension of the authors' flow work to waves, with a self-contained energy proof for the first scheme, but the paper's main selling point—contrast-independent stability—is asserted, not proved, and the numerical section is too thin to fill the gap.\n\nWhat's actually new: the three-layer schemes (3.8) and (3.9), the discrete energy (3.12), and the stability conditions (3.26) and (3.29). Theorem 1's energy conservation proof is algebraically consistent and, as far as I can tell, correct. Theorem 2 follows from it in a standard way. That is a solid core and worth building on.\n\nThe soft spots are all around the contrast-independence claim. After Theorem 3, the authors state that the m22, a22, and c22 norms are 'all independent of the contrast ... thanks to the construction of the cell problems [37]', and refer to [37] for the proofs. [37] is an overlapping-author preprint on flow problems; this paper gives no indication that it covers the wave-specific c22 term in (3.29). Nor is the coupling constant gamma in (3.24) shown to be contrast-independent, so (1−γ^2) might degrade even if the Rayleigh ratios are constant. The numerical examples do not close this: every test uses contrast 10^3 and τ=10^-3, with no maximum stable step as a function of contrast, no error magnitudes, and no runtimes. The optimized decomposition section similarly defers to [37] and relies on an unverified localization assumption.\n\nThese are addressable gaps, not fatal flaws. The scheme itself is coherent and the energy proof is real. But the advertised benefit is the contrast-independent stability, and that is currently a promise to the literature rather than a result. A serious referee should ask for a proof or a systematic numerical study of the norms and γ as functions of contrast, plus convergence rates and timings.\n\nWho should read this: anyone working on multiscale wave methods with high contrast. I would take it to a reading group and would cite the scheme and energy proof, with a caveat. It deserves peer review, not desk rejection.","headline":"A useful, mostly correct scheme paper whose headline claim of contrast-independent stability is deferred to an unpublished preprint and not closed by the numerics.","tokens_in":17907,"tokens_out":2914,"would_cite":true,"duration_ms":29126,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["65M12","65M60","35L05"],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper constructs two partially explicit time-stepping schemes for the high-contrast wave equation that split fast and slow dynamics so that the stability condition is claimed to be independent of the coefficient contrast.","keywords":["multicontinuum homogenization","partially explicit time discretization","high-contrast wave equation","discrete energy conservation","contrast-independent stability","multiscale basis functions","generalized eigenvalue decomposition","wave propagation in heterogeneous media"],"falsifier":"In the layered two-continuum example (5.3), compute the stability ratio $\\inf_W \\|W\\|_{m_{22}}^2 / \\|W\\|_{a_{22}}^2$ for contrast values $10$, $10^3$, and $10^6$ while keeping the coarse mesh fixed; if this ratio shrinks as the contrast grows, the stability condition is contrast-dependent and the central claim fails.","tokens_in":16929,"feed_emoji":"🌊","tokens_out":6773,"duration_ms":67308,"temperature":0.7,"pith_summary":"Multicontinuum homogenization represents the solution of the high-contrast wave equation through local averages on several continua, and the paper's central move is to split those continua into a fast, high-value group and a slow, low-value group. On top of that split the authors build two three-layer time integration schemes that solve the fast group implicitly and the slow group explicitly. The main claim is that the resulting stability conditions, (3.26) and (3.29), bound the time step only by the slow group's energy norms, and that with properly chosen continua those norms, and hence the time step, are independent of the coefficient contrast. The numerical examples show the schemes matching the accuracy of the fully implicit discretization while the ordinary explicit scheme diverges with the same step size; the practical stake is wave simulation in high-contrast media at explicit-like cost without contrast-driven time-step restrictions.","feed_headline":"Splitting fast and slow modes frees wave time steps from contrast","feed_subtitle":"A multicontinuum decomposition makes 1000:1 wave simulations nearly as cheap as explicit schemes.","key_machinery":"The machinery is the direct-sum decomposition $V_{\\mathrm{mc},H}=V_{\\mathrm{mc},1,H}\\oplus V_{\\mathrm{mc},2,H}$ of the discrete multicontinuum space, together with the bilinear forms $m_{ij}$, $a_{ij}$, and $c_{ij}$ defined through downscaling operators $T_j$. Each continuum is represented by multiscale basis functions obtained from local energy-minimizing cell problems with constraints on the averages of the solution; choosing indices $I_1$ for fast/high-value continua and $I_2$ for slow/low-value continua puts all contrast in the implicitly treated block. The stability proof runs through an exactly conserved discrete energy $E^{n+1/2}$; a strengthened Cauchy-Schwarz constant $\\gamma$ controls cross-coupling, and the final conditions bound $\\tau^2$ by $2(1-\\gamma^2)$ times the $m_{22}$-to-$a_{22}$ (or $m_{22}$-to-$(a_{22}+c_{22})$) ratio of the slow subspace, which is exactly the quantity the cell-problem construction is claimed to make contrast-independent.","core_discovery":"On the paper's own terms, the central discovery is that the standard multicontinuum homogenization expansion for the wave equation can be organized into two subspaces, one carrying the fast (high-contrast) dynamics and one the slow dynamics, so that the fully implicit semidiscrete problem can be replaced by the partially explicit three-layer schemes (3.8) and (3.9). The schemes conserve a discrete energy exactly (Theorem 1) and are stable under the conditions (3.26) and (3.29), which involve only the mass and stiffness norms of the explicitly treated slow subspace. Since, by the construction of the cell problems, those norms are claimed to be independent of the contrast for continua placed in low-value regions, the stability conditions are contrast-independent, and the method is presented as reaching essentially the accuracy of the fully implicit discretization at lower cost. The paper also proposes an optimized decomposition, through a tensor Rayleigh quotient problem and a cheaper generalized eigenvalue problem, to relax the stability restriction further.","pith_inferences":["If the contrast-independence proof transfers from the flow cell problems to the wave cell problems, the same splitting strategy should extend to other second-order hyperbolic systems, such as elastic waves or acoustics in fractured media, where high-contrast coefficients appear.","The tensor Rayleigh quotient formulation could be turned into an automatic continuum-selection tool for coefficient fields whose high- and low-value regions are not known in advance; the paper itself only treats predefined continua.","A natural testable extension is to push the contrast beyond $10^3$, the largest value reported here, and record whether the observed critical time step stays flat; under the paper's claim it should.","One could also test the optimized decomposition on coefficient fields with more than three continua, where the generalized eigenvalue problem would have to choose a genuinely mixed slow subspace rather than simply the low-value regions."],"forward_implications":["For two-value coefficients, choosing the high-value continuum as the implicit component and the low-value continuum as the explicit component makes the time-step restriction independent of the ratio $\\kappa_{\\max}/\\kappa_{\\min}$, so the same $\\tau$ is expected to work as the contrast is raised.","The partially explicit schemes cost less per step than the fully implicit scheme because only the fast block requires an implicit solve, while the numerical tests show errors nearly identical to the implicit reference solution.","When the two subspaces are $L^2$-orthogonal, scheme 2 decouples, so the explicit part can be updated independently of the implicit solve; the paper notes that mass lumping can remove the remaining coupling in scheme 1.","The optimized decomposition based on a generalized eigenvalue problem selects continua by taking linear combinations of the original ones, and the resulting stability restriction scales roughly like $H\\lambda_{i_0}^{-1/2}$, which relaxes the step-size bound without increasing computational cost."],"supporting_citations":[{"why":"defines the multicontinuum expansion and cell problems that produce the macroscopic variables and basis functions used throughout the paper.","marker":"[30]"},{"why":"provides the general multicontinuum homogenization theory underlying the space decomposition into continua.","marker":"[31]"},{"why":"introduces contrast-independent partially explicit time discretizations for flow problems, the template that the wave schemes extend.","marker":"[34]"},{"why":"gives the existing contrast-independent partially explicit scheme for wave problems whose stability strategy is adapted here.","marker":"[35]"},{"why":"is the companion paper whose cell-problem construction is cited for the key claim that the slow-continuum norms are contrast-independent.","marker":"[37]"},{"why":"supplies the strengthened Cauchy-Schwarz inequality used to define the coupling constant $\\gamma$ in the stability conditions.","marker":"[39]"}],"fun_headline_variants":["Fast/slow wave split removes contrast limit for explicit schemes","Splitting wave modes yields contrast-free stability conditions","Multicontinuum splitting breaks contrast barrier in wave simulations","Contrast-independent time steps from fast/slow wave decomposition","Wave solver: explicit slow, implicit fast, no contrast drag"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The headline result rests on an unproved assertion, deferred to a companion paper, that the cell-problem construction makes the energy stored in the low-coefficient continua insensitive to how large the high coefficient is.","fun_headline_variants_meta":{"raw":{"variants":["Fast/slow wave split removes contrast limit for explicit schemes","Splitting wave modes yields contrast-free stability conditions","Multicontinuum splitting breaks contrast barrier in wave simulations","Contrast-independent time steps from fast/slow wave decomposition","Wave solver: explicit slow, implicit fast, no contrast drag"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000904,"raw_usage":{"total_tokens":3884,"prompt_tokens":933,"completion_tokens":2951,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":549,"completion_tokens_details":{"reasoning_tokens":2869}},"tokens_in":549,"tokens_out":2951,"duration_ms":24639,"temperature":1.0,"reasoning_tokens":2869,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T11:37:28.500601+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"In the layered two-continuum example (5.3), compute the stability ratio $\\inf_W \\|W\\|_{m_{22}}^2 / \\|W\\|_{a_{22}}^2$ for contrast values $10$, $10^3$, and $10^6$ while keeping the coarse mesh fixed; if this ratio shrinks as the contrast grows, the stability condition is contrast-dependent and the central claim fails.","supporting_citations":[{"cited_title":"Multicontinuum homogenization and its relation to nonlocal multicontinuum theories","cited_arxiv_id":null,"evidence_quote":"defines the multicontinuum expansion and cell problems that produce the macroscopic variables and basis functions used throughout the paper."},{"cited_title":"Chung, Y","cited_arxiv_id":null,"evidence_quote":"provides the general multicontinuum homogenization theory underlying the space decomposition into continua."},{"cited_title":"Contrast-independent partially explicit time discretizations for multiscale flow problems.Journal of Computational Physics, 445:110578, 2021","cited_arxiv_id":null,"evidence_quote":"introduces contrast-independent partially explicit time discretizations for flow problems, the template that the wave schemes extend."},{"cited_title":"Contrast-independent partially explicit time discretizations for multiscale wave problems","cited_arxiv_id":null,"evidence_quote":"gives the existing contrast-independent partially explicit scheme for wave problems whose stability strategy is adapted here."},{"cited_title":"Strengthened cauchy-schwarz and h¨ older inequalities","cited_arxiv_id":null,"evidence_quote":"supplies the strengthened Cauchy-Schwarz inequality used to define the coupling constant $\\gamma$ in the stability conditions."}],"review_version":1}