alphaFeedback

WorksMathematics

Dissipative Hamiltonian representations of linear evolution equations

AI contributions
Not recorded for this paper
View details

Abstract

We investigate explicit and implicit dissipative Hamiltonian representations of generators of strongly continuous semigroups on real or complex Hilbert spaces. For a prescribed positive definite selfadjoint energy operator HH satisfying dom⁡(A)⊂dom⁡(H)\operatorname{dom}(A)\subset\operatorname{dom}(H), we show that the Lyapunov dissipation inequality is equivalent to an exact factorization A=DHA=DH, where DD can be chosen maximally dissipative. For bounded HH, the criterion also yields an implicit representation. We further study energy constructions based on closed intertwiners and output trajectories, distinguishing contractive energy-space dynamics from operator factorizations on the original state space. Examples illustrate the roles of noncoercivity, unboundedness, and domain compatibility.

Introduction

Dissipative Hamiltonian representations of dynamical systems express the evolution of a system in terms of an energy operator and a dissipation operator. Besides making energy balances transparent, these representations provide a structural starting point for investigating well-posedness, stability, and the properties of related operator pencils. For finite-dimensional linear differential equations and differential-algebraic equations, the relationships between various Hamiltonian, dissipative, and relation-based descriptions have been studied extensively; see, for instance, [4, 10, 20, 22, 23]. Infinite-dimensional counterparts arise in the analysis of partial differential equations (PDEs) and abstract port-Hamiltonian systems. The formulation of PDEs as abstract dissipative Hamiltonian differential-algebraic equations on extended state spaces and the spectral theory of the associated operator pencils have been investigated in [21, 32]; see also [27, 28]. Explicit and implicit port-Hamiltonian formulations of particular distributed-parameter models have been compared in [6]. The existing literature approaches dissipation in infinite-dimensional systems from several complementary perspectives. Passive system theory emphasizes input-output behavior and storage functions [3, 29], while port-Hamiltonian formulations describe partial differential equations directly in terms of energy and interconnection structures [14, 27]. We take a complementary viewpoint: starting from a given semigroup generator, we ask whether its dynamics admit a dissipative Hamiltonian representation and which properties the energy operator must satisfy. Related work in the finite-dimensional setting includes [1, 5]; see [2] for preliminary results in infinite dimensions.

More precisely, we consider the linear evolution equation

x˙(t)=Ax(t),(1)\dot{x}(t) = Ax(t), \tag*{(1)}

where A:dom⁡(A)⊂X→XA:\operatorname{dom}(A)\subset\mathcal{X}\to\mathcal{X} generates a strongly continuous semigroup T=(T(t))t≥0T=(T(t))_{t\geq0} on a Hilbert space X\mathcal{X}. Throughout the paper, Hilbert spaces may be real or complex, and all spaces appearing in the same statement are understood to be over the same scalar field. We denote the norm and inner product on X\mathcal{X} by ∥⋅∥X\|\cdot\|_{\mathcal{X}} and ⟨⋅,⋅⟩X\langle\cdot,\cdot\rangle_{\mathcal{X}}, respectively, and omit the subscripts whenever the underlying space is clear. We investigate when its dynamics can be represented in either of the two forms

x˙(t)=DeHx(t),(2)\dot{x}(t)=D_{\mathrm{e}}Hx(t), \tag*{(2)}
Mz˙(t)=Diz(t),(3)M\dot{z}(t)=D_{\mathrm{i}}z(t), \tag*{(3)}

with dissipative operators DeD_{\mathrm{e}} and DiD_{\mathrm{i}} and positive selfadjoint energy operators HH or MM, respectively. The dissipative operators DeD_{\mathrm{e}} and DiD_{\mathrm{i}} need not coincide. Throughout the paper, we assume that ker⁡H=ker⁡M={0}\ker H=\ker M=\{0\}, but do not require either energy operator to be coercive.

The implicit formulation (3) is reminiscent of the mass-operator formulations arising in weak formulations of partial differential equations and their finite element discretizations. We therefore also refer to it as a mass representation.

The case of energy operators with nontrivial kernels is also of interest, but lies outside the scope of the present paper; see [32]. For sufficiently regular solutions, the quadratic energies of the explicit and implicit formulations satisfy, respectively,

ddt12∥H1/2x(t)∥X2=Re⁡⟨DeHx(t),Hx(t)⟩X≤0,\frac{\mathrm{d}}{\mathrm{d}t}\frac{1}{2}\|H^{1/2}x(t)\|_{\mathcal{X}}^2=\operatorname{Re}\langle D_{\mathrm{e}}Hx(t),Hx(t)\rangle_{\mathcal{X}}\leq0,
ddt12∥M1/2z(t)∥X2=Re⁡⟨Diz(t),z(t)⟩X≤0.\frac{\mathrm{d}}{\mathrm{d}t}\frac{1}{2}\|M^{1/2}z(t)\|_{\mathcal{X}}^2=\operatorname{Re}\langle D_{\mathrm{i}}z(t),z(t)\rangle_{\mathcal{X}}\leq0.

In infinite dimensions, however, these formal energy balances do not automatically lead to well-defined operator representations. For the explicit equation (2), the product DeHD_{\mathrm{e}}H is defined only for states x∈dom⁡(H)x\in\operatorname{dom}(H) with Hx∈dom⁡(De)Hx\in\operatorname{dom}(D_{\mathrm{e}}). Thus, recovering the original evolution requires not merely a formal factorization, but equality of the operators AA and DeHD_{\mathrm{e}}H, including their domains.

The implicit formulation (3) involves a different difficulty. If MM is unbounded, then the term Mz˙M\dot{z} is meaningful only when z˙∈dom⁡(M)\dot{z}\in\operatorname{dom}(M), which need not hold for every classical solution of the original evolution equation. Even for bounded and injective MM, the dissipative operator Di=MAD_{\mathrm{i}}=MA need not be closed. If MM is not surjective, then λM−Di\lambda M-D_{\mathrm{i}} is not surjective for any λ\lambda.

A further distinction arises when the energy is defined by a closed quadratic form. Its dissipation inequality then may yield contractive dynamics on a suitable energy space without providing an operator factorization on the original state space X\mathcal{X}. We therefore distinguish exact operator representations, implicit equations with specified solution classes, and representations at the level of quadratic forms. In particular, we do not assume a decomposition of the dissipative operators into skew-adjoint and negative semidefinite selfadjoint parts. For unbounded operators, such a decomposition requires additional domain conditions.

A useful perspective comes from the renorming theory of semigroups. Every uniformly bounded semigroup on a Banach space is contractive with respect to the equivalent norm

∥x∥∗:=sup⁡t≥0∥T(t)x∥.\|x\|_*:=\sup_{t\geq0}\|T(t)x\|.

More generally, every strongly continuous semigroup becomes quasi-contractive after an equivalent Banach-space renorming; see [8]. Even when the original space is a Hilbert space, such a norm need not be induced by an inner product. In fact, there exist exponentially stable Hilbert-space semigroups that are not similar, through a boundedly invertible operator, to a contraction semigroup; see [17, 24]. By contrast, a contractive equivalent Hilbertian norm is induced by a bounded, coercive positive operator. This observation motivates allowing positive energy operators that are non-coercive or unbounded, while retaining a Hilbert-space framework.

A natural way of constructing such energy operators is through intertwining operators. Suppose that a closed, densely defined, injective operator E:dom⁡(E)⊂X→ZE : \operatorname{dom}(E) \subset\mathcal{X} \to\mathcal{Z} intertwines the strongly continuous semigroup TT with a contraction semigroup SS on a Hilbert space Z\mathcal{Z}, in the sense that dom⁡(E)\operatorname{dom}(E) is TT-invariant and

ET(t)x=S(t)Exfor all x∈dom⁡(E),t≥0.(4)ET(t)x = S(t)Ex \qquad\text{for all } x \in\operatorname{dom}(E),\quad t \ge0. \tag*{(4)}

The positive quadratic form h[x]=∥Ex∥Z2h[x] = \|Ex\|_{\mathcal{Z}}^2 is then closed and nonincreasing along trajectories in its domain, and its associated positive selfadjoint operator is H=E∗EH = E^*E. Neither EE nor its inverse on ran⁡E\operatorname{ran} E is required to be bounded. This form-theoretic construction is closely related to existing work on unbounded storage operators and pseudo-similarity, in particular [3]. It does not, without further domain assumptions, provide an operator identity of the form A=DHA = DH or an equivalent implicit equation in the original state space. Understanding this additional step is central to the representation problem that we discuss in this paper.

A particularly useful realization of the intertwiner is furnished by output trajectories. For a suitable observation operator CC with values in a Hilbert space Y\mathcal{Y}, consider

(OCx)(t)=CT(t)xfor all t≥0,(\mathcal{O}_{C}x)(t) = CT(t)x \qquad\text{for all } t \ge0,

where the domain consists of the initial states with square-integrable output trajectories. This mapping intertwines TT with the contractive left-shift semigroup on L2(0,∞;Y)L^2(0,\infty;\mathcal{Y}). The shift approach has a substantial history [18, 31], as do admissibility and observability criteria [9, 11, 13, 26, 30] [12]. Closed, potentially unbounded output-trajectory operators are also treated in [3]. In particular, the existence of an infinite-time admissible and exactly observable output operator is characterized by strong stability together with similarity to a contraction semigroup [24 Theorem 8.11]. Here we use observation maps primarily as a construction tool: we allow one-sided estimates, or neither upper nor lower estimates, and seek the resulting dissipative Hamiltonian representations rather than new admissibility criteria.

Our central result is a fixed-energy representation criterion: under the domain condition dom⁡(A)⊂dom⁡(H)\operatorname{dom}(A) \subset\operatorname{dom}(H), the Lyapunov dissipation inequality is equivalent to an exact factorization A=DHA = DH with a maximally dissipative operator DD. In particular, maximal dissipativity of the factor can always be achieved by extension, without any additional range condition. For bounded HH, the criterion is also equivalent to contractivity of the associated quadratic energy and yields an implicit representation. In particular, every exponentially stable semigroup has both representations with bounded energy. We then use pseudo-similarity and output-trajectory maps to construct energies that need not be bounded or coercive. For unbounded energy operators, a contractive quadratic form always determines a contraction semigroup on a canonical energy space, but a strong factorization on the original space demands additional domain compatibility. A dual construction and a variational interpretation complement the operator result. We stress that the aim of the paper is a unified representation framework, not a new general observability theorem.

The paper is organized as follows. § 2 introduces the representation concepts and pseudo-similarity. § 3 gives the fixed-energy criterion and bounded-energy consequences. § 4 treats closed intertwiners and the output-trajectory construction. § 5 discusses energy-space dynamics, duality, and variational realizations. § 6 provides examples and counterexamples, and § 7 a conclusion.

Preliminaries and representations

Throughout the paper, Hilbert-space inner products are linear in the first argument. All operators are linear, and equalities of operators are understood to include equality of their domains. For a Hilbert space X\mathcal{X}, L(X)\mathcal{L}(\mathcal{X}) denotes the bounded operators on X\mathcal{X}. If AA generates a strongly continuous semigroup TT, then AA is densely defined and closed, and T(t)dom⁡(A)⊂dom⁡(A)T(t)\operatorname{dom}(A) \subset\operatorname{dom}(A) with AT(t)x=T(t)AxAT(t)x = T(t)Ax for x∈dom⁡(A)x \in\operatorname{dom}(A); see [8 Lemma II.1.3 and Theorem II.1.4]. A subspace D⊂X×X\mathcal{D} \subset\mathcal{X} \times\mathcal{X} is called dissipative if

Re⁡⟨u,v⟩X≤0for all (u,v)∈D.\operatorname{Re}\langle u,v\rangle_{\mathcal{X}} \le0 \qquad\text{for all } (u,v) \in\mathcal{D}.

An operator is dissipative if its graph is dissipative, and is maximally dissipative if its graph has no proper dissipative linear extension. In particular, if a dissipative operator DD satisfies ran⁡(λI−D)=X\operatorname{ran}(\lambda I-D)=\mathcal{X} for some λ>0\lambda>0, then it is maximally dissipative by the injectivity property for dissipative extensions [8 Proposition II.3.14(i)]. By the Lumer–Phillips theorem [8 Corollary II.3.20], it also generates a strongly continuous contraction semigroup; see also [19]. More general linear relations are useful for singular Hamiltonian descriptions [10].

An operator HH will be called positive definite if it is selfadjoint, H≥0H \ge0, and ker⁡H={0}\operatorname{ker} H = \{0\}; note that this does not mean that H≥cIH \ge cI for some c>0c>0. The latter property will be called coercivity. If HH is unbounded, its positive quadratic form has domain dom⁡(H1/2)\operatorname{dom}(H^{1/2}).

Definition 2.1. Let AA be the generator of a strongly continuous semigroup TT on a Hilbert space X\mathcal{X}.

(i) An explicit dissipative Hamiltonian representation of AA is a factorization A=DHA=DH, where HH is a positive definite selfadjoint operator and DD is dissipative.

(ii) An implicit dissipative Hamiltonian representation of (1) consists of a positive definite selfadjoint operator MM and a dissipative operator DD such that the equation Mz˙=DzM\dot{z}=Dz is meaningful and equivalent to (1) for a specified class of solutions. If MM is bounded, a sufficient condition for equivalence of classical solutions is D=MAD=MA with dom⁡(D)=dom⁡(A)\operatorname{dom}(D)=\operatorname{dom}(A).

Note that the specification of the solution class in (ii) is essential. In particular, writing D=MAD=MA when MM is unbounded does not say that MAMA is defined on dom⁡(A)\operatorname{dom}(A). The same care is required for the product DHDH.

Definition 2.2. Let TjT_j be strongly continuous semigroups on Hilbert spaces Xj\mathcal{X}_j for j=1,2j=1,2. A pseudo-similarity from T1T_1 to T2T_2 is a closed, densely defined, injective operator E:dom⁡(E)⊂X1→X2E:\operatorname{dom}(E)\subset\mathcal{X}_1\to\mathcal{X}_2 such that dom⁡(E)\operatorname{dom}(E) is T1T_1-invariant and

ET1(t)x=T2(t)Exfor all x∈dom⁡(E),t≥0.ET_1(t)x=T_2(t)Ex \qquad\text{for all } x\in\operatorname{dom}(E),\quad t\ge0.

We may, and often will, replace X2\mathcal{X}_2 by ran⁡E‾\overline{\operatorname{ran} E}, which is invariant under T2T_2.

The word pseudo-similarity refers here to a directed intertwining relation; no boundedness of EE or E−1 ⁣:ran⁡E→X1E^{-1}\colon\operatorname{ran} E \to\mathcal{X}_1 is assumed. Various concepts of quasi- and pseudo-similarity are studied in operator theory; see, e.g., [3 Definition 1.8] and [24].

A fixed-energy representation criterion

We begin by studying the operator-theoretic core of the representation problem. Although the Hilbert-space Lyapunov inequality is classical, see e.g. [8, 9, 30], the purpose of the following construction is to identify the corresponding dissipative Hamiltonian factors, including their domains, within a single statement.

Theorem 3.1 (Fixed-energy representation criterion). Let AA be the generator of a strongly continuous semigroup TT on a Hilbert space X\mathcal{X}, and let HH be a positive definite selfadjoint operator on X\mathcal{X}. Suppose that

dom⁡(A)⊂dom⁡(H).(5)\operatorname{dom}(A) \subset\operatorname{dom}(H). \tag*{(5)}

Then the following assertions are equivalent:

(a) Re⁡⟨Ax,Hx⟩≤0\operatorname{Re}\langle Ax,Hx\rangle\le0 for every x∈dom⁡(A)x \in\operatorname{dom}(A);

(b) there exists a dissipative operator DD on X\mathcal{X} such that A=DHA = DH, including equality of domains;

(c) there exists a maximally dissipative operator DD on X\mathcal{X} such that A=DHA = DH, including equality of domains.

If these conditions hold, the canonical dissipative factor

dom⁡(De)=Hdom⁡(A),De(Hx)=Axfor all x∈dom⁡(A)(6)\operatorname{dom}(D_{\mathrm{e}})=H\operatorname{dom}(A), \qquad D_{\mathrm{e}}(Hx)=Ax \quad\text{for all }x\in\operatorname{dom}(A) \tag*{(6)}

is densely defined. Every dissipative extension of DeD_{\mathrm{e}} preserves the factorization A=DHA = DH. In particular, DeD_{\mathrm{e}} admits a maximally dissipative extension.

If HH is bounded, then the above assertions are also equivalent to

T(t)∗HT(t)≤Hfor all t≥0.(7)T(t)^*HT(t) \le H \qquad\text{for all }t\ge0. \tag*{(7)}

In this case the operator Di=HAD_{\mathrm{i}}=HA with dom⁡(Di)=dom⁡(A)\operatorname{dom}(D_{\mathrm{i}})=\operatorname{dom}(A) is dissipative, and the two equations

x˙=DeHx,Hz˙=Diz(8)\dot{x}=D_{\mathrm{e}}Hx,\qquad H\dot{z}=D_{\mathrm{i}}z \tag*{(8)}

both describe the dynamics generated by AA for classical solutions. Moreover, E=H1/2E=H^{1/2} intertwines TT with a contraction semigroup on ran⁡E=X\operatorname{ran} E=\mathcal{X}.

Proof. Assume (a). Injectivity of HH makes (3.2) well-defined, and

Re⁡⟨DeHx,Hx⟩=Re⁡⟨Ax,Hx⟩≤0for all x∈dom⁡(A).\operatorname{Re}\langle D_{\mathrm{e}}Hx,Hx\rangle=\operatorname{Re}\langle Ax,Hx\rangle\le0 \qquad\text{for all }x\in\operatorname{dom}(A).

Thus, DeD_{\mathrm{e}} is dissipative. Moreover, DeH=AD_{\mathrm{e}}H=A, including equality of domains, by the injectivity of HH.

We first show that DeD_{\mathrm{e}} is densely defined. Choose μ>0\mu>0 such that μI−A\mu I-A is bijective. Let z∈(Hdom⁡(A))⊥z\in(H\operatorname{dom}(A))^\perp and set x=(μI−A)−1zx=(\mu I-A)^{-1}z. Then

0=Re⁡⟨z,Hx⟩=μ⟨x,Hx⟩−Re⁡⟨Ax,Hx⟩≥μ⟨x,Hx⟩.0=\operatorname{Re}\langle z,Hx\rangle=\mu\langle x,Hx\rangle-\operatorname{Re}\langle Ax,Hx\rangle\ge\mu\langle x,Hx\rangle.

Positive definiteness of HH implies that x=0x=0, and hence z=0z=0. Therefore H‾ dom⁡(A)=X\overline{H}\,\operatorname{dom}(A)=\mathcal{X}.

Since DeD_e is densely defined and dissipative, its closure D0=De‾D_0=\overline{D_e} is dissipative. The range R=ran⁡(I−D0)\mathcal{R}=\operatorname{ran}(I-D_0) is closed, since D0D_0 is closed and

∥(I−D0)u∥≥∥u∥for all u∈dom⁡(D0).\lVert(I-D_0)u\rVert\ge\lVert u\rVert\qquad\text{for all }u\in\operatorname{dom}(D_0).

Define an extension DD by

dom⁡(D)=dom⁡(D0)∔R⊥,D(u+v)=D0u−v,\operatorname{dom}(D)=\operatorname{dom}(D_0)\mathbin{\dotplus}\mathcal{R}^{\perp}, \qquad D(u+v)=D_0u-v,

where u∈dom⁡(D0)u\in\operatorname{dom}(D_0) and v∈R⊥v\in\mathcal{R}^{\perp}. The sum is direct by the dissipativity of D0D_0. Since (I−D0)u⊥v(I-D_0)u\perp v, we obtain

Re⁡⟨D(u+v),u+v⟩=Re⁡⟨D0u,u⟩−∥v∥2≤0.\operatorname{Re}\langle D(u+v),u+v\rangle=\operatorname{Re}\langle D_0u,u\rangle-\lVert v\rVert^2\le0.

Furthermore,

ran⁡(I−D)=R⊕R⊥=X.\operatorname{ran}(I-D)=\mathcal{R}\oplus\mathcal{R}^{\perp}=\mathcal{X}.

Thus DD is maximally dissipative and generates a strongly continuous contraction semigroup. It remains to show that every dissipative extension D~\widetilde{D} of DeD_e satisfies A=D~HA=\widetilde{D}H. Since D~\widetilde{D} extends DeD_e, we have A⊆D~HA\subseteq\widetilde{D}H. Let x∈dom⁡(D~H)x\in\operatorname{dom}(\widetilde{D}H). By the surjectivity of μI−A\mu I-A, there exists y∈dom⁡(A)y\in\operatorname{dom}(A) such that

(μI−A)y=(μI−D~H)x.(\mu I-A)y=(\mu I-\widetilde{D}H)x.

With w=x−yw=x-y we obtain w∈dom⁡(D~H)w\in\operatorname{dom}(\widetilde{D}H) and D~Hw=μw\widetilde{D}Hw=\mu w. The dissipativity then yields that

0≥Re⁡⟨D~Hw,Hw⟩=μ⟨w,Hw⟩,0\ge\operatorname{Re}\langle\widetilde{D}Hw,Hw\rangle=\mu\langle w,Hw\rangle,

and the positive definiteness of HH implies that w=0w=0. Consequently, x=y∈dom⁡(A)x=y\in\operatorname{dom}(A), which proves that A=D~HA=\widetilde{D}H.

We have shown that (a) implies (c). The implication (c) ⇒\Rightarrow (b) is immediate, and (b) implies (a) by dissipativity.

Finally, suppose that HH is bounded. Differentiating ⟨HT(t)x,T(t)x⟩\langle HT(t)x,T(t)x\rangle for x∈dom⁡(A)x\in\operatorname{dom}(A), and using density, proves the equivalence of (a) and (3.3). The identity Di=HAD_i=HA is dissipative, and injectivity of HH gives the stated classical equivalence. On ran⁡H1/2\operatorname{ran}H^{1/2} define S(t)H1/2x=H1/2T(t)xS(t)H^{1/2}x=H^{1/2}T(t)x. The inequality (3.3) then makes S(t)S(t) contractive. It extends by density to a strongly continuous contraction semigroup; strong continuity follows from the boundedness of H1/2H^{1/2} and the strong continuity of TT.

Remark 3.2 (Maximality of the canonical factor). For every λ>0\lambda>0 one has

ran⁡(λI−De)=ran⁡(λH−A).\operatorname{ran}(\lambda I-D_e)=\operatorname{ran}(\lambda H-A).

Therefore, the canonical factor DeD_e is maximally dissipative if and only if ran⁡(λH−A)=X\operatorname{ran}(\lambda H-A)=\mathcal{X} for some λ>0\lambda>0.

If AA is bijective, then HA−1HA^{-1} is closed and everywhere defined, hence bounded. Therefore

λH−A=(λHA−1−I)A\lambda H-A=(\lambda HA^{-1}-I)A

is bijective for all sufficiently small λ>0\lambda>0. In this case DeD_e itself is maximally dissipative.

The range condition is not automatic, even when both AA and HH are bounded. The second construction in Subsection 6.2 shows that the canonical factor may have a proper dense, nonclosed domain, while its maximally dissipative extension still yields the same exact factorization.

Our next result discusses the implicit formulation.

Proposition 3.3 (The mass-range resolvent). Let M∈L(X)M \in\mathcal{L}(\mathcal{X}) be selfadjoint, positive definite, and suppose that Re⁡⟨MAx,x⟩≤0\operatorname{Re}\langle MAx,x\rangle\le0 for x∈dom⁡(A)x \in\operatorname{dom}(A). Set Di=MAD_{\mathrm{i}}=MA with dom⁡(Di)=dom⁡(A)\operatorname{dom}(D_{\mathrm{i}})=\operatorname{dom}(A). Then DiD_{\mathrm{i}} is dissipative, and for every scalar λ\lambda such that λI−A\lambda I-A is bijective, the pencil

L(λ)=λM−Di:dom⁡(A)⟶XL(\lambda)=\lambda M-D_{\mathrm{i}}:\operatorname{dom}(A)\longrightarrow\mathcal{X}

is injective, with ran⁡L(λ)=ran⁡M\operatorname{ran}L(\lambda)=\operatorname{ran}M. Thus L(λ)L(\lambda) is bijective as a map from dom⁡(A)\operatorname{dom}(A) onto ran⁡M\operatorname{ran}M, and

L(λ)−1M=(λI−A)−1on X.(9)L(\lambda)^{-1}M=(\lambda I-A)^{-1}\quad\text{on }\mathcal{X}. \tag*{(9)}

In particular, the implicit equation (3) describes the dynamics generated by AA for classical solutions, although, when MM is not surjective, no L(λ)L(\lambda) is surjective onto X\mathcal{X}.

Proof. The assumed inequality Re⁡⟨MAx,x⟩≤0\operatorname{Re}\langle MAx,x\rangle\le0 for x∈dom⁡(A)x \in\operatorname{dom}(A) shows that DiD_{\mathrm{i}} is dissipative. Fix λ\lambda such that λI−A\lambda I-A is bijective. Then L(λ)x=M(λI−A)xL(\lambda)x=M(\lambda I-A)x for every x∈dom⁡(A)x \in\operatorname{dom}(A). Since MM is injective and λI−A\lambda I-A is bijective, this identity proves injectivity, the range formula, and (9). Classical equivalence follows from the injectivity of MM.

Remark 3.4. The implicit factor Di=HAD_{\mathrm{i}}=HA need not be closed, let alone maximally dissipative, even though the two equations in (8) are classically equivalent. Theorem 3.3 shows why invertibility of the pencil L(λ)L(\lambda) on all of X\mathcal{X} is not the right criterion for the dynamics in this setting. For unbounded HH, condition (5) is an actual restriction, not a technical convenience. The first construction in Subsection 6.2 exhibits a contractive closed energy form whose associated energy operator does not satisfy this domain condition. If the condition holds, form contractivity implies (a) in Theorem 3.1; see § 5.

As a corollary we obtain that every exponentially stable semigroup admits both dissipative Hamiltonian representations.

Corollary 3.5. Every exponentially stable strongly continuous semigroup on a Hilbert space admits both an explicit and an implicit dissipative Hamiltonian representation with the same bounded, positive definite energy operator HH. For this choice of HH, the canonical explicit factor DeD_{\mathrm{e}} is maximally dissipative.

Proof. For the infinite-time Gramian, understood as a weak operator integral,

H=∫0∞T(t)∗T(t) dt,(10)H=\int_{0}^{\infty}T(t)^{*}T(t)\,\mathrm{d}t, \tag*{(10)}

exponential stability implies boundedness, since ∥T(t)∥≤Ke−ωt\lVert T(t)\rVert\le Ke^{-\omega t} gives ∥H∥≤K2/(2ω)\lVert H\rVert\le K^{2}/(2\omega), and ⟨Hx,x⟩=∫0∞∥T(t)x∥2 dt>0\langle Hx,x\rangle=\int_{0}^{\infty}\lVert T(t)x\rVert^{2}\,\mathrm{d}t>0 for x≠0x\ne0. Differentiating the energy for x∈dom⁡(A)x\in\operatorname{dom}(A) gives 2Re⁡⟨HAx,x⟩=−∥x∥22\operatorname{Re}\langle HAx,x\rangle=-\lVert x\rVert^{2}. The integral representation in [8 Theorem II.1.10] yields A−1=−∫0∞T(t) dtA^{-1}=-\int_{0}^{\infty}T(t)\,\mathrm{d}t. Hence Theorem 3.1 and Theorem 3.2 apply.

The construction of the Gramian (3.6) and its connection to Lyapunov equations are described in [30 Theorem 5.1.1]. Conversely, if the identity observation has square-integrable trajectories for all x∈Xx \in\mathcal{X}, the Datko–Pazy theorem yields exponential stability; see [7] and [8 Theorem V.1.8].

Remark 3.6 (Lyapunov equations and dissipative Hamiltonian representations). For bounded HH, the dissipation condition in Theorem 3.1 is precisely the Lyapunov inequality A∗H+HA≤0A^{*}H + HA \le0, understood in the quadratic-form sense on dom⁡(A)\operatorname{dom}(A). For the Gramian (3.6) of an exponentially stable semigroup, the following stronger operator identity holds: Hdom⁡(A)⊂dom⁡(A∗)H\operatorname{dom}(A) \subset\operatorname{dom}(A^{*}) and

(A∗H+HA)x=−xfor all x∈dom⁡(A).(A^{*}H + HA)x = -x \qquad\text{for all } x \in\operatorname{dom}(A).

Indeed, polarization of the Gramian identity yields

⟨HAx,y⟩+⟨Hx,Ay⟩=−⟨x,y⟩for all x,y∈dom⁡(A),\langle HAx,y\rangle+ \langle Hx,Ay\rangle= -\langle x,y\rangle\qquad\text{for all } x,y \in\operatorname{dom}(A),

which proves the assertion.

In finite dimensions, positive definiteness of HH implies coercivity, and hence the existence of an equivalent contractive Hilbertian norm. In infinite dimensions, however, even the strict Lyapunov equation need not yield a coercive solution. Thus exponential stability does not generally imply similarity to a contraction semigroup, as illustrated by the examples in § 6.

The significance of Theorem 3.1 is that even a non-coercive Lyapunov operator can yield an exact factorization A=DHA = DH with a maximally dissipative factor DD, while the canonical factor DeD_{\mathrm{e}} need not itself be maximally dissipative; see Theorem 3.2.

Pseudo-similarity and output trajectories

Closed energy forms from intertwiners

Theorem 4.1. Let X\mathcal{X} and Z\mathcal{Z} be Hilbert spaces, let E:dom⁡(E)⊂X→ZE : \operatorname{dom}(E) \subset\mathcal{X} \to\mathcal{Z} be a closed, densely defined, injective operator, and let SS be a contraction semigroup on Z\mathcal{Z}. Assume that EE intertwines TT with SS as in (1.4). Then

h(x,y)=⟨Ex,Ey⟩Z,dom⁡(h)=dom⁡(E),(11)\mathfrak{h}(x,y) = \langle Ex,Ey\rangle_{\mathcal{Z}}, \qquad\operatorname{dom}(\mathfrak{h}) = \operatorname{dom}(E), \tag*{(11)}

is a densely defined, closed, positive definite form, satisfying

h[T(t)x]≤h[x]for all x∈dom⁡(E),t≥0.(12)\mathfrak{h}[T(t)x] \le\mathfrak{h}[x] \qquad\text{for all } x \in\operatorname{dom}(E),\quad t \ge0. \tag*{(12)}

The associated positive selfadjoint operator is H=E∗EH = E^{*}E, with

dom⁡(H1/2)=dom⁡(E),∥H1/2x∥=∥Ex∥.\operatorname{dom}(H^{1/2}) = \operatorname{dom}(E), \qquad\lVert H^{1/2}x\rVert= \lVert Ex\rVert.

In particular, the semigroup TT is contractive on the completion of dom⁡(E)\operatorname{dom}(E) in the h1/2\mathfrak{h}^{1/2}-norm, which can be identified with ran⁡E\operatorname{ran} E.

Proof. The form in (4.1) is closed because EE is closed; it is strictly positive because EE is injective. The form representation theorem for positive selfadjoint operators yields that H=E∗EH = E^*E and the asserted square-root identity; see [15 Theorems VI.2.1 and VI.2.23]. By intertwining and the contractivity of SS, we have

∥ET(t)x∥=∥S(t)Ex∥≤∥Ex∥.\lVert ET(t)x\rVert= \lVert S(t)Ex\rVert\le\lVert Ex\rVert.

Finally the map x↦Exx \mapsto Ex identifies the completion in the energy norm with ran⁡E‾\overline{\operatorname{ran} E}, and the induced semigroup is precisely SS restricted to that closed invariant subspace.

The converse to the construction in Theorem 4.1 is also useful: a contractive closed form produces an intertwining operator without an output being specified beforehand.

Proposition 4.2. Let hh be a densely defined, closed, positive definite, nonnegative quadratic form on a Hilbert space X\mathcal{X}, whose domain is TT-invariant and satisfies (4.2). Let HH be its associated positive selfadjoint operator. Then E=H1/2E = H^{1/2} is a closed, densely defined, injective pseudo-similarity from TT to a contraction semigroup on ran⁡H1/2=X\operatorname{ran} H^{1/2} = \mathcal{X}.

Proof. The map S(t)S(t) defined initially by S(t)H1/2x=H1/2T(t)xS(t)H^{1/2}x = H^{1/2}T(t)x is well-defined and contractive on the dense subspace ran⁡H1/2\operatorname{ran} H^{1/2}; it extends uniquely to a contraction semigroup on X\mathcal{X}. It remains to check strong continuity. Fix x∈dom⁡(H1/2)x \in\operatorname{dom}(H^{1/2}) and let tn↓0t_n \downarrow0. The sequence H1/2T(tn)xH^{1/2}T(t_n)x is bounded by the form inequality (4.2). Every weakly convergent subsequence has weak limit H1/2xH^{1/2}x: indeed T(tn)x→xT(t_n)x \to x in X\mathcal{X} and the graph of the closed operator H1/2H^{1/2} is weakly closed. Moreover,

∥H1/2T(tn)x∥≤∥H1/2x∥.\lVert H^{1/2}T(t_n)x\rVert\le\lVert H^{1/2}x\rVert.

Weak lower semicontinuity forces convergence of these norms, and hence strong convergence to H1/2xH^{1/2}x. Density and contractivity give strong continuity on X\mathcal{X}.

The analysis in this subsection shows why closed energy forms provide a natural language for pseudo-similar contraction dynamics. However, they do not, by themselves, produce an everywhere valid operator product DHDH. We discuss this issue in § 5.

The output-trajectory transformation

An observation operator C:dom⁡(A)→YC : \operatorname{dom}(A) \to\mathcal{Y} that maps to a Hilbert space Y\mathcal{Y} is assumed to be bounded with respect to the graph norm ∥x∥dom⁡(A)2=∥x∥X2+∥Ax∥X2\lVert x\rVert_{\operatorname{dom}(A)}^2 = \lVert x\rVert_{\mathcal{X}}^2 + \lVert Ax\rVert_{\mathcal{X}}^2. It is called finite-time admissible if there exist τ>0\tau> 0 and βτ>0\beta_\tau> 0 such that

∫0τ∥CT(t)x∥Y2 dt≤βτ∥x∥X2for all x∈dom⁡(A).\int_0^\tau\lVert CT(t)x\rVert_{\mathcal{Y}}^2\,\mathrm{d}t \le\beta_\tau\lVert x\rVert_{\mathcal{X}}^2 \qquad\text{for all } x \in\operatorname{dom}(A).

It is infinite-time admissible if the same estimate holds with τ=∞\tau= \infty and a finite constant βτ\beta_\tau. An infinite-time admissible observation operator CC is called infinite-time exactly observable if there exists α>0\alpha> 0 such that

∫0∞∥CT(t)x∥Y2 dt≥α∥x∥X2for all x∈dom⁡(A).\int_0^\infty\lVert CT(t)x\rVert_{\mathcal{Y}}^2\,\mathrm{d}t \ge\alpha\lVert x\rVert_{\mathcal{X}}^2 \qquad\text{for all } x \in\operatorname{dom}(A).

See [26, 30] for these notions and related results.

We first consider C∈L(X,Y)C \in\mathcal{L}(\mathcal{X},\mathcal{Y}), for which finite-time admissibility is automatic. Define its possibly unbounded infinite-time output operator by

dom⁡(OC)={x∈X:∫0∞∥CT(t)x∥Y2 dt<∞},\operatorname{dom}(\mathcal{O}_C) = \left\{x \in\mathcal{X} : \int_0^\infty\lVert CT(t)x\rVert_{\mathcal{Y}}^2\,\mathrm{d}t < \infty\right\},
(OCx)(t)=CT(t)xfor all t≥0.(\mathcal{O}_C x)(t) = CT(t)x \qquad\text{for all } t \ge0.

The left-shift semigroup on L2(0,∞;Y)L^2(0,\infty;\mathcal{Y}) is

(S(t)f)(τ)=f(t+τ)for all t≥0 and a.e. τ≥0.(S(t)f)(\tau)=f(t+\tau)\qquad\text{for all }t\geq0\text{ and a.e. }\tau\geq0.

It is a strongly continuous contraction semigroup and it is strongly stable. The next proposition records elementary properties which also motivate more general output maps for graph-bounded CC; see [18], [3 Proposition 3.1], and [30 Secs. 4.3 and 4.6].

Proposition 4.3. For C∈L(X,Y)C\in\mathcal{L}(\mathcal{X},\mathcal{Y}), OC\mathcal{O}_C is closed, its domain is TT-invariant, and

OCT(t)x=S(t)OCx,for all x∈dom⁡(OC), t≥0.\mathcal{O}_CT(t)x=S(t)\mathcal{O}_Cx,\qquad\text{for all }x\in\operatorname{dom}(\mathcal{O}_C),\ t\geq0.

If, in addition, dom⁡(OC)\operatorname{dom}(\mathcal{O}_C) is dense and CT(⋅)x=0CT(\cdot)x=0 implies x=0x=0 for x∈dom⁡(OC)x\in\operatorname{dom}(\mathcal{O}_C), then OC\mathcal{O}_C is injective (the qualitative or approximate observability condition) and defines a pseudo-similarity to the shift semigroup. Moreover, every x∈dom⁡(OC)x\in\operatorname{dom}(\mathcal{O}_C) satisfies

∥OCT(t)x∥L22−∥OCx∥L22=−∫0t∥CT(τ)x∥Y2 dτfor all t≥0.\lVert\mathcal{O}_CT(t)x\rVert_{L^2}^2-\lVert\mathcal{O}_Cx\rVert_{L^2}^2=-\int_0^t\lVert CT(\tau)x\rVert_{\mathcal{Y}}^2\,\mathrm{d}\tau\qquad\text{for all }t\geq0.

Proof. Suppose that xn∈dom⁡(OC)x_n\in\operatorname{dom}(\mathcal{O}_C), xn→xx_n\to x in X\mathcal{X}, and OCxn→f\mathcal{O}_Cx_n\to f in L2(0,∞;Y)L^2(0,\infty;\mathcal{Y}). On every finite interval [0,τ][0,\tau], one has the convergence CT(⋅)xn→CT(⋅)xCT(\cdot)x_n\to CT(\cdot)x uniformly. Hence, f=CT(⋅)xf=CT(\cdot)x almost everywhere and x∈dom⁡(OC)x\in\operatorname{dom}(\mathcal{O}_C), proving closedness. The intertwining identity follows from the semigroup property, and invariance follows since the tail of an L2L^2-function is in L2L^2. The last identity is obtained by subtracting the integrals over [t,∞)[t,\infty) and [0,∞)[0,\infty).

The upper admissibility estimate and the lower observability estimate have a particularly simple interpretation in terms of OC\mathcal{O}_C. Assume that OC\mathcal{O}_C is densely defined and injective. Then

OC bounded⟺∫0∞∥CT(t)x∥Y2 dt≤β∥x∥2,\mathcal{O}_C\text{ bounded}\quad\Longleftrightarrow\quad\int_0^\infty\lVert CT(t)x\rVert_{\mathcal{Y}}^2\,\mathrm{d}t\leq\beta\lVert x\rVert^2,
OC−1 bounded on ran⁡OC⟺∫0∞∥CT(t)x∥Y2 dt≥α∥x∥2,\mathcal{O}_C^{-1}\text{ bounded on }\operatorname{ran}\mathcal{O}_C\quad\Longleftrightarrow\quad\int_0^\infty\lVert CT(t)x\rVert_{\mathcal{Y}}^2\,\mathrm{d}t\geq\alpha\lVert x\rVert^2,

where the respective estimates hold for some α,β>0\alpha,\beta>0 and all x∈dom⁡(OC)x\in\operatorname{dom}(\mathcal{O}_C). The upper estimate is equivalent to infinite-time admissibility of CC: by closedness and density it implies dom⁡(OC)=X\operatorname{dom}(\mathcal{O}_C)=\mathcal{X}. Under this assumption, the lower estimate is precisely the infinite-time exact-observability inequality for (A,C)(A,C).

Without admissibility, the lower estimate may still hold on the possibly proper domain dom⁡(OC)\operatorname{dom}(\mathcal{O}_C). This is a lower observability estimate on the natural domain of the output transformation, but does not imply infinite-time admissibility. In this case ran⁡OC\operatorname{ran}\mathcal{O}_C is closed. Thus, bounded invertibility on the range of the output transformation does not require admissibility of the underlying observation operator.

When both estimates hold, the associated semigroup is similar to a contraction semigroup. A precise existence characterization, including the required strong stability, is given in [24 Theorem 8.11]; compare the earlier results in [3, 11]. We quote these results rather than reproving their observation-theoretic content. The connection can also be formulated intrinsically in terms of Lyapunov operators. A bounded injective output map induces a bounded positive definite H=OC∗OCH = O_C^*O_C such that T(t)∗HT(t)≤HT(t)^*HT(t) \le H and ∥H1/2T(t)x∥→0\|H^{1/2}T(t)x\| \to0 as t→∞t \to\infty. Conversely, given such an HH, the nonnegative form

(x,y)⟼−⟨HAx,y⟩−⟨Hx,Ay⟩(x,y) \longmapsto-\langle HAx,y\rangle- \langle Hx,Ay\rangle

for all x,y∈dom⁡(A)x,y \in\operatorname{dom}(A), is bounded with respect to the graph norm on dom⁡(A)\operatorname{dom}(A). Factoring this nonnegative form yields a graph-norm bounded observation operator C:dom⁡(A)→YC : \operatorname{dom}(A) \to\mathcal{Y}, for a suitable Hilbert space Y\mathcal{Y}, such that

∫0∞∥CT(t)x∥2 dt=⟨Hx,x⟩\int_0^\infty\|CT(t)x\|^2\,dt = \langle Hx,x\rangle

for all x∈dom⁡(A)x \in\operatorname{dom}(A). In particular, CC is infinite-time admissible. If HH is additionally coercive, the pair (A,C)(A,C) is infinite-time exactly observable. These standard Lyapunov and admissibility constructions are discussed in [9, 24]; see also [30 Sec. 5.1].

Both one-sided output constructions for stable semigroups

There is a particularly simple way to construct the two one-sided cases without asserting any similarity to a contraction semigroup.

Theorem 4.4. Let TT be exponentially stable. Then:

(a) The identity observation C0=IC_0 = I produces a bounded, injective infinite-time output map

O0x=T(⋅)x.O_0x = T(\cdot)x.

(b) The graph-norm bounded observation C1=AC_1 = A, initially defined by

O1x=T(⋅)AxO_1x = T(\cdot)Ax

for all x∈dom⁡(A)x \in\operatorname{dom}(A), is closable as an operator from X\mathcal{X} into L2(0,∞;X)L^2(0,\infty;\mathcal{X}). Its closure O‾1\overline{O}_1 is densely defined, closed, injective, intertwines TT with the left shift, and has a bounded inverse on its range.

Consequently, every exponentially stable semigroup has both a bounded injective output transformation and a closed, bounded-below output transformation, even if it is not similar to a contraction semigroup.

Proof. Choose K≥1K \ge1 and ω>0\omega> 0 such that ∥T(t)∥≤Ke−ωt\|T(t)\| \le K e^{-\omega t} for all t≥0t \ge0. Part (a) follows immediately from

∫0∞∥T(t)x∥2 dt≤K22ω∥x∥2.\int_0^\infty\|T(t)x\|^2\,dt \le\frac{K^2}{2\omega}\|x\|^2.

For x∈dom⁡(A)x \in\operatorname{dom}(A), commutation with the semigroup and exponential stability show O1x∈L2O_1x \in L^2 and

∥O1x∥L2≤K2ω∥Ax∥.\|O_1x\|_{L^2} \le\frac{K}{\sqrt{2\omega}}\|Ax\|.

Moreover, differentiating the squared norm, integrating over [0,∞)[0,\infty), and applying the Cauchy–Schwarz inequality yield

∥x∥2=−2Re⁡∫0∞⟨AT(t)x,T(t)x⟩ dt\|x\|^2 = -2\operatorname{Re}\int_0^\infty\langle AT(t)x,T(t)x\rangle\,dt
≤2∥O1x∥L2∥O0x∥L2≤2∥O0∥∥O1x∥L2∥x∥.\le2\|O_1x\|_{L^2}\|O_0x\|_{L^2} \le2\|O_0\|\|O_1x\|_{L^2}\|x\|.

It follows that

∥O1x∥L2≥12∥O0∥∥x∥for all x∈dom⁡(A).(13)\lVert\mathcal{O}_1x\rVert_{L^2} \ge\frac{1}{2\lVert\mathcal{O}_0\rVert}\lVert x\rVert\qquad\text{for all } x \in\operatorname{dom}(A). \tag*{(13)}

To see closability, let BB be the generator of the left-shift semigroup, Bf=f′Bf=f' on H1(0,∞;X)H^1(0,\infty;\mathcal{X}). For x∈dom⁡(A)x\in\operatorname{dom}(A) we have O0x∈dom⁡(B)\mathcal{O}_0x\in\operatorname{dom}(B) and BO0x=O1xB\mathcal{O}_0x=\mathcal{O}_1x. Since O0\mathcal{O}_0 is bounded and BB is closed, BO0B\mathcal{O}_0 is closed. As O1\mathcal{O}_1 is its restriction to dom⁡(A)\operatorname{dom}(A), it is closable. The estimate (13) extends to its closure, proving injectivity and boundedness of the inverse. The intertwining identity holds on dom⁡(A)\operatorname{dom}(A) and passes to the closure using boundedness of T(t)T(t) and S(t)S(t).

Although C1=A:dom⁡(A)→XC_1=A:\operatorname{dom}(A)\to\mathcal{X} is graph-norm bounded, this does not ensure finite-time or infinite-time admissibility. For exponentially stable semigroups, these two notions of admissibility are equivalent; see [30 Sec. 4.6]. Its closed output transformation O‾1\overline{\mathcal{O}}_1 is bounded below. If it were also bounded, the associated semigroup would be similar to a contraction semigroup. Hence, for an exponentially stable semigroup not similar to a contraction, O‾1\overline{\mathcal{O}}_1 must be unbounded, while O‾0−1\overline{\mathcal{O}}_0^{-1} is unbounded on its range.

The existence of a closed, bounded-below pseudo-similarity in Theorem 4.4 cannot be extended to arbitrary uniformly bounded strongly continuous semigroups. Indeed, there exist semigroups TT for which no bounded, injective, dense-range intertwiner FF from a strongly continuous contraction semigroup SS satisfies T(t)F=FS(t)T(t)F=FS(t) for all t≥0t\ge0. The norm bound of TT can even be chosen arbitrarily close to one [25 Theorem 1.4 and Remark 6.3]. In particular, such a semigroup cannot admit a closed, densely defined, bounded-below pseudo-similarity to a contraction semigroup: its inverse on the closed range would provide precisely such an intertwiner.

Unbounded energies and operator representations

Let E:dom⁡(E)⊂X→ZE:\operatorname{dom}(E)\subset\mathcal{X}\to\mathcal{Z} be a closed, densely defined, injective intertwiner from TT to a contraction semigroup SS. Write

V=dom⁡(E),∥x∥V2=∥x∥X2+∥Ex∥Z2,K=ran⁡E‾,H=E∗E.\mathcal{V}=\operatorname{dom}(E), \qquad\lVert x\rVert_{\mathcal{V}}^2=\lVert x\rVert_{\mathcal{X}}^2+\lVert Ex\rVert_{\mathcal{Z}}^2, \qquad\mathcal{K}=\operatorname{ran}\overline{E}, \qquad H=E^*E.

The transition to the energy space draws on closed-form theory [15 Chap. VI, Sec. 2] and the theory of pseudo-similarity and storage functions [29 Sec. 9.2]; see also [3 Secs. 1 and 7]. We record it to specify which operator is maximally dissipative.

Proposition 5.1 (Energy-space realization and variational equation). The restricted semigroup TV=T∣VT_{\mathcal{V}}=T|_{\mathcal{V}} is strongly continuous on V\mathcal{V}, with generator AVA_{\mathcal{V}}, the part of AA in V\mathcal{V}:

dom⁡(AV)={x∈dom⁡(A)∩V:Ax∈V},AVx=Ax.\operatorname{dom}(A_{\mathcal{V}})=\{x\in\operatorname{dom}(A)\cap\mathcal{V}:Ax\in\mathcal{V}\}, \qquad A_{\mathcal{V}}x=Ax.

If BB generates S∣KS|_{\mathcal{K}}, then BB is maximally dissipative on K\mathcal{K},

BEx=EAVxfor all x∈dom⁡(AV),BEx=EA_{\mathcal{V}}x \qquad\text{for all }x\in\operatorname{dom}(A_{\mathcal{V}}),

and Edom⁡(AV)E\operatorname{dom}(A_{\mathcal{V}}) is a core for BB. Let V′\mathcal{V}' denote the anti-dual of V\mathcal{V}, that is, the space of continuous conjugate-linear functionals on V\mathcal{V}. Define

M:V→V′,⟨Mx,y⟩V′,V=⟨Ex,Ey⟩Z,\mathcal{M}:\mathcal{V}\to\mathcal{V}', \qquad\langle\mathcal{M}x,y\rangle_{\mathcal{V}',\mathcal{V}}=\langle Ex,Ey\rangle_{\mathcal{Z}},
D:dom⁡(AV)⊂V→V′,Dx=MAVx.\mathcal{D}:\operatorname{dom}(A_{\mathcal{V}})\subset\mathcal{V}\to\mathcal{V}', \qquad\mathcal{D}x=\mathcal{M}A_{\mathcal{V}}x.

Then M\mathcal{M} is bounded and injective, D\mathcal{D} is dissipative in the duality sense, and

Mz˙=Dzin V(14)\mathcal{M}\dot{z}=\mathcal{D}z \quad\text{in }\mathcal{V} \tag*{(14)}

is equivalent to z˙=AVz\dot{z}=A_{\mathcal{V}}z for V\mathcal{V}-classical solutions.

Proof. Intertwining gives ∥T(t)x∥V2≤∥T(t)∥2∥x∥22+∥Ex∥22\|T(t)x\|_{\mathcal{V}}^2\leq\|T(t)\|^2\|x\|_2^2+\|Ex\|_2^2. Since T(t)x→xT(t)x\to x and ET(t)x=S(t)Ex→ExET(t)x=S(t)Ex\to Ex, the restricted semigroup is strongly continuous on V\mathcal{V}. The generator formula follows from [30 Proposition 2.4.4]. Differentiating ETV(t)x=S(t)ExET_{\mathcal{V}}(t)x=S(t)Ex gives BEx=EAVxBEx=EA_{\mathcal{V}}x. For sufficiently large λ\lambda,

E(λ−AV)−1=(λ−B)−1Eon V.E(\lambda-A_{\mathcal{V}})^{-1}=(\lambda-B)^{-1}E\quad\text{on }\mathcal{V}.

As ran⁡E\operatorname{ran}E is dense in K\mathcal{K}, the resolvent maps it onto a graph-norm dense subset of dom⁡(B)\operatorname{dom}(B), giving the core property. Finally, M\mathcal{M} is bounded and injective since EE is, and

2Re⁡⟨Dx,x⟩V′,V=ddt∣t=0+∥S(t)Ex∥Z2≤0for all x∈dom⁡(AV).2\operatorname{Re}\langle\mathcal{D}x,x\rangle_{\mathcal{V}',\mathcal{V}}=\left.\frac{\mathrm{d}}{\mathrm{d}t}\right|_{t=0+}\|S(t)Ex\|_{\mathcal{Z}}^2\leq0\quad\text{for all }x\in\operatorname{dom}(A_{\mathcal{V}}).

Injectivity of M\mathcal{M} yields the equivalence.

There is also a useful dual consequence of the lower observation estimate. It concerns the adjoint dynamics and must not be confused with an operator identity for AA on X\mathcal{X}.

Proposition 5.2 (Duality for bounded-below intertwiners). In addition to the hypotheses in Theorem 5.1, assume that ∥Ex∥≥c∥x∥\|Ex\|\geq c\|x\| on dom⁡(E)\operatorname{dom}(E) for some c>0c>0. Then H≥c2IH\geq c^2I and P=H−1∈L(X)P=H^{-1}\in\mathcal{L}(\mathcal{X}) is positive definite, with

T(t)PT(t)∗≤Pfor all t≥0.T(t)PT(t)^*\leq P\quad\text{for all }t\geq0.

Consequently, the adjoint generator A∗A^* admits a bounded-energy factorization A∗=DaPA^*=D_aP with maximally dissipative DaD_a.

Proof. Closedness and the lower bound imply ran⁡E=K\operatorname{ran}E=\mathcal{K}. Thus F=E−1:K→XF=E^{-1}:\mathcal{K}\to\mathcal{X} is bounded, with dense range, and T(t)F=FS(t)∣KT(t)F=FS(t)|_{\mathcal{K}}. Moreover P=FF∗P=FF^*, by the polar decomposition of EE. Therefore,

T(t)PT(t)∗=FS(t)S(t)∗F∗≤FF∗=P.T(t)PT(t)^*=FS(t)S(t)^*F^*\leq FF^*=P.

The last assertion follows by applying Theorem 3.1 to A∗A^* and PP.

From energy forms to operators on X\mathcal{X}. Contractivity of the energy form alone does not imply an operator factorization on the original state space. Indeed, if H=E∗EH=E^*E and x∈dom⁡(A)∩dom⁡(H)x\in\operatorname{dom}(A)\cap\operatorname{dom}(H), then

0≥∥H1/2T(t)x∥2−∥H1/2x∥2≥2Re⁡⟨T(t)x−x,Hx⟩.0\geq\|H^{1/2}T(t)x\|^2-\|H^{1/2}x\|^2\geq2\operatorname{Re}\langle T(t)x-x,Hx\rangle.

Dividing by tt and taking the limit to zero gives Re⁡⟨Ax,Hx⟩≤0\operatorname{Re}\langle Ax,Hx\rangle\leq0. If dom⁡(A)⊂dom⁡(H)\operatorname{dom}(A)\subset\operatorname{dom}(H), Theorem 3.1 therefore provides an exact representation A=DHA=DH with maximally dissipative DD. The domain inclusion is also necessary for an exact product representation with this particular HH, as the theorem shows. The first construction in Subsection 6.2 shows that it may fail even for a closed, injective output-trajectory map.

Remark 5.3 (Strong implicit equations and duality). Suppose that MM is a closed operator and an X\mathcal{X}-valued equation Mz˙=DizM\dot{z}=D_i z is required to hold at t=0t=0 for every classical trajectory of z˙=Az\dot{z}=Az. This requires ran⁡A⊂dom⁡(M)\operatorname{ran} A \subset\operatorname{dom}(M). If AA is bijective, then dom⁡(M)=X\operatorname{dom}(M)=\mathcal{X}, so the closed graph theorem forces MM to be bounded. Thus, an unbounded mass operator requires a smaller solution class or an equation in a different space, such as (14). The latter is not an operator-pair representation on X\mathcal{X} in the sense of Theorem 2.1. When HH is bounded and coercive, the change of variables w=Hxw=Hx transforms x˙=DeHx\dot{x}=D_eHx into H−1w˙=DewH^{-1}\dot{w}=D_ew. Without coercivity the change of variables need not be a boundedly invertible coordinate transformation of X\mathcal{X}.

Examples and limitations

The following examples illustrate the scope and limitations of the representation results.

All four boundedness combinations

Let I⊂(0,∞)I\subset(0,\infty) be an interval, set X=L2(I)\mathcal{X}=L^2(I), and let

(Ax)(s)=−sx(s),dom⁡(A)={x∈L2(I):sx(s)∈L2(I)}.(Ax)(s)=-sx(s),\qquad\operatorname{dom}(A)=\{x\in L^2(I): sx(s)\in L^2(I)\}.

The semigroup T(t)x(s)=e−stx(s)T(t)x(s)=e^{-st}x(s) is selfadjoint, contractive, and strongly stable. For the identity observation,

∥O0x∥L2(0,∞;X)2=12∫I∣x(s)∣2s ds.\lVert O_0x\rVert_{L^2(0,\infty;\mathcal{X})}^2=\frac{1}{2}\int_I\frac{\lvert x(s)\rvert^2}{s}\,\mathrm{d}s.

The corresponding energy operator is multiplication by 1/(2s)1/(2s). Consequently each of the following cases is realized:

IIO0O_0O0−1O_0^{-1} on ran⁡O0\operatorname{ran}O_0
(1,2)(1,2)boundedbounded
(1,∞)(1,\infty)boundedunbounded
(0,1)(0,1)unboundedbounded
(0,∞)(0,\infty)unboundedunbounded

Table 1.

Thus, the first case satisfies both the upper and lower observation estimates, the second only the upper estimate, the third only the lower estimate on its natural domain, and the fourth neither. In particular, a lower observation estimate does not imply admissibility.

In every case, the output map is closed, densely defined, and injective. Notice that these four patterns can occur although the original semigroup is already contractive. The distinction concerns the particular energy representation, not just the contractivity of the semigroup itself.

Limitations of canonical factorizations

Let X=ℓ2(N)\mathcal{X}=\ell^2(\mathbb{N}) and A=diag⁡(−1/n)A=\operatorname{diag}(-1/n). Then AA is bounded and T(t)=diag⁡(e−t/n)T(t)=\operatorname{diag}(e^{-t/n}) is contractive and strongly stable, but not exponentially stable. For C=IC=I,

∥OCx∥L22=12∑n≥1n∣xn∣2,H=OC∗OC=diag⁡(n/2).\lVert O_Cx\rVert_{L^2}^2=\frac{1}{2}\sum_{n\geq1}n\lvert x_n\rvert^2,\qquad H=O_C^*O_C=\operatorname{diag}(n/2).

The output map is closed, densely defined, injective, and bounded below, but HH is unbounded. More precisely, dom⁡(OC)={x∈ℓ2:∑n≥1n∣xn∣2<∞}\operatorname{dom}(O_C)=\{x\in\ell^2:\sum_{n\geq1}n\lvert x_n\rvert^2<\infty\}, whereas dom⁡(H)={x∈ℓ2:∑n≥1n2∣xn∣2<∞}\operatorname{dom}(H)=\{x\in\ell^2:\sum_{n\geq1}n^2\lvert x_n\rvert^2<\infty\}. Since dom⁡(A)=X⊄dom⁡(H)\operatorname{dom}(A)=\mathcal{X}\not\subset\operatorname{dom}(H), it cannot yield any exact factorization A=DHA=DH on X\mathcal{X}, even though the energy is nonincreasing. The obstruction concerns this particular HH; trivially, A=AIA=AI is already a dissipative representation. The same semigroup also illustrates the distinction between the canonical dissipative factor and its maximally dissipative extensions. Indeed, Hb=diag⁡(1/n)H_b = \operatorname{diag}(1/n) is bounded, positive definite, and a contractive energy, but De=−ID_e = -I on ran⁡Hb\operatorname{ran} H_b is not closed and hence not maximally dissipative. Nevertheless, D=−ID = -I on X\mathcal{X} is a maximally dissipative extension of DeD_e, and A=DHbA = D H_b holds with equality of domains.

The Dirichlet heat equation

Let Ω⊂Rd\Omega\subset\mathbb{R}^{d} be a bounded domain and let A=ΔDA = \Delta_D be the Dirichlet Laplacian on L2(Ω)L^2(\Omega). Then AA is selfadjoint and there exists λ1>0\lambda_1 > 0 such that

⟨Ax,x⟩≤−λ1∥x∥2for all x∈dom⁡(A).\langle Ax, x\rangle\le-\lambda_1\lVert x\rVert^2 \qquad\text{for all } x \in\operatorname{dom}(A).

Consequently, the corresponding heat semigroup is exponentially stable. By the functional calculus for selfadjoint operators (see [15 Chap. VI, Sec. 5]), the identity observation yields

H0=−12A−1,De=−2A2on dom⁡(A2).H_0 = -\frac{1}{2}A^{-1}, \qquad D_e = -2A^2 \quad\text{on } \operatorname{dom}(A^2).

The bounded operator H0H_0 is positive definite but not coercive. The implicit factor is Di=−12ID_i = -\frac{1}{2}I restricted to dom⁡(A)\operatorname{dom}(A). With observation C1=AC_1 = A, we instead obtain

H1=−12A,De=−2Ion X.H_1 = -\frac{1}{2}A, \qquad D_e = -2I \quad\text{on } \mathcal{X}.

The closure of the output map initially defined by O1x=T(⋅)AxO_1x = T(\cdot)Ax on dom⁡(A)\operatorname{dom}(A) has domain dom⁡((−A)1/2)\operatorname{dom}((-A)^{1/2}), whereas dom⁡(H1)=dom⁡(A)\operatorname{dom}(H_1) = \operatorname{dom}(A). The energy operator H1H_1 is unbounded and coercive, but H1A=−12A2H_1A = -\frac{1}{2}A^2 is defined only on dom⁡(A2)\operatorname{dom}(A^2), not on all of dom⁡(A)\operatorname{dom}(A). On the other hand, for M=H0M = H_0 and Di=H0A=−12I∣dom⁡(A)D_i = H_0A = -\frac{1}{2}I|_{\operatorname{dom}(A)}, the operator pencil λH0−Di\lambda H_0 - D_i has range ran⁡H0=dom⁡(A)≠X\operatorname{ran} H_0 = \operatorname{dom}(A) \ne\mathcal{X} whenever λI−A\lambda I - A is bijective, while (9) recovers the heat resolvent. Thus, both the explicit and implicit representations can be meaningful without a boundedly invertible energy operator.

Exponential stability without Hilbertian similarity

Exponentially stable analytic Hilbert-space semigroups need not be similar to contraction semigroups [16, 17, 24]. More strikingly, there exist strongly continuous semigroups that are nilpotent and immediately compact, as well as strongly continuous semigroups that are quasi-nilpotent, immediately compact, and bounded holomorphic of angle π/2\pi/2. Neither class is similar even to a quasi-contraction semigroup. The construction uses a tensor-product splitting principle: similarity of a tensor product to a quasi-contraction semigroup forces the same property for its nonzero factors [25 Theorems 1.1 and 1.3].

Both classes of examples are exponentially stable. Theorem 3.5 therefore furnishes the bounded, positive definite Gramian H0H_0 from (10) and an exact factorization A=DeH0A = D_eH_0 with maximally dissipative DeD_e. For a nilpotent semigroup, the Gramian integral even terminates at a finite extinction time, and in both cases the strict identity 2Re⁡{H0Ax,x}=−∥x∥22\operatorname{Re}\{H_0Ax,x\} = -\lVert x\rVert^2 holds for all x∈dom⁡(A)x \in\operatorname{dom}(A). Nevertheless, H0H_0 cannot be coercive: any bounded, coercive energy operator satisfying (7) would yield a boundedly invertible similarity to a contraction semigroup, contradicting the examples. Thus strict Lyapunov dissipation, even together with very strong stability and compactness, does not imply an equivalent contractive Hilbertian norm.

The one-sided output constructions in Theorem 4.4 reflect the same obstruction: the identity output is bounded with unbounded inverse on its range, whereas the closure of the generator output is unbounded but bounded below.

Concluding remarks

We have characterized exact dissipative Hamiltonian factorizations of semigroup generators in terms of a Lyapunov inequality and a domain inclusion. Whenever such a factorization exists, its dissipative factor can be chosen maximally dissipative without additional assumptions. The canonical factor need not itself be maximally dissipative, but every dissipative extension preserves the exact operator factorization. For bounded energy operators, the same Lyapunov inequality also yields an implicit representation.

Closed intertwiners and output-trajectory maps provide energy operators and contractive dynamics on associated energy spaces. Bounded-below intertwiners additionally yield a dual Lyapunov inequality and a dissipative Hamiltonian factorization of the adjoint generator. These constructions do not, in general, guarantee an exact factorization of the original generator with the constructed energy operator. Domain compatibility remains essential.

Acknowledgments. Both authors gratefully acknowledge funding by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) – Project-ID 531152215 – CRC 1701 “Port-Hamiltonian Systems”.

References

  1. [1]F. Achleitner, A. Arnold, and E. A. Carlen, The hypocoercivity index for the short time behavior of linear time-invariant ODE systems, J. Differential Equations, 371 (2023), pp. 83–115, https://doi.org/10.1016/j.jde.2023.06.027.
  2. [2]A. Arnold, S. Egger, V. Mehrmann, and E. A. Nigsch, Connection of hypocoercivity and hypocontractivity via the Cayley transform. arXiv:2601.09656, 2026.arxiv.org/abs/2601.09656
  3. [3]J. A. Ball, S. ter Horst, and M. Kurula, The infinite-dimensional standard and strict bounded real lemmas in continuous time: The storage function approach, Complex Anal. Oper. Theory, 16 (2022), p. 84, https://doi.org/10.1007/s11785-022-01260-5.
  4. [4]C. Beattie, V. Mehrmann, H. Xu, and H. Zwart, Linear port-Hamiltonian descriptor systems, Math. Control Signals Syst., 30 (2018), p. 17, https://doi.org/10.1007/s00498-018-0223-3.
  5. [5]C. Beattie, V. Mehrmann, H. Xu, and H. Zwart, Optimally robust representations of linear time-invariant asymptotically stable and strictly passive systems. in preparation, 2026.
  6. [6]A. Bendimerad-Hohl, D. Matignon, G. Haine, and L. Lefèvre, On Stokes–Lagrange and Stokes–Dirac representations for 1D distributed port-Hamiltonian systems, IFAC-PapersOnLine, 58 (2024), pp. 238–243, https://doi.org/10.1016/j.ifacol.2024.10.174.
  7. [7]R. Datko, Extending a theorem of A. M. Liapunov to Hilbert space, J. Math. Anal. Appl., 32 (1970), pp. 610–616, https://doi.org/10.1016/0022-247X(70)90283-0.
  8. [8]K.-J. Engel and R. Nagel, One-Parameter Semigroups for Linear Evolution Equations, vol. 194 of Graduate Texts in Mathematics, Springer-Verlag, New York, 2000.
  9. [9]M.-C. Gao and J.-C. Hou, The infinite-time admissibility of observation operators and operator Lyapunov equations, Integral Equations Operator Theory, 35 (1999), pp. 53–64, https://doi.org/10.1007/BF01225527.
  10. [10]H. Gernandt, F. E. Haller, and T. Reis, A linear relation approach to port-Hamiltonian differential-algebraic equations, SIAM J. Matrix Anal. Appl., 42 (2021), pp. 1011–1044, https://doi.org/10.1137/20M1371166.
  11. [11]P. Grabowski and F. M. Callier, Admissible observation operators. Semigroup criteria of admissibility, Integral Equations Operator Theory, 25 (1996), pp. 182–198, https://doi.org/10.1007/ B. H. Haak and E. M. Ouhabaz, Exact observability, square functions and spectral theory, J. Funct. Anal., 262 (2012), pp. 2903–2927, https://doi.org/10.1016/j.jfa.2012.01.007.
  12. [13]B. Jacob and J. R. Partington, Admissibility of control and observation operators for semigroups: a survey, in Current Trends in Operator Theory and Its Applications, J. A. Ball, J. W. Helton, M. Klaus, and L. Rodman, eds., vol. 149 of Oper. Theory Adv. Appl., Birkhäuser, Basel, 2004, pp. 199–221, https://doi.org/10.1007/978-3-0348-7881-4_10.
  13. [14]B. Jacob and H. Zwart, Linear Port-Hamiltonian Systems on Infinite-Dimensional Spaces, vol. 223 of Operator Theory: Advances and Applications, Birkhäuser, Basel, 2012.DOI
  14. [15]T. Kato, Perturbation Theory for Linear Operators, Springer-Verlag, Heidelberg, 2nd ed., 1980.
  15. [16]C. Le Merdy, The similarity problem for bounded analytic semigroups on Hilbert space, Semigroup Forum, 56 (1998), pp. 205–224, https://doi.org/10.1007/PL00005942.
  16. [17]C. Le Merdy, A bounded compact semigroup on Hilbert space not similar to a contraction one, in Semigroups of Operators: Theory and Applications, A. V. Balakrishnan, ed., vol. 42 of Progr. Non-linear Differential Equations Appl., Birkhäuser, Basel, 2000, pp. 213–216, https://doi.org/10.1007/978-3-0348-8417-4_21.
  17. [18]N. Levan, The left shift semigroup approach to stability of distributed systems, J. Math. Anal. Appl., 152 (1990), pp. 354–367, https://doi.org/10.1016/0022-247X(90)90070-V.
  18. [19]G. Lumer and R. S. Phillips, Dissipative operators in a Banach space, Pacific J. Math., 11 (1961), pp. 679–698, https://doi.org/10.2140/pjm.1961.11.679.
  19. [20]C. Mehl, V. Mehrmann, and M. Wojtylak, Linear algebra properties of dissipative Hamiltonian descriptor systems, SIAM J. Matrix Anal. Appl., 39 (2018), pp. 1489–1519, https://doi.org/10.1137/18M1164275.
  20. [21]C. Mehl, V. Mehrmann, and M. Wojtylak, Spectral theory of infinite dimensional dissipative Hamiltonian systems, J. Dynam. Differential Equations, 38 (2026), pp. 551–577, https://doi.org/10.1007/s10888-025-010420-y.
  21. [22]V. Mehrmann and A. van der Schaft, Differential-algebraic systems with dissipative Hamiltonian structure, Math. Control Signals Syst., 35 (2023), pp. 541–584, https://doi.org/10.1007/s00498-023-00349-2.
  22. [23]V. Mehrmann and A. van der Schaft, Nonlinear differential-algebraic systems with dissipative port-Hamiltonian structure. arXiv:2610.11455, 2026.arxiv.org/abs/2610.11455
  23. [24]J. Oliva-Maza and Y. Tomilov, Similarity to contraction semigroups: structural properties, criteria, and applications to control theory. arXiv:2509.01007, 2025.
  24. [25]J. Oliva-Maza and Y. Tomilov, On similarity to contraction semigroups and tensor products, I, J. Funct. Anal., 291 (2026), p. 111492, https://doi.org/10.1016/j.jfa.2026.111492.
  25. [26]J. R. Partington and S. Pott, Admissibility and exact observability of observation operators for semigroups, Irish Math. Soc. Bull., 55 (2005), pp. 19–39, https://doi.org/10.33232/BIMS.0055.19.40.
  26. [27]F. M. Philipp, T. Reis, and M. Schaller, Infinite-dimensional port-Hamiltonian systems: a system node approach, Math. Control Signals Syst., 37 (2025), pp. 573–620, https://doi.org/10.1007/s00498-025-00412-0.
  27. [28]T. Reis, The infinite-dimensional dissipation inequality, ESAIM Control Optim. Calc. Var., 32 (2026), p. 70, https://doi.org/10.1051/cocv/2026053.
  28. [29]O. Staffans, Well-Posed Linear Systems, vol. 103 of Encyclopedia of Mathematics and its Applications, Cambridge University Press, Cambridge, 2005.
  29. [30]M. Tucsnak and G. Weiss, Observation and Control for Operator Semigroups, Birkhäuser Advanced Texts: Basler Lehrbücher, Birkhäuser, Basel, 2009, https://doi.org/10.1007/978-3-7643-8994-9.
  30. [31]G. Weiss, Admissible observation operators for linear semigroups, Israel J. Math., 65 (1989), pp. 17–43, https://doi.org/10.1007/BF02788172.
  31. [32]H. Zwart and V. Mehrmann, Abstract dissipative Hamiltonian differential-algebraic equations are everywhere, DAE Panel, 2 (2024), https://doi.org/10.52825/dae-p.v2i.957.

Paper details

vlnc:W_fxpdwvbt
v1 · latestHistory (1)

published · Fresh — added to the Hub within the last 30 days

added ·

Usage

views · 0

downloads · 0

A view counts once per browser per day, a download each time the PDF is requested.

CC BY 4.0

Original source

  • Mathematics

Not recorded for this paper · View statement →

Looking for related papers…
Contents