Local construction and the third-order input. Suppose the support is contained
in a ball of radius R R R centered at the origin. For every finite pair ( c , Q ) (c,Q) ( c , Q ) ,
the tilted law has the same support as μ \mu μ . Its covariance is positive definite
because isotropy makes that support full dimensional. Compact support permits
differentiation of the partition function to every order. Thus its mean and
covariance, the inverse covariance and its positive square root are smooth locally
in ( c , Q ) (c,Q) ( c , Q ) . The SDE consequently has a unique strong solution up to its explosion
time. Along this solution Q t Q_t Q t is positive semidefinite, so μ t \mu_t μ t is log-concave.
The stated Letwin inputs give, for any isotropic log-concave ν \nu ν and any v v v ,
∣ κ 3 ν ( v ) ∣ 2 ≤ 8 ∣ v ∣ 2 . |\kappa_3^\nu(v)|^2\le8|v|^2. ∣ κ 3 ν ( v ) ∣ 2 ≤ 8∣ v ∣ 2 . Indeed, for symmetric M M M , the covariance of ⟨ v , X ⟩ \langle v,X\rangle ⟨ v , X ⟩ with
X T M X X^TMX X T MX is ⟨ κ 3 ν ( v ) , M ⟩ \langle\kappa_3^\nu(v),M\rangle ⟨ κ 3 ν ( v ) , M ⟩ ; Cauchy–Schwarz and the
quadratic-chaos variance bound give its square at most 8 ∣ v ∣ 2 ∣ M ∣ H S 2 8|v|^2|M|_{\rm HS}^2 8∣ v ∣ 2 ∣ M ∣ HS 2 .
Take the supremum over symmetric matrices of Hilbert–Schmidt norm one.
The density and covariance equations. All calculations initially take place
before an exit time from a compact subset of the finite ( c , Q ) (c,Q) ( c , Q ) parameter space.
For g t ( x ) = exp ( ⟨ c t , x ⟩ − ⟨ Q t x , x ⟩ / 2 ) g_t(x)=\exp(\langle c_t,x\rangle-\langle Q_tx,x\rangle/2) g t ( x ) = exp (⟨ c t , x ⟩ − ⟨ Q t x , x ⟩ /2 ) , Itô’s formula gives
d g t ( x ) g t ( x ) = ⟨ x , A t − 1 a t ⟩ d t + ⟨ x , A t − 1 / 2 d B t ⟩ . \frac{dg_t(x)}{g_t(x)}=\langle x,A_t^{-1}a_t\rangle dt
+\langle x,A_t^{-1/2}dB_t\rangle. g t ( x ) d g t ( x ) = ⟨ x , A t − 1 a t ⟩ d t + ⟨ x , A t − 1/2 d B t ⟩ . Integrating gives
d Z t / Z t = ⟨ a t , A t − 1 a t ⟩ d t + ⟨ a t , A t − 1 / 2 d B t ⟩ dZ_t/Z_t=\langle a_t,A_t^{-1}a_t\rangle dt+
\langle a_t,A_t^{-1/2}dB_t\rangle d Z t / Z t = ⟨ a t , A t − 1 a t ⟩ d t + ⟨ a t , A t − 1/2 d B t ⟩ . The quotient rule, including its cross
variation, cancels the drift and yields
d p t ( x ) = p t ( x ) ⟨ A t − 1 / 2 ( x − a t ) , d B t ⟩ , p t = g t / Z t . dp_t(x)=p_t(x)\langle A_t^{-1/2}(x-a_t),dB_t\rangle,
\qquad p_t=g_t/Z_t. d p t ( x ) = p t ( x ) ⟨ A t − 1/2 ( x − a t ) , d B t ⟩ , p t = g t / Z t . Consequently, for any bounded measurable f f f on the support,
d ∫ f d μ t = ⟨ A t − 1 / 2 ∫ f ( x ) ( x − a t ) d μ t ( x ) , d B t ⟩ . d\int f\,d\mu_t=
\left\langle A_t^{-1/2}\int f(x)(x-a_t)\,d\mu_t(x),dB_t\right\rangle. d ∫ f d μ t = ⟨ A t − 1/2 ∫ f ( x ) ( x − a t ) d μ t ( x ) , d B t ⟩ . Apply this identity to x x x and x x T xx^T x x T , and subtract d ( a t a t T ) d(a_ta_t^T) d ( a t a t T ) .
The barycenter noise is A t 1 / 2 d B t A_t^{1/2}dB_t A t 1/2 d B t , whose quadratic variation is A t d t A_tdt A t d t .
The centered third moment is the remaining covariance noise. This proves the two
claimed SDEs, initially up to explosion.
Let ν t \nu_t ν t be the image of μ t \mu_t μ t under x ↦ A t − 1 / 2 ( x − a t ) x\mapsto A_t^{-1/2}(x-a_t) x ↦ A t − 1/2 ( x − a t ) .
For m ≥ 2 m\ge2 m ≥ 2 , affine covariance of the derivatives of the logarithmic Laplace
transform gives
∣ κ m ( t ) ( v ) ∣ t 2 = ∣ κ m ν t ( A t 1 / 2 v ) ∣ 2 . |\kappa_m(t)(v)|_t^2
=|\kappa_m^{\nu_t}(A_t^{1/2}v)|^2. ∣ κ m ( t ) ( v ) ∣ t 2 = ∣ κ m ν t ( A t 1/2 v ) ∣ 2 . In particular, the third-order estimate gives E 3 ≤ 8 ⟨ A t u , u ⟩ \mathcal E_3\le8\langle A_tu,u\rangle E 3 ≤ 8 ⟨ A t u , u ⟩ .
The quadratic form of ∑ i H i A t − 1 H i \sum_iH_iA_t^{-1}H_i ∑ i H i A t − 1 H i at v v v is exactly
∣ κ 3 ( t ) ( v ) ∣ t 2 |\kappa_3(t)(v)|_t^2 ∣ κ 3 ( t ) ( v ) ∣ t 2 . This proves the matrix inequality.
No finite-time explosion. Set G i = A t − 1 / 2 H i A t − 1 / 2 G_i=A_t^{-1/2}H_iA_t^{-1/2} G i = A t − 1/2 H i A t − 1/2 . The preceding
inequality implies ∑ i G i 2 ≤ 8 I \sum_iG_i^2\le8I ∑ i G i 2 ≤ 8 I . Applying Itô to log det A t \log\det A_t log det A t gives
d log det A t = ∑ i Tr ( G i ) d B t , i − ( n + 1 2 ∑ i Tr ( G i 2 ) ) d t . d\log\det A_t=\sum_i\operatorname{Tr}(G_i)dB_{t,i}
-\left(n+\frac12\sum_i\operatorname{Tr}(G_i^2)\right)dt. d log det A t = i ∑ Tr ( G i ) d B t , i − ( n + 2 1 i ∑ Tr ( G i 2 ) ) d t . Its drift lies between − 5 n -5n − 5 n and − n -n − n , and the bracket of its local martingale
is at most 8 n 2 d t 8n^2dt 8 n 2 d t , since ( Tr G i ) 2 ≤ n Tr ( G i 2 ) (\operatorname{Tr}G_i)^2\le n\operatorname{Tr}(G_i^2) ( Tr G i ) 2 ≤ n Tr ( G i 2 ) .
These bounds hold uniformly under stopping. On any deterministic interval [ 0 , T ] [0,T] [ 0 , T ]
the stopped stochastic integral extends continuously to the explosion time:
extend its integrand by zero after that time, and use its bounded quadratic
variation. Hence log det A t \log\det A_t log det A t is bounded below before T T T on each path.
Compact support gives A t ≤ R 2 I A_t\le R^2I A t ≤ R 2 I . The determinant lower bound therefore bounds
A t − 1 A_t^{-1} A t − 1 on that path before T T T . It also bounds the coefficients of c t , Q t c_t,Q_t c t , Q t ,
because ∣ a t ∣ ≤ R |a_t|\le R ∣ a t ∣ ≤ R . Their drift integrals and stochastic integrals have finite
limits. Local existence at these finite limiting parameters extends the solution,
contradicting a finite explosion time. The solution is therefore global.
Expectation of the covariance. The process e t A t e^tA_t e t A t is a local martingale by
its SDE, and all its entries are bounded by e T R 2 e^TR^2 e T R 2 on [ 0 , T ] [0,T] [ 0 , T ] . It is a true
martingale. Since A 0 = I A_0=I A 0 = I , its expectation is I I I . This proves
E A t = e − t I \mathbb EA_t=e^{-t}I E A t = e − t I and, in particular,
E E 2 ( t ) = e − t \mathbb E\mathcal E_2(t)=e^{-t} E E 2 ( t ) = e − t and
E E 3 ( t ) ≤ 8 e − t \mathbb E\mathcal E_3(t)\le8e^{-t} E E 3 ( t ) ≤ 8 e − t .
The cumulant equation. Write
M t ( z ) = ∫ e ⟨ z , x ⟩ d μ t ( x ) M_t(z)=\int e^{\langle z,x\rangle}\,d\mu_t(x) M t ( z ) = ∫ e ⟨ z , x ⟩ d μ t ( x ) ,
Λ t = log M t \Lambda_t=\log M_t Λ t = log M t , and v t ( z ) = ∇ Λ t ( z ) − a t v_t(z)=\nabla\Lambda_t(z)-a_t v t ( z ) = ∇ Λ t ( z ) − a t .
The bounded support and stopped parameter coefficients justify differentiation
and stochastic integration in either order, on compact sets of z z z . The density
equation and scalar Itô formula give
d Λ t ( z ) = ⟨ A t − 1 / 2 v t ( z ) , d B t ⟩ − 1 2 ⟨ v t ( z ) , A t − 1 v t ( z ) ⟩ d t . d\Lambda_t(z)=\langle A_t^{-1/2}v_t(z),dB_t\rangle
-\tfrac12\langle v_t(z),A_t^{-1}v_t(z)\rangle dt. d Λ t ( z ) = ⟨ A t − 1/2 v t ( z ) , d B t ⟩ − 2 1 ⟨ v t ( z ) , A t − 1 v t ( z )⟩ d t . At zero, v t ( 0 ) = 0 v_t(0)=0 v t ( 0 ) = 0 and its r r r th derivative, contracted against a further
vector, is κ r + 1 ( t ) \kappa_{r+1}(t) κ r + 1 ( t ) . The m m m th derivative of the noise is the asserted
( m + 1 ) (m+1) ( m + 1 ) st cumulant. Leibniz’s rule for the quadratic drift is the sum over
nonempty proper subsets J ⊂ [ m ] J\subset[m] J ⊂ [ m ] . Its 2 m 2m 2 m singleton or co-singleton terms
each equal κ m ( t ) \kappa_m(t) κ m ( t ) , because one factor is A t h j A_th_j A t h j and cancels A t − 1 A_t^{-1} A t − 1 .
The remaining terms occur in equal complementary pairs. Selecting the member
containing the first index and accounting for the prefactor 1 / 2 1/2 1/2 gives precisely
L m ( t ) L_m(t) L m ( t ) . This proves the equation, including L 3 = 0 L_3=0 L 3 = 0 .
Crude moment bounds, independent of the induction. For completeness there are
finite constants D n , m D_{n,m} D n , m such that every isotropic log-concave ν \nu ν on
R n \mathbb R^n R n satisfies ∣ κ m ν ( v ) ∣ 2 ≤ D n , m ∣ v ∣ 2 |\kappa_m^\nu(v)|^2\le D_{n,m}|v|^2 ∣ κ m ν ( v ) ∣ 2 ≤ D n , m ∣ v ∣ 2 . No dimension-free
claim is made here. To see this, a centered variance-one one-dimensional
log-concave marginal Y Y Y has a log-concave survival function S ( t ) S(t) S ( t ) .
This last fact follows either by integration of a log-concave density or directly
from its increasing hazard: for density f f f , the ratio
S ( t ) / f ( t ) = ∫ 0 ∞ f ( t + s ) / f ( t ) d s S(t)/f(t)=\int_0^\infty f(t+s)/f(t)\,ds S ( t ) / f ( t ) = ∫ 0 ∞ f ( t + s ) / f ( t ) d s is nonincreasing, by concavity of
log f \log f log f . Endpoints follow by limits. Chebyshev gives S ( 2 ) ≤ 1 / 4 S(2)\le1/4 S ( 2 ) ≤ 1/4 and
S ( − 2 ) ≥ 3 / 4 S(-2)\ge3/4 S ( − 2 ) ≥ 3/4 . Concavity of log S \log S log S therefore gives, for t ≥ 2 t\ge2 t ≥ 2 ,
S ( t ) ≤ 1 4 exp ( − ( log 3 ) ( t − 2 ) / 4 ) . S(t)\le\tfrac14\exp\bigl(-(\log3)(t-2)/4\bigr). S ( t ) ≤ 4 1 exp ( − ( log 3 ) ( t − 2 ) /4 ) . If S ( 2 ) = 0 S(2)=0 S ( 2 ) = 0 , its upper tail vanishes and the displayed bound is immediate.
The same argument applies to − Y -Y − Y . Integration of the two tails gives a finite
universal upper bound for E ∣ Y ∣ p \mathbb E|Y|^p E ∣ Y ∣ p at every fixed p p p .
For unit vectors v 1 , … , v m v_1,\ldots,v_m v 1 , … , v m , Hölder bounds every mixed moment
E ∏ j ∈ J ⟨ v j , X ⟩ \mathbb E\prod_{j\in J}\langle v_j,X\rangle E ∏ j ∈ J ⟨ v j , X ⟩ in terms of ∣ J ∣ |J| ∣ J ∣ alone.
The moment-cumulant formula
κ m ( v 1 , … , v m ) = ∑ π ∈ P ( [ m ] ) ( − 1 ) ∣ π ∣ − 1 ( ∣ π ∣ − 1 ) ! ∏ J ∈ π E ∏ j ∈ J ⟨ v j , X ⟩ \kappa_m(v_1,\ldots,v_m)=
\sum_{\pi\in\mathcal P([m])}(-1)^{|\pi|-1}(|\pi|-1)!
\prod_{J\in\pi}\mathbb E\prod_{j\in J}\langle v_j,X\rangle κ m ( v 1 , … , v m ) = π ∈ P ([ m ]) ∑ ( − 1 ) ∣ π ∣ − 1 ( ∣ π ∣ − 1 )! J ∈ π ∏ E j ∈ J ∏ ⟨ v j , X ⟩ then bounds each mixed cumulant by a finite C m C_m C m . Summing its square over the
n m − 1 n^{m-1} n m − 1 coordinate choices gives D n , m = n m − 1 C m 2 D_{n,m}=n^{m-1}C_m^2 D n , m = n m − 1 C m 2 .
Whitening the posterior now gives
E m ( t ) ≤ D n , m ⟨ A t u , u ⟩ ≤ D n , m R 2 . \mathcal E_m(t)\le D_{n,m}\langle A_tu,u\rangle\le D_{n,m}R^2. E m ( t ) ≤ D n , m ⟨ A t u , u ⟩ ≤ D n , m R 2 . For a summand of L m L_m L m with ∣ J ∣ = a |J|=a ∣ J ∣ = a , sum first over the whitened coordinates in
J c J^c J c , apply the same bound for κ m − a + 1 \kappa_{m-a+1} κ m − a + 1 , and then sum over those in
J ∖ { 1 } J\setminus\{1\} J ∖ { 1 } . Its squared weighted norm is at most
D n , m − a + 1 E a + 1 ( t ) D_{n,m-a+1}\mathcal E_{a+1}(t) D n , m − a + 1 E a + 1 ( t ) : the contracted vector is
A t − 1 κ a + 1 ( t ) ( u , h J ∖ { 1 } , ⋅ ) A_t^{-1}\kappa_{a+1}(t)(u,h_{J\setminus\{1\}},\cdot) A t − 1 κ a + 1 ( t ) ( u , h J ∖ { 1 } , ⋅ ) , whose A t A_t A t -norm
is exactly the sum over the final whitened contraction index. The finite sum
over J J J therefore bounds ∣ L m ( t ) ( u ) ∣ t 2 |L_m(t)(u)|_t^2 ∣ L m ( t ) ( u ) ∣ t 2 deterministically as well.
Finally write r = m − 1 r=m-1 r = m − 1 , T = κ m ( t ) ( u ) T=\kappa_m(t)(u) T = κ m ( t ) ( u ) ,
W i = κ m + 1 ( t ) ( u , ⋅ , … , ⋅ , A t − 1 / 2 e i ) W_i=\kappa_{m+1}(t)(u,\cdot,\ldots,\cdot,A_t^{-1/2}e_i) W i = κ m + 1 ( t ) ( u , ⋅ , … , ⋅ , A t − 1/2 e i ) ,
T ~ = ( A t − 1 / 2 ) ⊗ r T \widetilde T=(A_t^{-1/2})^{\otimes r}T T = ( A t − 1/2 ) ⊗ r T , and
W ~ i = ( A t − 1 / 2 ) ⊗ r W i \widetilde W_i=(A_t^{-1/2})^{\otimes r}W_i W i = ( A t − 1/2 ) ⊗ r W i .
The noise coefficient of E m \mathcal E_m E m is
β i = 2 ⟨ T ~ , W ~ i ⟩ − ⟨ T ~ , ∑ s = 1 r G i ( s ) T ~ ⟩ , \beta_i=2\langle\widetilde T,\widetilde W_i\rangle
-\left\langle\widetilde T,\sum_{s=1}^rG_i^{(s)}\widetilde T\right\rangle, β i = 2 ⟨ T , W i ⟩ − ⟨ T , s = 1 ∑ r G i ( s ) T ⟩ , where G i ( s ) G_i^{(s)} G i ( s ) acts on tensor slot s s s . Since ∑ i G i 2 ≤ 8 I \sum_iG_i^2\le8I ∑ i G i 2 ≤ 8 I ,
∑ i β i 2 ≤ 8 E m E m + 1 + 16 r 2 E m 2 . \sum_i\beta_i^2\le8\mathcal E_m\mathcal E_{m+1}
+16r^2\mathcal E_m^2. i ∑ β i 2 ≤ 8 E m E m + 1 + 16 r 2 E m 2 . This is bounded deterministically. The drift is integrable on finite intervals
as well: the second derivative of the inverse metric consists of finitely many
terms bounded by E m \mathcal E_m E m , E m + 1 \mathcal E_{m+1} E m + 1 and
E m ∣ L m ( u ) ∣ t \sqrt{\mathcal E_m}|L_m(u)|_t E m ∣ L m ( u ) ∣ t , using the same matrix inequality.
Thus localization may be removed in expectations. These observations establish
the final integrability assertion and supply the precise justification for the
energy argument of Lemma 8.4 . All assertions concerning a
fixed vector extend from unit u u u to arbitrary u u u by homogeneity; each energy
and squared norm acquires a factor ∣ u ∣ 2 |u|^2 ∣ u ∣ 2 . In particular,
∫ 0 ∞ E E 2 ( t ) d t = ∣ u ∣ 2 \int_0^\infty\mathbb E\mathcal E_2(t)\,dt=|u|^2 ∫ 0 ∞ E E 2 ( t ) d t = ∣ u ∣ 2 .