{"id":"02bf03da-1c6b-4e3a-a371-9dac738a523e","arxiv_id":"1908.01811","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Existence of weak solutions is proved for a large-strain elastodynamics model with repulsive long-range monopolar interactions, with and without diffusing charge, under a p-regularized electrostatic equation.","lead":"Charged elastic or poroelastic bodies can deform so much that parts of them overlap, which breaks the standard mathematics of long-range electric forces. This paper adds a nonlocal smoothing energy and uses an area formula to prove that the resulting dynamical equations have solutions, with the caveat that electrostatics is p-regularized rather than standard.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Equation (72) in Proposition 2 omits the ∇(φ(χ)) term when inverting ∂²_qqϕ, so the L²(H¹) estimate for q and the strong W^{1,1} convergence of q_k(t) are not proven as written.","rationale":"The paper's most valuable and least secure claim is the existence theorem for the diffusive system (Proposition 2). I agree with the reader that the Healey-Kroemer determinant bound is foundational and that the proof of Proposition 1 has terse passages; however, that bound is backed by a published theorem and the intended assumptions are recoverable from the text, so I do not see a concrete route by which it fails. The concrete defect I find is in the proof of Proposition 2: identity (72) is obtained by differentiating (63a), but the term ∇χ^T(∇φ∘χ) arising from the φ(χ) summand is dropped. That identity is the only mechanism offered to control ∇q_k in L² and to pass to the limit in terms involving ∇(qζ). As written, the estimate (70a) for q_k and the strong convergence q_k(t)→q(t) in W^{1,1}(Ω) do not follow. The missing term is bounded by the available estimates, so the gap is probably repairable, which is why I recommend keeping the reader's CONDITIONAL verdict rather than rejecting. A referee should require the corrected identity and a proof of strong convergence of the added term before the claim is accepted.","tokens_in":21550,"tokens_out":50015,"duration_ms":517065,"concrete_test":"Re-derive (72) from (63a) and the variational inequality (69), keeping the term ∇χ^T(∇φ∘χ). Then verify the two claims needed to repair the proof: (i) the added term is uniformly bounded in L²(Q) using the a-priori bounds on ∇φ_k and ∇χ_k; (ii) for a.a. t, ∇φ_k(t,χ_k(t))→∇φ(t,χ(t)) strongly in L¹(Ω), using φ_k→φ in C(B_R), ∇φ_k→∇φ in L^p(B_R), χ_k→χ in C¹(Ω), and a uniform bound on the multiplicity card χ_k^{-1}. If (i) or (ii) fails, Proposition 2's L²(H¹) estimate for q and the strong convergence of q_k(t) in W^{1,1} are not established.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The proof of Proposition 2 derives the estimate of ∇q_k from identity (72), obtained by applying ∇ to the inclusion (63a), µ∈∂_qϕ(∇χ,q)+φ(χ), and solving for ∇q. Correct differentiation gives ∇µ=∂²_Fqϕ∇²χ+∂²_qqϕ∇q+∇χ^T(∇φ∘χ), so (72) is missing the term ∇χ^T(∇φ∘χ). This is not cosmetic: (72) is the only step that yields the L²(I;H¹(Ω)) estimate for q_k in (70a), and the subsequent strong convergence q_k(t)→q(t) in W^{1,1}(Ω) for a.a. t is needed to pass to the limit in the ∇(qζ) term of (34a) and in (61b). Without a corrected identity, the estimate and the limit passage for the diffusive system are unsupported. The missing term is likely controllable, since ∇φ_k is bounded in L∞(I;L^p(R^d)) with p>d and ∇χ_k is bounded in L∞(I;L∞(Ω)), giving the product in L∞(I;L^p(Ω))⊂L²(Q); and ∇φ_k(χ_k)→∇φ(χ) strongly in L¹ should follow from the available strong convergences. But the paper does not provide this argument, and as written the proof contains a false identity.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a continuum model for the dynamics of charged elastic (and poroelastic) bodies at large strains, allowing for possible global self-interpenetration while retaining local injectivity. The key modeling ingredients are a nonlocal nonsimple free energy with singular kernel, a p-regularized nonlinear Poisson equation for the electrostatic potential, and Federer's area formula to define the Eulerian charge density through a set-valued inverse deformation. The analytical core consists of two Galerkin-based existence theorems: Proposition 1 establishes weak solutions to the elastodynamics system with fixed charge, and Proposition 2 extends this to the case of a diffusing charge driven by an electrochemical potential. The proofs rely crucially on the Healey-Kroemer determinant lower bound, on strong compactness of the deformation in C^1 obtained from the nonlocal second-gradient regularization, and on strong convergence of the scalar potential in C(B) for a.e. time.","tokens_in":21890,"tokens_out":12467,"duration_ms":129074,"significance":"If the technical gaps are repaired, the paper is a valuable contribution to the mathematical theory of electro-elastodynamics at large strains. Its use of Federer's change-of-variable formula to handle non-injective deformations is a genuine novelty with respect to earlier work, and the constructive Galerkin approach with explicit a priori estimates is appropriate for the hyperbolic-elliptic structure. The paper honestly states its limitations (monopolar interactions only, ideal-dielectric assumption, p-regularization instead of standard electrostatics), and it does not fit parameters or assume the desired conclusion. The main theorems are plausible and the overall framework is likely to be useful for further work on coupled electro-chemo-mechanics, but the proofs contain load-bearing gaps that need to be addressed.","major_comments":[{"comment":"Equation (72) is not the correct derivative of (63a). Since µ_k ∈ ∂_qϕ(∇χ_k,q_k)+φ_k(χ_k), spatial differentiation gives ∇µ_k = ∂²_Fqϕ(∇χ_k,q_k)∇²χ_k + ∂²_qqϕ(∇χ_k,q_k)∇q_k + ∇χ_k^T(∇φ_k∘χ_k). The last term is omitted in (72). This is load-bearing: (72) is the only step that yields the L²(I;H¹(Ω)) estimate for q_k in (70a), and the subsequent strong convergence ∇q_k(t)→∇q(t) in L¹ is used in the passage to the limit after (72). The missing term is likely controllable, since ∇φ_k is bounded in L∞(I;L^p(R^d)) with p>d and ∇χ_k is bounded in L∞(I;L∞(Ω)), so the product lies in L²(Q); however this argument is not given. As written, the proof of Proposition 2 is incomplete.","section":"§4, Eq. (72)"},{"comment":"The assertion that φ_k(t)→φ(t) strongly in C(B) for a.a. t∈(0,T) and every closed ball B is not justified in the text. The preceding compactness statements give only (51b) and the pointwise strong convergence (56) for χ_k(t). To obtain (57), one must show that a single subsequence works for all t outside a null set, using the continuity of the solution map of the monotone equation (58) as the right-hand side converges in W^{1,p}(R^d)^*, together with uniqueness of φ(t). This is a standard argument, but it should be written out explicitly because the limit passage in the composed terms (61a)–(61c) depends on it.","section":"§3, proof of Proposition 1, around (57)"},{"comment":"The successive-continuation argument leading to the uniform determinant lower bound det∇χ_k ≥ ε on the whole time interval is only sketched. The energy identity (43) and the Healey-Kroemer implication (40) are invoked before the global a priori estimate is established, so the proof should explain how the maximal time of existence is shown to be T: on any interval where the local Galerkin solution satisfies the energy bound, (40) prevents the determinant from approaching zero and hence allows continuation. Without this bootstrapping detail, the uniform bound (42) that underpins the change-of-variable formula and the control of ∇χ_k^{-⊤} is not fully justified.","section":"§3, Eq. (42)"}],"minor_comments":[{"comment":"In the displayed inequality following (48), the right-hand side is missing a dt: the term ∫_0^τ ‖χ̇_k(t)‖_{L²} should read ∫_0^τ ‖χ̇_k(t)‖_{L²} dt, and the displayed estimate should also account for the term ‖χ_k(t)‖²_{H^{2+γ}} dt that appears in the following line.","section":"§3, before Eq. (50)"},{"comment":"The notation ‖‖‖ 1/det(∇χ_k) ‖‖‖ contains a typographical artifact; it should presumably be the standard L∞(Q) norm.","section":"§3, Eq. (41a)"},{"comment":"The test function for the scalar mass-balance equation is written w∈H¹(Q;R^d), but w is scalar (it multiplies the scalar quantities M∇µ and q); it should likely be H¹(Q) or W¹,¹(Q) with w(T)=0.","section":"Definition 4.1, Eq. (68)"},{"comment":"There are several typographical and stylistic slips, such as \"Lagrangean\" instead of \"Lagrangian\", \"an result\", and inconsistent uses of upright and italic fonts for Eulerian and Lagrangian quantities. These should be corrected in the revision.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The false identity in Eq. (72) is the most serious technical defect, but the missing term appears controllable with the estimates already at hand, so the result is likely repairable rather than fundamentally flawed. The referees' main concern should be that the authors supply a complete argument for (57) and for the successive-continuation step in Proposition 1, and correct (72) together with the subsequent convergence proof. The paper's modeling framework and use of the Federer area formula are original and merit consideration after revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThis is a real working paper, not a gesture. What's new: it extends Roubíček and Tomassetti's previous work to dynamical problems with monopolar long-range self-interactions, where global invertibility of the deformation is not assumed. The trick is to use the nonlocal nonsimple regularization to get H^{2+γ} regularity, then Healey-Kroemer to keep det∇χ uniformly positive, and Federer's area formula to make sense of the pullback charge density under non-injective Lipschitz deformations. That combination is legitimate and the modeling discussion is honest about the ideal-dielectric assumption and p-regularization being limitations. I read the paper as doing what it claims: existence of weak solutions for the p-regularized system, with a serious Galerkin construction.\n\nWhere I'd push back. Proposition 1's proof is terse in places. The strong almost-everywhere convergence φ_k(t)→φ(t) in C(B) at (57) is asserted without a full argument; you only have φ_k bounded in L∞(I;W^{1,p}(R^d)) with p>d, and the equation gives some uniform continuity in time only in weaker norms. It's probably repairable by using the equation and the strong convergence of χ_k, but as printed it's a gap in the presentation.\n\nProposition 2 has a more concrete problem. Identity (72) is wrong: differentiating (63a) gives an extra ∇χ^T(∇φ∘χ) term when solving for ∇q. This is the step that yields the L²(I;H¹(Ω)) bound for q_k and the strong W^{1,1} convergence of q_k(t), both needed for the limit passage in the diffusive case. The missing term looks controllable — ∇φ_k is bounded in L∞(I;L^p) with p>d and ∇χ_k in L∞(I;L∞), so the product lands in L²(Q), and strong convergence should follow from the available convergences — but the paper doesn't give that argument. So as written, Proposition 2's proof has a load-bearing gap, not just a typo.\n\nThe negative claims are fine: no fitted parameters, no circularity, and the self-citations are to the authors' own earlier work that establishes the framework. The p-regularization is presented as a modeling compromise, which it is.\n\nBottom line: Proposition 1 is likely sound after some filling-in; Proposition 2 needs a corrected identity and the extra estimate. This is a subfield-level result but it's the kind of thing people doing numerical large-strain electromechanics will want to cite. I'd send it to a serious referee and ask them to check the diffusive section carefully, but I wouldn't desk-reject.","headline":"Solid extension of the authors' program to dynamic long-range self-interactions; Proposition 2 has a missing term in a key identity that needs fixing, but the paper deserves referee time.","tokens_in":22364,"tokens_out":1710,"would_cite":false,"duration_ms":16550,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35Q60","35Q74","65M60","74A30","74F15","76S99","78A30"],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper proves existence of weak solutions for the dynamics of charged elastic and poroelastic bodies at large strains, allowing global self-penetration.","keywords":["elastodynamics","poroelasticity","nonsimple materials","long-range interactions","electrostatic self-interactions","weak solutions","Galerkin approximation","large strains"],"falsifier":"A numerical or analytical example satisfying (29a)-(c), (30), and (31) in which the Galerkin solutions develop a point or time where det∇χ_k(t,x) tends to 0 while the kinetic plus stored energy remains bounded would refute the determinant bound (42) on which Proposition 1 rests; equivalently, exhibiting a sequence χ_k with bounded Φ(χ_k) but no uniform positive lower bound on det∇χ_k would falsify the Healey-Kroemer form (40) used here.","tokens_in":21356,"feed_emoji":"⚡","tokens_out":7072,"duration_ms":67176,"temperature":0.7,"pith_summary":"The paper proves existence of weak solutions for a dynamical model of an elastic body carrying electric charge and undergoing arbitrarily large strains. The central difficulty is that the electrostatic potential lives on the current, deformed configuration while the mechanics is formulated on the reference configuration, and the charge density involves the inverse deformation, which need not be globally injective in dynamics. The authors resolve this by admitting only local non-self-penetration, using Federer's change-of-variables formula to define the charge density through the preimage of the deformation, and adding a nonlocal second-gradient energy whose singular kernel forces deformations into a fractional Sobolev space. A theorem of Healey and Kroemer then gives a uniform positive lower bound on the determinant of the deformation gradient, which makes the change of variables legitimate and controls the pullback of the electric force. The main results, Propositions 1 and 2, show that Galerkin approximations exist on the whole time interval and their weak* limits solve the coupled system, including the diffusion of the charge in the poroelastic extension.","feed_headline":"Charged elastic bodies: existence of weak solutions at large strains","feed_subtitle":"Proof handles the electrostatic force even when the deformed body folds over itself.","key_machinery":"The central machinery is the nonlocal nonsimple energy (1)-(2) with a singular kernel satisfying (30), combined with the Healey-Kroemer determinant bound in the form (40): any deformation with bounded energy has det∇χ uniformly bounded below by a positive constant. This gives Lipschitz regularity and local invertibility, so Federer's change-of-variables formula (7) defines the Eulerian charge density (6) as a sum over the preimage points, and the force term q∇φ∘χ can be rewritten via (33) in a form that passes to the limit under only weak convergence of ∇²χ. The p-regularized Poisson equation (9) with p > d ensures the potential converges in the space of continuous functions on compact sets, which is needed for the composition φ∘χ.","core_discovery":"Proposition 1 establishes that the dynamical system (21)-(23), coupling elastodynamics with a p-regularized Poisson equation for the electrostatic potential, admits a weak solution in the sense of Definition 3.1. Proposition 2 extends this to the poroelastic case where the charge density q is an evolving diffusant, governed by the electrochemical potential inclusion (63a) and a Biot-type diffusion equation (63b,c), with the quadruple (χ, φ, q, μ) satisfying (34a), (34b), (68), and (69). The argument is constructive: nested Galerkin subspaces produce approximate solutions on the full time interval, uniform estimates (41a)-(41b) and (70) yield weakly* converging subsequences, and the determinant lower bound (42) supplies the strong convergence needed to pass to the limit in the terms containing ∇$χ^{{-T}}$ and φ∘χ.","pith_inferences":["If the determinant bound (40) holds beyond the assumptions made here, the same change-of-variables device could handle attractive monopolar interactions such as gravitation, but the paper explicitly leaves ε1 = 0 out, so Proposition 1 does not cover self-gravitating bodies.","The p-regularization of the Poisson equation is a technical device; the authors note as an open problem whether the standard electrostatic setting can be recovered, and a testable extension would be to attempt the same limit passage with p = 2 under stronger assumed regularity of the deformation.","The model's allowance of global self-penetration while retaining local injectivity suggests that charge conservation for overlapping deformation images is encoded in the multiplicity sum (6); one could test this by computing the electrostatic energy of a deformation that folds a charged body onto itself and comparing it with the unfolded configuration.","The multi-component version in Remark 3 merges fixed dopant charges with mobile diffusant charges, giving a direct route to polymer-electrolyte fuel cell and semiconductor device models, though the proof would need adaptation to the cross-diffusion mobility tensor."],"forward_implications":["For the purely elastic model of Section 2, a weak solution exists on any time interval for the initial-boundary-value problem, so the model can serve as a rigorous basis for engineering computations that ignore global self-penetration.","Including the charged diffusant gives weak solutions of the coupled elastodynamics-diffusion system, so Darcy, Fick, and drift responses such as (74)-(75) are obtained within a large-strain existence theory.","The Galerkin construction is explicit and the determinant bound rules out the Lavrentiev-type singularity on the approximation level, so the singular energy barrier at detF→0+ needs no additional regularization.","For the diffusion extension, the energy balance (66) holds, combining the conservative mechanical-electrostatic balance with the dissipative diffusion rate, which is exactly the structure the a priori estimates exploit."],"supporting_citations":[{"why":"Supplies the theorem, in the modified form (40), that bounded energy implies a uniform positive lower bound on det∇χ; this is the load-bearing determinant control.","marker":"[14]"},{"why":"Supplies the change-of-variables and area formula (7) that defines the Eulerian charge density without requiring global injectivity of the deformation.","marker":"[15]"},{"why":"Provides the classical Biot poroelastic consolidation model that underlies the large-strain diffusion energy (67).","marker":"[3]"},{"why":"Extends the Biot model to finite deformations of porous solids, framing the large-strain poroelastic setting used in Section 4.","marker":"[4]"},{"why":"Supplies the pullback formula (63d) for the mobility tensor under finite strain, which is used in the diffusion equation.","marker":"[8]"},{"why":"Provides the mathematical continuum-mechanics background and the modified Healey-Kroemer growth condition q > 2d/(2γ+2-d) used in (40).","marker":"[18]"},{"why":"The authors' previous analysis of magnetoelastic materials under diffusion at large strains, from which the Galerkin strategy and the nonsimple-material concept are carried over.","marker":"[27]"},{"why":"Justifies the ideal-dielectric assumption, which lets the free energy split additively into mechanical and electrostatic parts.","marker":"[31]"}],"fun_headline_variants":["Weak solutions exist for charged elastic bodies","Charged elastic bodies: weak solutions at large strains","Existence proof for charged elastic bodies with diffusion","Large-strain elasticity with charges: weak existence","Diffusion-driven charged elasticity: weak solutions"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"Everything rests on the assertion that any deformation with bounded energy keeps the determinant of its deformation gradient uniformly positive; if a bounded-energy motion could drive the determinant to zero, the charge-density formula and the entire weak formulation would lose meaning.","fun_headline_variants_meta":{"raw":{"variants":["Weak solutions exist for charged elastic bodies","Charged elastic bodies: weak solutions at large strains","Existence proof for charged elastic bodies with diffusion","Large-strain elasticity with charges: weak existence","Diffusion-driven charged elasticity: weak solutions"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000281,"raw_usage":{"total_tokens":1613,"prompt_tokens":843,"completion_tokens":770,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":459,"completion_tokens_details":{"reasoning_tokens":700}},"tokens_in":459,"tokens_out":770,"duration_ms":7794,"temperature":1.0,"reasoning_tokens":700,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T15:02:51.715791+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A numerical or analytical example satisfying (29a)-(c), (30), and (31) in which the Galerkin solutions develop a point or time where det∇χ_k(t,x) tends to 0 while the kinetic plus stored energy remains bounded would refute the determinant bound (42) on which Proposition 1 rests; equivalently, exhibiting a sequence χ_k with bounded Φ(χ_k) but no uniform positive lower bound on det∇χ_k would falsify the Healey-Kroemer form (40) used here.","supporting_citations":[{"cited_title":"Healey and S","cited_arxiv_id":null,"evidence_quote":"Supplies the theorem, in the modified form (40), that bounded energy implies a uniform positive lower bound on det∇χ; this is the load-bearing determinant control."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the change-of-variables and area formula (7) that defines the Eulerian charge density without requiring global injectivity of the deformation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the classical Biot poroelastic consolidation model that underlies the large-strain diffusion energy (67)."},{"cited_title":"Biot and G","cited_arxiv_id":null,"evidence_quote":"Extends the Biot model to finite deformations of porous solids, framing the large-strain poroelastic setting used in Section 4."},{"cited_title":"Duda, A.C","cited_arxiv_id":null,"evidence_quote":"Supplies the pullback formula (63d) for the mobility tensor under finite strain, which is used in the diffusion equation."},{"cited_title":"Kruˇ z´ ık and T","cited_arxiv_id":null,"evidence_quote":"Provides the mathematical continuum-mechanics background and the modified Healey-Kroemer growth condition q > 2d/(2γ+2-d) used in (40)."},{"cited_title":"Roub´ ıˇ cek and G","cited_arxiv_id":null,"evidence_quote":"The authors' previous analysis of magnetoelastic materials under diffusion at large strains, from which the Galerkin strategy and the nonsimple-material concept are carried over."},{"cited_title":"Zhao and Z","cited_arxiv_id":null,"evidence_quote":"Justifies the ideal-dielectric assumption, which lets the free energy split additively into mechanical and electrostatic parts."}],"review_version":1}