{"id":"72da93cb-9ab7-4727-a07a-652bdb3833ba","arxiv_id":"2508.18754","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"For a vector-valued Allen-Cahn equation with two sphere wells, the ε→0 limit is mean curvature flow coupled to harmonic map heat flow into S^{n-1} with a mixed boundary condition, rigorously justified by matched asymptotic expansions.","lead":"This paper rigorously derives the sharp-interface limit of a vector-valued Allen-Cahn equation whose potential vanishes on two concentric spheres. It shows the interface follows mean curvature flow while the order parameter evolves by harmonic map heat flow with a mixed jump condition, extending earlier radial-symmetric results.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The arbitrary-order approximate solution rests on an unproved existence/regularity theorem for the coupled mixed-boundary parabolic system (4.1)-(4.4); without it Theorems 1.1 and 1.3 do not follow.","rationale":"The reader's weakest_assumption identifies exactly the same structural gap: Section 4.2.2 asserts solvability of (4.1)-(4.4) with one sentence and no theorem, and the induction step (4.7) inherits that gap. My stress-test confirms this is the most load-bearing concern. The approximate solution u_K is the object on which both Theorem 1.1 and Theorem 1.3 rest; if the first-order coupled system cannot be solved, the construction stops and no high-order convergence can be claimed. I do not see a more fundamental objection. Concerns about the spectral lower bound are real but are addressed in detail in Section 5, with the scalar reductions and cross-term lemmas constituting a coherent (though compressed) argument; the global orthonormal frame E_i is a topological nuisance that can be localized. The mismatch between K=k+1 and the residual order is not fatal because the rescaled variable u_R only needs to remain O(1), which the O(1) forcing admits. The limit-system well-posedness is an explicit hypothesis, so it is not a flaw in the theorem as stated. Thus the reader's CONDITIONAL verdict is appropriate: the paper is a serious and plausible analysis, but the missing existence/regularity proof for the coupled mixed-boundary parabolic systems is a genuine, non-routine gap and should be supplied before full acceptance. I recommend no change to the reader's verdict; my agreement is full.","tokens_in":43329,"tokens_out":17406,"duration_ms":185264,"concrete_test":"Isolate the first induction step in a flat-interface half-space model: take Γ = R^{m-1}, d_0 = x_m, and choose smooth but nonconstant ω± on Γ with ω+=ω− and b²∂_νω+=a²∂_νω−. Linearize (4.1)-(4.4) around d_1=0, σ±_{1,β}=0, and compute the principal symbol of the resulting linearized boundary-value problem. Verify whether the surface operator in (4.1) is parabolic and whether the transmission boundary condition (4.3) satisfies the Lopatinski condition; then check whether the quadratic term in d_1 can be handled by the standard quasilinear theory (or by a contraction argument on a short time interval). If the symbol computation fails or the quasilinear step is not contractive, the concern lands and the proof of Theorem 1.1 is incomplete. If the model problem is well-posed and all terms in H are lower order, the gap is an omitted standard argument and the central claim is likely repairable.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Theorem 1.1 requires constructing u_K to arbitrary order, and the induction lives or dies on solving the coupled system for (d_1, σ±_{1,β}) in Section 4.2.2, then the analogous system for every k in (4.7). The only justification offered is: 'by trace theorem and basic parabolic estimate of existence, one can then determine d_1|Γ and σ±_{1,β}|Ω±' (Section 4.2.2). This is not a routine step. The d_1 equation (4.1) is a quasilinear parabolic equation on the moving interface: L_{1,0}|Γ depends linearly on d_1 (as the text itself claims), and equation (3.50) contains products like (d_1−z)L_{1,0}η_1ρ′_0, so d_1 enters quadratically. The text even states 'the equation of d_1 is nonlinear' before (3.51). Standard linear parabolic theory does not cover such an equation without smallness, quasilinear structural conditions, or a fixed-point argument, none of which is supplied. The σ±_{1,β} problem (4.2)-(4.4) is a transmission problem with a mixed Robin-type condition (4.3) coupling normal traces on the two sides and containing integrals of the moving-frame coefficients; no Lopatinski or Schauder analysis is given. The same gap recurs in the induction step (4.7). Consequently the existence of u_K—and therefore the error estimate of Theorem 1.3, which is stated for that u_K—is conditional on an unverified analytic input. This is the most load-bearing weakness: if the system (4.1)-(4.4) is not solvable with the required regularity, the expansion cannot be promoted past first order, and the high-order convergence claim collapses. The rest of the paper is a substantial, careful matched-expansion and spectral analysis, but this missing existence step is not a stylistic omission; it is a genuine gap in the proof as written.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the sharp-interface limit, as ε→0, of the vector-valued Allen-Cahn equation (1.5) for a radial potential F vanishing on two spheres |u|=a and |u|=b. The claimed limit is the system (1.8): harmonic map heat flow into S^{n-1} in the two bulk phases, mean curvature flow of the interface, and the mixed boundary condition b²∂_ν ω_+ = a²∂_ν ω_- across the interface. The proof is by matched asymptotic expansions: outer expansions in Ω±, an inner expansion with fast variable z=d/ε, an inductive construction of an arbitrary-order approximate solution u_K (Theorem 1.1), a spectral lower bound for the linearized operator around u_K via a reduction to scalar one-dimensional quadratic forms and cross-term/correction estimates (Theorem 1.2), and finally an energy estimate for u^ε - u_K (Theorem 1.3). The paper also records the formal derivation of the jump condition at (3.42).","tokens_in":43677,"tokens_out":12379,"duration_ms":136853,"significance":"If the main results were fully established, this would be a substantial contribution: it extends the Bronsard--Stoth analysis beyond n=2 and removes the radial symmetry assumption, and it gives a mechanism, based on the explicit construction of η_1(z), for handling unequal radii a≠b in the interfacial condition. The reduction of the vector-valued spectral problem to two scalar one-dimensional quadratic forms is an interesting and potentially reusable idea. However, the current manuscript has a load-bearing gap in the inductive construction of the approximate solutions, and the closing arguments of both the spectral proof and the error estimate are incomplete. The significance of the paper is therefore conditional on resolving these issues.","major_comments":[{"comment":"The induction for Theorem 1.1 hinges on solving the coupled parabolic system for (d_1, σ±_{1,β}) in §4.2.2 and then the analogous system (4.7) at every order. The text says only: “by trace theorem and basic parabolic estimate of existence, one can then determine d_1|Γ and σ±_{1,β}|Ω±.” This is not a routine step. As the authors themselves note, the d_1 equation (4.1) is nonlinear (see (3.50) and the discussion before (3.51)): L_{1,0} depends on d_1 linearly and the term L_{1,0}∫(d_1−z)η_1 ρ^1_0 dz is quadratic in d_1. The σ± equations (4.2)–(4.3) form a transmission problem with nonlocal mixed Robin-type boundary conditions and contain H(u_1,d_1). No existence theorem, fixed-point argument, smallness condition, or Lopatinski/regularity analysis is supplied, and the required regularity of the traces is not stated. Without a proof of solvability of (4.1)–(4.4) and of each step in (4.7), th","section":"§5.6, Eqs. (5.18)–(5.19)"},{"comment":"The proof of Theorem 1.2 ends with the sentence: “It is enough to estimate (5.18) and (5.19). Theorem 1.2 is derived directly.” The inequalities (5.18)–(5.19) are the claimed reductions, but the paper never carries out the final assembly. One must explicitly combine Lemma 5.7 (or its analogue for every j), Lemma 5.8, Lemma 5.9, and Lemmas 5.4–5.5, and choose the small parameter ν_0 strictly smaller than the coefficient 1/8 in the left-hand sides. The boundary-term estimate at the end of §5.6 also needs to be incorporated. As written, the proof verifies the ingredients but not the claimed scalar inequalities, so the spectral lower bound is not closed. This appears fixable by a careful bookkeeping argument, but it must be written out.","section":"§6, Eqs. (6.1), (6.6)"},{"comment":"In the error estimate, after defining u_R=(u^ε−u_K)/ε^k, the remainder term is bounded as R=O(ε^{K−1−k}) and then B^i R=O(ε^{K−1−k−i}) is stated after (6.1). With the stated choice K=k+1 this gives no positive power even for i=0 and negative powers for i≥1. The subsequent estimates (6.2)–(6.5) concern only the nonlinear terms ε^{k−2}B^i \\hat H; the term ∫ B^i R B^i u_R is not estimated and is absent from (6.6). Thus the differential inequality for E(u_R) is not justified. To close the continuation argument, either the approximate solution must be constructed with enough derivative estimates so that B^i R has a positive power of ε for all i≤⌊m/2⌋+1, or K must be chosen much larger than k plus the number of derivatives. This is a load-bearing issue for Theorem 1.3.","section":null}],"minor_comments":[{"comment":"The displayed formula f_B(ρ_0)=ρ_0²/ρ_0 should presumably be ρ_0''/ρ_0. Please correct the notation and verify the subsequent identities.","section":"Lemma 5.7"},{"comment":"Lemma 5.7 is stated only for i=0, j=1, but inequality (5.18) is needed for every j=1,…,n−1. The proof must state that the same estimate holds for all j with E_j replacing E_1, or provide the frame-independent argument.","section":"Theorem 1.1"},{"comment":"Theorem 1.1 only states R=O(ε^{K−1}) without specifying the norm or derivative bounds. Section 6 requires estimates on B^i R, so the construction should state explicitly what C^k/H^s controls are obtained for the remainder.","section":"Introduction"},{"comment":"The abstract and introduction invoke the “quasi-minimal connecting orbits” framework from [12], but the actual construction in §3.2 is explicit, using the heteroclinic ρ_0 and a geodesic interpolation. The relation to [12] should be clarified to avoid the impression that a black box is being used.","section":null}],"recommendation":"major_revision","confidential_remarks":"The central issue is §4.2.2: the inductive construction of u_K requires a real existence/regularity theorem for a coupled, nonlinear, transmission-type parabolic system with mixed boundary conditions, and the current one-sentence justification is far from sufficient. I would require a complete proof of this step before acceptance. The spectral and error-estimate gaps (§5.6 and §6) are less fundamental but also must be repaired. If the missing existence theorem can be supplied, the paper would be a strong contribution to the vector-valued Allen–Cahn asymptotic limit literature."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this is a serious new theorem — the first sharp-interface limit for a vector-valued Allen-Cahn equation with two concentric spherical wells in general dimension and without radial symmetry. It extends Bronsard-Stoth (n=2, radial) and uses the quasi-minimal connecting orbit framework of Fei-Lin-Wang-Zhang. If the proof fills in, it's a major advance in the KRS problem.\n\nWhat the paper does well: the matched asymptotic expansion is detailed and coherent. The derivation of the mixed boundary condition b^2∂ν ω+ = a^2 ∂ν ω- from compatibility conditions is careful. The choice of the interpolation η1(z) seems to effectively handle the unequal coefficients. The spectral reduction to scalar one-dimensional operators in Section 5 is a substantial technical contribution; Lemmas 5.2-5.5 look solid.\n\nNow the soft spots, in order of weight. First, the construction of the approximate solution u_K rests on solving the coupled parabolic system (4.1)-(4.4) for d1 and σ±1,β on the interface and in the bulk. The paper asserts this with 'by trace theorem and basic parabolic estimate of existence' and gives no reference or argument. That is not a routine step: the d1 equation is nonlinear — the text itself notes this before (3.51) — and the RHS contains quadratic dependence on d1. The same gap reappears at every induction step in (4.7). Without an existence/regularity theorem for these systems, Theorem 1.1 is not proved, and Theorem 1.3, which is stated for that approximate solution, does not follow. This is load-bearing.\n\nSecond, the proof of Theorem 1.2 jumps from (5.18)-(5.19) to 'Theorem 1.2 is derived directly.' The final assembly of the cross-term and correction-term estimates is sketched but not written out. I think it's repairable, but a referee will want the details.\n\nThird, Theorem 1.3's continuation argument is compressed into a few lines. Standard, but it should be explicit.\n\nNo circularity or parameter-fitting. The citation to [12] is appropriate.\n\nWho's this for: people working on phase transitions, singular perturbations, and the KRS problem. It deserves a serious referee; I'd send it out with expectation of major revision, asking for the missing solvability proof and a fuller spectral assembly.","headline":"Real new theorem for the KRS problem in general dimension, but the approximate-solution construction has a load-bearing regularity gap.","tokens_in":44260,"tokens_out":4387,"would_cite":false,"duration_ms":47739,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35B25","35K57","35K55"],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper establishes that, as ε→0, the two-sphere vector Allen-Cahn equation converges to mean curvature flow, harmonic map heat flow, and the weighted jump b²∂ν ω+ = a²∂ν ω−.","keywords":["vector-valued Allen-Cahn","sharp interface limit","mean curvature flow","harmonic map heat flow","matched asymptotic expansions","spectral lower bound","Keller-Rubinstein-Sternberg problem","two-well potential"],"falsifier":"Solve the first-order system (4.1)–(4.4) for a smooth non-planar limit solution, e.g., a shrinking sphere in $\\mathbb{R}^3$ with $n=3$; if $d_1|_\\Gamma$ and $\\sigma_{\\pm 1,\\beta}$ are not determined (or require additional compatibility data) by the asserted trace and parabolic-existence step, the arbitrary-order approximate solution of Theorem 1.1 does not exist and the convergence proof collapses.","tokens_in":43158,"feed_emoji":"","tokens_out":9853,"duration_ms":104319,"temperature":0.7,"texified_at":"2026-08-05T19:58:29.687051+00:00","pith_summary":"This paper establishes the sharp-interface limit of the vector-valued Allen-Cahn equation whose potential vanishes on two concentric spheres of radii a and b in $\\mathbb{R}^n$. Its main theorem says that, as the interface thickness $\\varepsilon \\to 0$, the interface itself moves by mean curvature flow, the vector field in the two bulk regions evolves by harmonic map heat flow into $S^{n-1}$, and across the interface the unit directions satisfy $\\omega_+ = \\omega_-$ together with the weighted normal-derivative jump $b^2 \\partial_\\nu \\omega_+ = a^2 \\partial_\\nu \\omega_-$. The proof constructs approximate solutions to arbitrarily high order in $\\varepsilon$ by matched asymptotic expansions, then proves a uniform spectral lower bound for the linearized operator at the approximate solution, which yields the quantitative error estimate $E(u_\\varepsilon(t) - u_K(t)) \\leq C_1 \\varepsilon^{2k}$. This removes the radial-symmetry assumption of the $n=2$ result of Bronsard-Stoth and provides the first rigorous derivation of the full system for general n, identifying the mixed boundary condition as a consequence of the compatibility conditions rather than an imposed ansatz.","texify_model":"deepseek-v4-flash","texify_usage":{"total_tokens":8307,"prompt_tokens":996,"completion_tokens":7311,"prompt_tokens_details":{"cached_tokens":0},"prompt_cache_hit_tokens":0,"prompt_cache_miss_tokens":996,"completion_tokens_details":{"reasoning_tokens":6180}},"feed_headline":"Sharp limit of two-sphere Allen-Cahn: curvature flow plus harmonic maps","feed_subtitle":"Rigorous ε→0 limit couples mean curvature flow, harmonic map heat flow, and weighted jump b²∂ν ω+=a²∂ν ω−.","key_machinery":"The load-bearing object is the inner profile $u_0(z,x,t) = \\rho_0(z)\\bar{\\omega}(\\eta_1(z),x,t)$, where $\\rho_0$ is the heteroclinic solution of $\\rho'' = f(\\rho)$ connecting $a$ to $b$, and $\\bar{\\omega}$ runs along a geodesic of $S^{n-1}$ from $\\omega_-$ to $\\omega_+$. The interpolation function $\\eta_1(z) = \\frac{b^2 - a^2 b^2 F(z)}{b^2 - a^2}$, with $F(z) = z^{-1} \\int_0^z \\rho_0(s)^{-2} \\, ds$, is chosen so that $B_\\nu \\bar{\\omega} = \\left[ \\frac{a^2 b^2 F}{b^2 - a^2} \\right] (B_\\nu \\omega_- - B_\\nu \\omega_+) + \\ldots$ ; combined with the jump condition $b^2 B_\\nu \\omega_+ = a^2 B_\\nu \\omega_-$ it makes $B_r(\\rho_0^2 B_r E_0) = O(1)$ at the interface. This reduces the spectral estimate for the linearized operator to two scalar one-dimensional quadratic forms $Q_0$ and $Q_1$ with potentials $\\rho_0^3 / \\rho_0'$ and $\\rho_0^2 / \\rho_0$, which are then estimated by integration by parts against the heteroclinic profil","core_discovery":"Dong and Wang claim that for any smooth solution ($\\Gamma_t, \\omega_\\pm$) of the sharp-interface system (1.8) on $[0,T]$, and for well-prepared initial data satisfying $E(u_\\varepsilon(0) - u_K(0)) \\leq C_0 \\varepsilon^{2k}$, the solution $u_\\varepsilon$ of (1.5) satisfies $E(u_\\varepsilon(t) - u_K(t)) \\leq C_1 \\varepsilon^{2k}$ for all $t \\in [0,T]$ (Theorem 1.3). The approximate solution $u_K$ is built in Theorem 1.1 by gluing outer expansions $u_{\\pm 0} = b\\omega_+$ in $\\Omega_+$ and $a\\omega_-$ in $\\Omega_-$ with inner expansions in the fast variable $z = d_\\varepsilon/\\varepsilon$; the leading inner state is $u_0 = \\rho_0(z)\\bar{\\omega}(\\eta_1(z),x,t)$, where $\\rho_0$ solves the heteroclinic ODE $\\rho'' = f(\\rho)$ with limits a and b. The key new point is that the first compatibility condition forces the interface speed to be mean curvature and the second forces the mi","pith_inferences":["The same construction should work for potentials vanishing on two concentric spheres with any dimensions and without the explicit quartic form, as long as the heteroclinic profile ρ0 and the geodesic \\bar ω exist; the key identity Br(ρ0² Br E0)=O(1) depends only on the ratio b²/a² entering the jump condition.","A direct testable extension: take the potential with unequal concentric sphere radii a≠b and n=2; the jump condition predicts the tangential component of ∂ν ω at the interface is weighted by b²/a², a quantitative signature that could be observed numerically in a shrinking sphere simulation.","The unproved regularity step in Section 4.2.2 suggests the full theorem may require a smallness or compatibility condition on the limit solution (Γt, ω±) for the coupled parabolic system to be solvable; supplying that existence proof would make the arbitrary-order construction unconditional.","If the spectral reduction to two scalar 1D operators is generic, the stability of the diffuse interface is governed by the heteroclinic profile alone; this might allow the same lower bound to be proved for non-smooth approximate profiles via a limiting argument, opening a route to the Keller-Rubinstein-Sternberg problem for general manifolds M1, M2."],"forward_implications":["For any smooth solution of (1.8), the constructed uK approximates the true Allen-Cahn solution in the energy E to order ε^{2k}; the convergence rate is limited only by how many expansion terms one is willing to compute.","The radial-symmetry restriction of Bronsard-Stoth is removed and the result holds for all n≥2 with the two-sphere potential, so the same sharp-interface system is available for vector order parameters in higher codimension.","The mixed boundary condition b²∂ν ω+ = a²∂ν ω− is derived rather than assumed; it shows that the ratio of the well radii controls how the S^{n-1}-valued field transmits normal derivatives across the interface.","The uniform spectral lower bound, once established at the approximate solution, gives the coercivity needed to run an energy method for the full nonlinear error, so the result transfers to perturbations of well-prepared initial data by standard continuation."],"supporting_citations":[{"why":"Supplies the n=2 radially symmetric case and the two-sphere potential whose limit this paper extends to general n.","marker":"[7]"},{"why":"Supplies the quasi-minimal connecting orbit method and arbitrary-order matched expansion used for the inner construction.","marker":"[12]"},{"why":"Introduces the vector-valued Allen-Cahn system and the formal limit to mean curvature flow plus harmonic map heat flow.","marker":"[31, 32]"},{"why":"Provides the static energy expansion for potentials with two high-dimensional wells that identifies the bulk and interface energy constants.","marker":"[23]"},{"why":"Provides a rigorous sharp-interface limit for a matrix-valued liquid crystal model, the closest prior dynamic result with a non-trivial interface condition.","marker":"[14]"},{"why":"Supplies the inner-outer gluing technique and matched asymptotic expansion framework used throughout the construction.","marker":"[3]"},{"why":"Supplies the spectral analysis of the linearized scalar Allen-Cahn operator, the scalar analogue whose vector version Theorem 1.2 extends.","marker":"[9]"}],"fun_headline_variants":["Two-sphere Allen-Cahn: sharp limit couples curvature flow and harmonic maps","Rigorous sharp limit for two-sphere Allen-Cahn: curvature + harmonic maps","Allen-Cahn sharp limit: mean curvature flow + harmonic heat","Two-sphere Allen-Cahn: sharp limit yields curvature flow and harmonic maps"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The load-bearing premise is that the coupled parabolic system (4.1)–(4.4) for the first-order corrections $d_1$ and $\\sigma_{\\pm 1,\\beta}$ has a solution with enough regularity, asserted in Section 4.2.2 without proof; every higher-order inductive step inherits the same unstated regularity requirement.","fun_headline_variants_meta":{"raw":{"variants":["Two-sphere Allen-Cahn: sharp limit couples curvature flow and harmonic maps","Rigorous sharp limit for two-sphere Allen-Cahn: curvature + harmonic maps","Allen-Cahn sharp limit: mean curvature flow + harmonic heat","Two-sphere Allen-Cahn: sharp limit yields curvature flow and harmonic maps"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001105,"raw_usage":{"total_tokens":4559,"prompt_tokens":973,"completion_tokens":3586,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":717,"completion_tokens_details":{"reasoning_tokens":3504}},"tokens_in":717,"tokens_out":3586,"duration_ms":29694,"temperature":1.0,"reasoning_tokens":3504,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T16:14:20.617244+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Solve the first-order system (4.1)–(4.4) for a smooth non-planar limit solution, e.g., a shrinking sphere in $\\mathbb{R}^3$ with $n=3$; if $d_1|_\\Gamma$ and $\\sigma_{\\pm 1,\\beta}$ are not determined (or require additional compatibility data) by the asserted trace and parabolic-existence step, the arbitrary-order approximate solution of Theorem 1.1 does not exist and the convergence proof collapses.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the n=2 radially symmetric case and the two-sphere potential whose limit this paper extends to general n."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the quasi-minimal connecting orbit method and arbitrary-order matched expansion used for the inner construction."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the static energy expansion for potentials with two high-dimensional wells that identifies the bulk and interface energy constants."},{"cited_title":"Peking Math","cited_arxiv_id":null,"evidence_quote":"Provides a rigorous sharp-interface limit for a matrix-valued liquid crystal model, the closest prior dynamic result with a non-trivial interface condition."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the inner-outer gluing technique and matched asymptotic expansion framework used throughout the construction."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the spectral analysis of the linearized scalar Allen-Cahn operator, the scalar analogue whose vector version Theorem 1.2 extends."}],"review_version":1}