Convergence of proximal solutions for evolution inclusions with time-dependent maximal monotone operators

TL;DR

This paper investigates the convergence of discretized solutions to evolution inclusions with time-dependent maximal monotone operators, establishing conditions for well-posedness and solution uniqueness using variational analysis.

Contribution

It introduces a novel convergence analysis for discretized solutions of time-dependent differential inclusions with maximal monotone operators, extending existing theory.

Findings

Discretized solutions converge to the unique continuous solution.

Provided conditions for well-posedness of related nonsmooth differential equations.

Applied passivity and linear matrix inequalities to establish solution properties.

Abstract

This article studies the solutions of time-dependent differential inclusions which is motivated by their utility in the modeling of certain physical systems. The differential inclusion is described by a time-dependent set-valued mapping having the property that, for a given time instant, the set-valued mapping describes a maximal monotone operator. Under certain mild assumptions on the regularity with respect to the time argument, we construct a sequence of functions parameterized by the sampling time that corresponds to the discretization of the continuous-time system. Using appropriate tools from functional and variational analysis, this sequence is then shown to converge to the unique solution of the original differential inclusion. The result is applied to develop conditions for well-posedness of differential equations interconnected with nonsmooth time-dependent complementarity relations, using passivity of underlying dynamics (equivalently expressed in terms of linear matrix inequalities).

Equations

$$\dot{x}\in -F(t,x),x(0)\in \mathrm{dom}F(0,\cdot ),$$

$$\displaystyle d_{\operatorname{H}}(S_{1},S_{2}):=\max{\big{\{}\sup_{z_{1}\in S_{1}}\operatorname{dist}(z_{1},S_{2}),\sup_{z_{2}\in S_{2}}\operatorname{dist}(z_{2},S_{1})\big{\}}}.$$

$$d_{\operatorname{H}}(S_{1},S_{2})=d_{\operatorname{H}}\big{(}\operatorname{cl}(S_{1}),\operatorname{cl}(S_{2})\big{)}.$$

$$|x-y|⩽\underset{z\in S_1}{\sup }\mathrm{dist}(z,S_2)\le d_H(S_1,S_2).$$

$$\mathrm{dom}F=\{x:F(x)\ne \varnothing \}$$

$$\mathrm{im}F=\{y:\text{there exists }x\text{ such that }y\in F(x)\}$$

$$\mathrm{graph}F=\{(x,y):y\in F(x)\}.$$

$$\underset{\ell \to \infty }{\mathrm{limsup}}{S}_{\ell }:=\left\{\xi |\exists N\in {\mathcal{N}}_{\infty }^{\#}\text{ and }{\xi }_{\ell }\in S_{\ell }\forall \ell \in N,\text{ s.t. }{\xi }_{\ell }\stackrel{N}{\to }\xi \right\}.$$

$$\underset{t\searrow t^{\ast }}{\mathrm{limsup}}G(t):=\bigcup _{t_{\ell }\searrow t^{\ast }}\underset{\ell \to \infty }{\mathrm{limsup}}G(t_{\ell }).$$

$$\underset{t\searrow t^{\ast }}{\mathrm{limsup}}G(t)=\left\{y\in {\mathbb{R}}^q|\begin{array}{rl}& \exists (t_{\ell },y_{\ell }{)}_{\ell \in \mathbb{N}}\subset [0,T]\times {\mathbb{R}}^q\text{ satisfying }{y}_{\ell }\in G(t_{\ell }),\\ & t_{\ell }\ge t,\text{ and }\underset{\ell \to \infty }{\lim }(t_{\ell },y_{\ell })=(t^{\ast },y)\end{array}\right\}.$$

$$\underset{t\searrow t^{\ast }}{\mathrm{limsup}}G(t)\subseteq G(t^{\ast }).$$

$$⟨x_1-x_2,y_1-y_2⟩⩾0$$

$$F^0(x):=\mathrm{proj}(0,F(x))$$

$$J_{\lambda }:=(I+\lambda F{)}^{-1}\text{and}{Y}_{\lambda }:=\frac{1}{\lambda }(I-J_{\lambda })$$

$$\operatorname{dis}(F_{1},F_{2}):=\sup_{\small\begin{array}[]{l}x_{1}\in\operatorname{dom}(F_{1}),y_{1}\in F_{1}(x_{1})\\ x_{2}\in\operatorname{dom}(F_{2}),y_{2}\in F_{2}(x_{2})\end{array}}{\dfrac{\langle y_{1}-y_{2},x_{2}-x_{1}\rangle}{1+|y_{1}|+|y_{2}|}}.$$

$$d_H(\mathrm{dom}(F_1),\mathrm{dom}(F_2))⩽\mathrm{dis}(F_1,F_2).$$

$$\displaystyle\operatorname{dis}\big{(}F(t,\cdot),F(s,\cdot)\big{)}\leqslant\varphi(t)-\varphi(s)\quad\forall\,s,t\text{ with }0\leqslant s\leqslant t\leqslant T.$$

$$\frac{x({\tau }_k)-x(t_k)}{{\tau }_k-t_k}=\frac{x({\tau }_k)-x(t^{\ast })}{{\tau }_k-t^{\ast }}\frac{{\tau }_k-t^{\ast }}{{\tau }_k-t_k}+\frac{x(t^{\ast })-x(t_k)}{t^{\ast }-t_k}\frac{t^{\ast }-t_k}{{\tau }_k-t_k}.$$

$$\xi (t)=x(0)+{\int }_0^{t}y(s)\dot{\psi }(s)ds$$

$$⟨\eta ,x_{\ell }(t)⟩=⟨\eta ,x_0⟩+{\int }_0^t⟨\eta ,y_{\ell }(s)⟩\dot{\psi }(s)ds$$

$$⟨\eta ,\xi (t)⟩=⟨\eta ,x_0⟩+{\int }_0^t⟨\eta ,y(s)⟩\dot{\psi }(s)ds.$$

$$x(t)=x_0+{\int }_0^{t}y(s)\dot{\psi }(s)ds.$$

$$\dot{x}(t)=\dot{\psi }(t)y(t)$$

$$f(t)\in [\text{smashoperator}[]\underset{\ell \to \infty }{\mathrm{liminf}}{f}_{\ell }(t),\text{smashoperator}[]\underset{\ell \to \infty }{\mathrm{limsup}}{f}_{\ell }(t)]$$

$$y(t)\in\operatorname{cl}\Big{(}\operatorname{conv}\big{(}\limsup_{\ell\rightarrow\infty}S_{\ell}(t)\big{)}\Big{)}$$

$$z(t)=\operatorname{proj}\big{(}y(t),S(t)\big{)}$$

$$a(t)=\frac{y(t)-z(t)}{|y(t)-z(t)|}\text{and}b(t)=⟨\frac{y(t)-z(t)}{|y(t)-z(t)|},\frac{y(t)+z(t)}{2}⟩$$

$$⟨a(t),y(t)⟩<b(t)<⟨a(t),z⟩$$

$$\zeta (t)\in [\text{smashoperator}[]\underset{\ell \to \infty }{\mathrm{liminf}}{\zeta }_{\ell }(t),\text{smashoperator}[]\underset{\ell \to \infty }{\mathrm{limsup}}{\zeta }_{\ell }(t)]$$

$$⟨a(t),y(t)⟩\in [⟨a(t),\underline{y}(t)⟩,⟨a(t),\overline{y}(t)⟩]$$

$$⟨a(t),y(t)⟩>b(t)$$

Full text

1 Introduction

The theory of monotone operators emerged as an important area of research within the field of nonlinear analysis in early 1960’s Kacu60 ; Mint62 ; Zara60 . Since then, we have seen applications of such operators in various disciplines, which include, but are not limited to, optimization algorithms, dynamical systems, and partial differential equations. Recent articles attouch2019 ; Comb18 ; RyuBoyd16 provide an overview of monotone operators appearing in optimization algorithms. The relevance of set-valued mappings in dynamical systems observed in Brez73 ; More77 , where the differential inclusions with maximal monotone operators are analyzed. Even in the systems of partial differential equations, the appearance of these operators brings tractability to proving existence of solutions Brow63 ; LeraLion65 ; Show93 ; Zeid90 . Applications of dynamical systems with maximal monotone operators range from modeling traffic equilibrium pang2012 to electronics addi11 . Relatively modern texts on analysis of monotone operators are BausComb17 ; Phel93 ; Simo08 .

This article is focused on studying maximal monotone operators in the context of mathematical models for dynamical systems, and the central object of our study is to investigate conditions for existence of solutions to the differential inclusion

[TABLE]

where $F:[0,\infty)\times\mathbb{R}^{n}\rightrightarrows\mathbb{R}^{n}$ has the property that, for each $t\geqslant 0$, $F(t,\cdot)$ is a maximal monotone operator. In studying this generic class of systems, we will refer to other types of nonsmooth dynamical systems which can be recast in the form (1). From a theoretical point of view, most of the earlier work had focused on differential inclusions with static maximal monotone operators, which is very elegantly collected in Brez73 , or see PeypSori10 for a recent overview on this subject. Common techniques used in analyzing such systems are either based on regularization, or discretization. For the former one, the so-called Yosida-Moreau approximations provide a single-valued Lipschitz function with a regularization parameter, and as this parameter converges to zero, it is shown that the corresponding solutions converge to the solution of the original differential inclusion. The discretization techniques rely on constructing piecewise constant interpolations of the sequence of points obtained from some discrete system with a sampling parameter. As the sampling parameter converges to zero, the corresponding sequence of solutions is shown to converge to the actual solution.

To the best of our knowledge, the first attempts for studying inclusion (1) with time-dependent operators $F(t,x)$, but with the domain of $F(t,\cdot)$ stationary for each $t\geqslant 0$, were carried out in Kato67 . Since then, several works have appeared which tackle dynamical systems with time-dependent multi-valued monotone operators. When $F(t,x)$ is the subdifferential of a time-dependent, proper, lower semicontinuous, and convex function $\varphi_{t}(\cdot)$, that is, $F(t,x)=\partial\varphi_{t}(x)$, then $F(t,\cdot)$ is a maximal monotone operator. Such systems, involving time-dependent subdifferentials, have been particularly studied in Arseni-Benou99 ; kandilakis96 ; kartsatos84 ; Kuttler00 ; otani94 ; yamazaki05 under varying degrees of regularity on the system data. Imposing further structure on the operator $F(t,\cdot)$, if we take $F(t,x)=\partial\psi_{S(t)}(x)$, where $S:[0,\infty)\rightrightarrows\mathbb{R}^{n}$ is closed and convex-valued mapping and $\psi$ is the indicator function associated with $S(t)$, then the resulting dynamics have been more commonly studied under the topic of sweeping processes. Starting from the seminal work of More77 , the research in this area has grown to study several generalizations of the fundamental model, see for example, the monographs Adly18 ; KunzMont00 ; Mont93 ; SiddManc02 for an overview, and the articles AdlyHadd18 ; AdlyHadd14 ; EdmoThib05 ; EdmoThib06 ; JourVilc18 ; KrejRoch11 ; Recu15 for more recent and focused expositions. Besides the cases where $F$ is expressed as a subdifferential of a convex function, certain classes of evolution variational inequalities BrogGoel11 ; PangStew08 ; TanwBrog18 can also be embedded in the framework of (1).

While all these aforementioned works can be represented by (1), they also rely on the particular structure of the set-valued mapping in their problem description for analysis of existence of solutions. Notable exceptions in the literature, which address directly the system class (1) are KunzMont97 ; Vlad91 . However, the regularity assumptions imposed in these works restrict the applicability of their results. Consequently, when applications of these dynamical models are studied, for example in control Briceno-Arias15 ; TanwBrog18 , the results that build on the works of KunzMont97 ; Vlad91 suffer similar limitations. Based on these observations, the motivation to study new set of conditions for existence of solutions to systems class (1) arises and our aim in this paper is to provide mild (read as mildest possible) conditions on regularity of system data, which allow us to cover a possibly larger class of systems. Moreover, we can recover most (if not all) of the results on time-dependent and static case with our approach.

Our approach builds on using the time-stepping algorithm pioneered in More77 , which was also used for studying existence of solutions for system (1) in KunzMont97 . This algorithm constructs a sequence of solutions, where each element of the sequence is an interpolation of points obtained by applying the proximal operator associated with domain of the multivalued function appearing in (1). With the help of an academic counterexample, we show how the assumptions imposed in KunzMont97 fail to hold for a dynamical system described by time-dependent complementarity relations. We study existence of solutions under conditions which overcome such restrictions. The basic idea is to construct a sequence of solutions. To construct an element of this approximate solution with a fixed sampling time, we first compute a set of points at sampled time instants by projecting the value of a certain function on the domain of the set-valued mapping. Using a novel interpolation technique among these discrete points, we obtain a sequence of absolutely continuous functions. Using the arguments based on Ascolà-Arzeli theorem, this sequence is shown to converge to an absolutely continuous function, which is then shown to be the unique solution of the original system. We generalize our result to the case where the right-hand side of (1) has a single-valued Lipschitz vector field in addition to the set-valued maximal monotone operator.

Moreover, because of the relaxed nature of assumptions, our results provide a constructive framework for studying differential equations with complementarity relations. Such nonsmooth relations form a particular subclass of maximal monotone operators, and have been useful in modeling systems with piecewise affine characteristics j06 ; Brog03 ; BrogThib10 . Earlier work on complementarity systems has focused on linear dynamics coupled with static complementarity relations CamlHeem02 ; Heem00 ; schumacher2004 ; j16 . Lately, it was shown in CamlSchu16 that an interconnection of static complementarity relation with ordinary differential equations yields a differential inclusion with static maximal monotone operator. However, time-dependence in complementarity relations has not been easy to treat with existing frameworks. Inspired by the result in CamlSchu16 , we provide conditions under which it is possible to recast the interconnection of an ordinary differential equation with time-dependent complementarity relation in the form of a differential inclusion with time-dependent maximal monotone operator, for which the existence of solutions is being studied in this article.

The remainder of the article is organized as follows: In Section 2, we provide appropriate background material from set-valued and functional analysis. In Section 3, a motivating example is provided to show how the current literature on differential inclusions with maximal monotone operators is inadequate for certain system classes. The main assumptions and the result is given in Section 4, followed by a detailed proof in Section 5. Section 6 deals with extensions of the main existence/uniqueness result towards non-autonomous case as well as Lipschitzian perturbations. The results are then studied in the context of linear ordinary differential equations coupled with time-dependent maximal monotone relations in Section 7. Finally, the paper closes with some concluding remarks in Section 8.

2 Preliminaries

In this section, we introduce notational conventions that will be in force throughout the paper as well as auxiliary results that will be employed later.

2.1 Vectors and matrices

We denote the set of real numbers by $\mathbb{R}$, nonnegative real numbers by $\mathbb{R}_{+}$, $n$-vectors of real numbers by $\mathbb{R}^{n}$, and $n\times m$ real-valued matrices by $\mathbb{R}^{n\times m}$.

To denote the scalar product of two vectors $x$, $y\in\mathbb{R}^{n}$, we use the notation $\left\langle x,y\right\rangle:=x^{\top}y$ where $x^{\top}$ denotes the transpose of $x$. The Euclidean norm of a vector $x$ is denoted by $|x|:=\left\langle x,x\right\rangle^{\frac{1}{2}}$. For a subspace of $\mathcal{W}$ of $\mathbb{R}^{n}$, $\mathcal{W}^{\perp}$ denotes the orthogonal subspace, that is $\{y\in\mathbb{R}^{n}:\left\langle x,y\right\rangle=0\text{ for all }x\in\mathcal{W}\}$.

We say that a (not necessarily symmetric) matrix $M\in\mathbb{R}^{n\times n}$ is positive semi-definite if $x^{\top}Mx\geqslant 0$ for all $x\in\mathbb{R}^{n}$. We sometimes write $M\geqslant 0$ meaning that $M$ is positive semi-definite. Also, we say that $M$ is positive definite if $M>0$ for all nonzero $x\in\mathbb{R}^{n}$.

2.2 Convex sets and related notions

The distance of a point $x$ to a set $S$ is defined by $\operatorname{dist}(x,S)=\inf\{|x-y|:y\in S\}$. If the set $S$ is closed and convex then for each $x\in\mathbb{R}^{n}$ there exists a unique point $y\in S$ such that $|x-y|=\operatorname{dist}(x,S)$. Such a point is called the projection of $x$ onto the set $S$ and will be denoted by $\operatorname{proj}(x,S)$.

The Hausdorff distance between two nonempty subsets of $\mathbb{R}^{m}$, say $S_{1}$ and $S_{2}$, is defined by:

[TABLE]

Since \operatorname{dist}(x,S)=\operatorname{dist}\big{(}x,\operatorname{cl}(S)\big{)} for any point $x$ and nonempty set $S$, the Hausdorff distance is invariant under closure, that is

[TABLE]

In addition, if $y=\operatorname{proj}(x,\operatorname{cl}(S_{2}))$ for some point $x\in\operatorname{cl}(S_{1})$, then we have

[TABLE]

2.3 Set-valued mappings

Let $F:\mathbb{R}^{m}\rightrightarrows\mathbb{R}^{n}$ be a set-valued mapping, that is $F(x)\subseteq\mathbb{R}^{n}$ for each $x\in\mathbb{R}^{m}$. We define its domain, image, and graph, respectively, as follows:

[TABLE]

The inverse mapping $F^{-1}:\mathbb{R}^{n}\rightrightarrows\mathbb{R}^{m}$ is defined by $F^{-1}(y)=\{x:y\in F(x)\}$.

In what follows we introduce a certain notion of continuity for set-valued mappings of a real variable. For a more detailed/general treatment we refer to (Rockafellar98, , Chp. 4 and 5).

Let $\mathcal{N}_{\infty}^{\#}$ denote the set of all subsequences of $\mathbb{N}$. For a sequence of sets $(S_{\ell})_{\ell\in\mathbb{N}}$ in $\mathbb{R}^{q}$, the outer limit is defined as the set

[TABLE]

For a given set-valued mapping $G:[0,T]\rightrightarrows\mathbb{R}^{q}$ for some $T>0$, we define

[TABLE]

It is known from (Rockafellar98, , p. 152) that

[TABLE]

We say that $G$ is outer semicontinuous from right at $t^{*}\in[0,T]$ if

[TABLE]

In case $G$ is outer semicontinuous from right at every $t^{*}\in[0,T]$, we say that $G$ is outer semicontinuous from right on $[0,T]$.

2.4 Maximal monotone operators

Throughout the paper, we are interested in maximal monotone set-valued mappings. A set valued-mapping $F:\mathbb{R}^{n}\rightrightarrows\mathbb{R}^{n}$ is said to be monotone if

[TABLE]

for all $(x_{i},y_{i})\in\operatorname{graph}(F)$. It is said to be maximal monotone if no enlargement of its graph is possible in $\mathbb{R}^{n}\times\mathbb{R}^{n}$ without destroying monotonicity. We refer to Brez73 and Rockafellar98 for detailed treatment of maximal monotone mappings.

If $F$ is maximal monotone, then it is closed and convex-valued, that is, $F(x)$ is a closed convex set for all $x\in\operatorname{dom}(F)$. This enables us to define the minimal section of a maximal monotone mapping $F$ by

[TABLE]

for $x\in\operatorname{dom}(F)$. Clearly, $F^{0}(x)$ is the least-norm element of the closed convex set $F(x)$, that is $|F^{0}(x)|\leqslant|y|$ for all $y\in F(x)$.

The *resolvent * $J_{\lambda}$ and Yosida approximation $Y_{\lambda}$ of $F$ are defined by

[TABLE]

for $\lambda>0$ where $I$ denotes the identity operator.

The following proposition collects some well-known facts (see e.g. Brez73 ) that will be employed in the sequel.

Proposition 1**.**

Suppose that $F:\mathbb{R}^{n}\rightrightarrows\mathbb{R}^{n}$ is a maximal monotone set-valued mapping. Then, the following statements hold for all $\lambda>0$:

1. i.

$\operatorname{dom}J_{\lambda}=\mathbb{R}^{n}$. 2. ii.

$J_{\lambda}$* is single-valued and non-expansive, that is $|J_{\lambda}(x_{1})-J_{\lambda}(x_{2})|\leqslant|x_{1}-x_{2}|$ for all $x_{1}$, $x_{2}\in\mathbb{R}^{n}$.* 3. iii.

$\lim_{\lambda\rightarrow 0}J_{\lambda}(x)=x$* for all $x\in\operatorname{dom}(F)$.* 4. iv.

$Y_{\lambda}$* is maximal monotone and *$\lambda^{-1}-$Lipschitzian. 5. v.

Y_{\lambda}(x)\in F\big{(}J_{\lambda}(x)\big{)}* for all $x\in\mathbb{R}^{n}$.* 6. vi.

For all $x\in\operatorname{dom}F$, $|Y_{\lambda}(x)|$ is nonincreasing in $\lambda$, $\lim_{\lambda\rightarrow 0}|Y_{\lambda}(x)|=|F^{0}(x)|$, and $|Y_{\lambda}(x)|\leqslant|F^{0}(x)|$.

Given two maximal monotone mappings, the pseudo-distance between them, introduced in Vlad91 , is defined as follows:

Definition 2**.**

The pseudo-distance between two maximal monotone mappings $F_{1}$ and $F_{2}$ is defined by

[TABLE]

The following lemma relates the Hausdorff distance between the domains of two maximal monotone operators with their pseudo-distance.

Lemma 3** (Vlad91 ).**

For any pair of maximal monotone mappings $F_{1}$ and $F_{2}$, it holds that

[TABLE]

Based on the pseduo-distance defined in Definition 2, one can introduce a notion of continuity for time-dependent maximal monotone operators as follows.

Definition 4** (Absolute continuity, Vlad91 ).**

Let $F:[0,T]\times\mathbb{R}^{n}\rightrightarrows\mathbb{R}^{n}$ be a time-dependent set-valued mapping such that $F(t,\cdot)$ is maximal monotone for each $t\in[0,\,T]$. We say that $t\mapsto F(t,\cdot)$ is absolutely continuous on $[0,T]$ if there exists a nondecreasing absolutely continuous function $\varphi\;:\;[0,\,T]\rightarrow\mathbb{R}$ such that

[TABLE]

2.5 Function spaces

The set of absolutely continuous, integrable, and square integrable functions defined from the interval $[t_{1},t_{2}]$ with $t_{1}<t_{2}$ to $\mathbb{R}^{n}$ are denoted, respectively, by $AC([t_{1},t_{2}],\mathbb{R}^{n})$, $L_{1}([t_{1},t_{2}],\mathbb{R}^{n})$, and $L_{2}([t_{1},t_{2}],\mathbb{R}^{n})$. Unless specified otherwise, we use the term almost everywhere with respect to Lebesgue measure, that is, a property holds almost everywhere on a set $X\subset\mathbb{R}^{n}$, if it holds on every subset of $X$ with nonzero Lebesgue measure.

Convergence of family of functions will play an important role in the sequel. For the sake of completeness, we state the well-known (see e.g. rudin ) Arzelá-Ascoli theorem for which we need some nomenclature.

Consider a collection $\mathcal{F}$ of functions $f:[0,T]\to\mathbb{R}^{n}$. We say that $\mathcal{F}$ is equicontinuous if for every $\varepsilon>0$, there exists a $\delta>0$ such that $|f(t)-f(s)|<\varepsilon$ for every $f\in\mathcal{F}$ and each $s,t$ satisfying $|t-s|<\delta$. We say that $\mathcal{F}$ is pointwise bounded if for every $t\in[0,T]$, there exists an $M_{t}<\infty$ such that $|f(t)|\leq M_{t}$ for every $f\in\mathcal{F}$.

Theorem 5** (Arzelá-Ascoli).**

Suppose that $\mathcal{F}$ is pointwise bounded equicontinuous collection of functions $f:[0,T]\to\mathbb{R}^{n}$. Every sequence $\{f_{n}\}$ in $\mathcal{F}$ has a subsequence that converges uniformly on every compact subset of $[0,T]$.

The following elementary results will be used later.

Lemma 6**.**

Let $x:[0,T]\rightarrow\mathbb{R}^{n}$ be a function and $t^{*}\in(0,T]$ be such that $\dot{x}(t^{*})$ exists. Suppose that $\{t_{k}\}$ and $\{\tau_{k}\}$ are two sequences such that $0\leqslant t_{k}\leqslant t^{*}\leqslant\tau_{k}\leqslant T$ and $t_{k}<\tau_{k}$ for all $k$ and $\lim_{k\uparrow\infty}t_{k}=\lim_{k\uparrow\infty}\tau_{k}=t^{*}$. Then, the sequence $\frac{x(\tau_{k})-x(t_{k})}{\tau_{k}-t_{k}}$ convergences to $\dot{x}(t^{*})$ on a subsequence.

**Proof. **Observe that

[TABLE]

Since $0\leqslant\frac{\tau_{k}-t^{*}}{\tau_{k}-t_{k}}\leqslant 1$ and $0\leqslant\frac{t^{*}-t_{k}}{\tau_{k}-t_{k}}\leqslant 1$, both must converge on a (common) subsequence. The rest follows from the hypothesis that $\dot{x}(t^{*})$ exists. $\blacksquare$

Lemma 7**.**

Suppose that a sequence of functions $(y_{\ell})_{\ell\in\mathbb{N}}$ weakly converges to $y$ in $L_{2}(d\psi,[0,T],\mathbb{R})$, for some $\psi\in AC([0,T],\mathbb{R})$. Let $(x_{\ell})_{\ell\in\mathbb{N}}$ be a sequence of absolutely continuous functions such that it converges uniformly to $x\in AC([0,T],\mathbb{R}^{n})$, and $\dot{x}_{\ell}(t)=\dot{\psi}(t)y_{\ell}(t)$, for each t\in\Gamma:=\bigl{\{}t\in[0,T]\,|\,x_{\ell},x and $\psi$ are differentiable at t\bigr{\}}. Then, it holds that $\dot{x}(t)=\dot{\psi}(t)y(t)$ for almost every $t\in\Gamma$.

**Proof. **Define the function $\xi:[0,T]\rightarrow\mathbb{R}^{n}$ by

[TABLE]

for $t\in[0,T]$. For every $\eta\in\mathbb{R}^{n}$, we have

[TABLE]

for all $\ell\in\mathbb{N}$ and

[TABLE]

Since $(y_{\ell})_{\ell\in\mathbb{N}}$ weakly converges to $y$, we have that \big{(}\left\langle\eta,x_{\ell}(t)\right\rangle\big{)}_{\ell\in\mathbb{N}} converges to $\left\langle\eta,\xi(t)\right\rangle$ for every $t\in[0,T]$ and every $\eta\in\mathbb{R}^{n}$. This means that \big{(}x_{\ell}(t)\big{)}_{\ell\in\mathbb{N}} converges to $\xi(t)$ for every $t\in[0,T]$. Hence, we see that $\xi(t)=x(t)$ for all $t\in[0,T]$ since \big{(}x_{\ell}\big{)}_{\ell\in\mathbb{N}} uniformly converges to $x$. Therefore, (4) yields

[TABLE]

In other words,

[TABLE]

for almost every $t\in\Gamma$. $\blacksquare$

For the next two statements, we recall that two measures are absolutely continuously equivalent if each one is absolutely continuous with respect to the other one.

Lemma 8**.**

Let $f_{\ell}:[0,T]\rightarrow\mathbb{R}$ be a sequence of functions with $\ell\in\mathbb{N}$ such that $|f_{\ell}(t)|\leqslant 1$ for all $\ell\in\mathbb{N}$ and $t\in[0,T]$. Suppose that the sequence $(f_{\ell})_{\ell\in\mathbb{N}}$ weakly converges to $f$ in $L_{2}(d\mu,[0,T],\mathbb{R})$ where $d\mu$ is absolutely continuously equivalent to Lebesgue measure. Then,

[TABLE]

for almost all $t\in[0,T]$.

**Proof. **Let $k\geqslant 1$ and define $g^{k}_{\ell}(t):=\sup_{q\geqslant k}f_{q}(t)-f_{\ell+k}(t)$. Note that $(g^{k}_{\ell})_{{\ell}\in\mathbb{N}}$ weakly converges in $L_{2}(d\mu,[0,T],\mathbb{R})$ to $g^{k}$ given by $g^{k}(t):=\sup_{q\geqslant k}f_{q}(t)-f(t)$. Since $g^{k}_{\ell}$ is nonnegative for all $\ell\in\mathbb{N}$ and $t\in[0,T]$, $g^{k}$ must be nonnegative for almost all $t\in[0,T]$. This means that $f(t)\leqslant\sup_{q\geqslant k}f_{q}(t)$ for almost all $t\in[0,T]$. Hence, $f(t)\leqslant\limsup_{\ell{\rightarrow}\infty}f_{\ell}(t)$ for almost all $t\in[0,T]$. Applying the same arguments to the sequence $(-f_{\ell})_{\ell\in\mathbb{N}}$, we can obtain $f(t)\geqslant\liminf_{\ell{\rightarrow}\infty}f_{\ell}(t)$ for almost all $t\in[0,T]$. $\blacksquare$

Lemma 9**.**

Let $y_{\ell}:[0,T]\rightarrow\mathbb{R}^{q}$ be a sequence of functions with $\ell\in\mathbb{N}$ such that $|y_{\ell}(t)|\leqslant 1$ for all $\ell\in\mathbb{N}$ and $t\in[0,T]$. Also let \big{(}S_{\ell}(t)\big{)}_{\ell\in\mathbb{N}} be a sequence of sets in $\mathbb{R}^{q}$ with $\ell\in\mathbb{N}$ and $t\in[0,T]$ such that $y_{\ell}(t)\in S_{\ell}(t)$ for all $\ell\in\mathbb{N}$ and $t\in[0,T]$. Suppose that $(y_{\ell})_{\ell\in\mathbb{N}}$ weakly converges to $y$ in $L_{2}(d\mu,[0,T],\mathbb{R}^{q})$ where $d\mu$ is absolutely continuously equivalent to Lebesgue measure. Then,

[TABLE]

for almost all $t\in[0,T]$.

**Proof. **Let S(t)=\operatorname{cl}\Big{(}\operatorname{conv}\big{(}\limsup_{\ell\rightarrow\infty}S_{\ell}(t)\big{)}\Big{)} for $t\in[0,T]$. It follows from (Rockafellar98, , Cor. 4.11) that $S(t)\neq\varnothing$ for each $t\in[0,T]$. Let $\Gamma=\{t\in[0,T]:y(t)\notin S(t)\}$. Define the function $z:[0,T]\rightarrow\mathbb{R}$ by

[TABLE]

for all $t\in\Gamma$ and $z(t)=0$ for all $t\in[0,T]\setminus\Gamma$. Note that $|y(t)-z(t)|>0$ for all $t\in\Gamma$. Also, we have $z\in L_{\infty}([0,T],\mathbb{R}^{n})$ since $S(t)$ contains an element in the unit ball of $\mathbb{R}^{n}$ for all $t\in[0,T]$. Now, define functions $a:[0,T]\rightarrow\mathbb{R}^{n}$ and $b:[0,T]\rightarrow\mathbb{R}$ by

[TABLE]

for all $t\in\Gamma$ and $a(t)=0$, $b(t)=0$ for all $t\in[0,T]\setminus\Gamma$. For all $t\in\Gamma$, the hyperplane $\mathcal{H}_{t}=\{\eta:\left\langle a(t),\eta\right\rangle=b(t)\}$ strictly separates the set $S(t)$ and the point $y(t)$, that is

[TABLE]

for all $z\in S(t)$ (see e.g. (bertsekas, , Prop. 1.5.3)). Note that $a\in L_{\infty}([0,T],\mathbb{R}^{n})$, and since $y\in L_{2}(d\mu,[0,T],\mathbb{R}^{n})$ and $z\in L_{\infty}([0,T],\mathbb{R}^{n})$, it follows that $b\in L_{2}(d\mu,[0,T],\mathbb{R})$. Therefore, the function $t\mapsto\left\langle a(t),w(t)\right\rangle$ belongs to $L_{2}(d\mu,[0,T],\mathbb{R})$ for every $w\in L_{2}(d\mu,[0,T],\mathbb{R}^{n})$. For each $\ell\in\mathbb{N}$, define $\zeta_{\ell}:[0,T]\rightarrow\mathbb{R}$ with $\zeta_{\ell}(t)=\left\langle a(t),y_{\ell}(t)\right\rangle$ for all $t\in\Gamma$ and $\zeta_{\ell}(t)=0$ for all $t\in[0,T]\setminus\Gamma$. Then, we see that $(\zeta_{\ell})_{\ell\in\mathbb{N}}$ weakly converges to $\zeta$ given by $\zeta(t)=\left\langle a(t),y(t)\right\rangle$ for all $t\in\Gamma$ and $\zeta(t)=0$ for all $t\in[0,T]\setminus\Gamma$. From Lemma 8, we see that

[TABLE]

for almost all $t\in[0,T]$. Since limit inferior (superior) can be obtained as the limit of a subsequence, (8) implies that for almost all $t\in\Gamma$

[TABLE]

where $\underline{y}(t)$ and $\overline{y}(t)$ belong to $S(t)$. Together with the second inequality in (7), this yields

[TABLE]

for almost all $t\in\Gamma$. In view of the first inequality in (7), this means that $\Gamma$ is a zero measure set. As such, we can conclude that

[TABLE]

for almost all $t\in[0,T]$. $\blacksquare$

3 Differential inclusions with maximal monotone mappings

Our goal is to study the existence of solutions to the differential inclusion

[TABLE]

where $F(t,\cdot):\mathbb{R}^{n}\rightrightarrows\mathbb{R}^{n}$ is maximal monotone for all $t\in[0,T]$. We say that $x\in AC([0,T],\mathbb{R}^{n})$ is a solution of (10) if $x(t)\in\operatorname{dom}F(t,\cdot)$ and $x$ satisfies (10) for almost all $t\in[0,T]$.

3.1 Related frameworks and their limitations

Historically, the evolution inclusions given in (10) have been a subject of research in mathematical community in different eras. However, the solutions to such equations have been proposed under rather strict conditions. Here, we provide a brief of list of the main results that exist concerning the existence and uniqueness of solutions for such systems. We also refer the reader to (HuPapa00, , Chapter II) for an overview of the special case when $F(t,\cdot)$ is the subdifferential of a convex function.

Single-valued operators with fixed domain Kato67 : The earliest results on solutions of dynamical systems (10) with time-dependent maximal monotone relations were proposed in Kato67 . The author focused on the case where $F(t,\cdot):\mathbb{R}^{n}\to\mathbb{R}^{n}$ is single-valued and $F(\cdot,x)$ is Lipschitz continuous, uniformly in $x$. The major restriction imposed here is that

[TABLE]

Under these conditions, there exists a Lipschitz continuous $x:\mathbb{R}_{+}\to\mathbb{R}^{n}$ such that (10) holds for Lebesgue almost every $t\geqslant 0$, and $x(t)\in\operatorname{dom}F(t,\cdot)$ for each $t\geqslant 0$.

Static maximal monotone operators (à la Brézis) with additive inputs Brez73 : In the classical book Brez73 dealing with differential inclusions with maximal monotone operators, we can find results dealing with inclusions of the form

[TABLE]

where $A$ is maximal monotone, $\gamma>0$ is a scalar, and $u:[0,\infty)\to\mathbb{R}^{n}$ is absolutely continuous. For such systems $\operatorname{dom}F(t,\cdot)=\operatorname{dom}A$ for each $t\geqslant 0$, that is, the domain of the multivalued operator is independent of time.

Dissipative operators Pave87 : Building up on the work of Kato Kato67 and Brézis Brez73 , we find results on evolution equations built on convergence of certain discrete approximations in Pave87 . When the results appearing in this line of work are applied to the maximal monotone case given in (10), it turns out that such results also require the strong assumption (11), where the domain of the operator does not change with time (Pave87, , Chap. 1, Sec. 4).

Moreau’s sweeping process More77 : The first real contribution in the literature with time-dependent domains is observed in the seminal work of More77 . The systems studied here within the umbrella of sweeping processes comprise differential inclusions with a special conic structure. We introduce a set-valued mapping $S:\mathbb{R}_{+}\rightrightarrows\mathbb{R}^{n}$, and let $\mathcal{N}_{S(t)}(x)$ denote the normal cone to the set $S(t)$ at a point $x\in S(t)$. The proposed system class is then described as:

[TABLE]

Thus, for each $t$ and $x$, $F(t,x)$ is a closed convex cone described by the subdifferential of the indicator function of $S(t)$, and hence $F(t,\cdot)$ is maximal monotone. Here, we see that $\operatorname{dom}(F(t,\cdot))=S(t)$ and since $S$ is time-dependent, the domain is allowed to vary with time. To describe the regularity imposed on $F(\cdot,x)$ with respect to time, we consider the Hausdorff distance, and in the simplest instances, it is assumed that, for every $t_{1},t_{2}\geqslant 0$

[TABLE]

that is the Hausdorff distance between the domains of $F(t,\cdot)$ is bounded by a Lipschitz continuous function of time. Under these assumptions, there exists a unique solution to (12) which is Lipschitz continuous. Different variants of this framework were then derived depending on how the Hausdorff distance varies with time, or whether we can relax the convexity assumption on $S(t)$ while preserving some nice properties of the subdifferential of the indicator function for that set. In short, sweeping processes provide the first instance in the literature on inclusions with a particular kind of maximal monotone operators which depend on time, and the corresponding domain may vary.

Maximal monotone operators with time-dependent domain Vlad91 ; KunzMont97 : As a generalization of the sweeping process, Vladimirov Vlad91 studied evolution inclusions where time-dependent domains were considered, with the hypothesis that the set-valued mapping $F(t,\cdot)$ is just maximal monotone for each $t\geqslant 0$, without any further structural or geometrical assumption. However, a very strong regularity assumption was imposed with respect to the pseudo-distance given in Definition 2. In particular, the mapping $F(t,\cdot)$ is required to be uniformly continuous, that is, there exists a sequence of piecewise constant operators $F_{i}:[0,T]\times\mathbb{R}^{n}\rightrightarrows\mathbb{R}^{n}$ such that for each $t\in[0,T]$

[TABLE]

Kunze and Monteiro-Marques KunzMont97 then generalized this line of work to consider systems where the regularity with respect to time can be relaxed, so that the pseudo-distance between $F(t_{1},\cdot)$ and $F(t_{2},\cdot)$ is bounded by $|\mu(t_{1})-\mu(t_{2})|$ for some function of bounded variation $\mu:[0,T]\to\mathbb{R}$. Certain results developed in the context of sweeping processes are thus covered within this framework. The work started by Vladimirov, and later generalized to some extent by Kunze and Monteiro-Marques, indeed is an attempt to deal with differential inclusions with most general time-dependent maximal monotone operators. However, they impose very strong assumptions in deriving their results which make their applicability somewhat restrictive. Indeed, as we show in the next section, strong continuity assumption is not necessary.

3.2 Motivation

A primary motivation for looking at inclusions of type (10) comes from differential equations where certain variables are related by a maximal monotone operator. In particular, consider systems described by

[TABLE]

where $x\in\mathbb{R}^{n}$, $(z,w)\in\mathbb{R}^{m}\times\mathbb{R}^{m}$, $v\in\mathbb{R}^{m}$, the matrices $(A,B,C,D)$ have appropriate dimensions and $\mathcal{M}:\mathbb{R}^{m}\rightrightarrows\mathbb{R}^{m}$ is a maximal monotone operator.

Systems of the form (13) can be alternatively described by (10) where

[TABLE]

By invoking (CamlSchu16, , Theorem 2), one can show that $F(t,\cdot)$ is maximal monotone for each $t\in\mathbb{R}_{+}$ under certain assumptions. Regularity with respect to time is critical here. In the works of Vlad91 ; KunzMont97 , existence and uniqueness of solutions to (10) is established under the assumption of absolute continuity in the sense of Definition 4. However, the mapping $t\mapsto F(t,\cdot)$ defined by (14) does not, in general, enjoy absolute continuity with respect to pseudo-distance even if $v$ is absolutely continuous. This is seen in the following example.

Example 10**.**

Consider a system of the form (13) where $n=1$, $m=2$,

[TABLE]

and $\mathcal{M}:\mathbb{R}^{2}\rightrightarrows\mathbb{R}^{2}$ is the set-valued mapping given by $\mathcal{M}(\zeta)=\{\eta:\eta\geqslant 0,\,\zeta\leqslant 0,\text{ and }\left\langle\eta,\zeta\right\rangle=0\}$. By invoking (CamlSchu16, , Theorem 2), it can be verified that the corresponding set-valued mapping $F(t,\cdot)$ as defined in (14) is maximal monotone for each $t$. Let $v:[0,T]\to\mathbb{R}^{2}$ be an absolutely continuous function such that for some $t_{1},t_{2}\in[0,T]$, we have

[TABLE]

Let $F_{i}:=F(t_{i},\cdot)$ with $i=1,2$. It can be verified that

[TABLE]

for any $\rho\geqslant 0$. From Definition 2, we get

[TABLE]

Since the righthand side is not bounded, we can conclude that set-valued mapping $F(t,\cdot)$ is not absolutely continuous in the sense of Definition 4. However, existence and uniqueness of solutions for this example would follow from our main results. Indeed, this example satisfies the hypothesis of Theorem 26.

4 Main results

The main goal of this paper is to investigate conditions (weaker than those of Vlad91 ; KunzMont97 ) that guarantee existence of solutions to (10). The uniqueness of solutions for a fixed initial condition follows easily from the maximal monotone property of the right-hand side of (10).

4.1 Existence and uniqueness of solutions to (10)

To state the main result of our paper, we introduce the following assumptions, where $T>0$ and $\mathbb{B}^{n}(r):=\{x\in\mathbb{R}^{n}:|x|\leqslant r\}$.

1. (A1)

For each $t\in[0,T]$, the operator $F(t,\cdot)$ is maximal monotone. 2. (A2)

The set-valued mapping $t\mapsto\operatorname{graph}F(t,\cdot)$ is outer semicontinuous from right on $[0,T]$. 3. (A3)

There exists a nondecreasing function $\varphi\in AC([0,T],\mathbb{R})$ such that

[TABLE]

1. (LG)

There exists $\sigma\in L_{1}([0,T],\mathbb{R}_{+})$ such that

[TABLE]

for all $x\in\operatorname{dom}{F(t,\cdot)}$ with $t\in[0,T]$.

1. (YB)

There exist a nondecreasing function $\theta:[0,T]\rightarrow\mathbb{R}$, and a scalar $\Lambda>0$, such that for every $0<\lambda<\Lambda$,

[TABLE]

for all $s\text{ with }0\leqslant s\leqslant T-\lambda$, and $x\in\operatorname{dom}F(s,\cdot)$ where $Y_{\lambda}(t,\cdot)$ denotes the Yosida approximation of $F(t,\cdot)$.

The result on existence and uniqueness of solutions now follows.

Theorem 11**.**

Consider the system (10) and assume that (A1)–(A3) hold. If the linear growth assumption (LG) holds, then there exists a unique solution of (10) for every x_{0}\in\operatorname{cl}\big{(}\operatorname{dom}F(0,\cdot)\big{)}. In case the Yosida approximation bound assumption (YB) holds, then there exists a unique solution of (10) for every $x_{0}\in\operatorname{dom}F(0,\cdot)$.

Remark 12**.**

The statement (YB) in the description of Theorem 11 can be checked without the explicit knowledge of the least norm element $|F^{0}(s,x)|$ for each $s\geq 0$. In particular, consider the following statement:

1. (YB’)

There exist a nondecreasing function $\theta:[0,T]\rightarrow\mathbb{R}$, and $\Lambda>0$, such that for every $0<\lambda<\Lambda$,

[TABLE]

for all $s\text{ with }0\leqslant s\leqslant T-\lambda$, and $x\in\operatorname{dom}F(s,\cdot)$.

Then, because of Proposition 1, item (vi), (YB’) implies (YB).

Remark 13**.**

If we consider the static maximal monotone operator $F:\mathbb{R}^{n}\rightrightarrows\mathbb{R}^{n}$, and the differential inclusion $\dot{x}\in-F(x)$, then the classical result on existence of solutions of such differential inclusions, such as (Brez73, , Theorem 3.1), follows directly from Theorem 11. In that case, conditions (A2) and (A3) hold trivially, because the operator is time-independent and one can check that (YB) holds with a constant function $\theta$, because of Proposition 1, item (vi).

4.2 Relevance of Theorem 11

In what follows we will show how the results of KunzMont97 as well as the results on sweeping processes can be recovered from Theorem 11.

4.2.1 Recovering the results of KunzMont97

As recalled in Sect. 3.1, the results in KunzMont97 imposed continuity with respect to the pseudo-distance introduced in Definition 2. We claim that if the mapping $t\mapsto F(t,\cdot):\mathbb{R}^{n}\rightrightarrows\mathbb{R}^{n}$ is absolutely continuous on $[0,T]$ then the set-valued mapping $t\mapsto\operatorname{graph}F(t,\cdot)$ is outer semicontinuous from right on $[0,T]$. To see this, let $(t_{\ell},x_{\ell},y_{\ell})_{\ell\in N}\subseteq[0,T]\times\mathbb{R}^{n}\times\mathbb{R}^{n}$ be a sequence such that $y_{\ell}\in F(t_{\ell},x_{\ell})$, $t_{\ell}\geq t$, and $\lim_{\ell\uparrow\infty}(t_{\ell},x_{\ell},y_{\ell})=(t,x,y)$ for some $t\in[0,T]$, $x,y\in\mathbb{R}^{n}$. What needs to be proven is that $y\in F(t,x)$. To see this, let $(\eta,\zeta)\in\mathbb{R}^{n+n}$ be such that $\zeta\in F(t,\eta)$. From absolute continuity of $t\mapsto F(t,\cdot)$, we have

[TABLE]

By letting $\ell$ tend to infinity, we obtain

[TABLE]

since $r$ is continuous. This means that $y\in F(t,x)$ as $F(t,\cdot)$ is maximal monotone. Another hypothesis required by the results of KunzMont97 is a linear growth condition that coincides with (LG).

4.2.2 Inclusions with subdifferentials of convex functions

Another closely related class of dynamical systems appearing in the literature is given by

[TABLE]

where $g(t,\cdot):\mathbb{R}^{n}\to\mathbb{R}\cup\{\infty\}$ is a convex, lower semicontinuous, and proper for each $t\geq 0$, and $\partial g(t,\cdot)$ denotes its subdifferential.

Corollary 14**.**

Consider the differential inclusion (15), and assume that the following hold:

• •

There exists a nondecreasing function $\varphi\in AC([0,T],\mathbb{R})$ such that

[TABLE]

• •

There exist a nondecreasing function $\theta:[0,T]\rightarrow\mathbb{R}$, and a scalar $\Lambda>0$, such that for every $0<\lambda<\Lambda$,

[TABLE]

for all $s\text{ with }0\leqslant s\leqslant T-\lambda$, and $x\in\operatorname{dom}\partial g(s,\cdot)$, where $g_{\lambda}(t,\cdot)$ denotes the continuously differentiable Yosida approximation of $g(t,\cdot)$, defined as

[TABLE]

Then, the differential inclusion (15) admits a unique solution.

The corollary, stated above, is a direct consequence of Theorem 11. Under the assumptions on $g(t,\cdot)$, (A1) holds. Moreover, Yosida approximation of the operator $\partial g$ can be expressed directly in terms of the Yosida approximation of the function $g$ itself. We also have the following result which ensures that, in such cases, (A3) implies (A2).

Lemma 15**.**

If, for each $t\geq 0$, $F(t,\cdot)=\partial g(t,\cdot)$ for some $g(t,\cdot)$ convex, lower semicontinuous and proper. Then, (A3) implies (A2).

**Proof. **Let $(t_{\ell},x_{\ell},y_{\ell})_{\ell\in\mathbb{N}}\subseteq[0,T]\times\mathbb{R}^{n}\times\mathbb{R}^{n}$ be a sequence such that $y_{\ell}\in\partial g(t_{\ell},x_{\ell})$ and moreover, $\lim_{\ell\uparrow\infty}(t_{\ell},x_{\ell},y_{\ell})=(t,x,y)$, $t_{\ell}\geqslant t$, for some $t\in[0,T]$, $x,y\in\mathbb{R}^{n}$. We now show that $y\in\partial g(t,x)$. For that, let $(\eta,\zeta)\in\mathbb{R}^{n+n}$ be such that $\zeta\in\partial g(t,\eta)$, and let $x_{\ell,\eta}^{\star}:=\operatorname{argmin}_{z\in\operatorname{dom}\partial g(t_{\ell},\cdot)}|z-\eta|$. Note that, $\operatorname{dom}\partial g(t,\cdot)$ is closed, convex and nonempty for each $t\geq 0$, and hence $x_{\ell,\eta}^{\star}$ is unique, and because of (A3), $x_{\ell,\eta}^{\star}\to\eta$ as $\ell\to\infty$. We, therefore, have

[TABLE]

By letting $\ell$ tend to infinity, and invoking (A3), we obtain

[TABLE]

*which implies $y\in\partial g(t,x)$, and hence assumption (A2) holds. *$\blacksquare$

Example 1.

As an illustration of Corollary 14, consider the function $g:\mathbb{R}_{+}\times\mathbb{R}\to\mathbb{R}_{+}$, defined by, $g(t,x)=t|x|$. The first item holds since the domain of $\partial g$ does not change with time. To check the second condition, Yosida approximation of the function $g(t,\cdot)$, for each $t\geq 0$, is described by:

[TABLE]

As a result, Yosida approximation of the operator $\partial g(t,\cdot)$, for each $t\geq 0$, is given by:

[TABLE]

and hence it follows that for each $t\geq s\geq 0$, and each $x\in\mathbb{R}$, we have

[TABLE]

and hence (YB’) holds with linear function $\theta(s)=s$.

4.2.3 Recovering the special case of sweeping processes

Sweeping process is a special case of (10), where $F(t,x)=\mathcal{N}_{S(t)}(x)$, with $S:[0,T]\rightrightarrows\mathbb{R}^{n}$ being a closed convex-valued mapping. The normal cone operator is by definition maximal monotone for each $t\in[0,T]$. Moreover, in this particular case, assumption (LG) is trivially satisfied since $\mathcal{N}_{S(t)}(x)$ is a cone and thus $0\in\mathcal{N}_{S(t)}(x)$, for each $t\in[0,T]$ and $x\in\operatorname{dom}\,\mathcal{N}_{S(t)}(\cdot)=S(t)$. Assumption (A3) can be written directly in terms of the set-valued mapping $S$ as follows: There exists a nondecreasing function $\varphi\in AC([0,T],\mathbb{R})$ such that

[TABLE]

The aforementioned condition is equivalent to saying that the retraction of the set-valued mapping $S(\cdot)$ is absolutely continuous, and it was first proposed in More77 for existence of solutions to first order sweeping processes. One can invoke Lemma 15 to check that assumption (A3) implies (A2). Indeed, by taking $g(t,x)=\begin{cases}0&\text{if }x\in S(t),\\ +\infty&\text{if }x\not\in S(t)\end{cases}$, and observing that $\mathcal{N}_{S(t)}(\cdot)=\partial g{(t,\cdot)}$.

Remark 16**.**

In Section 4.2.2, and Section 4.2.3, we considered a specific class of systems where the maximal monotone operator is defined by the subdifferential of the convex function. In the literature, we also find other results for differential inclusions with subdifferentials of time-dependent convex functions, see for example (HuPapa00, , Page 197), Kubo88 , Papa90 , and Tols17 . In some of these works, conditions are stated in terms of the function $g(\cdot,\cdot)$ itself, whereas we work with assumptions that result in the bounds on the corresponding Yosida approximations of the subdifferentials. Some results from Kubo88 can be used to write the conditions in terms of the function $g(\cdot,\cdot)$ itself, but in general, it is not clear how our conditions on the subdifferential operator can be translated to the conditions on the function $g(\cdot,\cdot)$.

5 Proof of Theorem 11

We are basically concerned with the existence of the solution in this proof, as the uniqueness readily follows from assumption (A1). The proof of existence is based on constructing a sequence of approximate solutions and showing that this sequence converges to a function which satisfies the differential inclusion (10). This is formally done in following main steps:

• •

Discretizing (10)

• •

Obtaining bounds on discrete values

• •

Construction of a sequence of approximate solutions

• •

Studying the limit of the sequence

Each of these steps is carried out as a subsection in the sequel, and it is shown that the limit we thus obtain is indeed a solution to (10).

5.1 Discretization of (10)

We first begin with discretizing (10). Let

[TABLE]

be a partition of the interval $[0,T]$. Define

[TABLE]

for $k\in\{1,2,\ldots,K\}$. Note that $\sum_{k=1}^{K}h_{k}=T$. We define the size of the partition $\Delta$ by $K(\Delta)$ and the granularity by $|\Delta|=\max_{k\in\{1,2,\ldots,K\}}h_{k}$. For simplicity, we write $K=K(\Delta)$ when $\Delta$ is clear from the context.

Next, consider the discretization of (10) based on the partition $\Delta$ given by

[TABLE]

for $k\in\{0,1,\ldots,K-1\}$. Alternatively, we have

[TABLE]

This resolvent based alternative form, together with assumption (A1) and Proposition 1, guarantees that the discretization (16) is well-defined in the sense that there exist $x_{0},x_{1},\ldots,x_{K}$ satisfying (16) (and hence (17)). If $\operatorname{dom}(F(t,\cdot))$ is closed for each $t\geqslant 0$, then it follows from (17) that $x_{k+1}$ is a projection of $x_{k}$ on the set $\operatorname{dom}(F(t_{k+1},\cdot))$; Put simply, in case the domain of the multivalued function $F(t,\cdot)$ is closed, the sequence of points $x_{k}$ is obtained by applying the proximal operator associated with the mapping $h_{k}F(t_{k},\cdot)$, as defined in Rock76 , PariBoyd13 . We will use a certain interpolation between the points $x_{k}$ to get an approximate solution for the differential inclusion (10).

5.2 Bounds on $x_{k}$ values

We aim at establishing bounds on $x_{k}$ that are independent of the underlying partition $\Delta$. The bounds obtained under the condition (LG) and (YB) are different and appear in Lemma 17 and Lemma 18, respectively. We prefer using the same notation for the bounds of our interest in these lemmas, and when we use these bounds later on, the reference to these bounds will be clear from the context.

Lemma 17**.**

Under assumption (LG), there exists a sequence of partitions $\{\Delta_{\ell}\}$ with $|\Delta_{\ell}|\to 0$, such that for each of these partitions, we have

[TABLE]

for each $k\in\{1,2,\ldots,K\}$; the constant $\beta\geq 0$, and the function $\psi:[0,T]\rightarrow\mathbb{R}_{+}$ are defined as

[TABLE]

with $\varphi$ satisfy (A3), $\sigma$ satisfying (LG), and

[TABLE]

**Proof. **To obtain the bounds given in (18) and (19), we start by analyzing the sequence (17) for a fixed partition and introduce some simplified notation for the corresponding operators:

[TABLE]

It then follows from (17) that

[TABLE]

for all $k\in\{0,1,\ldots,K-1\}$, where we recall that $K$ is the size of the chosen partition.

To establish (18), we first introduce auxiliary points $\bar{x}_{k}$ given by $\bar{x}_{0}=x_{0}$ and

[TABLE]

for $k\in\{0,1,\ldots,K-1\}$. Clearly, we have

[TABLE]

for all $k\in\{0,1,\ldots,K\}$. Then, it follows from assumption (A3) and (2) that

[TABLE]

for all $k\in\{1,2,\ldots,K\}$. We thus obtain, for each $k\in\{1,2,\ldots,K\}$,

[TABLE]

where (29b) follows from the repeated application of (29a), and (29c) uses the fact that $\varphi$ is nondecreasing. An immediate consequence of (29c) is that

[TABLE]

for each $k\in\{0,1,\ldots,K\}$ where $\alpha$ satisfies (22).

The bounds given in (18) and (19) are now obtained from assumption (LG). From (27), we have that, for every $\varepsilon>0$ and each $k\in\{1,2,\dots,K\}$, there exists a point $\bar{x}_{k}^{\varepsilon}$ satisfying

[TABLE]

and therefore

[TABLE]

for all $k\in\{0,1,\ldots,K\}$ in view of (30) and the triangle inequality. Next, we introduce the sequence of points $\bar{y}_{k}^{\varepsilon}$ based on $\bar{x}_{k}^{\varepsilon}$ by

[TABLE]

for all $k\in\{0,1,\ldots,K\}$. Note that

[TABLE]

Now, it follows from (31) and Proposition 1 that

[TABLE]

To obtain a bound on the right-hand side of (36), we employ assumption (LG). Take any partition $\Delta_{\ell}$, and assume momentarily that $\sigma$ is a constant function.

By using assumption (LG) and (33), we get

[TABLE]

for all $k\in\{0,1,\ldots,K\}$ and $\varepsilon>0$. In view of (26) and (34), Proposition 1 implies that

[TABLE]

Hence, using (37) and (38), we obtain

[TABLE]

for all $k\in\{1,2,\ldots,K\}$ and $\varepsilon>0$. Letting $\varepsilon$ tend to zero in (39) results in

[TABLE]

for all $k\in\{1,2,\ldots,K\}$. Therefore, we have

[TABLE]

for all $k\in\{1,2,\ldots,K\}$. Since the constant $\sigma$ is finite, we can conclude from (30) and (45) that

[TABLE]

for all $k\in\{1,2,\ldots,K\}$ where $\beta$ satisfies (20). This establishes (18).

Next, we proceed to establish (19) using the bound in assumption (LG) and (18). To this end, we continue using the notation introduced in (24), and introduce a sequence of auxiliary points $\xi_{k}$

[TABLE]

for all $k\in\{1,2,\ldots,K\}$. Clearly, we have

[TABLE]

for all $k\in\{1,2,\ldots,K\}$. Therefore, for all $\varepsilon>0$ there exist points $\xi_{k}^{\varepsilon}$ satisfying

[TABLE]

It follows from (2) and assumption (A3) that

[TABLE]

for all $k\in\{1,2,\ldots,K\}$. In view of (46) and the fact that $\varphi$ is nondecreasing, this means that

[TABLE]

for all $k\in\{1,2,\ldots,K\}$. As such, we have

[TABLE]

for all $k\in\{1,2,\ldots,K\}$ where $\gamma$ satisfies (23). Now, define

[TABLE]

for all $k\in\{1,2,\ldots,K\}$. Note that

[TABLE]

where we used the fact that $|x_{k}-\zeta_{k}^{\varepsilon}|\leqslant|x_{k-1}-\xi_{k}^{\varepsilon}|$ due to the resolvent being nonexpansive. Since $\xi_{k}^{\varepsilon}\in\operatorname{dom}F_{k}$ due to (49a), we can invoke the bound on the least norm element of $F(t_{k},\xi_{k}^{\varepsilon})$ in assumption (LG) to obtain

[TABLE]

for all $k\in\{1,\ldots,K\}$ and $\varepsilon>0$. Together with (49b) and (50), (53) leads to

[TABLE]

by taking the limit as $\varepsilon$ tends to zero. This establishes the bound given in (19) for constant $\sigma$ and arbitrary partition.

To obtain the bounds (18) and (19) for $\sigma$ locally integrable, and not necessarily constant, one can invoke (LaksLeel81, , Lemma 3.3.1). To use this result, one must choose the partitions $\Delta_{\ell}$ carefully, and work with a piecewise constant approximation (in $L_{1}$-norm) of the function $\sigma$. $\blacksquare$

Lemma 18**.**

Under the assumption (YB), for any partition $\Delta$, with $|\Delta|<\Lambda$ we have

[TABLE]

for each $k\in\{1,2,\ldots,K\}$; the constant $\beta\geq 0$, and the function $\psi:[0,T]\rightarrow\mathbb{R}_{+}$ are defined as

[TABLE]

with $\alpha:=2(\theta(T)-\Theta(0))+F^{0}(0,x_{0})$.

**Proof. **Let $x_{0}\in\operatorname{dom}F(0,\cdot)$. Since $x_{k+1}=J_{k+1}(x_{k})\in\operatorname{dom}F_{k+1}$, we observe that

[TABLE]

On the other hand,

[TABLE]

and hence,

[TABLE]

We therefore have

[TABLE]

and, in particular, for each $k\in\mathbb{N}$,

[TABLE]

Plugging this bound in (58), and using the fact that $\theta$ is nondecreasing, we obtain

[TABLE]

Letting $\alpha:=2(\theta(T)-\theta(0))+F_{0}^{0}(x_{0})$, we get

[TABLE]

so that (55) holds with $\psi(s)=\alpha\,s$. By repeated application of the trying inequality, $|x_{k}|\leq|x_{k-1}|+|x_{k}-x_{k-1}|$, we get

[TABLE]

so that (54) holds with $\beta:=|x_{0}|+\alpha T$. $\blacksquare$

In subsequent analysis, we refer to $\beta$ and $\psi$, either from Lemma 17 or Lemma 18, depending upon whether we assume (LG) or (YB).

5.3 Construction of a sequence of approximate solutions

Based on the $x_{k}$ values, we construct a sequence of absolutely continuous (in time) functions which approximate the actual solution of the system. To this end, note that the function $\psi$ defined above is strictly increasing and absolutely continuous. Now, define the piecewise continuous function $x_{\Delta}$ as

[TABLE]

where $t\in[t_{k},t_{k+1}]$ and $k\in\{0,1,\ldots,K-1\}$. By definition, $x_{\Delta}$ is a continuous function and

[TABLE]

for all $k\in\{0,1,\ldots,K\}$. We will show that

[TABLE]

is the desired solution to the inclusion (10). An important intermediate step in studying the convergence of the sequence $x_{\Delta}$ is to obtain the following uniform bound.

Lemma 19**.**

Let $\underline{\tau}$ and $\overline{\tau}$ be such that $0\leqslant\underline{\tau}<\overline{\tau}\leqslant T$. For any partition $\Delta$, it holds that

[TABLE]

**Proof. **From the definition of $x_{\Delta}$ (59) for a fixed partition $\Delta$, there exist integers $q$ and $r$ with $q+1\leqslant r$ such that $t_{q}\leqslant\underline{\tau}<t_{q+1}$ and $t_{r-1}<\overline{\tau}\leqslant t_{r}$. If $q+1=r$, then we have

[TABLE]

In a similar fashion, if $q+1<r$ then we have

[TABLE]

Hence, (61) is established. $\blacksquare$

5.4 Limit of the sequence

The bounds established in the previous section allow us to study the limiting behaviour of the sequence $(x_{\Delta_{\ell}})_{\ell\in\mathbb{N}}$.

Lemma 20**.**

Consider a sequence of partitions $(\Delta_{\ell})_{\ell\in\mathbb{N}}$ with $|\Delta_{\ell}|\rightarrow 0$ as $\ell$ tends to infinity. The sequence $(x_{\Delta_{\ell}})_{\ell\in\mathbb{N}}$ is equicontinuous.

**Proof. **Note first that $\psi$ introduced in (21) (or, (57)) is uniformly continuous on the compact interval $[0,T]$ as it is absolutely continuous on the same interval. Therefore, for any $\varepsilon>0$ there exists a positive number $\delta>0$ such that

[TABLE]

for all $\underline{\tau},\overline{\tau}\in[0,T]$ such that $|\overline{\tau}-\underline{\tau}|<\delta$. In view of (61), we have

[TABLE]

for all $\ell\in\mathbb{N}$ and $\underline{\tau},\overline{\tau}\in[0,T]$ such that $|\overline{\tau}-\underline{\tau}|<\delta$. Consequently, the sequence $(x_{\Delta_{\ell}})_{\ell\in\mathbb{N}}$ is equicontinuous. $\blacksquare$

Let $(\Delta_{\ell})_{\ell\in\mathbb{N}}$ be a sequence of partitions with $|\Delta_{\ell}|\rightarrow 0$ as $\ell$ tends to infinity. Since the sequence $(x_{\Delta_{\ell}})_{\ell\in\mathbb{N}}$ is also uniformly bounded in view of Lemma 17 (as well as Lemma 18), Theorem 5 (Arzelá-Ascoli theorem) implies that it converges uniformly to a continuous function $x$ on a subsequence, say $N\in\mathcal{N}_{\infty}^{\#}$. We claim that $x$ is absolutely continuous. To see this, let $\underline{\tau},\overline{\tau}\in[0,T]$ with $\underline{\tau}\leqslant\overline{\tau}$ and note that

[TABLE]

where the first inequality follows from the triangle inequality, the second from (61), and the third by taking the limit on the convergent subsequence $N$. Thus, absolute continuity of $x$ follows from absolute continuity of the function $\psi$.

Now, we want to show that $x$ is a solution of (10), that is

[TABLE]

for almost all $t\in[0,T]$.

Let $\Gamma\subseteq[0,T]$ be defined by $\Gamma=\{t\in(0,T):\psi\text{ and }x\text{ are both }\text{differentiable at }t\text{ and }t\not\in\cup_{\ell\in N}\Delta_{\ell}\}$. Since $\psi$ and $x$ are both absolutely continuous and $\cup_{\ell\in N}\Delta_{\ell}$ is countable, it is enough to show (65) for almost all $t\in\Gamma$.

For a partition $\Delta$, define

[TABLE]

for $t\in(t_{k},t_{k+1})$ and $y_{\Delta}(t_{k})=0$ for $t_{k}\in\Delta$.

From (59), we see that

[TABLE]

for all $t\in\Gamma$.

In view of (59) and Lemma 19, we see that $|y_{\Delta_{\ell}}|_{L_{\infty}}\leqslant 1$ for all $\ell$. Therefore, the sequence $(y_{\Delta_{\ell}})_{\ell\in N}$ is contained in the closed ball with radius $\sqrt{\psi(T)-\psi(0)}$ of the Hilbert space $L_{2}(d\psi,[0,T],\mathbb{R}^{n})$. As such, there exists a subsequence $N^{\prime}$ of $N$ such that $(y_{\Delta_{\ell}})_{\ell\in N^{\prime}}$ converges to $y$ weakly in $L_{2}(d\psi,[0,T],\mathbb{R}^{n})$. It then follows from Lemma 7 that

[TABLE]

for almost all $t\in\Gamma$.

Now, let $t^{*}\in\Gamma$. Then, for every $\ell\in N$, there must exist $k_{\ell}\in\{1,2,\ldots,K(\Delta_{\ell})\}$ with the property that $t_{k_{\ell}}<t^{*}<t_{k_{\ell}+1}$. Note that $\lim_{\ell\uparrow\infty}t_{k_{\ell}}=\lim_{\ell\uparrow\infty}t_{k_{\ell}+1}=t^{*}$ since $|\Delta_{\ell}|$ converges to zero as $\ell$ tends to infinity. By construction, we have

[TABLE]

Equivalently, we have

[TABLE]

Let $S_{\ell}(t^{*}):=-\frac{t_{k_{\ell}+1}-t_{k_{\ell}}}{\psi(t_{k_{\ell}+1})-\psi(t_{k_{\ell}})}F(t_{k_{\ell}+1},x_{\Delta_{\ell}}(t_{k_{\ell}+1}))$. From (69), we have that $y_{\ell}(t^{*})\in S_{\ell}(t^{*})$. We now invoke Lemma 9 and observe that $y(t^{*})\in\operatorname{cl}\left(\operatorname{conv}\left(\limsup_{\ell\to\infty}S_{\ell}(t^{*})\right)\right)$. Due to the outer-semicontinuity assumption, we have $\limsup_{\ell\to\infty}F(t_{k_{\ell}+1},x_{\Delta_{\ell}}(t_{k_{\ell}+1}))\subseteq F(t^{*},x(t^{*}))$. The set $F(t^{*},x(t^{*}))$ is closed and convex because of the maximal monotonicity property, and hence

[TABLE]

Since $\dot{\psi}(t^{*})\geqslant 1$, we get

[TABLE]

for each $t^{*}\in\Gamma$. $\blacksquare$

6 Extensions

In this section, we extend the results of Theorem 11. First, we consider non-autonomous differential inclusions of the form

[TABLE]

where $F(t,\cdot):\mathbb{R}^{n}\rightrightarrows\mathbb{R}^{n}$ is maximal monotone for all $t\geqslant 0$ and $u\in L_{1}([0,T],\mathbb{R}^{n})$. We begin with the following observation.

Lemma 21**.**

Consider the system (70). Let $u\in L_{1}([0,T],\mathbb{R}^{n})$ and $G$ be the set-valued mapping defined by G(t,\xi):=F\big{(}t,\xi+\int_{0}^{t}{u(\tau)\,d\tau}\big{)}. Then, the differential inclusion (70) admits a solution $x$ if and only if the differential inclusion

[TABLE]

admits a solution $\xi$.

**Proof. **For the ‘only if’ part, suppose that $x$ is a solution of (70). Define

[TABLE]

for all $t\geqslant 0$. Note that $\xi(0)=x_{0}$ and

[TABLE]

The ‘if’ part follows reversing the arguments. $\blacksquare$

Theorem 22**.**

Consider system (70) and suppose that $u\in L_{1}([0,T],\mathbb{R}^{n})$ and $F(t,\cdot)$ satisfies assumptions (A1), (A2), and (A3). If (LG) holds, then for every x_{0}\in\operatorname{cl}\big{(}\operatorname{dom}F(0,\cdot)\big{)}, there exists a unique solution $x\in AC([0,T],\mathbb{R}^{n})$. If (YB) holds, then for every $x_{0}\in\operatorname{dom}F(0,\cdot)$, there exists a unique solution $x\in AC([0,T],\mathbb{R}^{n})$.

**Proof. **In view of Lemma 21 and Theorem 11, it is enough to show that the time-dependent set-valued map $G$ defined by G(t,x)=F\big{(}t,x+\Phi(t)\big{)} with $\Phi(t)=\int_{0}^{t}{u(\tau)\,d\tau}$ satisfies assumptions (A1), (A2), (A3), and in addition, $G$ satisfies (LG) (resp. (YB)) if $F$ satisfies (LG) (resp. (YB)).

(A1): Since maximal monotonicity is preserved under translations (see e.g. (Rockafellar98, , Thm. 12.43)), $G$ satisfies assumption (A1) whenever $F$ satisfies it.

(A2): Note that $\operatorname{graph}G(t,\cdot)=\operatorname{graph}F(t,\cdot)-\{\Phi(t)\}\times\{0\}$. Since the set-valued mapping $t\mapsto\operatorname{graph}F(t,\cdot)$ is outer semicontinuous on $[0,T]$ by assumption and $\Phi$ is absolutely continuous, $t\mapsto\operatorname{graph}G(t,\cdot)$ is outer semicontinuous on $[0,T]$.

(A3): Note that $\operatorname{dom}G(t,\cdot)=\operatorname{dom}F(t,\cdot)-\Phi(t)$ for all $t\in[0,T]$. Therefore, we have

[TABLE]

for all $s,\,t$ with $0\leqslant s\leqslant t\leqslant T$ since $F$ satisfies assumption (A3). Note that

[TABLE]

for all $s,\,t$ with $0\leqslant s\leqslant t\leqslant T$ where

[TABLE]

Therefore, $G$ satisfies assumption (A3).

(LG): Since $F$ satisfies assumption (LG), we observe that

[TABLE]

for all $x\in\operatorname{dom}{G(t,\cdot)}$ with $t\in[0,T]$. Note that

[TABLE]

Together with (73), this results in

[TABLE]

for all $t\in[0,T]$ and for all $x\in\operatorname{dom}G(t,\cdot)$. Since the function $t\mapsto\sigma(t)$ is integrable and $t\mapsto 1+|\Phi(t)|$ is continuous, their product is integrable. Consequently, $G$ satisfies assumption (LG).

(YB): Let $Y_{\lambda}^{F}(t,\cdot)$ and $Y_{\lambda}^{G}(t,\cdot)$ denote the Yosida approximation of $F(t,\cdot)$ and $G(t,\cdot)$ respectively. With $F$ satisfying (YB), and observing that $x\in\operatorname{dom}G(s,\cdot)$ if and only if $x+\Phi(s)\in\operatorname{dom}F(s,\cdot)$, we get, for all $0\leq s<T$, $\lambda<\Lambda$,

[TABLE]

where $\overline{\theta}(s)=M_{u}+\theta(s)$ for some $M_{u}>0$. In writing down the second inequality, we used the fact that $Y_{\lambda}^{G}(t,x)=Y_{\lambda}^{F}(t,x+\Phi(t))$, and moreover, $\lambda^{-1}$-Lipschitzian property of $Y_{\lambda}^{G}(t,\cdot)$ along with integrability of $u$ yields

[TABLE]

for some $M_{u}>0$. $\blacksquare$

Now, we turn our attention to differential inclusions of the form

[TABLE]

where $F(t,\cdot):\mathbb{R}^{n}\rightrightarrows\mathbb{R}^{n}$ is maximal monotone for all $t\geqslant 0$, $f:\mathbb{R}^{n}\rightarrow\mathbb{R}^{n}$ is a function and $u\in L_{1}([0,T],\mathbb{R}^{n})$. Based on Theorem 22, we present the following existence and uniqueness result.

Theorem 23**.**

Consider system (76) and suppose that $u\in L_{1}([0,T],\mathbb{R}^{n})$, $f:\mathbb{R}^{n}\rightarrow\mathbb{R}^{n}$ is a Lipschitz continuous function, and $F(t,\cdot)$ satisfies assumptions (A1), (A2), and (A3). If (LG) holds, then for every x_{0}\in\operatorname{cl}\big{(}\operatorname{dom}F(0,\cdot)\big{)}, there exists a unique solution $x\in AC([0,T],\mathbb{R}^{n})$. If (YB) holds, then for every $x_{0}\in\operatorname{dom}F(0,\cdot)$, there exists a unique solution $x\in AC([0,T],\mathbb{R}^{n})$.

**Proof. **Let x_{0}\in\operatorname{cl}\big{(}\operatorname{dom}F(0,\cdot)\big{)} and let $x_{0}(t)=x_{0}$ for all $t\in[0,T]$. It follows from Theorem 22 that for each integer $\ell$ with $\ell\geqslant 1$ there exists a unique absolutely continuous function $x_{\ell+1}:[0,T]\rightarrow\mathbb{R}^{n}$ such that $x_{\ell+1}(0)=x_{0}$, $x_{\ell+1}(t)\in\operatorname{dom}F(t,\cdot)$ and the differential inclusion

[TABLE]

holds for almost all $t\in[0,T]$. In the rest of the proof, we will construct a solution of (76) by showing that the sequence $\{x_{\ell}(\tau)\}_{\ell\in\mathbb{N}}$ is a Cauchy sequence that converges to an absolutely continuous function which satisfies (76).

Step 1: The sequence $\{x_{\ell}(\tau)\}_{\ell\in\mathbb{N}}$ is Cauchy. By using (A1) and Lipschitzness of $f$, we see that

[TABLE]

for almost all $t\in[0,T]$ where $L_{f}$ is the Lipschitz constant of $f$. By integrating both sides from [math] to $\tau\in[0,T]$, we obtain

[TABLE]

Application of (Brez73, , Lemma A.5, p. 157) results in

[TABLE]

for all $\tau\in[0,T]$. Hence, we get

[TABLE]

for all $\tau\in[0,T]$ where $|\cdot|_{L_{\infty}}$ denotes the sup norm. Consequently, $(x_{\ell})_{\ell\in\mathbb{N}}$ converges uniformly on $[0,T]$ to a function $x$.

Step 2: The function $x$ belongs to $AC([0,T],\mathbb{R}^{n})$. For the moment, suppose that there exist an integer $L$ and a nondecreasing function $\widehat{\psi}:AC([0,T],\mathbb{R})$ such that

[TABLE]

for all $\ell\geqslant L$ and for all $\underline{\tau},\overline{\tau}$ with $0\leqslant\underline{\tau}<\overline{\tau}\leqslant T$. This would mean that we have

[TABLE]

for all $\underline{\tau},\overline{\tau}$ with $0\leqslant\underline{\tau}<\overline{\tau}\leqslant T$ where the first inequality follows from the triangle inequality, the second from (78) for all $\ell\geqslant L$, and the third by taking the limit as $\ell$ tends to infinity. Thus, absolute continuity of $x$ follows from absolute continuity of the function $\widehat{\psi}$.

To prove (78), we first observe that the triangle inequality and (77) result in

[TABLE]

for some $C_{1}>0$, and each $\ell\geqslant 1$, $\tau\in[0,T]$. Hence, in particular, $|x_{\ell}(\cdot)|$ is uniformly bounded for each $\ell\geqslant 1$. The bound in (78) is obtained differently in case if $F(\cdot,\cdot)$ satisfies (LG) or (YB).

If $F$ satisfies (LG): Similar to function $\overline{\varphi}$ in (72), let us introduce the function $\varphi_{\ell}$ as,

[TABLE]

so that, the Lipschitz continuity of $f$, $|f(z)|\leq|f(0)|+L_{f}|z|$, yields

[TABLE]

By introducing the function $\psi_{\ell}$, similar to (21), as

[TABLE]

and by letting $g(s):=s+(1+\gamma)\int_{0}^{s}\sigma(\tau)\,d\tau+2\int_{0}^{s}|u(\tau)|\,d\tau$, we get

[TABLE]

It follows from (64) that $|x_{\ell}(\overline{\tau})-x_{\ell}(\underline{\tau})|\leq\psi_{\ell}(\overline{\tau})-\psi_{\ell}(\underline{\tau})$. Thus we get

[TABLE]

for each $\ell\geqslant 1$, and hence (78) follows with $\widehat{\psi}(s):=g(s)+2\,L_{f}C_{1}s$, which is clearly an absolutely continuous function.

If $F$ satisfies (YB): In that case, the function $\psi_{\ell}$ is obtained from (57), so that $\psi_{\ell}(s):=\overline{\alpha}\,s$. The constant $\overline{\alpha}:=2(\overline{\theta}(T)-\overline{\theta}(0))$, where $\overline{\theta}$ is defined as in (75), that is, $\overline{\theta}(s):=M_{u}+C_{1}+\theta(s)$, for some $C_{1},M_{u}\geq 0$, which is independent of $\ell$. Hence, we take $\widehat{\psi}(s)=\overline{\alpha}\,s$, which is clearly absolutely continuous.

Step 3: The function $x$ satisfies (76). In view of Lemma 21, $x_{\ell}$ is a solution to the differential inclusion

[TABLE]

for all $\ell\geqslant 1$ where G_{\ell}(t,\xi)=F(t,\xi+\Phi_{\ell}(t)\big{)} with \Phi_{\ell}(t)=\int_{0}^{t}f\big{(}x_{\ell-1}(\tau)\big{)}+u(\tau)\,d\tau. Since $\operatorname{graph}G_{\ell}(t,\cdot)=\operatorname{graph}F(t,\cdot)-\{\Phi_{\ell}(t)\}\times\{0\}$.

Let us introduce the sequence $y_{\ell}:=\frac{dx_{\ell}}{d\widehat{\psi}}$, so that

[TABLE]

Because of the bound (78), $|y_{\ell}|_{L^{\infty}}\leq 1$, and there exists a subsequence $N$ such that $(y_{\ell})_{\ell\in N}$ converges to $y$ weakly in $L_{2}(d\widehat{\psi},[0,T],\mathbb{R}^{n})$. It then follows from Lemma 7 that

[TABLE]

for almost all $t\in\Gamma:=\{t\in[0,T]\,|\,x_{\ell},\ell\geqslant L,x,\text{ and }\psi\text{ are differentiable at }t\}$.

Let $t^{*}\in\Gamma$. By construction, we have

[TABLE]

or equivalently,

[TABLE]

In other words, $y_{\ell}(t^{*})$ belongs to the convex set $\frac{-1}{\dot{\widehat{\psi}}(t^{*})}G_{\ell}(t^{*},x_{\ell}(t^{*}))$. Using Lemma 9 with $S_{\ell}(t)=\frac{-1}{\dot{\widehat{\psi}}(t)}G_{\ell}(t,x_{\ell}(t))$, and recalling that $|\dot{\widehat{\psi}}(t)|\geqslant 1$ for each $t\in[0,T]$, we see that

[TABLE]

and the same holds for almost every $t^{*}\in[0,T]$. $\blacksquare$

7 Linear systems and maximal monotone relations

A particularly interesting class of time-dependent maximal monotone mappings arises from the interconnection of linear passive systems with maximal monotone relations. To formalize this class of systems, consider the linear system

[TABLE]

for some set-valued map $\mathcal{M}:\mathbb{R}^{m}\rightrightarrows\mathbb{R}^{m}$.

By solving $z$ from the relations (80b), (80c), and substituting in (80a), we obtain the differential inclusion

[TABLE]

where

[TABLE]

and

[TABLE]

The rest of this section is devoted to developing conditions under which the time-dependent set-valued mapping $H(t,\cdot)$ satisfies the hypotheses of Theorem 11. To establish such conditions, we first introduce the notion of passivity of a linear system.

A linear system $\Sigma(A,B,C,D)$

[TABLE]

is said to be passive, if there exists a nonnegative-valued storage function $V:\mathbb{R}^{n}\rightarrow\mathbb{R}_{+}$ such that the so-called dissipation inequality

[TABLE]

holds for all $t_{1},t_{2}$ with $t_{1}<t_{2}$ and for all trajectories $(z,x,w)\in L_{2}([t_{1},t_{2}],\mathbb{R}^{m})\times AC([t_{1},t_{2}],\mathbb{R}^{n})\times L_{2}([t_{1},t_{2}],\mathbb{R}^{m})$ of the system (84).

The classical Kalman-Yakubovich-Popov lemma states that the system (84) is passive if, and only if, the linear matrix inequalities

[TABLE]

admits a solution $K$. Moreover, $V(x)=\frac{1}{2}x^{\top}Kx$ defines a storage function in case $K$ is a solution of the linear matrix inequalities (86).

In the following proposition, we summarize some of the consequences of passivity that will be used later.

Proposition 1 ((camlibel:14, , Lem. 1)).

If $\Sigma(A,B,C,D)$ is passive with the storage function $x\mapsto\frac{1}{2}x^{\top}Kx$ then the following statements hold:

1. i.

$D$* is positive semi-definite.* 2. ii.

$(KB-C^{\top})\ker(D+D^{\top})=\{0\}$.

The following theorem states conditions that guarantee the hypotheses of Theorem 22 for the time-dependent set-valued mapping $H$ as defined in (82).

Theorem 24**.**

Let $T>0$. Suppose that

1. i.

$\Sigma(A,B,C,D)$* is passive with the storage function $x\mapsto\frac{1}{2}x^{\top}x$,* 2. ii.

$\mathcal{M}$* is maximal monotone,* 3. iii.

for all $t\in[0,T]$, we have111Here, $\operatorname{ri}(S)$ denotes the relative interior of $S$.* \operatorname{im}C\cap\operatorname{ri}\big{(}\operatorname{im}(\mathcal{M}+D)-v(t)\big{)}\neq\varnothing,* 4. iv.

$v$* is bounded on $[0,T]$,* 5. v.

there exists an absolutely continuous nondecreasing function $\theta:[0,T]\rightarrow\mathbb{R}$ such that

[TABLE]

for all $s,t$ with $0\leqslant s\leqslant t\leqslant T$. 6. vi.

There exists a positive number $\alpha$ such that

[TABLE]

for all $\eta\in\operatorname{dom}{(\mathcal{M}+D)^{-1}}$.

Then, $H$ satisfies assumptions (A1), (A2), (A3), and (LG).

**Proof. **For brevity, we define $W(t):=\operatorname{im}(\mathcal{M}+D)-v(t)$ for all $t\in[0,T]$. Note that $\operatorname{dom}H(t,\cdot)=C^{-1}W(t)$ for all $t\in[0,T]$.

(A1): It follows from (CamlSchu16, , Thm. 2) that the conditions (i.)-(iii.) imply that $H(t,\cdot)$ is a maximal monotone mapping for all $t\in[0,T]$. As such, $H$ satisfies assumption (A1).

(A3): Let $t$ and $s$ be such that $0\leqslant s\leqslant t\leqslant T$. Also, let $x\in C^{-1}W(s)$ and let y=\operatorname{proj}\big{(}x,C^{-1}W(t)\big{)}. Further, let \zeta=\operatorname{proj}\big{(}Cx,\operatorname{im}C\cap W(t)\big{)}. Therefore, there exists $\xi$ such that $\zeta=C\xi$. Without loss of generality, we can assume that $x-\xi\in\operatorname{im}C^{\top}$ since $\mathbb{R}^{n}=\operatorname{im}C^{\top}\oplus\ker C$. Now, we see that $Cx\in\operatorname{im}C\cap W(s)$ and \zeta=C\xi\in\operatorname{cl}\big{(}\operatorname{im}C\cap W(t)\big{)}. From (v.), we get

[TABLE]

Since $x-\xi\in\operatorname{im}C^{\top}$, there exists a positive number $\alpha$ such that

[TABLE]

Since $\xi\in C^{-1}W(t)$, we obtain

[TABLE]

Therefore, we see that

[TABLE]

This implies that

[TABLE]

Since $\operatorname{dom}H(t,\cdot)=C^{-1}W(t)$, we can conclude that $H$ satisfies assumption (A3).

(LG): From (vi.), we know that there exists a positive number $\alpha$ such that

[TABLE]

for all $\eta\in\operatorname{dom}{(\mathcal{M}+D)^{-1}}$. Let $x\in\operatorname{dom}{H(t,\cdot)}$ for some $t\in[0,T]$. Since $Ax-Bz^{0}\in H(t,x)$ where z^{0}\in\big{(}(\mathcal{M}+D)^{-1})^{0}\big{(}Cx+v(t)\big{)}, we have

[TABLE]

Since $x\in\operatorname{dom}{H(t,\cdot)}$, it follows that $Cx+v(t)\in\mathbb{B}^{m}(\rho)\cap\operatorname{dom}(\mathcal{M}+D)^{-1}$, and using (90), (91), and boundedness of $v$, we get

[TABLE]

for some positive number $\beta$ that does not depend on $t$. Therefore, we have

[TABLE]

for all $x\in\operatorname{dom}{H(t,\cdot)}$ with $t\in[0,T]$. In other words, $H$ satisfies assumption (LG).

(A2): Let $(t_{\ell},x_{\ell},y_{\ell})_{\ell\in\mathbb{N}}\subseteq[0,T]\times\mathbb{R}^{n}\times\mathbb{R}^{n}$ be a sequence such that $y_{\ell}\in H(t_{\ell},x_{\ell})$ and $\lim_{\ell\uparrow\infty}(t_{\ell},x_{\ell},y_{\ell})=(t,x,y)$ for some $t\in[0,T]$, $x,y\in\mathbb{R}^{n}$. What needs to be proven is that $y\in H(t,x)$.

Note that for each $\ell$ there exists z_{\ell}\in(\mathcal{M}+D)^{-1}\big{(}Cx_{\ell}+v(t_{\ell})\big{)} such that $y_{\ell}=-Ax_{\ell}+Bz_{\ell}$. Then, $(Bz_{\ell})_{\ell\in\mathbb{N}}$ converges. Let $\mathcal{W}$ be the subspace parallel to the affine hull of $\operatorname{im}(\mathcal{M}+D)=\operatorname{dom}(\mathcal{M}+D)^{-1}$. It follows from maximal monotonicity of $(\mathcal{M}+D)^{-1}$ that for each $\ell$

[TABLE]

holds for any $\zeta\in\mathcal{W}^{\perp}$. Now, let $z_{\ell}=z_{\ell}^{1}+z_{\ell}^{2}$ where $z_{\ell}^{1}\in\ker B\cap\mathcal{W}^{\perp}$ and

[TABLE]

Note that

[TABLE]

From (92) by taking $\zeta=-z_{\ell}^{1}$, we have

[TABLE]

Suppose that $(z^{2}_{\ell})_{\ell\in\mathbb{N}}$ is bounded. Hence, $(z^{2}_{\ell})_{\ell\in\mathbb{N}}$ converges on a subsequence $N$, say to $z$. Then, we see from $(x_{\ell},y_{\ell})=(x_{\ell},-Ax_{\ell}+Bz_{\ell})$ that $(x_{\ell},y_{\ell})_{\ell\in N}$ converges to $(x,y)=(x,-Ax+Bz)$. Since $H(t,\cdot)$ is maximal monotone and that implies the closedness of \operatorname{graph}\big{(}H(t,\cdot)\big{)}, we then can conclude that (x,y)\in\operatorname{graph}\big{(}H(t,\cdot)\big{)}, or equivalently $y\in H(t,x)$.

Therefore, it is enough to show that $(z^{2}_{\ell})_{\ell\in\mathbb{N}}$ is bounded. Suppose, on the contrary, that $z^{2}_{\ell}$ is unbounded. Without loss of generality, we can assume that the sequence $\frac{z^{2}_{\ell}}{|z^{2}_{\ell}|}$ converges. Define

[TABLE]

From (94) and the fact that $(Bz_{\ell})_{\ell\in\mathbb{N}}$ converges, we have

[TABLE]

Thus, we get

[TABLE]

Let $(\bar{x},\bar{y})\in\operatorname{graph}H(t,\cdot)$. Then, $\bar{y}=-A\bar{x}+B\bar{z}$ where

[TABLE]

Due to passivity with $K=I$ and monotonicity of $(\mathcal{M}+D)^{-1}$, it follows from (95) and (99) that

[TABLE]

By dividing by $|z^{2}_{\ell}|^{2}$, taking the limit as $\ell$ tends to infinity and using boundedness of $v$, we obtain

[TABLE]

Since $D$ is positive semi-definite due to the first statement of Proposition 1, this results in

[TABLE]

Then, it follows from (98), $K=I$, and the second statement of Proposition 1 that

[TABLE]

Let $\eta\in\operatorname{im}(\mathcal{M}+D)-v(t)$ and \zeta\in(\mathcal{M}+D)^{-1}\big{(}\eta+v(t)\big{)}. In view of monotonicity of $(\mathcal{M}+D)^{-1}$, we get

[TABLE]

from (95). Taking the limit as $\ell$ tends to infinity, employing boundedness of $v$, and using (102), we obtain

[TABLE]

From (102) and (104), we see that the hyperplane $\operatorname{span}(\{\zeta_{\infty}\})^{\perp}$ separates the sets $\operatorname{im}C$ and $\operatorname{im}(\mathcal{M}+D)-v(t)$. In view of $\operatorname{im}C=\operatorname{ri}(\operatorname{im}C)$ and (iii.), it follows from (rockafellar:70, , Thm. 11.3) that $\operatorname{im}C$ and $\operatorname{im}(\mathcal{M}+D)-v(t)$ cannot be properly separated. Therefore, both $\operatorname{im}C$ and $\operatorname{im}(\mathcal{M}+D)-v(t)$ must be contained in the hyperplane $\operatorname{span}(\{\zeta_{\infty}\})^{\perp}$. Thus, we see that $\operatorname{im}(\mathcal{M}+D)$ is contained in $v(t)+\operatorname{span}(\{\zeta_{\infty}\})^{\perp}$. Since $\mathcal{W}$ is the subspace parallel to the affine hull of $\operatorname{im}(\mathcal{M}+D)$, we get $\mathcal{W}\subseteq\operatorname{span}(\{\zeta_{\infty}\})^{\perp}$ which implies $\zeta_{\infty}\in\mathcal{W}^{\perp}$. Together with (98), we get

[TABLE]

In view of (93) and (96), this yields $\zeta_{\infty}=0$. This, however, clearly contradicts with (96). Therefore, $|z^{2}_{\ell}|$ must be bounded. $\blacksquare$

Next, we specialize the results of Theorem 24 to linear complementarity systems.

7.1 Linear complementarity systems

Linear complementarity systems are important instances of the differential inclusions of the form (81) with $\mathcal{M}$ described by so-called complementarity relations. In this section, we aim at presenting tailor-made conditions for existence and uniqueness of solutions to linear complementarity systems.

Consider a linear complementary system

[TABLE]

where $\mathcal{P}:\mathbb{R}^{m}\rightrightarrows\mathbb{R}^{m}$ is the maximal monotone set-valued mapping given by

[TABLE]

Next, we introduce the linear complementarity problem.

Given a vector $q\in\mathbb{R}^{m}$ and a matrix $M\in\mathbb{R}^{m\times m}$, the linear complementarity problem LCP($q,M$) is to find a vector $z\in\mathbb{R}^{m}$ such that

[TABLE]

We say that the LCP($q,M$) is feasible if there exists $z$ satisfying (106a) and (106b). If a vector $z$ is feasible and satisfies (106c) in addition, then we say that $z$ solves (is a solution of) LCP($q,M$). The set of all solutions of LCP($q,M$) will be denoted by SOL($q,M$).

A comprehensive study on LCPs can be found in cottle:92 . In the sequel, we will be interested in LCP($q,M$) where $M$ is a (not necessarily symmetric) positive semi-definite matrix.

Given a square matrix $M$, we define

[TABLE]

and its dual cone

[TABLE]

When $M$ is (not necessarily symmetric) a positive semi-definite matrix, the set $\mathcal{Q}_{M}$ is a convex cone and can be given by $\mathcal{Q}_{M}=\{z:z\geqslant 0,\,Mz\geqslant 0,\text{ and }(M+M^{\top})z=0\}$.

The following proposition characterizes the conditions under which an LCP with positive semi-definite $M$ matrix has solutions.

Proposition 25** (Cor. 3.8.10 of cottle:92 and Lem. 23 of camlibel:02a ).**

Let $M$ be a positive semi-definite matrix. Then, the following statements are equivalent:

1. i.

$q\in\mathcal{Q}_{M}^{+}$. 2. ii.

$\mathrm{LCP}(q,M)\text{ is feasible}$. 3. iii.

$\mathrm{LCP}(q,M)\text{ is solvable}$.

Moreover, the following statements hold:

1. iv.

For each $q\in\mathcal{Q}_{M}^{+}$, there exists a unique least-norm solution $z^{*}(q)\in\mathrm{SOL}(q,M)$ in the sense that $|z^{*}(q)|\leqslant|z|$ for all $z\in\mathrm{SOL}(q,M)$. 2. v.

There exists a positive number $\alpha$ such that

[TABLE]

Now, define

[TABLE]

Note that \operatorname{dom}H_{\mathcal{P}}(t,\cdot)=C^{-1}\big{(}\operatorname{im}(\mathcal{P}+D)-v(t)\big{)}. Moreover, $q\in(\mathcal{P}+D)(z)$ if and only if $-z\in\mathrm{SOL}(q,D)$. This means that $q\in(\mathcal{P}+D)(z)$ if and only if $q\in\mathcal{Q}_{D}^{+}$ in view of Proposition 25. In other words, \operatorname{dom}H_{\mathcal{P}}(t,\cdot)=C^{-1}\big{(}\mathcal{Q}_{D}^{+}-v(t)\big{)}.

The following theorem provides streamlined conditions that guarantee the hypotheses of Theorem 24 for the time-dependent set-valued mapping $H_{\mathcal{P}}$ as defined in (109).

Theorem 26**.**

Let $T>0$. Suppose that

1. i.

$\Sigma(A,B,C,D)$* is passive with the storage function $x\mapsto\frac{1}{2}x^{\top}x$,* 2. ii.

\operatorname{im}C\cap\operatorname{ri}\big{(}\operatorname{im}(\mathcal{P}+D)-v(t)\big{)}\neq\varnothing* for all $t\in[0,T]$,* 3. iii.

$v\in AC([0,T],\mathbb{R}^{m})$,

Then, $H_{\mathcal{P}}$ satisfies assumptions (A1), (A2), (A3), and (LG).

**Proof. **It is enough to show that $H_{\mathcal{P}}$ satisfies the hypotheses (i.)-(vi.) of Theorem 24. The first four hypotheses of Theorem 24 are readily satisfied. Therefore, we need to show that the remaining two also hold.

For the hypothesis (v.) of Theorem 24, we need a streamlined version of Hoffman’s bound on the polyhedral sets. To elaborate, let $\varnothing\neq\mathcal{R}\subseteq\mathbb{R}^{m}$ be a polyhedral set given by $\mathcal{R}=\{\zeta:R\zeta=0\text{ and }Q\zeta\leqslant q\}$ where $R,Q$ are matrices and $q$ is vector with appropriate sizes. Hoffman’s bound (see e.g. (pang, , Lemma 3.2.3)) asserts that there exists a positive number $\alpha$ that depend on $\mathcal{R}$ such that

[TABLE]

for all $x\in\mathbb{R}^{m}$ where $\max$ denotes componentwise maximum. By definition $\mathcal{Q}_{D}^{+}$ is a polyhedral cone. Therefore, we have $\mathcal{Q}_{D}^{+}=\{\eta\in\mathbb{R}^{m}:Q\eta\leqslant 0\}$ for some matrix $Q$. Let $E$ be a matrix such that $\operatorname{im}C=\ker E$. Then, we have

[TABLE]

for all $t\in[0,T]$.

Let $s,t$ be such that $0\leqslant s\leqslant t\leqslant T$ and w\in\operatorname{im}C\cap\big{(}\operatorname{im}(\mathcal{P}+D)-v(s)\big{)}. From (111), we see that

[TABLE]

Now, we have

[TABLE]

where $\beta$ is a positive number. Since $v$ is absolutely continuous, we have

[TABLE]

Then, (113) implies that

[TABLE]

for all $s,t$ with $0\leqslant s\leqslant t\leqslant T$ where $\theta(t)=\frac{1}{\beta}\int_{0}^{t}|\dot{v}(\tau)|\,d\tau$ for all $t\in[0,T]$. Clearly, $\theta$ is nondecreasing and absolutely continuous.

For the hypothesis (vi.) of Theorem 24, note that $\zeta\in(\mathcal{P}+D)^{-1}(\eta)$ if and only if $-\zeta\in\mathrm{SOL}(\eta,D)$. Due to Proposition 1, $D$ is positive semi-definite. Then, it follows from Proposition 25 that there exists a positive number $\alpha$ such that

[TABLE]

for all $\eta\in\operatorname{dom}{(\mathcal{P}+D)^{-1}}$. Therefore, $H_{\mathcal{P}}$ satisfies the hypothesis (vi.) of Theorem 24. $\blacksquare$

8 Conclusions

In this article, we have studied the existence of solutions to differential inclusions with time-dependent maximal monotone operators. With the help of an example, it is shown that our proposed conditions overcome the limitations of existing results. As a particular class of these inclusions, we consider differential equations coupled with time-dependent complementarity relations. For this system class, conditions for existence of solutions are derived explicitly in terms of system data. To build on these results, the conditions for existence of solutions can be relaxed for differential inclusions where the maximal monotone operators have a particular structure, for example TanwBrog18 .

Moving forward from the question of existence of solutions, it is also of interest to study the qualitative properties of the solutions of such systems, such as continuity with respect to initial data pang2009 . One can also investigate stability related problems for the generic class of dynamical systems, as has been done for some specific set-valued systems in TanwBrog14 . It also remains to be seen whether our proposed results provide any advantages in the study of optimal control problems, such as Briceno-Arias15 .