{"id":"31bc6697-e581-4aed-85cd-f8a58ef12172","arxiv_id":"2501.14930","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"The authors prove that linear boundary port-Hamiltonian systems on one-dimensional time-varying spatial domains admit a time-varying Stokes-Dirac structure, defining a new class of moving-boundary port-Hamiltonian systems.","lead":"This paper gives a port-Hamiltonian formulation for distributed systems living on a spatial interval whose boundaries move in time, adding moving-boundary power terms and a time-varying Dirac structure. It then uses this structure to build a moving-grid discretization for the telegrapher's equations and tests it numerically.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Theorem 1 fails when a boundary velocity vanishes: D̂(t)^⊥ strictly contains D̂(t), so the central Dirac-structure claim is false as stated.","rationale":"The reader's weakest assumption correctly identified Assumption 1(iii) as fragile, but the more serious and concrete defect is the zero-velocity degeneracy: the proof of Theorem 1 breaks down when √ȧ or √ḃ vanishes, and the theorem is actually false in that case. This is not a matter of restrictive scope; Assumption 1 explicitly allows zero velocities, and the numerical example includes instants with ȧ > 0, ḃ = 0 and with ȧ = ḃ = 0. Because Theorem 1 is the central mathematical result from which Proposition 3 and the claimed new class of moving-boundary port-Hamiltonian systems follow, a false theorem under stated assumptions invalidates the main contribution as written. The argument may be salvageable by requiring strictly nonzero boundary velocities (ȧ(t)ḃ(t) > 0) and by reconciling the Hermitian boundary pairing with the bilinear Dirac-structure definition, but the current version does not support the stated conclusion. Hence the verdict should move from CONDITIONAL to REJECT, with the possibility of resubmission after a corrected assumption set and a revised proof.","tokens_in":22039,"tokens_out":23141,"duration_ms":201726,"concrete_test":"Specialize the construction to the lossless transmission line with Q = I, L = 1, ȧ = 0, ḃ = 1, and check the element γ1 = (φ = 0, φ∂1 = 0, η∂1 = 0, φ∂2 = 1, η∂2 = −1/2, η = 0). Verify that the symmetric pairing of γ1 with every γ ∈ D̂(0) is zero, then verify that γ1 does not satisfy the second-block port equations in the definition of D̂(0). If both checks hold, D̂(0)^⊥ strictly contains D̂(0), contradicting Theorem 1.","verdict_should_be":"REJECT","load_bearing_attack":"Theorem 1 is false as stated. Assumption 1(iii) admits ȧ·ḃ ≥ 0, which permits one boundary velocity to be zero, but the proof's final inference in Appendix C requires both √ȧ and √ḃ to be nonzero. In Eq. (C.12), if α = √ȧ = 0 and β = √ḃ > 0, the equation reduces to φ∂2 = −2η∂2 (when η = 0), which does not imply the port equations φ∂2 = 0 and η∂2 = 0 required by (13). A concrete counterexample with n = 1, Q = 1, α = 0, β = 1 is the element γ1 = (φ = 0, φ∂1 = 0, η∂1 = 0, φ∂2 = 1, η∂2 = −1/2, η = 0). The pairing of γ1 with every element of D̂(t) vanishes, so γ1 ∈ D̂(t)^⊥, but γ1 violates the second-block port equations of (13), so γ1 ∉ D̂(t). Thus D̂(t)^⊥ ≠ D̂(t). This degeneracy occurs exactly at static-boundary instants, including t = 0 and t > 7.5 in the paper's own transmission-line numerical example, so the claimed time-varying Dirac structure is not a Dirac structure at those times.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a port-Hamiltonian formulation for linear boundary port-Hamiltonian distributed parameter systems on a one-dimensional time-varying interval. It introduces a normalized coordinate, derives transformed state dynamics (Lemma 2, Eq. (8)), a power balance (the second Lemma 2, Eq. (10)), and a time-varying Dirac structure (Theorem 1) with boundary ports (13), from which it defines moving-boundary port-Hamiltonian systems (Proposition 3). The framework is then applied to a dynamic-mesh discretization of the telegrapher's equations, with numerical validation. The central contribution is the Dirac-structure characterization; the discretization is presented as a consequence of that structure.","tokens_in":22309,"tokens_out":14521,"duration_ms":162982,"significance":"If the Dirac-structure result were correct under the stated assumptions, this would be a genuinely useful contribution: it gives a concrete, physically interpretable way to preserve power balance under moving boundaries and a principled basis for dynamic-mesh discretization. The chain-rule derivation of the dynamics and the energy-balance computation are transparent, and the telegrapher example is well chosen. The authors also deserve credit for explicitly constructing time-varying ports and for reporting the non-preservation of the discretized power balance rather than hiding it. However, the main theorem is not correct as stated: the proof degenerates at instants where one boundary velocity vanishes, and Assumption 1 as written also allows cases not covered by the square-root identities used in the proof. The significance of the paper therefore depends on repairing this degeneracy and clarifying the admissible velocity regimes.","major_comments":[{"comment":"The Dirac-structure claim fails when one boundary velocity vanishes, which is allowed by Assumption 1(iii). The final step of the proof infers the second-block port equations from Eq. (C.12), but if, say, a_dot(t)=0 and b_dot(t)>0, Eq. (C.12) yields only a one-dimensional relation and does not imply the two port equalities in (13). A concrete counterexample is n=1, Q=1, a_dot=0, b_dot=1, and the element gamma1 = (phi=0, phi_d1=0, eta_d1=0, phi_d2=1, eta_d2=-1/2, eta=0). For every element gamma=(f, f_d, e, e_d) in D_hat(t), the second-block port equations give e_d2 = (1/2) f_d2, so the pairing of gamma1 with gamma vanishes; hence gamma1 is in D_hat(t)^perp. But gamma1 does not satisfy the port equations of (13), so gamma1 is not in D_hat(t). Thus D_hat(t)^perp is strictly larger than D_hat(t). This is not a remote corner case: in the numerical example of Section 4.2, b_dot(0)=0 and both velocities are zero for t>7.5. The theorem should either assume strictly nonzero velocities of the same sign and treat static instants separately, or the port construction must be modified to remove the degeneracy. As stated, the central claim is false.","section":"Theorem 1 and Appendix C, Eq. (C.12)"},{"comment":"Assumption 1(i) states only that a,b are continuous, but the entire derivation differentiates them: Eq. (8), Theorem 1, and the boundary-port definition (13) all use a_dot and b_dot. The assumption should be C^1, or piecewise C^1 with the theorem stated on each differentiability interval. Moreover, the numerical example in Section 4.2 defines b(t) piecewise with a jump in b_dot at t=7.5; as written, the theorem is asserted at that instant even though the boundary velocity is not defined there. This is a load-bearing regularity gap, not merely a presentation issue.","section":"Assumption 1(i), Lemma 2, Theorem 1"},{"comment":"The proof uses the identity sqrt(a_dot) sqrt(b_dot) = sqrt(a_dot b_dot). With the standard principal branch this identity holds only when both velocities are nonnegative. Assumption 1(iii) as written also permits a_dot<0 and b_dot<0, for which sqrt(a_dot) sqrt(b_dot) = -sqrt(a_dot b_dot). If the intended scope is expanding or contracting domains with nonnegative velocities, Assumption 1(iii) must be replaced by a_dot, b_dot >= 0; otherwise the proof must specify a branch convention and handle the sign. As stated, the Dirac-structure proof does not cover the full range of Assumption 1.","section":"Appendix C, Eqs. (C.3) and (C.11)"}],"minor_comments":[{"comment":"The power-balance result is also labelled Lemma 2, duplicating the dynamics lemma of Section 3.1; the second one should be renumbered and cross-references updated.","section":"Section 3.3"},{"comment":"The definition of phi_hat(t, s_hat) appears to omit the factor sqrt(b(t)-a(t)): it is written as phi(t, h(t, s_hat)), while Eq. (4) applies the factor to all components of x. This is inconsistent with the Hamiltonian scaling and with the energy-balance expressions that follow.","section":"Section 3.2"},{"comment":"In Eq. (13), the sentence about z denoting the complex conjugate of z is unclear because no overlines appear in the displayed matrix. It should be stated explicitly which entries involving sqrt(a_dot) and sqrt(b_dot) are conjugated, and for negative velocities the branch of the square root must be specified.","section":"Section 3.5.2 and Theorem 1"},{"comment":"The set comprehension in Theorem 1 has a typo: it should read (f_hat, f_d_hat, e_hat, e_d_hat) in B_hat(t), not (f_hat, f_d_hat, e_hat, e_d_hat in B_hat(t)).","section":"Theorem 1 statement"}],"recommendation":"major_revision","confidential_remarks":"The paper is in scope and the core idea is good, but the main theorem is overstated and the stress-test concern is real. I would ask the authors to strengthen Assumption 1 or add a separate treatment of zero-velocity instants, align the numerical example with the corrected assumptions, and re-check the negative-velocity case. If those points are addressed, the paper could be suitable for publication; in the current form it should not be accepted as is."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Take a look at arXiv:2501.14930. The authors extend port-Hamiltonian systems to one-dimensional domains with moving boundaries, and the first half is genuinely useful: the change of variables, the state dynamics in (8), and the power balance in (10) are clean, and they correctly identify the earlier moving-interface papers as special cases. The Dirac structure idea is the right one.\n\nThe problem is Theorem 1. As stated, it is false. Assumption 1(iii) allows one boundary velocity to be zero, and the proof of D̂⊥ ⊂ D̂ in Appendix C breaks exactly there. The step from (C.12) to the port equations requires both √ȧ and √ḃ to be nonzero; if ȧ=0, the bilinear pairing only forces a relation between φ∂2 and η∂2, not the individual port equations. Concretely, with n=1, Q=1, ȧ=0, ḃ=1, the element γ1=(φ=0, φ∂1=0, η∂1=0, φ∂2=1, η∂2=-1/2, η=0) is in D̂⊥ but violates the port equations, so D̂⊥ ≠ D̂. This is not a pathological corner: their own transmission-line example has ḃ(0)=0 and both velocities zero for t>7.5.\n\nThe fix is probably to require ȧḃ>0 strictly, or to treat zero-velocity instants with a degenerate structure separately. Either way, the paper overclaims as written. The theorem is likely correct on the open set where both boundaries move in the same direction at nonzero speed, and the construction is a real contribution. The minor issues—Assumption 1(i) says continuity but uses derivatives, and the numerical section is a single illustrative run with no convergence study or code—are secondary. The discretization honesty is a plus: they state clearly it is not structure-preserving.\n\nI'd send this to referees; the flaw is concrete and fixable, and the result is significant enough to warrant the effort. The verdict should be revise before accept, not reject.","headline":"The moving-boundary port-Hamiltonian construction is a real contribution and likely correct when both boundary velocities are nonzero, but Theorem 1 as stated is false at instants where one velocity vanishes.","tokens_in":22853,"tokens_out":6051,"would_cite":false,"duration_ms":46929,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["37K05","93C20","65M50","35Q93"],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper proves that a linear port-Hamiltonian system on a time-varying spatial interval admits a time-varying Dirac structure, formally establishing moving-boundary port-Hamiltonian systems as a generalization of the fixed-domain theory.","keywords":["port-Hamiltonian systems","Dirac structure","moving boundary","free boundary problems","distributed parameter systems","transport equation","spatial discretization","telegrapher's equations"],"falsifier":"Evaluate the Dirac orthogonality condition $\\hat D(t) = \\hat D^\\perp(t)$ symbolically for the ports (13) under $\\dot a(t) < 0 < \\dot b(t)$: the substitution $\\sqrt{\\dot a}\\sqrt{\\dot b} = \\sqrt{\\dot a\\dot b}$ used in Appendix C ceases to hold over the reals, so a direct check of whether the subspace is self-orthogonal under the complex pairing settles whether Theorem 1 survives outside Assumption 1(iii).","tokens_in":21821,"feed_emoji":"📐","tokens_out":9122,"duration_ms":69831,"temperature":0.7,"pith_summary":"This paper aims to extend boundary port-Hamiltonian systems — PDE models whose energy exchange with the environment is encoded in a Dirac structure — from fixed spatial domains to domains whose endpoints move over time. The authors prove that a linear first-order port-Hamiltonian PDE on a time-varying interval $[a(t), b(t)]$ still carries a Dirac structure, provided the domain is re-parameterized onto a fixed reference interval and the boundary motion is folded into the port variables. This yields the new class of moving-boundary port-Hamiltonian systems, which reduces to the classical definition when the boundaries are stationary. The result matters because moving domains arise in phase-transition problems, deforming mechanical systems, and lithography equipment, where ignoring conservation structure is known to damage controller design and simulation. The paper also turns the structure into a dynamic-mesh spatial discretization for the telegrapher's equations and checks it numerically.","feed_headline":"Moving-boundary PDEs get a port-Hamiltonian structure","feed_subtitle":"A time-varying Dirac structure captures energy flow through moving endpoints, enabling dynamic-mesh simulation.","key_machinery":"The load-bearing object is the time-varying Stokes–Dirac structure $\\hat D(t)$ on the fixed reference domain $\\hat s \\in [0,1]$, obtained by the change of spatial coordinates $s = a(t) + (b(t)-a(t))\\hat s$ and the state rescaling $\\hat x = \\sqrt{b-a}\\,x$. The coordinate change converts the moving endpoints into two new flow terms: a compression term $-\\frac{1}{2}(\\dot b-\\dot a)/(b-a)\\hat x$ and a translation term $\\partial_{\\hat s}\\big((\\dot a + (\\dot b-\\dot a)\\hat s)/(b-a)\\hat x\\big)$, which the paper interprets physically and packages into the port equations (13). The proof of the Dirac property combines Stokes' theorem, which handles the internal Hamiltonian operator, with the Leibniz integral rule, which handles the moving integration bounds; the price of this packaging is that the boundary ports take complex values whenever a boundary velocity is negative, with the pairing $\\hat f_\\partial^H \\hat e_\\partial$ remaining real under Assumption 1(iii).","core_discovery":"On the paper's own terms, the central discovery is Theorem 1: for a linear boundary port-Hamiltonian system on a time-varying interval $[a(t),b(t)]$ whose endpoints satisfy Assumption 1, the time-dependent subspace $\\hat D(t)$ defined by the flow equation $\\hat f = -\\frac{1}{2}\\frac{\\dot b-\\dot a}{b-a}Q^{-1}\\hat e + \\partial_{\\hat s}\\big(\\frac{\\dot a+(\\dot b-\\dot a)\\hat s}{b-a}Q^{-1}\\hat e\\big) + \\hat J(t)\\hat e$ together with the port equations (13) is a Dirac structure for all $t$. From this, Proposition 3 concludes that the moving-boundary system is exactly the dynamical system $(\\partial_t\\hat x, \\hat f_\\partial, Q\\hat x, \\hat e_\\partial) \\in \\hat D(t)$, a genuine generalization of the static-domain definition whose two extra flow terms correspond to compression of the domain and translation of material points through it. The power balance acquires boundary terms equal to the boundary velocity times the energy density at the boundary, and when $\\dot a = \\dot b = 0$ all formulas collapse to the classical Stokes-Dirac structure.","pith_inferences":["A natural next step is a structure-preserving dynamic-mesh discretization: the scheme in Section 4 deliberately drops exact power balance when boundaries move, while Theorem 1 provides the discrete Dirac structure such a scheme would need to respect exactly.","The domain of validity is genuinely fixed by Assumption 1(iii): opposite-sign boundary motion (one end expanding while the other contracts) falls outside the theorem, and a different port assignment or a realification of the structure would be needed to cover it — for instance breathing domains or two-sided Stefan-like fronts moving against each other.","Because the construction is generic in the Hamiltonian operator $J = J_0 + J_1\\partial_s$, the same moving-boundary treatment should apply to transport, shallow-water, and beam-type models, not only the lossless transmission line shown here.","The complex-valued ports suggest that the moving-boundary structure may be viewed as a complexified Dirac structure whose real power pairing is recovered from the boundary-velocity sign condition; making that interpretation explicit could connect the framework to para-complex or phase-field formulations of moving interfaces."],"forward_implications":["When the boundary velocities vanish, Proposition 3 reduces to the standard static-domain port-Hamiltonian definition, so the new class strictly generalizes the old one.","The power balance of a moving-domain system gains the boundary power flow $\\frac{1}{2}(\\dot b\\, \\hat e^\\top(1)\\hat x(1) - \\dot a\\, \\hat e^\\top(0)\\hat x(0))$: boundary motion itself transports energy in or out, beyond the usual port power.","For the lossless transmission line, total charge and flux in the moving segment change only through boundary terms; boundary motion acts like an extra current equal to the charge density times the boundary velocity.","The Dirac structure yields a dynamic-grid discretization of the telegrapher's equations whose mesh tracks the moving boundaries, with numerical error that decreases as the number of elements grows.","The discretization recovers the structure-preserving scheme of [13] whenever the domain is static, and otherwise approximately preserves power balance with bounded error."],"supporting_citations":[{"why":"Supplies Lemma 1, the Stokes-Dirac structure for fixed domains that the paper generalizes.","marker":"[33]"},{"why":"Gives the construction $M = I_n$, $S_1 = -\\tfrac{1}{2}J_1$ for full-rank $J_1$, used for the transmission-line ports.","marker":"[20]"},{"why":"Provides the mixed-element Hamiltonian discretization that Section 4 modifies for dynamic meshing.","marker":"[13]"},{"why":"Justifies Dirac structures over the complex field, needed for the complex boundary ports.","marker":"[18]"},{"why":"The Leibniz integral rule used to differentiate the moving integration bounds.","marker":"[26]"},{"why":"The original boundary port-Hamiltonian formulation that this paper extends.","marker":"[31]"},{"why":"Background reference for the boundary port-Hamiltonian distributed parameter systems framework.","marker":"[11]"}],"fun_headline_variants":["Moving boundaries gain a port-Hamiltonian structure","Time-varying Dirac structure for moving-boundary systems","Dynamic meshing for moving-boundary port-Hamiltonian systems","Port-Hamiltonian form for PDEs with moving boundaries"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole construction assumes the two boundary velocities always have the same sign (both boundaries move in the same direction), because the proof needs $\\sqrt{\\dot a}\\sqrt{\\dot b} = \\sqrt{\\dot a \\dot b}$ to keep the port pairings real.","fun_headline_variants_meta":{"raw":{"variants":["Moving boundaries gain a port-Hamiltonian structure","Time-varying Dirac structure for moving-boundary systems","Dynamic meshing for moving-boundary port-Hamiltonian systems","Port-Hamiltonian form for PDEs with moving boundaries"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000551,"raw_usage":{"total_tokens":2599,"prompt_tokens":888,"completion_tokens":1711,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":504,"completion_tokens_details":{"reasoning_tokens":1647}},"tokens_in":504,"tokens_out":1711,"duration_ms":11968,"temperature":1.0,"reasoning_tokens":1647,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T14:46:55.917318+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Evaluate the Dirac orthogonality condition $\\hat D(t) = \\hat D^\\perp(t)$ symbolically for the ports (13) under $\\dot a(t) < 0 < \\dot b(t)$: the substitution $\\sqrt{\\dot a}\\sqrt{\\dot b} = \\sqrt{\\dot a\\dot b}$ used in Appendix C ceases to hold over the reals, so a direct check of whether the subspace is self-orthogonal under the complex pairing settles whether Theorem 1 survives outside Assumption 1(iii).","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies Lemma 1, the Stokes-Dirac structure for fixed domains that the paper generalizes."},{"cited_title":"Le Gorrec, H","cited_arxiv_id":null,"evidence_quote":"Gives the construction $M = I_n$, $S_1 = -\\tfrac{1}{2}J_1$ for full-rank $J_1$, used for the transmission-line ports."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the mixed-element Hamiltonian discretization that Section 4 modifies for dynamic meshing."},{"cited_title":"Jeltsema and A","cited_arxiv_id":null,"evidence_quote":"Justifies Dirac structures over the complex field, needed for the complex boundary ports."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The Leibniz integral rule used to differentiate the moving integration bounds."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The original boundary port-Hamiltonian formulation that this paper extends."},{"cited_title":"Duindam, A","cited_arxiv_id":null,"evidence_quote":"Background reference for the boundary port-Hamiltonian distributed parameter systems framework."}],"review_version":1}