{"id":"a5a0069b-83e6-431b-a634-dafa4968dd3d","arxiv_id":"2607.27119","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.5,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"Unified Transform Method yields novel contour-integral solution formulas for nonhomogeneous IBVPs and interface problems for the Barenblatt–Sobolev–Galpern pseudoparabolic equation.","lead":"The paper derives explicit contour-integral formulas for solutions of the Barenblatt pseudoparabolic PDE on the half-line, finite interval, and across interfaces. These closed-form representations give analysts a direct handle on regularity, asymptotics, and free-boundary extensions in models used for porous media, heat conduction, and solid-state batteries.","discovery_kind":"extension","skeptic_critique":null,"referee_report":{"model":"moonshotai/kimi-k3","summary":"The manuscript applies the Fokas unified transform method (UTM) to the nonhomogeneous Barenblatt-Zheltov-Kochina pseudoparabolic equation u_t - αu_xxt - βu_xx = f. Four problems are treated: (1) a half-line IBVP with the non-standard boundary condition (1.7) mixing u(0,t), u_t(0,t), u_x(0,t), u_xt(0,t); (2) a two-medium interface problem on the line with two interface conditions coupling the half-line solutions, which is shown to reduce to two decoupled half-line problems under the algebraic restriction det Γ12 = det Γ34 = 0 in (1.13); (3) a general interface problem with six interface conditions, solved via an 8×8 Cramer system subject to a list of vanishing 6×6 minors of the coefficient matrix; (4) a finite-interval IBVP with Robin-type analogues (1.15). For Problem 1 the authors derive formula (1.18) from the global relation (2.1) and verify in detail (Theorems 1-3, Section 3) that the contour-integral formula is jointly C∞, solves the PDE, and attains the initial and boundary data uniformly. Problems 2-4 are treated by derivation only, via contour deformations around the pole i/α of the dispersion function ω(λ) = λ²/(1+αλ²) and biholomorphic changes of spectral variable.","tokens_in":47282,"tokens_out":4427,"duration_ms":89281,"significance":"If the results hold, the paper provides the first explicit contour-integral solution formulas for interface problems for the Barenblatt-Sobolev-Galpern pseudoparabolic equation, an equation with a long modeling pedigree (Milne 1926, Barenblatt-Zheltov-Kochina 1960, two-temperature heat conduction) and current applied interest. The half-line result (Problem 1) is not merely derived but verified: Theorems 1-3 establish joint C^∞ regularity, uniform attainment of the initial condition, uniform attainment of the non-standard boundary condition, and spatial decay, by a careful six-step argument involving contour deformation, Jordan's lemma, and the symbol identity (1.16). The derivations are self-contained and parameter-free — no fitted quantities enter. The identification of the non-standard boundary/interface conditions dictated by the global relation, and of the separable class (1.13), are useful structural observations. The finite-interval formula (6.16) is a nontrivial extension of UTM interval technology to a third-order mixed-derivative symbol. These are solid, citable contributions to the UTM literature, conditional on the verification gaps for Problems 3 and 4 being closed.","major_comments":[{"comment":"The manuscript states that 'if this problem is solvable, we have to assume (1.13)' (det Γ12 = det Γ34 = 0). What is actually shown is narrower: if (1.13) fails, the parts of the integrals (4.13)-(4.14) containing the opposite half-line's unknown spectral functions J*_L, J*_R do not cancel, so the UTM elimination step fails. Unsolvability of the interface problem itself does not follow — other methods (e.g., direct spectral or energy methods for the coupled system) are not ruled out. Please rephrase as a sufficient condition for the present method. More importantly for scope: standard physical interface conditions (continuity of u, continuity of flux v = -α∂x u + β∂x∂t u) generically give det Γ12 ≠ 0. The paper should state explicitly which physically motivated interface laws, if any, satisfy (1.13); otherwise the class of interface problems solved here is considerably narrower than the引言","section":"§4, display (1.13)"},{"comment":"The Proposition lists twelve vanishing conditions on 6×6 minors Γ_{j1...j6} (together with the α_Rβ_R, α_Lβ_L combinations) under which (5.7)-(5.8) reduce to the effective formulas (5.9)-(5.10). No example is given of an interface matrix γ_ij satisfying these conditions beyond the separable cases already covered by §4/Problem 2. If the only admissible matrices are (equivalent to) the separable ones, §5 adds no genuinely new solvable class. Please either exhibit a concrete non-separable 6×6 interface system satisfying all twelve conditions — ideally one with a physical interpretation — or show the conditions define a strictly larger algebraic set than the separable locus.","section":"§5, Proposition and displays (5.9)-(5.10)"},{"comment":"Sections 2-3 give a careful derivation and a full six-step verification (Theorems 1-3) for Problem 1, but no analogous verification is provided for the interface and finite-interval formulas (5.9), (5.10) and (6.16). Problem 2 inherits the verification via reduction to Problem 1, but Problems 3 and 4 do not. In §6 the situation is genuinely more delicate than in §3: after substituting (6.10)-(6.11) into (6.3), the integrands acquire denominators Δ(λ) whose zeros (roots of (κ0² - κl²) sin(lλ)-type expressions) must be controlled relative to the contours γ(±i/α), and the boundary-attainment argument of Steps 4-6 needs to be redone for the Robin-type conditions (1.15). A theorem for (6.16) (PDE, initial condition, boundary conditions, with uniformity) and at least a statement for (5.9)-(5.10), with an indication of how the §3 machinery adapts, seems necessary for the title claim about ' ' ","section":"§§5-6, formulas (5.9), (5.10), (6.16)"},{"comment":"The vanishing of the J* and U0, Ul terms in (6.12)-(6.15) is justified only by 'the integrands are bounded in the vicinities of the points ±i/α'. Boundedness alone does not make a contour integral vanish; presumably the intended argument is that γ(±i/α) can be shrunk to the point, provided the integrand is analytic in a punctured neighborhood with at most integrable singularity — but the denominators Δ(λ) and κ0, κl have their own zeros/singularities, and the condition Δ(λ) ≠ 0 is only imposed at generic λ. Please make the argument precise, including the local behavior of Δ near ±i/α (the stated asymptotics Δ(λ) ~ (β/α²)(1+λ²/α²)^{-2} sin(lλ)·(2i)^{-1} as λ→±i/α suggest a blow-up, not boundedness — the two estimates given seem to need reconciliation).","section":"§6, displays (6.12)-(6.15)"}],"minor_comments":[{"comment":"Cross-referencing errors in the theorem statements: Theorem 1 cites assumption '(1.4)' where (1.8) is meant, and refers to the solution formula as '(1.14)' where (1.18) is meant (also in Theorems 2 and 3). After the Notation paragraph, 'the boundary condition (1.3)' should be (1.7). On p. 5, 'assumption (1.9) imply' should be (1.13).","section":"§1, Theorems 1-3"},{"comment":"Reference [38] appears twice with different content (Chatziafratis-Kamvissis-Stratis, Stud. Appl. Math. 2023 and Chatziafratis-Karali-Synolakis, Stud. Appl. Math. 2026); renumber the list.","section":"References"},{"comment":"No uniqueness statement accompanies Theorem 1: the formula is shown to produce a solution, but it is not stated within which class the solution of Problem 1 is unique. Even a brief remark (or a reference to [17]) would help, since the abstract advertises future well-posedness work.","section":"§3"},{"comment":"In (3.5) and (3.9), the interchange of ∂^{m+n}/∂x^m∂t^n with the integrals is stated to hold 'provided M is sufficiently large'; please state the explicit threshold (e.g., M > n + m + 1 in terms of the O(λ^{-(M+1)}) decay in (3.2), (3.6)) so the reader can verify the dominated-convergence step.","section":"§3, displays (3.5), (3.9)"},{"comment":"Assorted typos: 'eqution was obtaind' (p. 2), 'sourse' (§1, Problem 4 discussion), 'orientaion' (Fig. 1 caption), 'procced', 'deterninants' (§4). The definition of σ_M below (3.2) should be displayed as a sum over k of u^{(k-1)}(0)/(iλ)^k — as typeset it is hard to parse.","section":"throughout"},{"comment":"The contour C(i,ε) in (4.3), (4.6) and §5 is used without a definition in the text; only Figures 4-5 hint at it. Please define it explicitly (a small circle around μ = i in the μ-plane, with ε chosen via the biholomorphic maps φ_R, φ_L) and state the orientation convention.","section":"§4, display (4.3)"}],"recommendation":"major_revision","confidential_remarks":"The bibliography is heavily weighted toward the first author's recent UTM output (refs. 15-20, 22, 32, 34-41). These supply technique background rather than the results themselves, so this is not a circularity concern, but the editor may wish to confirm that the present manuscript is adequately delineated from [17] (half-line Barenblatt via UTM) and [22] (general interface problems for advection-diffusion-reaction), so that the novel content here is precisely the mixed-derivative Barenblatt symbol, the non-standard boundary/interface conditions it forces, and the finite-interval formula. The abstract's breadth-of-applications paragraph (batteries, nanotechnology) is disproportionate to what is analyzed; a modest trim would improve the paper's reception."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The one thing worth knowing is that this paper delivers explicit, verified contour-integral representations for the Barenblatt pseudoparabolic equation under fully nonhomogeneous data and nonstandard mixed boundary/interface operators that the global relation itself forces. That is new relative to the authors’ own half-line work in [17] and to the broader UTM literature on polynomial-dispersion equations.\n\nWhat they do well is the classical three-step Fokas program: global relation, algebraic elimination of unknown boundary transforms, then direct verification that the resulting formula (1.18) is C∞, solves the PDE via the symbol identity, recovers the initial datum by Fourier inversion plus Jordan, and attains the mixed boundary condition uniformly on compact time intervals (Theorems 1–3). The finite-interval Robin-type formula (6.16) and the general six-condition interface formulas (5.9)–(5.10) are written out cleanly. Contour deformations and growth estimates are handled carefully; the arguments are long but conventional and self-contained.\n\nThe soft spot is real but limited: for the two-sided interface problem to separate, the 2×2 blocks of interface coefficients must have vanishing determinants. That condition is imposed so the opposite-side spectral transforms cancel inside the deformed contours; it is not derived from the underlying continuum mechanics or battery models the introduction invokes. The applied paragraphs on solid-state batteries and free boundaries are forward-looking only—no analysis is supplied. Low-regularity well-posedness is left open, which is normal for a pure representation paper.\n\nThis is for people who already work with the Fokas method or Sobolev-type equations and need closed-form starting points for asymptotics or nonlinear extensions. Citation pattern is heavy on the group’s prior UTM papers, but those supply technique, not the target formulas. I would send it to referees; the derivations are honest and the formulas are usable. Engage if you need the representations; skip if you only care about the physical interface conditions without the algebraic restriction.","headline":"Solid, incremental UTM formulas for nonstandard Barenblatt IBVPs and interfaces; the math checks out, the physics framing is aspirational, and the det=0 restriction is the main caveat.","tokens_in":40808,"tokens_out":521,"would_cite":false,"duration_ms":19498,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35K70","35A22","35C15","35G16","44A15"],"pacs":[],"model":"grok-4.5","headline":"Contour-integral formulas solve nonstandard interface and boundary problems for the Barenblatt pseudo-parabolic equation.","keywords":["Barenblatt pseudo-parabolic equation","Sobolev-Galpern equations","unified transform method","interface problems","contour integrals","initial-boundary value problems","porous media","solid-state batteries"],"falsifier":"Direct numerical evaluation of the contour integrals against a high-resolution finite-difference or spectral solution of the same initial-boundary-value problem on a half-line or finite interval; systematic discrepancy in the recovered boundary traces or in the PDE residual would refute the claimed formulas.","tokens_in":41073,"feed_emoji":"∫","tokens_out":905,"duration_ms":18971,"temperature":0.7,"pith_summary":"The paper derives explicit solution formulas, written as contour integrals in the complex Fourier plane, for the Barenblatt–Zheltov–Kochina pseudo-parabolic equation on the half-line, the whole line with an interface, and a finite interval. The equation models seepage in fissured rock, non-Newtonian flow, two-temperature heat conduction and transport in batteries; the boundary and interface conditions it forces are nonstandard combinations of the unknown and its mixed derivatives. Using the unified transform method, the authors eliminate unknown boundary transforms via global relations and obtain closed-form integral representations that recover the initial data and the prescribed boundary/interface data. These formulas open the way to asymptotic analysis, regularity theory and free-boundary problems in which an interface moves according to energy balance.","feed_headline":"Contour integrals solve Barenblatt interface problems","feed_subtitle":"Explicit formulas recover nonstandard boundary data for a key pseudo-parabolic model in porous media and batteries","key_machinery":"The unified transform method: a global relation that couples the Fourier transforms of the solution and of the unknown boundary values is combined with algebraic invariances and contour deformation around the poles of the symbol, allowing those unknown transforms to be eliminated and yielding an explicit contour-integral representation.","core_discovery":"Under standard smoothness assumptions on the data, the contour-integral formula (1.18) defines a jointly smooth solution of the nonhomogeneous Barenblatt equation on the half-line that attains the initial condition uniformly and the natural nonstandard boundary condition uniformly on compact time intervals; analogous integral representations are obtained for separable interface problems on the line and for Robin-type problems on a finite interval.","pith_inferences":["Because the admissible interface conditions are constrained by algebraic determinants rather than by physics alone, many physically natural jump relations may lie outside the solvable class and will require a different spectral treatment.","The same global-relation machinery should apply, with only notational changes, to systems of coupled Barenblatt-type equations that appear in multi-temperature or multi-porosity models.","Once the linear formulas are in hand, a fixed-point argument in appropriate function spaces becomes feasible for mild nonlinearities such as power-law source terms or concentration-dependent mobility."],"forward_implications":["The explicit formulas supply the linear building blocks needed for asymptotic and well-posedness studies of the same equation.","They furnish exact benchmarks for numerical schemes applied to pseudo-parabolic models arising in battery and porous-media transport.","The same contour technique extends, in principle, to free-boundary problems in which the interface location itself evolves by an energy-balance law.","The half-line and finite-interval representations reduce, when the higher-order coefficient vanishes, to the classical Robin problem for the heat equation, recovering known formulas as a special case."],"fun_headline_variants":["Contour integrals solve Barenblatt interface problems","Explicit contour formulas for Barenblatt pseudo-parabolic interfaces","Integral reps recover nonstandard Barenblatt boundary data","Unified transform yields contour solutions on line and interval","Contour integrals handle nonhomogeneous Barenblatt interface PDEs"],"cache_read_input_tokens":32896,"weakest_assumption_plain":"For the two-sided interface problem to separate into independent half-line problems, two 2-by-2 blocks of interface coefficients must have vanishing determinants; otherwise the unknown spectral data from the opposite side do not cancel and the elimination step fails.","fun_headline_variants_meta":{"raw":{"variants":["Contour integrals solve Barenblatt interface problems","Explicit contour formulas for Barenblatt pseudo-parabolic interfaces","Integral reps recover nonstandard Barenblatt boundary data","Unified transform yields contour solutions on line and interval","Contour integrals handle nonhomogeneous Barenblatt interface PDEs"]},"model":"grok-4.5","effort":"low","cost_usd":0.005221,"raw_usage":{"total_tokens":1474,"prompt_tokens":797,"num_sources_used":0,"completion_tokens":63,"cost_in_usd_ticks":52208000,"prompt_tokens_details":{"text_tokens":797,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":614,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":797,"tokens_out":63,"duration_ms":11663,"temperature":1.0,"reasoning_tokens":614,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-30T12:13:46.826904+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"Direct numerical evaluation of the contour integrals against a high-resolution finite-difference or spectral solution of the same initial-boundary-value problem on a half-line or finite interval; systematic discrepancy in the recovered boundary traces or in the PDE residual would refute the claimed formulas.","supporting_citations":[],"review_version":1}