{"id":"6f08f390-3cd0-4443-b31c-f6109f422fc0","arxiv_id":"1908.06873","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"For several classes of cross-diffusion systems, normal ellipticity of the diffusion matrix together with symmetry of the transformed matrix is shown to be sufficient, and for constant or nearly constant matrices equivalent, to the existence of an entropy structure.","lead":"This paper gives conditions under which cross-diffusion systems with nonsymmetric diffusion matrices can be rewritten in a form where the transformed diffusion matrix is positive definite, the entropy structure needed for global existence proofs. It proves normal ellipticity for several concrete models and constructs new entropy densities, including a quadratic entropy for fluid mixtures.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The sufficiency result in Theorem 7(ii) rests on symmetry of h''A (detailed balance); without it, the gap between normal ellipticity and entropy structure remains open.","rationale":"The reader's conditional verdict is appropriate. The matrix-factorization arguments underlying Theorem 7 are sound: the necessity of normal ellipticity, the sufficiency under symmetry, and the diagonalizability consequence all follow from Propositions 2, 3, and 6, which are essentially Sylvester inertia and Lyapunov arguments. The local error in Remark 12 (the determinant identity) is real but does not enter the proofs of the main theorems. The most load-bearing issue for the paper's central claim is the reliance of the general sufficiency criterion on the symmetry of h''A, which in practice means the detailed-balance condition. This is not a flaw in the proof but a genuine gap in the characterization: the paper does not resolve whether normal ellipticity alone, or some weaker condition, can suffice for an entropy structure with nonsymmetric Onsager matrix. The paper itself acknowledges this open question, and the reader's weakest_assumption identifies exactly this point. No adjustment to the conditional verdict is needed; the paper should be accepted or conditionally accepted with the requested clarifications and correction of Remark 12.","tokens_in":20229,"tokens_out":39815,"duration_ms":365528,"concrete_test":"Use the open 3x3 example after Lemma 21, A(u) = [[u3,0,u1],[u2,u1,0],[0,u3,u2]] for u > 0, which is normally elliptic and does not satisfy detailed balance. Attempt to construct a strictly convex h explicitly, for example by solving the Hessian integrability equations with a polynomial ansatz of low degree and checking positive definiteness of h''(u)A(u) symbolically or via semidefinite programming over a grid of u. If such an h exists, the symmetry condition in Theorem 7(ii) is not necessary and the gap is real; if a proof of nonexistence emerges, the detailed-balance route may be unavoidable.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Theorem 7(ii) is the paper's main general tool for constructing entropies: if A(u) is normally elliptic and there exists a strictly convex h with h''(u)A(u) symmetric, then h''(u)A(u) is positive definite. In every application (SKT, volume-filling, fluid mixtures) this symmetry is supplied by the detailed-balance condition (10), which is a strong structural assumption on the coefficients and is not derived from the equations. For a normally elliptic diffusion matrix that fails such a symmetry, the paper provides no sufficient condition and explicitly leaves open whether an entropy structure with a nonsymmetric Onsager matrix can exist (see the example after Lemma 21). Thus the central characterization is conditional rather than complete: the logical gap between normal ellipticity and existence of an entropy is not closed. This is a limitation of scope, not an internal inconsistency, but it is load-bearing because Theorem 7(ii) is the statement that would justify applying the entropy method to a new model.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies when the diffusion matrix A(u) of a quasilinear cross-diffusion system (1) admits an entropy structure, i.e., a strictly convex h with h''(u)A(u) positive definite. The main results are: an entropy structure implies normal ellipticity of A(u) (Theorem 7(i)); if A(u) is normally elliptic and h''(u)A(u) is symmetric, then h''(u)A(u) is positive definite (Theorem 7(ii)); for constant or nearly constant diffusion matrices, normal ellipticity is equivalent to existence of an entropy structure (Propositions 11 and 13); and for entropies that are sums of single-variable functions, positive definiteness reduces to positivity of all leading principal minors of A(u) (Proposition 14). The paper also constructs entropies for SKT, volume-filling, and fluid-mixture models, and proves normal ellipticity of the SKT diffusion matrix without the detailed-balance condition (Lemma 20), together with related results.","tokens_in":20404,"tokens_out":27281,"duration_ms":250787,"significance":"The matrix-factorization approach is attractive and, for the most part, cleanly executed. Theorems 7, Propositions 2, 3, 6, 11, 13, 14, and Lemmas 20-22 appear correct and provide useful, checkable criteria. The new normal-ellipticity result for the n-species SKT model without detailed balance is a genuine contribution, and the paper is honest about the remaining open gap between normal ellipticity and entropy structure (see the example after Lemma 21). However, the 'second entropy' results in Section 6 are not justified as stated: the proofs establish positivity of h''(u)A(u) but do not verify that the constructed h is strictly convex, and the stated assumptions do not imply convexity.","major_comments":[{"comment":"The proofs show only that z^T h''(u)A(u)z > 0 for z != 0, but the definition of entropy structure also requires h to be strictly convex. The assumptions in Proposition 17 do not imply that h'' = diag(pi) Q is positive definite. For example, take n = 2, pi = (1,1), p1 = u1 + 2u2, p2 = 2u1 + u2. The detailed-balance condition (10) holds, Q = [[1,2],[2,1]] is invertible, and h''(u)A(u) = Q diag(u) Q is positive definite for u in R_+^2. However, the Poincare potential h = (1/2)u1^2 + 2u1u2 + (1/2)u2^2 has Hessian Q with eigenvalues 3 and -1, so h is not convex on any convex open domain. Proposition 19 has the same defect. The propositions need an additional assumption that ensures h'' is positive definite, such as normal ellipticity of Q as in Proposition 15, or a modified conclusion.","section":"Section 6, Propositions 17 and 19"},{"comment":"The displayed identity det H = integral over [0,infinity) of det(e^{-A^T t}) det(e^{-A t}) dt is false, because the determinant does not commute with integration. For A = diag(lambda1, lambda2), the Lyapunov solution is H = diag(1/(2 lambda1), 1/(2 lambda2)), so det H = 1/(4 lambda1 lambda2), whereas the formula in the remark gives 1/(2(lambda1 + lambda2)). The formula H = integral_0^infty e^{-A^T t} e^{-A t} dt is correct, but the determinant consequence should be removed or corrected.","section":"Remark 12"}],"minor_comments":[{"comment":"For the system as written, the diffusion matrix is A = [[1, -u1], [1, delta]], but the displayed factorization A1 A2 equals [[1, -u1], [delta, 1]]. The reported eigenvalues 1 ± i sqrt(delta u1) correspond to the factorized matrix, not to the stated system. Please correct either the system or the factorization and entropy computation.","section":"Section 3, Case 1.1, Eq. (7)"},{"comment":"The example matrix [[u3,0,u1],[u2,u1,0],[0,u3,u2]] satisfies the conditions of Lemma 21 with the triple (3,1,2), not with (1,2,3), since b11 = b22 = b33 = 0 in this example.","section":"Section 7, after Lemma 21"},{"comment":"There are several typographical errors: 'elipticity' in the abstract should be 'ellipticity'; 'symmetrix' in Case 2.1 should be 'symmetric'; 'tripel' in Lemma 21 should be 'triple'.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The main body of the paper (Sections 2-5 and 7) is largely sound, but the Section 6 'second entropy' propositions assert false statements as written. If the authors add a suitable convexity/normal-ellipticity assumption or weaken the conclusion appropriately, the revision should be straightforward. I would also ask them to fix the false determinant identity in Remark 12 and the inconsistent Keller-Segel example before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a solid, useful paper. The two results I didn't know were in the literature are the normal ellipticity of the n-species SKT diffusion matrix without detailed balance (Lemma 20) and the existence of a second, quadratic entropy for fluid mixtures under a curl-free condition (Proposition 17). The matrix-factorization framework around Theorem 7 is clean, and the central proofs are correct. The paper is honest about its main limitation: Theorem 7(ii) requires h''(u)A(u) to be symmetric, and in every application that symmetry comes from the detailed-balance condition, not from the equations themselves. The gap between normal ellipticity and entropy structure for nonsymmetric Onsager matrices is left open, and the paper says so explicitly in the discussion after Lemma 21. That is a limitation of scope, not a hidden flaw.\n\nThe one genuine error is in Remark 12. The formula det H = ∫ det(e^{-A^T t}e^{-A t})dt is not valid; determinants do not pass through integrals. For a concrete check, take A = diag(λ1, λ2), which gives det H = 1/(4λ1λ2) while the integral equals 1/(2(λ1+λ2)). This is a localized mistake in a remark, and nothing downstream depends on it, but it should be corrected or deleted.\n\nThe weakest part of the paper is the sufficiency side. The general sufficient condition rests on a structural symmetry assumption, so the title's \"when\" is answered only within the symmetric-Onsager class. The authors are upfront about this, and it doesn't undermine the main theorems. The citation pattern is appropriate: the factorization results are partly known, but the paper credits them and applies them to genuinely new cases.\n\nMy recommendation: send it to peer review. A competent referee can verify the main theorems quickly, and the paper will be useful to anyone working on entropy methods for cross-diffusion systems. I'd ask for a correction to Remark 12 and maybe a slightly more modest title, but the body holds up.","headline":"A clean, correct paper on entropy structure for cross-diffusion; the n-species SKT normal ellipticity result is the real news, but the general characterization only runs inside the detailed-balance class and Remark 12 has a false determinant identity.","tokens_in":20938,"tokens_out":3362,"would_cite":true,"duration_ms":33597,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35K40","35K55","35Q92","35Q79","15A23","15A24"],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper proves that normal ellipticity plus a symmetric Onsager matrix guarantees an entropy structure.","keywords":["Cross diffusion","entropy method","normal ellipticity","matrix factorization","Lyapunov equation","population model","volume-filling model","fluid mixture model"],"falsifier":"Take the 3-species SKT-type matrix from Lemma 21, A(u) = [[u3, 0, u1], [u2, u1, 0], [0, u3, u2]], and determine whether the system has an entropy structure. A proof that an entropy exists would show that Theorem 7(ii)'s symmetry assumption is not necessary; a proof that none exists would confirm the gap between normal ellipticity and entropy structure.","tokens_in":20004,"feed_emoji":"","tokens_out":9619,"duration_ms":84998,"temperature":0.7,"pith_summary":"Cross-diffusion systems couple the gradient of one species into the flux of another, and their diffusion matrices are typically neither symmetric nor positive definite, which blocks the usual energy estimates. This paper asks when such a system can be rewritten in entropy variables so that the transformed diffusion matrix becomes positive definite, yielding a Lyapunov functional and the gradient estimates needed for global existence. The answer is a matrix-factorization criterion: normal ellipticity, meaning all eigenvalues of the diffusion matrix have positive real part, is necessary for an entropy structure, and it becomes sufficient once there is a strictly convex entropy density making the Onsager matrix symmetric. For constant or nearly constant diffusion matrices, the paper shows that entropy structure and normal ellipticity are equivalent. This matters because the entropy structure is the main analytical tool for proving global-in-time existence of solutions to these systems.","feed_headline":"A symmetric Onsager matrix turns normal ellipticity into an entropy","feed_subtitle":"This condition supplies the Lyapunov functional and gradient estimates needed for global existence.","key_machinery":"The central mechanism is the factorization A(u) = A1A2 or A(u) = A2A1, with A1 symmetric positive definite. A matrix is normally elliptic when all eigenvalues have positive real part, and Lyapunov's matrix theorem says such matrices are exactly those admitting a factorization with A2 positive definite. If A2 can be chosen symmetric, the factorization forces A2 to be positive definite whenever A is normally elliptic, and taking A1 = h''(u)^{-1} turns A2 into h''(u)A(u). This factorization is the bridge between normal ellipticity and entropy structure; the detailed-balance condition pi_i partial p_i/partial u_j = pi_j partial p_j/partial u_i supplies the symmetry of the Onsager matrix in the model applications.","core_discovery":"The central claim is Theorem 7: if the system has an entropy structure with strictly convex entropy density h, then the diffusion matrix A(u) is normally elliptic for every u; conversely, if A(u) is normally elliptic and there exists a strictly convex h in $C^{2}$(D) with h''(u)A(u) symmetric for all u, then h''(u)A(u) is positive definite, so the system has an entropy structure. If in addition h''(u)A(u) is symmetric, A(u) is diagonalizable with positive eigenvalues. The proof factorizes A(u) = A1A2 with A1 = h''(u)^{-1} symmetric positive definite and A2 = h''(u)A(u); Lyapunov's matrix theorem and inertia arguments convert normal ellipticity of A into positive definiteness of A2. Applications show that the criterion reproduces the known Boltzmann entropies for the SKT population model and the volume-filling model, and it produces two entropies, one Boltzmann-type and one quadratic, for fluid mixtures with detailed balance. The paper also proves normal ellipticity of the n-species SKT diffusion matrix without detailed balance, a property not previously established.","pith_inferences":["The determinant formula det H = 1/(2 tr A) for the constant-matrix Lyapunov function suggests that entropy production rates for nearly constant diffusion matrices could be estimated from the trace alone; a testable extension would compare such bounds with numerical solutions for slowly varying A(u).","The Poincare-lemma construction behind the second entropy indicates that any flux function F with invertible Jacobian and curl-free pi_i F_i field admits an entropy, which could be checked for reaction-cross-diffusion systems with lower-order terms.","The paper leaves open whether a normally elliptic matrix that is not diagonalizable with positive eigenvalues can still support an entropy structure with a nonsymmetric Onsager matrix; the 3-species SKT-type example from Lemma 21 is a concrete test case."],"forward_implications":["For any system satisfying Theorem 7(ii), the function t maps to integral of h(u(t)) is nonincreasing and the production term supplies gradient estimates, opening the boundedness-by-entropy route to global existence.","For constant diffusion matrices, normal ellipticity alone is equivalent to an entropy structure, with an explicit Lyapunov function h(u) = (1/2) u^T H u and H given by an integral of exp(-A^T t) exp(-A t).","A constant normally elliptic matrix perturbed by a small bounded nonlinear term still admits an entropy structure, so the criterion is stable under perturbation.","When the entropy density is a sum of single-variable functions and h''A is symmetric, positive definiteness is equivalent to positivity of all leading principal minors of A, avoiding eigenvalue computations.","For fluid mixtures with detailed balance and an invertible pressure Jacobian, two distinct entropy densities exist, one Boltzmann-type and one quadratic, giving complementary gradient estimates."],"supporting_citations":[{"why":"Supplies the baseline result that normal ellipticity is the minimal condition for local existence of smooth solutions.","marker":"[1]"},{"why":"Provides the theorem that a matrix is diagonalizable iff it factorizes as a symmetric positive definite times a symmetric matrix, used in Proposition 3 and Theorem 7.","marker":"[5]"},{"why":"Contains the Lyapunov matrix theorems that underlie the factorization and normal-ellipticity characterizations in Section 2.","marker":"[15]"},{"why":"Establishes global existence and the Boltzmann entropy for the n-species SKT model under detailed balance, the main example whose entropy is reconstructed.","marker":"[9]"},{"why":"Derives the volume-filling model and proves its entropy, the template for the factorization in Case 2.2.","marker":"[28]"},{"why":"Derives the fluid-mixture model with partial pressure gradients from a particle system, the setting for the two new entropies in Section 6.","marker":"[8]"},{"why":"Introduces the SKT population model whose n-species diffusion matrix is proved normally elliptic in Section 7.","marker":"[27]"}],"fun_headline_variants":["Normal ellipticity plus Onsager symmetry gives entropy structure","Entropy structure iff normal elliptic for nearly constant diffusion","Matrix factorization yields new entropy criteria for cross-diffusion","SKT population model: normal ellipticity proven without detailed balance"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The sufficient direction of the main theorem assumes that a strictly convex entropy density exists making the Onsager matrix symmetric on the whole domain; in the model examples this symmetry is the detailed-balance condition on the coefficients, a structural restriction not derived from the equations.","fun_headline_variants_meta":{"raw":{"variants":["Normal ellipticity plus Onsager symmetry gives entropy structure","Entropy structure iff normal elliptic for nearly constant diffusion","Matrix factorization yields new entropy criteria for cross-diffusion","SKT population model: normal ellipticity proven without detailed balance"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000523,"raw_usage":{"total_tokens":2527,"prompt_tokens":944,"completion_tokens":1583,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":560,"completion_tokens_details":{"reasoning_tokens":1516}},"tokens_in":560,"tokens_out":1583,"duration_ms":15144,"temperature":1.0,"reasoning_tokens":1516,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:32:27.674887+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take the 3-species SKT-type matrix from Lemma 21, A(u) = [[u3, 0, u1], [u2, u1, 0], [0, u3, u2]], and determine whether the system has an entropy structure. A proof that an entropy exists would show that Theorem 7(ii)'s symmetry assumption is not necessary; a proof that none exists would confirm the gap between normal ellipticity and entropy structure.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the baseline result that normal ellipticity is the minimal condition for local existence of smooth solutions."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the theorem that a matrix is diagonalizable iff it factorizes as a symmetric positive definite times a symmetric matrix, used in Proposition 3 and Theorem 7."},{"cited_title":"Horn and C","cited_arxiv_id":null,"evidence_quote":"Contains the Lyapunov matrix theorems that underlie the factorization and normal-ellipticity characterizations in Section 2."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes global existence and the Boltzmann entropy for the n-species SKT model under detailed balance, the main example whose entropy is reconstructed."},{"cited_title":"Zamponi and A","cited_arxiv_id":null,"evidence_quote":"Derives the volume-filling model and proves its entropy, the template for the factorization in Case 2.2."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Derives the fluid-mixture model with partial pressure gradients from a particle system, the setting for the two new entropies in Section 6."},{"cited_title":"Shigesada, K","cited_arxiv_id":null,"evidence_quote":"Introduces the SKT population model whose n-species diffusion matrix is proved normally elliptic in Section 7."}],"review_version":1}