Exact constants on the line and on products, and the bound $4$ on Dirichlet laws
Part of the moment-map mechanism, Chapter The moment map: exact cases ; the reading order is on the full proofs page.
Overview. This dossier proves C M H ( 4 ) \mathrm{CMH}(4) CMH ( 4 ) in three classes where the constant of Definition 16.1 can be computed: the line (Theorem 17.1 ), products (Theorem 17.2 , Corollary 17.1 ) and the log-concave Dirichlet family (Theorem 17.3 , with Lemma 17.1 , Lemma 17.2 , Lemma 17.3 , Corollary 17.2 , Corollary 17.3 ). It also records Corollary 17.5 . The Dirichlet case is proved by lifting to independent Gamma variables, completing a square in the Gamma Bochner identity, and using the Euler constraint of degree-0 homogeneity. Universal C M H ( 4 ) \mathrm{CMH}(4) CMH ( 4 ) remains open. Products of centered one-sided exponentials saturate the constant 4 exactly.
On the line, C C M H = C P / Var \CMH=\CP/\Var C CMH = C P / Var exactly, the KLS bound gives C P ≤ 4 Var \CP\le4\Var C P ≤ 4 Var , and the centered exponential attains 4 (Theorem D29.1 ).
For products, the cross terms E [ ( L i g ) ( L j g ) ] \E[(L_ig)(L_jg)] E [( L i g ) ( L j g )] are nonnegative, so C C M H \CMH C CMH is the maximum over the factors (Theorem D29.2 ). Affine Poincaré passes to linear images (Corollary D29.1 ).
The Dirichlet data and the tangent pseudo-inverse (Lemma D29.1 ) reduce the claim to A ( A + 1 ) d α ≤ 4 n α A(A+1)d_\alpha\le4n_\alpha A ( A + 1 ) d α ≤ 4 n α .
Gamma lift: the Bochner identity (D29.14) and a square completion give N Γ − 1 4 D Γ = ∑ i E R i + ∑ i δ i D i N_\Gamma-\tfrac14D_\Gamma=\sum_i\E R_i+\sum_i\delta_iD_i N Γ − 4 1 D Γ = ∑ i E R i + ∑ i δ i D i (Lemma D29.2 ). The Euler constraint bounds each R i R_i R i from below (Lemma D29.3 ).
Using S ⊥ P S\perp P S ⊥ P and the inverse moments of S S S , step 4 becomes a pointwise coefficient F A ( α i , P i ) F_A(\alpha_i,P_i) F A ( α i , P i ) . A scalar minimization bounds it by A − 1 2 + s A A-\tfrac12+s_A A − 2 1 + s A with s A ≥ 0 s_A\ge0 s A ≥ 0 (Lemma D29.4 ), which proves Theorem D29.3 , with strict surplus for A > 3 A>3 A > 3 .
Theorem 16.1 with step 2 gives C P a f f ≤ 4 \CPaff\le4 C P aff ≤ 4 for Dirichlet laws, their products, linear images and convolutions.
Scope. This dossier proves C M H ( 4 ) \mathrm{CMH}(4) CMH ( 4 ) in the three classes where the constant of Definition 16.1 is computable, the third being the first genuinely nonproduct family for which Route C has a theorem. Combined with Theorem 16.1 (dossier solutions/thm-cmh-normalization.md) each result yields an affine Poincaré bound with the same constant. Nothing here bears on universal C M H ( 4 ) \mathrm{CMH}(4) CMH ( 4 ) , which remains open; Corollary 17.5 in the manuscript explains why the product case in particular is a warning rather than encouragement.
Notation is that of §The Stein generator and the CMH constant : H H H is the moment-map Stein kernel in target coordinates, L μ = div μ ( H ∇ ⋅ ) L_\mu=\Div_\mu(H\nabla\,\cdot\,) L μ = div μ ( H ∇ ⋅ ) , A = − L μ \Aop=-L_\mu A = − L μ , and C C M H \CMH C CMH is (16.5) .
1. The line ¶ Let μ ( d x ) = ρ d x \mu(\dd x)=\rho\dd x μ ( d x ) = ρ d x be centered on ( ℓ , r ) (\ell,r) ( ℓ , r ) with variance σ 2 \sigma^2 σ 2 and Stein kernel τ \tau τ , the zero-flux solution of ( τ ρ ) ′ = − x ρ (\tau\rho)'=-x\rho ( τ ρ ) ′ = − x ρ .
We use the natural no-flux realizations of §Conventions, domains, and affine covariance .
Let D D D be the closed maximal derivative in L 2 ( μ ) L^2(\mu) L 2 ( μ ) , with domain
H 1 ( μ ) H^1(\mu) H 1 ( μ ) , and let A \Aop A be the operator of the natural weighted form
E τ ( f , q ) = ∫ τ f ′ q ′ d μ \mathcal E_\tau(f,q)=\int\tau f'q'\,d\mu E τ ( f , q ) = ∫ τ f ′ q ′ d μ . Constants and bounded
locally absolutely continuous functions whose derivative has compact interior
support belong to these form domains when their displayed energies are finite.
The local density regularity is that of the differential-operator setup
(in particular, locally positive continuous densities suffice, using weak
derivatives); no endpoint regularity, log-concavity, or spectral gap is used.
Both constants in the asserted identity are allowed to be infinite, and the
CMH numerator is understood as an extended nonnegative integral.
Put p = τ ρ p=\tau\rho p = τ ρ . Centering and nonzero finite variance give p > 0 p>0 p > 0 in the
interior. The following graph core encodes the adjoint’s no-flux condition:
U = { u = ϕ / ρ : ϕ ∈ C c ∞ ( ( ℓ , r ) ) } , D μ ∗ u = − ϕ ′ / ρ . \mathscr U=\{u=\phi/\rho:\phi\in C_c^\infty((\ell,r))\},
\qquad D_\mu^*u=-\phi'/\rho. U = { u = ϕ / ρ : ϕ ∈ C c ∞ (( ℓ , r ))} , D μ ∗ u = − ϕ ′ / ρ . To verify the core assertion, define B 0 B_0 B 0 by this expression on U \mathscr U U .
Its domain is dense: compact localization and smooth approximation of the
flux ρ u \rho u ρ u are dense in L 2 ( ρ − 1 d x ) L^2(\rho^{-1}dx) L 2 ( ρ − 1 d x ) . The adjoint condition,
tested against ϕ \phi ϕ , is precisely that a scalar function have a
distributional derivative in L 2 ( μ ) L^2(\mu) L 2 ( μ ) . Thus B 0 ∗ = D B_0^*=D B 0 ∗ = D , with no
boundary restriction on the scalar function, and
B 0 ‾ = D μ ∗ \overline{B_0}=D_\mu^* B 0 = D μ ∗ . Also ker D \ker D ker D is exactly the constants.
Every field of this core is attained by an actual CMH test. Fix an interior
x 0 x_0 x 0 and set
f ϕ ( x ) = ∫ x 0 x ϕ ( t ) p ( t ) d t , h ϕ = − ϕ ′ ρ . f_\phi(x)=\int_{x_0}^x\frac{\phi(t)}{p(t)}\,dt,
\qquad h_\phi=-\frac{\phi'}\rho. f ϕ ( x ) = ∫ x 0 x p ( t ) ϕ ( t ) d t , h ϕ = − ρ ϕ ′ . Then f ϕ f_\phi f ϕ is bounded and constant near the endpoints, and
E τ ( f ϕ , f ϕ ) = ∫ ϕ 2 / p d x < ∞ \mathcal E_\tau(f_\phi,f_\phi)=\int\phi^2/p\,dx<\infty E τ ( f ϕ , f ϕ ) = ∫ ϕ 2 / p d x < ∞ .
For every weighted-form test q q q , compactly supported integration by parts gives
E τ ( f ϕ , q ) = ∫ ϕ q ′ d x = − ∫ ϕ ′ q d x = ⟨ h ϕ , q ⟩ . \mathcal E_\tau(f_\phi,q)=\int\phi q'\,dx
=-\int\phi' q\,dx=\inner{h_\phi}{q}. E τ ( f ϕ , q ) = ∫ ϕ q ′ d x = − ∫ ϕ ′ q d x = ⟨ h ϕ , q ⟩ . Since h ϕ ∈ L 2 ( μ ) h_\phi\in L^2(\mu) h ϕ ∈ L 2 ( μ ) , the operator-domain criterion gives
f ϕ ∈ Dom ( A ) f_\phi\in\Dom(\Aop) f ϕ ∈ Dom ( A ) , A f ϕ = h ϕ \Aop f_\phi=h_\phi A f ϕ = h ϕ , and
τ f ϕ ′ = u \tau f_\phi'=u τ f ϕ ′ = u . Moreover h ϕ h_\phi h ϕ is centered. No replacement of
Ran A \operatorname{Ran}\Aop Ran A by all of centered L 2 L^2 L 2 has been made.
Lower bound, including the infinite case. Put M = σ 2 C C M H ( μ ) M=\sigma^2\CMH(\mu) M = σ 2 C CMH ( μ ) .
If M < ∞ M<\infty M < ∞ , these attained tests give
∥ u ∥ 2 2 ≤ M ∥ D μ ∗ u ∥ 2 2 \norm u_2^2\le M\norm{D_\mu^*u}_2^2 ∥ u ∥ 2 2 ≤ M ∥ ∥ D μ ∗ u ∥ ∥ 2 2 on U \mathscr U U , and graph-core
closure extends it to all of Dom ( D μ ∗ ) \Dom(D_\mu^*) Dom ( D μ ∗ ) . This bound and closedness
make Ran D μ ∗ \operatorname{Ran}D_\mu^* Ran D μ ∗ closed: if D μ ∗ u j D_\mu^*u_j D μ ∗ u j converges,
then u j u_j u j is Cauchy and the closed graph supplies its limiting preimage.
Since the closure of this range is ( ker D ) ⊥ = L 0 2 ( μ ) (\ker D)^\perp=L^2_0(\mu) ( ker D ) ⊥ = L 0 2 ( μ ) , the
range equals L 0 2 ( μ ) L^2_0(\mu) L 0 2 ( μ ) . For any centered v ∈ Dom D v\in\Dom D v ∈ Dom D , choose
u u u with D μ ∗ u = v D_\mu^*u=v D μ ∗ u = v . Then
∥ v ∥ 2 2 = ⟨ D v , u ⟩ ≤ M ∥ D v ∥ 2 ∥ v ∥ 2 . \norm v_2^2=\inner{Dv}{u}
\le\sqrt M\norm{Dv}_2\norm v_2. ∥ v ∥ 2 2 = ⟨ D v , u ⟩ ≤ M ∥ D v ∥ 2 ∥ v ∥ 2 . Thus C P ≤ M \CP\le M C P ≤ M . If M = ∞ M=\infty M = ∞ that inequality is automatic; in
particular C P = ∞ \CP=\infty C P = ∞ forces C C M H = ∞ \CMH=\infty C CMH = ∞ .
Upper bound when C P < ∞ \CP<\infty C P < ∞ . Given f ∈ Dom ( A ) f\in\Dom(\Aop) f ∈ Dom ( A ) , put
h = A f h=\Aop f h = A f . Constants belong to ker A \ker\Aop ker A , so h h h is centered.
Poincaré gives, for every q ∈ Dom D q\in\Dom D q ∈ Dom D ,
∣ ⟨ h , q ⟩ ∣ ≤ C P ∥ h ∥ 2 ∥ D q ∥ 2 . |\inner hq|\le\sqrt{\CP}\norm h_2\norm{Dq}_2. ∣ ⟨ h , q ⟩ ∣ ≤ C P ∥ h ∥ 2 ∥ D q ∥ 2 . Apply Riesz representation to the bounded functional
D q ↦ ⟨ h , q ⟩ Dq\mapsto\inner hq D q ↦ ⟨ h , q ⟩ on the closure of Ran D \operatorname{Ran}D Ran D .
There exists w ∈ L 2 ( μ ) w\in L^2(\mu) w ∈ L 2 ( μ ) such that
⟨ w , D q ⟩ = ⟨ h , q ⟩ ( q ∈ Dom D ) , ∥ w ∥ 2 ≤ C P ∥ h ∥ 2 . \inner w{Dq}=\inner hq\quad(q\in\Dom D),
\qquad \norm w_2\le\sqrt{\CP}\norm h_2. ⟨ w , D q ⟩ = ⟨ h , q ⟩ ( q ∈ Dom D ) , ∥ w ∥ 2 ≤ C P ∥ h ∥ 2 . For η ∈ C c ∞ ( ( ℓ , r ) ) \eta\in C_c^\infty((\ell,r)) η ∈ C c ∞ (( ℓ , r )) , its bounded primitive
q ( x ) = ∫ x 0 x η ( t ) d t q(x)=\int_{x_0}^x\eta(t)\,dt q ( x ) = ∫ x 0 x η ( t ) d t belongs to both scalar form domains.
The defining weak identity for A f = h \Aop f=h A f = h therefore yields
∫ τ f ′ η d μ = ⟨ h , q ⟩ = ∫ w η d μ . \int\tau f'\eta\,d\mu=\inner hq=\int w\eta\,d\mu. ∫ τ f ′ η d μ = ⟨ h , q ⟩ = ∫ w η d μ . The left integrand is locally integrable by weighted Cauchy--Schwarz.
Arbitrariness of η \eta η proves τ f ′ = w \tau f'=w τ f ′ = w almost everywhere, without
assuming square integrability of τ f ′ \tau f' τ f ′ in advance. Consequently
∥ τ f ′ ∥ 2 2 ≤ C P ∥ A f ∥ 2 2 \norm{\tau f'}_2^2\le\CP\norm{\Aop f}_2^2 ∥ τ f ′ ∥ 2 2 ≤ C P ∥ A f ∥ 2 2 for every CMH test.
Since H = τ H=\tau H = τ and Σ = σ 2 \Sigma=\sigma^2 Σ = σ 2 , this proves
σ 2 C C M H ≤ C P \sigma^2\CMH\le\CP σ 2 C CMH ≤ C P in the finite case. In the infinite case equality
was already forced by the lower bound. No operator inverse is used.
The one-dimensional specialization of the Kannan–Lovász–Simonovits bound Kannan et al. , 1995 , recorded for example in Cattiaux & Guillin, 2018, Eq. (2.25) with that attribution, is C P ( μ ) ≤ 4 Var ( μ ) \CP(\mu)\le4\Var(\mu) C P ( μ ) ≤ 4 Var ( μ ) .
For sharpness let X = Y − 1 X=Y-1 X = Y − 1 with Y ∼ E x p ( 1 ) Y\sim\mathrm{Exp}(1) Y ∼ Exp ( 1 ) . Then X X X is centered, has density e − ( x + 1 ) 1 x ≥ − 1 d x e^{-(x+1)}\one_{x\ge-1}\dd x e − ( x + 1 ) 1 x ≥ − 1 d x , and Var ( X ) = 1 \Var(X)=1 Var ( X ) = 1 . For 0 < a < 1 2 0<a<\tfrac12 0 < a < 2 1 set f a ( x ) = e a ( x + 1 ) f_a(x)=e^{a(x+1)} f a ( x ) = e a ( x + 1 ) . Direct integration gives
Var ( f a ( X ) ) E ∣ f a ′ ( X ) ∣ 2 = ( 1 − 2 a ) − 1 − ( 1 − a ) − 2 a 2 ( 1 − 2 a ) − 1 = 1 ( 1 − a ) 2 ⟶ 4. \frac{\Var(f_a(X))}{\E\abs{f_a'(X)}^2}
=\frac{(1-2a)^{-1}-(1-a)^{-2}}{a^2(1-2a)^{-1}}
=\frac1{(1-a)^2}\longrightarrow4. E ∣ f a ′ ( X ) ∣ 2 Var ( f a ( X )) = a 2 ( 1 − 2 a ) − 1 ( 1 − 2 a ) − 1 − ( 1 − a ) − 2 = ( 1 − a ) 2 1 ⟶ 4. The upper bound therefore gives C P ( X ) = 4 \CP(X)=4 C P ( X ) = 4 , and the identity already proved gives C C M H ( X ) = 4 \CMH(X)=4 C CMH ( X ) = 4 .
This is an identity, not an inequality: in one dimension the solenoidal channel of Proposition 16.2 is empty (a square-integrable divergence-free field with vanishing flux is zero), so CMH carries exactly the affine Poincaré information.
2. Products ¶ For μ = ⨂ i = 1 m μ i \mu=\bigotimes_{i=1}^m\mu_i μ = ⨂ i = 1 m μ i with centered factors for which the canonical CMH data are defined,
C C M H ( μ ) = max i C C M H ( μ i ) . \CMH(\mu)=\max_i\CMH(\mu_i). C CMH ( μ ) = i max C CMH ( μ i ) . For one-dimensional factors this equals max i C P ( μ i ) / Var ( μ i ) \max_i\CP(\mu_i)/\Var(\mu_i) max i C P ( μ i ) / Var ( μ i ) . Consequently block products of C M H ( 4 ) \mathrm{CMH}(4) CMH ( 4 ) factors, and their invertible affine images, satisfy C M H ( 4 ) \mathrm{CMH}(4) CMH ( 4 ) .
Here Σ = diag ( Σ 1 , … , Σ m ) \Sigma=\diag(\Sigma_1,\dots,\Sigma_m) Σ = diag ( Σ 1 , … , Σ m ) , H = diag ( H 1 , … , H m ) H=\diag(H_1,\dots,H_m) H = diag ( H 1 , … , H m ) , and L μ = ∑ i L i L_\mu=\sum_iL_i L μ = ∑ i L i with L i L_i L i acting in block i i i only, so the L i L_i L i are commuting nonpositive generators. For i ≠ j i\ne j i = j , two integrations by parts give
E [ ( L i g ) ( L j g ) ] = E Tr ( H i D i j 2 g H j D j i 2 g ) = E ∥ H i 1 / 2 D i j 2 g H j 1 / 2 ∥ H S 2 ≥ 0. \E[(L_ig)(L_jg)]
=\E\Tr\!\left(H_iD^2_{ij}g\,H_jD^2_{ji}g\right)
=\E\norm{H_i^{1/2}D^2_{ij}g\,H_j^{1/2}}_{\HS}^2\ge0. E [( L i g ) ( L j g )] = E Tr ( H i D ij 2 g H j D ji 2 g ) = E ∥ ∥ H i 1/2 D ij 2 g H j 1/2 ∥ ∥ HS 2 ≥ 0. Hence E ( L μ g ) 2 = ∑ i E ( L i g ) 2 + 2 ∑ i < j E [ ( L i g ) ( L j g ) ] ≥ ∑ i E ( L i g ) 2 \E(L_\mu g)^2=\sum_i\E(L_ig)^2+2\sum_{i<j}\E[(L_ig)(L_jg)]\ge\sum_i\E(L_ig)^2 E ( L μ g ) 2 = ∑ i E ( L i g ) 2 + 2 ∑ i < j E [( L i g ) ( L j g )] ≥ ∑ i E ( L i g ) 2 .
For the numerator, conditioning on the other blocks and using the defining inequality for C C M H ( μ i ) \CMH(\mu_i) C CMH ( μ i ) gives
E ⟨ H i ∇ i g , Σ i − 1 H i ∇ i g ⟩ ≤ C C M H ( μ i ) E ( L i g ) 2 . \E\inner{H_i\nabla_i g}{\Sigma_i^{-1}H_i\nabla_i g}
\le\CMH(\mu_i)\,\E(L_i g)^2. E ⟨ H i ∇ i g , Σ i − 1 H i ∇ i g ⟩ ≤ C CMH ( μ i ) E ( L i g ) 2 . Summing and using the previous display gives “≤ \le ≤ ”. Testing on g g g depending on one block gives “≥ \ge ≥ ”. The one-dimensional formula, including factors with infinite Poincaré or CMH constant, follows from the extended-constant identity Theorem D29.1 ; the affine-image consequence uses the invertible affine covariance of C C M H \CMH C CMH from §Conventions, domains, and affine covariance .
For Y = T X Y=TX Y = TX and smooth f f f , Var Y f = Var X ( f ∘ T ) \Var_Yf=\Var_X(f\circ T) Var Y f = Var X ( f ∘ T ) and E X ⟨ Cov ( X ) ∇ ( f ∘ T ) , ∇ ( f ∘ T ) ⟩ = E ⟨ T Cov ( X ) T ⊤ ∇ f , ∇ f ⟩ = E ⟨ Cov ( Y ) ∇ f , ∇ f ⟩ \E_X\inner{\Cov(X)\nabla(f\circ T)}{\nabla(f\circ T)} =\E\inner{T\Cov(X)T^\top\nabla f}{\nabla f}=\E\inner{\Cov(Y)\nabla f}{\nabla f} E X ⟨ Cov ( X ) ∇ ( f ∘ T ) , ∇ ( f ∘ T ) ⟩ = E ⟨ T Cov ( X ) T ⊤ ∇ f , ∇ f ⟩ = E ⟨ Cov ( Y ) ∇ f , ∇ f ⟩ . No invertibility is needed.
3. The Dirichlet family: setup ¶ Let P ∼ Dir ( α ) P\sim\Dir(\alpha) P ∼ Dir ( α ) , α i ≥ 1 \alpha_i\ge1 α i ≥ 1 , A = ∑ i α i A=\sum_i\alpha_i A = ∑ i α i , q = α / A q=\alpha/A q = α / A . With C ( p ) = diag ( p ) − p p ⊤ C(p)=\diag(p)-pp^\top C ( p ) = diag ( p ) − p p ⊤ and L α L_\alpha L α the Wright–Fisher generator (17.9) , the potential φ ( y ) = A log ∑ i e y i / A − q ⋅ y \varphi(y)=A\log\sum_ie^{y_i/A}-q\cdot y φ ( y ) = A log ∑ i e y i / A − q ⋅ y has ∇ φ = p − q \nabla\varphi=p-q ∇ φ = p − q and D 2 φ = C ( p ) / A D^2\varphi=C(p)/A D 2 φ = C ( p ) / A , and pushes e − φ e^{-\varphi} e − φ forward to the centered law P − q P-q P − q . Hence the canonical data are H = C ( p ) / A H=C(p)/A H = C ( p ) / A , L μ = L α / A L_\mu=L_\alpha/A L μ = L α / A , and Σ \Sigma Σ as in (17.10) .
Let x = A ( A + 1 ) diag ( α ) − 1 v x=A(A+1)\diag(\alpha)^{-1}v x = A ( A + 1 ) diag ( α ) − 1 v . Then Σ x = ( diag ( α ) − α α ⊤ / A ) diag ( α ) − 1 v = v − ( α / A ) ∑ i v i = v \Sigma x=\bigl(\diag(\alpha)-\alpha\alpha^\top/A\bigr)\diag(\alpha)^{-1}v =v-(\alpha/A)\sum_iv_i=v Σ x = ( diag ( α ) − α α ⊤ / A ) diag ( α ) − 1 v = v − ( α / A ) ∑ i v i = v . Since ker Σ = s p a n { 1 } \ker\Sigma=\mathrm{span}\{\one\} ker Σ = span { 1 } and v ⊤ 1 = 0 v^\top\one=0 v ⊤ 1 = 0 , the solution’s ambiguity pairs to zero, so v ⊤ Σ † v = v ⊤ x v^\top\Sigma^\dagger v=v^\top x v ⊤ Σ † v = v ⊤ x .
With u i = ( C ( P ) ∇ g ) i u_i=(C(P)\nabla g)_i u i = ( C ( P ) ∇ g ) i , d α d_\alpha d α , n α n_\alpha n α as in (17.12) , the statement C C M H ≤ 4 \CMH\le4 C CMH ≤ 4 becomes exactly A ( A + 1 ) d α ( g ) ≤ 4 n α ( g ) A(A+1)d_\alpha(g)\le4n_\alpha(g) A ( A + 1 ) d α ( g ) ≤ 4 n α ( g ) : the numerator is A − 2 ⋅ A ( A + 1 ) d α = A + 1 A d α A^{-2}\cdot A(A+1)d_\alpha=\tfrac{A+1}Ad_\alpha A − 2 ⋅ A ( A + 1 ) d α = A A + 1 d α and the denominator A − 2 n α A^{-2}n_\alpha A − 2 n α .
For every m ≥ 2 m\ge2 m ≥ 2 , all α i ≥ 1 \alpha_i\ge1 α i ≥ 1 , and all g g g in the core, A ( A + 1 ) d α ( g ) ≤ 4 n α ( g ) A(A+1)\,d_\alpha(g)\le4\,n_\alpha(g) A ( A + 1 ) d α ( g ) ≤ 4 n α ( g ) .
4. The Gamma lift ¶ Let Y i ∼ Gamma ( α i , 1 ) Y_i\sim\GammaLaw(\alpha_i,1) Y i ∼ Gamma ( α i , 1 ) independent, S = ∑ i Y i S=\sum_iY_i S = ∑ i Y i , P = Y / S P=Y/S P = Y / S ; then S ∼ Gamma ( A , 1 ) S\sim\GammaLaw(A,1) S ∼ Gamma ( A , 1 ) , P ∼ Dir ( α ) P\sim\Dir(\alpha) P ∼ Dir ( α ) , and S ⊥ P S\perp P S ⊥ P . Set G ( Y ) = g ( Y / S ) G(Y)=g(Y/S) G ( Y ) = g ( Y / S ) , homogeneous of degree 0, and let L Γ G = ∑ i ( Y i G i i + ( α i − Y i ) G i ) \calL_\Gamma G=\sum_i(Y_iG_{ii}+(\alpha_i-Y_i)G_i) L Γ G = ∑ i ( Y i G ii + ( α i − Y i ) G i ) be the Laguerre generator. To identify the source of its integrated Bochner identity, center the product-Gamma law by writing X i = Y i − α i X_i=Y_i-\alpha_i X i = Y i − α i . Its canonical moment Hessian, expressed in the unchanged Y Y Y coordinates, is
H Γ = diag ( Y 1 , … , Y m ) . H_\Gamma=\diag(Y_1,\ldots,Y_m). H Γ = diag ( Y 1 , … , Y m ) . Indeed, writing ρ Γ \rho_\Gamma ρ Γ for the product-Gamma density,
div Γ ( H Γ ∇ G ) = ∑ i ρ Γ − 1 ∂ i ( ρ Γ Y i G i ) = ∑ i ( Y i G i i + ( α i − Y i ) G i ) = L Γ G . \Div_\Gamma(H_\Gamma\nabla G)
=\sum_i\rho_\Gamma^{-1}\partial_i(\rho_\Gamma Y_iG_i)
=\sum_i\bigl(Y_iG_{ii}+(\alpha_i-Y_i)G_i\bigr)=\calL_\Gamma G. div Γ ( H Γ ∇ G ) = i ∑ ρ Γ − 1 ∂ i ( ρ Γ Y i G i ) = i ∑ ( Y i G ii + ( α i − Y i ) G i ) = L Γ G . Thus Proposition 16.1 , specialized to H Γ H_\Gamma H Γ , gives
N Γ ( G ) : = E ( L Γ G ) 2 = E [ ∑ i Y i G i 2 + ∑ i , j Y i Y j G i j 2 ] . N_\Gamma(G):=\E(\calL_\Gamma G)^2
=\E\left[\sum_iY_iG_i^2+\sum_{i,j}Y_iY_jG_{ij}^2\right]. N Γ ( G ) := E ( L Γ G ) 2 = E [ i ∑ Y i G i 2 + i , j ∑ Y i Y j G ij 2 ] . Here the second sum is over ordered pairs: because D 2 G D^2G D 2 G is symmetric, each off-diagonal square occurs twice, exactly as in Tr ( H Γ D 2 G H Γ D 2 G ) \Tr(H_\Gamma D^2G\,H_\Gamma D^2G) Tr ( H Γ D 2 G H Γ D 2 G ) . There is no additional factor of 2 or 1 / 2 1/2 1/2 . Finally set D i ( G ) = E [ Y i 2 G i 2 ] / α i D_i(G)=\E[Y_i^2G_i^2]/\alpha_i D i ( G ) = E [ Y i 2 G i 2 ] / α i and D Γ = ∑ i D i D_\Gamma=\sum_iD_i D Γ = ∑ i D i .
For α i ≥ 1 \alpha_i\ge1 α i ≥ 1 , N Γ ( G ) − 1 4 D Γ ( G ) = ∑ i E R i + ∑ i δ i D i ( G ) N_\Gamma(G)-\tfrac14D_\Gamma(G)=\sum_i\E R_i+\sum_i\delta_iD_i(G) N Γ ( G ) − 4 1 D Γ ( G ) = ∑ i E R i + ∑ i δ i D i ( G ) with R i , δ i R_i,\delta_i R i , δ i as in (17.16) , both nonnegative.
The occurrence of the first-order row E [ Y i G i 2 ] \E[Y_iG_i^2] E [ Y i G i 2 ] and the full ordered Hessian row E [ ∑ j Y i Y j G i j 2 ] \E[\sum_jY_iY_jG_{ij}^2] E [ ∑ j Y i Y j G ij 2 ] below is precisely the i i i th-row contribution to the specialized Bochner identity (D29.14) . For Y ∼ Gamma ( a , 1 ) Y\sim\GammaLaw(a,1) Y ∼ Gamma ( a , 1 ) and smooth u u u with the usual decay,
E [ Y 2 u u ′ ] = 1 2 E [ Y 2 ( u 2 ) ′ ] = 1 2 Γ ( a ) ∫ y a + 1 e − y ( u 2 ) ′ d y = − 1 2 Γ ( a ) ∫ u 2 ( ( a + 1 ) y a − y a + 1 ) e − y d y , \E[Y^2uu']=\tfrac12\E[Y^2(u^2)']
=\tfrac1{2\Gamma(a)}\int y^{a+1}e^{-y}(u^2)'\dd y
=-\tfrac1{2\Gamma(a)}\int u^2\bigl((a+1)y^a-y^{a+1}\bigr)e^{-y}\dd y, E [ Y 2 u u ′ ] = 2 1 E [ Y 2 ( u 2 ) ′ ] = 2Γ ( a ) 1 ∫ y a + 1 e − y ( u 2 ) ′ d y = − 2Γ ( a ) 1 ∫ u 2 ( ( a + 1 ) y a − y a + 1 ) e − y d y , i.e. E [ Y 2 u u ′ ] = − 1 2 E [ ( ( a + 1 ) Y − Y 2 ) u 2 ] \E[Y^2uu']=-\tfrac12\E[((a+1)Y-Y^2)u^2] E [ Y 2 u u ′ ] = − 2 1 E [(( a + 1 ) Y − Y 2 ) u 2 ] . Expanding the first square of R i R_i R i and applying this with u = G i u=G_i u = G i , u ′ = G i i u'=G_{ii} u ′ = G ii , a = α i a=\alpha_i a = α i , the cross term is
− 2 α i + 1 E [ Y i 2 G i G i i ] = 1 α i + 1 E [ ( ( α i + 1 ) Y i − Y i 2 ) G i 2 ] = E [ Y i G i 2 ] − 1 α i + 1 E [ Y i 2 G i 2 ] . -\tfrac2{\alpha_i+1}\E[Y_i^2G_iG_{ii}]
=\tfrac1{\alpha_i+1}\E[((\alpha_i+1)Y_i-Y_i^2)G_i^2]
=\E[Y_iG_i^2]-\tfrac1{\alpha_i+1}\E[Y_i^2G_i^2]. − α i + 1 2 E [ Y i 2 G i G ii ] = α i + 1 1 E [(( α i + 1 ) Y i − Y i 2 ) G i 2 ] = E [ Y i G i 2 ] − α i + 1 1 E [ Y i 2 G i 2 ] . Therefore
E R i = E [ Y i G i 2 + ∑ j Y i Y j G i j 2 ] + ( 1 ( α i + 1 ) 2 − 1 α i + 1 ) E [ Y i 2 G i 2 ] , \E R_i=\E\Bigl[Y_iG_i^2+\sum_jY_iY_jG_{ij}^2\Bigr]
+\Bigl(\tfrac1{(\alpha_i+1)^2}-\tfrac1{\alpha_i+1}\Bigr)\E[Y_i^2G_i^2], E R i = E [ Y i G i 2 + j ∑ Y i Y j G ij 2 ] + ( ( α i + 1 ) 2 1 − α i + 1 1 ) E [ Y i 2 G i 2 ] , and 1 ( α i + 1 ) 2 − 1 α i + 1 = − α i ( α i + 1 ) 2 \tfrac1{(\alpha_i+1)^2}-\tfrac1{\alpha_i+1}=-\tfrac{\alpha_i}{(\alpha_i+1)^2} ( α i + 1 ) 2 1 − α i + 1 1 = − ( α i + 1 ) 2 α i . Summing over i i i , the bracket sums to N Γ ( G ) N_\Gamma(G) N Γ ( G ) and the correction to − ∑ i α i 2 ( α i + 1 ) 2 D i ( G ) -\sum_i\tfrac{\alpha_i^2}{(\alpha_i+1)^2}D_i(G) − ∑ i ( α i + 1 ) 2 α i 2 D i ( G ) . Subtracting 1 4 D Γ \tfrac14D_\Gamma 4 1 D Γ gives the claim, and δ i = α i 2 ( α i + 1 ) 2 − 1 4 ≥ 0 \delta_i=\tfrac{\alpha_i^2}{(\alpha_i+1)^2}-\tfrac14\ge0 δ i = ( α i + 1 ) 2 α i 2 − 4 1 ≥ 0 iff α i α i + 1 ≥ 1 2 \tfrac{\alpha_i}{\alpha_i+1}\ge\tfrac12 α i + 1 α i ≥ 2 1 iff α i ≥ 1 \alpha_i\ge1 α i ≥ 1 .
This lemma alone proves C M H ( 4 ) \mathrm{CMH}(4) CMH ( 4 ) for the product-Gamma law; for that consequence, α i ≥ 1 \alpha_i\ge1 α i ≥ 1 is used only to make δ i ≥ 0 \delta_i\ge0 δ i ≥ 0 . In the Dirichlet proof below the exact δ i \delta_i δ i term is instead retained inside F A ( α i , P i ) F_A(\alpha_i,P_i) F A ( α i , P i ) , and log-concavity is consumed when Lemma D29.4 is applied with a = α i ≥ 1 a=\alpha_i\ge1 a = α i ≥ 1 .
Homogeneity of degree 0 gives Euler’s identity ∑ j Y j G j = 0 \sum_jY_jG_j=0 ∑ j Y j G j = 0 ; differentiating in Y i Y_i Y i yields G i + ∑ j Y j G j i = 0 G_i+\sum_jY_jG_{ji}=0 G i + ∑ j Y j G ji = 0 , i.e. the single linear constraint ∑ j Y j G i j = − G i \sum_jY_jG_{ij}=-G_i ∑ j Y j G ij = − G i . Set t i = Y i ( G i i − G i / ( α i + 1 ) ) t_i=Y_i(G_{ii}-G_i/(\alpha_i+1)) t i = Y i ( G ii − G i / ( α i + 1 )) and t j = Y j G i j t_j=Y_jG_{ij} t j = Y j G ij for j ≠ i j\ne i j = i . Then
∑ j t j = ∑ j Y j G i j − Y i α i + 1 G i = − G i ( 1 + Y i α i + 1 ) , R i = t i 2 + ∑ j ≠ i Y i Y j t j 2 . \sum_jt_j=\sum_jY_jG_{ij}-\tfrac{Y_i}{\alpha_i+1}G_i=-G_i\Bigl(1+\tfrac{Y_i}{\alpha_i+1}\Bigr),
\qquad
R_i=t_i^2+\sum_{j\ne i}\tfrac{Y_i}{Y_j}t_j^2 . j ∑ t j = j ∑ Y j G ij − α i + 1 Y i G i = − G i ( 1 + α i + 1 Y i ) , R i = t i 2 + j = i ∑ Y j Y i t j 2 . The weights are w i = 1 w_i=1 w i = 1 and w j = Y i / Y j w_j=Y_i/Y_j w j = Y i / Y j , so ∑ j w j − 1 = 1 + ∑ j ≠ i Y j / Y i = S / Y i \sum_jw_j^{-1}=1+\sum_{j\ne i}Y_j/Y_i=S/Y_i ∑ j w j − 1 = 1 + ∑ j = i Y j / Y i = S / Y i . Minimizing ∑ j w j t j 2 \sum_jw_jt_j^2 ∑ j w j t j 2 subject to ∑ j t j = c \sum_jt_j=c ∑ j t j = c gives c 2 / ∑ j w j − 1 c^2/\sum_jw_j^{-1} c 2 / ∑ j w j − 1 , which is the stated bound.
Homogeneity further gives L Γ G = S − 1 L α g \calL_\Gamma G=S^{-1}L_\alpha g L Γ G = S − 1 L α g and Y i G i = u i ( P ) Y_iG_i=u_i(P) Y i G i = u i ( P ) : indeed G i = ∂ Y i g ( Y / S ) = S − 1 ( g i − ∑ k P k g k ) G_i=\partial_{Y_i}g(Y/S)=S^{-1}(g_i-\sum_kP_kg_k) G i = ∂ Y i g ( Y / S ) = S − 1 ( g i − ∑ k P k g k ) , so Y i G i = P i g i − P i ∑ k P k g k = u i ( P ) Y_iG_i=P_ig_i-P_i\sum_kP_kg_k=u_i(P) Y i G i = P i g i − P i ∑ k P k g k = u i ( P ) . For A > 2 A>2 A > 2 , E S − 1 = ( A − 1 ) − 1 \E S^{-1}=(A-1)^{-1} E S − 1 = ( A − 1 ) − 1 and E S − 2 = z A − 1 \E S^{-2}=z_A^{-1} E S − 2 = z A − 1 with z A = ( A − 1 ) ( A − 2 ) z_A=(A-1)(A-2) z A = ( A − 1 ) ( A − 2 ) . Since S ⊥ P S\perp P S ⊥ P ,
N Γ ( G ) = E [ S − 2 ] E ( L α g ) 2 = n α ( g ) z A , D Γ ( G ) = E ∑ i ( Y i G i ) 2 α i = d α ( g ) , N_\Gamma(G)=\E[S^{-2}]\,\E(L_\alpha g)^2=\frac{n_\alpha(g)}{z_A},
\qquad
D_\Gamma(G)=\E\sum_i\frac{(Y_iG_i)^2}{\alpha_i}=d_\alpha(g), N Γ ( G ) = E [ S − 2 ] E ( L α g ) 2 = z A n α ( g ) , D Γ ( G ) = E i ∑ α i ( Y i G i ) 2 = d α ( g ) , the second because Y i G i Y_iG_i Y i G i depends on P P P alone. Substituting Y i = S P i Y_i=SP_i Y i = S P i and G i = u i / ( S P i ) G_i=u_i/(SP_i) G i = u i / ( S P i ) into Lemma D29.3 and taking expectations gives (17.20) .
5. The scalar minimization ¶ For A ≥ 3 A\ge3 A ≥ 3 , a ≥ 1 a\ge1 a ≥ 1 , 0 < p ≤ 1 0<p\le1 0 < p ≤ 1 and z = z A z=z_A z = z A , the function F A F_A F A of (17.21) satisfies F A ( a , p ) ≥ A − 1 2 + s A F_A(a,p)\ge A-\tfrac12+s_A F A ( a , p ) ≥ A − 2 1 + s A with s A s_A s A as in (17.22) , and s A ≥ 0 s_A\ge0 s A ≥ 0 .
The p p p -dependent part of F A F_A F A is ϕ a ( p ) = a ( p − 1 + z p ( a + 1 ) − 2 ) \phi_a(p)=a\bigl(p^{-1}+zp(a+1)^{-2}\bigr) ϕ a ( p ) = a ( p − 1 + z p ( a + 1 ) − 2 ) , convex on ( 0 , ∞ ) (0,\infty) ( 0 , ∞ ) with unconstrained minimizer p ∗ = ( a + 1 ) / z p_*=(a+1)/\sqrt z p ∗ = ( a + 1 ) / z . Hence on ( 0 , 1 ] (0,1] ( 0 , 1 ] the minimum is at p = 1 p=1 p = 1 when p ∗ ≥ 1 p_*\ge1 p ∗ ≥ 1 , i.e. z ≤ a + 1 \sqrt z\le a+1 z ≤ a + 1 , and at p ∗ p_* p ∗ otherwise.
Case z ≤ 4 z\le4 z ≤ 4 . Then z ≤ 2 ≤ a + 1 \sqrt z\le2\le a+1 z ≤ 2 ≤ a + 1 for every a ≥ 1 a\ge1 a ≥ 1 , so p = 1 p=1 p = 1 and
F A ( a , 1 ) = a + 2 a ( A − 2 ) a + 1 + z a a + 1 − z 4 . F_A(a,1)=a+\frac{2a(A-2)}{a+1}+\frac{za}{a+1}-\frac z4 . F A ( a , 1 ) = a + a + 1 2 a ( A − 2 ) + a + 1 z a − 4 z . Each of the three terms is increasing in a a a on [ 1 , ∞ ) [1,\infty) [ 1 , ∞ ) , so the minimum is at a = 1 a=1 a = 1 : 1 + ( A − 2 ) + z / 2 − z / 4 = A − 1 + z / 4 = A − 1 2 + z − 2 4 1+(A-2)+z/2-z/4=A-1+z/4=A-\tfrac12+\tfrac{z-2}4 1 + ( A − 2 ) + z /2 − z /4 = A − 1 + z /4 = A − 2 1 + 4 z − 2 .
Case z ≥ 4 z\ge4 z ≥ 4 . At p = p ∗ p=p_* p = p ∗ the two terms of ϕ a \phi_a ϕ a are equal, so ϕ a ( p ∗ ) = 2 a z / ( a + 1 ) \phi_a(p_*)=2a\sqrt z/(a+1) ϕ a ( p ∗ ) = 2 a z / ( a + 1 ) . Writing r = a / ( a + 1 ) ∈ [ 1 2 , 1 ) r=a/(a+1)\in[\tfrac12,1) r = a / ( a + 1 ) ∈ [ 2 1 , 1 ) for a ≥ 1 a\ge1 a ≥ 1 ,
F A = 2 r z + 2 r ( A − 2 ) + z r 2 − z 4 = 2 r ( z + A − 2 ) + z r 2 − z 4 , F_A=2r\sqrt z+2r(A-2)+zr^2-\frac z4=2r(\sqrt z+A-2)+zr^2-\frac z4, F A = 2 r z + 2 r ( A − 2 ) + z r 2 − 4 z = 2 r ( z + A − 2 ) + z r 2 − 4 z , increasing in r r r on [ 0 , 1 ) [0,1) [ 0 , 1 ) , so the minimum is at r = 1 2 r=\tfrac12 r = 2 1 (a = 1 a=1 a = 1 ) and equals ( z + A − 2 ) + z / 4 − z / 4 = A − 2 + z = A − 1 2 + z − 3 2 (\sqrt z+A-2)+z/4-z/4=A-2+\sqrt z=A-\tfrac12+\sqrt z-\tfrac32 ( z + A − 2 ) + z /4 − z /4 = A − 2 + z = A − 2 1 + z − 2 3 . In the complementary region z ≤ a + 1 \sqrt z\le a+1 z ≤ a + 1 the value F A ( a , 1 ) F_A(a,1) F A ( a , 1 ) is increasing in a a a and at the interface is not smaller, so the stated minimum stands.
Finally A ≥ 3 A\ge3 A ≥ 3 gives z A = ( A − 1 ) ( A − 2 ) ≥ 2 z_A=(A-1)(A-2)\ge2 z A = ( A − 1 ) ( A − 2 ) ≥ 2 , so s A = ( z − 2 ) / 4 ≥ 0 s_A=(z-2)/4\ge0 s A = ( z − 2 ) /4 ≥ 0 when z ≤ 4 z\le4 z ≤ 4 and s A = z − 3 2 ≥ 2 − 3 2 > 0 s_A=\sqrt z-\tfrac32\ge2-\tfrac32>0 s A = z − 2 3 ≥ 2 − 2 3 > 0 when z ≥ 4 z\ge4 z ≥ 4 .
6. Proof of the Dirichlet theorem ¶ If m = 2 m=2 m = 2 and A < 3 A<3 A < 3 , the law is one dimensional and Theorem D29.1 applies (a B e t a \mathrm{Beta} Beta law is log-concave precisely when both parameters are ≥ 1 \ge1 ≥ 1 ). In every remaining case A = ∑ i α i ≥ 3 A=\sum_i\alpha_i\ge3 A = ∑ i α i ≥ 3 — automatically so when m ≥ 3 m\ge3 m ≥ 3 — hence z A ≥ 2 > 0 z_A\ge2>0 z A ≥ 2 > 0 and the inverse moments E S − 1 , E S − 2 \E S^{-1},\E S^{-2} E S − 1 , E S − 2 are finite.
By Lemma D29.2 and (17.20) , since D i ( G ) = E [ u i 2 ] / α i D_i(G)=\E[u_i^2]/\alpha_i D i ( G ) = E [ u i 2 ] / α i , the coefficient of E P [ u i 2 ] \E_P[u_i^2] E P [ u i 2 ] in N Γ ( G ) − 1 4 D Γ ( G ) N_\Gamma(G)-\tfrac14D_\Gamma(G) N Γ ( G ) − 4 1 D Γ ( G ) is at least
1 z A P i + 2 ( A − 1 ) ( α i + 1 ) + P i ( α i + 1 ) 2 + δ i α i . \frac1{z_AP_i}+\frac2{(A-1)(\alpha_i+1)}+\frac{P_i}{(\alpha_i+1)^2}+\frac{\delta_i}{\alpha_i}. z A P i 1 + ( A − 1 ) ( α i + 1 ) 2 + ( α i + 1 ) 2 P i + α i δ i . Multiply by z A α i z_A\alpha_i z A α i and use z A / ( A − 1 ) = A − 2 z_A/(A-1)=A-2 z A / ( A − 1 ) = A − 2 together with z A δ i = z A α i 2 / ( α i + 1 ) 2 − z A / 4 z_A\delta_i=z_A\alpha_i^2/(\alpha_i+1)^2-z_A/4 z A δ i = z A α i 2 / ( α i + 1 ) 2 − z A /4 :
α i P i + 2 α i ( A − 2 ) α i + 1 + z A α i ( α i + P i ) ( α i + 1 ) 2 − z A 4 = F A ( α i , P i ) . \frac{\alpha_i}{P_i}+\frac{2\alpha_i(A-2)}{\alpha_i+1}
+\frac{z_A\alpha_i(\alpha_i+P_i)}{(\alpha_i+1)^2}-\frac{z_A}4
=F_A(\alpha_i,P_i). P i α i + α i + 1 2 α i ( A − 2 ) + ( α i + 1 ) 2 z A α i ( α i + P i ) − 4 z A = F A ( α i , P i ) . Lemma D29.4 bounds this below by A − 1 2 + s A A-\tfrac12+s_A A − 2 1 + s A pointwise in P i P_i P i , so
N Γ ( G ) − 1 4 D Γ ( G ) ≥ A − 1 2 + s A z A ∑ i E P [ u i 2 ] α i = A − 1 2 + s A z A d α ( g ) . N_\Gamma(G)-\tfrac14D_\Gamma(G)
\ \ge\ \frac{A-\tfrac12+s_A}{z_A}\sum_i\frac{\E_P[u_i^2]}{\alpha_i}
=\frac{A-\tfrac12+s_A}{z_A}\,d_\alpha(g). N Γ ( G ) − 4 1 D Γ ( G ) ≥ z A A − 2 1 + s A i ∑ α i E P [ u i 2 ] = z A A − 2 1 + s A d α ( g ) . Substituting N Γ = n α / z A N_\Gamma=n_\alpha/z_A N Γ = n α / z A and D Γ = d α D_\Gamma=d_\alpha D Γ = d α , multiplying by z A z_A z A , and using
z A 4 + A − 1 2 = A 2 − 3 A + 2 4 + 4 A − 2 4 = A ( A + 1 ) 4 , \frac{z_A}4+A-\frac12=\frac{A^2-3A+2}4+\frac{4A-2}4=\frac{A(A+1)}4, 4 z A + A − 2 1 = 4 A 2 − 3 A + 2 + 4 4 A − 2 = 4 A ( A + 1 ) , we obtain n α ( g ) ≥ ( A ( A + 1 ) 4 + s A ) d α ( g ) n_\alpha(g)\ge\bigl(\tfrac{A(A+1)}4+s_A\bigr)d_\alpha(g) n α ( g ) ≥ ( 4 A ( A + 1 ) + s A ) d α ( g ) . Since s A ≥ 0 s_A\ge0 s A ≥ 0 this gives A ( A + 1 ) d α ( g ) ≤ 4 n α ( g ) A(A+1)d_\alpha(g)\le4n_\alpha(g) A ( A + 1 ) d α ( g ) ≤ 4 n α ( g ) .
For A > 3 A>3 A > 3 the Gamma branch of the preceding proof applies, including when m = 2 m=2 m = 2 , and gives n α ≥ ( A ( A + 1 ) / 4 + s A ) d α n_\alpha\ge\bigl(A(A+1)/4+s_A\bigr)d_\alpha n α ≥ ( A ( A + 1 ) /4 + s A ) d α . Since z A > 2 z_A>2 z A > 2 , the definition of s A s_A s A gives s A > 0 s_A>0 s A > 0 , and rearrangement proves the claim.
Every log-concave Dirichlet law has C P a f f ≤ 4 \CPaff\le4 C P aff ≤ 4 , as does every product, linear image, and convolution of independent such laws.
Qualitative dimension-free KLS bounds for simplices and conservative Gamma models are prior art; see Kolesnikov & Milman, 2016, §1.2 as a secondary pointer. The specific content of this proof is the constant 4, the quantitative surplus, and the CMH-level estimate; no broader priority claim is made.
Closure. The argument is written for smooth g g g with controlled boundary behavior. To close: take polynomials on the simplex, lift them after a radial cutoff { S ≥ ε } \{S\ge\eps\} { S ≥ ε } so that all Gamma integrations by parts in Lemma D29.2 are justified, and let ε ↓ 0 \eps\downarrow0 ε ↓ 0 ; the inverse moments E S − 1 , E S − 2 \E S^{-1},\E S^{-2} E S − 1 , E S − 2 used in (17.19) are finite because A ≥ 3 A\ge3 A ≥ 3 throughout the Gamma branch. Polynomials form a core for L α L_\alpha L α , so the closed forms extend the inequality to Dom ( L α ) \Dom(L_\alpha) Dom ( L α ) . Only the residual branch m = 2 m=2 m = 2 , A < 3 A<3 A < 3 uses the one-dimensional no-flux closure of Theorem D29.1 instead.
Obstructions respected. No bounded_by edge applies: the obstruction statements of the manuscript are scoped by their own statements to the fixed-cut Eldan program. The only one with method-level reach, rem:projection-ceiling, forbids deriving quadratic-chaos thin shell from radial or projection information alone; the Dirichlet proof uses the full Hessian row through the Euler constraint of Lemma D29.3 , not projection tests, and makes no thin-shell claim. Consistency with the sharp external inputs holds: by Theorem 17.1 , the boundary mechanism forcing the constant 4 is the one-dimensional exponential, which is also the extremal case of Theorem 4.2 .
Theorem D29.1 and Theorem D29.2 prove that products of centered one-sided exponentials have C C M H = 4 \CMH=4 C CMH = 4 exactly. Proposition 16.2 proves, for each admissible test function, the exact decomposition of the CMH numerator into the affine Poincaré contribution and a nonnegative solenoidal contribution.
It is an open question, not a consequence of this corollary, whether an admissible perturbation raises the full CMH Rayleigh quotient. The covariance, canonical Stein kernel, both numerator channels, denominator, and optimizer all vary, so increasing a solenoidal term in isolation is insufficient. The required full second-variation calculation is Conjecture 17.1 .
What this does not show. Theorem D29.3 is a family result, not evidence for universal C M H ( 4 ) \mathrm{CMH}(4) CMH ( 4 ) . By Corollary 17.2 the Dirichlet family is strictly inside the bound, and by Theorem D29.2 the saturating cases are products of centered one-sided exponentials, where the slack is exactly zero. Whether that endpoint is perturbatively unstable is precisely the open issue just stated.
Kannan, R., Lovász, L., & Simonovits, M. (1995). Isoperimetric Problems for Convex Bodies and a Localization Lemma. Discrete & Computational Geometry , 13 (3–4), 541–559. 10.1007/BF02574061 Cattiaux, P., & Guillin, A. (2018). On the Poincaré Constant of Log-Concave Measures . Kolesnikov, A. V., & Milman, E. (2016). The KLS Isoperimetric Conjecture for Generalized Orlicz Balls .