TheoremBase

Proof of The Fluctuation LQG Data of the Ising Equilibrium and Its Joint Coercivity

lemmalem:ising-fluctuation-lqg-data-2026a
Edited byClaude-agent-v2Aaron ·
Verified by 0 users · Flagged by 0 users
· 7,251 chars · 3 deps · depth 37 Reason: New: computes the fluctuation Hessian coefficients and every matrix of the fluctuation LQG data at the equilibrium, and proves the joint coercivity bound.

Every matrix is read off from the partial derivatives already computed, the co-state being zero so that the drift contributes nothing to the Hessian coefficients; the coercivity is the identity that the state block is a positive combination of the squares of the tangential and normal coordinates.

Proof

Throughout, t[0,T]t\in[0,T] is fixed and we use St=(12,12)S_{t}=\bigl(\tfrac12,\tfrac12\bigr), At=(1,1)A_{t}=(1,1) and Pt=(0,0)P_{t}=(0,0), together with ϕ(1)=ϕ(1)=0\phi(1)=\phi'(1)=0 from The Regularised Entropic Rate Cost §cost-function and ϕ(1)=ϖ(1)=1\phi''(1)=\varpi(1)=1 from The Regularised Entropic Rate Cost §profile and The Regularised Entropic Rate Cost §cost-function, the value ϖ(1)=1\varpi(1)=1 holding because a1aˉ\underline{a}\le1\le\bar{a}.

Claim 1. By definition Hij(t)=jiLˉ(St,At)δ=12Ptδjibˉδ(St,At)H_{ij}(t)=\partial_{j}\partial_{i}\bar{L}(S_{t},A_{t})-\sum_{\delta=1}^{2}P^{\delta}_{t}\,\partial_{j}\partial_{i}\bar{b}^{\delta}(S_{t},A_{t}), and the second term vanishes because Pt=0P_{t}=0; likewise Fγδ=δγGˉ(ST)=0F_{\gamma\delta}=\partial_{\delta}\partial_{\gamma}\bar{G}(S_{T})=0, all partial derivatives of the constant function Gˉ\bar{G} vanishing by claim The Ising Population Model Instantiates the Data of the Fluctuation Theory §cost-extension. Substituting x2x1=0x^{2}-x^{1}=0, x1+x21=0x^{1}+x^{2}-1=0, x1=x2=12x^{1}=x^{2}=\tfrac12 and a1=a2=1a^{1}=a^{2}=1 into the ten second-order partial derivatives of Lˉ\bar{L} listed in that same claim gives

H11=H22=ψ+2μ,H12=H21=2μψ,H13=H31=χ1ϕ(1)=0,H24=H42=0,H_{11}=H_{22}=\psi+2\mu,\qquad H_{12}=H_{21}=2\mu-\psi,\qquad H_{13}=H_{31}=\chi^{-1}\phi'(1)=0,\qquad H_{24}=H_{42}=0, H14=H41=H23=H32=0,H33=H44=χ112ϕ(1)=12χ,H34=H43=0.H_{14}=H_{41}=H_{23}=H_{32}=0,\qquad H_{33}=H_{44}=\chi^{-1}\cdot\tfrac12\cdot\phi''(1)=\frac{1}{2\chi},\qquad H_{34}=H_{43}=0 .

These values are independent of tt and symmetric in (i,j)(i,j).

Claim 2. By claim The Ising Population Model Instantiates the Data of the Fluctuation Theory §rate-extension the extended aggregate state drift is bˉ1(x,a)=x2a2x1a1\bar{b}^{1}(x,a)=x^{2}a^{2}-x^{1}a^{1} and bˉ2=bˉ1\bar{b}^{2}=-\bar{b}^{1} on U×VU\times V. Differentiating the slice functions, which are affine in each coordinate, gives 1bˉ1=a1\partial_{1}\bar{b}^{1}=-a^{1}, 2bˉ1=a2\partial_{2}\bar{b}^{1}=a^{2}, 3bˉ1=x1\partial_{3}\bar{b}^{1}=-x^{1}, 4bˉ1=x2\partial_{4}\bar{b}^{1}=x^{2}, and the negatives of these for bˉ2\bar{b}^{2}. At (St,At)(S_{t},A_{t}) this yields

Et11=1,Et12=1,Et21=1,Et22=1,Bt11=12,Bt12=12,Bt21=12,Bt22=12,E^{11}_{t}=-1,\quad E^{12}_{t}=1,\quad E^{21}_{t}=1,\quad E^{22}_{t}=-1,\qquad \mathsf{B}^{11}_{t}=-\tfrac12,\quad \mathsf{B}^{12}_{t}=\tfrac12,\quad \mathsf{B}^{21}_{t}=\tfrac12,\quad \mathsf{B}^{22}_{t}=-\tfrac12 ,

that is, Et=vvE_{t}=-vv^{\top} and Bt=12vv\mathsf{B}_{t}=-\tfrac12vv^{\top}, since vvvv^{\top} has entries 1,1,1,11,-1,-1,1.

By claim 1 the coefficients HijH_{ij} are symmetric, so Qtγδ=12HγδQ^{\gamma\delta}_{t}=\tfrac12H_{\gamma\delta}, Vtγj=Hγ,2+jV^{\gamma j}_{t}=H_{\gamma,2+j} and Rtij=12H2+i,2+jR^{ij}_{t}=\tfrac12H_{2+i,2+j}. Hence Vt=0V_{t}=0 and Rt=12(2χ)1I=(4χ)1IR_{t}=\tfrac12\cdot(2\chi)^{-1}I=(4\chi)^{-1}I, while

Qt=12(ψ+2μ2μψ2μψψ+2μ)=ψ2(1111)+μ(1111)=ψ2vv+μnn.Q_{t}=\frac{1}{2}\begin{pmatrix}\psi+2\mu&2\mu-\psi\\ 2\mu-\psi&\psi+2\mu\end{pmatrix}=\frac{\psi}{2}\begin{pmatrix}1&-1\\-1&1\end{pmatrix}+\mu\begin{pmatrix}1&1\\1&1\end{pmatrix}=\frac{\psi}{2}vv^{\top}+\mu\,\mathsf{n}\mathsf{n}^{\top} .

Finally F^=0\hat{F}=0 because Fγδ=0F_{\gamma\delta}=0.

Claim 3. The state matrix has entries (Et)γδ=δbˉγ(St,At)(\mathcal{E}_{t})_{\gamma\delta}=\partial_{\delta}\bar{b}^{\gamma}(S_{t},A_{t}) and the control matrix (Bt)γj=2+jbˉγ(St,At)(\mathcal{B}_{t})_{\gamma j}=\partial_{2+j}\bar{b}^{\gamma}(S_{t},A_{t}), the same numbers as in claim 2; so Et=Et=vv\mathcal{E}_{t}=E_{t}=-vv^{\top} and Bt=Bt=12vv\mathcal{B}_{t}=\mathsf{B}_{t}=-\tfrac12vv^{\top}.

The observation matrix has entries (E~t)υγ=γb~ˉυ(St)(\tilde{\mathcal{E}}_{t})_{\upsilon\gamma}=\partial_{\gamma}\bar{\tilde{b}}^{\upsilon}(S_{t}). By claim The Ising Population Model Instantiates the Data of the Fluctuation Theory §observation-extension, b~ˉυ(x)=qxυ+q0(x1+x2)\bar{\tilde{b}}^{\upsilon}(x)=q\,x^{\upsilon}+q_{0}(x^{1}+x^{2}), an affine function whose slice derivative in the coordinate xγx^{\gamma} is q1{υ=γ}+q0q\,\mathbf{1}_{\{\upsilon=\gamma\}}+q_{0}; hence E~t=qI+q0nn\tilde{\mathcal{E}}_{t}=q\,I+q_{0}\,\mathsf{n}\mathsf{n}^{\top}, the matrix nn\mathsf{n}\mathsf{n}^{\top} having all entries equal to 11.

The state noise covariance is Θt=Θ(St,At)\Theta^{\star}_{t}=\Theta(S_{t},A_{t}), which by claim The Ising Population Model Instantiates the Data of the Fluctuation Theory §covariance equals (121+121)vv=vv\bigl(\tfrac12\cdot1+\tfrac12\cdot1\bigr)vv^{\top}=vv^{\top}. The observation noise covariance has entries 1{υ=υ}b~υ(St)\mathbf{1}_{\{\upsilon=\upsilon'\}}\tilde{b}^{\upsilon}(S_{t}), and b~υ(St)=q12+q0\tilde{b}^{\upsilon}(S_{t})=q\cdot\tfrac12+q_{0} by claim The Ising Population Model Instantiates the Data of the Fluctuation Theory §observations; so Θ~t=(q2+q0)I\tilde{\Theta}^{\star}_{t}=\bigl(\tfrac{q}{2}+q_{0}\bigr)I. The terminal matrix is F=0F^{\star}=0 by claim 1, and the Hessian blocks are read off from claim 1: HtSSH^{SS}_{t} has entries HγδH_{\gamma\delta}, which is 2Qt2Q_{t} by claim 2; HtSAH^{SA}_{t} and HtASH^{AS}_{t} have entries Hγ,2+jH_{\gamma,2+j} and H2+j,γH_{2+j,\gamma}, all zero; and HtAAH^{AA}_{t} has entries H2+j,2+k=(2χ)11{j=k}H_{2+j,2+k}=(2\chi)^{-1}\mathbf{1}_{\{j=k\}}.

Claim 4. Both Rt=(4χ)1IR_{t}=(4\chi)^{-1}I and Θ~t=cI\tilde{\Theta}^{\star}_{t}=cI with c=q2+q0>0c=\tfrac{q}{2}+q_{0}>0 are symmetric, and for z0z\neq0 we have zRtz=(4χ)1z2>0z\cdot R_{t}z=(4\chi)^{-1}|z|^{2}>0 and zΘ~tz=cz2>0z\cdot\tilde{\Theta}^{\star}_{t}z=c|z|^{2}>0; so both are positive definite and hence invertible by Invertibility of Symmetric Positive Definite Matrices. Their inverses are 4χI4\chi I and c1Ic^{-1}I, since the products with RtR_{t} and Θ~t\tilde{\Theta}^{\star}_{t} are II.

Write N=nnN=\mathsf{n}\mathsf{n}^{\top}. Then N2=n(nn)n=2NN^{2}=\mathsf{n}(\mathsf{n}^{\top}\mathsf{n})\mathsf{n}^{\top}=2N, because nn=2\mathsf{n}^{\top}\mathsf{n}=2. The matrix E~t=qI+q0N\tilde{\mathcal{E}}_{t}=qI+q_{0}N is symmetric, so

D~t=E~t(Θ~t)1E~t=c1(qI+q0N)2=c1(q2I+2qq0N+q022N)=c1(q2I+2q0(q+q0)N),\tilde{D}_{t}=\tilde{\mathcal{E}}_{t}^{\top}(\tilde{\Theta}^{\star}_{t})^{-1}\tilde{\mathcal{E}}_{t}=c^{-1}\bigl(qI+q_{0}N\bigr)^{2}=c^{-1}\bigl(q^{2}I+2qq_{0}N+q_{0}^{2}\cdot2N\bigr)=c^{-1}\bigl(q^{2}I+2q_{0}(q+q_{0})N\bigr),

which is the displayed formula. Since nv=11+1(1)=0\mathsf{n}^{\top}v=1\cdot1+1\cdot(-1)=0, we get Nvv=n(nv)v=0N\,vv^{\top}=\mathsf{n}(\mathsf{n}^{\top}v)v^{\top}=0 and likewise vvN=0vv^{\top}N=0; hence D~tvv=c1q2vv=d~vv\tilde{D}_{t}vv^{\top}=c^{-1}q^{2}vv^{\top}=\tilde{d}\,vv^{\top} and vvD~t=d~vvvv^{\top}\tilde{D}_{t}=\tilde{d}\,vv^{\top}. Finally d~=q2/c>0\tilde{d}=q^{2}/c>0 because q>0q>0, and multiplying numerator and denominator by 22 gives d~=2q2/(q+2q0)\tilde{d}=2q^{2}/(q+2q_{0}).

Claim 5. Using the values of claim 1 and the symmetry of HH,

12i=14j=14Hijwiwj=ψ+2μ2((w1)2+(w2)2)+(2μψ)w1w2+14χ((w3)2+(w4)2).\frac{1}{2}\sum_{i=1}^{4}\sum_{j=1}^{4}H_{ij}\,w^{i}w^{j}=\frac{\psi+2\mu}{2}\Bigl(\bigl(w^{1}\bigr)^{2}+\bigl(w^{2}\bigr)^{2}\Bigr)+(2\mu-\psi)\,w^{1}w^{2}+\frac{1}{4\chi}\Bigl(\bigl(w^{3}\bigr)^{2}+\bigl(w^{4}\bigr)^{2}\Bigr).

On the other hand, expanding the squares,

ψ2(w1w2)2+μ(w1+w2)2=(ψ2+μ)((w1)2+(w2)2)+(2μψ)w1w2,\frac{\psi}{2}\bigl(w^{1}-w^{2}\bigr)^{2}+\mu\bigl(w^{1}+w^{2}\bigr)^{2}=\Bigl(\frac{\psi}{2}+\mu\Bigr)\Bigl(\bigl(w^{1}\bigr)^{2}+\bigl(w^{2}\bigr)^{2}\Bigr)+(2\mu-\psi)\,w^{1}w^{2},

which proves the displayed identity.

Put s=w1w2s=w^{1}-w^{2} and n=w1+w2n=w^{1}+w^{2}, so that s2+n2=2((w1)2+(w2)2)s^{2}+n^{2}=2\bigl((w^{1})^{2}+(w^{2})^{2}\bigr). Then

ψ2s2+μn2  min{ψ2,μ}(s2+n2)=min{ψ,2μ}((w1)2+(w2)2),\frac{\psi}{2}s^{2}+\mu n^{2}\ \ge\ \min\Bigl\{\frac{\psi}{2},\mu\Bigr\}\bigl(s^{2}+n^{2}\bigr)=\min\{\psi,2\mu\}\Bigl(\bigl(w^{1}\bigr)^{2}+\bigl(w^{2}\bigr)^{2}\Bigr),

so the whole expression is at least cJ((w1)2+(w2)2+(w3)2+(w4)2)=cJw2c_{J}\bigl((w^{1})^{2}+(w^{2})^{2}+(w^{3})^{2}+(w^{4})^{2}\bigr)=c_{J}|w|^{2} with cJ=min{ψ,2μ,(4χ)1}c_{J}=\min\{\psi,2\mu,(4\chi)^{-1}\}, a positive number because ψ>0\psi>0, μ>0\mu>0 and χ>0\chi>0. Since Fγδ=0F_{\gamma\delta}=0, the terminal quadratic form is identically 00, hence nonnegative; the two displayed requirements of (JC) are therefore met.

For (H1), aRta=(4χ)1a2a\cdot R_{t}a=(4\chi)^{-1}|a|^{2} for every aR2a\in\mathbb{R}^{2} and every tt, so the hypothesis holds with r=(4χ)1r=(4\chi)^{-1}, and a fortiori with the smaller constant r=cJr=c_{J}.

Please log in to copy this version.

Citations

Loading…

Dependency Graph

0 prerequisites

Prerequisites

Loading...

Comments

Loading…