Each result cited below is universally quantified over the data in its own statement.
The notation of The Wick-Square Problem on the Sobolev Triple of Order Two: Standing Notation , Hilbert Triples: Standing Notation and Background and The Galerkin Problems of the Wick-Square Problem: Standing Notation for Cubes of Modes, Coordinates, the Free Galerkin Quadratic and the Riccati Potential is in force, and elementary facts about real numbers, finite sums, reindexing of finite sums along bijections, absolute values and squares are those carried by The Real Numbers: Standing Notation and Background . The items Viscosity Subsolution, Supersolution and Solution of a Second-Order Equation on a Hilbert Triple , The δ \delta δ -Envelopes u δ − u^-_\delta u δ − and u δ + u^+_\delta u δ + of a Function on an Open Subset of a Hilbert Triple , Second-Order Equation Operator on an Open Subset of a Hilbert Triple and Its δ \delta δ -Shifts , Basic Properties of the δ \delta δ -Envelopes: Semicontinuity, Duality, Closed Superlevel Sets, Bounds and Monotonicity , First- and Second-Order Conditions at a Local Extremum of a C 2 C^2 C 2 Function Penalised by h h h on the Small Space and The Penalty Function h = 1 2 ∣ ⋅ ∣ V 2 h=\tfrac12|\cdot|_V^2 h = 2 1 ∣ ⋅ ∣ V 2 of a Hilbert Triple: Expansion Identities, Closed Sublevel Sets, Density and Local Bounds are stated in the setting of Hilbert Triples: Standing Notation and Background ; their notation (H H H , V V V , A A A , D ( A ) D(A) D ( A ) , h h h , C 2 ( H ) C^{2}(H) C 2 ( H ) , S y m ( H ) \mathrm{Sym}(H) Sym ( H ) , S y m ( V ) \mathrm{Sym}(V) Sym ( V ) , Y ∣ V Y|_{V} Y ∣ V , I V I_{V} I V ) is the notation fixed for the present triple by Hilbert Triples: Standing Notation and Background , and we apply them with U = H U=H U = H , so that V ∩ U = V V\cap U=V V ∩ U = V and W = D ( A ) W=D(A) W = D ( A ) , as recorded in Hilbert Triples: Standing Notation and Background §open-sets .
Since Γ N \Gamma_{N} Γ N is nonempty, m ≥ 1 m\ge1 m ≥ 1 . For i ∈ [ m ] i\in[m] i ∈ [ m ] we write k i = κ N ( i ) k_{i}=\kappa_{N}(i) k i = κ N ( i ) and E i = e k i E_{i}=e_{k_{i}} E i = e k i . Since κ N \kappa_{N} κ N is a bijection from [ m ] [m] [ m ] onto Γ N \Gamma_{N} Γ N (The Galerkin Problems of the Wick-Square Problem: Standing Notation for Cubes of Modes, Coordinates, the Free Galerkin Quadratic and the Riccati Potential §cutoff ), ∑ k ∈ Γ N f ( k ) = ∑ i = 1 m f ( k i ) \sum_{k\in\Gamma_{N}}f(k)=\sum_{i=1}^{m}f(k_{i}) ∑ k ∈ Γ N f ( k ) = ∑ i = 1 m f ( k i ) for every real-valued f f f on Γ N \Gamma_{N} Γ N , and k i ≠ k j k_{i}\neq k_{j} k i = k j for i ≠ j i\neq j i = j . Recall 1 ≤ μ k 1\le\mu_{k} 1 ≤ μ k for every mode k k k (The Wick-Square Problem on the Torus: Standing Notation §modes ).
Step 1 (Coordinate identities). Every unit family e k e_{k} e k lies in H − 1 = D ( A ) ⊆ V ⊆ H H^{-1}=D(A)\subseteq V\subseteq H H − 1 = D ( A ) ⊆ V ⊆ H (The Wick-Square Problem on the Torus: Standing Notation §units , The Sobolev Triple of Order Two for the Wick-Square Problem: the Domain Is the Sobolev Space of Order -1, Unit Families Are Square-Summable, the Gradient Form and the Riccati Drift §domain , Hilbert Triples: Standing Notation and Background §triple ). By The Sobolev Triple of Order Two for the Wick-Square Problem: the Domain Is the Sobolev Space of Order -1, Unit Families Are Square-Summable, the Gradient Form and the Riccati Drift §pairings , for every z ∈ H z\in H z ∈ H and every mode k k k ,
⟨ z , e k ⟩ H = z ( k ) μ k 3 , and, if z ∈ V , ⟨ z , e k ⟩ V = z ( k ) μ k 2 . ( 1 ) \langle z,e_{k}\rangle_{H}=\frac{z(k)}{\mu_{k}^{3}},\qquad\text{and, if }z\in V,\qquad \langle z,e_{k}\rangle_{V}=\frac{z(k)}{\mu_{k}^{2}}.\qquad(1) ⟨ z , e k ⟩ H = μ k 3 z ( k ) , and, if z ∈ V , ⟨ z , e k ⟩ V = μ k 2 z ( k ) . ( 1 )
Taking z = e k z=e_{k} z = e k gives ∣ e k ∣ H 2 = μ k − 3 ≤ 1 |e_{k}|_{H}^{2}=\mu_{k}^{-3}\le1 ∣ e k ∣ H 2 = μ k − 3 ≤ 1 , hence ∣ e k ∣ H ≤ 1 |e_{k}|_{H}\le1 ∣ e k ∣ H ≤ 1 ; and ⟨ e l , e k ⟩ V = e l ( k ) / μ k 2 \langle e_{l},e_{k}\rangle_{V}=e_{l}(k)/\mu_{k}^{2} ⟨ e l , e k ⟩ V = e l ( k ) / μ k 2 , which is μ k − 2 \mu_{k}^{-2} μ k − 2 for l = k l=k l = k and 0 0 0 for l ≠ k l\neq k l = k .
(a) By The Galerkin Problems of the Wick-Square Problem: Standing Notation for Cubes of Modes, Coordinates, the Free Galerkin Quadratic and the Riccati Potential §coordinates , ι N ζ = ∑ i = 1 m ζ i E i \iota_{N}\zeta=\sum_{i=1}^{m}\zeta_{i}E_{i} ι N ζ = ∑ i = 1 m ζ i E i has value ζ i \zeta_{i} ζ i at k i k_{i} k i and 0 0 0 off Γ N \Gamma_{N} Γ N , so p r N ι N ζ = ζ \mathrm{pr}_{N}\iota_{N}\zeta=\zeta pr N ι N ζ = ζ for ζ ∈ R m \zeta\in\mathbb{R}^{m} ζ ∈ R m ; and for every family x x x , Π N x = ι N ( p r N x ) = ∑ i = 1 m x ( k i ) E i \Pi_{N}x=\iota_{N}(\mathrm{pr}_{N}x)=\sum_{i=1}^{m}x(k_{i})E_{i} Π N x = ι N ( pr N x ) = ∑ i = 1 m x ( k i ) E i agrees with x x x at each k i k_{i} k i and vanishes off Γ N \Gamma_{N} Γ N . Hence p r N Π N = p r N \mathrm{pr}_{N}\Pi_{N}=\mathrm{pr}_{N} pr N Π N = pr N and Π N Π N = Π N \Pi_{N}\Pi_{N}=\Pi_{N} Π N Π N = Π N , and p r N \mathrm{pr}_{N} pr N , ι N \iota_{N} ι N , Π N \Pi_{N} Π N are linear, the operations on families being pointwise. For x ∈ H x\in H x ∈ H , Π N x ∈ H − 1 ⊆ V \Pi_{N}x\in H^{-1}\subseteq V Π N x ∈ H − 1 ⊆ V by The Shifted Hamilton-Jacobi-Bellman Equation of the Wick-Square Problem with Noise, Control and Cost Cut Off to a Cube of Modes §data . For a mode k ∉ Γ N k\notin\Gamma_{N} k ∈ / Γ N , e k e_{k} e k vanishes at every k i k_{i} k i , so Π N e k = 0 H \Pi_{N}e_{k}=0_{H} Π N e k = 0 H .
(b) Let ζ , ζ ′ ∈ R m \zeta,\zeta'\in\mathbb{R}^{m} ζ , ζ ′ ∈ R m and let c ∈ H c\in H c ∈ H vanish at every k i k_{i} k i . By Inner Products Against Finite Sums, and Orthonormal Families, in a Real Inner Product Space §combinations and (1),
⟨ ι N ζ , c ⟩ H = ∑ i = 1 m ζ i c ( k i ) μ k i 3 = 0 , ⟨ ι N ζ , ι N ζ ′ ⟩ V = ∑ i = 1 m ζ i ζ i ′ μ k i 2 , \langle\iota_{N}\zeta,c\rangle_{H}=\sum_{i=1}^{m}\zeta_{i}\frac{c(k_{i})}{\mu_{k_{i}}^{3}}=0,\qquad
\langle\iota_{N}\zeta,\iota_{N}\zeta'\rangle_{V}=\sum_{i=1}^{m}\frac{\zeta_{i}\zeta'_{i}}{\mu_{k_{i}}^{2}}, ⟨ ι N ζ , c ⟩ H = i = 1 ∑ m ζ i μ k i 3 c ( k i ) = 0 , ⟨ ι N ζ , ι N ζ ′ ⟩ V = i = 1 ∑ m μ k i 2 ζ i ζ i ′ ,
and, if moreover c ∈ V c\in V c ∈ V , ⟨ ι N ζ , c ⟩ V = ∑ i ζ i c ( k i ) / μ k i 2 = 0 \langle\iota_{N}\zeta,c\rangle_{V}=\sum_{i}\zeta_{i}c(k_{i})/\mu_{k_{i}}^{2}=0 ⟨ ι N ζ , c ⟩ V = ∑ i ζ i c ( k i ) / μ k i 2 = 0 . With The Penalty Function h = 1 2 ∣ ⋅ ∣ V 2 h=\tfrac12|\cdot|_V^2 h = 2 1 ∣ ⋅ ∣ V 2 of a Hilbert Triple: Expansion Identities, Closed Sublevel Sets, Density and Local Bounds §expansion-v , applied to the points ι N ζ \iota_{N}\zeta ι N ζ and ι N ζ + c \iota_{N}\zeta+c ι N ζ + c of V V V , this gives, for c ∈ V c\in V c ∈ V vanishing at every k i k_{i} k i ,
h ( ι N ζ + c ) = h ( ι N ζ ) + h ( c ) , h ( ι N ζ ) = 1 2 ∑ i = 1 m ζ i 2 μ k i 2 . ( 2 ) h(\iota_{N}\zeta+c)=h(\iota_{N}\zeta)+h(c),\qquad h(\iota_{N}\zeta)=\tfrac12\sum_{i=1}^{m}\frac{\zeta_{i}^{2}}{\mu_{k_{i}}^{2}}.\qquad(2) h ( ι N ζ + c ) = h ( ι N ζ ) + h ( c ) , h ( ι N ζ ) = 2 1 i = 1 ∑ m μ k i 2 ζ i 2 . ( 2 )
(c) Let a , a ′ ∈ H a,a'\in H a , a ′ ∈ H . By Inner Products Against Finite Sums, and Orthonormal Families, in a Real Inner Product Space §combinations and (1), ⟨ Π N a , a ′ ⟩ H = ∑ i = 1 m a ( k i ) a ′ ( k i ) / μ k i 3 \langle\Pi_{N}a,a'\rangle_{H}=\sum_{i=1}^{m}a(k_{i})a'(k_{i})/\mu_{k_{i}}^{3} ⟨ Π N a , a ′ ⟩ H = ∑ i = 1 m a ( k i ) a ′ ( k i ) / μ k i 3 , which is symmetric in a , a ′ a,a' a , a ′ ; so ⟨ Π N a , a ′ ⟩ H = ⟨ a , Π N a ′ ⟩ H \langle\Pi_{N}a,a'\rangle_{H}=\langle a,\Pi_{N}a'\rangle_{H} ⟨ Π N a , a ′ ⟩ H = ⟨ a , Π N a ′ ⟩ H . The family a − Π N a a-\Pi_{N}a a − Π N a vanishes at every k i k_{i} k i , so by (b) ⟨ Π N a , a − Π N a ⟩ H = 0 \langle\Pi_{N}a,a-\Pi_{N}a\rangle_{H}=0 ⟨ Π N a , a − Π N a ⟩ H = 0 , and Elementary Identities in a Real Inner Product Space §expansion gives ∣ a ∣ H 2 = ∣ Π N a ∣ H 2 + ∣ a − Π N a ∣ H 2 ≥ ∣ a − Π N a ∣ H 2 |a|_{H}^{2}=|\Pi_{N}a|_{H}^{2}+|a-\Pi_{N}a|_{H}^{2}\ge|a-\Pi_{N}a|_{H}^{2} ∣ a ∣ H 2 = ∣ Π N a ∣ H 2 + ∣ a − Π N a ∣ H 2 ≥ ∣ a − Π N a ∣ H 2 , whence ∣ a − Π N a ∣ H ≤ ∣ a ∣ H |a-\Pi_{N}a|_{H}\le|a|_{H} ∣ a − Π N a ∣ H ≤ ∣ a ∣ H , both sides being nonnegative. For y ∈ V y\in V y ∈ V , (2) with ζ = p r N y \zeta=\mathrm{pr}_{N}y ζ = pr N y and c = y − Π N y c=y-\Pi_{N}y c = y − Π N y gives
h ( y ) = h ( Π N y ) + h ( y − Π N y ) . ( 3 ) h(y)=h(\Pi_{N}y)+h(y-\Pi_{N}y).\qquad(3) h ( y ) = h ( Π N y ) + h ( y − Π N y ) . ( 3 )
(d) For a ∈ H a\in H a ∈ H and a mode k k k , (1) and The Cauchy-Schwarz Inequality in a Real Inner Product Space give ∣ a ( k ) ∣ = μ k 3 ∣ ⟨ a , e k ⟩ H ∣ ≤ μ k 3 ∣ a ∣ H ∣ e k ∣ H ≤ μ k 3 ∣ a ∣ H |a(k)|=\mu_{k}^{3}|\langle a,e_{k}\rangle_{H}|\le\mu_{k}^{3}|a|_{H}|e_{k}|_{H}\le\mu_{k}^{3}|a|_{H} ∣ a ( k ) ∣ = μ k 3 ∣ ⟨ a , e k ⟩ H ∣ ≤ μ k 3 ∣ a ∣ H ∣ e k ∣ H ≤ μ k 3 ∣ a ∣ H . For ζ ∈ R m \zeta\in\mathbb{R}^{m} ζ ∈ R m the family ι N ζ \iota_{N}\zeta ι N ζ has value ζ i \zeta_{i} ζ i at k i k_{i} k i by (a), so Inner Products Against Finite Sums, and Orthonormal Families, in a Real Inner Product Space §combinations and (1) give, as in (b),
∣ ι N ζ ∣ H 2 = ⟨ ι N ζ , ι N ζ ⟩ H = ∑ i = 1 m ζ i ζ i μ k i 3 = ∑ i = 1 m ζ i 2 μ k i 3 ≤ ∑ i = 1 m ζ i 2 = ∥ ζ ∥ 2 , |\iota_{N}\zeta|_{H}^{2}=\langle\iota_{N}\zeta,\iota_{N}\zeta\rangle_{H}=\sum_{i=1}^{m}\zeta_{i}\frac{\zeta_{i}}{\mu_{k_{i}}^{3}}=\sum_{i=1}^{m}\frac{\zeta_{i}^{2}}{\mu_{k_{i}}^{3}}\le\sum_{i=1}^{m}\zeta_{i}^{2}=\lVert\zeta\rVert^{2}, ∣ ι N ζ ∣ H 2 = ⟨ ι N ζ , ι N ζ ⟩ H = i = 1 ∑ m ζ i μ k i 3 ζ i = i = 1 ∑ m μ k i 3 ζ i 2 ≤ i = 1 ∑ m ζ i 2 = ∥ ζ ∥ 2 ,
the inequality termwise because μ k i − 3 ≤ 1 \mu_{k_{i}}^{-3}\le1 μ k i − 3 ≤ 1 , and the last equality by Euclidean Norm on R n \mathbb{R}^n R n . Both ∣ ι N ζ ∣ H |\iota_{N}\zeta|_{H} ∣ ι N ζ ∣ H and ∥ ζ ∥ \lVert\zeta\rVert ∥ ζ ∥ being nonnegative,
∣ ι N ζ ∣ H ≤ ∥ ζ ∥ . ( 4 ) |\iota_{N}\zeta|_{H}\le\lVert\zeta\rVert .\qquad(4) ∣ ι N ζ ∣ H ≤ ∥ ζ ∥ . ( 4 )
Put C N = ∑ i = 1 m μ k i 3 C_{N}=\sum_{i=1}^{m}\mu_{k_{i}}^{3} C N = ∑ i = 1 m μ k i 3 ; its terms are at least 1 1 1 since 1 ≤ μ k 1\le\mu_{k} 1 ≤ μ k , so μ k i 3 ≤ C N \mu_{k_{i}}^{3}\le C_{N} μ k i 3 ≤ C N for each i i i , and 1 ≤ μ k 1 3 ≤ C N 1\le\mu_{k_{1}}^{3}\le C_{N} 1 ≤ μ k 1 3 ≤ C N because 1 ≤ m 1\le m 1 ≤ m . For a ∈ H a\in H a ∈ H , ι N ( p r N a ) = Π N a \iota_{N}(\mathrm{pr}_{N}a)=\Pi_{N}a ι N ( pr N a ) = Π N a , and the identity ∣ a ∣ H 2 = ∣ Π N a ∣ H 2 + ∣ a − Π N a ∣ H 2 |a|_{H}^{2}=|\Pi_{N}a|_{H}^{2}+|a-\Pi_{N}a|_{H}^{2} ∣ a ∣ H 2 = ∣ Π N a ∣ H 2 + ∣ a − Π N a ∣ H 2 of (c) gives ∣ Π N a ∣ H 2 ≤ ∣ a ∣ H 2 |\Pi_{N}a|_{H}^{2}\le|a|_{H}^{2} ∣ Π N a ∣ H 2 ≤ ∣ a ∣ H 2 ; so, by Euclidean Norm on R n \mathbb{R}^n R n and the computation above with ζ = p r N a \zeta=\mathrm{pr}_{N}a ζ = pr N a ,
∥ p r N a ∥ 2 = ∑ i = 1 m μ k i 3 a ( k i ) 2 μ k i 3 ≤ C N ∑ i = 1 m a ( k i ) 2 μ k i 3 = C N ∣ Π N a ∣ H 2 ≤ C N ∣ a ∣ H 2 ≤ C N 2 ∣ a ∣ H 2 , \lVert\mathrm{pr}_{N}a\rVert^{2}=\sum_{i=1}^{m}\mu_{k_{i}}^{3}\,\frac{a(k_{i})^{2}}{\mu_{k_{i}}^{3}}\le C_{N}\sum_{i=1}^{m}\frac{a(k_{i})^{2}}{\mu_{k_{i}}^{3}}=C_{N}|\Pi_{N}a|_{H}^{2}\le C_{N}|a|_{H}^{2}\le C_{N}^{2}|a|_{H}^{2}, ∥ pr N a ∥ 2 = i = 1 ∑ m μ k i 3 μ k i 3 a ( k i ) 2 ≤ C N i = 1 ∑ m μ k i 3 a ( k i ) 2 = C N ∣ Π N a ∣ H 2 ≤ C N ∣ a ∣ H 2 ≤ C N 2 ∣ a ∣ H 2 ,
whence, both sides being nonnegative,
∥ p r N a ∥ ≤ C N ∣ a ∣ H ( a ∈ H ) . ( 4 ′ ) \lVert\mathrm{pr}_{N}a\rVert\le C_{N}|a|_{H}\qquad(a\in H).\qquad(4') ∥ pr N a ∥ ≤ C N ∣ a ∣ H ( a ∈ H ) . ( 4 ′ )
Step 2 (Clause 1, and semicontinuity of w w w ). For a , a ′ ∈ H a,a'\in H a , a ′ ∈ H , p r N a − p r N a ′ = p r N ( a − a ′ ) \mathrm{pr}_{N}a-\mathrm{pr}_{N}a'=\mathrm{pr}_{N}(a-a') pr N a − pr N a ′ = pr N ( a − a ′ ) by linearity of p r N \mathrm{pr}_{N} pr N (Step 1(a)), so (4') gives ∥ p r N a − p r N a ′ ∥ ≤ C N ∣ a − a ′ ∣ H \lVert\mathrm{pr}_{N}a-\mathrm{pr}_{N}a'\rVert\le C_{N}|a-a'|_{H} ∥ pr N a − pr N a ′ ∥ ≤ C N ∣ a − a ′ ∣ H . Given a real ε > 0 \varepsilon>0 ε > 0 , every a ′ ∈ H a'\in H a ′ ∈ H with ∣ a − a ′ ∣ H < ε / C N |a-a'|_{H}<\varepsilon/C_{N} ∣ a − a ′ ∣ H < ε / C N thus satisfies d E ( p r N a , p r N a ′ ) < ε d_{E}(\mathrm{pr}_{N}a,\mathrm{pr}_{N}a')<\varepsilon d E ( pr N a , pr N a ′ ) < ε , so p r N : H → R m \mathrm{pr}_{N}:H\to\mathbb{R}^{m} pr N : H → R m is continuous on H H H relative to H H H in the sense of Continuous Map Between Metric Spaces . By Semicontinuity and Continuity Under Composition with a Continuous Map , applied with the metric spaces ( R m , d E ) (\mathbb{R}^{m},d_{E}) ( R m , d E ) and ( H , d H ) (H,d_{H}) ( H , d H ) , with its map taken to be p r N \mathrm{pr}_{N} pr N on H H H and its function taken to be Z Z Z on R m \mathbb{R}^{m} R m : if Z Z Z is continuous on ( R m , d E ) (\mathbb{R}^{m},d_{E}) ( R m , d E ) then w = Z ∘ p r N w=Z\circ\mathrm{pr}_{N} w = Z ∘ pr N is continuous on H H H (claim 3), which is continuity on H H H in the sense of Real Hilbert Spaces: Standing Notation and Background §topology ; this proves clause 1 . Likewise, by claims 1 and 2 there, w w w is upper semicontinuous on H H H when Z Z Z is upper semicontinuous on R m \mathbb{R}^{m} R m , and lower semicontinuous on H H H when Z Z Z is lower semicontinuous on R m \mathbb{R}^{m} R m .
Step 3 (Composition with an affine contraction). For a ∈ H a\in H a ∈ H put L a = a − Π N a La=a-\Pi_{N}a L a = a − Π N a . By Step 1(a),(c), L : H → H L:H\to H L : H → H is linear, ∣ L a ∣ H ≤ ∣ a ∣ H |La|_{H}\le|a|_{H} ∣ L a ∣ H ≤ ∣ a ∣ H , and ⟨ L a , a ′ ⟩ H = ⟨ a , a ′ ⟩ H − ⟨ Π N a , a ′ ⟩ H = ⟨ a , L a ′ ⟩ H \langle La,a'\rangle_{H}=\langle a,a'\rangle_{H}-\langle\Pi_{N}a,a'\rangle_{H}=\langle a,La'\rangle_{H} ⟨ L a , a ′ ⟩ H = ⟨ a , a ′ ⟩ H − ⟨ Π N a , a ′ ⟩ H = ⟨ a , L a ′ ⟩ H for a , a ′ ∈ H a,a'\in H a , a ′ ∈ H . Let c ∈ H c\in H c ∈ H , let χ ∈ C 2 ( H ) \chi\in C^{2}(H) χ ∈ C 2 ( H ) , and define χ T : H → R \chi^{T}:H\to\mathbb{R} χ T : H → R by χ T ( y ) = χ ( c + L y ) \chi^{T}(y)=\chi(c+Ly) χ T ( y ) = χ ( c + L y ) . We claim that χ T ∈ C 2 ( H ) \chi^{T}\in C^{2}(H) χ T ∈ C 2 ( H ) , with
D χ T ( y ) = L D χ ( c + L y ) , D 2 χ T ( y ) = B y , B y ( a , a ′ ) = D 2 χ ( c + L y ) ( L a , L a ′ ) . D\chi^{T}(y)=L\,D\chi(c+Ly),\qquad D^{2}\chi^{T}(y)=B_{y},\quad B_{y}(a,a')=D^{2}\chi(c+Ly)(La,La'). D χ T ( y ) = L Dχ ( c + L y ) , D 2 χ T ( y ) = B y , B y ( a , a ′ ) = D 2 χ ( c + L y ) ( L a , L a ′ ) .
Fix y ∈ H y\in H y ∈ H and write x = c + L y x=c+Ly x = c + L y , so c + L ( y + v ) = x + L v c+L(y+v)=x+Lv c + L ( y + v ) = x + Lv for v ∈ H v\in H v ∈ H . Differentiability, in the sense of Fréchet Differentiability and the Gradient on an Open Subset of a Real Inner Product Space §differentiable : given a real ε > 0 \varepsilon>0 ε > 0 , let η > 0 \eta>0 η > 0 be a radius for χ \chi χ at x x x and ε \varepsilon ε ; for v ∈ H v\in H v ∈ H with ∣ v ∣ H < η |v|_{H}<\eta ∣ v ∣ H < η one has ∣ L v ∣ H ≤ ∣ v ∣ H < η |Lv|_{H}\le|v|_{H}<\eta ∣ Lv ∣ H ≤ ∣ v ∣ H < η , hence
∣ χ T ( y + v ) − χ T ( y ) − ⟨ L D χ ( x ) , v ⟩ H ∣ = ∣ χ ( x + L v ) − χ ( x ) − ⟨ D χ ( x ) , L v ⟩ H ∣ ≤ ε ∣ L v ∣ H ≤ ε ∣ v ∣ H . \bigl|\chi^{T}(y+v)-\chi^{T}(y)-\langle LD\chi(x),v\rangle_{H}\bigr|=\bigl|\chi(x+Lv)-\chi(x)-\langle D\chi(x),Lv\rangle_{H}\bigr|\le\varepsilon|Lv|_{H}\le\varepsilon|v|_{H}. χ T ( y + v ) − χ T ( y ) − ⟨ L Dχ ( x ) , v ⟩ H = χ ( x + Lv ) − χ ( x ) − ⟨ Dχ ( x ) , Lv ⟩ H ≤ ε ∣ Lv ∣ H ≤ ε ∣ v ∣ H .
So χ T \chi^{T} χ T is differentiable at y y y with gradient L D χ ( x ) LD\chi(x) L Dχ ( x ) , and by uniqueness of the gradient (Fréchet Differentiability and the Gradient on an Open Subset of a Real Inner Product Space §gradient ) D χ T ( y ) = L D χ ( x ) D\chi^{T}(y)=LD\chi(x) D χ T ( y ) = L Dχ ( x ) ; as y y y was arbitrary, χ T \chi^{T} χ T is differentiable on H H H and D χ T ( y ′ ) = L D χ ( c + L y ′ ) D\chi^{T}(y')=LD\chi(c+Ly') D χ T ( y ′ ) = L Dχ ( c + L y ′ ) for every y ′ ∈ H y'\in H y ′ ∈ H . Continuity of the gradient map: ∣ L D χ ( c + L y ′ ) − L D χ ( c + L y ) ∣ H ≤ ∣ D χ ( c + L y ′ ) − D χ ( c + L y ) ∣ H |LD\chi(c+Ly')-LD\chi(c+Ly)|_{H}\le|D\chi(c+Ly')-D\chi(c+Ly)|_{H} ∣ L Dχ ( c + L y ′ ) − L Dχ ( c + L y ) ∣ H ≤ ∣ Dχ ( c + L y ′ ) − Dχ ( c + L y ) ∣ H and ∣ ( c + L y ′ ) − ( c + L y ) ∣ H ≤ ∣ y ′ − y ∣ H |(c+Ly')-(c+Ly)|_{H}\le|y'-y|_{H} ∣ ( c + L y ′ ) − ( c + L y ) ∣ H ≤ ∣ y ′ − y ∣ H , so continuity of D χ D\chi Dχ at x x x (part of The Classes C 1 C^1 C 1 and C 2 C^2 C 2 on an Open Subset of a Real Inner Product Space §c1 ) yields continuity of D χ T D\chi^{T} D χ T at y y y with the same radius. The form B y B_{y} B y : symmetry, additivity and homogeneity follow from those of D 2 χ ( x ) D^{2}\chi(x) D 2 χ ( x ) and linearity of L L L , and ∣ B y ( a , a ′ ) ∣ ≤ ∥ D 2 χ ( x ) ∥ ∣ L a ∣ H ∣ L a ′ ∣ H ≤ ∥ D 2 χ ( x ) ∥ ∣ a ∣ H ∣ a ′ ∣ H |B_{y}(a,a')|\le\lVert D^{2}\chi(x)\rVert|La|_{H}|La'|_{H}\le\lVert D^{2}\chi(x)\rVert|a|_{H}|a'|_{H} ∣ B y ( a , a ′ ) ∣ ≤ ∥ D 2 χ ( x )∥ ∣ L a ∣ H ∣ L a ′ ∣ H ≤ ∥ D 2 χ ( x )∥ ∣ a ∣ H ∣ a ′ ∣ H by Elementary Properties of Bounded Symmetric Bilinear Forms: Norm, Quadratic Form, Order and Continuity §bound ; so B y ∈ S y m ( H ) B_{y}\in\mathrm{Sym}(H) B y ∈ Sym ( H ) by Bounded Symmetric Bilinear Forms on a Real Inner Product Space: Norm, Order, Identity Form and Restriction §form . Second derivative, in the sense of The Second Derivative and the Hessian on an Open Subset of a Real Inner Product Space §expansion with ρ = 1 \rho=1 ρ = 1 (every ball lies in H H H and χ T \chi^{T} χ T is differentiable on H H H ): given ε > 0 \varepsilon>0 ε > 0 , let δ ′ ′ ≤ 1 \delta''\le1 δ ′′ ≤ 1 be a radius witnessing that D 2 χ ( x ) D^{2}\chi(x) D 2 χ ( x ) is a second derivative of χ \chi χ at x x x (with ρ = 1 \rho=1 ρ = 1 , admissible since the condition does not depend on ρ \rho ρ ). For v , a ∈ H v,a\in H v , a ∈ H with ∣ v ∣ H < δ ′ ′ |v|_{H}<\delta'' ∣ v ∣ H < δ ′′ , ∣ L v ∣ H < δ ′ ′ |Lv|_{H}<\delta'' ∣ Lv ∣ H < δ ′′ , and by self-adjointness of L L L ,
∣ ⟨ D χ T ( y + v ) − D χ T ( y ) , a ⟩ H − B y ( v , a ) ∣ = ∣ ⟨ D χ ( x + L v ) − D χ ( x ) , L a ⟩ H − D 2 χ ( x ) ( L v , L a ) ∣ ≤ ε ∣ L v ∣ H ∣ L a ∣ H ≤ ε ∣ v ∣ H ∣ a ∣ H . \bigl|\langle D\chi^{T}(y+v)-D\chi^{T}(y),a\rangle_{H}-B_{y}(v,a)\bigr|=\bigl|\langle D\chi(x+Lv)-D\chi(x),La\rangle_{H}-D^{2}\chi(x)(Lv,La)\bigr|\le\varepsilon|Lv|_{H}|La|_{H}\le\varepsilon|v|_{H}|a|_{H}. ⟨ D χ T ( y + v ) − D χ T ( y ) , a ⟩ H − B y ( v , a ) = ⟨ Dχ ( x + Lv ) − Dχ ( x ) , L a ⟩ H − D 2 χ ( x ) ( Lv , L a ) ≤ ε ∣ Lv ∣ H ∣ L a ∣ H ≤ ε ∣ v ∣ H ∣ a ∣ H .
So B y B_{y} B y is a second derivative of χ T \chi^{T} χ T at y y y , and by uniqueness of the Hessian (The Second Derivative and the Hessian on an Open Subset of a Real Inner Product Space §hessian ) D 2 χ T ( y ) = B y D^{2}\chi^{T}(y)=B_{y} D 2 χ T ( y ) = B y . Continuity of the Hessian map: for y ′ ∈ H y'\in H y ′ ∈ H and a , a ′ ∈ H a,a'\in H a , a ′ ∈ H , ∣ ( B y ′ − B y ) ( a , a ′ ) ∣ = ∣ ( D 2 χ ( c + L y ′ ) − D 2 χ ( x ) ) ( L a , L a ′ ) ∣ ≤ ∥ D 2 χ ( c + L y ′ ) − D 2 χ ( x ) ∥ ∣ a ∣ H ∣ a ′ ∣ H |(B_{y'}-B_{y})(a,a')|=|(D^{2}\chi(c+Ly')-D^{2}\chi(x))(La,La')|\le\lVert D^{2}\chi(c+Ly')-D^{2}\chi(x)\rVert\,|a|_{H}|a'|_{H} ∣ ( B y ′ − B y ) ( a , a ′ ) ∣ = ∣ ( D 2 χ ( c + L y ′ ) − D 2 χ ( x )) ( L a , L a ′ ) ∣ ≤ ∥ D 2 χ ( c + L y ′ ) − D 2 χ ( x )∥ ∣ a ∣ H ∣ a ′ ∣ H by Elementary Properties of Bounded Symmetric Bilinear Forms: Norm, Quadratic Form, Order and Continuity §bound , so by the second part of that claim d S y m ( B y ′ , B y ) ≤ d S y m ( D 2 χ ( c + L y ′ ) , D 2 χ ( x ) ) d_{\mathrm{Sym}}(B_{y'},B_{y})\le d_{\mathrm{Sym}}(D^{2}\chi(c+Ly'),D^{2}\chi(x)) d Sym ( B y ′ , B y ) ≤ d Sym ( D 2 χ ( c + L y ′ ) , D 2 χ ( x )) , and continuity of D 2 χ D^{2}\chi D 2 χ at x x x (part of The Classes C 1 C^1 C 1 and C 2 C^2 C 2 on an Open Subset of a Real Inner Product Space §c2 ) gives continuity of y ↦ B y y\mapsto B_{y} y ↦ B y . This proves the claim.
Step 4 (Restriction of a C 2 C^{2} C 2 function to a coordinate slice). Let y 0 ∈ H y_{0}\in H y 0 ∈ H and χ ∈ C 2 ( H ) \chi\in C^{2}(H) χ ∈ C 2 ( H ) , put S ( z ) = ι N z + y 0 S(z)=\iota_{N}z+y_{0} S ( z ) = ι N z + y 0 for z ∈ R m z\in\mathbb{R}^{m} z ∈ R m and χ 0 ( z ) = χ ( S ( z ) ) \chi_{0}(z)=\chi(S(z)) χ 0 ( z ) = χ ( S ( z )) . By (4) and linearity of ι N \iota_{N} ι N , ∣ S ( z ′ ) − S ( z ) ∣ H ≤ ∥ z ′ − z ∥ |S(z')-S(z)|_{H}\le\lVert z'-z\rVert ∣ S ( z ′ ) − S ( z ) ∣ H ≤ ∥ z ′ − z ∥ . We claim that χ 0 \chi_{0} χ 0 is of class C 2 C^{2} C 2 on R m \mathbb{R}^{m} R m (C^k Maps on a Euclidean Open Set ), with
∂ i χ 0 ( z ) = ⟨ D χ ( S ( z ) ) , E i ⟩ H , ∂ l ∂ i χ 0 ( z ) = D 2 χ ( S ( z ) ) ( E l , E i ) ( i , l ∈ [ m ] ) . \partial_{i}\chi_{0}(z)=\langle D\chi(S(z)),E_{i}\rangle_{H},\qquad \partial_{l}\partial_{i}\chi_{0}(z)=D^{2}\chi(S(z))(E_{l},E_{i})\qquad(i,l\in[m]). ∂ i χ 0 ( z ) = ⟨ Dχ ( S ( z )) , E i ⟩ H , ∂ l ∂ i χ 0 ( z ) = D 2 χ ( S ( z )) ( E l , E i ) ( i , l ∈ [ m ]) .
Continuity: if Θ : H → R \Theta:H\to\mathbb{R} Θ : H → R is continuous at S ( z ) S(z) S ( z ) with radius η \eta η for a given ε \varepsilon ε , then Θ ∘ S \Theta\circ S Θ ∘ S is continuous at z z z with the same radius η \eta η . This applies to χ \chi χ (continuous by Basic Properties of Differentiability on an Open Subset of a Real Inner Product Space §continuous ), to x ↦ ⟨ D χ ( x ) , E i ⟩ H x\mapsto\langle D\chi(x),E_{i}\rangle_{H} x ↦ ⟨ Dχ ( x ) , E i ⟩ H (as ∣ ⟨ D χ ( x ′ ) − D χ ( x ) , E i ⟩ H ∣ ≤ ∣ D χ ( x ′ ) − D χ ( x ) ∣ H |\langle D\chi(x')-D\chi(x),E_{i}\rangle_{H}|\le|D\chi(x')-D\chi(x)|_{H} ∣ ⟨ Dχ ( x ′ ) − Dχ ( x ) , E i ⟩ H ∣ ≤ ∣ Dχ ( x ′ ) − Dχ ( x ) ∣ H by The Cauchy-Schwarz Inequality in a Real Inner Product Space and ∣ E i ∣ H ≤ 1 |E_{i}|_{H}\le1 ∣ E i ∣ H ≤ 1 , and D χ D\chi Dχ is continuous) and to x ↦ D 2 χ ( x ) ( E l , E i ) x\mapsto D^{2}\chi(x)(E_{l},E_{i}) x ↦ D 2 χ ( x ) ( E l , E i ) (as ∣ ( D 2 χ ( x ′ ) − D 2 χ ( x ) ) ( E l , E i ) ∣ ≤ ∥ D 2 χ ( x ′ ) − D 2 χ ( x ) ∥ |(D^{2}\chi(x')-D^{2}\chi(x))(E_{l},E_{i})|\le\lVert D^{2}\chi(x')-D^{2}\chi(x)\rVert ∣ ( D 2 χ ( x ′ ) − D 2 χ ( x )) ( E l , E i ) ∣ ≤ ∥ D 2 χ ( x ′ ) − D 2 χ ( x )∥ by Elementary Properties of Bounded Symmetric Bilinear Forms: Norm, Quadratic Form, Order and Continuity §continuity , and D 2 χ D^{2}\chi D 2 χ is continuous). First partials, in the sense of Partial Derivative on a Euclidean Open Set : moving the i i i th coordinate of z z z by t t t moves S ( z ) S(z) S ( z ) to S ( z ) + t E i S(z)+tE_{i} S ( z ) + t E i . Given ε > 0 \varepsilon>0 ε > 0 , let η > 0 \eta>0 η > 0 be a radius for differentiability of χ \chi χ at S ( z ) S(z) S ( z ) with gradient D χ ( S ( z ) ) D\chi(S(z)) Dχ ( S ( z )) and ε / 2 \varepsilon/2 ε /2 , the gradient D χ ( S ( z ) ) D\chi(S(z)) Dχ ( S ( z )) being by definition the unique vector with this property (Fréchet Differentiability and the Gradient on an Open Subset of a Real Inner Product Space §gradient ); for 0 < ∣ t ∣ < η 0<|t|<\eta 0 < ∣ t ∣ < η , ∣ t E i ∣ H ≤ ∣ t ∣ < η |tE_{i}|_{H}\le|t|<\eta ∣ t E i ∣ H ≤ ∣ t ∣ < η , so ∣ χ ( S ( z ) + t E i ) − χ ( S ( z ) ) − t ⟨ D χ ( S ( z ) ) , E i ⟩ H ∣ ≤ ε 2 ∣ t ∣ |\chi(S(z)+tE_{i})-\chi(S(z))-t\langle D\chi(S(z)),E_{i}\rangle_{H}|\le\tfrac{\varepsilon}{2}|t| ∣ χ ( S ( z ) + t E i ) − χ ( S ( z )) − t ⟨ Dχ ( S ( z )) , E i ⟩ H ∣ ≤ 2 ε ∣ t ∣ , and dividing by ∣ t ∣ |t| ∣ t ∣ the difference quotient of χ 0 \chi_{0} χ 0 is within ε 2 < ε \tfrac{\varepsilon}{2}<\varepsilon 2 ε < ε of ⟨ D χ ( S ( z ) ) , E i ⟩ H \langle D\chi(S(z)),E_{i}\rangle_{H} ⟨ Dχ ( S ( z )) , E i ⟩ H . Second partials: with ω i ( z ) = ⟨ D χ ( S ( z ) ) , E i ⟩ H \omega_{i}(z)=\langle D\chi(S(z)),E_{i}\rangle_{H} ω i ( z ) = ⟨ Dχ ( S ( z )) , E i ⟩ H one has ω i ( z ′ ) − ω i ( z ) = ⟨ D χ ( S ( z ) + t E l ) − D χ ( S ( z ) ) , E i ⟩ H \omega_{i}(z')-\omega_{i}(z)=\langle D\chi(S(z)+tE_{l})-D\chi(S(z)),E_{i}\rangle_{H} ω i ( z ′ ) − ω i ( z ) = ⟨ Dχ ( S ( z ) + t E l ) − Dχ ( S ( z )) , E i ⟩ H when z ′ z' z ′ is z z z with l l l th coordinate moved by t t t ; given ε > 0 \varepsilon>0 ε > 0 , a radius δ ′ ′ \delta'' δ ′′ witnessing The Second Derivative and the Hessian on an Open Subset of a Real Inner Product Space §expansion for χ \chi χ at S ( z ) S(z) S ( z ) with ε / 2 \varepsilon/2 ε /2 and the form D 2 χ ( S ( z ) ) D^{2}\chi(S(z)) D 2 χ ( S ( z )) , which is by definition the unique second derivative of χ \chi χ at S ( z ) S(z) S ( z ) (The Second Derivative and the Hessian on an Open Subset of a Real Inner Product Space §hessian ), used with its increment t E l tE_{l} t E l and direction E i E_{i} E i , gives for 0 < ∣ t ∣ < δ ′ ′ 0<|t|<\delta'' 0 < ∣ t ∣ < δ ′′ the bound ∣ ω i ( z ′ ) − ω i ( z ) − t D 2 χ ( S ( z ) ) ( E l , E i ) ∣ ≤ ε 2 ∣ t ∣ |\omega_{i}(z')-\omega_{i}(z)-t\,D^{2}\chi(S(z))(E_{l},E_{i})|\le\tfrac{\varepsilon}{2}|t| ∣ ω i ( z ′ ) − ω i ( z ) − t D 2 χ ( S ( z )) ( E l , E i ) ∣ ≤ 2 ε ∣ t ∣ , so the difference quotient of ω i \omega_{i} ω i is within ε 2 < ε \tfrac{\varepsilon}{2}<\varepsilon 2 ε < ε of D 2 χ ( S ( z ) ) ( E l , E i ) D^{2}\chi(S(z))(E_{l},E_{i}) D 2 χ ( S ( z )) ( E l , E i ) . Thus χ 0 \chi_{0} χ 0 is continuous with continuous first partials, i.e. of class C 1 C^{1} C 1 , and each ∂ i χ 0 = ω i \partial_{i}\chi_{0}=\omega_{i} ∂ i χ 0 = ω i is continuous with continuous first partials, i.e. of class C 1 C^{1} C 1 ; by clause 2 of C^k Maps on a Euclidean Open Set (with k = 1 k=1 k = 1 ) and its scalar convention, χ 0 \chi_{0} χ 0 is of class C 2 C^{2} C 2 on R m \mathbb{R}^{m} R m , and the formulas hold.
Step 5 (Set-up for clauses 2 and 3). Let σ ∈ { 1 , − 1 } \sigma\in\{1,-1\} σ ∈ { 1 , − 1 } : in case σ = 1 \sigma=1 σ = 1 assume that Z Z Z is a viscosity subsolution of F N ♯ \mathcal{F}^{\sharp}_{N} F N ♯ on R m \mathbb{R}^{m} R m , and in case σ = − 1 \sigma=-1 σ = − 1 that it is a viscosity supersolution, both in the sense of The Galerkin Hamilton-Jacobi-Bellman Equations of the Wick-Square Problem in the Coordinates of a Cube of Modes §equation , that is of Viscosity Subsolution and Supersolution of a Second-Order Equation for the open set R m \mathbb{R}^{m} R m . Then Z Z Z is upper semicontinuous on R m \mathbb{R}^{m} R m if σ = 1 \sigma=1 σ = 1 and lower semicontinuous if σ = − 1 \sigma=-1 σ = − 1 , so by Step 2 w w w is upper, respectively lower, semicontinuous on H H H . By Basic Properties of the δ \delta δ -Envelopes: Semicontinuity, Duality, Closed Superlevel Sets, Bounds and Monotonicity §semicontinuous-case , w w w is bounded above (if σ = 1 \sigma=1 σ = 1 ), respectively below (if σ = − 1 \sigma=-1 σ = − 1 ), near each point of H H H ; this one-sided local bound comes from the semicontinuity of w w w alone, and no bound on Z Z Z is assumed or used anywhere in this proof. So the standing hypothesis of Viscosity Subsolution, Supersolution and Solution of a Second-Order Equation on a Hilbert Triple §subsolution , respectively Viscosity Subsolution, Supersolution and Solution of a Second-Order Equation on a Hilbert Triple §supersolution , holds; and for every real δ > 0 \delta>0 δ > 0 the envelope w δ σ w^{\sigma}_{\delta} w δ σ , meaning w δ − w^{-}_{\delta} w δ − if σ = 1 \sigma=1 σ = 1 and w δ + w^{+}_{\delta} w δ + if σ = − 1 \sigma=-1 σ = − 1 (The δ \delta δ -Envelopes u δ − u^-_\delta u δ − and u δ + u^+_\delta u δ + of a Function on an Open Subset of a Hilbert Triple §minus , The δ \delta δ -Envelopes u δ − u^-_\delta u δ − and u δ + u^+_\delta u δ + of a Function on an Open Subset of a Hilbert Triple §plus ), satisfies
w δ σ ( x ) = w ( x ) − σ δ h ( x ) ( x ∈ V ) . w^{\sigma}_{\delta}(x)=w(x)-\sigma\delta h(x)\qquad(x\in V). w δ σ ( x ) = w ( x ) − σ δ h ( x ) ( x ∈ V ) .
The operator F N ♯ F^{\sharp}_{N} F N ♯ is a second-order equation operator on H H H relative to ( H , V , A ) (H,V,A) ( H , V , A ) by The Shifted Hamilton-Jacobi-Bellman Equation of the Wick-Square Problem with Noise, Control and Cost Cut Off to a Cube of Modes §operator . Write F σ F^{\sigma} F σ for its shift ( F N ♯ ) δ − (F^{\sharp}_{N})^{-}_{\delta} ( F N ♯ ) δ − if σ = 1 \sigma=1 σ = 1 and ( F N ♯ ) δ + (F^{\sharp}_{N})^{+}_{\delta} ( F N ♯ ) δ + if σ = − 1 \sigma=-1 σ = − 1 ; by Second-Order Equation Operator on an Open Subset of a Hilbert Triple and Its δ \delta δ -Shifts §shifted , in both cases
F σ ( x , r , p , Y ) = F N ♯ ( x , r + σ δ h ( x ) , p + σ δ A x , Y ∣ V + σ δ I V ) F^{\sigma}(x,r,p,Y)=F^{\sharp}_{N}\bigl(x,\ r+\sigma\delta h(x),\ p+\sigma\delta Ax,\ Y|_{V}+\sigma\delta I_{V}\bigr) F σ ( x , r , p , Y ) = F N ♯ ( x , r + σ δ h ( x ) , p + σ δ A x , Y ∣ V + σ δ I V )
for x ∈ D ( A ) x\in D(A) x ∈ D ( A ) , r ∈ R r\in\mathbb{R} r ∈ R , p ∈ H p\in H p ∈ H and Y ∈ S y m ( H ) Y\in\mathrm{Sym}(H) Y ∈ Sym ( H ) , where Y ∣ V + ( − δ ) I V Y|_{V}+(-\delta)I_{V} Y ∣ V + ( − δ ) I V is the form Y ∣ V − δ I V Y|_{V}-\delta I_{V} Y ∣ V − δ I V .
Fix a real δ > 0 \delta>0 δ > 0 , φ ∈ C 2 ( H ) \varphi\in C^{2}(H) φ ∈ C 2 ( H ) , a point x ^ ∈ V \hat{x}\in V x ^ ∈ V at which w δ σ − φ w^{\sigma}_{\delta}-\varphi w δ σ − φ has a local maximum (if σ = 1 \sigma=1 σ = 1 ), respectively a local minimum (if σ = − 1 \sigma=-1 σ = − 1 ), relative to V V V , and a real ε > 0 \varepsilon>0 ε > 0 . By Local Maximum of a Function Relative to a Subset of a Metric Space and Local Minimum of a Function Relative to a Subset of a Metric Space , with the metric d H d_{H} d H (Hilbert Triples: Standing Notation and Background §open-sets ), there is a real r > 0 r>0 r > 0 such that every x ∈ V x\in V x ∈ V with ∣ x − x ^ ∣ H < r |x-\hat{x}|_{H}<r ∣ x − x ^ ∣ H < r satisfies
σ [ ( w ( x ) − σ δ h ( x ) − φ ( x ) ) − ( w ( x ^ ) − σ δ h ( x ^ ) − φ ( x ^ ) ) ] ≤ 0. ( 5 ) \sigma\Bigl[\bigl(w(x)-\sigma\delta h(x)-\varphi(x)\bigr)-\bigl(w(\hat{x})-\sigma\delta h(\hat{x})-\varphi(\hat{x})\bigr)\Bigr]\le0.\qquad(5) σ [ ( w ( x ) − σ δ h ( x ) − φ ( x ) ) − ( w ( x ^ ) − σ δ h ( x ^ ) − φ ( x ^ ) ) ] ≤ 0. ( 5 )
Put z ^ = p r N x ^ ∈ R m \hat{z}=\mathrm{pr}_{N}\hat{x}\in\mathbb{R}^{m} z ^ = pr N x ^ ∈ R m and y 0 = x ^ − Π N x ^ = L x ^ ∈ V y_{0}=\hat{x}-\Pi_{N}\hat{x}=L\hat{x}\in V y 0 = x ^ − Π N x ^ = L x ^ ∈ V ; then x ^ = ι N z ^ + y 0 \hat{x}=\iota_{N}\hat{z}+y_{0} x ^ = ι N z ^ + y 0 , y 0 y_{0} y 0 vanishes at every k i k_{i} k i , and z ^ i = x ^ ( k i ) \hat{z}_{i}=\hat{x}(k_{i}) z ^ i = x ^ ( k i ) .
Step 6 (The point x ^ \hat{x} x ^ lies in D ( A ) D(A) D ( A ) , and the shifted gradient has no high modes). Define T : H → H T:H\to H T : H → H by T y = Π N x ^ + L y Ty=\Pi_{N}\hat{x}+Ly T y = Π N x ^ + L y . For y ∈ V y\in V y ∈ V : T y ∈ V Ty\in V T y ∈ V ; p r N T y = p r N x ^ + p r N y − p r N y = z ^ \mathrm{pr}_{N}Ty=\mathrm{pr}_{N}\hat{x}+\mathrm{pr}_{N}y-\mathrm{pr}_{N}y=\hat{z} pr N T y = pr N x ^ + pr N y − pr N y = z ^ by Step 1(a), so w ( T y ) = Z ( z ^ ) = w ( x ^ ) w(Ty)=Z(\hat{z})=w(\hat{x}) w ( T y ) = Z ( z ^ ) = w ( x ^ ) ; T y = ι N z ^ + L y Ty=\iota_{N}\hat{z}+Ly T y = ι N z ^ + L y with L y ∈ V Ly\in V L y ∈ V vanishing at every k i k_{i} k i , so (2) gives h ( T y ) = h ( Π N x ^ ) + h ( L y ) h(Ty)=h(\Pi_{N}\hat{x})+h(Ly) h ( T y ) = h ( Π N x ^ ) + h ( L y ) ; T x ^ = x ^ T\hat{x}=\hat{x} T x ^ = x ^ ; and T y − x ^ = L ( y − x ^ ) Ty-\hat{x}=L(y-\hat{x}) T y − x ^ = L ( y − x ^ ) , so ∣ T y − x ^ ∣ H ≤ ∣ y − x ^ ∣ H |Ty-\hat{x}|_{H}\le|y-\hat{x}|_{H} ∣ T y − x ^ ∣ H ≤ ∣ y − x ^ ∣ H .
Let b N b_{N} b N be the map assigning to y , y ′ ∈ H y,y'\in H y , y ′ ∈ H the number b N ( y , y ′ ) = δ ∑ i = 1 m y ( k i ) y ′ ( k i ) / μ k i 2 b_{N}(y,y')=\delta\sum_{i=1}^{m}y(k_{i})y'(k_{i})/\mu_{k_{i}}^{2} b N ( y , y ′ ) = δ ∑ i = 1 m y ( k i ) y ′ ( k i ) / μ k i 2 . It is symmetric, additive and homogeneous in each argument (the operations being pointwise), and by Step 1(d) ∣ b N ( y , y ′ ) ∣ ≤ ( δ ∑ i = 1 m μ k i 4 ) ∣ y ∣ H ∣ y ′ ∣ H |b_{N}(y,y')|\le\bigl(\delta\sum_{i=1}^{m}\mu_{k_{i}}^{4}\bigr)|y|_{H}|y'|_{H} ∣ b N ( y , y ′ ) ∣ ≤ ( δ ∑ i = 1 m μ k i 4 ) ∣ y ∣ H ∣ y ′ ∣ H ; so b N ∈ S y m ( H ) b_{N}\in\mathrm{Sym}(H) b N ∈ Sym ( H ) . Let Q ( y ) = 1 2 b N ( y , y ) Q(y)=\tfrac12 b_{N}(y,y) Q ( y ) = 2 1 b N ( y , y ) . By Affine and Quadratic Functions on a Real Hilbert Space are of Class C 2 C^2 C 2 §form , Q ∈ C 2 ( H ) Q\in C^{2}(H) Q ∈ C 2 ( H ) with D Q ( y ) = T b N y DQ(y)=T_{b_{N}}y D Q ( y ) = T b N y , where ⟨ T b N y , a ⟩ H = b N ( y , a ) \langle T_{b_{N}}y,a\rangle_{H}=b_{N}(y,a) ⟨ T b N y , a ⟩ H = b N ( y , a ) for a ∈ H a\in H a ∈ H (The Bounded Symmetric Operator Represented by a Bounded Symmetric Bilinear Form on a Real Hilbert Space §representation ). By (2) with ζ = p r N y \zeta=\mathrm{pr}_{N}y ζ = pr N y , Q ( y ) = δ h ( Π N y ) Q(y)=\delta h(\Pi_{N}y) Q ( y ) = δ h ( Π N y ) for y ∈ H y\in H y ∈ H . Put
ψ σ ( y ) = − φ ( T y ) + σ Q ( y ) ( y ∈ H ) . \psi_{\sigma}(y)=-\varphi(Ty)+\sigma Q(y)\qquad(y\in H). ψ σ ( y ) = − φ ( T y ) + σ Q ( y ) ( y ∈ H ) .
By Step 3 (with c = Π N x ^ c=\Pi_{N}\hat{x} c = Π N x ^ , χ = φ \chi=\varphi χ = φ ) and Constants, Sums, Scalar Multiples and Differences of Differentiable Functions on an Open Subset of a Real Inner Product Space §scalar , Constants, Sums, Scalar Multiples and Differences of Differentiable Functions on an Open Subset of a Real Inner Product Space §sum , ψ σ ∈ C 2 ( H ) \psi_{\sigma}\in C^{2}(H) ψ σ ∈ C 2 ( H ) and, as T x ^ = x ^ T\hat{x}=\hat{x} T x ^ = x ^ ,
D ψ σ ( x ^ ) = − L D φ ( x ^ ) + σ T b N x ^ . ( 6 ) D\psi_{\sigma}(\hat{x})=-L\,D\varphi(\hat{x})+\sigma T_{b_{N}}\hat{x}.\qquad(6) D ψ σ ( x ^ ) = − L D φ ( x ^ ) + σ T b N x ^ . ( 6 )
For y ∈ V y\in V y ∈ V , Q ( y ) = δ h ( Π N y ) Q(y)=\delta h(\Pi_{N}y) Q ( y ) = δ h ( Π N y ) and (3) give ψ σ ( y ) − σ δ h ( y ) = − φ ( T y ) − σ δ h ( L y ) \psi_{\sigma}(y)-\sigma\delta h(y)=-\varphi(Ty)-\sigma\delta h(Ly) ψ σ ( y ) − σ δ h ( y ) = − φ ( T y ) − σ δ h ( L y ) , and by the properties of T T T above this equals ( w ( T y ) − σ δ h ( T y ) − φ ( T y ) ) − w ( x ^ ) + σ δ h ( Π N x ^ ) \bigl(w(Ty)-\sigma\delta h(Ty)-\varphi(Ty)\bigr)-w(\hat{x})+\sigma\delta h(\Pi_{N}\hat{x}) ( w ( T y ) − σ δ h ( T y ) − φ ( T y ) ) − w ( x ^ ) + σ δ h ( Π N x ^ ) . Subtracting the same identity at y = x ^ y=\hat{x} y = x ^ (where T x ^ = x ^ T\hat{x}=\hat{x} T x ^ = x ^ ),
( ψ σ ( y ) − σ δ h ( y ) ) − ( ψ σ ( x ^ ) − σ δ h ( x ^ ) ) = ( w − σ δ h − φ ) ( T y ) − ( w − σ δ h − φ ) ( x ^ ) . \bigl(\psi_{\sigma}(y)-\sigma\delta h(y)\bigr)-\bigl(\psi_{\sigma}(\hat{x})-\sigma\delta h(\hat{x})\bigr)=\bigl(w-\sigma\delta h-\varphi\bigr)(Ty)-\bigl(w-\sigma\delta h-\varphi\bigr)(\hat{x}). ( ψ σ ( y ) − σ δ h ( y ) ) − ( ψ σ ( x ^ ) − σ δ h ( x ^ ) ) = ( w − σ δ h − φ ) ( T y ) − ( w − σ δ h − φ ) ( x ^ ) .
If y ∈ V y\in V y ∈ V and ∣ y − x ^ ∣ H < r |y-\hat{x}|_{H}<r ∣ y − x ^ ∣ H < r , then T y ∈ V Ty\in V T y ∈ V and ∣ T y − x ^ ∣ H < r |Ty-\hat{x}|_{H}<r ∣ T y − x ^ ∣ H < r , so (5) at x = T y x=Ty x = T y shows that σ \sigma σ times the left side is at most 0 0 0 . Thus for σ = 1 \sigma=1 σ = 1 the function ψ 1 − δ h \psi_{1}-\delta h ψ 1 − δ h has a local maximum at x ^ \hat{x} x ^ relative to V V V , and First- and Second-Order Conditions at a Local Extremum of a C 2 C^2 C 2 Function Penalised by h h h on the Small Space §maximum (with λ = δ \lambda=\delta λ = δ ) gives x ^ ∈ D ( A ) \hat{x}\in D(A) x ^ ∈ D ( A ) and δ A x ^ = D ψ 1 ( x ^ ) \delta A\hat{x}=D\psi_{1}(\hat{x}) δ A x ^ = D ψ 1 ( x ^ ) ; for σ = − 1 \sigma=-1 σ = − 1 the function ψ − 1 + δ h \psi_{-1}+\delta h ψ − 1 + δ h has a local minimum at x ^ \hat{x} x ^ relative to V V V , and First- and Second-Order Conditions at a Local Extremum of a C 2 C^2 C 2 Function Penalised by h h h on the Small Space §minimum gives x ^ ∈ D ( A ) \hat{x}\in D(A) x ^ ∈ D ( A ) and δ A x ^ = − D ψ − 1 ( x ^ ) \delta A\hat{x}=-D\psi_{-1}(\hat{x}) δ A x ^ = − D ψ − 1 ( x ^ ) . In both cases x ^ ∈ D ( A ) \hat{x}\in D(A) x ^ ∈ D ( A ) and σ δ A x ^ = D ψ σ ( x ^ ) \sigma\delta A\hat{x}=D\psi_{\sigma}(\hat{x}) σ δ A x ^ = D ψ σ ( x ^ ) . Put
p σ = D φ ( x ^ ) + σ δ A x ^ ∈ H . p_{\sigma}=D\varphi(\hat{x})+\sigma\delta A\hat{x}\in H . p σ = D φ ( x ^ ) + σ δ A x ^ ∈ H .
Let k ∉ Γ N k\notin\Gamma_{N} k ∈ / Γ N . Then L e k = e k Le_{k}=e_{k} L e k = e k (Step 1(a)), so ⟨ L D φ ( x ^ ) , e k ⟩ H = ⟨ D φ ( x ^ ) , e k ⟩ H \langle LD\varphi(\hat{x}),e_{k}\rangle_{H}=\langle D\varphi(\hat{x}),e_{k}\rangle_{H} ⟨ L D φ ( x ^ ) , e k ⟩ H = ⟨ D φ ( x ^ ) , e k ⟩ H by self-adjointness of L L L , and ⟨ T b N x ^ , e k ⟩ H = b N ( x ^ , e k ) = 0 \langle T_{b_{N}}\hat{x},e_{k}\rangle_{H}=b_{N}(\hat{x},e_{k})=0 ⟨ T b N x ^ , e k ⟩ H = b N ( x ^ , e k ) = 0 because e k e_{k} e k vanishes at every k i k_{i} k i . Pairing σ δ A x ^ = D ψ σ ( x ^ ) \sigma\delta A\hat{x}=D\psi_{\sigma}(\hat{x}) σ δ A x ^ = D ψ σ ( x ^ ) with e k e_{k} e k and using (6) gives σ δ ⟨ A x ^ , e k ⟩ H = − ⟨ D φ ( x ^ ) , e k ⟩ H \sigma\delta\langle A\hat{x},e_{k}\rangle_{H}=-\langle D\varphi(\hat{x}),e_{k}\rangle_{H} σ δ ⟨ A x ^ , e k ⟩ H = − ⟨ D φ ( x ^ ) , e k ⟩ H , i.e. ⟨ p σ , e k ⟩ H = 0 \langle p_{\sigma},e_{k}\rangle_{H}=0 ⟨ p σ , e k ⟩ H = 0 , so by (1)
p σ ( k ) = 0 ( k ∉ Γ N ) . ( 7 ) p_{\sigma}(k)=0\qquad(k\notin\Gamma_{N}).\qquad(7) p σ ( k ) = 0 ( k ∈ / Γ N ) . ( 7 )
Moreover, by The Sobolev Triple of Order Two for the Wick-Square Problem: the Domain Is the Sobolev Space of Order -1, Unit Families Are Square-Summable, the Gradient Form and the Riccati Drift §domain and (1), ⟨ A x ^ , E i ⟩ H = μ k i x ^ ( k i ) / μ k i 3 = z ^ i / μ k i 2 \langle A\hat{x},E_{i}\rangle_{H}=\mu_{k_{i}}\hat{x}(k_{i})/\mu_{k_{i}}^{3}=\hat{z}_{i}/\mu_{k_{i}}^{2} ⟨ A x ^ , E i ⟩ H = μ k i x ^ ( k i ) / μ k i 3 = z ^ i / μ k i 2 for i ∈ [ m ] i\in[m] i ∈ [ m ] , hence
⟨ D φ ( x ^ ) , E i ⟩ H + σ δ z ^ i μ k i 2 = ⟨ p σ , E i ⟩ H = p σ ( k i ) μ k i 3 . ( 8 ) \langle D\varphi(\hat{x}),E_{i}\rangle_{H}+\sigma\delta\frac{\hat{z}_{i}}{\mu_{k_{i}}^{2}}=\langle p_{\sigma},E_{i}\rangle_{H}=\frac{p_{\sigma}(k_{i})}{\mu_{k_{i}}^{3}}.\qquad(8) ⟨ D φ ( x ^ ) , E i ⟩ H + σ δ μ k i 2 z ^ i = ⟨ p σ , E i ⟩ H = μ k i 3 p σ ( k i ) . ( 8 )
Step 7 (The Euclidean test function). With S ( z ) = ι N z + y 0 ∈ V S(z)=\iota_{N}z+y_{0}\in V S ( z ) = ι N z + y 0 ∈ V as in Step 4, define Ψ : R m → R \Psi:\mathbb{R}^{m}\to\mathbb{R} Ψ : R m → R by Ψ ( z ) = φ ( S ( z ) ) + σ δ h ( S ( z ) ) \Psi(z)=\varphi(S(z))+\sigma\delta h(S(z)) Ψ ( z ) = φ ( S ( z )) + σ δ h ( S ( z )) . Since y 0 ∈ V y_{0}\in V y 0 ∈ V vanishes at every k i k_{i} k i , (2) gives
Ψ ( z ) = φ 0 ( z ) + σ δ h ( y 0 ) + σ δ 2 ∑ i = 1 m z i 2 μ k i 2 , φ 0 = φ ∘ S . \Psi(z)=\varphi_{0}(z)+\sigma\delta h(y_{0})+\frac{\sigma\delta}{2}\sum_{i=1}^{m}\frac{z_{i}^{2}}{\mu_{k_{i}}^{2}},\qquad \varphi_{0}=\varphi\circ S . Ψ ( z ) = φ 0 ( z ) + σ δ h ( y 0 ) + 2 σ δ i = 1 ∑ m μ k i 2 z i 2 , φ 0 = φ ∘ S .
By Step 4, φ 0 \varphi_{0} φ 0 is of class C 2 C^{2} C 2 on R m \mathbb{R}^{m} R m ; the remaining terms form a function built from constants and coordinate functions by sums, scalar multiples and products, which is smooth on R m \mathbb{R}^{m} R m by claims 2 and 3 of Constants, Coordinate Functions, Sums and Products of C k C^k C k Functions on a Euclidean Open Set , and hence of class C 2 C^{2} C 2 by Smooth Map on a Euclidean Open Set ; so Ψ \Psi Ψ is of class C 2 C^{2} C 2 on R m \mathbb{R}^{m} R m by claim 3 of Constants, Coordinate Functions, Sums and Products of C k C^k C k Functions on a Euclidean Open Set . By claim 1 there (the i i i th partial of the l l l th coordinate function being 1 1 1 if l = i l=i l = i and 0 0 0 otherwise, directly from Partial Derivative on a Euclidean Open Set ) and Step 4, for z ∈ R m z\in\mathbb{R}^{m} z ∈ R m and i , l ∈ [ m ] i,l\in[m] i , l ∈ [ m ] ,
∂ i Ψ ( z ) = ⟨ D φ ( S ( z ) ) , E i ⟩ H + σ δ z i μ k i 2 , ∂ l ∂ i Ψ ( z ) = D 2 φ ( S ( z ) ) ( E l , E i ) + { σ δ μ k i − 2 ( l = i ) , 0 ( l ≠ i ) . \partial_{i}\Psi(z)=\langle D\varphi(S(z)),E_{i}\rangle_{H}+\sigma\delta\frac{z_{i}}{\mu_{k_{i}}^{2}},\qquad
\partial_{l}\partial_{i}\Psi(z)=D^{2}\varphi(S(z))(E_{l},E_{i})+\begin{cases}\sigma\delta\mu_{k_{i}}^{-2}&(l=i),\\0&(l\neq i).\end{cases} ∂ i Ψ ( z ) = ⟨ D φ ( S ( z )) , E i ⟩ H + σ δ μ k i 2 z i , ∂ l ∂ i Ψ ( z ) = D 2 φ ( S ( z )) ( E l , E i ) + { σ δ μ k i − 2 0 ( l = i ) , ( l = i ) .
For z ∈ R m z\in\mathbb{R}^{m} z ∈ R m , p r N S ( z ) = z + p r N y 0 = z \mathrm{pr}_{N}S(z)=z+\mathrm{pr}_{N}y_{0}=z pr N S ( z ) = z + pr N y 0 = z , so w ( S ( z ) ) = Z ( z ) w(S(z))=Z(z) w ( S ( z )) = Z ( z ) and
Z ( z ) − Ψ ( z ) = ( w − σ δ h − φ ) ( S ( z ) ) , S ( z ^ ) = x ^ , ∣ S ( z ) − x ^ ∣ H ≤ ∥ z − z ^ ∥ Z(z)-\Psi(z)=\bigl(w-\sigma\delta h-\varphi\bigr)(S(z)),\qquad S(\hat{z})=\hat{x},\qquad |S(z)-\hat{x}|_{H}\le\lVert z-\hat{z}\rVert Z ( z ) − Ψ ( z ) = ( w − σ δ h − φ ) ( S ( z )) , S ( z ^ ) = x ^ , ∣ S ( z ) − x ^ ∣ H ≤ ∥ z − z ^ ∥
by (4). Hence every z z z with d E ( z , z ^ ) < r d_{E}(z,\hat{z})<r d E ( z , z ^ ) < r satisfies σ [ ( Z − Ψ ) ( z ) − ( Z − Ψ ) ( z ^ ) ] ≤ 0 \sigma[(Z-\Psi)(z)-(Z-\Psi)(\hat{z})]\le0 σ [( Z − Ψ ) ( z ) − ( Z − Ψ ) ( z ^ )] ≤ 0 by (5): Z − Ψ Z-\Psi Z − Ψ has a local maximum (if σ = 1 \sigma=1 σ = 1 ), respectively minimum (if σ = − 1 \sigma=-1 σ = − 1 ), at z ^ \hat{z} z ^ relative to R m \mathbb{R}^{m} R m . The defining property in Viscosity Subsolution and Supersolution of a Second-Order Equation , applied to the test function Ψ \Psi Ψ , gives
σ F N ♯ ( z ^ , Z ( z ^ ) , D Ψ ( z ^ ) , D 2 Ψ ( z ^ ) ) ≤ 0. ( 9 ) \sigma\,\mathcal{F}^{\sharp}_{N}\bigl(\hat{z},Z(\hat{z}),D\Psi(\hat{z}),D^{2}\Psi(\hat{z})\bigr)\le0.\qquad(9) σ F N ♯ ( z ^ , Z ( z ^ ) , D Ψ ( z ^ ) , D 2 Ψ ( z ^ ) ) ≤ 0. ( 9 )
Step 8 (Translation term by term). Let s = w δ σ ( x ^ ) = w ( x ^ ) − σ δ h ( x ^ ) s=w^{\sigma}_{\delta}(\hat{x})=w(\hat{x})-\sigma\delta h(\hat{x}) s = w δ σ ( x ^ ) = w ( x ^ ) − σ δ h ( x ^ ) and X σ = D 2 φ ( x ^ ) ∣ V + σ δ I V X_{\sigma}=D^{2}\varphi(\hat{x})|_{V}+\sigma\delta I_{V} X σ = D 2 φ ( x ^ ) ∣ V + σ δ I V . Since x ^ ∈ D ( A ) \hat{x}\in D(A) x ^ ∈ D ( A ) (Step 6), Step 5 gives
F σ ( x ^ , s , D φ ( x ^ ) , D 2 φ ( x ^ ) ) = F N ♯ ( x ^ , w ( x ^ ) , p σ , X σ ) , F^{\sigma}\bigl(\hat{x},s,D\varphi(\hat{x}),D^{2}\varphi(\hat{x})\bigr)=F^{\sharp}_{N}\bigl(\hat{x},w(\hat{x}),p_{\sigma},X_{\sigma}\bigr), F σ ( x ^ , s , D φ ( x ^ ) , D 2 φ ( x ^ ) ) = F N ♯ ( x ^ , w ( x ^ ) , p σ , X σ ) ,
and w ( x ^ ) = Z ( z ^ ) w(\hat{x})=Z(\hat{z}) w ( x ^ ) = Z ( z ^ ) . We compare the terms of The Shifted Hamilton-Jacobi-Bellman Equation of the Wick-Square Problem with Noise, Control and Cost Cut Off to a Cube of Modes §operator with those of The Galerkin Hamilton-Jacobi-Bellman Equations of the Wick-Square Problem in the Coordinates of a Cube of Modes §shifted , using S ( z ^ ) = x ^ S(\hat{z})=\hat{x} S ( z ^ ) = x ^ .
Gradient. By Step 7 and (8), ∂ i Ψ ( z ^ ) = ⟨ p σ , E i ⟩ H \partial_{i}\Psi(\hat{z})=\langle p_{\sigma},E_{i}\rangle_{H} ∂ i Ψ ( z ^ ) = ⟨ p σ , E i ⟩ H for every i i i . The gradient D Ψ ( z ^ ) D\Psi(\hat{z}) D Ψ ( z ^ ) is the vector of the partials ∂ i Ψ ( z ^ ) \partial_{i}\Psi(\hat{z}) ∂ i Ψ ( z ^ ) (Differential Calculus and Convexity on Euclidean Open Sets: Standing Notation §derivatives , Gradient of a Real-Valued Function on a Euclidean Open Set ), and the squared Euclidean norm of a vector is the sum of the squares of its coordinates (Euclidean Norm on R n \mathbb{R}^n R n ). Hence 1 2 ∥ D Ψ ( z ^ ) ∥ 2 = 1 2 ∑ i = 1 m ⟨ p σ , E i ⟩ H 2 = 1 2 ∑ k ∈ Γ N ⟨ p σ , e k ⟩ H 2 = 1 2 G N ( p σ , p σ ) \tfrac12\lVert D\Psi(\hat{z})\rVert^{2}=\tfrac12\sum_{i=1}^{m}\langle p_{\sigma},E_{i}\rangle_{H}^{2}=\tfrac12\sum_{k\in\Gamma_{N}}\langle p_{\sigma},e_{k}\rangle_{H}^{2}=\tfrac12 G_{N}(p_{\sigma},p_{\sigma}) 2 1 ∥ D Ψ ( z ^ ) ∥ 2 = 2 1 ∑ i = 1 m ⟨ p σ , E i ⟩ H 2 = 2 1 ∑ k ∈ Γ N ⟨ p σ , e k ⟩ H 2 = 2 1 G N ( p σ , p σ ) , with G N G_{N} G N as in The Shifted Hamilton-Jacobi-Bellman Equation of the Wick-Square Problem with Noise, Control and Cost Cut Off to a Cube of Modes §data .
Drift. By The Galerkin Problems of the Wick-Square Problem: Standing Notation for Cubes of Modes, Coordinates, the Free Galerkin Quadratic and the Riccati Potential §functions and claim 1 of Constants, Coordinate Functions, Sums and Products of C k C^k C k Functions on a Euclidean Open Set , ∂ i Φ N ( z ^ ) = ( μ k i + 2 q k i ) z ^ i \partial_{i}\Phi_{N}(\hat{z})=(\mu_{k_{i}}+2q_{k_{i}})\hat{z}_{i} ∂ i Φ N ( z ^ ) = ( μ k i + 2 q k i ) z ^ i ; since D Φ N ( z ^ ) D\Phi_{N}(\hat{z}) D Φ N ( z ^ ) and D Ψ ( z ^ ) D\Psi(\hat{z}) D Ψ ( z ^ ) are the vectors of partials (Gradient of a Real-Valued Function on a Euclidean Open Set ) and the dot product is the sum of the coordinatewise products (Difference, Dot Product, and Orthogonality in R n \mathbb{R}^n R n ), and z ^ i = x ^ ( k i ) \hat{z}_{i}=\hat{x}(k_{i}) z ^ i = x ^ ( k i ) ,
D Φ N ( z ^ ) ⋅ D Ψ ( z ^ ) = ∑ i = 1 m ( μ k i + 2 q k i ) x ^ ( k i ) ⟨ p σ , E i ⟩ H . D\Phi_{N}(\hat{z})\cdot D\Psi(\hat{z})=\sum_{i=1}^{m}\bigl(\mu_{k_{i}}+2q_{k_{i}}\bigr)\hat{x}(k_{i})\,\langle p_{\sigma},E_{i}\rangle_{H}. D Φ N ( z ^ ) ⋅ D Ψ ( z ^ ) = i = 1 ∑ m ( μ k i + 2 q k i ) x ^ ( k i ) ⟨ p σ , E i ⟩ H .
On the other side, x ^ ∈ H − 1 \hat{x}\in H^{-1} x ^ ∈ H − 1 and p σ ∈ H p_{\sigma}\in H p σ ∈ H , so by The Sobolev Triple of Order Two for the Wick-Square Problem: the Domain Is the Sobolev Space of Order -1, Unit Families Are Square-Summable, the Gradient Form and the Riccati Drift §pairings ⟨ A x ^ , p σ ⟩ H \langle A\hat{x},p_{\sigma}\rangle_{H} ⟨ A x ^ , p σ ⟩ H is the lattice sum of k ↦ x ^ ( k ) p σ ( k ) / μ k 2 k\mapsto\hat{x}(k)p_{\sigma}(k)/\mu_{k}^{2} k ↦ x ^ ( k ) p σ ( k ) / μ k 2 . By (7) this family vanishes off Γ N \Gamma_{N} Γ N ; for M ≥ N M\ge N M ≥ N one has Γ N ⊆ Γ M \Gamma_{N}\subseteq\Gamma_{M} Γ N ⊆ Γ M (The Wick-Square Problem on the Torus: Standing Notation §cubes ), so the M M M th cube sum equals ∑ k ∈ Γ N x ^ ( k ) p σ ( k ) / μ k 2 \sum_{k\in\Gamma_{N}}\hat{x}(k)p_{\sigma}(k)/\mu_{k}^{2} ∑ k ∈ Γ N x ^ ( k ) p σ ( k ) / μ k 2 , and by Cube Sums of Families on the Integer Lattice §lattice-sum so does the lattice sum. With (8),
⟨ A x ^ , p σ ⟩ H = ∑ i = 1 m x ^ ( k i ) p σ ( k i ) μ k i 2 = ∑ i = 1 m μ k i x ^ ( k i ) ⟨ p σ , E i ⟩ H . \langle A\hat{x},p_{\sigma}\rangle_{H}=\sum_{i=1}^{m}\frac{\hat{x}(k_{i})p_{\sigma}(k_{i})}{\mu_{k_{i}}^{2}}=\sum_{i=1}^{m}\mu_{k_{i}}\hat{x}(k_{i})\langle p_{\sigma},E_{i}\rangle_{H}. ⟨ A x ^ , p σ ⟩ H = i = 1 ∑ m μ k i 2 x ^ ( k i ) p σ ( k i ) = i = 1 ∑ m μ k i x ^ ( k i ) ⟨ p σ , E i ⟩ H .
By The Shifted Hamilton-Jacobi-Bellman Equation of the Wick-Square Problem with Noise, Control and Cost Cut Off to a Cube of Modes §data and Step 1(a), B q N x ^ = Π N ( B q x ^ ) = ∑ i = 1 m 2 q k i x ^ ( k i ) E i B^{N}_{q}\hat{x}=\Pi_{N}(B_{q}\hat{x})=\sum_{i=1}^{m}2q_{k_{i}}\hat{x}(k_{i})E_{i} B q N x ^ = Π N ( B q x ^ ) = ∑ i = 1 m 2 q k i x ^ ( k i ) E i , so by Inner Products Against Finite Sums, and Orthonormal Families, in a Real Inner Product Space §combinations ⟨ B q N x ^ , p σ ⟩ H = ∑ i = 1 m 2 q k i x ^ ( k i ) ⟨ E i , p σ ⟩ H \langle B^{N}_{q}\hat{x},p_{\sigma}\rangle_{H}=\sum_{i=1}^{m}2q_{k_{i}}\hat{x}(k_{i})\langle E_{i},p_{\sigma}\rangle_{H} ⟨ B q N x ^ , p σ ⟩ H = ∑ i = 1 m 2 q k i x ^ ( k i ) ⟨ E i , p σ ⟩ H . Adding, by Elementary Identities in a Real Inner Product Space §bilinear ,
⟨ A x ^ + B q N x ^ , p σ ⟩ H = D Φ N ( z ^ ) ⋅ D Ψ ( z ^ ) . \langle A\hat{x}+B^{N}_{q}\hat{x},p_{\sigma}\rangle_{H}=D\Phi_{N}(\hat{z})\cdot D\Psi(\hat{z}). ⟨ A x ^ + B q N x ^ , p σ ⟩ H = D Φ N ( z ^ ) ⋅ D Ψ ( z ^ ) .
Trace. The Hessian D 2 Ψ ( z ^ ) D^{2}\Psi(\hat{z}) D 2 Ψ ( z ^ ) is the matrix with entries ∂ i ∂ l Ψ ( z ^ ) \partial_{i}\partial_{l}\Psi(\hat{z}) ∂ i ∂ l Ψ ( z ^ ) (Differential Calculus and Convexity on Euclidean Open Sets: Standing Notation §derivatives , Hessian Matrix of a C^2 Function ), so by Trace of a Real Square Matrix its trace is ∑ i = 1 m ∂ i ∂ i Ψ ( z ^ ) \sum_{i=1}^{m}\partial_{i}\partial_{i}\Psi(\hat{z}) ∑ i = 1 m ∂ i ∂ i Ψ ( z ^ ) , and by Step 7, tr D 2 Ψ ( z ^ ) = ∑ i = 1 m ( D 2 φ ( x ^ ) ( E i , E i ) + σ δ μ k i − 2 ) = ∑ i = 1 m X σ ( E i , E i ) \operatorname{tr}D^{2}\Psi(\hat{z})=\sum_{i=1}^{m}\bigl(D^{2}\varphi(\hat{x})(E_{i},E_{i})+\sigma\delta\mu_{k_{i}}^{-2}\bigr)=\sum_{i=1}^{m}X_{\sigma}(E_{i},E_{i}) tr D 2 Ψ ( z ^ ) = ∑ i = 1 m ( D 2 φ ( x ^ ) ( E i , E i ) + σ δ μ k i − 2 ) = ∑ i = 1 m X σ ( E i , E i ) , since X σ ( E i , E i ) = D 2 φ ( x ^ ) ( E i , E i ) + σ δ ⟨ E i , E i ⟩ V X_{\sigma}(E_{i},E_{i})=D^{2}\varphi(\hat{x})(E_{i},E_{i})+\sigma\delta\langle E_{i},E_{i}\rangle_{V} X σ ( E i , E i ) = D 2 φ ( x ^ ) ( E i , E i ) + σ δ ⟨ E i , E i ⟩ V by Bounded Symmetric Bilinear Forms on a Real Inner Product Space: Norm, Order, Identity Form and Restriction §identity and Bounded Symmetric Bilinear Forms on a Real Inner Product Space: Norm, Order, Identity Form and Restriction §restriction , and ⟨ E i , E i ⟩ V = μ k i − 2 \langle E_{i},E_{i}\rangle_{V}=\mu_{k_{i}}^{-2} ⟨ E i , E i ⟩ V = μ k i − 2 by Step 1. On the other side, with κ \kappa κ the enumeration of The Wick-Square Problem on the Sobolev Triple of Order Two: Standing Notation §noise , a bijection of N \mathbb{N} N onto Z n \mathbb{Z}^{n} Z n (an enumeration being by definition such a bijection, Properties of the Fourier Coefficients on the Torus, and the Realisation of Weighted Coefficient Families ), the set J = { j ∈ N : κ ( j ) ∈ Γ N } J=\{j\in\mathbb{N}:\kappa(j)\in\Gamma_{N}\} J = { j ∈ N : κ ( j ) ∈ Γ N } is finite and nonempty and κ \kappa κ maps it bijectively onto Γ N \Gamma_{N} Γ N ; the terms X σ ( f j N , f j N ) X_{\sigma}(f^{N}_{j},f^{N}_{j}) X σ ( f j N , f j N ) vanish for j ∉ J j\notin J j ∈ / J , since then f j N = 0 H f^{N}_{j}=0_{H} f j N = 0 H and X σ X_{\sigma} X σ is bilinear. So every partial sum of index at least max J \max J max J of the series defining T r f N X σ \mathrm{Tr}_{f^{N}}X_{\sigma} Tr f N X σ in Square-Summable Sequences in the Form Space of a Hilbert Triple and the Trace of a Form along Them §trace equals ∑ j ∈ J X σ ( e κ ( j ) , e κ ( j ) ) = ∑ k ∈ Γ N X σ ( e k , e k ) = ∑ i = 1 m X σ ( E i , E i ) \sum_{j\in J}X_{\sigma}(e_{\kappa(j)},e_{\kappa(j)})=\sum_{k\in\Gamma_{N}}X_{\sigma}(e_{k},e_{k})=\sum_{i=1}^{m}X_{\sigma}(E_{i},E_{i}) ∑ j ∈ J X σ ( e κ ( j ) , e κ ( j ) ) = ∑ k ∈ Γ N X σ ( e k , e k ) = ∑ i = 1 m X σ ( E i , E i ) , and the series converges to this number. Hence tr D 2 Ψ ( z ^ ) = T r f N X σ \operatorname{tr}D^{2}\Psi(\hat{z})=\mathrm{Tr}_{f^{N}}X_{\sigma} tr D 2 Ψ ( z ^ ) = Tr f N X σ .
Running cost. ι N z ^ = Π N x ^ ∈ H − 1 \iota_{N}\hat{z}=\Pi_{N}\hat{x}\in H^{-1} ι N z ^ = Π N x ^ ∈ H − 1 , so the running cost g ∣ H − 1 g|_{H^{-1}} g ∣ H − 1 of F N ♯ \mathcal{F}^{\sharp}_{N} F N ♯ takes the value g ( Π N x ^ ) g(\Pi_{N}\hat{x}) g ( Π N x ^ ) at ι N z ^ \iota_{N}\hat{z} ι N z ^ .
Combining these identities,
F σ ( x ^ , s , D φ ( x ^ ) , D 2 φ ( x ^ ) ) = γ Z ( z ^ ) − ν 2 T r f N X σ + 1 2 G N ( p σ , p σ ) + ⟨ A x ^ + B q N x ^ , p σ ⟩ H − g ( Π N x ^ ) = F N ♯ ( z ^ , Z ( z ^ ) , D Ψ ( z ^ ) , D 2 Ψ ( z ^ ) ) , F^{\sigma}\bigl(\hat{x},s,D\varphi(\hat{x}),D^{2}\varphi(\hat{x})\bigr)=\gamma Z(\hat{z})-\tfrac{\nu}{2}\mathrm{Tr}_{f^{N}}X_{\sigma}+\tfrac12G_{N}(p_{\sigma},p_{\sigma})+\langle A\hat{x}+B^{N}_{q}\hat{x},p_{\sigma}\rangle_{H}-g(\Pi_{N}\hat{x})=\mathcal{F}^{\sharp}_{N}\bigl(\hat{z},Z(\hat{z}),D\Psi(\hat{z}),D^{2}\Psi(\hat{z})\bigr), F σ ( x ^ , s , D φ ( x ^ ) , D 2 φ ( x ^ ) ) = γ Z ( z ^ ) − 2 ν Tr f N X σ + 2 1 G N ( p σ , p σ ) + ⟨ A x ^ + B q N x ^ , p σ ⟩ H − g ( Π N x ^ ) = F N ♯ ( z ^ , Z ( z ^ ) , D Ψ ( z ^ ) , D 2 Ψ ( z ^ ) ) ,
and therefore, by (9), σ F σ ( x ^ , s , D φ ( x ^ ) , D 2 φ ( x ^ ) ) ≤ 0 \sigma F^{\sigma}(\hat{x},s,D\varphi(\hat{x}),D^{2}\varphi(\hat{x}))\le0 σ F σ ( x ^ , s , D φ ( x ^ ) , D 2 φ ( x ^ )) ≤ 0 .
Step 9 (Conclusion). Take y = x ^ y=\hat{x} y = x ^ , which lies in D ( A ) = W D(A)=W D ( A ) = W by Step 6, the number s = w δ σ ( x ^ ) s=w^{\sigma}_{\delta}(\hat{x}) s = w δ σ ( x ^ ) of Step 8, the vector D φ ( x ^ ) ∈ H D\varphi(\hat{x})\in H D φ ( x ^ ) ∈ H in the gradient slot of Viscosity Subsolution, Supersolution and Solution of a Second-Order Equation on a Hilbert Triple and the form Y = D 2 φ ( x ^ ) ∈ S y m ( H ) Y=D^{2}\varphi(\hat{x})\in\mathrm{Sym}(H) Y = D 2 φ ( x ^ ) ∈ Sym ( H ) . Then ∣ y − x ^ ∣ H = 0 |y-\hat{x}|_{H}=0 ∣ y − x ^ ∣ H = 0 , ∣ w δ σ ( y ) − w δ σ ( x ^ ) ∣ = 0 |w^{\sigma}_{\delta}(y)-w^{\sigma}_{\delta}(\hat{x})|=0 ∣ w δ σ ( y ) − w δ σ ( x ^ ) ∣ = 0 , ∣ s − w δ σ ( x ^ ) ∣ = 0 |s-w^{\sigma}_{\delta}(\hat{x})|=0 ∣ s − w δ σ ( x ^ ) ∣ = 0 , ∣ D φ ( x ^ ) − D φ ( x ^ ) ∣ H = 0 |D\varphi(\hat{x})-D\varphi(\hat{x})|_{H}=0 ∣ D φ ( x ^ ) − D φ ( x ^ ) ∣ H = 0 and ∥ Y − D 2 φ ( x ^ ) ∥ = ∥ 0 S y m ∥ = 0 \lVert Y-D^{2}\varphi(\hat{x})\rVert=\lVert0_{\mathrm{Sym}}\rVert=0 ∥ Y − D 2 φ ( x ^ )∥ = ∥ 0 Sym ∥ = 0 (Elementary Properties of Bounded Symmetric Bilinear Forms: Norm, Quadratic Form, Order and Continuity §norm-axioms ), all less than ε \varepsilon ε . If σ = 1 \sigma=1 σ = 1 , Step 8 gives ( F N ♯ ) δ − ( y , s , D φ ( x ^ ) , Y ) ≤ 0 < ε (F^{\sharp}_{N})^{-}_{\delta}(y,s,D\varphi(\hat{x}),Y)\le0<\varepsilon ( F N ♯ ) δ − ( y , s , D φ ( x ^ ) , Y ) ≤ 0 < ε ; if σ = − 1 \sigma=-1 σ = − 1 , it gives − ε < 0 ≤ ( F N ♯ ) δ + ( y , s , D φ ( x ^ ) , Y ) -\varepsilon<0\le(F^{\sharp}_{N})^{+}_{\delta}(y,s,D\varphi(\hat{x}),Y) − ε < 0 ≤ ( F N ♯ ) δ + ( y , s , D φ ( x ^ ) , Y ) . Since δ \delta δ , φ \varphi φ , x ^ \hat{x} x ^ and ε \varepsilon ε were arbitrary, w w w is a viscosity subsolution of F N ♯ F^{\sharp}_{N} F N ♯ on H H H when Z Z Z is a viscosity subsolution of F N ♯ \mathcal{F}^{\sharp}_{N} F N ♯ on R m \mathbb{R}^{m} R m , by Viscosity Subsolution, Supersolution and Solution of a Second-Order Equation on a Hilbert Triple §subsolution , which is clause 2 ; and w w w is a viscosity supersolution of F N ♯ F^{\sharp}_{N} F N ♯ on H H H when Z Z Z is a viscosity supersolution of F N ♯ \mathcal{F}^{\sharp}_{N} F N ♯ on R m \mathbb{R}^{m} R m , by Viscosity Subsolution, Supersolution and Solution of a Second-Order Equation on a Hilbert Triple §supersolution , which is clause 3 .