{"id":"c559b7fa-5f2a-4822-8f54-b23a8e09cdc5","arxiv_id":"2608.11586","paper_version":1,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"For p-Laplace type equations, Phi(Du) is C^k whenever Phi is positively homogeneous of degree m with m>k/alpha and Du is C^{0,alpha}.","lead":"This paper proves that well-chosen powers and other homogeneous functions of the gradient of p-Laplace solutions become classically differentiable to any desired order, even where the gradient itself is not smooth. The result gives a clean sufficient condition, degree greater than wanted derivatives divided by the Hölder exponent, and applies to elliptic, parabolic, and porous-medium variants.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No significant objection identified","rationale":"The reader identified the a priori Hölder regularity (1.2) and the analogous parabolic and porous-medium hypotheses as the weakest assumptions. I partially agree: those are the only genuinely conditional inputs, but they are explicit hypotheses of the theorems rather than hidden gaps. The central elliptic theorem is internally consistent: the intrinsic radius r≃|Du|^{1/α} makes the rescaled equation uniformly elliptic, the condition α≤1/(p-1) exactly controls the rescaled right-hand side, the derivative decay estimates are dimensionally consistent, and the flat extension lemma supplies the C^k conclusion with vanishing derivatives on the critical set. I found no algebraic error, no missing factor in the scaling identities (3.14) and (4.2), and no circular step. The parabolic and porous-medium sections are compressed but honestly conditional, and their intrinsic time scales ϑ=r^2 q^{2-p} and ϑ=r^2 q^{1-µ} are correct. Therefore the reader's ACCEPT verdict remains appropriate.","tokens_in":18843,"tokens_out":30911,"duration_ms":336012,"concrete_test":"Check the sharp one-dimensional model u(x)=(p-1)/p |x|^{p/(p-1)} with f=1: the exact C^k threshold for |Du|^m=|x|^{m/(p-1)} is m>k(p-1). Verify that Theorem 1.1 with α=1/(p-1) reproduces this threshold, and check for p=3, k=1,2 that derivatives of |Du|^m indeed vanish on {0} exactly when the strict inequality m>k/α holds, confirming that Lemma 2.1's strict condition µ>k is not off by one.","verdict_should_be":"UNCHANGED","load_bearing_attack":"After careful review of Theorem 1.1 and its proof, I find no internal inconsistency or unsupported step in the central elliptic argument. The intrinsic scaling (3.11), the normalized equation (3.16), the uniform bounds on the rescaled right-hand side (3.17)-(3.19), the chain-rule decay (4.1), and the extension Lemma 2.1 all check out. The genuinely conditional input is the a priori Hölder bound (1.2), with analogous parabolic and porous-medium hypotheses (7.4) and (8.2) explicitly assumed and honestly labeled. For the elliptic p-Laplace equation, (1.2) is supplied by classical regularity theory, so this is not a flaw in the claimed implication. The parabolic and porous-medium sections are genuinely conditional but do not bear on the correctness of Theorem 1.1.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the classical higher regularity of positively homogeneous functions of the gradient of solutions to the inhomogeneous p-Laplace equation and several related problems. Theorem 1.1 states that if Du is locally C^{0,\\alpha} with \\alpha \\le \\min\\{\\alpha_0, 1/(p-1)\\}, and \\Phi is smooth away from zero and positively homogeneous of degree m, then \\Phi(Du) is locally C^k whenever m>k/\\alpha, with all derivatives up to order k vanishing on the critical set Z(u). The proof combines an intrinsic rescaling r\\simeq |Du|^{1/\\alpha}, uniform Schauder estimates on a normalized uniformly elliptic equation, and a flat extension lemma across Z(u). Analogous results are proved for autonomous anisotropic operators (Theorem 5.1), elliptic Uhlenbeck systems (Theorem 6.1), parabolic p-Laplace systems (Theorem 7.1), and powers of nonnegative solutions to the porous medium equation (Proposition 8.1).","tokens_in":18972,"tokens_out":41375,"duration_ms":424958,"significance":"If correct, the main theorem provides a clean, explicit sufficient condition for the C^k regularity of nonlinear gradient quantities across the critical set, complementing existing results on stress-field regularity and second-order estimates. The proof is transparent and self-contained for the scalar elliptic case, with no fitted parameters: the threshold is expressed purely in terms of the available H\\\"older exponent and the homogeneity degree. The conditional hypotheses (1.2), (7.4), and (8.2) are honestly stated, and for the scalar elliptic equation (1.2) is supplied by classical regularity theory. The extensions to vectorial, parabolic, and porous-medium settings enlarge the paper's scope. The main limitation is that the systems and parabolic results inherit the a priori H\\\"older assumption on the gradient or solution; this is clearly disclosed in the text. I found no internal inconsistency or circularity in the central derivation, and the constants are explicit or standard.","major_comments":[],"minor_comments":[{"comment":"The uniform parabolic bootstrap leading to (7.20) is asserted with a bare citation to [51]. I recommend adding two sentences explaining that the linearized coefficients DA(D_yV) have a uniform parabolic C^{0,\\bar\\alpha} norm on Q_2, using (7.4) and the definitions of r and \\vartheta, and that iterative application of the linear parabolic Schauder theory gives (7.20). This will make the argument easier for readers to verify.","section":"Section 7, proof of Theorem 7.1"},{"comment":"The notation (p-2)_+ and (2-p)_+ should be explicitly defined (for instance, as max{p-2,0} and max{2-p,0}) to avoid ambiguity.","section":"Remark 7.3"},{"comment":"The proof of (8.5)-(8.6) is summarized as \"the same coordinatewise extension argument used in Theorem 7.1.\" Since the threshold in (8.6) involves \\Lambda_\\mu = \\max\\{\\gamma_\\mu,\\eta_\\mu\\}, I suggest spelling out the margin calculation, e.g., that m-a\\gamma_\\mu-b\\eta_\\mu>1/\\alpha for the spatial step follows from m>k\\Lambda_\\mu and a+b<k, in parallel with the presentation in Theorem 7.1.","section":"Section 8, Proposition 8.1"},{"comment":"The notation \\partial^P_i A^\\beta_j(DV) is not defined. Please define the derivatives of A with respect to the matrix entries of P so that the pointwise form of the differentiated system is unambiguous.","section":"Section 6, equation (6.16)"},{"comment":"The passage from (3.7) to the energy estimate for \\delta_h v is terse. Adding one sentence clarifying that (3.7) is tested with \\zeta^2\\delta_h v before passing to the limit would improve readability.","section":"Section 3, Lemma 3.1 proof"},{"comment":"The statement that the constants in (3.19) may depend on H is slightly imprecise because the factor 1/8 in (3.15) is independent of H. A short clarification that H enters only through the choice of q_0 and the bound on the normalized right-hand side would be helpful.","section":"Section 3, Proposition 3.2"}],"recommendation":"minor_revision","confidential_remarks":"I concur with the reader's assessment that the central elliptic argument is sound and that no load-bearing mathematical error is present. The parabolic and porous-medium sections are more compressed but their arguments are recoverable from standard linear Schauder theory once the uniform estimates on the normalized coefficients are written out. The AI-assistance disclosure is transparent and does not affect my evaluation. The paper is a solid contribution to the regularity theory of p-Laplace-type equations and is well within the scope of math.AP."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Main take: this paper gives a genuinely useful criterion. If you have a p-Laplace-type solution whose gradient is C^{0,alpha}, then any smooth positively homogeneous function of the gradient gains k classical derivatives as soon as the homogeneity degree exceeds k/alpha, and all derivatives up to order k vanish on the critical set. The intrinsic scaling argument is standard, but the statement is new and clean: it applies uniformly to scalar, anisotropic, vectorial, parabolic, and porous-medium settings, and the critical-set vanishing is a nice extra.\n\nThe scalar elliptic proof is complete and honest. The extension lemma (Lemma 2.1) is elementary and correct; the rescaled equation, the bound on the normalized right-hand side via alpha <= 1/(p-1), and the Schauder bootstrap all check out. I verified the scaling identities and the Faà di Bruno expansion; no hidden circularity. The Hölder hypothesis (1.2) is external, from classical theory, so it is not a flaw.\n\nWhere are the soft spots? The parabolic section is the most compressed. The bootstrap leading to (7.20) is asserted with a reference to Ladyzhenskaya–Solonnikov–Ural'tseva and not fully detailed, and the definition of gamma and eta, while consistent, is a bit opaque at first read. The porous-medium section is even shorter; the uniform parabolicity after normalization is plausible, but I would want the details spelled out for the mixed derivative bootstrap. These are matters of exposition, not apparent errors. The vectorial section is also a bit terse in the iteration, but the structure is clear.\n\nThe citation pattern looks fine; the related work is discussed, and the authors are explicit that they assume rather than prove the appropriate gradient Hölder bounds in the anisotropic, vectorial, and parabolic cases. The sharpness remarks are honest, including the constant-source example and the free-boundary obstruction for the porous medium.\n\nBottom line: this is a solid paper. The central theorem is worth knowing, the proof is reproducible in the scalar case, and the conditional extensions are clearly labeled. For a PDE regularity reader it deserves a serious referee. I would suggest a minor revision asking for a more detailed parabolic bootstrap and a short discussion of why the thresholds are natural. I would send it to a good journal.","headline":"A clean, correct sufficient criterion for higher differentiability of homogeneous gradient compositions in p-Laplace-type problems; the scalar elliptic part is solid, the parabolic and porous-medium extensions are more compressed but honest.","tokens_in":19470,"tokens_out":1762,"would_cite":true,"duration_ms":19417,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35J92","35K92","35B65","35J60","35K65"],"pacs":[],"model":"deepseek-v4-flash","headline":"For p-Laplace solutions, any smooth homogeneous function of the gradient of degree m > k/α is classically k-times differentiable, even across the critical set where the gradient vanishes.","keywords":["p-Laplace equation","gradient regularity","homogeneous composition","critical set","Schauder estimate","p-Laplace system","porous medium equation","parabolic p-Laplace system"],"falsifier":"Find a weak solution of $\\operatorname{div}(|Du|^{p-2}Du)=f$ with smooth $f$, whose gradient is Hölder continuous of exponent $\\alpha\\le 1/(p-1)$, and a smooth positively homogeneous $\\Phi$ of degree $m>k/\\alpha$ for which the continuous extension of $\\Phi(Du)$ fails to be $C^k$ at a critical point. The theorem asserts no such pair exists; the one-dimensional profile (1.12) shows the strict inequality is necessary, so any violation would settle the claim.","tokens_in":18664,"feed_emoji":"📐","tokens_out":12525,"duration_ms":127168,"temperature":0.7,"pith_summary":"This paper establishes a higher-regularity principle for the inhomogeneous $p$-Laplace equation: although a solution's gradient need not be $C^2$ near points where it vanishes, sufficiently high powers (or any smooth positively homogeneous function) of the gradient are classically $C^k$ across the entire critical set. The sufficient threshold is $m > k/\\alpha$, where $\\alpha$ is the gradient's Hölder exponent, capped at $\\alpha \\le 1/(p-1)$; the result therefore applies to every $p \\in (1,\\infty)$ once gradient Hölder regularity is known. The proof exploits an intrinsic scale $r \\simeq |Du|^{1/\\alpha}$ on which the equation becomes uniformly elliptic, so Schauder estimates control all derivatives, and an elementary extension lemma carries the zero values across the critical set. The same mechanism yields analogous statements for anisotropic equations, elliptic and parabolic $p$-Laplace systems, and high powers of nonnegative solutions of the porous medium equation.","feed_headline":"Gradient powers gain C^k regularity once degree exceeds k/alpha","feed_subtitle":"Even where the gradient is not C^2, its powers are classically smooth once m > k/alpha.","key_machinery":"The central object is the intrinsic noncritical scale $r\\simeq |Du|^{1/\\alpha}$: at a point where $|Du|=q>0$, this is the radius on which the gradient stays comparable to $q$ by the Hölder assumption. After the normalization $v(y)=(u(x+ry)-u(x))/(rq)$, the equation becomes uniformly elliptic on a fixed ball, so classical Schauder estimates bound all derivatives of $v$ uniformly; the rescaled right-hand side is controlled exactly because $\\alpha\\le 1/(p-1)$. Rescaling back yields derivative decay $|D^\\ell u(x)|\\lesssim q^{1-(\\ell-1)/\\alpha}$ and, via homogeneity and the multivariate Faà di Bruno chain rule, $|D^j(\\Phi(Du))|\\lesssim q^{m-j/\\alpha}$. The flat extension lemma (Lemma 2.1) then converts the decay estimate $|D^j F|\\le C\\,\\operatorname{dist}(x,Z)^{\\mu-j}$ with $\\mu=\\alpha m>k$ into genuine $C^k$ differentiability across the critical set, with all derivatives equal to zero on $Z$.","core_discovery":"On the paper's own terms: Let $u$ solve $\\operatorname{div}(|Du|^{p-2}Du)=f$ with smooth $f$, and assume the gradient is locally Hölder continuous with exponent $\\alpha \\le \\min\\{\\alpha_0, 1/(p-1)\\}$, where $\\alpha_0$ is the classical exponent from gradient regularity theory. Then for any smooth $\\Phi:\\mathbb{R}^d\\setminus\\{0\\}\\to\\mathbb{R}^N$ positively homogeneous of degree $m$, extended by $\\Phi(0)=0$, the composition $\\Phi(Du)$ lies in $C^k_{\\rm loc}(\\Omega;\\mathbb{R}^N)$ whenever $m>k/\\alpha$, and all derivatives of order at most $k$ vanish on the critical set $Z(u)=\\{Du=0\\}$. Quantitatively, near noncritical points the estimates read $|D^j(\\Phi(Du))|\\le C\\,|Du|^{m-j/\\alpha}$ for $0\\le j\\le k$. The same statement holds for autonomous anisotropic operators and for elliptic and parabolic $p$-Laplace systems; in the parabolic case spatial and temporal derivatives carry different homogeneity costs, and the argument gives $C^k$ criteria for high powers $u^m$ of nonnegative solutions of the porous medium equation.","pith_inferences":["If better control of the forcing term (for instance, vanishing or smallness near the critical set) were available, the cap $\\alpha\\le 1/(p-1)$ might be relaxed and the homogeneity threshold lowered; the paper's normalization shows the cap is used only to keep the rescaled right-hand side bounded.","The argument suggests a general 'hidden smoothness' principle for degenerate operators whose ellipticity is homogeneous of degree $p-2$: whenever the gradient is Hölder continuous, homogeneous powers of the gradient become classically differentiable with a threshold tied to the ratio $m\\alpha$; this could be probed for double-phase or variable-exponent operators as soon as gradient Hölder estimate","Because all derivatives up to order $k$ vanish on the critical set, $\\Phi(Du)$ is flat to order $k$ near critical points; this flatness could serve as a quantitative tool to locate the critical set through the leading-order term $|Du|^{m-k/\\alpha}$ in the derivative estimates, for instance in symmetry or nodal-set arguments."],"forward_implications":["For every $1<p<\\infty$, once the gradient Hölder regularity (1.2) is in hand, $|Du|^m\\in C^k_{\\rm loc}$ for all $m>k/\\alpha$; in particular, $|Du|^m\\in C^1_{\\rm loc}$ whenever $m>\\max\\{1/\\alpha_0, p-1\\}$.","The sufficient degree depends on the gradient Hölder exponent and the desired order, but not on the size of the datum $f$: a larger force changes only the constants and the neighborhood on which the intrinsic estimate applies, not the homogeneity threshold.","For vectorial $p$-Laplace systems, the same conclusion holds for $\\Phi(DU)$ provided the full matrix gradient is Hölder continuous; the linearized operator is uniformly strongly elliptic on matrix annuli, which is all the Schauder bootstrap needs.","In the parabolic system, the regularity splits: spatial differentiability up to order $k$ requires $m>k\\gamma$, while full space-time $C^k$ requires $m>k\\eta$, where $\\gamma$ and $\\eta$ encode the separate spatial and temporal homogeneity costs.","For the porous medium equation, $u^m$ gains spatial $C^k$ and space-time $C^k$ regularity once $m$ exceeds the corresponding thresholds $\\gamma_\\mu$ or $\\Lambda_\\mu$, consistent with the free-boundary obstruction that the pressure power $u^{\\mu-1}$ is typically Lipschitz but not $C^1$."],"supporting_citations":[{"why":"Supplies the classical interior Schauder estimates for elliptic equations used to bootstrap regularity of the normalized solution on the intrinsic scale.","marker":"[39]"},{"why":"One of the foundational proofs of local C^{1,\\alpha_0} regularity for p-Laplace solutions, the hypothesis (1.2) on which Theorem 1.1 rests.","marker":"[37]"},{"why":"Same role as [37]: establishes the C^{1+\\alpha} local regularity of weak solutions of degenerate elliptic equations that gives the starting Hölder exponent.","marker":"[30]"},{"why":"Provides the interior Schauder estimates for elliptic systems used in the vectorial p-Laplace theorem.","marker":"[38]"},{"why":"Supplies the parabolic Schauder estimates used to bootstrap regularity in the parabolic system theorem.","marker":"[51]"},{"why":"Extends the gradient Hölder regularity to a broader class of quasilinear elliptic equations, supporting the assumed regularity in the anisotropic theorem.","marker":"[70]"},{"why":"Establishes interior Hölder continuity of the spatial gradient for parabolic p-Laplace systems, the input assumption (7.4) for the parabolic result.","marker":"[31]"}],"fun_headline_variants":["Gradient compositions gain C^k when m exceeds k/alpha","Homogeneous gradient powers are C^k past m > k/alpha","When degree beats ratio, gradient compositions turn C^k","Classical smoothness emerges for gradient powers above threshold","C^k gradient compositions from Holder gradients"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument presupposes that the gradient is already Hölder continuous with some exponent $\\alpha \\le \\min\\{\\alpha_0, 1/(p-1)\\}$; it does not prove that regularity, and the required degree $m>k/\\alpha$ becomes more demanding if only a smaller exponent is available.","fun_headline_variants_meta":{"raw":{"variants":["Gradient compositions gain C^k when m exceeds k/alpha","Homogeneous gradient powers are C^k past m > k/alpha","When degree beats ratio, gradient compositions turn C^k","Classical smoothness emerges for gradient powers above threshold","C^k gradient compositions from Holder gradients"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000876,"raw_usage":{"total_tokens":3837,"prompt_tokens":1040,"completion_tokens":2797,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":656,"completion_tokens_details":{"reasoning_tokens":2719}},"tokens_in":656,"tokens_out":2797,"duration_ms":23654,"temperature":1.0,"reasoning_tokens":2719,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T00:37:17.050017+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Find a weak solution of $\\operatorname{div}(|Du|^{p-2}Du)=f$ with smooth $f$, whose gradient is Hölder continuous of exponent $\\alpha\\le 1/(p-1)$, and a smooth positively homogeneous $\\Phi$ of degree $m>k/\\alpha$ for which the continuous extension of $\\Phi(Du)$ fails to be $C^k$ at a critical point. The theorem asserts no such pair exists; the one-dimensional profile (1.12) shows the strict inequality is necessary, so any violation would settle the claim.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"One of the foundational proofs of local C^{1,\\alpha_0} regularity for p-Laplace solutions, the hypothesis (1.2) on which Theorem 1.1 rests."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the parabolic Schauder estimates used to bootstrap regularity in the parabolic system theorem."},{"cited_title":"DiBenedetto and A","cited_arxiv_id":null,"evidence_quote":"Establishes interior Hölder continuity of the spatial gradient for parabolic p-Laplace systems, the input assumption (7.4) for the parabolic result."}],"review_version":1}