Put r = m − 1 r=m-1 r = m − 1 , T t = κ m ( t ) ( u ) T_t=\kappa_m(t)(u) T t = κ m ( t ) ( u ) and
W i , t = κ m + 1 ( t ) ( u , ⋅ , … , ⋅ , A t − 1 / 2 e i ) W_{i,t}=\kappa_{m+1}(t)(u,\cdot,\ldots,\cdot,A_t^{-1/2}e_i) W i , t = κ m + 1 ( t ) ( u , ⋅ , … , ⋅ , A t − 1/2 e i ) .
The dynamics give
d A = ∑ i H i d B i − A d t , d T = ∑ i W i d B i − ( m T + L m ( u ) ) d t . dA=\sum_iH_i\,dB_i-A\,dt,\qquad
dT=\sum_iW_i\,dB_i-(mT+L_m(u))\,dt. d A = i ∑ H i d B i − A d t , d T = i ∑ W i d B i − ( m T + L m ( u )) d t . For positive definite A A A define
F ( A , T ) = ⟨ T , ( A − 1 ) ⊗ r T ⟩ F(A,T)=\langle T,(A^{-1})^{\otimes r}T\rangle F ( A , T ) = ⟨ T , ( A − 1 ) ⊗ r T ⟩ . For any fixed symmetric
positive definite P P P ,
F ( A , T ) = F ( P A P , P ⊗ r T ) . F(A,T)=F(PAP,P^{\otimes r}T). F ( A , T ) = F ( P A P , P ⊗ r T ) . This identity follows from P ( P A P ) − 1 P = A − 1 P(PAP)^{-1}P=A^{-1} P ( P A P ) − 1 P = A − 1 and also holds for the
directional derivatives of F F F , with the corresponding congruence applied to
directions. At the particular time where the generator is evaluated choose
P = A − 1 / 2 P=A^{-1/2} P = A − 1/2 . This is a change of variables in a derivative at a fixed point,
not a stochastic change of coordinates; it adds no drift.
Use tildes for the transformed tensors, and let G i = P H i P G_i=PH_iP G i = P H i P .
The first derivative is evaluated in the direction
( − I , − m T ~ − L ~ ) (-I,-m\widetilde T-\widetilde L) ( − I , − m T − L ) and the quadratic directions are
( G i , W ~ i ) (G_i,\widetilde W_i) ( G i , W i ) .
For symmetric H H H , write H ( s ) H^{(s)} H ( s ) for its action on slot s s s of an order-r r r
tensor, and put S H = ∑ s = 1 r H ( s ) S_H=\sum_{s=1}^r H^{(s)} S H = ∑ s = 1 r H ( s ) . Different slots commute. Multiplying
the expansions of ( I + ε H ) − 1 (I+\varepsilon H)^{-1} ( I + ε H ) − 1 in the r r r slots gives
( ( I + ε H ) − 1 ) ⊗ r = I − ε S H + ε 2 2 ( S H 2 + ∑ s = 1 r ( H ( s ) ) 2 ) + O ( ε 3 ) . ((I+\varepsilon H)^{-1})^{\otimes r}
=I-\varepsilon S_H+\frac{\varepsilon^2}{2}
\left(S_H^2+\sum_{s=1}^r(H^{(s)})^2\right)+O(\varepsilon^3). (( I + ε H ) − 1 ) ⊗ r = I − ε S H + 2 ε 2 ( S H 2 + s = 1 ∑ r ( H ( s ) ) 2 ) + O ( ε 3 ) . It follows by multiplying on each side by T ~ + ε W \widetilde T+\varepsilon W T + ε W that
the coefficient of ε \varepsilon ε in F ( I + ε H , T ~ + ε W ) F(I+\varepsilon H,\widetilde T+\varepsilon W) F ( I + ε H , T + ε W )
is
2 ⟨ T ~ , W ⟩ − ⟨ T ~ , S H T ~ ⟩ , 2\langle\widetilde T,W\rangle-\langle\widetilde T,S_H\widetilde T\rangle, 2 ⟨ T , W ⟩ − ⟨ T , S H T ⟩ , and its coefficient of ε 2 \varepsilon^2 ε 2 is
∣ W ∣ 2 − 2 ⟨ W , S H T ~ ⟩ + 1 2 ∣ S H T ~ ∣ 2 + 1 2 ∑ s = 1 r ∣ H ( s ) T ~ ∣ 2 . |W|^2-2\langle W,S_H\widetilde T\rangle
+\tfrac12|S_H\widetilde T|^2
+\tfrac12\sum_{s=1}^r|H^{(s)}\widetilde T|^2. ∣ W ∣ 2 − 2 ⟨ W , S H T ⟩ + 2 1 ∣ S H T ∣ 2 + 2 1 s = 1 ∑ r ∣ H ( s ) T ∣ 2 . These are the first derivative and half the second derivative respectively.
The deterministic direction has S − I = − r I S_{-I}=-rI S − I = − r I and therefore contributes
− ( m + 1 ) ∣ T ~ ∣ 2 − 2 ⟨ T ~ , L ~ ⟩ -(m+1)|\widetilde T|^2-2\langle\widetilde T,\widetilde L\rangle − ( m + 1 ) ∣ T ∣ 2 − 2 ⟨ T , L ⟩ .
Consequently the full drift equals
− ( m + 1 ) ∣ T ~ ∣ 2 − 2 ⟨ T ~ , L ~ ⟩ + ∑ i ∣ W ~ i ∣ 2 − 2 ∑ i ⟨ W ~ i , S G i T ~ ⟩ + 1 2 ∑ i ∣ S G i T ~ ∣ 2 + 1 2 ∑ i , s ∣ G i ( s ) T ~ ∣ 2 . \begin{split}
&-(m+1)|\widetilde T|^2-2\langle\widetilde T,\widetilde L\rangle
+\sum_i|\widetilde W_i|^2
-2\sum_i\langle\widetilde W_i,S_{G_i}\widetilde T\rangle\\
&\hspace{12mm}+\tfrac12\sum_i|S_{G_i}\widetilde T|^2
+\tfrac12\sum_{i,s}|G_i^{(s)}\widetilde T|^2.
\end{split} − ( m + 1 ) ∣ T ∣ 2 − 2 ⟨ T , L ⟩ + i ∑ ∣ W i ∣ 2 − 2 i ∑ ⟨ W i , S G i T ⟩ + 2 1 i ∑ ∣ S G i T ∣ 2 + 2 1 i , s ∑ ∣ G i ( s ) T ∣ 2 . Discard the last two nonnegative sums. Young’s inequality bounds the absolute
value of the cross term by
1 2 ∑ i ∣ W ~ i ∣ 2 + 2 ∑ i ∣ S G i T ~ ∣ 2 \frac12\sum_i|\widetilde W_i|^2+2\sum_i|S_{G_i}\widetilde T|^2 2 1 ∑ i ∣ W i ∣ 2 + 2 ∑ i ∣ S G i T ∣ 2 .
The third-order matrix estimate is ∑ i G i 2 ≤ 8 I \sum_iG_i^2\le8I ∑ i G i 2 ≤ 8 I . Thus
∑ i ∣ S G i T ~ ∣ 2 ≤ r ∑ s = 1 r ∑ i ∣ G i ( s ) T ~ ∣ 2 ≤ 8 r 2 ∣ T ~ ∣ 2 . \sum_i|S_{G_i}\widetilde T|^2
\le r\sum_{s=1}^r\sum_i|G_i^{(s)}\widetilde T|^2
\le8r^2|\widetilde T|^2. i ∑ ∣ S G i T ∣ 2 ≤ r s = 1 ∑ r i ∑ ∣ G i ( s ) T ∣ 2 ≤ 8 r 2 ∣ T ∣ 2 . Also ∣ T ~ ∣ 2 = E m |\widetilde T|^2=\mathcal E_m ∣ T ∣ 2 = E m ,
∑ i ∣ W ~ i ∣ 2 = E m + 1 \sum_i|\widetilde W_i|^2=\mathcal E_{m+1} ∑ i ∣ W i ∣ 2 = E m + 1 and
⟨ T ~ , L ~ ⟩ = ⟨ κ m ( u ) , L m ( u ) ⟩ t \langle\widetilde T,\widetilde L\rangle=
\langle\kappa_m(u),L_m(u)\rangle_t ⟨ T , L ⟩ = ⟨ κ m ( u ) , L m ( u ) ⟩ t . The drift lower bound follows because
m + 1 + 16 ( m − 1 ) 2 ≤ 17 m 2 m+1+16(m-1)^2\le17m^2 m + 1 + 16 ( m − 1 ) 2 ≤ 17 m 2 for m ≥ 3 m\ge3 m ≥ 3 .
By the finite-time integrability established in Lemma 8.3 ,
the expectation of the Itô stochastic integral on [ 0 , T ] [0,T] [ 0 , T ] is zero. Equivalently,
one may first stop in the parameter space and then remove the stop using its
deterministic finite-order posterior energy bounds. Integrating the drift gives
E m ( 0 ) + 1 2 ∫ 0 T E E m + 1 ( t ) d t ≤ E E m ( T ) + C m 2 ∫ 0 T E E m ( t ) d t + 2 ∫ 0 T E ⟨ κ m ( t ) ( u ) , L m ( t ) ( u ) ⟩ t d t . \begin{split}
\mathcal E_m(0)+\tfrac12\int_0^T\mathbb E\mathcal E_{m+1}(t)\,dt
\le{}&\mathbb E\mathcal E_m(T)+Cm^2\int_0^T\mathbb E\mathcal E_m(t)\,dt\\
&+2\int_0^T\mathbb E\langle\kappa_m(t)(u),L_m(t)(u)\rangle_t\,dt.
\end{split} E m ( 0 ) + 2 1 ∫ 0 T E E m + 1 ( t ) d t ≤ E E m ( T ) + C m 2 ∫ 0 T E E m ( t ) d t + 2 ∫ 0 T E ⟨ κ m ( t ) ( u ) , L m ( t ) ( u ) ⟩ t d t . The absolute value of the last integral is at most I m J m \sqrt{I_mJ_m} I m J m by
Cauchy–Schwarz on the product of time and probability spaces. Since the
nonnegative function t ↦ E E m ( t ) t\mapsto\mathbb E\mathcal E_m(t) t ↦ E E m ( t ) is integrable, there
exists an increasing deterministic sequence T j → ∞ T_j\to\infty T j → ∞ for which
E E m ( T j ) → 0 \mathbb E\mathcal E_m(T_j)\to0 E E m ( T j ) → 0 . In fact one can select T j T_j T j in successive
unit intervals whose energy integral tends to zero. Monotone convergence on
the two nonnegative energy integrals along this sequence proves the claimed
inequality, and proves finiteness of the next-order integral at the same time.
Integrability alone is not being used to assert terminal decay at every time.
For arbitrary fixed u u u , all terms are homogeneous of degree two in u u u , so the
same conclusions follow by scaling; the case u = 0 u=0 u = 0 has all terms zero.