## Anti-selfdual Lagrangians: Variational resolutions of non self-adjoint equations and dissipative evolutions

### January 20, 2005

Abstract We develop the concept and the calculus of anti-self dual (ASD) Lagrangians which seems inherent to many questions in mathematical physics, geometry, and differential equations. They are natural extensions of gradients of convex functions –hence of self-adjoint positive operators– which usually drive dissipative systems, but also rich enough to provide representations for the superposition of such gradients with skew-symmetric operators which normally generate unitary flows. They yield variational formulations and resolutions for large classes of non-potential boundary value problems and initial-value parabolic equations. Solutions are minima of functionals of the form $I\left(u\right)=L\left(u,\Lambda u\right)$  (resp. $I\left(u\right)={\int }_{0}^{T}L\left(t,u\left(t\right),\stackrel{˙}{u}\left(t\right)+{\Lambda }_{t}u\left(t\right)\right)dt$  ) where $L$  is an anti-self dual Lagrangian and where ${\Lambda }_{t}$  are essentially skew-adjoint operators. However, and just like the self (and antiself ) dual equations of quantum field theory (e.g. Yang-Mills) the equations associated to such minima are not derived from the fact they are critical points of the functional $I$  , but because they are also zeroes of the Lagrangian $L$  itself. The approach has many advantages: It solves variationally many equations and systems that cannot be obtained as Euler-Lagrange equations of action functionals, since they can involve non self-adjoint or other non-potential operators; It also associates variational principles to variational inequalities, and to various dissipative initial-value first order parabolic problems.
These equations can therefore now be analyzed with the full range of methods –computational or not– that are available for variational settings. Most remarkable are the permanence properties that ASD Lagrangians possess making them more pervasive than expected and quite easy to construct.
In this first of a series of papers, we describe the basic theory of anti-self dual Lagrangians and some of its first applications involving mostly bounded linear operators. In forthcoming papers, we extend the theory to deal with non bounded operators ([20) and with non-linear maps ([16) .
Contents

1 Introduction

Non self-adjoint problems such as the transport equation:
 $\begin{array}{c}\left\{\begin{array}{ccc}-{\Sigma }_{i=1}^{n}{a}_{i}\frac{\partial u}{\partial {x}_{i}}+{a}_{0}u& =& \beta \left(u\right)+fon\Omega \subset I{R}^{n}\\ u\left(x\right)& =& 0on{\Sigma }_{-}.\end{array}\end{array}$ (1)
where $\mathbf{a}=\left({a}_{i}{\right)}_{i}:\Omega \to {\mathbf{R}}^{n}$  is a smooth vector field, $\beta$  is a convex function, $f\in {L}^{2}\left(\Omega \right)$  , and where ${\Sigma }_{-}=\left\{x\in \partial \Omega ;\mathbf{a}\left(x\right)\cdot \mathbf{n}\left(x\right)<0\right\}$  , $\mathbf{n}$  being the outer normal vector, are not of Euler-Lagrange type and their solutions are not normally obtained as critical points of functionals of the form ${\int }_{\Omega }F\left(x,u\left(x\right),\nabla u\left(x\right)\right)dx$  . Similarly, dissipative initial value problems such as the heat equation or those describing porous media:
 $\begin{array}{c}\left\{\begin{array}{ccc}\frac{\partial u}{\partial t}& =& \Delta {u}^{m}+fon\Omega ×\left[0,T\right]\\ u\left(0,x\right)& =& {u}_{0}\left(x\right)on\Omega ,\\ \end{array}\end{array}$ (2)
are not normally solved by the methods of the calculus of variations since they do not correspond to Euler-Lagrange equations of action functionals of the form ${\int }_{0}^{T}L\left(t,x\left(t\right),\stackrel{˙}{x}\left(t\right)dt$  .
However, physicists have managed to obtain variationally many of the basic first order equations of quantum field theory by minimizing their associated action functionals. These are the celebrated self (antiself ) dual equations of Yang-Mills, Seiberg-Witten and Ginzburg-Landau which are not derived from the fact they are critical points (i.e., from the corresponding Euler-Lagrange equations) but from the fact that they are zeros of the Lagrangian itself, which is the case as long as the action functional attains a natural and –a priori– known minimum (See for example [22).
From a totally different perspective, Brezis and Ekeland formulated about 30 years ago in [8an intriguing minimization principle which can be associated to the heat equation and other gradient flows of convex energy functionals. Again the applicability of their principle was conditional on identifying the minimum value of the functional. Later, Auchmuty ([1, [2) proposed a framework in which he formalizes and generalizes the Brezis-Ekeland procedure in order to apply it to operator equations of non-potential type. However, the applicability of this variational principle remained conditional on evaluating the minimum value and in most cases could not be used to establish existence and uniqueness of solutions.
In this paper, we develop a general framework where such variational principles are applicable. It is based on the concept of anti-selfdual (ASD) Lagrangians which seems inherent to many important boundary value problems as well as parabolic evolution equations. For such Lagrangians $L$  and for skew-adjoint operators ${\Lambda }_{t}$  , solutions are obtained as minima of functionals of the form $I\left(u\right)=L\left(u,\Lambda u\right)orI\left(u\right)={\int }_{0}^{T}L\left(t,u\left(t\right),\stackrel{˙}{u}\left(t\right)+{\Lambda }_{t}u\left(t\right)\right)dt.$  The minimal value will always be zero, and the equations associated to such minima are not derived from the fact they are critical points of the functional $I$  , but because they are also zeroes of the Lagrangian $L$  itself. More specifically, the solutions will satisfy $L\left(u,\Lambda u\right)+〈u,\Lambda u〉=0andL\left(t,u\left(t\right),\stackrel{˙}{u}\left(t\right)+{\Lambda }_{t}u\left(t\right)\right)+〈u\left(t\right),\stackrel{˙}{u}\left(t\right)〉=0,$  for all time, which is reminiscent of the conservation laws enjoyed by Hamiltonians. This provides variational formulations and complete proofs for the conditional results of Brezis-Ekeland, Auchmuty and others.
As importantly, we show that ASD Lagrangians possess remarkable permanence properties making them more prevalent than expected and quite easy to construct and/or identify. The variational game changes from the analytical proofs of existence of extremals for general Lagrangians, to a more algebraic search of an appropriate ASD Lagrangian for which the minimization problem is remarkably simple. This makes them efficient new tools for proving existence and uniqueness results for a large array of differential equations.The basic idea is simple and is an elaboration on our work in [19where we gave complete variational proofs of the existence and uniqueness of gradient flows of convex energy functionals, and the one in [15, where we give a variational proof for the existence and uniqueness of solutions of certain non-linear transport equations. Starting with an equation of the form
 $\begin{array}{c}-Au\in \partial \phi \left(u\right)\end{array}$ (3)
it is well known that it can be formulated –and sometimes solved– variationally whenever $A:X\to {X}^{*}$  is a selfadjoint bounded linear operator and $\phi$  is an appropriate functional on $X$  . Indeed, in this case it can be reduced to the equation $0\in \partial \psi \left(u\right)$  , where $\psi$  is the functional
 $\begin{array}{c}\psi \left(u\right)=\phi \left(u\right)+\frac{1}{2}〈Au,u〉.\end{array}$ (4)
A solution can then be obtained by minimization whenever $\phi$  is convex and lower semi-continuous and whenever $A$  is positive (i.e., $〈u,Au〉\ge 0$  ) or better if $A$  is coercive (i.e., if for some $c>0$  , $〈u,Au〉\ge c\parallel u{\parallel }^{2}$  for all $u\in X$  ).
But this variational procedure fails when $A$  is not self-adjoint, or when $A$  is a non-potential operator (i.e., when $A$  is not a gradient vector field), and definitely when $A$  is not linear. In this case, the Brezis-Ekeland procedure –as formalized by Auchmuty– consists of simply minimizing the functional
 $\begin{array}{c}I\left(u\right)=\phi \left(u\right)+{\phi }^{*}\left(-Au\right)+〈u,Au〉\end{array}$ (5)
where ${\phi }^{*}$  is the Legendre dual of $\phi$  defined on ${X}^{*}$  by ${\phi }^{*}\left(p\right)=sup\left\{〈x,p〉-\phi \left(x\right);x\in X\right\}.$  Legendre duality yields that $\alpha :={inf}_{u\in X}I\left(u\right)\ge 0$  , and the key observation made by several authors is the following simple Fact: If the infimum $\alpha =0$  and if it is attained at $\overline{u}\in X$  then we are in the limiting case of the Fenchel-Legendre duality, $\phi \left(\overline{u}\right)+{\phi }^{*}\left(-A\overline{u}\right)=〈\overline{u},-A\overline{u}〉$  and therefore $-A\overline{u}\in \partial \phi \left(\overline{u}\right)$  .
Note that the procedure does not require any assumption on $A$  , and very general coercivity assumptions on $\phi$  often ensure the existence of a minimum. However, the difficulty here is different from standard minimization problems in that besides the problem of existence of a minimum, one has to insure that the infimum is actually zero. This is obviously not the case for general operators $A$  , though one can always write (and many authors did) the variational principle (5) for the operator equation (3).
In this paper, we tackle the real difficulty of when the infimum $\alpha$  is actually zero and we try to identify a class of nonpotential operators $F\left(u\right)$  for which the equation and the initial-value problem
 $\begin{array}{c}0\in F\left(u\right)and\left\{\begin{array}{ccc}-\stackrel{˙}{u}\left(t\right)& \in & F\left(u\left(t\right)\right)\\ u\left(0\right)& =& {u}_{0}\\ \end{array}\end{array}$ (6)
can be solved by the above variational procedure. We show here that this is essentially the case whenever $F\left(u\right)=Bu+\partial \phi \left(u\right)$  , where $\phi$  is a convex lower semicontinuous function and when $B$  is a skew-adjoint operator. We note that –when $\Lambda$  is linear– such operators form a very important subset of the class of maximal monotone operators for which there is already an extensive theory ([7, [4). The interest here is in the new variational approach based on the concept of anti-selfdual Lagrangians which possesses remarkable permanence properties that maximal monotone operators either do not satisfy or do so via substantially more elaborate methods. In a forthcoming paper ([16) we establish similar results for operators of the form $F\left(u\right)=\Lambda u+Bu+\partial \phi \left(u\right)$  where $\Lambda$  is an appropriate non-linear conservative operator, $B$  is linear and positive, and $\phi$  is convex, the superposition of which is not normally covered by the theory of maximal monotone operators.
In this paper, we establish the algebraic structure of ASD Lagrangians, emphasizing issues on how to build and identify complex ASD Lagrangians from the more basic ones. To keep the key ideas transparent, we chose to deal with the case when the operators are bounded and linear, leaving the more analytically involved cases of unbounded and nonlinear operators to forthcoming papers.
This –bounded linear– case already has many interesting features, especially in boundary value problems of the form:
 $\begin{array}{c}\left\{\begin{array}{ccc}-Ax+f& \in & \partial \phi \left(x\right)\\ B\left(x\right)& =& a\\ \end{array}\end{array}$ (7)
where $B$  is a boundary operator on $X$  (related to the positive operator $A$  ), as well as parabolic evolution equations of the form:
 $\begin{array}{c}\left\{\begin{array}{ccc}-{A}_{t}x\left(t\right)-\stackrel{˙}{x}\left(t\right)& \in & \partial \phi \left(t,x\left(t\right)\right)a.e.t\in \left[0,T\right]\\ {B}_{t}\left(x\left(t\right)\right)& =& a\left(t\right)a.et\in \left[0,T\right]\\ x\left(0\right)& =& {x}_{0}\end{array}\end{array}$ (8)
where ${x}_{0}$  is a given initial value and where $a\left(t\right)$  is a prescribed boundary value.
We start by presenting – in section 2– the special variational properties of the class of $R$  -Antiselfdual Lagrangians, where $R$  is any automorphism of the state space. This should already give an idea of their relevance in the existence theory of certain PDEs, and will hopefully motivate the in-depth study of their permanence properties. Beyond this first section, we will only deal with the anti-symmetric case, i.e., when $R\left(x\right)=-x$  , in which case $R$  -antiselfdual Lagrangians will be called anti-selfdual Lagrangians (ASD). We shall see that this class of Lagrangians already covers a great deal of applications which warranted that this paper as well as ([20, [16) be solely devoted to this case. However, the theory involving other automorphisms $R$  will also be very useful, especially in applications to Hamiltonian systems and this will be developed in [21.
In section 3, we establish the basic permanence properties of anti-self dual Lagrangians as well as their special variational features while focussing on stationary equations and systems. This restrictive looking class turns out to be quite rich. In section (3) we deal with boundary value problems where appropriate selfdual boundary Lagrangians are appropriately added to the “interior Lagrangian” to make it anti-selfdual allowing us to solve problems with prescribed boundary terms.
In section (4), we show how ASD Lagrangians “lift” to path spaces allowing us to solve with the same variational approach several parabolic equations –including gradient flows. In section (5), we associate to each autonomous Lagrangian, a semi-group of contractions which emphasizes again that such Lagrangians are natural extensions of gradients of convex functions, of positive operators as well as of the ”superpositions” of the two actions. In section 6, we give a glimpse on how the theory can help in solving variationally certain implicit PDEs, a project for future investigation.
As mentioned above, in this paper we describe the basics of the ASD theory emphasizing its stability under various operations and its rich structure. So we stuck with the simplest of examples leaving more complicated PDE settings to forthcoming papers. In ([20), we extend the theory to deal with linear but unbounded operators, and in ([16) we tackle various non-linear but appropriately defined “skew-adjoint” operators such as those appearing in the Navier-Stokes and other equations of hydrodynamics. Finally, I would like to thank Yann Brenier, Eric Séré, Leo Tzou and Abbas Moameni for the many extremely fruitful discussions and their valuable input into this project.

2 Basic variational properties of R-antiselfdual Lagrangians

We consider the class $\mathcal{ℒ}\left(X\right)$  of convex Lagrangians $L$  on a reflexive Banach space $X$  : these are all functions $L:X×{X}^{*}\to IR\cup \left\{+\infty \right\}$  which are convex and lower semi-continuous (in both variables) and which are not identically $+\infty$  . The Legendre-Fenchel dual (in both variables) of $L$  is defined at any pair $\left(q,y\right)\in {X}^{*}×X$  by:
${L}^{*}\left(q,y\right)=sup\left\{〈q,x〉+〈y,p〉-L\left(x,p\right);x\in X,p\in {X}^{*}\right\}$
Definition 2.1 Given a bounded linear operator $R:X\to X$  , say that:
(1) $L$  is an $R$  -antiselfdual Lagrangian on $X×{X}^{*}$  , if
 $\begin{array}{c}{L}^{*}\left(p,x\right)=L\left(-Rx,-{R}^{*}p\right)forall\left(p,x\right)\in {X}^{*}×X.\end{array}$ (9)
(2) $L$  is partially $R$  -antiselfdual, if
 $\begin{array}{c}{L}^{*}\left(0,x\right)=L\left(-Rx,0\right)forallx\in X.\end{array}$ (10)
(3) $L$  is $R$  -antiselfdual on the graph of $\Lambda$  , the latter being a map from $X$  into ${X}^{*}$  , if
 $\begin{array}{c}{L}^{*}\left(\Lambda x,x\right)=L\left(-Rx,-{R}^{*}\circ \Lambda x\right)forallx\in X.\end{array}$ (11)
(4) More generally, if $Y×Z$  is any subset of $X×{X}^{*}$  , we shall say that $L$  is $R$  -antiself dual on the elements of $Y×Z$  if ${L}^{*}\left(p,x\right)=L\left(-Rx,-{R}^{*}p\right)$  for all $\left(p,x\right)\in Y×Z$  .
A typical example of an $R$  -antiselfdual Lagrangian is $L\left(x,p\right)=\phi \left({R}^{-1}x\right)+{\phi }^{*}\left(-p\right)$  and $M\left(x,p\right)=\phi \left(-x\right)+{\phi }^{*}\left(\left({R}^{*}{\right)}^{-1}p\right)$  where $\phi$  is a convex lower semi-continuous function and $R$  is an invertible operator on $X$  . More generally, $L\left(x,p\right)=\phi \left(Rx\right)+{\phi }^{*}\left(-{S}^{*}p\right)$  is an $\left(S\circ R{\right)}^{-1}$  -antiselfdual Lagrangian. Moreover, if $\Lambda :X\to {X}^{*}$  is such that $\Lambda \circ \left(S\circ R{\right)}^{-1}$  is skew-adjoint, then $L\left(x,p\right)=\phi \left(Rx\right)+{\phi }^{*}\left(-{S}^{*}\Lambda x-{S}^{*}p\right)$  is also an $\left(S\circ R{\right)}^{-1}$  -antiselfdual Lagrangian. Our basic premise in this paper is that many boundary value problems can be solved by minimizing functionals of the form $I\left(x\right)=L\left(x,\Lambda x\right)$  where $L$  is a $R$  -Antiselfdual Lagrangian and provided $\Lambda \circ R$  is a skew-adjoint operator. However, their main relevance to our study stems from the fact that –generically– the infimum is actually equal to $0$  . It is this latter property that allows for novel variational formulations and resolutions of several basic PDEs and evolution equations, which –often because of lack of self-adjointness– do not normally fit the Euler-Lagrange framework.
As mentioned above, if $L$  is a $R$  -Antiselfdual Lagrangian and if $\Lambda :X\to {X}^{*}$  is an operator such that $\Lambda \circ R$  is skew adjoint, then the Lagrangian ${L}_{\Lambda }\left(x,p\right)=L\left(x,\Lambda x+p\right)$  is again $R$  -Antiselfdual. In other words, Minimizing $L\left(x,\Lambda x\right)$  amounts to minimizing ${L}_{\Lambda }\left(x,0\right)$  which is covered by the following very simple –yet far reaching– proposition. Again, its relevance comes from the evaluation of the minimum and not from the –more standard– question about its attainability.
We start by noticing that for a $R$  -Antiselfdual Lagrangian, we readily have:
 $\begin{array}{c}L\left(Rx,{R}^{*}p\right)\ge -〈Rx,p〉forevery\left(x,p\right)\in X×{X}^{*},\end{array}$ (12)
and if $L$  is partially anti-selfdual, then
 $\begin{array}{c}I\left(x\right)=L\left(Rx,0\right)\ge 0foreveryx\in X,\end{array}$ (13)
So, we are looking into an interesting variational situation, where the minima can also be zeros of the functionals. Here are some necessary conditions for the existence of such minima.
Proposition 2.1 Let $L$  be a convex lower-semi continuous functional on a reflexive Banach space $X×{X}^{*}$  . Assume that $L$  is a partially $R$  -Antiselfdual Lagrangian and that for some ${x}_{0}\in X$  , the function $p\to L\left({x}_{0},p\right)$  is bounded above on a neighborhood of the origin in ${X}^{*}$  . Then there exists $\overline{x}\in X$  , such that:
 $\begin{array}{c}\left\{\begin{array}{ccc}L\left(-R\overline{x},0\right)& =& {inf}_{x\in X}L\left(x,0\right)=0.\\ \left(0,\overline{x}\right)& \in & \partial L\left(-R\overline{x},0\right).\end{array}\end{array}$ (14)
Proof: This follows from the basic duality theory in convex optimization. Indeed, if $\left({\mathcal{P}}_{p}\right)$  is the primal minimization problem $h\left(p\right)={inf}_{x\in X}L\left(x,p\right)$  in such a way that $\left({\mathcal{P}}_{0}\right)$  is the initial problem $h\left(0\right)={inf}_{x\in X}L\left(x,0\right)$  , then the dual problem $\left({\mathcal{P}}^{*}\right)$  is ${sup}_{y\in X}-{L}^{*}\left(0,y\right)$  , and we have the weak duality formula $inf{\mathcal{P}}_{0}:={inf}_{x\in X}L\left(x,0\right)\ge {sup}_{y\in X}-{L}^{*}\left(0,y\right):=sup{\mathcal{P}}^{*}.$  The “partial $R$  -Antiselfdual ity” of $L$  gives that
 $\begin{array}{c}{inf}_{x\in X}L\left(x,0\right)\ge {sup}_{y\in X}-{L}^{*}\left(0,y\right)={sup}_{y\in X}-L\left(-Ry,0\right).\end{array}$ (15)
Note now that $h$  is convex on ${X}^{*}$  and that its Legendre conjugate satisfies ${h}^{*}\left(y\right)={L}^{*}\left(0,y\right)=L\left(-Ry,0\right)$  on $X$  . If now $h$  is subdifferentiable at $0$  (i.e., if the problem $\left({\mathcal{P}}_{0}\right)$  is stable), then for any $\overline{x}\in \partial h\left(0\right)$  , we have $h\left(0\right)+{h}^{*}\left(\overline{x}\right)=0$  , which means that $-{inf}_{x\in X}L\left(x,0\right)=-h\left(0\right)={h}^{*}\left(\overline{x}\right)={L}^{*}\left(0,\overline{x}\right)=L\left(-R\overline{x},0\right)\ge {inf}_{x\in X}L\left(x,0\right).$  It follows that ${inf}_{x\in X}L\left(x,0\right)=L\left(-R\overline{x},0\right)\le 0$  and in view of ( 13 ), we get that the infimum of $\left(\mathcal{P}\right)$  is zero and attained at $-R\overline{x}$  , while the supremum of $\left({\mathcal{P}}^{*}\right)$  is attained at $\overline{x}$  . In this case we can write $L\left(-R\overline{x},0\right)+{L}^{*}\left(0,\overline{x}\right)=0$  which yields that $\left(0,\overline{x}\right)\in \partial L\left(-R\overline{x},0\right)$  .
If now for some ${x}_{0}\in X$  , the function $p\to L\left({x}_{0},p\right)$  is bounded above on a neighborhood of the origin in ${X}^{*}$  , then $h\left(p\right)\le {inf}_{x\in X}L\left(x,p\right)\le L\left({x}_{0},p\right)$  and therefore $h$  is subdifferentiable at $0$  and we are done.
Remark 2.2 The above holds under the condition that $x\to L\left(Rx,0\right)$  is coercive in the following sense:
 $\begin{array}{c}{lim}_{\parallel x\parallel \to \infty }\frac{L\left(Rx,0\right)}{\parallel x\parallel }=+\infty .\end{array}$ (16)
Indeed since ${h}^{*}\left(y\right)={L}^{*}\left(0,y\right)=L\left(Ry,0\right)$  on $X$  , we get that that ${h}^{*}$  is coercive on $X$  , which means that $h$  is bounded above on neighborhoods of zero in ${X}^{*}$  .
Remark 2.3 The proof above requires only that $L$  is a Lagrangian satisfying
 $\begin{array}{c}{L}^{*}\left(0,x\right)\ge L\left(-Rx,0\right)\ge 0forallx\in X.\end{array}$ (17)
Now we can deduce the following
Theorem 2.4 Let $R:X\to X$  be a bounded linear operator on a reflexive Banach space $X$  and let $\Lambda :X\to {X}^{*}$  be another operator such that $\Lambda \circ R$  is skew adjoint. Let $L$  be a Lagrangian on $X$  that is $R$  -antiselfdual on the graph of $-{\Lambda }^{*}$  , and assume that ${lim}_{\parallel x\parallel \to \infty }\frac{L\left(Rx,\Lambda Rx\right)}{\parallel x\parallel }=+\infty$  . Then there exists $\overline{x}\in X$  , such that:
 $\begin{array}{c}\left\{\begin{array}{ccc}L\left(-R\overline{x},-\Lambda R\overline{x}\right)& =& {inf}_{x\in X}L\left(x,\Lambda x\right)=0.\\ \left(-{\Lambda }^{*}\overline{x},\overline{x}\right)& \in & \partial L\left(-R\overline{x},-\Lambda R\overline{x}\right).\end{array}\end{array}$ (18)
Proof: We first prove that the Lagrangian defined as $M\left(x,p\right)=L\left(x,\Lambda x+p\right)$  is partially $R$  -Antiselfdual. Indeed fix $\left(q,y\right)\in {X}^{*}×X$  , set $r=\Lambda x+p$  and write:
 $\begin{array}{ccc}{M}^{*}\left(q,y\right)& =& sup\left\{〈q,x〉+〈y,p〉-L\left(x,\Lambda x+p\right);\left(x,p\right)\in X×{X}^{*}\right\}\end{array}$
 $\begin{array}{ccc}& =& sup\left\{〈q,x〉+〈y,r-\Lambda x〉-L\left(x,r\right);\left(x,r\right)\in X×{X}^{*}\right\}\end{array}$
 $\begin{array}{ccc}& =& sup\left\{〈q-{\Lambda }^{*}y,x〉+〈y,r〉-L\left(x,r\right);\left(x,r\right)\in X×{X}^{*}\right\}\end{array}$
 $\begin{array}{ccc}& =& {L}^{*}\left(q-{\Lambda }^{*}y,y\right).\end{array}$
If $q=0$  , then ${M}^{*}\left(0,y\right)={L}^{*}\left(-{\Lambda }^{*}y,y\right)=L\left(-Ry,{R}^{*}{\Lambda }^{*}y\right)=L\left(-Ry,-\Lambda Ry\right)=M\left(-Ry,0\right)$  , and $M$  is therefore partially $R$  -Antiselfdual.
It follows from the previous proposition applied to $M$  , that there exists $\overline{x}\in X$  such that:
$L\left(-R\overline{x},-\Lambda R\overline{x}\right)=M\left(-R\overline{x},0\right)={inf}_{x\in X}M\left(x,0\right)={inf}_{x\in X}L\left(x,\Lambda x\right)=0.$  Now note that $L\left(-R\overline{x},-\Lambda R\overline{x}\right)=L\left(-R\overline{x},{R}^{*}{\Lambda }^{*}\overline{x}\right)={L}^{*}\left(-{\Lambda }^{*}\overline{x},\overline{x}\right),$  hence $L\left(-R\overline{x},-\Lambda R\overline{x}\right)+{L}^{*}\left(-{\Lambda }^{*}\overline{x},\overline{x}\right)=0=〈\left(-R\overline{x},-\Lambda R\overline{x}\right),\left(-{\Lambda }^{*}\overline{x},\overline{x}\right)〉.$  It follows from the limiting case of Legendre duality that $\left(-{\Lambda }^{*}\overline{x},\overline{x}\right)\in \partial L\left(-R\overline{x},-\Lambda R\overline{x}\right)$  .

3 Permanence properties of Anti-selfdual Lagrangians

The concept of $R$  -Antiselfduality for a general automorphism $R$  is relevant for dealing with certain Hamiltonian systems [21and will be pursued in full generality in a forthcoming paper [17. We shall however concentrate in the sequel on the class of anti-selfdual Lagrangians (ASD), meaning those $R$  -Antiselfdual Lagrangians corresponding to the identity operator $R\left(x\right)=x$  . In other words, (1) $L$  is said to be an anti-selfdual Lagrangian on $X×{X}^{*}$  , if
 $\begin{array}{c}{L}^{*}\left(p,x\right)=L\left(-x,-p\right)forall\left(p,x\right)\in {X}^{*}×X.\end{array}$ (19)
(2) $L$  is partially anti-self dual, if
 $\begin{array}{c}{L}^{*}\left(0,x\right)=L\left(-x,0\right)forallx\in X.\end{array}$ (20)
(3) $L$  is anti-self dual on the graph of $\Lambda$  , the latter being a map from $X$  into ${X}^{*}$  , if
 $\begin{array}{c}{L}^{*}\left(\Lambda x,x\right)=L\left(-x,-\Lambda x\right)forallx\in X.\end{array}$ (21)
(4) More generally, if $Y×Z$  is any subset of $X×{X}^{*}$  , we shall say that $L$  is anti-self dual on the elements of $Y×Z$  if ${L}^{*}\left(p,x\right)=L\left(-x,-p\right)$  for all $\left(p,x\right)\in Y×Z$  .
Denote by ${\mathcal{ℒ}}_{AD}\left(X\right)$  the class of anti-selfdual (ASD) Lagrangians on a given Banach space $X$  . We shall see that this is already a very interesting and natural class of Lagrangians as they appear in several basic PDEs and evolution equations. The basic example of an anti-selfdual Lagrangian is given by a function $L$  on $X×{X}^{*}$  , of the form
 $\begin{array}{c}L\left(x,p\right)=\phi \left(x\right)+{\phi }^{*}\left(-p\right)\end{array}$ (22)
where $\phi$  is a convex and lower semi-continuous function on $X$  and ${\phi }^{*}$  is its Legendre conjugate on ${X}^{*}$  . We shall call them the Basic ASD-Lagrangians. A key element of this theory is that the family of ASD Lagrangians is much richer and goes well beyond convex functions and their conjugates, since they are naturally compatible with skew-symmetric operators. Indeed if $\Lambda :X\to {X}^{*}$  is skew-symmetric (i.e., ${\Lambda }^{*}=-\Lambda$  ), the Lagrangian
 $\begin{array}{c}M\left(x,p\right)=\phi \left(x\right)+{\phi }^{*}\left(-\Lambda x-p\right)\end{array}$ (23)
is also anti-self dual, and if in addition $\Lambda$  is invertible then the same holds true for
 $\begin{array}{c}N\left(x,p\right)=\phi \left(x+{\Lambda }^{-1}p\right)+{\phi }^{*}\left(\Lambda x\right).\end{array}$ (24)
Basic properties of ASD Lagrangians The class ${\mathcal{ℒ}}_{AD}\left(X\right)$  enjoys a remarkable number of permanence properties. Indeed, we define on the class of Lagrangians $\mathcal{ℒ}\left(X\right)$  the following operations:
• Scalar multiplication: If $\lambda >0$  and $L\in \mathcal{ℒ}\left(X\right)$  , define the Lagrangian $\lambda \cdot L$  on $X×{X}^{*}$  by:
$\left(\lambda \cdot L\right)\left(x,p\right)={\lambda }^{2}L\left(\frac{x}{\lambda },\frac{p}{\lambda }\right).$
• Addition: If $L,M\in \mathcal{ℒ}\left(X\right)$  , define the Lagrangian $L+M$  on $X×{X}^{*}$  by:
$\left(L\oplus M\right)\left(x,p\right)=inf\left\{L\left(x,r\right)+M\left(x,p-r\right);r\in {X}^{*}\right\}$
• Convolution: If $L,M\in \mathcal{ℒ}\left(X\right)$  , define the Lagrangian $L◆M$  on $X×{X}^{*}$  by:
$\left(L◆M\right)\left(x,p\right)=inf\left\{L\left(z,p\right)+M\left(x-z,p\right);z\in X\right\}$
• Right operator shift: If $L\in \mathcal{ℒ}\left(X\right)$  and $\Lambda :X\to {X}^{*}$  is a bounded linear operator, define the Lagrangian ${L}_{\Lambda }$  on $X×{X}^{*}$  by ${L}_{\Lambda }\left(x,p\right):=L\left(x,\Lambda x+p\right).$
• Left operator shift: If $L\in \mathcal{ℒ}\left(X\right)$  and if $\Lambda :X\to {X}^{*}$  is an invertible operator, define the Lagrangian ${}_{\Lambda }L$  on $X×{X}^{*}$  by:
${}_{\Lambda }L\left(x,p\right):=L\left(x+{\Lambda }^{-1}p,\Lambda x\right).$
• Free product: If $\left\{{L}_{i};i\in I\right\}$  is a finite family of Lagrangians on reflexive Banach spaces $\left\{{X}_{i};i\in I\right\}$  , define the Lagrangian $L:={\Sigma }_{i\in I}{L}_{i}$  on $\left({\Pi }_{i\in I}{X}_{i}\right)×\left({\Pi }_{i\in I}{X}_{i}^{*}\right)$  by $L\left(\left({x}_{i}{\right)}_{i},\left({p}_{i}{\right)}_{i}\right)={\Sigma }_{i\in I}{L}_{i}\left({x}_{i},{p}_{i}\right).$
• Twisted A-product: If $L\in \mathcal{ℒ}\left(X\right)$  and $M\in \mathcal{ℒ}\left(Y\right)$  where $X$  and $Y$  are two reflexive spaces, then for any bounded linear operator $A:X\to {Y}^{*}$  , define the Lagrangian $L{\oplus }_{A}M$  on $\left(X×Y\right)×\left({X}^{*}×{Y}^{*}\right)$  by $\left(L{\oplus }_{A}M\right)\left(\left(x,y\right),\left(p.q\right)\right):=L\left(x,{A}^{*}y+p\right)+M\left(y,-Ax+q\right).$
• A-antidualisation: If $\phi$  is any convex function on $X×Y$  and $A$  is any bounded linear operator $A:X\to {Y}^{*}$  , define the Lagrangian $L{\oplus }_{as}A$  on $\left(X×Y\right)×\left({X}^{*}×{Y}^{*}\right)$  by $\phi {\oplus }_{as}A\left(\left(x,y\right),\left(p.q\right)\right)=\phi \left(x,y\right)+{\phi }^{*}\left(-{A}^{*}y-p,Ax-q\right).$
The above defined convolution operation should not be confused with the standard convolution for $L$  and $M$  as convex functions in both variables. It is easy to see that in the case where $L\left(x,p\right)=\phi \left(x\right)+{\phi }^{*}\left(-p\right)$  and $M\left(x,p\right)=\psi \left(x\right)+{\psi }^{*}\left(-p\right)$  , addition corresponds to taking $\left(L\oplus M\right)\left(x,p\right)=\left(\phi +\psi \right)\left(x\right)+{\phi }^{*}◆{\psi }^{*}\left(-p\right)$  while convolution reduces to:
$\left(L◆M\right)\left(x,p\right)=\left(\phi ◆\psi \right)\left(x\right)+\left({\phi }^{*}+{\psi }^{*}\right)\left(-p\right).$  which also means that they are dual operations. We do not know whether this is true in general, but for the sequel we shall only need the following:
Lemma 3.1 Let $X$  be a reflexive Banach space and consider two Lagrangians $L$  and $M$  in $\mathcal{ℒ}\left(X\right)$  .
Then the following hold:
• 1. If $\lambda >0$  , then $\left(\lambda \cdot L\right)◆M=\lambda \cdot \left(L◆M\right)$  .
• 2. $\left(L\oplus M{\right)}^{*}\le {L}^{*}◆{M}^{*}$  and $\left(L◆M{\right)}^{*}\le {L}^{*}\oplus {M}^{*}$  .
• 3. If $L$  or $M$  is a basic ASD Lagrangian, then $\left(L\oplus M{\right)}^{*}={L}^{*}◆{M}^{*}$  and $\left(L◆M{\right)}^{*}={L}^{*}\oplus {M}^{*}$  .
• 4. If $L$  and $M$  are in ${\mathcal{ℒ}}_{AD}\left(X\right)$  , then ${L}^{*}\oplus {M}^{*}\left(q,y\right)=L◆M\left(-y,-q\right)$  for every $\left(y,q\right)\in X×{X}^{*}$  .
Proof: (1) is straightforward. To prove (2), fix $\left(q,y\right)\in {X}^{*}×X$  and write:
 $\begin{array}{ccc}& & \left(L◆M{\right)}^{*}\left(q,y\right)\end{array}$
 $\begin{array}{ccc}& =& sup\left\{〈q,x〉+〈y,p〉-L\left(z,p\right)-M\left(x-z,p\right);\left(z,x,p\right)\in X×X×{X}^{*}\right\}\end{array}$
 $\begin{array}{ccc}& =& sup\left\{〈q,v+z〉+〈y,p〉-L\left(z,p\right)-M\left(v,p\right);\left(z,v,p\right)\in X×X×{X}^{*}\right\}\end{array}$
 $\begin{array}{ccc}& =& sup\left\{〈q,v+z〉+sup\left\{〈y,p〉-L\left(z,p\right)-M\left(v,p\right);p\in {X}^{*}\right\};\left(z,v\right)\in X×X\right\}\end{array}$
 $\begin{array}{ccc}& =& {sup}_{\left(z,v\right)\in X×X}\left\{〈q,v+z〉+{inf}_{w\in X}\left\{{sup}_{{p}_{1}\in {X}^{*}}\left(〈w,{p}_{1}〉-L\left(z,{p}_{1}\right)\right)+{sup}_{{p}_{2}\in {X}^{*}}\left(〈y-w,{p}_{2}〉-M\left(v,{p}_{2}\right)\right)\right\}\right\}\end{array}$
 $\begin{array}{ccc}& \le & {inf}_{w\in X}\left\{{sup}_{\left(z,{p}_{1}\right)\in X×{X}^{*}}\left\{〈q,z〉+〈w,{p}_{1}〉-L\left(z,{p}_{1}\right)\right)\right\}+{sup}_{\left(v,{p}_{2}\right)\in X×{X}^{*}}\left\{〈q,v〉+〈y-w,{p}_{2}〉-M\left(v,{p}_{2}\right)\right\}\end{array}$
 $\begin{array}{ccc}& =& {inf}_{w\in X}\left\{{L}^{*}\left(q,w\right)+{M}^{*}\left(q,y-w\right)\right\}\end{array}$
 $\begin{array}{ccc}& =& \left({L}^{*}\oplus {M}^{*}\right)\left(q,y\right).\end{array}$
For (3) assume that $M\left(x,p\right)=\phi \left(x\right)+{\phi }^{*}\left(-p\right)$  where $\phi$  is a convex lower semi-continuous function.
Fix $\left(q,y\right)\in {X}^{*}×X$  and write:
 $\begin{array}{ccc}\left(L◆M{\right)}^{*}\left(q,y\right)& =& sup\left\{〈q,x〉+〈y,p〉-L\left(z,p\right)-M\left(x-z,p\right);\left(z,x,p\right)\in X×X×{X}^{*}\right\}\end{array}$
 $\begin{array}{ccc}& =& sup\left\{〈q,v+z〉+〈y,p〉-L\left(z,p\right)-M\left(v,p\right);\left(z,v,p\right)\in X×X×{X}^{*}\right\}\end{array}$
 $\begin{array}{ccc}& =& {sup}_{p\in {X}^{*}}\left\{〈y,p〉+{sup}_{\left(z,v\right)\in X×X}\left\{〈q,v+z〉-L\left(z,p\right)-\phi \left(v\right)\right\}\right\}-{\phi }^{*}\left(-p\right)\right\}\end{array}$
 $\begin{array}{ccc}& =& {sup}_{p\in {X}^{*}}\left\{〈y,p〉+{sup}_{z\in X}\left\{〈q,z〉-L\left(z,p\right)\right\}+{sup}_{v\in X}\left\{〈q,v〉-\phi \left(v\right)\right\}-{\phi }^{*}\left(-p\right)\right\}\end{array}$
 $\begin{array}{ccc}& =& {sup}_{p\in {X}^{*}}\left\{〈y,p〉+{sup}_{z\in X}\left\{〈q,z〉-L\left(z,p\right)\right\}+{\phi }^{*}\left(q\right)-{\phi }^{*}\left(-p\right)\right\}\end{array}$
 $\begin{array}{ccc}& =& {sup}_{p\in {X}^{*}}{sup}_{z\in X}\left\{〈y,p〉+〈q,z〉-L\left(z,p\right)-{\phi }^{*}\left(-p\right)\right\}+{\phi }^{*}\left(q\right)\end{array}$
 $\begin{array}{ccc}& =& \left(L+T{\right)}^{*}\left(q,y\right)+{\phi }^{*}\left(q\right)\end{array}$
where $T\left(z,p\right):={\phi }^{*}\left(-p\right)$  for all $\left(z,p\right)\in X×{X}^{*}$  . Note now that
 $\begin{array}{ccc}{T}^{*}\left(q,y\right)={sup}_{z,p}\left\{〈q,z〉+〈y,p〉-{\phi }^{*}\left(-p\right)\right\}=\left\{\begin{array}{ccc}+\infty & if& q\ne 0\\ \phi \left(-y\right)& if& q=0\end{array}& & \end{array}$
in such a way that by using the duality between sums and convolutions in both variables, we get
 $\begin{array}{ccc}\left(L+T{\right)}^{*}\left(q,y\right)& =& conv\left({L}^{*},{T}^{*}\right)\left(q,y\right)\end{array}$
 $\begin{array}{ccc}& =& {inf}_{r\in {X}^{*},z\in X}\left\{{L}^{*}\left(r,z\right)+{T}^{*}\left(-r+q,-z+y\right)\right\}\end{array}$
 $\begin{array}{ccc}& =& {inf}_{z\in X}\left\{{L}^{*}\left(q,z\right)+\phi \left(z-y\right)\right\}\end{array}$
 $\begin{array}{ccc}& & \end{array}$
and finally
 $\begin{array}{ccc}\left(L◆M{\right)}^{*}\left(q,y\right)& =& \left(L+T{\right)}^{*}\left(q,y\right)+{\phi }^{*}\left(q\right)\end{array}$
 $\begin{array}{ccc}& =& {inf}_{z\in X}\left\{{L}^{*}\left(q,z\right)+\phi \left(z-y\right)\right\}+{\phi }^{*}\left(q\right)\end{array}$
 $\begin{array}{ccc}& =& {inf}_{z\in X}\left\{{L}^{*}\left(q,z\right)+{\phi }^{*}\left(q\right)+\phi \left(z-y\right)\right\}\end{array}$
 $\begin{array}{ccc}& =& \left({L}^{*}\oplus {M}^{*}\right)\left(q,y\right).\end{array}$
The rest follows in the same way. For (4) write
 $\begin{array}{ccc}\left({L}^{*}\oplus {M}^{*}\right)\left(q,y\right)& =& {inf}_{w\in X}\left\{{L}^{*}\left(q,w\right)+{M}^{*}\left(q,y-w\right)\right\}\end{array}$
 $\begin{array}{ccc}& =& {inf}_{w\in X}\left\{L\left(-w,-q\right)+M\left(w-y,-q\right)\right\}\end{array}$
 $\begin{array}{ccc}& =& \left(L◆M\right)\left(-y,-q\right).\end{array}$
The following proposition summarizes some of the remarkable permanence properties of ASD Lagrangians.
Proposition 3.1 Let $X$  be a reflexive Banach space, then the following holds:
• 1. If $L$  is in ${\mathcal{ℒ}}_{AD}\left(X\right)$  , then ${L}^{*}\in {\mathcal{ℒ}}_{AD}\left({X}^{*}\right)$  , and if $\lambda >0$  , then $\lambda \cdot L$  also belong to ${\mathcal{ℒ}}_{AD}\left(X\right)$  .
• 2. If $L$  and $M$  are in ${\mathcal{ℒ}}_{AD}\left(X\right)$  and one of them is basic, then the Lagrangians $L\oplus M$  , and $L◆M$  also belong to ${\mathcal{ℒ}}_{AD}\left(X\right)$  .
• 3. If ${L}_{i}\in {\mathcal{ℒ}}_{AD}\left({X}_{i}\right)$  where ${X}_{i}$  is a reflexive Banach space for each $i\in I$  , then ${\Sigma }_{i\in I}{L}_{i}$  is in ${\mathcal{ℒ}}_{AD}\left({\Pi }_{i\in I}{X}_{i}\right)$  .
• 4. If $L\in {\mathcal{ℒ}}_{AD}\left(X\right)$  and $\Lambda :X\to {X}^{*}$  is a skew-adjoint bounded linear operator (i.e., ${\Lambda }^{*}=-\Lambda$  ), then the Lagrangian ${L}_{\Lambda }$  is also in ${\mathcal{ℒ}}_{AD}\left(X\right)$  .
• 5. If $L\in {\mathcal{ℒ}}_{AD}\left(X\right)$  and if $\Lambda :X\to {X}^{*}$  is an invertible skew-adjoint operator, then the Lagrangian ${}_{\Lambda }L$  is also in ${\mathcal{ℒ}}_{AD}\left(X\right)$  .
• 6. If $L\in {\mathcal{ℒ}}_{AD}\left(X\right)$  and $M\in {\mathcal{ℒ}}_{AD}\left(Y\right)$  , then for any bounded linear operator $A:X\to {Y}^{*}$  , the Lagrangian $L{\oplus }_{A}M$  belongs to ${\mathcal{ℒ}}_{AD}\left(X×Y\right)$
• 7. If $\phi$  is a proper convex lower semi-continuous function on $X×Y$  and $A$  is any bounded linear operator $A:X\to {Y}^{*}$  , then $\phi {\oplus }_{as}A$  belongs to ${\mathcal{ℒ}}_{AD}\left(X×Y\right)$
Proof: (1) and the stability by multiplication with a scalar is straightforward. (2) follows from the above lemma and (3) is obvious. To show (4) fix $\left(q,y\right)\in {X}^{*}×X$  , set $r=\Lambda x+p$  and write:
 $\begin{array}{ccc}{L}_{\Lambda }^{*}\left(q,y\right)& =& sup\left\{〈q,x〉+〈y,p〉-L\left(x,\Lambda x+p\right);\left(x,p\right)\in X×{X}^{*}\right\}\end{array}$
 $\begin{array}{ccc}& =& sup\left\{〈q,x〉+〈y,r-\Lambda x〉-L\left(x,r\right);\left(x,r\right)\in X×{X}^{*}\right\}\end{array}$
 $\begin{array}{ccc}& =& sup\left\{〈q+\Lambda y,x〉+〈y,r〉-L\left(x,r\right);\left(x,r\right)\in X×{X}^{*}\right\}\end{array}$
 $\begin{array}{ccc}& =& {L}^{*}\left(q+\Lambda y,y\right)=L\left(-y,-\Lambda y-q\right)\end{array}$
 $\begin{array}{ccc}& =& {L}_{\Lambda }\left(-y,-q\right).\end{array}$
For (5) let $r=x-{\Lambda }^{-1}p$  and $s=\Lambda x$  and write
 $\begin{array}{ccc}{}_{\Lambda }{L}^{*}\left(q,y\right)& =& sup\left\{〈q,x〉+〈y,p〉-L\left(x-{\Lambda }^{-1}p,\Lambda x\right);\left(x,p\right)\in X×{X}^{*}\right\}\end{array}$
 $\begin{array}{ccc}& =& sup\left\{〈q,{\Lambda }^{-1}s〉+〈y,s-\Lambda r〉-L\left(r,s\right);\left(r,s\right)\in X×{X}^{*}\right\}\end{array}$
 $\begin{array}{ccc}& =& sup\left\{〈-{\Lambda }^{-1}q+y,s〉+〈\Lambda y,r〉-L\left(r,s\right);\left(r,s\right)\in X×{X}^{*}\right\}\end{array}$
 $\begin{array}{ccc}& =& {L}^{*}\left(\Lambda y,-{\Lambda }^{-1}q+y\right)=L\left(-y+{\Lambda }^{-1}q,-\Lambda y\right)\end{array}$
 $\begin{array}{ccc}& =& {}_{\Lambda }L\left(-y,-q\right).\end{array}$
For (6), it is enough to notice that for $\left(\stackrel{~}{x},\stackrel{~}{p}\right)\in \left(X×Y\right)×\left({X}^{*}×{Y}^{*}\right)$  , we can write $L{\oplus }_{A}M\left(\stackrel{~}{x},\stackrel{~}{p}\right)=\left(L+M\right)\left(\stackrel{~}{x},\stackrel{~}{A}\stackrel{~}{x}+\stackrel{~}{p}\right)$  where $\stackrel{~}{A}:X×Y\to {X}^{*}×{Y}^{*}$  is the skew-adjoint operator defined by $\stackrel{~}{A}\left(\stackrel{~}{x}\right)=\stackrel{~}{A}\left(\left(x,y\right)\right)=\left({A}^{*}y,-Ax\right).$  Assertion (7) follows from (4) since $\phi {\oplus }_{as}A\left(\left(x,y\right),\left(p,q\right)\right)=\phi \left(x,y\right)+{\phi }^{*}\left(-{A}^{*}y-p,Ax-q\right)={L}_{\stackrel{~}{A}}\left(\left(x,y\right),\left(p,q\right)\right)$  where $L\left(\left(x,y\right),\left(p.q\right)\right):=\phi \left(x,y\right)+{\phi }^{*}\left(-p,-q\right)$  is obviously in ${\mathcal{ℒ}}_{AD}\left(X×Y\right)$  and where $\stackrel{~}{A}:X×Y\to {X}^{*}×{Y}^{*}$  is again the skew-adjoint operator defined by $\stackrel{~}{A}\left(\left(x,y\right)\right)=\left({A}^{*}y,-Ax\right).$
Remark 3.2 The proof of (4) and (5) above clearly shows that ${L}_{\Lambda }$  (resp., ${}_{\Lambda }L$  ) is partially anto-selfdual if and only if $L$  is anti-selfdual on the graph of $\Lambda$  .
Remark 3.3 An important use of the above proposition is when ${M}_{\lambda }\left(x,p\right)=\frac{\parallel x{\parallel }^{2}}{2{\lambda }^{2}}+\frac{{\lambda }^{2}\parallel p{\parallel }^{2}}{2}$  , then ${L}_{\lambda }=L◆{M}_{\lambda }$  is a $\lambda$  -regularization of the Lagrangian $L$  , which is reminescent of the Yosida theory for operators and for convex functions.
This will be most useful in [16and [20.
Remark 3.4 Denote by ${\mathcal{ℒ}}_{AD}^{+}\left(X\right)$  the cone of sub-ASD Lagrangians: i.e., those $L$  in $\mathcal{ℒ}\left(X\right)$  such that
 $\begin{array}{c}{L}^{*}\left(p,x\right)\ge L\left(-x,-p\right)forall\left(x,p\right)\in X×{X}^{*}.\end{array}$ (25)
A typical example is a Lagrangian of the form $L\left(x,p\right)=\phi \left(Cx\right)+{\phi }^{*}\left(-p\right)$  where either $C$  is a surjective operator from $X$  onto itself, or when $C$  has a dense range and $\phi$  is continuous. It is easy to see that ${\mathcal{ℒ}}_{AD}^{+}\left(X\right)$  also satisfies the following permanence properties:
• 1. If $L$  is in ${\mathcal{ℒ}}_{AD}^{+}\left(X\right)$  , $M$  is a basic $ASD$  -Lagrangian and $\lambda >0$  , then the Lagrangians $L+M$  , $L◆M$  and $\lambda \cdot L$  also belong to ${\mathcal{ℒ}}_{AD}^{+}\left(X\right)$  .
• 2. If ${L}_{i}\in {\mathcal{ℒ}}_{AD}^{+}\left({X}_{i}\right)$  where ${X}_{i}$  is a reflexive Banach space for each $i\in I$  , then ${\oplus }_{i\in I}{L}_{i}$  is in ${\mathcal{ℒ}}_{AD}^{+}\left({\Pi }_{i\in I}{X}_{i}\right)$  .
• 3. If $L\in {\mathcal{ℒ}}_{AD}^{+}\left(X\right)$  and $\Lambda :X\to {X}^{*}$  is skew-adjoint then ${L}_{\Lambda }$  is also in ${\mathcal{ℒ}}_{AD}^{+}\left(X\right)$  .
• 4. If $L\in {\mathcal{ℒ}}_{AD}^{+}\left(X\right)$  and if $\Lambda :X\to {X}^{*}$  is an invertible skew-adjoint operator, then ${}_{\Lambda }L$  is also in ${\mathcal{ℒ}}_{AD}^{+}\left(X\right)$  .
• 5. If $L\in {\mathcal{ℒ}}_{AD}^{+}\left(X\right)$  and $M\in {\mathcal{ℒ}}_{AD}^{+}\left(Y\right)$  , then for any bounded linear operator $A:X\to {Y}^{*}$  , the Lagrangian $L{\oplus }_{A}M$  belongs to ${\mathcal{ℒ}}_{AD}^{+}\left(X×Y\right)$

4 ASD Lagrangians in variational problems with no boundary constraint

An immediare corollary of Theorem 2.4 in the special case of ASD Lagrangians is the following result which will be used repeatedly in the sequel.
Theorem 4.1 Let $\Lambda :X\to {X}^{*}$  be a bounded linear skew-adjoint operator on a reflexive Banach space $X$  , and let $L$  be an anti-self dual Lagrangian on the graph of $\Lambda$  . Assume one of the following hypothesis:
• 1. (A) ${lim}_{\parallel x\parallel \to \infty }\frac{L\left(x,\Lambda x\right)}{\parallel x\parallel }=+\infty$  , or
• 2. (B) The operator $\Lambda$  is invertible and the map $x\to L\left(x,0\right)$  is bounded above on a neighborhood of the origin of $X$  .
Then there exists $\overline{x}\in X$  , such that:
 $\begin{array}{c}\left\{\begin{array}{ccc}L\left(\overline{x},\Lambda \overline{x}\right)& =& {inf}_{x\in X}L\left(x,\Lambda x\right)=0.\\ \left(-\Lambda \overline{x},-\overline{x}\right)& \in & \partial L\left(\overline{x},\Lambda \overline{x}\right).\end{array}\end{array}$ (26)
Proof: It suffices to apply Theorem 2.4 in the case where $R\left(x\right)=-x$  . In the case where $\Lambda$  is also invertible, then we directly apply Proposition 2.1 to the Lagrangian ${}_{\Lambda }L\left(x,p\right)=L\left(x+{\Lambda }^{-1}p,\Lambda x\right)$  which is partially anti-selfdual.
We note that in view of Remark 2.7, it is sufficient to have a Lagrangian $L$  in ${\mathcal{ℒ}}_{AD}^{+}\left(X\right)$  that is non-negative on the graph of $\Lambda$  , that is if $L\left(x,\Lambda x\right)\ge 0$  for all $x\in X$  .
Example 1: A variational formulation for the Lax-Milgram theorem Given a bilinear continuous functional $a$  on a Banach space $X$  , and assuming that $a$  is coercive:
i.e., for some $\lambda >0$  , we have that $a\left(v,v\right)\ge \lambda \parallel v{\parallel }^{2}$  for every $v\in X$  . It is well known that if $a$  is symmetric, then for any $f\in {X}^{*}$  , we can use a variational approach to find $u\in X$  , such that for every $v\in X$  , we have $a\left(u,v\right)=〈v,f〉$  . The procedure amounts to minimize on $H$  the convex functional $\psi \left(u\right)=\frac{1}{2}a\left(u,u\right)-〈u,f〉.$  The theorem of Lax-Milgram deals with the case when $a$  is not symmetric, for which the above variational argument does not work. Theorem 4.1 however yields the following variational formulation and proof of the original Lax-Milgram theorem.
Corollary 4.2 Let $a$  be a coercive continuous bilinear form on $X×X$  . For any $f\in {X}^{*}$  , consider the functional $I\left(v\right)=\psi \left(v\right)+{\psi }^{*}\left(-\Lambda v\right)$  where $\psi \left(v\right)=\frac{1}{2}a\left(v,v\right)-〈v,f〉$  , ${\psi }^{*}$  its Legendre conjugate and where $\Lambda :X\to {X}^{*}$  is the skew-adjoint operator defined by $〈\Lambda v,w〉=\frac{1}{2}\left(a\left(v,w\right)-a\left(w,v\right)\right).$  Then, there exists $u\in X$  , such that $I\left(u\right)={inf}_{v\in H}I\left(v\right)=0anda\left(u,v\right)=〈v,f〉foreveryv\in X.$
Proof: Consider the Lagrangian $L\left(x,p\right)=\psi \left(x\right)+{\psi }^{*}\left(-p\right)$  which is clearly anti-self dual. Apply Theorem 4.1 and note that:
$L\left(u,\Lambda u\right)=0ifandonlyif\psi \left(u\right)+{\psi }^{*}\left(-\Lambda u\right)=0=-〈\Lambda u,u〉,$  which means that $-\Lambda u\in \partial \psi \left(u\right)$  . In other words, we have for every $v\in X$  $-\frac{1}{2}\left(a\left(u,v\right)-a\left(v,u\right)\right)=\frac{1}{2}\left(a\left(u,v\right)+a\left(v,u\right)\right)-〈v,f〉$  which yields our claim.Example 2: Inverting variationally a non-selfadjoint matrix An immediate finite dimensional application of the above corollary is the following variational solution for the linear equation $Ax=y$  where $A$  is an $n×n$  -matrix and $y\in I{R}^{n}$  . It then suffices to minimize $I\left(x\right)=\frac{1}{2}〈Ax,x〉+\frac{1}{2}〈{A}_{s}^{-1}\left(y-{A}_{a}x\right),y-{A}_{a}x〉-〈y,x〉.$  on $I{R}^{n}$  , where ${A}_{a}$  is the anti-symmetric part of $A$  and ${A}_{s}^{-1}$  is the inverse of the symmetric part. If $A$  is coercive, i.e., $〈Ax,x〉\ge c|x{|}^{2}$  for all $x\in {\mathbf{R}}^{n}$  , then there is a solution $\overline{x}\in {\mathbf{R}}^{n}$  to the equation obtained as $I\left(\overline{x}\right)={inf}_{x\in {\mathbf{R}}^{n}}I\left(x\right)=0.$

4.1 ASD Lagrangians as representations of certain maximal monotone operators

As noted above, the basic examples of anti-selfdual Lagrangians are of the form
 $\begin{array}{c}L\left(x,p\right)=\phi \left(x\right)+{\phi }^{*}\left(-Bx-p\right)\end{array}$ (27)
where $\phi$  is a convex and lower semi-continuous function on $X$  , ${\phi }^{*}$  is its Legendre conjugate on ${X}^{*}$  and where $B:X\to {X}^{*}$  is skew-symmetric. This suggests that ASD Lagrangians are natural extensions of operators of the form $A+\partial \phi$  , where $A$  is positive and $\phi$  is convex. This is an important subclass of maximal monotone operators which can now be resolved variationally.
Indeed, first consider the cone $\mathcal{C}\left(X\right)$  of all bounded below, proper convex lower semi-continuous functions on $X$  , and let $\mathcal{A}\left(X\right)$  be the cone of all positive bounded linear operators from $X$  into ${X}^{*}$  (i.e., $〈Ax,x〉\ge 0$  for all $x\in X$  ). Consider also the subclasses ${\mathcal{C}}_{0}\left(X\right)=\left\{\phi \in \mathcal{C}\left(X\right);{inf}_{x\in X}\phi \left(x\right)=0\right\}and{\mathcal{A}}_{0}\left(X\right)=\left\{A\in \mathcal{A}\left(X\right);{A}^{*}=-A\right\}.$
Proposition 4.1 (1) There is a projection $\Pi :\left(\mathcal{C}\left(X\right),\mathcal{A}\left(X\right)\right)\to \left({\mathcal{C}}_{0}\left(X\right),{\mathcal{A}}_{0}\left(X\right)\right)$  such that if $\left({\phi }_{0},{A}_{0}\right)$  is the image of $\left(\phi ,A\right)$  by $\Pi$  , then a pair $\left(x,f\right)\in X×{X}^{*}$  satisfies $\left(A+\partial \phi \right)\left(x\right)=f$  if and only if $\left({A}_{0}+\partial {\phi }_{0}\right)\left(x\right)=f$  .
(2) For any pair $\left(\phi ,A\right)\in \mathcal{C}\left(X\right)×\mathcal{A}\left(X\right)$  there exists a Lagrangian ${L}_{{}_{\left(\phi ,A\right)}}\in {\mathcal{ℒ}}_{AD}\left(X\right)$  such that the equation $\left(A+\partial \phi \right)\left(x\right)=0$  has a solution $\overline{x}\in X$  if and only if the functional $I\left(x\right)={L}_{{}_{\left(\phi ,A\right)}}\left(x,0\right)$  attains its infimum.
Proof: (1) Define the projection as follows: For $\left(\phi ,A\right)\in \left(\mathcal{C}\left(X\right),\mathcal{A}\left(X\right)\right)$  , decompose $A$  into a symmetric ${A}^{s}$  and an anti-symmetric part ${A}^{a}$  , by simply writing ${A}^{s}=\frac{1}{2}\left(A+{A}^{*}\right)$  and ${A}^{a}=\frac{1}{2}\left(A-{A}^{*}\right)$  . Let ${\phi }_{0}$  be the convex functional $\psi +{\psi }^{*}\left(0\right)$  , where $\psi \left(x\right)=\frac{1}{2}〈Ax,x〉+\phi \left(x\right)$  , and define the projection as $\Pi \left(\phi ,A\right)=\left({\phi }_{0},{A}^{a}\right).$  (2) Associate to each pair $\left(\phi ,A\right)\in \mathcal{C}\left(X\right)×\mathcal{A}\left(X\right)$  , the anti-selfdual Lagrangian ${L}_{\left(\phi ,A\right)}\left(x,p\right)={L}_{\left({\phi }_{0},{A}_{a}\right)}\left(x,p\right)={\phi }_{0}\left(x\right)+{\phi }_{0}^{*}\left(-{A}^{a}x-p\right)forany\left(x,p\right)\in X×{X}^{*},$  where $\left({\phi }_{0},{A}^{a}\right)$  is the projection of $\left(\phi ,A\right)$  . The fact that the minimum of $I\left(x\right)={\phi }_{0}\left(x\right)+{\phi }_{0}^{*}\left(-{A}^{a}x\right)$  is equal to $0$  and is attained at some $\overline{x}\in X$  means that ${\phi }_{0}\left(\overline{x}\right)+{\phi }_{0}^{*}\left(-{A}^{a}\overline{x}\right)=0=-〈{A}^{a}\overline{x},\overline{x}〉$  which yields, in view of Legendre-Fenchel duality that $-{A}^{a}\overline{x}\in \partial {\phi }_{0}\left(\overline{x}\right)={A}^{s}\overline{x}+\partial \phi \left(\overline{x}\right)$  , hence $\overline{x}$  satisfies $-Ax\in \partial \phi \left(x\right)$  .
Remark 4.3 We note the following relations between classical operations on functions and operators and the operations on ASD Lagrangians.
• For $\lambda >0$  and $\phi \in \mathcal{C}\left(X\right)$  , we have $\lambda \cdot {L}_{\left(\phi ,A\right)}={L}_{\left({\lambda }^{2}\phi \left(\frac{\cdot }{\lambda }\right),A\right)}$  .
• ${L}_{\left({\phi }_{1},{A}_{1}\right)}+{L}_{\left({\phi }_{2},{A}_{2}\right)}={L}_{\left({\phi }_{1}+{\phi }_{2},{A}_{1}+{A}_{2}\right)}$  .
• ${L}_{\left({\phi }_{\lambda },0\right)}={L}_{\left(\phi ,0\right)}◆{L}_{\left(\frac{\parallel x{\parallel }^{2}}{2{\lambda }^{2}},0\right)}$  where ${\phi }_{\lambda }$  is the Yosida regularization of $\phi$  .
• More generally, ${L}_{\left({\phi }_{1},0\right)}◆{L}_{\left({\phi }_{2},0\right)}={L}_{\left({\phi }_{1}◆{\phi }_{2},0\right)}$
In the sequel, whenever $\phi$  is a functional on $X$  and $f\in {X}^{*}$  , we shall denote by $\phi +f$  the functional defined for $x\in X$  by $\phi \left(x\right)+〈f,x〉$  . Now we can a variational resolution to the following nonlinear Lax-Milgram type result.
Corollary 4.4 Assume one of the following conditions on a pair $\left(\phi ,A\right)\in \mathcal{C}\left(X\right)×\mathcal{A}\left(X\right)$  :
• 1. (A) ${lim}_{\parallel x\parallel \to \infty }\parallel x{\parallel }^{-1}\left(\phi \left(x\right)+\frac{1}{2}〈Ax,x〉\right)=+\infty$  , or
• 2. (B) The operator ${A}^{a}=\frac{1}{2}\left(A-{A}^{*}\right):X\to {X}^{*}$  is onto and $\phi$  is bounded above on the bounded sets of $X$  .
Then, there exists for any $f\in {X}^{*}$  , a solution $\overline{x}\in X$  to the equation $-Ax+f\in \partial \phi \left(x\right)$  that can be obtained as a minimizer of the problem:
 $\begin{array}{c}{inf}_{x\in X}\left\{\psi \left(x\right)+{\psi }^{*}\left(-{A}^{a}x\right)\right\}=0\end{array}$ (28)
where $\psi$  is the convex functional $\psi \left(x\right)=\frac{1}{2}〈Ax,x〉+\phi \left(x\right)-〈f,x〉$  .
Proof: Associate to $\left(\psi ,A\right)\in \mathcal{C}\left(X\right)×\mathcal{A}\left(X\right)$  , the anti-selfdual Lagrangian ${L}_{\left(\psi ,A\right)}\left(x,p\right)={L}_{\left(\psi ,{A}^{a}\right)}\left(x,p\right)=\psi \left(x\right)+{\psi }^{*}\left(-{A}^{a}x-p\right),for\left(x,p\right)\in X×{X}^{*},$  The fact that the minimum in ( 28 ) is attained at some $\overline{x}\in X$  , follows from Theorem 4.1. It means that $\psi \left(\overline{x}\right)+{\psi }^{*}\left(-{A}^{a}\overline{x}\right)=-〈{A}^{a}\overline{x},\overline{x}〉$  which yields, in view of Legendre-Fenchel duality that $-{A}^{a}\overline{x}\in \partial \psi \left(\overline{x}\right)={A}^{s}\overline{x}+\partial \phi \left(\overline{x}\right)-f$  , hence $\overline{x}$  satisfies $-Ax+f\in \partial \phi \left(x\right)$  .
Remark 4.5 All what is needed in the above proposition is that the function $\phi$  be $A$  -convex for some operator $A$  , meaning that $\psi \left(x\right)=\frac{1}{2}〈Ax,x〉+\phi \left(x\right)$  is convex and lower semi-continuous.
Example 3: A variational principle for a non-symmetric Dirichlet problem Let $\mathbf{a}:\Omega \to {\mathbf{R}}^{\mathbf{n}}$  be a smooth function on a bounded domain $\Omega$  of ${\mathbf{R}}^{\mathbf{n}}$  , and consider the first order linear operator $Av=\mathbf{a}\cdot \nabla v={\Sigma }_{i=1}^{n}{a}_{i}\frac{\partial v}{\partial {x}_{i}}$  Assume that the vector field ${\Sigma }_{i=1}^{n}{a}_{i}\frac{\partial v}{\partial {x}_{i}}$  is actually the restriction of a smooth vector field ${\Sigma }_{i=1}^{n}{\overline{a}}_{i}\frac{\partial v}{\partial {x}_{i}}$  defined on an open neighborhood $X$  of $\overline{\Omega }$  and that each ${\overline{a}}_{i}$  is a ${C}^{1,1}$  function on $X$  . Consider the Dirichlet problem:
 $\begin{array}{c}\left\{\begin{array}{ccc}\Delta u+{\Sigma }_{i=1}^{n}{a}_{i}\frac{\partial u}{\partial {x}_{i}}& =& |u{|}^{p-2}u+fon\Omega \\ u& =& 0on\partial \Omega .\\ \end{array}\end{array}$ (29)
If ${a}_{i}=0$  , then to find a solution, it is sufficient to minimize the functional $\Phi \left(u\right)=\frac{1}{2}{\int }_{\Omega }|\nabla u{|}^{2}dx+\frac{1}{p}{\int }_{\Omega }|u{|}^{p}dx+{\int }_{\Omega }fudx$  and get the solution of $\partial \Phi \left(u\right)=0$  .
However, if the non self-adjoint term $a$  is not zero, we can use the above to get
Theorem 4.6 Assume $div\left(\mathbf{a}\right)\ge 0$  on $\Omega$  , and consider on ${H}_{0}^{1}\left(\Omega \right)$  , the functional $I\left(u\right)=\Psi \left(u\right)+{\Psi }^{*}\left(\mathbf{a}.\nabla u+\frac{1}{2}div\left(\mathbf{a}\right)u\right)$  where $\Psi \left(u\right)=\frac{1}{2}{\int }_{\Omega }|\nabla u{|}^{2}dx+\frac{1}{p}{\int }_{\Omega }|u{|}^{p}dx+{\int }_{\Omega }fudx+\frac{1}{4}{\int }_{\Omega }div\left(\mathbf{a}\right)|u{|}^{2}dx$  and ${\Psi }^{*}$  is its Legendre transform. Then, there exists $\overline{u}\in {H}_{0}^{1}\left(\Omega \right)$  such that:
$I\left(\overline{u}\right)=inf\left\{I\left(u\right);u\in {H}_{0}^{1}\left(\Omega \right)\right\}=0,$  and $\overline{u}$  is a solution of  29  $\left(\right)$  .
Proof: Indeed, $\Psi$  is clearly convex and lower semi-continuous on ${H}_{0}^{1}\left(\Omega \right)$  while the operator $\Lambda u=-\mathbf{a}.\nabla u-\frac{1}{2}div\left(\mathbf{a}\right)u$  is skew-adjoint, since ${\int }_{\Omega }\left(\mathbf{a}.\nabla u\right)u+\frac{1}{2}div\left(\mathbf{a}\right){u}^{2}dx=0.$  Again the functional $I\left(u\right)=\Psi \left(u\right)+{\Psi }^{*}\left(\mathbf{a}.\nabla u+\frac{1}{2}div\left(\mathbf{a}\right)u$  is given by a self-dual Lagrangian $L\left(u,\Lambda u\right)$  where $L\left(u,v\right)=\Psi \left(u\right)+{\Psi }^{*}\left(v\right)$  is defined on ${H}_{0}^{1}\left(\Omega \right)×{H}^{-1}\left(\Omega \right)$  . The existence follows from Theorem 4.1, since $\Psi$  is clearly coercive. Note that $\overline{u}$  then satisfies $\mathbf{a}.\nabla \overline{u}+\frac{1}{2}div\left(\mathbf{a}\right)\overline{u}=\partial \Psi \left(\overline{u}\right)=-\Delta \overline{u}+{\overline{u}}^{p-1}+f+\frac{1}{2}div\left(\mathbf{a}\right)\overline{u}$  and therefore $\overline{u}$  is a solution for  29  $\left(\right)$  .Example 4: A variational solution for variational inequalities Given again a bilinear continuous functional $a$  on $X×X$  , and $\phi :X\to \mathbf{R}$  a convex l.s.c, then solving the corresponding variational inequality amounts to constructing for any $f\in {X}^{*}$  , a point $y\in X$  such that for all $z\in X$  ,
 $\begin{array}{c}a\left(y,y-z\right)+\phi \left(y\right)-\phi \left(z\right)\le 〈y-z,f〉.\end{array}$ (30)
It is well known that this problem can be rewritten as $f\in Ay+\partial \phi \left(y\right)$  where $A$  is the bounded linear operator from $X$  into ${X}^{*}$  defined by $a\left(u,v\right)=〈Au,v〉$  . This means that the variational inequality ( 30 ) can be rewritten and solved using the variational principle (21).
For example, we can solve variationally the following “obstacle” problem.
Corollary 4.7 Let $a$  be bilinear continuous functional $a$  on a reflexive Banach space $X×X$  so that $a\left(v,v\right)\ge \lambda \parallel v{\parallel }^{2}$  , and let $K$  be a convex closed subset of $X$  . Then, for any $f\in {X}^{*}$  , there is $\overline{x}\in K$  such that
 $\begin{array}{c}a\left(\overline{x},\overline{x}-z\right)\le 〈\overline{x}-z,f〉forallz\in K.\end{array}$ (31)
The point $\overline{x}$  can be obtained as a minimizer of the following problem:
${inf}_{x\in X}\left\{\phi \left(x\right)+\left(\phi +{\psi }_{K}{\right)}^{*}\left(-\Lambda x\right)\right\}=0$  where $\phi \left(u\right)=\frac{1}{2}a\left(u,u\right)-〈f,x〉$  , $\Lambda :X\to {X}^{*}$  is the skew-adjoint operator defined by $〈\Lambda u,v〉=\frac{1}{2}\left(a\left(u,v\right)-a\left(v,u\right)\right).$  and where ${\psi }_{K}\left(x\right)=0$  on $K$  and $+\infty$  elsewhere.

4.2 ASD Lagrangians and anti-Hamiltonian systems

Recall that an important class of Hamiltonian systems can be written as $\left({A}^{*}y,Ax\right)\in \partial H\left(x,y\right)$  where $A:X\to {Y}^{*}$  is a –normally symmetric– operator and $H$  is a convex (Hamiltonian) on $X×Y$  . The next proposition show however that the theory of ASD-Lagrangians is more suited for “Anti-Hamiltonian” systems of the form $\left(-{A}^{*}y,Ax\right)\in \partial H\left(x,y\right).$
Proposition 4.2 Let $\phi$  be any coercive and proper convex lower semi-continuous function on $X×Y$  with $\left(0,0\right)\in dom\left(\phi \right)$  , and let $A:X\to {Y}^{*}$  be any bounded linear operator. Assume ${B}_{1}:X\to {X}^{*}$  (resp., ${B}_{2}:Y\to {Y}^{*}$  ) are skew-adjoint operators, then there exists $\left(\overline{x},\overline{y}\right)\in X×Y$  such that
 $\begin{array}{c}\left(-{A}^{*}\overline{y}+{B}_{1}\overline{x},A\overline{x}+{B}_{2}\overline{y}\right)\in \partial \phi \left(\overline{x},\overline{y}\right).\end{array}$ (32)
The solution is obtained as a minimizer on $X×Y$  of the functional $I\left(x,y\right)=\phi \left(x,y\right)+{\phi }^{*}\left(-{A}^{*}y+{B}_{1}x,Ax+{B}_{2}y\right).$
Proof: It is enough to apply Theorem 4.1 to the ASD Lagrangian $L\left(\left(x,y\right),\left(p,q\right)\right)=\phi \left(x,y\right)+{\phi }^{*}\left(-{A}^{*}y+{B}_{1}x-p,Ax+{B}_{2}y-q\right).$  obtained by shifting to the right the ASD Lagrangian $\phi {\oplus }_{as}A$  by the skew-adjoint operator $\left(-{B}_{1},-{B}_{2}\right)$  . This yields that $I\left(x,y\right)=L\left(\left(x,y\right),\left(0,0\right)\right)$  attains its minimum at some $\left(\overline{x},\overline{y}\right)\in X×Y$  and that the minimum is actually $0$  . In other words,
 $\begin{array}{ccc}0& =& I\left(\overline{x},\overline{y}\right)=\phi \left(\overline{x},\overline{y}\right)+{\phi }^{*}\left(-{A}^{*}\overline{y}+{B}_{1}\overline{x},A\overline{x}+{B}_{2}\overline{y}\right)\end{array}$
 $\begin{array}{ccc}& =& \phi \left(\overline{x},\overline{y}\right)+{\phi }^{*}\left(-{A}^{*}\overline{y}+{B}_{1}\overline{x},A\overline{x}+{B}_{2}\overline{y}\right)-〈\left(\overline{x},\overline{y}\right),\left(-{A}^{*}\overline{y}+{B}_{1}\overline{x},A\overline{x}+{B}_{2}\overline{y}\right)〉\end{array}$
from which the equation follows.
Corollary 4.8 Given positive operators ${B}_{1}:X\to {X}^{*}$  , ${B}_{2}:Y\to {Y}^{*}$  and convex functions ${\phi }_{1}$  in $\mathcal{C}\left(X\right)$  and ${\phi }_{2}$  in $\mathcal{C}\left(Y\right)$  having $0$  in their respective domains, we consider the convex functionals ${\psi }_{1}\left(x\right)=\frac{1}{2}〈{B}_{1}x,x〉+{\phi }_{1}\left(x\right)$  and ${\psi }_{2}\left(x\right)=\frac{1}{2}〈{B}_{2}x,x〉+{\phi }_{2}\left(x\right)$  . Let $A:X\to {Y}^{*}$  be a bounded linear operator such that ${lim}_{\parallel x\parallel +\parallel y\parallel \to \infty }\frac{{\psi }_{1}\left(x\right)+{\psi }_{2}\left(y\right)}{\parallel x\parallel +\parallel y\parallel }=+\infty ,$  Then, for any $\left(f,g\right)\in {X}^{*}×{Y}^{*}$  and any $c\in IR$  , there exists a solution $\left(\overline{x},\overline{y}\right)\in X×Y$  to the system of equations
 $\begin{array}{c}\left\{\begin{array}{ccc}-{A}^{*}y-{B}_{1}x+f& \in & \partial {\phi }_{1}\left(x\right)\\ {c}^{2}Ax-{B}_{2}y+g& \in & \partial {\phi }_{2}\left(y\right).\end{array}\end{array}$ (33)
It can be obtained as a minimizer of the problem:
 $\begin{array}{c}{inf}_{x,y\in X×Y}\left\{{\chi }_{1}\left(x\right)+{\chi }_{1}^{*}\left(-{B}_{1}^{a}x-{A}^{*}y\right)+{\chi }_{2}\left(y\right)+{\chi }_{2}^{*}\left(-{B}_{2}^{a}y+{c}^{2}Ax\right)\right\}=0\end{array}$ (34)
where ${B}_{1}^{a}$  (resp., ${B}_{2}^{a}$  ) are the skew-symmetric parts of ${B}_{1}$  and ${B}_{2}$  and where ${\chi }_{1}\left(x\right)={\psi }_{1}\left(x\right)-〈f,x〉$  and ${\chi }_{2}\left(x\right)={\psi }_{2}\left(x\right)-〈g,x〉$
Proof: This follows by applying the above proposition to the convex function $\phi \left(x,y\right)={\chi }_{1}\left(x\right)+{\chi }_{2}\left(y\right)$  and the skew-symmetric operators $-{B}_{1}^{a}$  and $-{B}_{2}^{a}$  . Note that the operator $\stackrel{~}{A}:X×Y\to {X}^{*}×{Y}^{*}$  defined by $\stackrel{~}{A}\left(x,y\right)=\left({A}^{*}y,-{c}^{2}Ax\right)$  is skew adjoint once we equip $X×Y$  with the scalar product $〈\left(x,y\right),\left(p,q\right)〉=〈x,p〉+{c}^{-2}〈y,q〉.$  We then get
 $\begin{array}{c}\left\{\begin{array}{ccc}-{A}^{*}y-{B}_{1}^{a}x+f& \in & \partial {\phi }_{1}\left(x\right)+{B}_{1}^{s}\left(x\right)\\ {c}^{2}Ax-{B}_{2}^{a}y+g& \in & \partial {\phi }_{2}\left(y\right)+{B}_{2}^{s}\left(y\right).\end{array}\end{array}$ (35)
which gives the result.
Another approach consists of associating to the pairs $\left({\phi }_{1},{B}_{1}\right)$  and $\left({\phi }_{2},{B}_{1}\right)$  the anti-selfdual Lagrangians $L\left(x,p\right)={\psi }_{1}\left(x\right)-〈f,x〉+{\psi }_{1}^{*}\left(-{B}_{1}^{a}x+f-p\right),for\left(x,p\right)\in X×{X}^{*},$  and $M\left(y,q\right)={\psi }_{2}\left(y\right)-〈g,y〉+{\psi }_{2}^{*}\left(-{B}_{2}^{a}y+g-q\right),for\left(y,q\right)\in Y×{Y}^{*},$  Now apply Theorem 4.1 to the twisted-sum Lagrangian $L{\oplus }_{A}M$  .
Example 5: A variational principle for coupled equations Let ${\mathbf{b}}_{1}:\Omega \to {\mathbf{R}}^{\mathbf{n}}$  and ${\mathbf{b}}_{2}:\Omega \to {\mathbf{R}}^{\mathbf{n}}$  be two smooth vector fields on a bounded domain $\Omega$  of ${\mathbf{R}}^{\mathbf{n}}$  , verifying the conditions in example 3 and let ${B}_{1}v={\mathbf{b}}_{1}\cdot \nabla v$  and ${B}_{2}v={\mathbf{b}}_{2}\cdot \nabla v$  be the corresponding first order linear operators. Consider the Dirichlet problem:
 $\begin{array}{c}\left\{\begin{array}{ccc}\Delta \left(v+u\right)+{\mathbf{b}}_{1}\cdot \nabla u& =& {u}^{p-1}+fon\Omega \\ \Delta \left(v-{c}^{2}u\right)+{\mathbf{b}}_{2}\cdot \nabla v& =& {v}^{q-1}+gon\Omega \\ u=v& =& 0on\partial \Omega .\\ \end{array}\end{array}$ (36)
We can use the above to get
Theorem 4.9 Assume $div\left({\mathbf{b}}_{1}\right)\ge 0$  and $div\left({\mathbf{b}}_{2}\right)\ge 0$  on $\Omega$  , $1  and consider on ${H}_{0}^{1}\left(\Omega \right)×{H}_{0}^{1}\left(\Omega \right)$  the functional $I\left(u,v\right)=\Psi \left(u\right)+{\Psi }^{*}\left({\mathbf{b}}_{1}.\nabla u+\frac{1}{2}div\left({\mathbf{b}}_{1}\right)u+\Delta v\right)+\Phi \left(v\right)+{\Phi }^{*}\left({\mathbf{b}}_{2}.\nabla v+\frac{1}{2}div\left({\mathbf{b}}_{2}\right)v-{c}^{2}\Delta u\right)$  where $\Psi \left(u\right)=\frac{1}{2}{\int }_{\Omega }|\nabla u{|}^{2}dx+\frac{1}{p}{\int }_{\Omega }|u{|}^{p}dx+{\int }_{\Omega }fudx+\frac{1}{4}{\int }_{\Omega }div\left({\mathbf{b}}_{1}\right)|u{|}^{2}dx,$  $\Phi \left(v\right)=\frac{1}{2}{\int }_{\Omega }|\nabla v{|}^{2}dx+\frac{1}{q}{\int }_{\Omega }|v{|}^{q}dx+{\int }_{\Omega }gvdx+\frac{1}{4}{\int }_{\Omega }div\left({\mathbf{b}}_{2}\right)|v{|}^{2}dx$  and ${\Psi }^{*}$  and ${\Phi }^{*}$  are their Legendre transforms. Then there exists $\left(\overline{u},\overline{v}\right)\in {H}_{0}^{1}\left(\Omega \right)×{H}_{0}^{1}\left(\Omega \right)$  such that:
$I\left(\overline{u},\overline{v}\right)=inf\left\{I\left(u,v\right);\left(u,v\right)\in {H}_{0}^{1}\left(\Omega \right)×{H}_{0}^{1}\left(\Omega \right)\right\}=0,$  and $\left(\overline{u},\overline{v}\right)$  is a solution of  36  $\left(\right)$  .
We can also reduce general minimization problems of functionals of the form $I\left(x\right)=\phi \left(x\right)+\psi \left(Ax\right)$  to the much easier problem of minimizing ASD Lagrangians. Indeed we have
Proposition 4.3 Let $\phi$  (resp., $\psi$  ) be a convex lower semi-continuous function on a reflexive Banach space $X$  (resp. ${Y}^{*}\right)$  and let $A:X\to {Y}^{*}$  be a bounded linear operator. To minimize the functional $I\left(x\right)=\phi \left(x\right)+\psi \left(Ax\right)$  on $X$  , we consider on $X×Y$  the functional $I\left(x,y\right)=\phi \left(x\right)+{\psi }^{*}\left(y\right)+{\phi }^{*}\left(-{A}^{*}y\right)+\psi \left(Ax\right).$  Assuming ${lim}_{\parallel x\parallel +\parallel y\parallel \to \infty }I\left(x,y\right)=+\infty$  , then the infimum of $I$  is zero and is attained at a point $\left(\overline{x},\overline{y}\right)$  which determines the extremals of the min-max problem:
$sup\left\{-{\psi }^{*}\left(y\right)-{\phi }^{*}\left(-{A}^{*}y\right)=inf\left\{\phi \left(x\right)+\psi \left(Ax\right);x\in X\right\}.$
They also satisfy the system:
 $\begin{array}{c}\left\{\begin{array}{ccc}-{A}^{*}y& \in & \partial \phi \left(x\right)\\ Ax& \in & \partial {\psi }^{*}\left(y\right).\end{array}\end{array}$ (37)
Proof: It is sufficient to note that $I\left(x,y\right)=L\left(\left(x,y\right),\left(0,0\right)$  where $L$  is an anti-self dual Lagrangian defined on $X×Y$  by:
$L\left(\left(x,y\right),\left(p,q\right)\right)=\phi \left(x\right)+{\psi }^{*}\left(y\right)+{\phi }^{*}\left(-{A}^{*}y-p\right)+\psi \left(Ax-q\right).$  By considering more general twisted sum Lagrangians, we obtain the following application
Theorem 4.10 Let $X$  and $Y$  be two reflexive Banach spaces and let $A:X\to {Y}^{*}$  be any bounded linear operator. Assume $L\in {\mathcal{ℒ}}_{AD}\left(X\right)$  and $M\in {\mathcal{ℒ}}_{AD}\left(Y\right)$  are such that ${lim}_{\parallel x\parallel +\parallel y\parallel \to \infty }\frac{L\left(x,{A}^{*}y\right)+M\left(y,-Ax\right)}{\parallel x\parallel +\parallel y\parallel }=+\infty ,$  Then there exists $\left(\overline{x},\overline{y}\right)\in X×Y$  , such that:
 $\begin{array}{c}L\left(\overline{x},{A}^{*}\overline{y}\right)+M\left(\overline{y},-A\overline{x}\right)={inf}_{\left(x,y\right)\in X×Y\right)}L\left(x,{A}^{*}y\right)+M\left(y,-Ax\right)=0.\end{array}$ (38)
Moreover, we have
 $\begin{array}{c}\left\{\begin{array}{ccc}L\left(\overline{x},{A}^{*}\overline{y}\right)+〈\overline{x},{A}^{*}\overline{y}〉& =& 0\\ M\left(\overline{y},-A\overline{x}\right)+〈\overline{y},-Ax〉& =& 0\\ \left(-{A}^{*}\overline{y},-\overline{x}\right)& \in & \partial L\left(\overline{x},{A}^{*}\overline{x}\right)\\ \left(A\overline{x},-\overline{y}\right)& \in & \partial M\left(\overline{y},-A\overline{x}\right)\end{array}\end{array}$ (39)
Proof: It is sufficient to apply Theorem 4.1 to the ASD Lagrangian $L{\oplus }_{A}M$  .

5 ASD Lagrangians associated to boundary value problems

For problems involving boundaries, we may start with an ASD Lagrangian $L$  , but the operator $\Lambda$  may be skew-adjoint modulo a term involving the boundary. Assuming we can represent this term by a pair of operators $\left({b}_{1},{b}_{2}\right)$  from $X$  into a Hilbert space ${H}_{1}×{H}_{2}$  which correspond to an adequate splitting of the boundary, then we may try to recover anti-selfduality by adding a correcting term via a boundary Lagrangian $\ell$  . In this section, we look into frameworks where Lagrangians of the form $M\left(x,p\right)=L\left(x,\Lambda x+p\right)+\ell \left({b}_{1}\left(x\right),{b}_{2}\left(x\right)\right)$  can be made anti-selfdual.

5.1 Anti-selfduality involving boundary Lagrangians

Definition 5.1 (1) A boundary operator will be any surjective continuous linear map $\left({b}_{1},{b}_{2}\right):X\to {H}_{1}×{H}_{2}$  from $X$  onto the product of Hilbert spaces ${H}_{1}×{H}_{2}$  .
(2) An operator $\Lambda :X\to {X}^{*}$  is said to be skew-symmetric modulo the boundary operator $\left({b}_{1},{b}_{2}\right)$  , if for every $x,y\in X$  ,
 $\begin{array}{c}〈\Lambda x,y{〉}_{{}_{\left(X,{X}^{*}\right)}}=-〈\Lambda y,x{〉}_{{}_{\left(X,{X}^{*}\right)}}+〈{b}_{2}\left(x\right),{b}_{2}\left(y\right){〉}_{{}_{{H}_{2}}}-〈{b}_{1}\left(x\right),{b}_{1}\left(y\right){〉}_{{}_{{H}_{1}}}\end{array}$ (40)
That is, ${\Lambda }^{*}=-\Lambda +{b}_{2}^{*}{b}_{2}-{b}_{1}^{*}{b}_{1}$  which means that the operator $\Lambda -\frac{1}{2}\left({b}_{2}^{*}{b}_{2}-{b}_{1}^{*}{b}_{1}\right)$  is skew symmetric.
We shall then say that we have a skew symmetric triplet $\left(\Lambda ,{b}_{1},{b}_{2}\right)$  .
We also consider a Boundary Lagrangian $\ell :{H}_{1}×{H}_{2}\to IR\cup \left\{+\infty \right\}$  which is also proper convex and lower semi-continuous, and its Legendre transform in both variable, ${\ell }^{*}\left({h}_{1},{h}_{2}\right)=sup\left\{〈{k}_{1},{h}_{1}〉+〈{k}_{2},{h}_{2}〉-\ell \left({k}_{1},{k}_{2}\right);{k}_{1}\in {H}_{1},{k}_{2}\in {H}_{2}\right\}$
Definition 5.2 We say that $\ell$  is a self-dual boundary Lagrangian if
 $\begin{array}{c}{\ell }^{*}\left(-{h}_{1},{h}_{2}\right)=\ell \left({h}_{1},{h}_{2}\right)forall\left({h}_{1},{h}_{2}\right)\in {H}_{1}×{H}_{2}.\end{array}$ (41)
It is easy to see that such a boundary Lagrangian will always satisfy the inequality
 $\begin{array}{c}\ell \left(r,s\right)\ge \frac{1}{2}\left(\parallel s{\parallel }^{2}-\parallel r{\parallel }^{2}\right)forall\left(r,s\right)\in {H}_{1}×{H}_{2}.\end{array}$ (42)
The basic example of a self dual boundary Lagrangian is given by a function $\ell$  on ${H}_{1}×{H}_{2}$  , of the form $\ell \left(r,s\right)={\psi }_{1}\left(r\right)+{\psi }_{2}\left(s\right)$  , with ${\psi }_{1}^{*}\left(r\right)={\psi }_{1}\left(-r\right)$  and ${\psi }_{2}^{*}\left(s\right)={\psi }_{2}\left(s\right)$  . Here the choices for ${\psi }_{1}$  and ${\psi }_{2}$  are rather limited and the typical sample is:
${\psi }_{1}\left(r\right)=\frac{1}{2}\parallel r{\parallel }^{2}-2〈a,r〉+\parallel a{\parallel }^{2},and{\psi }_{2}\left(s\right)=\frac{1}{2}\parallel s{\parallel }^{2}.$  where $a$  is given in ${H}_{1}$  .
Boundary operators allow us to build new ASD Lagrangians. We shall present several ways to do so, which correspond to various conditions that $\Lambda$  and the Lagrangian $L$  may or may not satisfy in applications.
Proposition 5.1 Let $\ell$  be a self dual boundary Lagrangian on the Hilbertian product ${H}_{1}×{H}_{2}$  , and let $\left(\Lambda ,{b}_{1},{b}_{2}\right):X\to {X}^{*}×{H}_{1}×{H}_{2}$  be a skew symmetric triplet where $X$  is a reflexive Banach space.
Suppose $H$  is a linear subspace of ${X}^{*}$  containing $Range\left(\Lambda \right)$  such that ${X}_{0}=Ker\left({b}_{1},{b}_{2}\right)$  is dense in $X$  for the $\sigma \left(X,H\right)$  -topology, and consider $L$  to be a Lagrangian on $X$  such that for each $p\in {X}^{*}$  , the map $x\to L\left(x,p\right)$  is continuous for the $\sigma \left(X,H\right)$  -topology.
• 1. If $L$  is anti-self dual on the graph of $\Lambda$  , then the Lagrangian $M\left(x,p\right)=L\left(x,\Lambda x+p\right)+\ell \left({b}_{1}\left(x\right),{b}_{2}\left(x\right)\right)$  is partially anti-self dual.
• 2. If $L$  is anti-selfdual on the elements of $X×H$  , then $M$  is also anti-self dual on the elements of $X×H$  .
Proof: Proof: Fix $\left(q,y\right)\in {X}^{*}×X$  and calculate
 $\begin{array}{ccc}{M}^{*}\left(q,y\right)& =& sup\left\{〈q,x{〉}_{X}+〈y,p{〉}_{X}-M\left(x,p\right);\left(x,p\right)\in X×{X}^{*}\right\}\end{array}$
 $\begin{array}{ccc}& =& sup\left\{〈q,x{〉}_{X}+〈y,p{〉}_{X}-L\left(x,\Lambda x+p\right)-\ell \left({b}_{1}\left(x\right),{b}_{2}\left(x\right)\right);\left(x,p\right)\in X×{X}^{*}\right\}\end{array}$
Setting $r=\Lambda x+p$  , we obtain
 $\begin{array}{ccc}{M}^{*}\left(q,y\right)& =& sup\left\{〈x,q〉+〈y,r-\Lambda x〉-L\left(x,r\right)-\ell \left({b}_{1}\left(x\right),{b}_{2}\left(x\right)\right);\left(x,r\right)\in X×{X}^{*}\right\}\end{array}$
 $\begin{array}{ccc}& =& sup\left\{〈x,q〉+〈{b}_{1}\left(y\right),{b}_{1}\left(x\right)〉-〈{b}_{2}\left(y\right),{b}_{2}\left(x\right)〉+〈\Lambda y,x〉+〈y,r〉\end{array}$
 $\begin{array}{ccc}& & -L\left(x,r\right)-\ell \left({b}_{1}\left(x\right),{b}_{2}\left(x\right)\right);\left(x,r\right)\in X×{X}^{*}\right\}\end{array}$
 $\begin{array}{ccc}& =& sup\left\{〈x,q+\Lambda y〉+〈y,r〉-L\left(x,r\right)\end{array}$
 $\begin{array}{ccc}& & +〈{b}_{1}\left(y\right),{b}_{1}\left(x\right)〉-〈{b}_{2}\left(y\right),{b}_{2}\left(x\right)〉-\ell \left({b}_{1}\left(x\right),{b}_{2}\left(x\right)\right);\left(x,r\right)\in X×{X}^{*}\right\}\end{array}$
 $\begin{array}{ccc}& =& sup\left\{〈x,q+\Lambda y〉+〈y,r〉-L\left(x,r\right)+〈{b}_{1}\left(y\right),{b}_{1}\left(x+{x}_{0}\right)〉-〈{b}_{2}\left(y\right),{b}_{2}\left(x+{x}_{0}\right)〉\end{array}$
 $\begin{array}{ccc}& & -\ell \left({b}_{1}\left(x+{x}_{0}\right),{b}_{2}\left(x\right)+{x}_{0}\right);\left({x}_{0},x,r\right)\in {X}_{0}×X×{X}^{*}\right\}\end{array}$
 $\begin{array}{ccc}& =& sup\left\{〈w-{x}_{0},q+\Lambda y〉+〈y,r〉-L\left(w-{x}_{0},r\right)+〈{b}_{1}\left(y\right),{b}_{1}\left(w\right)〉-〈{b}_{2}\left(y\right),{b}_{2}\left(w\right)〉\end{array}$
 $\begin{array}{ccc}& & -\ell \left({b}_{1}\left(w\right),{b}_{2}\left(w\right);\left({x}_{0},w,r\right)\in {X}_{0}×X×{X}^{*}\right\}\end{array}$
Now suppose $\left(q,y\right)\in H×X$  , and use the fact that ${X}_{0}$  is $\sigma \left(X,H\right)$  dense in $X$  , that $Range\Lambda \subset H$  and the continuity of $x\to L\left(x,p\right)$  in that topology to obtain
 $\begin{array}{ccc}{M}^{*}\left(q,y\right)& =& sup\left\{〈z,q+\Lambda y〉+〈y,r〉-L\left(z,r\right)+〈{b}_{1}\left(y\right),{b}_{1}\left(w\right)〉-〈{b}_{2}\left(y\right),{b}_{2}\left(w\right)〉\end{array}$
 $\begin{array}{ccc}& & -\ell \left({b}_{1}\left(w\right),{b}_{2}\left(w\right);\left(z,w,r\right)\in X×X×{X}^{*}\right\}\end{array}$
 $\begin{array}{ccc}& =& sup\left\{〈z,q+\Lambda y〉+〈y,r〉-L\left(z,r\right);\left(z,r\right)\in X×{X}^{*}\right\}\end{array}$
 $\begin{array}{ccc}& & +sup\left\{〈{b}_{1}\left(y\right),{b}_{1}\left(w\right)〉-〈{b}_{2}\left(y\right),{b}_{2}\left(w\right)〉-\ell \left({b}_{1}\left(w\right),{b}_{2}\left(w\right)\right);w\in X\right\}\end{array}$
 $\begin{array}{ccc}& =& sup\left\{〈z,q+\Lambda y〉+〈y,r〉-L\left(z,r\right);\left(z,r\right)\in X×{X}^{*}\right\}\end{array}$
 $\begin{array}{ccc}& & +sup\left\{〈{b}_{1}\left(y\right),a〉-〈{b}_{2}\left(y\right),b〉-\ell \left(a,b\right);\left(a,b\right)\in {H}_{1}×{H}_{2}\right\}\end{array}$
 $\begin{array}{ccc}& =& {L}^{*}\left(q+\Lambda y,y\right)+{\ell }^{*}\left({b}_{1}\left(y\right),-{b}_{2}\left(y\right)\right)\end{array}$
 $\begin{array}{ccc}& =& L\left(-y,-q-\Lambda y\right)+\ell \left(-{b}_{1}\left(y\right),-{b}_{2}\left(y\right)\right)\end{array}$
 $\begin{array}{ccc}& =& M\left(-y,-q\right)\end{array}$
Here is another situation that occurs in certain applications.
Definition 5.3 Say that $\left({b}_{1},{b}_{2}\right)$  is a regular boundary operator if there is a projection $\Pi :X\to {X}_{0}:=Ker\left({b}_{1},{b}_{2}\right)$  so that the bounded linear map $\left(\Pi ,{b}_{1},{b}_{2}\right):X\to Ker\left({b}_{1},{b}_{2}\right)\oplus {H}_{1}\oplus {H}_{2}$  is an isomorphism.
Denote by $K:X\to {X}_{0}$  the projection in such a way that the bounded linear map $\left(K,{b}_{1},{b}_{2}\right):X\to {X}_{0}\oplus {H}_{1}\oplus {H}_{2}$  is an isomorphism. We can identify ${X}^{*}$  with the space ${X}_{0}^{*}\oplus {H}_{1}\oplus {H}_{2}$  in such a way that the duality between $X$  and ${X}^{*}$  is given by:
$〈x,p〉=〈x,\left({p}_{0},{p}_{1},{p}_{2}\right)〉=〈x,{K}^{*}{p}_{0}{〉}_{X,{X}^{*}}+〈{b}_{1}\left(x\right),{p}_{1}{〉}_{{H}_{1}}+〈{b}_{2}\left(x\right),{p}_{2}{〉}_{{H}_{2}}.$
Proposition 5.2 Let $\ell$  be a self dual boundary Lagrangian on the Hilbertian product ${H}_{1}×{H}_{2}$  , and let $\left(\Lambda ,{b}_{1},{b}_{2}\right):X\to {X}^{*}×{H}_{1}×{H}_{2}$  be a regular skew symmetric triplet where $X$  is a reflexive Banach space. Consider $L$  to be a Lagrangian on $X$  such that for each $x\in X$  , the map $p\to L\left(x,p\right)$  is continuous on ${X}^{*}$  .
• 1. If $L$  is a Lagrangian on $X×{X}^{*}$  that is anti-self dual on the graph of $\Lambda$  , then the Lagrangian $M\left(x,p\right)=L\left(x,\Lambda x+{K}^{*}{p}_{0}\right)+\ell \left({b}_{1}\left(x\right)+{p}_{1},{b}_{2}\left(x\right)-{p}_{2}\right)$  is partially anti-self dual. Here ${X}^{*}$  is identified with ${X}_{0}^{*}\oplus {H}_{1}\oplus {H}_{2}$  .
• 2. If $L$  is anti-self dual on $X$  , then $M$  is anti-self dual on $X×{X}_{0}^{*}$  .
Proof: Fix $\left(q,y\right)\in {X}^{*}×X$  , with $q=\left({q}_{0},0,0\right)$  and calculate
 $\begin{array}{ccc}{M}^{*}\left(q,y\right)& =& sup\left\{〈q,x{〉}_{X}+〈y,p{〉}_{X}-M\left(x,p\right);\left(x,p\right)\in X×{X}^{*}\right\}\end{array}$
 $\begin{array}{ccc}& =& sup\left\{〈x,{K}^{*}{q}_{0}〉+〈{b}_{1}\left(y\right),{p}_{1}〉+〈{b}_{2}\left(y\right),{p}_{2}〉+〈y,{K}^{*}{p}_{0}〉\end{array}$
 $\begin{array}{ccc}& & -L\left(x,\Lambda x+{K}^{*}{p}_{0}\right)-\ell \left({b}_{1}\left(x\right)+{p}_{1},{b}_{2}\left(x\right)-{p}_{2}\right);x\in X,{p}_{0}\in {X}_{0}^{*},{p}_{1}\in {H}_{1},{p}_{2}\in {H}_{2}\right\}\end{array}$
where $\left({p}_{0},{p}_{1},{p}_{2}\right)$  represent $p\in {X}^{*}$  .
Since the operator $A=\Lambda -\frac{1}{2}\left({b}_{2}^{*}{b}_{2}-{b}_{1}^{*}{b}_{1}\right)$  is skew-adjoint on $X$  , we can apply the results of the last section to get that $\left(A+\epsilon I\right)$  is onto for each $\epsilon >0$  . In other words, $A$  has dense range in ${X}^{*}$  , which yields that $Range\left(\Lambda \right)+{K}^{*}\left({X}_{0}^{*}\right)=Range\left(\Lambda \right)+Ker\left({b}_{1},{b}_{2}{\right)}^{\perp }\supset Range\left(\Lambda -\frac{1}{2}\left({b}_{2}^{*}{b}_{2}-{b}_{1}^{*}{b}_{1}\right)=Range\left(A\right)$  is also dense in ${X}^{*}$  .
Setting $r=\Lambda x+{K}^{*}{p}_{0}$  , ${f}_{1}={b}_{1}\left(x\right)+{p}_{1}$  and ${f}_{2}={b}_{2}\left(x\right)-{p}_{2}$  , we obtain that
 $\begin{array}{ccc}{M}^{*}\left(q,y\right)& =& sup\left\{〈x,{K}^{*}{q}_{0}〉+〈y,r-\Lambda x〉+〈{b}_{1}\left(y\right),{f}_{1}-{b}_{1}\left(x\right)〉+〈{b}_{2}\left(y\right),{b}_{2}\left(x\right)-{f}_{2}〉.\end{array}$
 $\begin{array}{ccc}& & -L\left(x,r\right)-\ell \left({f}_{1},{f}_{2}\right);x\in X,r\in Range\left(\Lambda \right)+{K}^{*}\left({X}_{0}^{*}\right),{f}_{1}\in {H}_{1},{f}_{2}\in {H}_{2}\right\}\end{array}$
 $\begin{array}{ccc}& =& sup\left\{〈x,{K}^{*}{q}_{0}〉+〈y,r-\Lambda x〉+〈{b}_{1}\left(y\right),{f}_{1}-{b}_{1}\left(x\right)〉+〈{b}_{2}\left(y\right),{b}_{2}\left(x\right)-{f}_{2}〉.\end{array}$
 $\begin{array}{ccc}& & -L\left(x,r\right)-\ell \left({f}_{1},{f}_{2}\right);x\in X,r\in {X}^{*},{f}_{1}\in {H}_{1},{f}_{2}\in {H}_{2}\right\}\end{array}$
 $\begin{array}{ccc}& =& sup\left\{〈x,{K}^{*}{q}_{0}〉+〈{b}_{1}\left(y\right),{f}_{1}〉-〈{b}_{2}\left(y\right),{f}_{2}〉+〈y,r〉+〈\Lambda y,x〉\end{array}$
 $\begin{array}{ccc}& & -L\left(x,r\right)-\ell \left({f}_{1},{f}_{2}\right);x\in X,r\in {X}^{*},{f}_{1}\in {H}_{1},{f}_{2}\in {H}_{2}\right\}\end{array}$
 $\begin{array}{ccc}& =& {L}^{*}\left({K}^{*}{q}_{0}+\Lambda y,y\right)+{\ell }^{*}\left({b}_{1}\left(y\right),-{b}_{2}\left(y\right)\right)\end{array}$
 $\begin{array}{ccc}& =& L\left(-y,-{K}^{*}{q}_{0}-\Lambda y\right)+\ell \left(-{b}_{1}\left(y\right),-{b}_{2}\left(y\right)\right)\end{array}$
 $\begin{array}{ccc}& =& M\left(-y,-q\right)\end{array}$
since $\left({q}_{0},0,0\right)$  represents $q$  in ${X}_{0}^{*}×{H}_{1}×{H}_{2}$  .
In the case where $\Lambda$  is essentially onto (modulo the boundary) we have yet another useful setting.
Definition 5.4 Say that a skew symmetric triplet $\left(\Lambda ,{b}_{1},{b}_{2}\right):X\to {X}^{*}×{H}_{1}×{H}_{2}$  is a nice boundary operator if the map $\left(\Lambda ,{b}_{1}\right):X\to Range\left(\Lambda \right)\oplus {H}_{1}$  is an isomorphism.
In this case, we identify ${X}^{*}$  with the space ${X}_{0}\oplus {H}_{1}$  where ${X}_{0}=X/Ker\left(\Lambda \right)\right)\sim Range\left(\Lambda \right)$  in such a way that the duality between $X$  and ${X}^{*}$  is given by:
$〈x,p〉=〈x,\left({p}_{0},{p}_{1}\right)〉=〈\Lambda x,{p}_{0}〉+〈{b}_{1}\left(x\right),{p}_{1}〉.$
Proposition 5.3 Let $\ell$  be a self dual boundary Lagrangian on the Hilbertian product ${H}_{1}×{H}_{2}$  , and let $\left(\Lambda ,{b}_{1},{b}_{2}\right):X\to {X}^{*}×{H}_{1}×{H}_{2}$  be a skew symmetric triplet on a reflexive Banach space $X$  that is a nice boundary operator. Then
• 1. If $L$  is a Lagrangian on $X×{X}^{*}$  that is anti-self dual on the graph of $\Lambda$  , then the Lagrangian $N\left(x,p\right)=L\left(x+{p}_{0},\Lambda x\right)+\ell \left({b}_{1}\left(x\right)+{p}_{1},{b}_{2}\left(x\right)\right)$  is partially anti-self dual on $X$  .
• 2. If $L$  is anti-self dual on $X$  , then $N$  is anti-self dual on the elements of $X×\left({X}_{0}\oplus \left\{0\right\}\right)$
Proof: Indeed, fix $\left(\left({q}_{0},0\right),y\right)\in \left({X}_{0}×{H}_{1}\right)×X$  and calculate
 $\begin{array}{ccc}{N}^{*}\left(q,y\right)& =& sup\left\{〈\Lambda x,{q}_{0}〉+〈\Lambda y,{p}_{0}〉+〈{b}_{1}\left(y\right),{p}_{1}〉\end{array}$
 $\begin{array}{ccc}& & -L\left(x+{p}_{0},\Lambda x\right)+\ell \left({b}_{1}\left(x\right)+{p}_{1},{b}_{2}\left(x\right)\right);x\in X,{p}_{0}\in {X}_{0},{p}_{1}\in {H}_{1}\right\}\end{array}$
Setting $r=x+{p}_{0}$  , ${f}_{1}={b}_{1}\left(x\right)+{p}_{1}$  , we obtain that
 $\begin{array}{ccc}{N}^{*}\left(q,y\right)& =& sup\left\{〈\Lambda x,{q}_{0}〉+〈\Lambda y,r-x〉+〈{b}_{1}\left(y\right),{f}_{1}-{b}_{1}\left(x\right)〉\end{array}$
 $\begin{array}{ccc}& & -L\left(r,\Lambda x\right)-\ell \left({f}_{1},{b}_{2}\left(x\right)\right);x\in X,r\in X,{f}_{1}\in {H}_{1}\right\}\end{array}$
 $\begin{array}{ccc}& =& sup\left\{〈\Lambda x,{q}_{0}〉+〈\Lambda y,r〉+〈y,\Lambda x〉+〈{b}_{1}\left(y\right),{f}_{1}〉-〈{b}_{2}\left(y\right),{b}_{2}\left(x\right)〉\end{array}$
 $\begin{array}{ccc}& & -L\left(r,\Lambda x\right)-\ell \left({f}_{1},{b}_{2}\left(x\right)\right);x\in X,r\in X,{f}_{1}\in {H}_{1}\right\}\end{array}$
Since $X$  can be identified with ${X}_{0}\oplus {H}_{2}$  via the correspondence $x\to \left(\Lambda x,{b}_{2}\left(x\right)\right)$  , we obtain:
 $\begin{array}{ccc}{N}^{*}\left(q,y\right)& =& sup\left\{〈s,{q}_{0}+y〉+〈\Lambda y,r〉+〈{b}_{1}\left(y\right),{f}_{1}〉-〈{b}_{2}\left(y\right),{f}_{2}〉\end{array}$
 $\begin{array}{ccc}& & -L\left(r,s\right)-\ell \left({f}_{1},{f}_{2}\right);s\in {X}_{0},r\in X,{f}_{1}\in {H}_{1},{f}_{2}\in {H}_{2}\right\}\end{array}$
 $\begin{array}{ccc}& =& sup\left\{〈s,{q}_{0}+y〉+〈\Lambda y,r〉-L\left(r,s\right);s\in {X}_{0},r\in X\right\}\end{array}$
 $\begin{array}{ccc}& & +sup\left\{〈{b}_{1}\left(y\right),{f}_{1}〉-〈{b}_{2}\left(y\right),{f}_{2}〉-\ell \left({f}_{1},{f}_{2}\right);{f}_{1}\in {H}_{1},{f}_{2}\in {H}_{2}\right\}\end{array}$
 $\begin{array}{ccc}& =& {L}^{*}\left(\Lambda y,{q}_{0}+y\right)+{\ell }^{*}\left({b}_{1}\left(y\right),-{b}_{2}\left(y\right)\right)\end{array}$
 $\begin{array}{ccc}& =& L\left(-{q}_{0}-y,-\Lambda y\right)+\ell \left(-{b}_{1}\left(y\right),-{b}_{2}\left(y\right)\right)\end{array}$
 $\begin{array}{ccc}& =& N\left(-y,-q\right).\end{array}$

5.2 Variational properties of ASD Lagrangians with boundary terms

One can now deduce the following
Theorem 5.5 Let $\left(\Lambda ,{b}_{1},{b}_{2}\right):X\to {X}^{*}×{H}_{1}×{H}_{2}$  be a skew symmetric triplet, $\ell$  a self dual boundary Lagrangian on ${H}_{1}×{H}_{2}$  and let $L:X×{X}^{*}\to \mathbf{R}\cup \left\{+\infty \right\}$  be anti-self dual on the graph of $\Lambda$  . Assume one of the following hypothesis:
• 1. (A) The boundary operator $\left({b}_{1},{b}_{2}\right)$  is regular and for every $x\in X$  , the map $p\to L\left(x,p\right)$  is bounded on the bounded sets of ${X}^{*}$  .
• 2. (B) The triplet $\left(\Lambda ,{b}_{1},{b}_{2}\right)$  is a nice boundary operator and the map $x\to L\left(x,0\right)$  is bounded on the bounded sets of $X$  .
Then, there exists $\overline{x}\in X$  such that:
 $\begin{array}{c}L\left(\overline{x},\Lambda \overline{x}\right)+\ell \left({b}_{1}\overline{x},{b}_{2}\overline{x}\right)={inf}_{x\in X}\left\{L\left(x,\Lambda x\right)+\ell \left({b}_{1}x,{b}_{2}x\right)\right\}=0.\end{array}$ (43)
Moreover, we have
 $\begin{array}{c}\left\{\begin{array}{ccc}L\left(\overline{x},\Lambda \overline{x}\right)+〈\overline{x},\Lambda \overline{x}〉& =& 0.\\ \left(-\Lambda \overline{x},-\overline{x}\right)& \in & \partial L\left(\overline{x},\Lambda \overline{x}\right)\\ \ell \left({b}_{1}\left(\overline{x}\right),{b}_{2}\left(\overline{x}\right)\right)& =& \frac{1}{2}\left(\parallel {b}_{2}\overline{x}{\parallel }^{2}-\parallel {b}_{1}\overline{x}{\parallel }^{2}\right).\end{array}\end{array}$ (44)
In particular, for any $a\in {H}_{1}$  there exists $\overline{x}\in X$  such that ${b}_{1}\left(\overline{x}\right)=a$  and satisfying ( 44 ). It is obtained as a minimizer on $X$  of the functional $I\left(x\right)=L\left(x,\Lambda x\right)+\frac{1}{2}\parallel {b}_{1}\left(x\right){\parallel }^{2}-2〈a,{b}_{1}\left(x\right)〉+\parallel a{\parallel }^{2}+\frac{1}{2}\parallel {b}_{2}\left(x\right){\parallel }^{2}.$
Proof: Under case (A), we use proposition 3.2 to get that the Lagrangian $M\left(x,p\right)=L\left(x,\Lambda x+{p}_{0}\right)+\ell \left({b}_{1}\left(x\right)+{p}_{1},{b}_{2}\left(x\right)-{p}_{2}\right)$  is partially anti-self dual. In case (B), we use Proposition 5.3 to conclude that the Lagrangian $N\left(x,p\right)=L\left(x+{p}_{0},\Lambda x\right)+\ell \left({b}_{1}\left(x\right)+{p}_{1},{b}_{2}\left(x\right)\right)$  is partially anti-self dual on $X$  .
In both cases, the hypothesis implies that $M\left(0,p\right)$  (resp., $N\left(0,p\right)$  ) is bounded above on the bounded sets of ${X}^{*}$  . Theorem 4.1 then applies to yield $\overline{x}\in X$  such that ( 43 ) is satisfied.
To establish ( 44 ), write
 $\begin{array}{ccc}L\left(x,\Lambda x\right)+\ell \left({b}_{1}x,{b}_{2}x\right)& =& L\left(x,\Lambda x\right)+〈x,\Lambda x〉-〈x,\Lambda x〉+\ell \left({b}_{1}x,{b}_{2}x\right)\end{array}$
 $\begin{array}{ccc}& =& L\left(x,\Lambda x\right)+〈x,\Lambda x〉-\frac{1}{2}\left(\parallel {b}_{2}x{\parallel }^{2}-\parallel {b}_{1}x{\parallel }^{2}\right)+\ell \left({b}_{1}x,{b}_{2}x\right).\end{array}$
Since $L\left(x,p\right)\ge -〈x,p〉$  and $\ell \left(r,s\right)\ge \frac{1}{2}\left(\parallel s{\parallel }^{2}-\parallel r{\parallel }^{2}\right)$  , we immediately obtain ( 44 ).
In particular, for any $a\in {H}_{1}$  , consider the boundary Lagrangian, $\ell \left(r,s\right)=\frac{1}{2}\parallel r{\parallel }^{2}-2〈a,r〉+\parallel a{\parallel }^{2}+\frac{1}{2}\parallel s{\parallel }^{2}.$  which is clearly self-dual. We then get
 $\begin{array}{ccc}L\left(x,\Lambda x\right)+\ell \left({b}_{1}x,{b}_{2}x\right)& =& L\left(x,\Lambda x\right)+〈x,\Lambda x〉-\frac{1}{2}\left(\parallel {b}_{2}x{\parallel }^{2}-\parallel {b}_{1}x{\parallel }^{2}\right)+\ell \left({b}_{1}x,{b}_{2}x\right)\end{array}$
 $\begin{array}{ccc}& =& L\left(x,\Lambda x\right)+〈x,\Lambda x〉+\parallel {b}_{1}\left(x\right)-a{\parallel }^{2}.\end{array}$
In other words, $\overline{x}\in X$  is a solution of ${inf}_{x\in X}\left\{L\left(x,\Lambda x\right)+〈x,\Lambda x〉+\parallel {b}_{1}\left(x\right)-a{\parallel }^{2}\right\}=0$  , and since $L\left(x,p\right)\ge -〈x,p〉$  , we obtain:
 $\begin{array}{c}\left\{\begin{array}{ccc}L\left(\overline{x},\Lambda \overline{x}\right)+〈\overline{x},\Lambda \overline{x}〉& =& 0.\\ {b}_{1}\left(\overline{x}\right)& =& a.\end{array}\end{array}$ (45)

5.3 Variational principle for operators which are positive modulo a boundary

Consider again $\left({b}_{1},{b}_{2}\right):X\to {H}_{1}×{H}_{2}$  to be a regular boundary operator.
Definition 5.6 Say that $A:X\to {X}^{*}$  is positive modulo the boundary operator $\left({b}_{1},{b}_{2}\right)$  if the operator $A-\frac{1}{2}\left({b}_{2}^{*}{b}_{2}-{b}_{1}^{*}{b}_{1}\right)$  is positive.
Corollary 5.7 Let $A:X\to {X}^{*}$  be positive modulo a boundary operator $\left({b}_{1},{b}_{2}\right)$  and set $\Lambda ={A}^{a}+\frac{1}{2}\left({b}_{2}^{*}{b}_{2}-{b}_{1}^{*}{b}_{1}\right)$  . Let $\phi$  be a convex function in $\mathcal{C}\left(X\right)$  with $0$  in its domain and such that one of the following conditions holds:
• 1. (A) The boundary operator $\left({b}_{1},{b}_{2}\right)$  is regular and ${lim}_{\parallel x\parallel \to \infty }\parallel x{\parallel }^{-1}\left\{\phi \left(x\right)+\frac{1}{2}〈Ax,x〉-\frac{1}{4}\left(\parallel {b}_{2}x{\parallel }^{2}-\parallel {b}_{1}x{\parallel }^{2}\right)\right\}=+\infty .$
• 2. (B) The triplet $\left(\Lambda ,{b}_{1},{b}_{2}\right)$  is a nice boundary operator and $\phi$  is bounded on the bounded sets of $X$  .
Then for any $a\in {H}_{1}$  and any $f\in {X}^{*}$  , the equation
 $\begin{array}{c}\left\{\begin{array}{ccc}-Ax& \in & \partial \phi \left(x\right)+f\\ {b}_{1}\left(x\right)& =& a\\ \end{array}\end{array}$ (46)
has a solution $\overline{x}\in X$  that is a minimizer of the problem:
$I\left(x\right)=\psi \left(x\right)+{\psi }^{*}\left(-{A}^{a}x-\frac{1}{2}{b}_{2}^{*}{b}_{2}x+\frac{1}{2}{b}_{1}^{*}{b}_{1}x\right)+\frac{1}{2}\left(\parallel {b}_{1}\left(x\right){\parallel }^{2}+\parallel {b}_{2}\left(x\right){\parallel }^{2}\right)-2〈a,{b}_{1}\left(x\right)〉+\parallel a{\parallel }^{2}$  where $\psi \left(x\right)=\phi \left(x\right)+\frac{1}{2}〈Ax,x〉-\frac{1}{4}\left(\parallel {b}_{2}x{\parallel }^{2}-\parallel {b}_{1}x{\parallel }^{2}\right)+〈f,x〉$  .
Proof: Let $B=A-\frac{1}{2}\left({b}_{2}^{*}{b}_{2}-{b}_{1}^{*}{b}_{1}\right)$  and decompose it into its symmetric ${B}^{s}={A}^{s}-\frac{1}{2}\left({b}_{2}^{*}{b}_{2}-{b}_{1}^{*}{b}_{1}\right)$  and its anti-symmetric part ${B}^{a}={A}^{a}$  , by simply writing ${B}^{s}=\frac{1}{2}\left(B+{B}^{*}\right)$  and ${B}^{a}=\frac{1}{2}\left(B-{B}^{*}\right)$  . We can then write $A={B}^{s}+\Lambda$  where $\Lambda :={A}^{a}+\frac{1}{2}\left({b}_{2}^{*}{b}_{2}-{b}_{1}^{*}{b}_{1}\right)={B}^{a}+\frac{1}{2}\left({b}_{2}^{*}{b}_{2}-{b}_{1}^{*}{b}_{1}\right)$  which means that $\Lambda$  is skew-symmetric modulo the boundary operator $\left({b}_{1},{b}_{2}\right)$  . For $f\in {X}^{*}$  , consider the convex functional $\psi \left(x\right)=\phi \left(x\right)+\frac{1}{2}〈{B}^{s}x,x〉+〈f,x〉=\phi \left(x\right)+\frac{1}{2}〈Ax,x〉-\frac{1}{4}\left(\parallel {b}_{2}x{\parallel }^{2}-\parallel {b}_{1}x{\parallel }^{2}\right)+〈f,x〉.$  The proposed minimization problems amounts to applying Theorem 5.5 to the anti-self dual Lagrangian $L\left(x,p\right)=\psi \left(x\right)+{\psi }^{*}\left(-p\right)$  , the operator $\Lambda$  and the boundary Lagrangian $\ell \left(r,s\right)=\frac{1}{2}\parallel r{\parallel }^{2}-2〈a,r〉+\parallel a{\parallel }^{2}+\frac{1}{2}\parallel s{\parallel }^{2}$  . Note that $I\left(x\right)=\psi \left(x\right)+{\psi }^{*}\left(-\Lambda x\right)+〈x,\Lambda x〉+\parallel {b}_{1}\left(x\right)-a{\parallel }^{2}.$  The fact that the minimum is attained at some $\overline{x}$  and is equal to $0$  , implies that ${b}_{1}\left(\overline{x}\right)=a$  and that $\psi \left(\overline{x}\right)+{\psi }^{*}\left(-\Lambda \overline{x}\right)=-〈\Lambda \overline{x},\overline{x}〉$  which means that $-{A}^{a}\left(\overline{x}\right)-\frac{1}{2}\left({b}_{2}^{*}{b}_{2}-{b}_{1}^{*}{b}_{1}\right)\left(\overline{x}\right)=-\Lambda \overline{x}\in \partial \psi \left(\overline{x}\right)=\partial \phi \left(\overline{x}\right)+{A}^{s}\left(\overline{x}\right)-\frac{1}{2}\left({b}_{2}^{*}{b}_{2}-{b}_{1}^{*}{b}_{1}\right)\left(\overline{x}\right)+f$  and therefore $-A\overline{x}\in \partial \phi \left(\overline{x}\right)+f$  .
Remark 5.8 Again the above applies to functions $\phi$  that are A-convex modulo a boundary $\left({b}_{1},{b}_{2}\right)$  meaning those functions $\phi$  such that there exists an operator $A$  such that $\psi \left(x\right)=\phi \left(x\right)+\frac{1}{2}〈Ax,x〉-\frac{1}{4}\left(\parallel {b}_{2}x{\parallel }^{2}-\parallel {b}_{1}x{\parallel }^{2}\right)$  is convex and lower semi-continuous.
Example 6: A variational principle for non-linear transport equations As in example 3, Let $\mathbf{a}:\Omega \to {\mathbf{R}}^{\mathbf{n}}$  and ${a}_{0}:\Omega \to \mathbf{R}$  be two smooth functions on a bounded domain $\Omega$  of ${\mathbf{R}}^{\mathbf{n}}$  , and consider the first order linear operator $Av=\mathbf{a}\cdot \nabla v={\Sigma }_{i=1}^{n}{a}_{i}\frac{\partial v}{\partial {x}_{i}}and\Lambda v=\mathbf{a}\cdot \nabla v+{a}_{0}v.$  As in [5, we shall assume throughout that the vector field ${\Sigma }_{i=1}^{n}{a}_{i}\frac{\partial v}{\partial {x}_{i}}$  is actually the restriction of a smooth vector field ${\Sigma }_{i=1}^{n}{\overline{a}}_{i}\frac{\partial v}{\partial {x}_{i}}$  defined on an open neighborhood $X$  of $\overline{\Omega }$  and that each ${\overline{a}}_{i}$  is a ${C}^{1,1}$  function on $X$  . We also assume that the boundary of $\Omega$  is piecewise ${C}^{1}$  , in such a way that the outer normal $\mathbf{n}$  is defined almost everywhere on $\partial \Omega$  . In this case, if we denote by ${\Sigma }_{-}=\left\{x\in \partial \Omega ;\mathbf{n}\left(x\right)\cdot \mathbf{a}\left(x\right)<0\right\}and{\Sigma }_{+}=\partial \Omega \{\Sigma }_{-}=\left\{x\in \partial \Omega ;\mathbf{n}\left(x\right)\cdot \mathbf{a}\left(x\right)\ge 0\right\},$  then a trace ${u}_{{|}_{{\Sigma }_{-}}}$  makes sense in ${L}_{loc}^{2}\left({\Sigma }_{-}\right)$  as soon as $u\in {L}^{2}\left(\Omega \right)$  and $\Lambda u\in {L}^{2}\left(\Omega \right)$  .
Let now $\beta :\mathbf{R}\to \mathbf{R}$  be a continuous nondecreasing function so that its antiderivative $j$  is convex, and let $f\in {L}^{2}\left(\Omega \right)$  . We are interested in finding variationally solutions for the nonlinear transport equation:
 $\begin{array}{c}\left\{\begin{array}{ccc}-\mathbf{a}\cdot \nabla u-{a}_{0}u& =& \beta \left(u\right)+fon\Omega \\ u\left(x\right)& =& {u}_{0}on{\Sigma }_{-}.\\ \end{array}\end{array}$ (47)
First, we identify the appropriate underlying space. Consider the space ${H}^{1}\left(\Omega \right)=\left\{u\in {L}^{2}\left(\Omega \right);Au\in {L}^{2}\left(\Omega \right)\right\}.$  equipped with the norm $\parallel u{\parallel }_{{H}^{1}}=\parallel u{\parallel }_{2}+\parallel Au{\parallel }_{2}$  . As noticed in [5, that a function $u$  belongs to ${H}^{1}\left(\Omega \right)$  does not necessarily guarantee that its trace ${u}_{{|}_{{\Sigma }_{-}}}$  is in the space ${L}_{A}^{2}\left({\Sigma }_{-}\right)=\left\{u\in {L}_{loc}^{2}\left({\Sigma }_{-}\right);{\int }_{{\Sigma }_{-}}|u\left(x\right){|}^{2}|\mathbf{n}\left(x\right)\cdot \mathbf{a}\left(x\right)|d\sigma <+\infty \right\}.$  However, if $u\in {H}^{1}\left(\Omega \right)$  and ${u}_{{\Sigma }_{-}}\in {L}_{A}^{2}\left({\Sigma }_{-}\right)$  , then necessarily ${u}_{{\Sigma }_{+}}\in {L}_{A}^{2}\left({\Sigma }_{+}\right)$  . The appropriate space for our setting is therefore ${H}_{A}^{1}\left(\Omega \right)=\left\{u\in {H}^{1}\left(\Omega \right);{u}_{{|}_{{\Sigma }_{-}}}\in {L}_{A}^{2}\left({\Sigma }_{-}\right)\right\}.$  equipped with the norm $\parallel u{\parallel }_{{H}_{A}^{1}}=\parallel u{\parallel }_{2}+\parallel Au{\parallel }_{2}+\parallel {u}_{{|}_{{\Sigma }_{-}}}{\parallel }_{{}_{{L}_{A}^{2}\left({\Sigma }_{-}\right)}}$  .
To define appropriate boundary spaces, we follow [11and consider for each open subset $\Gamma$  of $\partial \Omega$  , the space ${H}_{{}_{00}}^{{}^{1/2}}\left(\Gamma \right)=\left\{v\in {L}_{A}^{2}\left(\Gamma \right);\exists w\in {H}^{1}\left(\Omega \right),w=0on\partial \Omega \\Gamma ,andw=von\Gamma \right\}$  A trace theorem ([23, Vol III. p. 307) or [3) yields that the restriction mapping $u\to {u}_{\Gamma }$  is a continuous surjective map from $V=\left\{v\in {H}^{1}\left(\Omega \right);{v}_{{|}_{\partial \Omega \\Gamma }}=0\right\}$  onto ${H}_{{}_{00}}^{{}^{1/2}}\left(\Gamma \right)$  . It follows that there is a continuous surjection from ${H}_{A}^{1}\left(\Omega \right)$  onto ${H}_{0}^{1}\left(\Omega \right)\oplus {H}_{{}_{00}}^{{}^{1/2}}\left({\Sigma }_{-}\right)\oplus {H}_{{}_{00}}^{{}^{1/2}}\left({\Sigma }_{+}\right)$  via the map $T:{H}_{A}^{1}\left(\Omega \right)\to {H}_{0}^{1}\left(\Omega \right)\oplus {H}_{{}_{00}}^{{}^{1/2}}\left({\Sigma }_{-}\right)\oplus {H}_{{}_{00}}^{{}^{1/2}}\left({\Sigma }_{+}\right),$  given by $Tu=\left(Ku,{u}_{|{\Sigma }_{-}},{u}_{{|}_{{\Sigma }_{+}}}\right)$  , where $K:{H}_{A}^{1}\left(\Omega \right)\to {H}_{0}^{1}\left(\Omega \right)$  is the operator that associates to $u\in {H}_{A}^{1}\left(\Omega \right)$  the unique function $w\in {H}_{0}^{1}\left(\Omega \right)$  such that $\Delta w=\Delta u$  and $w=0$  on $\partial \Omega$  .
If now ${a}_{0}\left(x\right)-\frac{1}{2}diva\left(x\right)\ge 0$  on $\Omega$  , then $\Lambda$  is positive modulo the boundary operators $u\to \left({u}_{{|}_{{\Sigma }_{-}}},u{|}_{{\Sigma }_{+}}\right)\in {L}^{2}\left({\Sigma }_{-}\right)×{L}^{2}\left({\Sigma }_{+}\right)$  , since ${\int }_{\Omega }u\Lambda udx={\int }_{\Omega }\left(\left(\mathbf{a}\cdot \nabla u\right)u+{a}_{0}|u{|}^{2}\right)dx={\int }_{\Omega }\left({a}_{0}-\frac{1}{2}div\mathbf{a}\right)|u{|}^{2}dx+{\int }_{\partial \Omega }|u{|}^{2}\mathbf{n}\cdot \mathbf{a}d\sigma .$  and the operator ${\Lambda }_{1}\left(u\right):=\mathbf{a}\cdot \nabla u+\frac{1}{2}div\left(\mathbf{a}\right)u=\Lambda \left(u\right)-\left({a}_{0}-\frac{1}{2}div\mathbf{a}\right)u$  is therefore skew-adjoint modulo that boundary since then
 $\begin{array}{c}{\int }_{\Omega }v{\Lambda }_{1}udx=-{\int }_{\Omega }u{\Lambda }_{1}vdx+{\int }_{\partial \Omega }uv\mathbf{n}\cdot \mathbf{a}d\sigma .\end{array}$ (48)
We can now state:
Theorem 5.9 Assume the coercivity condition ${a}_{0}\left(x\right)-\frac{1}{2}diva\left(x\right)\ge \alpha >0$  on $\Omega$  .
For any $f\in {L}^{2}\left(\Omega \right)$  and ${u}_{0}\in {L}_{A}^{2}\left({\Sigma }_{-}\right)$  , consider the following functional on the space ${H}_{A}^{1}\left(\Omega \right)$
 $\begin{array}{ccc}I\left(u\right)& =& \psi \left(u\right)+{\psi }^{*}\left(-\mathbf{a}\cdot \nabla u-\frac{1}{2}\left(div\mathbf{a}\right)u\right)+\frac{1}{2}{\int }_{{\Sigma }_{+}}|u\left(x\right){|}^{2}\mathbf{n}\left(x\right)\cdot \mathbf{a}\left(x\right)d\sigma -\frac{1}{2}{\int }_{{\Sigma }_{-}}|u\left(x\right){|}^{2}\mathbf{n}\left(x\right)\cdot \mathbf{a}\left(x\right)d\sigma \end{array}$
 $\begin{array}{ccc}& & +2{\int }_{{\Sigma }_{-}}u\left(x\right){u}_{0}\left(x\right)\mathbf{n}\left(x\right)\cdot \mathbf{a}\left(x\right)d\sigma -{\int }_{{\Sigma }_{-}}|{u}_{0}\left(x\right){|}^{2}\mathbf{n}\left(x\right)\cdot \mathbf{a}\left(x\right)d\sigma .\end{array}$
and where $\psi$  is the convex functional on ${L}^{2}\left(\Omega \right)$  defined by $\psi \left(u\right)={\int }_{\Omega }\left\{j\left(u\left(x\right)\right)+f\left(x\right)u\left(x\right)+\frac{1}{2}\left({a}_{0}-\frac{1}{2}diva\right){u}^{2}\right)\right\}dx$  and where ${\psi }^{*}$  is its Legendre conjugate.
Then there exists a solution $\overline{u}$  for ( 47 ) that is obtained as a minimizer of the problem:
$I\left(\overline{u}\right)=inf\left\{I\left(u\right);u\in {H}_{A}^{1}\left(\Omega \right)\right\}=0.$
Proof: The only problem remaining is the fact that the convex functional $\psi$  defined by:
$\psi \left(u\right)={\int }_{\Omega }\left\{j\left(u\left(x\right)\right)+f\left(x\right)u\left(x\right)+\frac{1}{2}\left({a}_{0}-\frac{1}{2}diva\right){u}^{2}\right)\right\}dx$  is not necessarily coercive on ${H}_{A}^{1}\left(\Omega \right)$  , so we consider instead for each $\epsilon >0$  , the functional ${\phi }_{\epsilon }\left(u\right)=\psi \left(u\right)+\frac{\epsilon }{2}{\int }_{\Omega }|\nabla u{|}^{2}dx$  which obviously is. Assuming without loss that ${u}_{0}=0$  and setting
 $\begin{array}{c}{I}_{\epsilon }\left(u\right)={\phi }_{\epsilon }\left(u\right)+{\phi }_{\epsilon }^{*}\left(-{\Lambda }_{1}u\right)+\frac{1}{2}{\int }_{{\Sigma }_{+}}|u\left(x\right){|}^{2}\mathbf{n}\cdot \mathbf{a}d\sigma -\frac{1}{2}{\int }_{{\Sigma }_{-}}|u\left(x\right){|}^{2}\mathbf{n}\cdot \mathbf{a}d\sigma .\end{array}$ (49)
The above lemma now applies and we get ${u}_{\epsilon }\in {H}_{A}^{1}\left(\Omega \right)$  such that ${inf}_{u\in {H}_{A}^{1}\left(\Omega \right)}{I}_{\epsilon }\left(u\right)={I}_{\epsilon }\left({u}_{\epsilon }\right)=0,$  This means that ${u}_{\epsilon }$  belongs to $Dom\left(\partial {\phi }_{\epsilon }\right)$  and satisfies $-{\Lambda }_{1}{u}_{\epsilon }\in \partial {\phi }_{\epsilon }\left({u}_{\epsilon }\right)$  , which implies $-{\Lambda }_{1}{u}_{\epsilon }=\beta \left({u}_{\epsilon }\right)+f+\left({a}_{0}-\frac{1}{2}diva\right){u}_{\epsilon }-\epsilon \Delta {u}_{\epsilon }.$  In other words, we have for each $\epsilon >0$  ,
 $\begin{array}{c}\left\{\begin{array}{ccc}\epsilon \Delta {u}_{\epsilon }-{a}_{0}{u}_{\epsilon }-{\Sigma }_{i=1}^{n}{a}_{i}\frac{\partial {u}_{\epsilon }}{\partial {x}_{i}}& =& \beta \left({u}_{\epsilon }\right)+fon\Omega \\ {u}_{\epsilon }& =& 0on{\Sigma }_{-},\\ \frac{\partial {u}_{\epsilon }}{\partial n}& =& 0on\partial \Omega \{\Sigma }_{-}.\end{array}\end{array}$ (50)
It is now standard to show that, as $\epsilon \to 0$  , ${u}_{\epsilon }$  converges in ${L}^{2}\left(\Omega \right)$  to a solution $u$  of ( 47 ). For details, see Bardos [5.

5.4 ASD Lagrangians on intermediate Hilbert spaces

As one can see in the previous example, it is more desirable to have coercivity on the space ${L}^{2}\left(\Omega \right)$  and therefore we need to “extend” anti-selfduality from the Banach space ${H}_{A}^{1}\left(\Omega \right)$  to the ambient Hilbert space ${L}^{2}\left(\Omega \right)$  .
This situation is common in applications to partial differential equations, where an ambient Hilbert space $H$  is usually present in such a way that $X$  is a dense subset of $H$  , and the identity injection $i:X\to H$  is continuous. The scalar product and the norm of $H$  are denoted by $\left(u,v\right)$  and $||$  respectively. By duality, the adjoint ${i}^{*}:{H}^{*}\to {X}^{*}$  is also one-to-one with dense range. One often identifies $H$  with its dual ${H}^{*}$  , in such a way that we have a representation of the form:
$X\subset H\equiv {H}^{*}\subset {X}^{*}$  . In this representation, we have $〈h,x〉=\left(h,x\right)$  whenever $h\in H$  and $x\in X$  .
The pair $\left(X,H\right)$  is sometimes called an evolution pair. We shall need the following notion.
Definition 5.10 Let $\left(\Lambda ,{b}_{1},{b}_{2}\right):X\to {X}^{*}×{H}_{1}×{H}_{2}$  be a skew symmetric triplet on a reflexive Banach space, and let $H$  be a Hilbert space so that $\left(X,H\right)$  is an evolution pair. We say that $\left(X,H,\Lambda \right)$  is a maximal evolution triple if ${X}_{0}=Ker\left({b}_{1},{b}_{2}\right)$  is dense in $H$  , $\Lambda$  maps $X$  into $H$  and if $X=\left\{x\in H;sup\left\{〈x,\Lambda y{〉}_{H}-\frac{1}{2}\left(\parallel {b}_{1}\left(x\right){\parallel }_{{H}_{1}}^{2}+\parallel {b}_{2}\left(x\right){\parallel }_{{H}_{2}}^{2}\right);y\in X,\parallel y{\parallel }_{H}\le 1\right\}<+\infty \right\}.$
Proposition 5.4 Let $\ell$  be a self dual boundary Lagrangian on the Hilbertian product ${H}_{1}×{H}_{2}$  , let $\left(\Lambda ,{b}_{1},{b}_{2}\right):X\to {X}^{*}×{H}_{1}×{H}_{2}$  be a regular skew symmetric triplet on a reflexive Banach space $X$  , and let $H$  be a Hilbert space such that $\left(X,H,\Lambda \right)$  is a maximal evolution triple.
If $L$  is anti-self dual on $H$  such that for each $p\in H$  , the map $x\to L\left(x,p\right)$  is continuous on $H$  , then the Lagrangian $M\left(x,p\right)=\left\{\begin{array}{ccc}L\left(x,\Lambda x+p\right)+\ell \left({b}_{1}\left(x\right),{b}_{2}\left(x\right)\right)& if& x\in X\\ +\infty & & otherwise\\ \end{array}$  is also anti-self dual on $H×H$  .
Proof: Fix $\left(q,y\right)\in H×X$  and calculate
 $\begin{array}{ccc}{M}^{*}\left(q,y\right)& =& sup\left\{〈q,x{〉}_{H}+〈y,p{〉}_{H}-M\left(x,p\right);\left(x,p\right)\in H×H\right\}\end{array}$
 $\begin{array}{ccc}& =& sup\left\{〈q,x{〉}_{X}+〈y,p{〉}_{X}-L\left(x,\Lambda x+p\right)-\ell \left({b}_{1}\left(x\right),{b}_{2}\left(x\right)\right);\left(x,p\right)\in X×H\right\}\end{array}$
Setting $r=\Lambda x+p$  , we obtain since $y\in X$  , that
 $\begin{array}{ccc}{M}^{*}\left(q,y\right)& =& sup\left\{〈x,q〉+〈y,r-\Lambda x〉-L\left(x,r\right)-\ell \left({b}_{1}\left(x\right),{b}_{2}\left(x\right)\right);\left(x,r\right)\in X×H\right\}\end{array}$
 $\begin{array}{ccc}& =& sup\left\{〈x,q〉+〈{b}_{1}\left(y\right),{b}_{1}\left(x\right)〉-〈{b}_{2}\left(y\right),{b}_{2}\left(x\right)〉+〈\Lambda y,x〉+〈y,r〉\end{array}$
 $\begin{array}{ccc}& & -L\left(x,r\right)-\ell \left({b}_{1}\left(x\right),{b}_{2}\left(x\right)\right);\left(x,r\right)\in X×H\right\}\end{array}$
 $\begin{array}{ccc}& =& sup\left\{〈x,q+\Lambda y〉+〈y,r〉-L\left(x,r\right)\end{array}$
 $\begin{array}{ccc}& & +〈{b}_{1}\left(y\right),{b}_{1}\left(x\right)〉-〈{b}_{2}\left(y\right),{b}_{2}\left(x\right)〉-\ell \left({b}_{1}\left(x\right),{b}_{2}\left(x\right)\right);\left(x,r\right)\in X×H\right\}\end{array}$
 $\begin{array}{ccc}& =& sup\left\{〈x,q+\Lambda y〉+〈y,r〉-L\left(x,r\right)+〈{b}_{1}\left(y\right),{b}_{1}\left(x+{x}_{0}\right)〉-〈{b}_{2}\left(y\right),{b}_{2}\left(x+{x}_{0}\right)〉\end{array}$
 $\begin{array}{ccc}& & -\ell \left({b}_{1}\left(x+{x}_{0}\right),{b}_{2}\left(x\right)+{x}_{0}\right);\left({x}_{0},x,r\right)\in {X}_{0}×X×H\right\}\end{array}$
 $\begin{array}{ccc}& =& sup\left\{〈w-{x}_{0},q+\Lambda y〉+〈y,r〉-L\left(w-{x}_{0},r\right)+〈{b}_{1}\left(y\right),{b}_{1}\left(w\right)〉-〈{b}_{2}\left(y\right),{b}_{2}\left(w\right)〉\end{array}$
 $\begin{array}{ccc}& & -\ell \left({b}_{1}\left(w\right),{b}_{2}\left(w\right);\left({x}_{0},w,r\right)\in {X}_{0}×X×H\right\}\end{array}$
Now use the fact that ${X}_{0}$  is dense in $H$  , and the continuity of $x\to L\left(x,p\right)$  on $H$  , to obtain
 $\begin{array}{ccc}{M}^{*}\left(q,y\right)& =& sup\left\{〈z,q+\Lambda y〉+〈y,r〉-L\left(z,r\right)+〈{b}_{1}\left(y\right),{b}_{1}\left(w\right)〉-〈{b}_{2}\left(y\right),{b}_{2}\left(w\right)〉\end{array}$
 $\begin{array}{ccc}& & -\ell \left({b}_{1}\left(w\right),{b}_{2}\left(w\right);\left(z,w,r\right)\in H×X×H\right\}\end{array}$
 $\begin{array}{ccc}& =& sup\left\{〈z,q+\Lambda y〉+〈y,r〉-L\left(z,r\right);\left(z,r\right)\in H×H\right\}\end{array}$
 $\begin{array}{ccc}& & +sup\left\{〈{b}_{1}\left(y\right),{b}_{1}\left(w\right)〉-〈{b}_{2}\left(y\right),{b}_{2}\left(w\right)〉-\ell \left({b}_{1}\left(w\right),{b}_{2}\left(w\right)\right);w\in X\right\}\end{array}$
 $\begin{array}{ccc}& =& sup\left\{〈z,q+\Lambda y〉+〈y,r〉-L\left(z,r\right);\left(z,r\right)\in H×H\right\}\end{array}$
 $\begin{array}{ccc}& & +sup\left\{〈{b}_{1}\left(y\right),a〉-〈{b}_{2}\left(y\right),b〉-\ell \left(a,b\right);\left(a,b\right)\in {H}_{1}×{H}_{2}\right\}\end{array}$
 $\begin{array}{ccc}& =& {L}^{*}\left(q+\Lambda y,y\right)+{\ell }^{*}\left({b}_{1}\left(y\right),-{b}_{2}\left(y\right)\right)\end{array}$
 $\begin{array}{ccc}& =& L\left(-y,-q-\Lambda y\right)+\ell \left(-{b}_{1}\left(y\right),-{b}_{2}\left(y\right)\right)\end{array}$
 $\begin{array}{ccc}& =& M\left(-y,-q\right)\end{array}$
If now $\left(q,y\right)\in H×\left(H\X\right)$  , then
 $\begin{array}{ccc}{M}^{*}\left(q,y\right)& \ge & sup\left\{〈x,q〉+〈y,r-\Lambda x〉-L\left(x,r\right)-\ell \left({b}_{1}\left(x\right),{b}_{2}\left(x\right)\right);\left(x,r\right)\in X×H\right\}\end{array}$
 $\begin{array}{ccc}& \ge & sup\left\{-\parallel x{\parallel }_{H}\parallel q{\parallel }_{H}+〈y,\Lambda x〉-L\left(x,0\right)-\ell \left({b}_{1}\left(x\right),{b}_{2}\left(x\right)\right);x\in X\right\}\end{array}$
 $\begin{array}{ccc}& \ge & sup\left\{-\parallel q{\parallel }_{H}+〈y,\Lambda x〉-C-\frac{1}{2}\left(\parallel {b}_{1}\left(x\right){\parallel }_{{H}_{1}}^{2}+\parallel {b}_{2}\left(x\right){\parallel }_{{H}_{2}}^{2}\right);x\in X,\parallel x{\parallel }_{H}\le 1\right\}\end{array}$
 $\begin{array}{ccc}& =& +\infty ,\end{array}$
since otherwise $y\in X$  which is a contradiction.
Remark 5.11 In practice, the Hilbert space is usually given and $X$  is usually obtained from the domain of some unbounded operator on $H$  . The full scope of this setting is developed in [20. For now, we give the following illustrative example
Example 7: More general transport equations Consider the following general transport equation
 $\begin{array}{c}\left\{\begin{array}{ccc}-\Lambda u& =& |u{|}^{p-2}u+Bu+fon\Omega \\ u\left(x\right)& =& {u}_{0}\left(x\right)on{\Sigma }_{-}.\\ \end{array}\end{array}$ (51)
where $B:{L}^{2}\left(\Omega \right)\to {L}^{2}\left(\Omega \right)$  is a positive bounded linear operator, $f\in {L}^{2}\left(\Omega \right)$  and ${u}_{0}\in {L}_{A}^{2}\left({\Sigma }_{-}\right)$  .
We again decompose $B$  into a symmetric and an anti-symmetric part, ${B}_{s}$  and ${B}_{a}$  , by writing ${B}_{s}=\frac{1}{2}\left(B+{B}^{*}\right)$  and ${B}_{a}=\frac{1}{2}\left(B-{B}^{*}\right)$  , and we consider the convex functional defined on ${L}^{2}$  by:
$\psi \left(u\right)=\frac{1}{2}{\int }_{\Omega }\left(\frac{1}{p}|u\left(x\right){|}^{p}+u\left(x\right)\left(Bu\right)\left(x\right)+f\left(x\right)u\left(x\right)+\frac{1}{2}\left({a}_{0}-\frac{1}{2}div\mathbf{a}\right){u}^{2}\right)dx$  and its conjugate ${\psi }^{*}$  . Let again ${\Lambda }_{1}$  be the operator ${\Lambda }_{1}\left(u\right)={\Sigma }_{i=1}^{n}{a}_{i}\frac{\partial u}{\partial {x}_{i}}+\frac{1}{2}div\left(\mathbf{a}\right)u=\Lambda \left(u\right)-\left({a}_{0}-\frac{1}{2}div\mathbf{a}\right)u.$  The functional on ${L}^{2}\left(\Omega \right)$  , is now defined as
 $\begin{array}{ccc}\stackrel{~}{I}\left(u\right)& =& \psi \left(u\right)+{\psi }^{*}\left(-{\Lambda }_{1}u-{B}_{a}u\right)\end{array}$
 $\begin{array}{ccc}& & +\frac{1}{2}{\int }_{{\Sigma }_{+}}|u\left(x\right){|}^{2}\mathbf{n}\left(x\right)\cdot \mathbf{a}\left(x\right)d\sigma -\frac{1}{2}{\int }_{{\Sigma }_{-}}|u\left(x\right){|}^{2}\mathbf{n}\left(x\right)\cdot \mathbf{a}\left(x\right)d\sigma \end{array}$
 $\begin{array}{ccc}& & +2{\int }_{{\Sigma }_{-}}u\left(x\right){u}_{0}\left(x\right)\mathbf{n}\left(x\right)\cdot \mathbf{a}\left(x\right)d\sigma -{\int }_{{\Sigma }_{-}}|{u}_{0}\left(x\right){|}^{2}\mathbf{n}\left(x\right)\cdot \mathbf{a}\left(x\right)d\sigma .\end{array}$
if $u\in {H}_{A}^{1}\left(\Omega \right)$  and $+\infty$  elsewhere. On can then verify the following.
Theorem 5.12 If $1  , then there exists $\overline{u}\in {H}_{A}^{1}\left(\Omega \right)$  such that
 $\begin{array}{c}\stackrel{~}{I}\left(\overline{u}\right)=inf\left\{\stackrel{~}{I}\left(u\right);u\in {H}_{A}^{1}\left(\Omega \right)\right\}=0.\end{array}$ (52)
and $\overline{u}$  solves equation ( 51 ).

5.5 ASD Lagrangians for coupled equations with prescribed boundaries

Assume $L\in {\mathcal{ℒ}}_{AD}\left(X\right)$  and $M\in {\mathcal{ℒ}}_{AD}\left(Y\right)$  where $X$  and $Y$  are two reflexive Banach spaces and let $A:X\to {Y}^{*}$  be any bounded linear operator. Let $\left(\Lambda ,{b}_{1},{b}_{2}\right):X\to {X}^{*}×{H}_{1}×{H}_{2}$  (resp.
$\left(\Gamma ,{c}_{1},{c}_{2}\right):Y\to {Y}^{*}×{K}_{1}×{K}_{2}$  ) be skew symmetric triplets, and let $\ell$  (resp., $m$  be a self dual boundary Lagrangian on ${H}_{1}×{H}_{2}$  (resp., ${K}_{1}×{K}_{2}$  ), in such a way that the Lagrangians ${L}_{\Lambda }\left(x,p\right)=L\left(x,\Lambda x+p\right)+\ell \left({b}_{1}\left(x\right),{b}_{2}\left(x\right)\right)and{M}_{\Gamma }\left(y,q\right)=L\left(y,\Gamma y+q\right)+\ell \left({c}_{1}\left(y\right),{b}_{2}\left(y\right)\right)$  are ASD and therefore the Lagrangian ${L}_{\Lambda }{\oplus }_{{}_{A}}{M}_{\Gamma }\left(\left(x,y\right),\left(p,q\right)\right)=L\left(x,{A}^{*}y+\Lambda x+p\right)+M\left(y,-Ax+\Gamma y+q\right)+\ell \left({b}_{1}x,{b}_{2}x\right)+m\left({c}_{1}y,{c}_{2}y\right)$  is also anti-selfdual. Consider the functional $I\left(x,y\right)={L}_{\Lambda }{\oplus }_{{}_{A}}{M}_{\Gamma }\left(\left(x,y\right),\left(0,0\right)\right)$  , that is $I\left(x,y\right):=L\left(x,{A}^{*}y+\Lambda x\right)+M\left(y,-Ax+\Gamma y\right)+\ell \left({b}_{1}x,{b}_{2}x\right)+m\left({c}_{1}y,{c}_{2}y\right).$  We can now state
Theorem 5.13 Assume that ${lim}_{\parallel x\parallel +\parallel y\parallel \to \infty }\frac{I\left(x,y\right)}{\parallel x\parallel +\parallel y\parallel }=+\infty .$  Then there exists $\left(\overline{x},\overline{y}\right)\in X×Y$  such that:
 $\begin{array}{c}I\left(\overline{x},\overline{y}\right)={inf}_{\left(x,y\right)\in X×Y}I\left(x,y\right)=0.\end{array}$ (53)
In particular, for any $a\in {H}_{1}$  and $b\in {K}_{1}$  , there exists $\left(\overline{x},\overline{y}\right)\in X×Y$  such that:
 $\begin{array}{c}\left\{\begin{array}{ccc}L\left(\overline{x},{A}^{*}\overline{y}+\Lambda \overline{x}\right)+〈\overline{x},{A}^{*}\overline{y}+\Lambda \overline{x}〉& =& 0\\ M\left(\overline{y},-A\overline{x}+\Gamma \overline{y}\right)+〈\overline{y},-A\overline{x}+\Gamma \overline{y}〉& =& 0\\ \left(-{A}^{*}\overline{y}-\Lambda \overline{x},-\overline{x}\right)& \in & \partial L\left(\overline{x},{A}^{*}\overline{x}+\Lambda \overline{x}\right)\\ \left(A\overline{x}-\Gamma \overline{y},-\overline{y}\right)& \in & \partial M\left(\overline{y},-A\overline{x}+\Gamma \overline{y}\right)\\ {b}_{1}\left(\overline{x}\right)& =& a\\ {c}_{1}\left(\overline{y}\right)& =& b\end{array}\end{array}$ (54)
It is obtained as a minimizer on $X×Y$  of the functional
 $\begin{array}{ccc}I\left(x,y\right)& =& L\left(x,{A}^{*}y+\Lambda x\right)+\frac{1}{2}\parallel {b}_{1}\left(x\right){\parallel }^{2}-2〈a,{b}_{1}\left(x\right)〉+\parallel a{\parallel }^{2}+\frac{1}{2}\parallel {b}_{2}\left(x\right){\parallel }^{2}\end{array}$
 $\begin{array}{ccc}& & +M\left(y,-Ax+\Gamma y\right)+\frac{1}{2}\parallel {c}_{1}\left(y\right){\parallel }^{2}-2〈b,{c}_{1}\left(y\right)〉+\parallel b{\parallel }^{2}+\frac{1}{2}\parallel {c}_{2}\left(y\right){\parallel }^{2}.\end{array}$
Proof: Note that we can rewrite
 $\begin{array}{ccc}I\left(x,y\right)& =& L\left(x,{A}^{*}y+\Lambda x\right)+〈x,{A}^{*}y+\Lambda x〉\end{array}$
 $\begin{array}{ccc}& & +M\left(y,-Ax+\Gamma y\right)+〈y,-Ax+\Gamma y〉\end{array}$
 $\begin{array}{ccc}& & +\parallel {b}_{1}\left(x\right)-a{\parallel }^{2}+\parallel {c}_{1}\left(x\right)-b{\parallel }^{2},\end{array}$
in such a way that if $I\left(\overline{x},\overline{y}\right)=0$  , then the fact that the sum of each two consecutive terms constituting $I$  above is non-negative, prove our claim ( 54 ).
Corollary 5.14 Let ${B}_{1}:X\to {X}^{*}$  (resp., ${B}_{2}:Y\to {Y}^{*}$  ) be positive operators modulo a regular boundary $\left({b}_{1},{b}_{2}\right):X\to {H}_{1}×{H}_{2}$  (resp., $\left({c}_{1},{c}_{2}\right):Y\to {K}_{1}×{K}_{2}$  ), let ${\phi }_{1}$  (resp ${\phi }_{2}$  ) be a convex function in $\mathcal{C}\left(X\right)$  (resp. in $\mathcal{C}\left(Y\right)$  ) and consider the convex functions ${\psi }_{1}\left(x\right)=\frac{1}{2}〈{B}_{1}x,x〉+{\phi }_{1}\left(x\right)-\frac{1}{4}\left(\parallel {b}_{2}x{\parallel }^{2}-\parallel {b}_{1}x{\parallel }^{2}\right)$  ${\psi }_{2}\left(x\right)=\frac{1}{2}〈{B}_{2}x,x〉+{\phi }_{2}\left(x\right)-\frac{1}{4}\left(\parallel {c}_{2}x{\parallel }^{2}-\parallel {c}_{1}x{\parallel }^{2}\right)$  Let $A:X\to {Y}^{*}$  be a bounded linear operator such that ${lim}_{\parallel x\parallel +\parallel y\parallel \to \infty }\frac{{\psi }_{1}\left(x\right)+{\psi }_{2}\left(y\right)}{\parallel x\parallel +\parallel y\parallel }=+\infty .$  Then, for any $\left(a,b\right)\in {H}_{1}×{K}_{1}$  , any $\left(f,g\right)\in {X}^{*}×{Y}^{*}$  and any $\alpha \in R$  , there exists a solution $\left(\overline{x},\overline{y}\right)\in X×Y$  to the system of equations
 $\begin{array}{c}\left\{\begin{array}{ccc}-{A}^{*}y-{B}_{1}x+f& \in & \partial {\phi }_{1}\left(x\right)\\ {\alpha }^{2}Ax-{B}_{2}y+g& \in & \partial {\phi }_{2}\left(y\right)\\ {b}_{1}\left(\overline{x}\right)& =& a\\ {b}_{2}\left(\overline{y}\right)& =& b\end{array}\end{array}$ (55)
It is obtained as a minimizer on $X×Y$  of the functional:
 $\begin{array}{ccc}I\left(x,y\right)& =& {\chi }_{1}\left(x\right)+{\chi }_{1}^{*}\left(-{B}_{1}^{a}x-\frac{1}{2}{b}_{2}^{*}{b}_{2}x+\frac{1}{2}{b}_{1}^{*}{b}_{1}x-{A}^{*}y\right)+\frac{1}{2}\left(\parallel {b}_{1}x{\parallel }^{2}+\parallel {b}_{2}x{\parallel }^{2}\right)-2〈a,{b}_{1}\left(x\right)〉+\parallel a{\parallel }^{2}\end{array}$
 $\begin{array}{ccc}& & +{\chi }_{2}\left(y\right)+{\chi }_{2}^{*}\left(-{B}_{2}^{a}y-\frac{1}{2}{c}_{2}^{*}{c}_{2}y+\frac{1}{2}{c}_{1}^{*}{c}_{1}y+{\alpha }^{2}Ax\right)+\frac{1}{2}\left(\parallel {c}_{1}y{\parallel }^{2}+\parallel {c}_{2}y{\parallel }^{2}\right)-2〈b,{c}_{1}y〉+\parallel b{\parallel }^{2}\end{array}$
where ${\chi }_{1}\left(x\right)={\psi }_{1}\left(x\right)-〈f,x〉$  and ${\chi }_{2}\left(y\right)={\psi }_{2}\left(y\right)-〈g,y〉$  .
Proof: Associate the following anti-selfdual Lagrangians on $X×{X}^{*}$  and $Y×{Y}^{*}$  respectively, $L\left(x,p\right)={\chi }_{1}\left(x\right)+{\chi }_{1}^{*}\left(-{B}_{1}^{a}x-\frac{1}{2}{b}_{2}^{*}{b}_{2}x+\frac{1}{2}{b}_{1}^{*}{b}_{1}x-p\right)+\frac{1}{2}\left(\parallel {b}_{1}x{\parallel }^{2}+\parallel {b}_{2}x{\parallel }^{2}\right)-2〈a,{b}_{1}\left(x\right)〉+\parallel a{\parallel }^{2}$  $M\left(y,q\right)={\chi }_{2}\left(y\right)+{\psi }_{2}^{*}\left(-{B}_{2}^{a}y-\frac{1}{2}{c}_{2}^{*}{c}_{2}y+\frac{1}{2}{c}_{1}^{*}{c}_{1}y-q\right)+\frac{1}{2}\left(\parallel {c}_{1}y{\parallel }^{2}+\parallel {c}_{2}y{\parallel }^{2}\right)-2〈b,{c}_{1}y〉+\parallel b{\parallel }^{2}$  Now apply the preceeding corollary to these two ASD Lagrangians, to the operator ${\alpha }^{2}A:X\to {Y}^{*}$  and to the product $X×Y$  equipped with the scalar product $〈\left(x,y\right),\left(p,q\right)〉=〈x,p〉+{\alpha }^{-2}〈y,q〉.$  We get
 $\begin{array}{c}\left\{\begin{array}{ccc}-{A}^{*}y-{B}_{1}^{a}x+f& \in & \partial {\phi }_{1}\left(x\right)+{B}_{1}^{s}\left(x\right)\\ {\alpha }^{2}Ax-{B}_{2}^{a}y+g& \in & \partial {\phi }_{2}\left(y\right)+{B}_{2}^{s}\left(y\right)\\ {b}_{1}\left(\overline{x}\right)& =& a\\ {b}_{2}\left(\overline{y}\right)& =& b.\end{array}\end{array}$ (56)
which gives the result.
Example 8: A variational principle for a coupled system with prescribed boundary conditions Let $\mathbf{a}:\Omega \to {\mathbf{R}}^{\mathbf{n}}$  and $\mathbf{b}:\Omega \to {\mathbf{R}}^{\mathbf{n}}$  be two smooth vector fields on a bounded domain $\Omega$  of ${\mathbf{R}}^{\mathbf{n}}$  , verifying the conditions in example 6 and consider their corresponding first order linear operator ${B}_{1}u=\mathbf{a}\cdot \nabla u$  and ${B}_{2}v=\mathbf{b}\cdot \nabla v$  . Let ${\Sigma }_{-}^{1}=\left\{x\in \partial \Omega ;\mathbf{a}\cdot \mathbf{n}\left(x\right)<0\right\}and{\Sigma }_{-}^{2}=\left\{x\in \partial \Omega ;\mathbf{b}\cdot \mathbf{n}\left(x\right)<0\right\}.$  For ${u}_{0}\in {L}_{{B}_{1}}^{2}\left({\Sigma }_{-}^{1}\right)$  and ${v}_{0}\in {L}_{{B}_{2}}^{2}\left({\Sigma }_{-}^{2}\right)$  , consider the Dirichlet problem:
 $\begin{array}{c}\left\{\begin{array}{ccc}\Delta v-\mathbf{a}\cdot \nabla u-{a}_{0}u& =& |u{|}^{p-2}u+fon\Omega \\ -{\alpha }^{2}\Delta u-\mathbf{b}\cdot \nabla v-{b}_{0}v& =& |v{|}^{q-2}v+gon\Omega \\ u& =& {u}_{0}on{\Sigma }_{-}^{1}\\ v& =& {v}_{0}on{\Sigma }_{-}^{2}.\\ \end{array}\end{array}$ (57)
We can use the above to get
Theorem 5.15 Assume ${a}_{0}\left(x\right)-\frac{1}{2}div\mathbf{a}\left(x\right)\ge \alpha >0$  and ${b}_{0}\left(x\right)-\frac{1}{2}div\mathbf{b}\left(x\right)\ge \alpha >0$  on $\Omega$  , $2  . For any $f,g\in {L}^{2}\left(\Omega \right)$  and $\left({u}_{0},{v}_{0}\right)\in {L}_{A}^{2}\left({\Sigma }_{-}^{1}\right)×{L}_{B}^{2}\left({\Sigma }_{-}^{2}\right)$  , consider on ${H}_{A}^{1}\left(\Omega \right)×{H}_{B}^{1}\left(\Omega \right)$  the functional
 $\begin{array}{ccc}I\left(u,v\right)& =& \Psi \left(u\right)+{\Psi }^{*}\left(-\mathbf{a}.\nabla u-\frac{1}{2}div\left(\mathbf{a}\right)u+\Delta v\right)\end{array}$
 $\begin{array}{ccc}& & +\frac{1}{2}{\int }_{{\Sigma }_{+}^{1}}|u\left(x\right){|}^{2}\mathbf{n}\left(x\right)\cdot \mathbf{a}\left(x\right)d\sigma -\frac{1}{2}{\int }_{{\Sigma }_{-}^{1}}|u\left(x\right){|}^{2}\mathbf{n}\left(x\right)\cdot \mathbf{a}\left(x\right)d\sigma \end{array}$
 $\begin{array}{ccc}& & +2{\int }_{{\Sigma }_{-}^{1}}u\left(x\right){u}_{0}\left(x\right)\mathbf{n}\left(x\right)\cdot \mathbf{a}\left(x\right)d\sigma -{\int }_{{\Sigma }_{-}^{1}}|{u}_{0}\left(x\right){|}^{2}\mathbf{n}\left(x\right)\cdot \mathbf{a}\left(x\right)d\sigma \end{array}$
 $\begin{array}{ccc}& & +\Phi \left(v\right)+{\Phi }^{*}\left(-\mathbf{b}.\nabla v-\frac{1}{2}div\left(\mathbf{b}\right)v-{\alpha }^{2}\Delta u\right)\end{array}$ (58)
 $\begin{array}{ccc}& & +\frac{1}{2}{\int }_{{\Sigma }_{+}^{2}}|v\left(x\right){|}^{2}\mathbf{n}\left(x\right)\cdot \mathbf{b}\left(x\right)d\sigma -\frac{1}{2}{\int }_{{\Sigma }_{-}^{2}}|v\left(x\right){|}^{2}\mathbf{n}\left(x\right)\cdot \mathbf{b}\left(x\right)d\sigma \end{array}$
 $\begin{array}{ccc}& & +2{\int }_{{\Sigma }_{-}^{2}}v\left(x\right){u}_{0}\left(x\right)\mathbf{n}\left(x\right)\cdot \mathbf{b}\left(x\right)d\sigma -{\int }_{{\Sigma }_{-}^{2}}|{v}_{0}\left(x\right){|}^{2}\mathbf{n}\left(x\right)\cdot \mathbf{b}\left(x\right)d\sigma ,\end{array}$
where $\Psi \left(u\right)=\frac{1}{p}{\int }_{\Omega }|u{|}^{p}dx+{\int }_{\Omega }fudx+\frac{1}{2}{\int }_{\Omega }\left({a}_{0}-\frac{1}{2}div\left(\mathbf{a}\right)\right)|u{|}^{2}dx,$  $\Phi \left(v\right)=\frac{1}{q}{\int }_{\Omega }|v{|}^{q}dx+{\int }_{\Omega }gvdx+\frac{1}{2}{\int }_{\Omega }\left({b}_{0}-\frac{1}{2}div\left(\mathbf{b}\right)\right)|v{|}^{2}dx,$  and ${\Psi }^{*}$  and ${\Phi }^{*}$  are their Legendre transforms. The infimum is zero and there exists a minimizer $\left(\overline{u},\overline{v}\right)\in {H}_{A}^{1}\left(\Omega \right)×{H}_{B}^{1}\left(\Omega \right)$  that is a solution of  57  $\left(\right)$  .
The conditions on $\mathbf{a}$  and $\mathbf{b}$  insure that the first order linear operators ${B}_{1}u:=\mathbf{a}\cdot \nabla u+{a}_{0}u$  (resp., ${B}_{2}v:=\mathbf{b}\cdot \nabla v+{b}_{0}v$  ) are positive modulo the boundary operators $u\to \left({u}_{{|}_{{\Sigma }_{-}^{1}}},{u}_{{|}_{{\Sigma }_{+}^{1}}}\right)\in {L}_{A}^{2}\left({\Sigma }_{-}^{1}\right)×{L}_{A}^{2}\left({\Sigma }_{+}^{1}\right)$  (resp., $v\to \left({v}_{{|}_{{\Sigma }_{-}^{2}}},{v}_{{|}_{{\Sigma }_{+}^{2}}}\right)\in {L}_{B}^{2}\left({\Sigma }_{-}^{2}\right)×{L}_{B}^{2}\left({\Sigma }_{+}^{2}\right)$  ). Apply now the above with $A={\alpha }^{2}\Delta$  .

6 Time dependent anti-self dual Lagrangians

Let $H$  be a Hilbert space with $〈,〉$  as scalar product and let $\left[0,T\right]$  be a fixed real interval where $\left(0  . Consider the classical space ${L}_{H}^{2}$  of Bochner integrable functions from $\left[0,T\right]$  into $H$  with norm denoted by $\parallel \cdot {\parallel }_{2}$  , as well as the Hilbert space ${A}_{H}^{2}=\left\{u:\left[0,T\right]\to H;\stackrel{˙}{u}\in {L}_{H}^{2}\right\}$  consisting of all absolutely continuous arcs $u:\left[0,T\right]\to H$  , equipped with the norm $\parallel u{\parallel }_{{A}_{H}^{{}^{2}}}=\left(\parallel u\left(0\right){\parallel }_{H}^{2}+{\int }_{0}^{T}\parallel \stackrel{˙}{u}{\parallel }^{2}dt{\right)}^{\frac{1}{2}}.$
Definition 6.1 Let $L:\left[0,T\right]×H×H\to IR\cup \left\{+\infty \right\}$  be measurable with respect to the $\sigma$  -field generated by the products of Lebesgue sets in $\left[0,T\right]$  and Borel sets in $H×H$  . We say that $L$  is an anti-self dual Lagrangian (ASD) on $\left[0,T\right]×H×H$  if for any $t\in \left[0,T\right]$  , the map ${L}_{t}:\left(x,p\right)\to L\left(t,x,p\right)$  is in ${\mathcal{ℒ}}_{AD}\left(H\right)$  : that is if ${L}^{*}\left(t,p,x\right)=L\left(t,-x,-p\right)forall\left(x,p\right)\in H×H.$  where here ${L}^{*}$  is the Legendre transform in the last two variables.
The most basic time-dependent $ASD$  -Lagrangians are again of the form $L\left(t,x,p\right)=\phi \left(t,x\right)+{\phi }^{*}\left(t,-p\right)$  where for each $t$  , the function $x\to \phi \left(t,x\right)$  is convex and lower semi-continuous. We now show how this property naturally “lifts” to path space.

6.1 ASD Lagrangians on path spaces

Proposition 6.1 Suppose that $L$  is an anti-self dual Lagrangian on $\left[0,T\right]×H×H$  , then for each $\omega \in \mathbf{R}$  , the Lagrangian $M\left(u,p\right):={\int }_{0}^{T}{e}^{2wt}L\left(t,{e}^{-wt}u\left(t\right),{e}^{-wt}p\left(t\right)\right)dt$  is anti-self dual on ${L}_{H}^{2}$  .
Proof: It is sufficient to show that for any Lagrangian $L\left(t,x,p\right)$  , we have the formula:
 $\begin{array}{ccc}{M}^{*}\left(p,s\right)& :=& sup\left\{{\int }_{0}^{T}\left(〈p\left(t\right),u\left(t\right)〉+〈s\left(t\right),v\left(t\right)〉-L\left(t,u\left(t\right),v\left(t\right)\right)\right)dt;\left(u,v\right)\in {L}_{H}^{2}×{L}_{H}^{2}\right\}\end{array}$
 $\begin{array}{ccc}& =& {\int }_{0}^{T}{L}^{*}\left(t,p\left(t\right),s\left(t\right)\right)dt.\end{array}$
For that, first note that for all $u,v\in {L}_{H}^{2}$  and $p,s\in {L}_{H}^{{}^{2}}$  , we have:
${\int }_{0}^{T}L\left(t,u\left(t\right),v\left(t\right)\right)dt+{\int }_{0}^{T}{L}^{*}\left(t,p\left(t\right),s\left(t\right)\right)dt\ge {\int }_{0}^{T}\left(〈p\left(t\right),u\left(t\right)〉+〈s\left(t\right),v\left(t\right)〉dt,$  which implies
 $\begin{array}{ccc}{\int }_{0}^{T}{L}^{*}\left(t,p\left(t\right),s\left(t\right)\right)dt& \ge & sup\left\{{\int }_{0}^{T}\left(〈p\left(t\right),u\left(t\right)〉+〈s\left(t\right),v\left(t\right)〉-L\left(t,u\left(t\right),v\left(t\right)\right)dt;\left(u,v\right)\in {L}_{H}^{2}×{L}_{H}^{2}\right\}\end{array}$
 $\begin{array}{ccc}& =& {M}^{*}\left(p,s\right)\end{array}$
For the reverse inequality, assume ${M}^{*}\left(p,s\right)<{\int }_{0}^{T}{L}^{*}\left(t,p\left(t\right),s\left(t\right)\right)dt$  for some $\left(p,s\right)$  in ${L}_{H}^{2}×{L}_{H}^{2}$  , and let $\mu \left(t\right)$  be such that $\mu \left(t\right)<{L}^{*}\left(t,p\left(t\right),s\left(t\right)\right)$  for all $t$  while ${\int }_{0}^{T}\mu \left(t\right)dt>{M}^{*}\left(p,s\right)$  . We then have for all $t$  , $-\mu \left(t\right)>-{L}^{*}\left(t,p\left(t\right),s\left(t\right)\right)=inf\left\{L\left(t,u,v\right)-〈u,p\left(t\right)〉-〈v,s\left(t\right)〉;\left(u,v\right)\in H×H\right\}.$  By a standard measurable selection theorem (see[10), there exists a measurable pair $\left({u}_{1},{u}_{2}\right)\in {L}_{H}^{2}×{L}_{H}^{2}$  such that $-\mu \left(t\right)\ge L\left(t,{u}_{1}\left(t\right),{u}_{2}\left(t\right)\right)-〈{u}_{1}\left(t\right),p\left(t\right)〉-〈{u}_{2}\left(t\right),s\left(t\right)〉$  . Therefore
 $\begin{array}{ccc}{M}^{*}\left(p,s\right)& <& {\int }_{0}^{T}\mu \left(t\right)dt\le {\int }_{0}^{T}-L\left(t,{u}_{1}\left(t\right),{u}_{2}\left(t\right)\right)+〈{u}_{1}\left(t\right),s\left(t\right)〉+〈{u}_{2}\left(t\right),p\left(t\right)〉dt\end{array}$
 $\begin{array}{ccc}& \le & sup\left\{{\int }_{0}^{T}\left(〈p\left(t\right),u\left(t\right)〉+〈s\left(t\right),v\left(t\right)〉-L\left(t,u\left(t\right),v\left(t\right)\right)\right)dt;\left(u,v\right)\in {L}_{H}^{2}×{L}_{H}^{2}\right\}\end{array}$
 $\begin{array}{ccc}& =& {M}^{*}\left(p,s\right)\end{array}$
A representation of $\left({A}_{H}^{{}^{2}}\right)$  : One way to represent the space ${A}_{H}^{{}^{2}}$  is to identify it with the product space $H×{L}_{H}^{2}$  , in such a way that its dual $\left({A}_{H}^{2}{\right)}^{*}$  can also be identified with $H×{L}_{H}^{2}$  via the formula:
$〈u,\left({p}_{1},{p}_{0}\right){〉}_{{}_{{A}_{H}^{2},H×{L}_{H}^{2}}}=〈u\left(0\right),{p}_{1}{〉}_{{}_{H}}+{\int }_{0}^{T}〈\stackrel{˙}{u}\left(t\right),{p}_{0}\left(t\right)〉dt.$  where $u\in {A}_{H}^{2}$  and $\left({p}_{1},{p}_{0}\right)\in H×{L}_{H}^{2}$  .
Proposition 6.2 Suppose $L$  is an anti-self dual Lagrangian on $\left[0,T\right]×H×H$  and that $\ell$  is a self-dual boundary Lagrangian on $H×H$  , then the Lagrangian defined on ${A}_{H}^{2}×\left({A}_{H}^{2}{\right)}^{*}={A}_{H}^{2}×\left(H×{L}_{H}^{2}\right)$  by $N\left(u,p\right)={\int }_{0}^{T}L\left(t,u\left(t\right)+{p}_{0}\left(t\right),\stackrel{˙}{u}\left(t\right)\right)dt+\ell \left(u\left(0\right)+{p}_{1},u\left(T\right)\right)$  is anti-selfdual on ${A}_{H}^{2}×\left({L}_{H}^{2}×\left\{0\right\}\right)$  .
Proof: For $\left(v,q\right)\in {A}_{H}^{2}×\left({A}_{H}^{2}{\right)}^{*}$  with $q$  represented by $\left({q}_{0}\left(t\right),0\right)$  write:
 $\begin{array}{ccc}{N}^{*}\left(q,v\right)& =& {sup}_{{p}_{1}\in H}{sup}_{{p}_{0}\in {L}_{H}^{2}}{sup}_{u\in {A}_{H}^{2}}\left\{〈{p}_{1},v\left(0\right)〉+{\int }_{0}^{T}\left[〈{p}_{0}\left(t\right),\stackrel{˙}{v}\left(t\right)〉+〈{q}_{0}\left(t\right),\stackrel{˙}{u}\left(t\right)〉-L\left(t,u\left(t\right)+{p}_{0}\left(t\right),\stackrel{˙}{u}\left(t\right)\right)\right]dt\end{array}$
 $\begin{array}{ccc}& & -\ell \left(u\left(0\right)+{p}_{1},u\left(T\right)\right)\right\}.\end{array}$
Making a substitution $u\left(0\right)+{p}_{1}=a\in H$  and $u\left(t\right)+{p}_{0}\left(t\right)=y\left(t\right)\in {L}_{H}^{2}$  , we obtain
 $\begin{array}{ccc}{N}^{*}\left(q,v\right)& =& {sup}_{a\in H}{sup}_{y\in {L}_{H}^{2}}{sup}_{u\in {A}_{H}^{2}}\left\{〈a-u\left(0\right),v\left(0\right)〉-\ell \left(a,u\left(T\right)\right)\end{array}$
 $\begin{array}{ccc}& & +{\int }_{0}^{T}\left[〈y\left(t\right)-u\left(t\right),\stackrel{˙}{v}\left(t\right)〉+〈{q}_{0}\left(t\right),\stackrel{˙}{u}\left(t\right)〉-L\left(t,y\left(t\right),\stackrel{˙}{u}\left(t\right)\right)\right]dt\end{array}$
Since $\stackrel{˙}{u}$  and $\stackrel{˙}{v}\in {L}_{H}^{2}$  , we have:
${\int }_{0}^{T}〈u,\stackrel{˙}{v}〉=-{\int }_{0}^{T}〈\stackrel{˙}{u},v〉+〈v\left(T\right),u\left(T\right)〉-〈v\left(0\right),u\left(0\right)〉,$  which implies
 $\begin{array}{ccc}{N}^{*}\left(q,v\right)& =& {sup}_{a\in H}{sup}_{y\in {L}_{H}^{2}}{sup}_{u\in {A}_{H}^{2}}\left\{〈a,v\left(0\right)〉-〈v\left(T\right),u\left(T\right)〉-\ell \left(a,u\left(T\right)\right)\end{array}$
 $\begin{array}{ccc}& & +{\int }_{0}^{T}\left[〈y\left(t\right),\stackrel{˙}{v}\left(t\right)〉+〈v\left(t\right)+{q}_{0}\left(t\right),\stackrel{˙}{u}\left(t\right)〉-L\left(t,y\left(t\right),\stackrel{˙}{u}\left(t\right)\right)\right]dt\right\}.\end{array}$
Identify now ${A}_{H}^{{}^{2}}$  with $H×{L}_{H}^{2}$  via the correspondence:
 $\begin{array}{ccc}\left(b,r\right)\in H×{L}_{H}^{2}& ↦& b+{\int }_{t}^{T}r\left(s\right)ds\in {A}_{H}^{{}^{2}}\end{array}$
 $\begin{array}{ccc}u\in {A}_{H}^{{}^{2}}& ↦& \left(u\left(T\right),-\stackrel{˙}{u}\left(t\right)\right)\in H×{L}_{H}^{2}.\end{array}$
We finally obtain
 $\begin{array}{ccc}{N}^{*}\left(q,v\right)& =& {sup}_{a\in H}{sup}_{b\in H}\left\{〈a,v\left(0\right)〉-〈v\left(T\right),b〉-\ell \left(a,b\right)\end{array}$
 $\begin{array}{ccc}& +& {sup}_{y\in {L}_{H}^{2}}{sup}_{r\in {L}_{H}^{2}}{\int }_{0}^{T}\left[〈y\left(t\right),\stackrel{˙}{v}\left(t\right)〉+〈v\left(t\right)+{q}_{0}\left(t\right),r\left(t\right)〉-L\left(t,y\left(t\right),r\left(t\right)\right)\right]dt\end{array}$
 $\begin{array}{ccc}& =& {\int }_{0}^{T}{L}^{*}\left(t,\stackrel{˙}{v}\left(t\right),v\left(t\right)+{q}_{0}\left(t\right)\right)dt+{\ell }^{*}\left(v\left(0\right),-v\left(T\right)\right)\end{array}$
 $\begin{array}{ccc}& =& {\int }_{0}^{T}L\left(t,-v\left(t\right)-{q}_{0}\left(t\right),-\stackrel{˙}{v}\left(t\right)\right)dt+\ell \left(-v\left(0\right),-v\left(T\right)\right)\end{array}$
 $\begin{array}{ccc}& =& N\left(-v,-q\right).\end{array}$

6.2 ASD Lagrangians in the calculus of variations

Theorem 6.2 Suppose $L$  is an anti-self dual Lagrangian on $\left[0,T\right]×H×H$  and $\ell$  is a self-dual boundary Lagrangian on $H×H$  , and consider the following functional ${I}_{\ell ,L}\left(u\right)={\int }_{0}^{T}L\left(t,u\left(t\right),\stackrel{˙}{u}\left(t\right)\right)dt+\ell \left(u\left(0\right),u\left(T\right)\right).$  Suppose there exists $C>0$  such that for all $x\in {L}_{H}^{2}$  , ${\int }_{0}^{T}L\left(t,x\left(t\right),0\right)dt\le C\left(1+\parallel x{\parallel }_{{L}_{H}^{2}}^{2}\right).$  Then there exists $v\in {A}_{H}^{2}$  such that $\left(v\left(t\right),\stackrel{˙}{v}\left(t\right)\right)\in \text{Dom}\left(L\right)$  for almost all $t\in \left[0,T\right]$  and ${I}_{\ell ,L}\left(v\right)={inf}_{u\in {A}_{H}^{2}}{I}_{\ell ,L}\left(u\right)=0.$  In particular, for every ${v}_{0}\in H$  the following functional ${I}_{\ell ,L}\left(u\right)={\int }_{0}^{T}L\left(t,u\left(t\right),\stackrel{˙}{u}\left(t\right)\right)dt+\frac{1}{2}\parallel u\left(0\right){\parallel }^{2}-2〈{v}_{0},u\left(0\right)〉+\parallel {v}_{0}{\parallel }^{2}+\frac{1}{2}\parallel u\left(T\right){\parallel }^{2}$  has minimum equal to zero on ${A}_{H}^{2}$  . It is attained at a unique path $v$  which then satisfies:
 $\begin{array}{c}v\left(0\right)={v}_{0}and\left(v\left(t\right),\stackrel{˙}{v}\left(t\right)\right)\in \text{Dom}\left(L\right)foralmostallt\in \left[0,T\right],\end{array}$ (59)
 $\begin{array}{c}\frac{d}{dt}{\partial }_{p}L\left(t,v\left(t\right),\stackrel{˙}{v}\left(t\right)\right)={\partial }_{x}L\left(t,v\left(t\right),\stackrel{˙}{v}\left(t\right)\right)\end{array}$ (60)
 $\begin{array}{c}\left(-\stackrel{˙}{v}\left(t\right),-v\left(t\right)\right)\in \partial L\left(t,v\left(t\right),\stackrel{˙}{v}\left(t\right)\right),\end{array}$ (61)
 $\begin{array}{c}\parallel v\left(t\right){\parallel }_{H}^{2}=\parallel {v}_{0}{\parallel }^{2}-2{\int }_{0}^{t}L\left(s,v\left(s\right),\stackrel{˙}{v}\left(s\right)\right)dsforeveryt\in \left[0,T\right].\end{array}$ (62)
If $L$  is autonomous and $v\in {C}^{1}\left(\left[0,T\right],H\right)$  , then for all $t\in \left[0,T\right]$  , we have:
 $\begin{array}{c}\parallel \stackrel{˙}{v}\left(t\right)\parallel \le \parallel \stackrel{˙}{v}\left(0\right)\parallel .\end{array}$ (63)
Proof: Apply Proposition 6.2 to get that $N\left(u,p\right)={\int }_{0}^{T}L\left(t,u\left(t\right)+{p}_{0}\left(t\right),\stackrel{˙}{u}\left(t\right)\right)dt+\ell \left(u\left(0\right)+{p}_{1},u\left(T\right)\right)$  is partially anti-self dual on ${A}_{H}^{2}$  . It now suffices to apply Theorem 4.1 since in this case $N\left(0,p\right)={\int }_{0}^{T}L\left(t,{p}_{0}\left(t\right),0\right)dt+\ell \left({p}_{1},0\right)\le {C}_{2}\left(1+\parallel {p}_{0}{\parallel }_{{L}_{H}^{2}}^{2}\right)+\parallel {p}_{1}{\parallel }_{H}^{2}$  , which means that $N\left(0,p\right)$  is bounded on the bounded sets of $\left({A}_{H}^{2}{\right)}^{*}$  .
For a given ${v}_{0}\in H$  , use the boundary Lagrangian $\ell \left(r,s\right)=\frac{1}{2}\parallel r{\parallel }^{2}-2〈{v}_{0},r〉+\parallel {v}_{0}{\parallel }^{2}+\frac{1}{2}\parallel s{\parallel }^{2}.$  which is clearly self-dual. We then get ${I}_{\ell ,L}\left(u\right)={\int }_{0}^{T}\left[L\left(t,u\left(t\right),\stackrel{˙}{u}\left(t\right)\right)+〈u\left(t\right),\stackrel{˙}{u}\left(t\right)〉\right]dt+\parallel u\left(0\right)-{v}_{0}{\parallel }^{2}.$  Since $L\left(t,x,p\right)\ge -〈x,p〉$  for all $\left(t,x,p\right)\in \left[0,T\right]×H×H$  , the fact that ${I}_{\ell ,L}\left(v\right)={inf}_{u\in {A}_{H}^{2}}{I}_{\ell ,L}\left(u\right)=0$  , then yields $v\left(0\right)={v}_{0}$  and that
 $\begin{array}{c}L\left(s,v\left(s\right),\stackrel{˙}{v}\left(s\right)+〈v\left(s\right),\stackrel{˙}{v}\left(s\right)〉=0foralmostalls\in \left[0,T\right].\end{array}$ (64)
This clearly yields ( 62 ), since we then have:
$\frac{d\left(|v\left(s\right){|}^{2}\right)}{ds}=-2L\left(s,v\left(s\right),\stackrel{˙}{v}\left(s\right)\right).$  To prove ( 61 ), use ( 64 ) and the fact that $L$  is anti-selfdual to write:
$L\left(s,v\left(s\right),\stackrel{˙}{v}\left(s\right)+{L}^{*}\left(s,-\stackrel{˙}{v}\left(s\right),-v\left(s\right)\right)+〈\left(v\left(s\right),\stackrel{˙}{v}\left(s\right)\right),\left(\stackrel{˙}{v}\left(s\right),v\left(s\right)\right)〉=0.$  Now apply Legendre-Fenchel duality in the space $H×H$  . The uniqueness and ( 63 ) follow from the following observation.
Lemma 6.3 Suppose $L\left(t,,\right)$  is convex on $H×H$  for each $t\in \left[0,T\right]$  , and that $x\left(t\right)$  and $v\left(t\right)$  are two paths in ${C}^{1}\left(\left[0,T\right],H\right)$  satisfying $x\left(0\right)={x}_{0}$  , $v\left(0\right)={v}_{0}$  , $-\left(\stackrel{˙}{x},x\right)\in \partial L\left(t,x,\stackrel{˙}{x}\right)$  and $-\left(\stackrel{˙}{v},v\right)\in \partial L\left(t,v,\stackrel{˙}{v}\right)$  . Then $\parallel x\left(t\right)-v\left(t\right)\parallel \le \parallel x\left(0\right)-v\left(0\right)\parallel$  for each $t\in \left[0,T\right]$  .
Proof: Estimate $\alpha \left(t\right)=\frac{d}{dt}\frac{\parallel x\left(t\right)-v\left(t\right){\parallel }^{2}}{2}$  as follows:
 $\begin{array}{ccc}\alpha \left(t\right)& =& 〈v\left(t\right)-x\left(t\right),\stackrel{˙}{v}\left(t\right)-\stackrel{˙}{x}\left(t\right)〉\end{array}$
 $\begin{array}{ccc}& =& \frac{1}{2}〈v\left(t\right)-x\left(t\right),\stackrel{˙}{v}\left(t\right)-\stackrel{˙}{x}\left(t\right)〉+\frac{1}{2}〈\stackrel{˙}{v}\left(t\right)-\stackrel{˙}{x}\left(t\right),v\left(t\right)-x\left(t\right)〉\end{array}$
 $\begin{array}{ccc}& =& \frac{1}{2}\left(〈\left(v\left(t\right)-x\left(t\right),\stackrel{˙}{v}\left(t\right)-\stackrel{˙}{x}\left(t\right)\right),\left(\stackrel{˙}{v}\left(t\right)-\stackrel{˙}{x}\left(t\right),v\left(t\right)-x\left(t\right)\right){〉}_{H×H}\right)\end{array}$
 $\begin{array}{ccc}& =& \frac{1}{2}\left(〈\left(v\left(t\right)-x\left(t\right),\stackrel{˙}{v}\left(t\right)-\stackrel{˙}{x}\left(t\right)\right),\left({L}_{x}\left(t,x\left(t\right),\stackrel{˙}{x}\left(t\right)\right)-{L}_{x}\left(t,v\left(t\right),\stackrel{˙}{v}\left(t\right)\right),{L}_{y}\left(t,x\left(t\right),\stackrel{˙}{x}\left(t\right)\right)-{L}_{y}\left(t,v\left(t\right),\stackrel{˙}{v}\left(t\right)\right)〉\right)\end{array}$
 $\begin{array}{ccc}& =& \frac{1}{2}\left(〈\left(v\left(t\right),\stackrel{˙}{v}\left(t\right)\right)-\left(x\left(t\right),\stackrel{˙}{x}\left(t\right)\right),\left({L}_{x}\left(t,x\left(t\right),\stackrel{˙}{x}\left(t\right)\right),{L}_{y}\left(t,x\left(t\right),\stackrel{˙}{x}\left(t\right)\right)\right)-\left({L}_{x}\left(t,v\left(t\right),\stackrel{˙}{v}\left(t\right)\right),{L}_{y}\left(t,v\left(t\right),\stackrel{˙}{v}\left(t\right)\right)\right)〉\right)\end{array}$
 $\begin{array}{ccc}& =& \frac{1}{2}\left(〈\left(v\left(t\right),\stackrel{˙}{v}\left(t\right)\right)-\left(x\left(t\right),\stackrel{˙}{x}\left(t\right)\right),\partial L\left(t,x\left(t\right),\stackrel{˙}{x}\left(t\right)\right)-\partial L\left(t,v\left(t\right),\stackrel{˙}{v}\left(t\right)\right)〉\right)\end{array}$
 $\begin{array}{ccc}& \le & 0\end{array}$
in view of the convexity of $L$  .
It then follows that $\parallel x\left(t\right)-v\left(t\right)\parallel \le \parallel x\left(0\right)-v\left(0\right)\parallel$  for all $t>0$  . Now if $L$  is autonomous, $v\left(t\right)$  and $x\left(t\right)=v\left(t+h\right)$  are solutions for any $h>0$  , so that ( 63 ) follows from the above.

6.3 ASD Lagrangians associated to gradient flows

The most basic example of a self-dual Lagrangian already provides a variational formulation and proof of existence for gradient flows. The following extends some of the resuts in [19.
Theorem 6.4 Let $\phi :\left[0,T\right]×H\to \mathbf{R}\cup \left\{+\infty \right\}$  be a measurable function with respect to the $\sigma$  -field in $\left[0,T\right]×H$  generated by the products of Lebesgue sets in $\left[0,T\right]$  and Borel sets in $H$  . Assume that for every $t\in \left[0,T\right]$  , the function $\phi \left(t,\cdot \right)$  is convex and lower semicontinuous on $H$  , and ${A}_{t}$  is a bounded linear positive operator on $H$  such that for some positive functions $\gamma ,{\beta }^{-1}\in {L}^{\infty }\left[0,T\right]$  , we have
 $\begin{array}{c}\beta \left(t\right)\parallel x{\parallel }^{p}\le \phi \left(t,x\right)+\frac{1}{2}〈{A}_{t}x,x〉\le \gamma \left(t\right)\parallel x{\parallel }^{q}.\end{array}$ (65)
Then, for any ${u}_{0}\in H$  , the functional
 $\begin{array}{c}I\left(u\right)=\frac{1}{2}\left(|u\left(0\right){|}^{2}+|u\left(T\right){|}^{2}\right)-2〈u\left(0\right),{u}_{0}〉+|{u}_{0}{|}^{2}+{\int }_{0}^{T}\left[\psi \left(t,u\left(t\right)\right)+{\psi }^{*}\left(t,-{A}_{t}^{a}u\left(t\right)-\stackrel{˙}{u}\left(t\right)\right)\right]dt\end{array}$ (66)
where $\psi$  is the convex functiona $\psi \left(t,x\right)=\phi \left(t,x\right)+\frac{1}{2}〈{A}_{t}x,x〉$  has a unique minimizer $v$  in ${A}_{H}^{2}$  such that:
 $\begin{array}{c}I\left(v\right)={inf}_{u\in {A}_{H}^{2}}I\left(u\right)=0.\end{array}$ (67)
Among the paths in ${A}_{H}^{2}$  , $v$  is the unique solution to
 $\begin{array}{c}\left\{\begin{array}{ccc}-{A}_{t}u\left(t\right)-\stackrel{˙}{v}\left(t\right)& \in & \partial \phi \left(t,v\left(t\right)\right)a.e.on\left[0,T\right]\\ v\left(0\right)& =& {u}_{0}.\end{array}\end{array}$ (68)
Proof: This follows directly from Theorem 6.2 applied to the anti-selfdual Lagrangian $L\left(t,x,p\right)=\psi \left(t,x\right)+{\psi }^{*}\left(t,-{A}_{t}^{a}x-p\right)$  . Note that the conditions ( 65 ) yield that ${\int }_{0}^{T}L\left(t,x\left(t\right),0\right)dt={\int }_{0}^{T}\psi \left(t,x\left(t\right)\right)+{\psi }^{*}\left(t,{A}_{t}^{a}x\left(t\right)\right)dt$  is bounded on the bounded sets of ${L}_{H}^{2}$  .

6.4 Variational resolution of parabolic equations with prescribed boundaries

Suppose now that for each $t\in \left[0,T\right]$  , $\left({b}_{1}^{t},{b}_{2}^{t}\right):{X}_{t}\to {H}_{1}^{t}×{H}_{2}^{t}$  are regular boundary operators from a reflexive Banach space ${X}_{t}$  into Hilbert spaces ${H}_{1}^{t},{H}_{2}^{t}$  , and that there are operators ${\Lambda }_{t}:{X}_{t}\to {X}_{t}^{*}$  which are skew-adjoint modulo the boundary $\left({b}_{1}^{t},{b}_{2}^{t}\right)$  , that is for every $x,y\in {X}_{t}$  , we have:
$〈{\Lambda }_{t}x,y{〉}_{{X}_{t}}=-〈{\Lambda }_{t}y,x{〉}_{{X}_{t}}+〈{b}_{2}^{t}\left(x\right),{b}_{2}^{t}\left(y\right){〉}_{{}_{{H}_{2}^{t}}}-〈{b}_{1}^{t}\left(x\right),{b}_{1}^{t}\left(y\right){〉}_{{}_{{H}_{1}^{t}}}.$  Suppose $H$  is a Hilbert space such that for each $t$  , $\left({X}_{t},{\Lambda }_{t},H\right)$  is a maximal evolution triple, in particular ${X}_{t}\subset H\subset {X}_{t}^{*}$  and ${\Lambda }_{t}:{X}_{t}\to H$  . Now starting with a time-dependent ASD Lagrangian $L$  on $H$  , and self-dual state-boundary Lagrangians ${m}_{t}:{H}_{1}^{t}×{H}_{2}^{t}\to IR\cup \left\{+\infty \right\}$  , we get by Proposition 5.4 that $M\left(t,x,p\right)=\left\{\begin{array}{ccc}L\left(t,x,{\Lambda }_{t}x+p\right)+{m}_{t}\left({b}_{1}^{t}\left(x\right),{b}_{2}^{t}\left(x\right)\right)& if& x\in {X}_{t}\\ +\infty & & otherwise\\ \end{array}$  is also anti-self dual on $H×H$  for each $t\in \left[0,T\right]$  .
If now $\ell$  is a self-dual time-boundary Lagrangian on $H$  , then $\stackrel{~}{M}\left(u,p\right)={\int }_{0}^{T}\left\{M\left(t,u\left(t\right),p\left(t\right)+\stackrel{˙}{u}\left(t\right)\right)\right\}dt+\ell \left(u\left(0\right),u\left(T\right)\right)$  is partially anti-self dual Lagrangian on ${A}_{H}^{2}$  , and Theorem 6.2 then applies to get that $I\left(u\right)=\stackrel{~}{M}\left(u,0\right)={\int }_{0}^{T}\left\{L\left(t,u\left(t\right),{\Lambda }_{t}u\left(t\right)+\stackrel{˙}{u}\left(t\right)\right)+{m}_{t}\left({b}_{1}^{t}u\left(t\right),{b}_{2}^{t}u\left(t\right)\right)\right\}dt+\ell \left(u\left(0\right),u\left(T\right)\right)$  has a minimum at $\overline{v}\left(t\right)$  , and that the minimal value is zero. Applying the theorem with the time boundary Lagrangian on $H$  , $\ell \left(x,p\right)=\frac{1}{2}\parallel x{\parallel }^{2}-2〈{v}_{0},x〉+\parallel {v}_{0}{\parallel }^{2}+\frac{1}{2}\parallel p{\parallel }^{2}$  where ${v}_{0}$  is a given initial value in $H$  , and with a state boundary Lagrangian ${m}_{t}\left(x,p\right)=\frac{1}{2}\parallel x{\parallel }^{2}-2〈\gamma \left(t\right),x〉+\parallel \gamma \left(t\right){\parallel }^{2}+\frac{1}{2}\parallel p{\parallel }^{2},$  where $\gamma \left(t\right)$  is prescribed in ${H}_{1}^{t}$  for each $t$  , we get that $\overline{v}\left(t\right)$  satisfies:
 $\begin{array}{c}\left\{\begin{array}{ccc}L\left(t,v\left(t\right),{\Lambda }_{t}v\left(t\right)+\stackrel{˙}{v}\left(t\right)\right)+〈v\left(t\right),{\Lambda }_{t}v\left(t\right)+\stackrel{˙}{v}\left(t\right)〉& =& 0a.e.t\in \left[0,T\right]\\ \left(-{\Lambda }_{t}v\left(t\right)-\stackrel{˙}{v}\left(t\right),-v\left(t\right)\right)& \in & \partial L\left(t,v\left(t\right),\stackrel{˙}{v}\left(t\right)\right)\\ {b}_{1}^{t}\left(v\left(t\right)\right)& =& \gamma \left(t\right)a.et\in \left[0,T\right]\\ v\left(0\right)& =& {v}_{0}\end{array}\end{array}$ (69)
By starting with the most basic Lagrangian $L\left(t,x,p\right)=\phi \left(t,x\right)+{\phi }^{*}\left(t,-p\right)$  we get
Theorem 6.5 Under the above conditions on $\left({X}_{t},H,{H}_{1}^{t},{H}_{2}^{t},{b}_{1}^{t},{b}_{2}^{t}\right)$  , consider bounded linear operators ${A}_{t}:{X}_{t}\to {X}_{t}^{*}$  such that ${A}_{t}-\frac{1}{2}\left(\left({b}_{2}^{t}{\right)}^{*}{b}_{2}^{t}-\left({b}_{1}^{t}{\right)}^{*}{b}_{1}^{t}\right)$  is positive and denote by ${\Lambda }_{t}$  the operator ${\Lambda }_{t}=\frac{1}{2}\left({A}_{t}-{A}_{t}^{*}\right)+\frac{1}{2}\left(\left({b}_{2}^{t}{\right)}^{*}{b}_{2}^{t}-\left({b}_{1}^{t}{\right)}^{*}{b}_{1}^{t}\right)$  which is skew-adjoint modulo the boundary.
For each $t\in \left[0,T\right]$  , suppose $\left({X}_{t},H,{\Lambda }_{t}\right)$  is a maximal evolution triple and that $\phi \left(t,\cdot \right)$  is a convex continuous function on $H$  . For $f\in {L}^{2}\left(\left[0,T\right];H\right)$  , ${v}_{0}\in H$  and $\gamma \left(t\right)\in {H}_{t}^{1}$  consider the following functional on ${A}_{H}^{2}$  ,
 $\begin{array}{ccc}I\left(u\right)& =& {\int }_{0}^{T}\left\{\psi \left(t,u\left(t\right)\right)+{\psi }^{*}\left(t,-{\Lambda }_{t}u\left(t\right)-\stackrel{˙}{u}\left(t\right)\right)+\frac{1}{2}\left(|{b}_{1}^{t}u\left(t\right){|}^{2}+|{b}_{2}^{t}u\left(t\right){|}^{2}\right)-2〈\gamma \left(t\right),u\left(t\right)〉+|\gamma \left(t\right){|}^{2}\right)\right\}dt\end{array}$
 $\begin{array}{ccc}& & +\frac{1}{2}\left(|u\left(0\right){|}^{2}+|u\left(T\right){|}^{2}\right)-2〈u\left(0\right),{v}_{0}〉+|{v}_{0}{|}^{2},\end{array}$
where $\psi \left(t,x\right)=\phi \left(t,x\right)+\frac{1}{2}〈{A}_{t}x,x〉-\frac{1}{4}\left(\parallel {b}_{2}^{t}x{\parallel }^{2}-\parallel {b}_{1}^{t}x{\parallel }^{2}\right)+〈f\left(t\right),x〉$  . Suppose there is $C>0$  so that for every $x\in {A}_{H}^{2}$  , ${\int }_{0}^{T}\psi \left(t,x\left(t\right)\right)+{\psi }^{*}\left(t,-{\Lambda }_{t}x\left(t\right)dt\le C\left(1+\parallel x{\parallel }_{{L}_{H}^{2}}^{2}\right).$  Then there exists $v\in {A}_{H}^{2}$  such that $I\left(v\right)={inf}_{u\in {A}_{H}^{2}}I\left(u\right)=0.$  Moreover, $v$  solves
 $\begin{array}{c}\left\{\begin{array}{ccc}-{A}_{t}v\left(t\right)-\stackrel{˙}{v}\left(t\right)& \in & \partial \phi \left(t,v\left(t\right)\right)+f\left(t\right)a.e.t\in \left[0,T\right]\\ {b}_{1}^{t}\left(v\left(t\right)\right)& =& \gamma \left(t\right)a.et\in \left[0,T\right]\\ v\left(0\right)& =& {v}_{0}.\end{array}\end{array}$ (70)
Example 9: Non linear Transport evolutions With the notation of Example 6, we consider the equation
 $\begin{array}{c}\left\{\begin{array}{ccc}-\frac{\partial u}{\partial t}-{\Sigma }_{i=1}^{n}{a}_{i}\frac{\partial u}{\partial {x}_{i}}-{a}_{0}u& =& \beta \left(u\right)+fon\left[0,T\right]×\Omega \\ u\left(t,x\right)& =& \gamma \left(t,x\right)on\left[0,T\right]×{\Sigma }_{-}.\\ u\left(0,x\right)& =& {u}_{0}\left(x\right)on\Omega \end{array}\end{array}$ (71)
where ${u}_{0}\in {H}_{A}^{1}\left(\Omega \right)$  , $f\in {L}^{2}\left(\Omega \right)$  and where $\gamma \left(t\right)\in {L}^{2}\left({\Sigma }_{-},|\mathbf{n}\cdot \mathbf{a}|dx\right)$  for each $t\in \left[0,T\right]$  . Let $\psi \left(u\right)={\int }_{\Omega }\left\{j\left(u\left(x\right)\right)+f\left(x\right)u\left(x\right)+\frac{1}{2}\left({a}_{0}-\frac{1}{2}\mathbf{d}\mathbf{i}\mathbf{v}\mathbf{a}\right)|u{|}^{2}\right)\right\}dx$
Theorem 6.6 Assume ${a}_{0}\left(x\right)-\frac{1}{2}\mathbf{d}\mathbf{i}\mathbf{v}a\left(x\right)\ge \alpha >0$  on $\Omega$  , and consider the following functional on the space $X:={A}^{2}\left(\left[0,T\right];{H}_{A}^{1}\left(\Omega \right)\right)$  .
 $\begin{array}{ccc}I\left(u\right)& =& {\int }_{0}^{T}\left\{\psi \left(u\left(t\right)\right)+{\psi }^{*}\left(-\mathbf{a}\cdot {\nabla }_{x}u\left(t\right)-\frac{1}{2}\left(diva\right)u\left(t\right)-\stackrel{˙}{u}\left(t\right)\right)\right\}dt\end{array}$
 $\begin{array}{ccc}& & +{\int }_{0}^{T}\left\{\frac{1}{2}{\int }_{{\Sigma }_{+}}|u\left(t,x\right){|}^{2}\mathbf{n}\cdot \mathbf{a}d\sigma -\frac{1}{2}{\int }_{{\Sigma }_{-}}|u\left(t,x\right){|}^{2}\mathbf{n}\cdot \mathbf{a}d\sigma +{\int }_{{\Sigma }_{-}}\left(|\gamma \left(t,x\right){|}^{2}-2\gamma \left(t,x\right)u\left(t,x\right)\right)|\mathbf{n}\cdot \mathbf{a}|d\sigma \right\}dt\end{array}$
 $\begin{array}{ccc}& & +{\int }_{\Omega }\left\{\frac{1}{2}\left(|u\left(0,x\right){|}^{2}+|u\left(x,T\right){|}^{2}\right)-2〈u\left(0,x\right),{u}_{0}\left(x\right)〉+|{u}_{0}\left(x\right){|}^{2}\right\}dx.\end{array}$
There exists $\overline{u}\in X$  such that $I\left(\overline{u}\right)={inf}_{u\in X}I\left(u\right)=0$  and which solves equation ( 71 ).
These results will be improved in [20.
Variational resolution for parabolic-elliptic variational inequalities Consider for each time $t$  , a bilinear continuous functional ${a}_{t}$  on a Hilbert space $H×H$  and a convex l.s.c function $\phi \left(t,\cdot \right):H\to \mathbf{R}\cup \left\{+\infty \right\}$  . Solving the corresponding parabolic variational inequality amounts to constructing for a given $f\in {L}^{2}\left(\left[0,T\right];H\right)$  and ${x}_{0}\in H$  , a path $x\left(t\right)\in {A}_{H}^{2}\left(\left[0,T\right]\right)$  such that for all $z\in H$  ,
 $\begin{array}{c}〈\stackrel{˙}{x}\left(t\right),x\left(t\right)-z\right)+{a}_{t}\left(x\left(t\right),x\left(t\right)-z\right)+\phi \left(t,x\left(t\right)\right)-\phi \left(t,z\right)\le 〈x\left(t\right)-z,f\left(t\right)〉.\end{array}$ (72)
for almost all $t\in \left[0,T\right]$  . This problem can be rewritten as: $f\left(t\right)\in \stackrel{˙}{y}\left(t\right)+{A}_{t}y\left(t\right)+\partial \phi \left(t,y\right)$  , where ${A}_{t}$  is the bounded linear operator on $H$  defined by ${a}_{t}\left(u,v\right)=〈{A}_{t}u,v〉$  . This means that the variational inequality ( 72 ) can be rewritten and solved using the variational principle in Theorem 6.5 For example, one can then solve variationally the following ”obstacle ” problem.
Corollary 6.7 Let $\left({a}_{t}{\right)}_{t}$  be bilinear continuous functionals on $H×H$  satisfying:
• For some $\lambda >0$  , ${a}_{t}\left(v,v\right)\ge \lambda \parallel v{\parallel }^{2}$  on $H$  for every $t\in \left[0,T\right]$  .
• The map $u\to {\int }_{0}^{T}{a}_{t}\left(u\left(t\right),u\left(t\right)\right)dt$  is continuous on ${L}_{H}^{2}$  .
If $K$  is a convex closed subset of $H$  , then for any $f\in {L}^{2}\left(\left[0,T\right];H\right)$  and any ${x}_{0}\in K$  , there exists a path $x\in {A}_{H}^{2}\left(\left[0,T\right]\right)$  such that $x\left(0\right)={x}_{0},x\left(t\right)\in Kforalmostallt\in \left[0,T\right]and$  $〈\stackrel{˙}{x}\left(t\right),x\left(t\right)-z〉+{a}_{t}\left(x\left(t\right),x\left(t\right)-z\right)\le 〈x\left(t\right)-z,f〉forallz\in K.$  The path $x\left(t\right)$  is obtained as a minimizer of the following functional on ${A}_{H}^{2}\left(\left[0,T\right]\right)$  :
$I\left(y\right)={\int }_{0}^{T}\left\{\phi \left(t,y\left(t\right)\right)+\left(\phi \left(t,\cdot \right)+{\psi }_{K}{\right)}^{*}\left(-\stackrel{˙}{y}\left(t\right)-{\Lambda }_{t}y\left(t\right)\right)\right\}dt+\frac{1}{2}\left(|y\left(0\right){|}^{2}+|y\left(T\right){|}^{2}\right)-2〈y\left(0\right),{x}_{0}〉+|{x}_{0}{|}^{2}.$  Here $\phi \left(t,y\right)=\frac{1}{2}{a}_{t}\left(y,y\right)-〈f\left(t\right),y〉$  and ${\psi }_{K}\left(y\right)=0$  on $K$  and $+\infty$  elsewhere, while ${\Lambda }_{t}:H\to H$  is the skew-adjoint operator defined by $〈{\Lambda }_{t}u,v〉=\frac{1}{2}\left({a}_{t}\left(u,v\right)-{a}_{t}\left(v,u\right)\right)$  .

7 Semi-groups associated to autonomous anti-selfdual Lagrangians

When the Lagrangian $L\left(x,p\right)$  is autonomous, the situation is much nicer since we can associate a flow without stringent boundedness or coercivity conditions. Indeed, we can then use a Yosida-type regularization of ASD-Lagrangian reminiscent of the standard theory for operators and for convex functions. Let us define the Partial Domain of $\partial L$  to be the set:
$Do{m}_{1}\left(\partial L\right)=\left\{x\in X;thereexistsp,q\in {X}^{*}suchthat\left(p,0\right)\in \partial L\left(x,q\right)\right\}.$  Note that if $L\left(x,p\right)=\phi \left(x\right)+{\phi }^{*}\left(-p\right)$  with $0$  assumed to be in the domain of $\partial \phi$  , then ${x}_{0}$  belongs to $Do{m}_{1}\left(\partial L\right)$  if and only if it belongs to the domain of $\partial \phi$  .
We then obtain the following result.
Theorem 7.1 Let $L$  be an anti-selfdual Lagrangian on a Hilbert space $H$  that is uniformly convex in the first variable. Assuming $Do{m}_{1}\left(\partial L\right)$  is non-empty, then there exists a semi-group of 1-Lipschitz maps $\left({T}_{t}{\right)}_{t\in {\mathbf{R}}^{+}}$  on $H$  such that ${T}_{0}=Id$  and for any ${x}_{0}\in Do{m}_{1}\left(\partial L\right)$  , the path $x\left(t\right)={T}_{t}x$  satisfies the following:
 $\begin{array}{c}\frac{d}{dt}{\partial }_{p}L\left(x\left(t\right),\stackrel{˙}{x}\left(t\right)\right)={\partial }_{x}L\left(x\left(t\right),\stackrel{˙}{x}\left(t\right)\right)\end{array}$ (73)
 $\begin{array}{c}\left(-\stackrel{˙}{x}\left(t\right),-x\left(t\right)\right)\in \partial L\left(x\left(t\right),\stackrel{˙}{x}\left(t\right)\right)\end{array}$ (74)
and
 $\begin{array}{c}\parallel x\left(t\right){\parallel }_{H}^{2}=\parallel x{\parallel }^{2}-2{\int }_{0}^{t}L\left(x\left(s\right),\stackrel{˙}{x}\left(s\right)\right)dsforeveryt\in \left[0,T\right].\end{array}$ (75)
The path $x=\left(x\left(t\right){\right)}_{t}=\left({T}_{t}x{\right)}_{t}$  is obtained as a minimizer on ${A}_{H}^{2}$  of the functional $I\left(u\right)={\int }_{0}^{T}L\left(u\left(t\right),\stackrel{˙}{u}\left(t\right)\right)dt+\frac{1}{2}\parallel u\left(0\right){\parallel }^{2}-2〈x,u\left(0\right)〉+\parallel x{\parallel }^{2}+\frac{1}{2}\parallel u\left(T\right){\parallel }^{2},$  where $I\left(x\right)={inf}_{u\in {A}_{H}^{2}}I\left(u\right)$
As mentioned above, we can associate to the Lagrangian $L\left(x,p\right)$  its $\lambda$  -regularization by considering ${L}_{\lambda }=L◆{T}_{\lambda }$  where ${T}_{\lambda }\left(x,p\right)=\frac{\parallel x{\parallel }^{2}}{2{\lambda }^{2}}+\frac{{\lambda }^{2}\parallel p{\parallel }^{2}}{2}$  . Then ${L}_{\lambda }$  satisfies the hypothesis of Theorem 6.2, and we can then find for each initial point $v\in H$  , a path ${v}_{\lambda }\in {A}_{H}^{2}$  , with ${v}_{\lambda }\left(0\right)=v$  , which verify the above properties.
The uniform convexity of $L$  in the first variable insures that the regularization ${L}_{\lambda }$  is uniformly convex in both variables which then yield ${C}^{1}$  -solutions. The 1-Lipschitz property follows from Lemma 4.3, since in the autonomous case, we can apply it to a solution $u\left(t\right)$  and its translate $v\left(t\right)=u\left(t+h\right)$  to get $\parallel u\left(t+h\right)-u\left(t\right)\parallel \le \parallel u\left(h\right)-u\left(0\right)\parallel$  for all $t$  , which yields that ${lim}_{h\to 0}\frac{\parallel u\left(t+h\right)-u\left(t\right)\parallel }{h}\le {lim}_{h\to 0}\frac{\parallel u\left(h\right)-u\left(0\right)\parallel }{h}$  . The rest of the argument amounts to analyzing what happens when $\lambda \to 0$  .
The details will be given in [20.
We can also deal with the following situation which can sometimes do away with coercivity assumptions and to also cover the case of semi-convex potentials.
Theorem 7.2 Let $L$  be an anti-selfdual Lagrangian on a Hilbert space $H$  that is uniformly convex in the first variable. Assuming $Do{m}_{1}\left(\partial L\right)$  is non-empty, then for any $\omega \in \mathbf{R}$  there exists a semi-group of maps $\left({T}_{t}{\right)}_{t\in {\mathbf{R}}^{+}}$  on $H$  such that: ${T}_{0}=Id$  and $\parallel {T}_{t}x-{T}_{t}y\parallel \le {e}^{-\omega t}\parallel x-y\parallel$  for any $x,y\in H$  .
Moreover, for any ${x}_{0}\in Do{m}_{1}\left(\partial L\right)$  the path $x\left(t\right)={T}_{t}{x}_{0}$  satisfies the following:
 $\begin{array}{ccc}-\left(\stackrel{˙}{x}\left(t\right)+\omega x\left(t\right),x\left(t\right)\right)& \in & \partial L\left(x\left(t\right),\stackrel{˙}{x}\left(t\right)+\omega x\left(t\right)\right)\end{array}$ (76)
 $\begin{array}{ccc}x\left(0\right)& =& {x}_{0}.\end{array}$
The path $x\left(t\right)$  is obtained as a minimizer on ${A}_{H}^{2}$  of the functional $\stackrel{~}{I}\left(u\right)={\int }_{0}^{T}{e}^{2\omega t}L\left(u\left(t\right),\omega u\left(t\right)+\stackrel{˙}{u}\left(t\right)\right)dt+\frac{1}{2}\parallel u\left(0\right){\parallel }^{2}-2〈{x}_{0},u\left(0\right)〉+\parallel {x}_{0}{\parallel }^{2}+\frac{1}{2}\parallel {e}^{\omega T}u\left(T\right){\parallel }^{2}$  in such a way that $\stackrel{~}{I}\left(x\right)={inf}_{u\in {A}_{H}^{2}}\stackrel{~}{I}\left(u\right)=0.$
Proof: We associate to $L$  , the anti-selfdual Lagrangian ${L}_{\omega }\left(t,x,p\right):=\left({e}^{\omega t}\cdot L\right)\left(x,p\right)={e}^{2\omega t}L\left({e}^{-\omega t}x,{e}^{-\omega t}p\right).$  Note that if $y\left(t\right)$  satisfies:
 $\begin{array}{c}\left(-\stackrel{˙}{y}\left(t\right),-y\left(t\right)\right)\in \partial {L}_{\omega }\left(t,y\left(t\right),\stackrel{˙}{y}\left(t\right)\right)\end{array}$ (77)
then $x\left(t\right)={e}^{-\omega t}y\left(t\right)$  satisfies
 $\begin{array}{c}-\left(\stackrel{˙}{x}\left(t\right)+\omega x\left(t\right),x\left(t\right)\right)\in \partial L\left(x\left(t\right),\stackrel{˙}{x}\left(t\right)+\omega x\left(t\right)\right)\end{array}$ (78)
However, we cannot apply Theorem 7.1 directly to the Lagrangian ${L}_{\omega }$  because the latter is not autonomous. However, we shall see in [20that the Yosida regularization argument still works in this case, since we have the following property:
$\left({e}^{\omega t}\cdot L\right)◆{M}_{\lambda }={e}^{\omega t}\cdot \left(L◆{M}_{\lambda }\right).$  Now we can deduce the following which was established in [19in the case of gradient flows of convex potentials (i.e., when $A=0$  and $\omega =0$  ), and in [18in the case of gradient flows of semi-convex functions (i.e., when $A=0$  and $\omega >0$  ).
Theorem 7.3 Let $\phi$  be a proper, bounded below, convex lower semi-continuous functional on $H$  such that $0\in Dom\partial \phi$  and let $A$  be a positive bounded linear operator on $H$  . For any $\omega \in \mathbf{R}$  and ${x}_{0}\in Dom\partial \phi$  , consider the following functional on ${A}_{H}^{2}$  :
 $\begin{array}{ccc}I\left(u\right)& =& {\int }_{0}^{T}{e}^{2\omega t}\left\{\psi \left(u\left(t\right)\right)+{\psi }^{*}\left(-{A}^{a}u\left(t\right)-\omega u\left(t\right)-\stackrel{˙}{u}\left(t\right)\right)\right\}dt\end{array}$
 $\begin{array}{ccc}& & +\frac{1}{2}\parallel u\left(0\right){\parallel }^{2}-2〈{x}_{0},u\left(0\right)〉+\parallel {x}_{0}{\parallel }^{2}+\frac{1}{2}\parallel {e}^{\omega T}u\left(T\right){\parallel }^{2}\end{array}$
where ${A}^{a}$  is the anti-symmetric part of $A$  , and $\psi \left(u\right)=\phi \left(u\right)+\frac{1}{2}〈Au,u〉$  . The minimum of $I$  is then zero and is attained at a path $x\left(t\right)$  which is a solution of
 $\begin{array}{c}\left\{\begin{array}{ccc}-Ax\left(t\right)-\omega x\left(t\right)-\stackrel{˙}{x}\left(t\right)& \in & \partial \phi \left(x\left(t\right)\right)a.e.t\in \left[0,T\right]\\ x\left(0\right)& =& {x}_{0}.\end{array}\end{array}$ (79)

7.1 Nonlinear parabolic equations

Example 10: Quasi-linear parabolic equations Let $\Omega$  be a smooth bounded domain in $I{R}^{n}$  . For $p\ge \frac{n-2}{n+2}$  , the Sobolev space ${W}_{0}^{1,p+1}\left(\Omega \right)\subset H:={L}^{2}\left(\Omega \right)$  , and so we define on ${L}^{2}\left(\Omega \right)$  the functional
 $\begin{array}{c}\phi \left(u\right)=\left\{\begin{array}{cc}\frac{1}{p+1}{\int }_{\Omega }|\nabla u{|}^{p+1}& on{W}_{0}^{1,p+1}\left(\Omega \right)\\ +\infty & elsewhere\end{array}\end{array}$ (80)
Its conjugate is then
 $\begin{array}{c}{\phi }^{*}\left(v\right)=\frac{p}{p+1}{\int }_{\Omega }|\nabla {\Delta }_{p}^{-1}v{|}^{\frac{p+1}{p}}dx.\end{array}$ (81)
then for any $\omega \in IR$  , any ${u}_{0}\in {W}_{0}^{1,p+1}\left(\Omega \right)$  and any $f\in {W}^{-1,\frac{p+1}{p}}\left(\Omega \right)$  , that the infimum of the functional
 $\begin{array}{ccc}I\left(u\right)& =& \frac{1}{p+1}{\int }_{0}^{T}{e}^{2\omega t}{\int }_{\Omega }\left(|\nabla u\left(t,x\right){|}^{p+1}-\left(p+1\right)f\left(x\right)u\left(x,t\right)\right)dxdt\end{array}$
 $\begin{array}{ccc}& & +p{\int }_{0}^{T}{e}^{2\omega t}{\int }_{\Omega }\left(|\nabla {\Delta }_{p}^{-1}\left(f\left(x\right)-\omega u\left(t,x\right)-\frac{\partial u}{\partial t}\left(t,x\right)\right)\right){|}^{\frac{p+1}{p}}\right)dxdt\end{array}$
 $\begin{array}{ccc}& & -2{\int }_{\Omega }u\left(0,x\right){u}_{0}\left(x\right)dx+{\int }_{\Omega }|{u}_{0}\left(x\right){|}^{2}dx+\frac{1}{2}{\int }_{\Omega }\left(|u\left(0,x\right){|}^{2}+{e}^{2T}|u\left(T,x\right){|}^{2}\right)dx\end{array}$
on the space ${A}_{H}^{2}$  is equal to zero and is attained uniquely at an ${W}_{0}^{1,p+1}\left(\Omega \right)$  -valued path $u$  such that ${\int }_{0}^{T}\parallel \stackrel{˙}{u}\left(t\right){\parallel }_{2}^{2}dt<+\infty$  and which is a solution of the equation:
 $\begin{array}{c}\left\{\begin{array}{ccc}\frac{\partial u}{\partial t}& =& {\Delta }_{p}u+\omega u+fon\Omega ×\left[0,T\right]\\ u\left(0,x\right)& =& {u}_{0}on\Omega \\ u\left(t,0\right)& =& 0on\partial \Omega .\end{array}\end{array}$ (82)
Similarly, we can deal with the equation
 $\begin{array}{c}\left\{\begin{array}{ccc}\frac{\partial u}{\partial t}\left(t,x\right)& =& {\Delta }_{p}u-Au+\omega u\left(t,x\right)+fon\Omega ×\left[0,T\right]\\ u\left(0,x\right)& =& {u}_{0}on\Omega \\ u\left(t,0\right)& =& 0on\partial \Omega .\end{array}\end{array}$ (83)
whenever $A$  is a positive operator on ${L}^{2}\left(\Omega \right)$  , Example 11: Porous media equations Let $H={H}^{-1}\left(\Omega \right)$  equipped with the norm induced by the scalar product $〈u,v〉={\int }_{\Omega }u\left(-\Delta {\right)}^{-1}vdx=〈u,v{〉}_{{H}^{-1}\left(\Omega \right)}.$  For $m\ge \frac{n-2}{n+2}$  , we have ${L}^{m+1}\left(\Omega \right)\subset {H}^{-1}$  , and so we can consider the functional
 $\begin{array}{c}\phi \left(u\right)=\left\{\begin{array}{cc}\frac{1}{m+1}{\int }_{\Omega }|u{|}^{m+1}& on{L}^{m+1}\left(\Omega \right)\\ +\infty & elsewhere\end{array}\end{array}$ (84)
and its conjugate
 $\begin{array}{c}{\phi }^{*}\left(v\right)=\frac{m}{m+1}{\int }_{\Omega }|{\Delta }^{-1}v{|}^{\frac{m+1}{m}}dx.\end{array}$ (85)
Then, for any $\omega \in IR$  , ${u}_{0}\in {H}^{-1}\left(\Omega \right)$  and $f\in {L}^{2}\left(\Omega \right)$  , the infimum of the functional
 $\begin{array}{ccc}I\left(u\right)& =& \frac{1}{m+1}{\int }_{0}^{T}{e}^{2\omega t}{\int }_{\Omega }\left(|u\left(t,x\right){|}^{m+1}dx+m|\left(-\Delta {\right)}^{-1}\left(f\left(x\right)-\omega u\left(t,x\right)-\frac{\partial u}{\partial t}\left(t,x\right)\right){|}^{\frac{m+1}{m}}\right)dxdt\end{array}$
 $\begin{array}{ccc}& & -{\int }_{0}^{T}{e}^{2\omega t}{\int }_{\Omega }u\left(x,t\right)\left(\left(-\Delta {\right)}^{-1}f\right)\left(x\right)dxdt+{\int }_{\Omega }|\nabla \left(-\Delta {\right)}^{-1}{u}_{0}\left(x\right){|}^{2}dx\end{array}$
 $\begin{array}{ccc}& & -2{\int }_{\Omega }{u}_{0}\left(x\right)\left(-\Delta {\right)}^{-1}u\left(0,x\right)dx+\frac{1}{2}\left(\parallel u\left(0\right){\parallel }_{{}_{{H}^{-1}}}^{2}+{e}^{2T}\parallel u\left(T\right){\parallel }_{{}_{{H}^{-1}}}^{2}\right)\end{array}$
on the space ${A}_{H}^{2}$  is equal to zero and is attained uniquely at an ${L}^{m+1}\left(\Omega \right)$  -valued path $u$  such that ${\int }_{0}^{T}\parallel \stackrel{˙}{u}\left(t\right){\parallel }_{H}^{2}dt<+\infty$  and which is a solution of the equation:
 $\begin{array}{c}\left\{\begin{array}{ccc}\frac{\partial u}{\partial t}\left(t,x\right)& =& \Delta {u}^{m}+\omega u\left(t,x\right)+fon\Omega ×\left[0,T\right]\\ u\left(0,x\right)& =& {u}_{0}on\Omega .\\ \end{array}\end{array}$ (86)

7.2 Variational resolution for coupled flows and wave-type equations

Again, ASD Lagrangians are suited to treat variationally coupled evolution equations.
Proposition 7.1 Let $\phi$  be a proper convex lower semi-continuous function on $X×Y$  and let $A:X\to {Y}^{*}$  be any bounded linear operator. Assume ${B}_{1}:X\to X$  (resp., ${B}_{2}:Y\to Y$  ) are positive operators, then for any $\left({x}_{0},{y}_{0}\right)\in dom\left(\partial \phi \right)$  and any $\left(f,g\right)\in {X}^{*}×{Y}^{*}$  , there exists a path $\left(x\left(t\right),y\left(t\right)\right)\in {A}_{X}^{2}×{A}_{Y}^{2}$  such that
 $\begin{array}{ccc}-\stackrel{˙}{x}\left(t\right)-{A}^{*}y\left(t\right)-{B}_{1}x\left(t\right)+f& \in & {\partial }_{1}\phi \left(x\left(t\right),y\left(t\right)\right)\end{array}$
 $\begin{array}{ccc}-\stackrel{˙}{y}\left(t\right)+Ax\left(t\right)-{B}_{2}y\left(t\right)+g& \in & {\partial }_{2}\phi \left(x\left(t\right),y\left(t\right)\right)\end{array}$
 $\begin{array}{ccc}x\left(0\right)& =& {x}_{0}\end{array}$
 $\begin{array}{ccc}y\left(0\right)& =& {y}_{0}.\end{array}$
The solution is obtained as a minimizer on ${A}_{X}^{2}×{A}_{Y}^{2}$  of the following functional
 $\begin{array}{ccc}I\left(x,y\right)& =& {\int }_{0}^{T}\left\{\psi \left(x\left(t\right),y\left(t\right)\right)+{\psi }^{*}\left(-{A}^{*}y\left(t\right)-{B}_{1}^{a}x\left(t\right)-\stackrel{˙}{x}\left(t\right),Ax\left(t\right)-{B}_{2}^{a}y\left(t\right)-\stackrel{˙}{y}\left(t\right)\right)\right\}dt\end{array}$
 $\begin{array}{ccc}& & +\frac{1}{2}\parallel x\left(0\right){\parallel }^{2}-2〈{x}_{0},x\left(0\right)〉+\parallel {x}_{0}{\parallel }^{2}+\frac{1}{2}\parallel x\left(T\right){\parallel }^{2}\end{array}$
 $\begin{array}{ccc}& & +\frac{1}{2}\parallel y\left(0\right){\parallel }^{2}-2〈{y}_{0},y\left(0\right)〉+\parallel {y}_{0}{\parallel }^{2}+\frac{1}{2}\parallel y\left(T\right){\parallel }^{2}.\end{array}$
whose infimum is zero. Here ${B}_{1}^{a}$  (resp., ${B}_{2}^{a}$  ) are the skew-symmetric parts of ${B}_{1}$  and ${B}_{2}$  and $\psi \left(x,y\right)=\phi \left(x,y\right)+\frac{1}{2}〈{B}_{1}x,x〉-〈f,x〉+\frac{1}{2}〈{B}_{2}y,y〉-〈g,x〉$
Proof: It is enough to apply Theorem 7.1 to the ASD Lagrangian $L\left(\left(x,y\right),\left(p,q\right)\right)=\psi \left(x,y\right)+{\psi }^{*}\left(-{A}^{*}y-{B}_{1}^{a}x-p,Ax-{B}_{2}^{a}y-q\right).$  obtained by shifting to the right the ASD Lagrangian $\phi {\oplus }_{as}A$  by the skew-adjoint operator $\left({B}_{1}^{a},{B}_{2}^{a}\right)$  .
If $\left(\overline{x}\left(t\right),\overline{y}\left(t\right)\right)$  is where the infimum is attained, then we get
 $\begin{array}{ccc}0& =& I\left(\overline{x},\overline{y}\right)\end{array}$
 $\begin{array}{ccc}& =& {\int }_{0}^{T}\left\{\psi \left(\overline{x}\left(t\right),\overline{y}\left(t\right)\right)+{\psi }^{*}\left(-{A}^{*}\overline{y}\left(t\right)-{B}_{1}^{a}\overline{x}\left(t\right)-\stackrel{˙}{\overline{x}}\left(t\right),A\overline{x}\left(t\right)-{B}_{2}^{a}\overline{y}\left(t\right)-\stackrel{˙}{\overline{y}}\left(t\right)\right)\end{array}$
 $\begin{array}{ccc}& & -〈\left(\overline{x}\left(t\right),\overline{y}\left(t\right)\right),\left(-{A}^{*}\overline{y}\left(t\right)-{B}_{1}^{a}\overline{x}\left(t\right)-\stackrel{˙}{\overline{x}}\left(t\right),A\overline{x}\left(t\right)-{B}_{2}^{a}\overline{y}\left(t\right)-\stackrel{˙}{\overline{y}}\left(t\right)\right)〉\right\}dt\end{array}$
 $\begin{array}{ccc}& & +\parallel x\left(0\right)-{x}_{0}{\parallel }^{2}++\parallel y\left(0\right)-{y}_{0}{\parallel }^{2}\end{array}$
It follows that $\overline{x}\left(0\right)={x}_{0}$  , $\overline{y}\left(0\right)=0$  and the integrand is zero for almost all $t$  which yields
 $\begin{array}{ccc}-\stackrel{˙}{x}\left(t\right)-{A}^{*}y\left(t\right)-{B}_{1}^{a}x\left(t\right)& \in & {\partial }_{1}\psi \left(x\left(t\right),y\left(t\right)\right)={\partial }_{1}\phi \left(x\left(t\right),y\left(t\right)\right)+{B}_{1}^{s}x\left(t\right)-f\end{array}$
 $\begin{array}{ccc}-\stackrel{˙}{y}\left(t\right)+Ax\left(t\right)-{B}_{2}^{a}y\left(t\right)& \in & {\partial }_{2}\psi \left(x\left(t\right),y\left(t\right)\right)={\partial }_{2}\phi \left(x\left(t\right),y\left(t\right)\right)+{B}_{2}^{s}y\left(t\right)-g\end{array}$
 $\begin{array}{ccc}x\left(0\right)& =& {x}_{0}\end{array}$
 $\begin{array}{ccc}y\left(0\right)& =& {y}_{0}.\end{array}$
Consider now two convex lower semi-continuous ${\phi }_{1}$  and ${\phi }_{2}$  on Hilbert spaces $X$  and $Y$  respectively, as well as two positive operators ${B}_{1}$  on $X$  and ${B}_{2}$  on $Y$  . For any $\left(f,g\right)\in X×Y$  , consider the convex functionals ${\psi }_{1}\left(x\right)=\frac{1}{2}〈{B}_{1}x,x〉+{\phi }_{1}\left(x\right)$  and ${\psi }_{2}\left(x\right)=\frac{1}{2}〈{B}_{2}x,x〉+{\phi }_{2}\left(x\right)$  , and the anti-selfdual Lagrangians $L\left(x,p\right)={\psi }_{1}\left(x\right)-〈f,x〉+{\psi }_{1}^{*}\left(-{B}_{1}^{a}x+f-p\right),for\left(x,p\right)\in X×X,$  and $M\left(y,q\right)={\psi }_{2}\left(y\right)-〈g,y〉+{\psi }_{2}^{*}\left(-{B}_{2}^{a}y+g-q\right),for\left(y,q\right)\in Y×Y,$  For $w,{w}^{\prime }\in \mathbf{R}$  , we associate the following time-dependent ASD Lagrangian: ${L}_{\omega }\left(t,x,p\right):={e}^{-2wt}L\left({e}^{wt}x,{e}^{wt}p\right)and{M}_{{\omega }^{\prime }}\left(t,y,q\right)={e}^{-2{w}^{\prime }t}M\left({e}^{{w}^{\prime }t}y,{e}^{{w}^{\prime }t}q\right).$  Let $A:X\to Y$  be any bounded linear operator and consider for any $c\in IR$  the following twisted ASD Lagrangian on $X×Y$  $\left({L}_{\omega }{\oplus }_{{c}^{2}A}{M}_{{\omega }^{\prime }}\right)\left(t,\left(x,y\right),\left(p.q\right)\right):={L}_{\omega }\left(t,x,{A}^{*}y+p\right)+{M}_{{\omega }^{\prime }}\left(t,y,-{c}^{2}Ax+q\right).$  where the duality in $X×Y$  is given by $〈\left(x,y\right),\left(p,q\right)〉=〈x,p〉+{c}^{-2}〈y,q〉$  . Applying Theorem 7.2, we obtain
Proposition 7.2 Assume $0\in Dom\left(\partial {\phi }_{1}\right)$  and $0\in Dom\left(\partial {\phi }_{2}\right)$  , and consider the following functional on ${A}_{X}^{2}×{A}_{Y}^{2}$  :
 $\begin{array}{ccc}I\left(u,v\right)& =& {\int }_{0}^{T}{e}^{-2\omega t}\left\{{\psi }_{1}\left({e}^{\omega t}u\left(t\right)\right)+{\psi }_{1}^{*}\left({e}^{\omega t}\left(-{A}^{*}v\left(t\right)-{B}_{1}^{a}u\left(t\right)-\stackrel{˙}{u}\left(t\right)\right)\right\}dt\end{array}$
 $\begin{array}{ccc}& & +{\int }_{0}^{T}{e}^{-2{\omega }^{\prime }t}\left\{{\psi }_{2}\left({e}^{{\omega }^{\prime }t}v\left(t\right)\right)+{\psi }_{2}^{*}\left({e}^{{\omega }^{\prime }t}\left({c}^{2}Au\left(t\right)-{B}_{2}^{a}v\left(t\right)-\stackrel{˙}{v}\left(t\right)\right)\right\}dt\end{array}$
 $\begin{array}{ccc}& & +\frac{1}{2}\parallel u\left(0\right){\parallel }^{2}-2〈{x}_{0},u\left(0\right)〉+\parallel {x}_{0}{\parallel }^{2}+\frac{1}{2}\parallel u\left(T\right){\parallel }^{2}\end{array}$
 $\begin{array}{ccc}& & +\frac{1}{2}\parallel v\left(0\right){\parallel }^{2}-2〈{y}_{0},v\left(0\right)〉+\parallel {y}_{0}{\parallel }^{2}+\frac{1}{2}\parallel v\left(T\right){\parallel }^{2}.\end{array}$
The minimum of $I$  is then zero and is attained at a path $\left(\overline{x}\left(t\right),\overline{y}\left(t\right)$  , in such a way that $x\left(t\right)={e}^{\omega t}\overline{x}\left(t\right)$  and $y\left(t\right)={e}^{{\omega }^{\prime }t}\overline{y}\left(t\right)$  form a solution of the system of equations
 $\begin{array}{c}\left\{\begin{array}{ccc}-\stackrel{˙}{x}\left(t\right)+\omega x\left(t\right)-{A}^{*}y\left(t\right)-{B}_{1}x\left(t\right)+f& \in & \partial {\phi }_{1}\left(x\left(t\right)\right)\\ -\stackrel{˙}{y}\left(t\right)+{\omega }^{\prime }y\left(t\right)+{c}^{2}Ax\left(t\right)-{B}_{2}y\left(t\right)+g& \in & \partial {\phi }_{2}\left(y\left(t\right)\right)\\ x\left(0\right)& =& {x}_{0}\\ y\left(0\right)& =& {y}_{0}.\end{array}\end{array}$ (87)
Example 12: A variational principle for coupled equations Let ${\mathbf{b}}_{1}:\Omega \to {\mathbf{R}}^{\mathbf{n}}$  and ${\mathbf{b}}_{2}:\Omega \to {\mathbf{R}}^{\mathbf{n}}$  be two smooth vector fields on a bounded domain $\Omega$  of ${\mathbf{R}}^{\mathbf{n}}$  , verifying the conditions in example 3 and consider their corresponding first order linear operator ${B}_{1}v={\mathbf{b}}_{1}\cdot \nabla vand{B}_{2}v={\mathbf{b}}_{2}\cdot \nabla v.$  Consider the Dirichlet problem:
 $\begin{array}{c}\left\{\begin{array}{ccc}-\frac{\partial u}{\partial t}-\Delta \left(v-u\right)+{\mathbf{b}}_{1}\cdot \nabla u& =& |{u}^{p-2}|u+fon\left(0,T\right]×\Omega \\ -\frac{\partial v}{\partial t}+\Delta \left(v+{c}^{2}u\right)+{\mathbf{b}}_{2}\cdot \nabla v& =& |{v}^{q-2}|v+gon\left(0,T\right]×\Omega \\ u\left(t,x\right)=v\left(t,x\right)& =& 0on\left(0,T\right]×\partial \Omega .\\ u\left(0,x\right)& =& {u}_{0}\left(x\right)forx\in \Omega .\\ v\left(0,x\right)& =& {v}_{0}\left(x\right)forx\in \Omega .\end{array}\end{array}$ (88)
We can use the above to get
Theorem 7.4 Assume $div\left({\mathbf{b}}_{1}\right)\ge 0$  and $div\left({\mathbf{b}}_{2}\right)\ge 0$  on $\Omega$  , $1  and consider on ${A}_{{}_{{H}_{0}^{1}\left(\Omega \right)\right)}}^{2}×{A}_{{}_{{H}_{0}^{1}\left(\Omega \right)}}^{2}$  the functional
 $\begin{array}{ccc}I\left(u,v\right)& =& {\int }_{0}^{T}\left\{\Psi \left(u\left(t\right)\right)+{\Psi }^{*}\left({\mathbf{b}}_{1}.\nabla u\left(t\right)+\frac{1}{2}div\left({\mathbf{b}}_{1}\right)u\left(t\right)-\Delta v\left(t\right)-\stackrel{˙}{u}\left(t\right)\right)\right\}dt\end{array}$
 $\begin{array}{ccc}& & +{\int }_{0}^{T}\left\{\Phi \left(v\left(t\right)\right)+{\Phi }^{*}\left({\mathbf{b}}_{2}.\nabla v\left(t\right)+\frac{1}{2}div\left({\mathbf{b}}_{2}\right)v\left(t\right)+{c}^{2}\Delta u\left(t\right)-\stackrel{˙}{v}\left(t\right)\right)\right\}dt\end{array}$
 $\begin{array}{ccc}& & +{\int }_{\Omega }\left\{\frac{1}{2}\left(|u\left(0,x\right){|}^{2}+|u\left(T,x\right){|}^{2}\right)-2u\left(0,x\right){u}_{0}\left(x\right)+|{u}_{0}\left(x\right){|}^{2}\right\}dx\end{array}$
 $\begin{array}{ccc}& & +{\int }_{\Omega }\left\{\frac{1}{2}\left(|v\left(0,x\right){|}^{2}+|v\left(T,x\right){|}^{2}\right)-2v\left(0,x\right){v}_{0}\left(x\right)+|{v}_{0}\left(x\right){|}^{2}\right\}dx\end{array}$
where $\Psi \left(u\right)=\frac{1}{2}{\int }_{\Omega }|\nabla u{|}^{2}dx+\frac{1}{p}{\int }_{\Omega }|u{|}^{p}dx+{\int }_{\Omega }fudx+\frac{1}{4}{\int }_{\Omega }div\left({\mathbf{b}}_{1}\right)|u{|}^{2}dx,$  $\Phi \left(v\right)=\frac{1}{2}{\int }_{\Omega }|\nabla v{|}^{2}dx+\frac{1}{q}{\int }_{\Omega }|v{|}^{q}dx+{\int }_{\Omega }gvdx+\frac{1}{4}{\int }_{\Omega }div\left({\mathbf{b}}_{2}\right)|v{|}^{2}dx$  and ${\Psi }^{*}$  and ${\Phi }^{*}$  are their Legendre transforms. Then there exists $\left(\overline{u},\overline{v}\right)\in {A}_{{}_{{H}_{0}^{1}\left(\Omega \right)\right)}}^{2}×{A}_{{}_{{H}_{0}^{1}\left(\Omega \right)}}^{2}$  such that:
$I\left(\overline{u},\overline{v}\right)=inf\left\{I\left(u,v\right);\left(u,v\right)\in {A}_{{}_{{H}_{0}^{1}\left(\Omega \right)\right)}}^{2}×{A}_{{}_{{H}_{0}^{1}\left(\Omega \right)}}^{2}\right\}=0,$  and $\left(\overline{u},\overline{v}\right)$  is a solution of  88  $\left(\right)$  .
Example 13: Pressureless gaz of sticky particles Motivated by the recent work of Brenier [7we consider equations of the form
 $\begin{array}{c}{\partial }_{tt}X={c}^{2}{\partial }_{yy}X-{\partial }_{t}{\partial }_{a}\mu ,{\partial }_{a}X\ge 0,\mu \ge 0.\end{array}$ (89)
where here $X\left(t\right):=X\left(t,a,y\right)$  is a function on $K=\left[0,1\right]×\mathbf{R}/\mathbf{Z}$  , and $\mu \left(t,a,y\right)$  is a nonnegative measure that plays the role of a Lagrange multiplier for the constraint ${\partial }_{a}X\ge 0$  . Following Brenier, we reformulate the problem with the following system:
 $\begin{array}{c}\left\{\begin{array}{ccc}-\stackrel{˙}{X}\left(t\right)-\frac{\partial U}{\partial y}\left(t\right)& \in & \partial {\phi }_{1}\left(X\left(t\right)\right)\\ -\stackrel{˙}{U}\left(t\right)+\frac{\partial X}{\partial y}\left(t\right)& =& 0\\ X\left(0\right)& =& {X}_{0}\\ U\left(0\right)& =& {U}_{0}\end{array}\end{array}$ (90)
where ${\phi }_{1}$  is the convex function defined on ${L}^{2}\left(K\right)$  by
 $\begin{array}{c}{\phi }_{1}\left(X\right)=\left\{\begin{array}{cc}0& if{\partial }_{a}X\ge 0\\ +\infty & elsewhere\end{array}\end{array}$ (91)
We can solve this system with the above method by first setting ${\phi }_{2}\left(U\right)=0$  for every $U\in {L}^{2}\left(K\right)$  and by considering the Hilbert spaces $X=Y={H}_{per}^{2}\left(K\right)$  to be the subspace of ${A}_{K}^{2}$  consisting of functions that are periodic in $y$  . Define on this space the operator $AX=\frac{\partial X}{\partial y}$  in such a way that ${A}^{*}=-A$  . We consider now the functional
 $\begin{array}{ccc}I\left(X,U\right)& =& {\int }_{0}^{T}\left\{{\phi }_{1}\left(X\left(t\right)\right)+{\phi }_{1}^{*}\left(-\frac{\partial U}{\partial y}\left(t\right)-\stackrel{˙}{X}\left(t\right)\right)\right\}dt\end{array}$
 $\begin{array}{ccc}& & +{\int }_{0}^{T}\left\{{\phi }_{2}^{*}\left(\frac{\partial X}{\partial y}\left(t\right)-\stackrel{˙}{U}\left(t\right)\right)\right\}dt\end{array}$
 $\begin{array}{ccc}& & +\frac{1}{2}\parallel X\left(0\right){\parallel }^{2}-2〈{X}_{0},X\left(0\right)〉+\parallel {X}_{0}{\parallel }^{2}+\frac{1}{2}\parallel X\left(T\right){\parallel }^{2}\end{array}$
 $\begin{array}{ccc}& & +\frac{1}{2}\parallel U\left(0\right){\parallel }^{2}-2〈{Y}_{0},U\left(0\right)〉+\parallel {Y}_{0}{\parallel }^{2}+\frac{1}{2}\parallel U\left(T\right){\parallel }^{2}\end{array}$
If $\left({X}_{0},{U}_{0}\right)$  are such that ${\partial }_{a}{X}_{0}\ge 0$  , then the minimum of $I$  is then zero and is attained at a path $\left(\overline{X}\left(t\right),\overline{U}\left(t\right)$  , which solves the above system of equations.

8 Variational resolution of certain implicit PDEs

Motivated by the time dependent case, we briefly describe in this section how ASD Lagrangians can be used to solve variationally certain types of implicit differential equations. Indeed, letting $\Omega$  be a bounded domain in $I{R}^{n}$  , we let $\mu ={\mu }^{+}-{\mu }^{-}$  be a signed finite measure on the boundary $\partial \Omega$  .
We consider an equation of the form
 $\begin{array}{ccc}L\left(x,u\left(x\right),\Lambda u\left(x\right)\right)+〈u\left(x\right),\left(\Lambda u\right)\left(x\right){〉}_{{}_{X,{X}^{*}}}& =& 0on\Omega \end{array}$ (92)
 $\begin{array}{ccc}u& =& 0{\mu }^{-}a.e.on\Sigma .\end{array}$
where $L:\Omega ×X×{X}^{*}\to IR$  for some Banach space $X$  and where $\Lambda$  is an operator from a space $H\left(\Omega ;X\right)$  of $X$  -valued functions on $\Omega$  to a space $K\left(\Omega ;{X}^{*}\right)$  of ${X}^{*}$  -valued functions on $\Omega$  .
Suppose now that for each $x\in \Omega$  , the Lagrangian $L\left(x,\cdot ,\cdot \right)$  is anti-selfdual on $X×{X}^{*}$  and suppose that $\Lambda$  is skew-adjoint modulo the boundary $\Sigma$  in the following sense:
There exists a coercive map $J:X\to {X}^{*}$  (normally a duality map) such that for any $u,v$  in $H\left(\Omega ;X\right)$  we have:
${\int }_{\Omega }〈\left(\Lambda u\right)\left(x\right),v\left(x\right){〉}_{{}_{X,{X}^{*}}}dx=-{\int }_{\Omega }〈\left(\Lambda v\right)\left(x\right),u\left(x\right){〉}_{{}_{X,{X}^{*}}}dx+{\int }_{\partial \Omega }〈Jv\left(x\right),u\left(x\right){〉}_{{}_{X,{X}^{*}}}d\mu .$  On can then consider the Lagrangian $\mathcal{ℒ}\left(u,p\right)={\int }_{\Omega }L\left(x,u\left(x\right),\Lambda u\left(x\right)+p\left(x\right)\right)dx+\frac{1}{2}{\int }_{\partial \Omega }〈Ju\left(x\right),u\left(x\right){〉}_{{}_{X,{X}^{*}}}d{\mu }^{+}+\frac{1}{2}{\int }_{\partial \Omega }〈Ju\left(x\right),u\left(x\right){〉}_{{}_{X,{X}^{*}}}d{\mu }^{-}$  on $H\left(\Omega ;X\right)×K\left(\Omega ;{X}^{*}\right)$  . Under the right conditions, the Lagrangian $\mathcal{ℒ}$  has every chance to be anti-selfdual and therefore Theorem 4.1 can apply to yield that the infimum of the functional $I\left(u\right):=\mathcal{ℒ}\left(u,0\right)={\int }_{\Omega }L\left(x,u\left(x\right),\Lambda u\left(x\right)\right)dx+\frac{1}{2}{\int }_{\partial \Omega }〈Ju\left(x\right),u\left(x\right){〉}_{{}_{X,{X}^{*}}}d{\mu }^{+}+\frac{1}{2}{\int }_{\partial \Omega }〈Ju\left(x\right),u\left(x\right){〉}_{{}_{X,{X}^{*}}}d{\mu }^{-}.$  is zero and that there exists $\overline{u}\in H\left(\Omega ;X\right)$  such that $I\left(\overline{u}\right)=inf\left\{I\left(u\right);u\in {H}_{A}^{1}\left(\Omega \right)\right\}=0.$  It follows that
 $\begin{array}{ccc}0& =& {\int }_{\Omega }L\left(x,u\left(x\right),\Lambda u\left(x\right)\right)dx+\frac{1}{2}{\int }_{\partial \Omega }〈Ju\left(x\right),u\left(x\right){〉}_{{}_{X,{X}^{*}}}d{\mu }^{+}+\frac{1}{2}{\int }_{\partial \Omega }〈Ju\left(x\right),u\left(x\right){〉}_{{}_{X,{X}^{*}}}d{\mu }^{-}\end{array}$
 $\begin{array}{ccc}& =& {\int }_{\Omega }\left(L\left(x,u\left(x\right),\Lambda u\left(x\right)\right)+〈u\left(x\right),\left(\Lambda u\right)\left(x\right){〉}_{{}_{X,{X}^{*}}}\right)dx+\frac{1}{2}{\int }_{\partial \Omega }〈Ju\left(x\right),u\left(x\right){〉}_{{}_{X,{X}^{*}}}d{\mu }^{-}\end{array}$
Since both integrands are non-negative and $J$  is coercive we deduce that $\overline{u}$  satisfies equation ( 92 ).
The key here is that “pointwise” the operator $\Lambda$  may not have any particular property but on the “average” it is skew-adjoint on the function space and this allows for the variational approach to apply. We illustrate the method on the following example.
Example 14: An implicit transport equation Let $\Omega$  be a bounded domain in $I{R}^{n}$  and let $\mathbf{a}$  be a vector field as in Example (6) as well as its entrance set ${\Sigma }_{-}$  . Consider the following “toy” equation for a function $u:\Omega \to IR$  :
 $\begin{array}{ccc}\frac{1}{2}{u}^{2}+u+\frac{1}{2}\left(-\mathbf{a}\cdot \nabla u-1{\right)}^{2}-\left(\mathbf{a}\cdot \nabla u\right)u& =& 0on\Omega \end{array}$ (93)
 $\begin{array}{ccc}u& =& 0on{\Sigma }_{-}.\end{array}$
Assuming that $div\mathbf{a}=0$  , we can solve the above equation by minimizing the following functional $I\left(u\right)={\int }_{\Omega }\left(\frac{1}{2}{u}^{2}+u+\frac{1}{2}\left(-\mathbf{a}\cdot \nabla u-1{\right)}^{2}\right)dx+\frac{1}{2}{\int }_{{\Sigma }_{+}}|u\left(x\right){|}^{2}\mathbf{n}\left(x\right)\cdot \mathbf{a}\left(x\right)d\sigma -\frac{1}{2}{\int }_{{\Sigma }_{-}}|u\left(x\right){|}^{2}\mathbf{n}\left(x\right)\cdot \mathbf{a}\left(x\right)d\sigma$  on the space ${H}_{A}^{1}\left(\Omega \right)$  . Indeed, we can rewrite
 $\begin{array}{ccc}I\left(u\right)& =& \psi \left(u\right)+{\psi }^{*}\left(-\Lambda u\right)+\frac{1}{2}{\int }_{{\Sigma }_{+}}|u\left(x\right){|}^{2}\mathbf{n}\left(x\right)\cdot \mathbf{a}\left(x\right)d\sigma -\frac{1}{2}{\int }_{{\Sigma }_{-}}|u\left(x\right){|}^{2}\mathbf{n}\left(x\right)\cdot \mathbf{a}\left(x\right)d\sigma \end{array}$
where $\psi$  is the convex functional on ${L}^{2}\left(\Omega \right)$  defined by:
$\psi \left(u\right)=\frac{1}{2}{\int }_{\Omega }\left(|u\left(x\right){|}^{2}+u\left(x\right)\right)dx$  and ${\psi }^{*}$  is its Legendre conjugate ${\psi }^{*}\left(v\right)=\frac{1}{2}{\int }_{\Omega }\left(v\left(x\right)-1{\right)}^{2}dx.$  Note that $\Lambda u=\mathbf{a}\cdot \nabla u$  is here skew-adjoint modulo the boundary. There exists then a solution $\overline{u}$  for ( 93 ) that is obtained as a minimizer of the problem:
$I\left(\overline{u}\right)=inf\left\{I\left(u\right);u\in {H}_{A}^{1}\left(\Omega \right)\right\}=0.$  Rewrite now
 $\begin{array}{ccc}0& =& I\left(\overline{u}\right)={\int }_{\Omega }\left(\frac{1}{2}{\overline{u}}^{2}+\overline{u}+\frac{1}{2}\left(-\mathbf{a}\cdot \nabla \overline{u}-1{\right)}^{2}\right)dx+\frac{1}{2}{\int }_{{\Sigma }_{+}}|\overline{u}{|}^{2}\mathbf{n}\cdot \mathbf{a}d\sigma -\frac{1}{2}{\int }_{{\Sigma }_{-}}|\overline{u}{|}^{2}\mathbf{n}\cdot \mathbf{a}d\sigma \end{array}$
 $\begin{array}{ccc}& =& {\int }_{\Omega }\left(\frac{1}{2}{\overline{u}}^{2}+\overline{u}+\frac{1}{2}\left(-\mathbf{a}\cdot \nabla \overline{u}-1{\right)}^{2}-\left(\mathbf{a}\cdot \nabla \overline{u}\right)\overline{u}\right)dx+{\int }_{{\Sigma }_{-}}|\overline{u}{|}^{2}|\mathbf{n}\cdot \mathbf{a}|d\sigma .\end{array}$
The equation follows because each integrand is non-negative.
References

1. G. Auchmuty. Saddle points and existence-uniqueness for evolution equations, Differential Integral Equations, 6 (1993), 1161–1171.
2. G. Auchmuty. Variational principles for operator equations and initial value problems, Nonlinear Analysis, Theory, Methods and Applications Vol. 12, No.5, pp. 531-564 (1988).
3. C. Baiocchi, A. Capelo. Variational and quasivariational inequalities: Applications to free boundary value problems. Wiley, New York (1984).
4. V. Barbu, Optimal control of variational inequalities, Research Notes in Mathematics 100, Pitman (1984).
5. C. Bardos. Problèmes aux limites pour les equations aux dérivées partielles du premier ordre a coefficients réels; Théorèmes d'approximation; Application à l'équation de transport, Ann. scient. Ec. Norm. Sup., ${4}^{e}$  série, t. 3, (1970), p. 185-233.
6. Y. Brenier, Order preserving vibrating strings and applications to Electrodynamics and Magnetohydrodynamics, Preprint (2004)
7. H. Brezis. Operateurs maximaux monotones et semi-groupes de contractions dans les espaces de Hilbert, North Holland, Amsterdam-London, 1973.
8. H. Brezis, I. Ekeland, Un principe variationnel associé à certaines equations paraboliques. Le cas independant du temps, C.R. Acad. Sci. Paris Sér. A 282 (1976), 971–974.
9. H. Brezis, I. Ekeland, Un principe variationnel associé à certaines equations paraboliques. Le cas dependant du temps, C.R. Acad. Sci. Paris Sér. A 282 (1976), 1197–1198.
10. C. Castaing, M. Valadier, Convex Analysis and Measurable Multifunctions, Springer-Verlag, New York, 1977.
11. R. Dautray, J.L. Lions. Mathematical Analysis and Numerical Methods for Science and Technology, Vol 2, Functional and Variational Methods, Springer-Verlag, (1988).
12. L.C. Evans, Partial Differential Equations, Graduate Studies in Mathematics, vol. 19, Amer. Math. Soc., Providence, 1998.
13. N. Ghoussoub, A theory of anti-selfdual Lagrangians: Stationary case, C.R. Acad. Sci., Paris, Ser. I 340 (2005) 245-250
14. N. Ghoussoub, A theory of anti-selfdual Lagrangians: Dynamical case, C.R. Acad. Sci., Paris, Ser. I 340 (2005) p. 325-330
15. N. Ghoussoub, A variational principle for non-linear transport equations, Communications on Pure and Applied Analysis, To appear (2004) 10 p.
16. N. Ghoussoub, Anti-self dual Hamiltonians: Variational resolution for Navier-Stokes equations and other nonlinear evolutions, submitted (2005) 30 pages
17. N. Ghoussoub, $R$  -antiselfduality and the variational resolution of partial differential systems, In preparation (2005)
18. N. Ghoussoub, R. McCann A least action principle for steepest descent in a non-convex landscape, Contemporary Math. Vol 362 (2004) p. 177-187.
19. N. Ghoussoub, L. Tzou. A variational principle for gradient flows, Math. Annalen, Vol 30, 3 (2004) p. 519-549.
20. N. Ghoussoub, L. Tzou. Anti-selfdual Lagrangians II: Unbounded non self-adjoint operators and evolution equations, (2005). Submitted 30 pages
21. N. Ghoussoub and L. Tzou: Iterations of Anti-selfdual Lagrangians, Hamiltonian systems and multiparameter gradient flows, Submitted (2005) 20 pages
22. J. Jost, Riemannian Geometry and Geometric Analysis, Springer Universitytext (2002).
23. J.L. Lions, E. Magenes Non-homogeneous boundary value problems and applications. Vol.3, Springer-Verlag, berlin (1973).
24. R. T. Rockafellar. Existence and duality theorems for convex problems of Bolza. Trans. Amer. Math. Soc. 159 (1971), 1–40.