Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

Letwin: matrix and quadratic estimates

Part of the results of the literature written out here, Chapter Family 4: moment maps, Monge–Ampère, and Stein kernels; the reading order is on the full proofs page.

Overview. This dossier gives a proof of the precise matrix and quadratic estimates imported as Theorem 4.2 and Theorem 25.1. It follows the mechanism of Letwin, 2026, pinned to arXiv:2607.24164v1, and expands the coordinate change and the analytic domains needed in that mechanism. This is an author’s dossier, not its independent certification.

  1. On a bounded regular target, construct cutoffs in the source coordinates and justify integrating the differentiated Monge–Ampère identity before using its positive terms.

  2. Transform the third-order tensor and the fixed matrix together. The pointwise comparison then combines with Brascamp–Lieb to give the matrix estimate.

  3. Transfer the Stein kernel by congruence, apply the published H−1H^{-1} inequality on a generally non-isotropic law, and remove singular matrices and target regularity.

The source locations are Theorem 2.5 (PDF pp. 8–13), Theorem 1.2 (p. 3, proof on pp. 14–16), and Appendix A (pp. 16–18). Page numbers refer to the PDF, with its first page numbered one. The versioned HTML also carries the theorem and lemma numbers. The proof of the general KLS bound, Theorem 1.1 of the source, is outside this dossier.

Statements and published inputs

The established inputs used below are the compact-target regularity theorem Theorem 4.1, the pointwise estimate of Klartag, 2014, Theorem 1.1, the convex-boundary Brascamp–Lieb inequality Kolesnikov & Milman, 2017, Theorem 1.2(1), and Barthe & Klartag, 2019, Propositions 10 and 27. We use these as published results, not the conclusions of either recent preprint. The application of each is specified below. Gaussian convolution preserves log-concavity, and finite-dimensional log-concave probabilities have finite polynomial moments; these standard facts are used only in the last approximation step.

Regularity, symmetry and integration

The regularity node applies with P=K‾P=\overline K. Indeed, multiplying VV by a smooth cutoff equal to one near PP gives a globally smooth extension of V∣PV|_P; its exponential is positive and globally smooth, as required by that node. The extension need not be convex outside PP. Thus φ\varphi is smooth, H≻0H\succ0, and ∇φ\nabla\varphi maps Rn\mathbb R^n diffeomorphically onto KK. The change-of-variables equation is

log⁡det⁡H=−φ+V(∇φ).\log\det H=-\varphi+V(\nabla\varphi).

The hypotheses of Klartag’s Theorem 1.1 are satisfied: VV and each of its derivatives are bounded on KK, by smoothness near the compact set K‾\overline K. It gives Tr⁡H≤2RK2\operatorname{Tr}H\le2R_K^2, where RK=sup⁡x∈K∣x∣R_K=\sup_{x\in K}|x|. In particular HH and ∇φ\nabla\varphi are bounded; no positive uniform lower bound on HH is asserted or used.

Use HabH^{ab} for the entries of H−1H^{-1} and define

Lf=Hab∂abf−Va(∇φ)∂af,Γ(f,g)=Hab∂af∂bg.Lf=H^{ab}\partial_{ab}f-V_a(\nabla\varphi)\partial_af, \qquad \Gamma(f,g)=H^{ab}\partial_af\partial_bg.

Repeated indices are summed. Differentiating the inverse Hessian and using symmetry of third derivatives gives ∂aHab=−Hbi∂ilog⁡det⁡H\partial_a H^{ab}=-H^{bi}\partial_i\log\det H. The logarithmic equation therefore implies the divergence formula

Lf=eφ∂a(e−φHab∂bf).Lf=e^{\varphi}\partial_a(e^{-\varphi}H^{ab}\partial_bf).

Consequently ∫(Lf)g dν=−∫Γ(f,g) dν\int (Lf)g\,d\nu=-\int\Gamma(f,g)\,d\nu when either smooth function has compact support. This assertion is local and requires no tail estimate on H−1H^{-1}.

Ordinary Euclidean integration by parts with a cutoff ζ(y/R)\zeta(y/R) gives E∂if=Ef∂iφ\mathbb E\partial_i f=\mathbb E f\partial_i\varphi whenever all displayed terms and ff are integrable: the extra term is bounded by CR−1E∣f∣C R^{-1}\mathbb E|f|. Applying it to f=∂jφf=\partial_j\varphi, whose derivative and all other terms are bounded, gives

EνHij=Eν(∂iφ∂jφ)=EμXiXj=δij.\mathbb E_\nu H_{ij}=\mathbb E_\nu(\partial_i\varphi\partial_j\varphi) =\mathbb E_\mu X_iX_j=\delta_{ij}.

For later use, every bounded smooth uu satisfies Var⁡νu≤∫Γ(u) dν\operatorname{Var}_\nu u\le\int\Gamma(u)\,d\nu, allowing infinity on the right. Here the precise boundary input is Kolesnikov–Milman, Theorem 1.2(1) (DOI). Set MR={φ≤R}M_R=\{\varphi\le R\} for R>min⁡φR>\min\varphi, with its Euclidean metric and conditional probability νR=ν∣MR/ν(MR)\nu_R=\nu|_{M_R}/\nu(M_R). These closed sublevels are compact, connected and convex. Since H≻0H\succ0, their boundaries are smooth and ∇φ\nabla\varphi does not vanish there. With outward normal ∇φ/∣∇φ∣\nabla\varphi/|\nabla\varphi|, their second fundamental form on tangent vectors is H/∣∇φ∣⪰0H/|\nabla\varphi|\succeq0 (vacuously on the zero-dimensional boundary when n=1n=1). The conditional potential is φ+log⁡ν(MR)\varphi+\log\nu(M_R), so the Euclidean N=∞N=\infty curvature tensor in the cited theorem is exactly D2φ=H≻0D^2\varphi=H\succ0; its prefactor N/(N−1)N/(N-1) equals one. The theorem applies to arbitrary C1(MR)C^1(M_R) tests, with no Neumann condition on uu, and yields

Var⁡νRu≤1ν(MR)∫MR∇uTH−1∇u dν.\operatorname{Var}_{\nu_R}u \le\frac{1}{\nu(M_R)}\int_{M_R}\nabla u^TH^{-1}\nabla u\,d\nu.

The domains exhaust Rn\mathbb R^n and ν(MR)→1\nu(M_R)\to1. Boundedness of uu passes its conditional first and second moments; monotone convergence passes the unnormalized nonnegative energy integrals. This proves the claimed inequality, including the case of infinite energy.

The full tensor comparison

Differentiating the logarithmic Monge–Ampère identity twice gives

LH+H=A+Q,A=H(D2V∘∇φ)H,Qij=HacHbdφabiφcdj.LH+H=A+Q,\qquad A=H(D^2V\circ\nabla\varphi)H,\qquad Q_{ij}=H^{ac}H^{bd}\varphi_{abi}\varphi_{cdj}.

For clarity, its first derivative is Habφabi=−φi+VaφaiH^{ab}\varphi_{abi}=-\varphi_i+V_a\varphi_{ai}, and differentiating it uses ∂jHab=−HacφcdjHdb\partial_jH^{ab}=-H^{ac}\varphi_{cdj}H^{db}. The drift term VaφaijV_a\varphi_{aij} moves to the left and gives LHijLH_{ij}. Both AA and QQ are positive semidefinite: convexity gives the first, and the second is the Gram matrix of H−1/2(∂iH)H−1/2H^{-1/2}(\partial_iH)H^{-1/2}.

Fix B⪰0B\succeq0 and introduce

S=Tr⁡(BHBH),DB=HabTr⁡(B∂aHB∂bH),CB=Tr⁡(BHBQ),AB=Tr⁡(BHBA).S=\operatorname{Tr}(BHBH),\quad D_B=H^{ab}\operatorname{Tr}(B\partial_aH B\partial_bH),\quad C_B=\operatorname{Tr}(BHBQ),\quad A_B=\operatorname{Tr}(BHBA).

All four are nonnegative. For DBD_B this follows by contracting the Gram matrix of B1/2(∂aH)B1/2B^{1/2}(\partial_aH)B^{1/2} with H−1H^{-1}; for CB,ABC_B,A_B it follows from the trace of a product of positive semidefinite matrices. Product differentiation gives

LS=−2S+2(DB+AB+CB).LS=-2S+2(D_B+A_B+C_B).

SS is bounded. The cutoff lemma with c=2c=2 first gives LS∈L1LS\in L^1 and ELS=0\mathbb E LS=0. The nonnegative sum is therefore integrable, and each summand is integrable separately. Only at this point do we obtain

ES=EDB+EAB+ECB.\mathbb E S=\mathbb E D_B+\mathbb E A_B+\mathbb E C_B.

The pointwise comparison is CB≥DBC_B\ge D_B. Here is the coordinate calculation including the invariance that permits normalization. At the point in question take an invertible constant matrix PP, put y=Pzy=Pz, and transform

H~=PTHP,B~=P−1BP−T,T~abc=PiaPjbPkcTijk,Tijk=φijk.\widetilde H=P^THP,\quad \widetilde B=P^{-1}BP^{-T},\quad \widetilde T_{abc}=P_{ia}P_{jb}P_{kc}T_{ijk},\qquad T_{ijk}=\varphi_{ijk}.

The inverse Hessian transforms as P−1H−1P−TP^{-1}H^{-1}P^{-T}. Contracting two copies of T~\widetilde T with two inverse Hessians gives Q~=PTQP\widetilde Q=P^TQP; hence Tr⁡(B~H~B~Q~)=CB\operatorname{Tr}(\widetilde B\widetilde H\widetilde B\widetilde Q)=C_B. Moreover

∂zaH~=∑iPiaPT(∂yiH)P,\partial_{z_a}\widetilde H =\sum_i P_{ia}P^T(\partial_{y_i}H)P,

so cyclicity of the trace and ∑ab(H~−1)abPiaPjb=(H−1)ij\sum_{ab}(\widetilde H^{-1})_{ab}P_{ia}P_{jb}=(H^{-1})_{ij} show D~B~=DB\widetilde D_{\widetilde B}=D_B. Thus each scalar in the comparison is invariant. Choose PP so that PTHP=IP^THP=I, followed by an orthogonal change that diagonalizes B~=diag⁡(b1,…,bn)\widetilde B=\operatorname{diag}(b_1,\ldots,b_n). It remains symmetric and positive semidefinite by congruence. In these coordinates

DB=∑ijkbibjTijk2,CB=∑ijkbi2Tijk2,D_B=\sum_{ijk}b_i b_j T_{ijk}^2,\qquad C_B=\sum_{ijk}b_i^2T_{ijk}^2,

and full symmetry of TT implies

CB−DB=12∑ijk(bi−bj)2Tijk2≥0.C_B-D_B=\frac12\sum_{ijk}(b_i-b_j)^2T_{ijk}^2\ge0.

This is a change at a single point with a constant PP; no derivative of a moving normalizing frame is taken. The global mean identity and Brascamp–Lieb step remain in the original coordinates. We conclude that EDB≤12ES\mathbb E D_B\le\tfrac12\mathbb E S.

Proof of the matrix assertion

Quadratic functions, congruence, and singular matrices

For a centered log-concave probability λ\lambda, use the homogeneous dual norm

∥h∥H−1(λ)=sup⁡{∫hg dλ: g locally Lipschitz,g∈L2(λ),∫∣∇g∣2 dλ≤1},∫h dλ=0.\|h\|_{H^{-1}(\lambda)} =\sup\left\{\int hg\,d\lambda:\ g\text{ locally Lipschitz},\quad g\in L^2(\lambda),\quad\int|\nabla g|^2\,d\lambda\le1\right\}, \qquad \int h\,d\lambda=0.

The published Barthe–Klartag Proposition 10 states that Var⁡λf≤∑i∥∂if∥H−1(λ)2\operatorname{Var}_\lambda f\le\sum_i\|\partial_i f\|_{H^{-1}(\lambda)}^2 provided f,∂if∈L2f,\partial_i f\in L^2 and ∫∂if dλ=0\int\partial_i f\,d\lambda=0 for every ii. It requires log-concavity, not isotropy. Their Proposition 27 supplies density of Cc∞(Rn)C_c^\infty(\mathbb R^n) in the weighted H1H^1 norm.

In target coordinates let τ(x)=H((∇φ)−1(x))\tau(x)=H((\nabla\varphi)^{-1}(x)). For g∈Cc∞(Rn)g\in C_c^\infty(\mathbb R^n) the bounded functions g(∇φ)g(\nabla\varphi) and their derivatives permit ordinary source integration by parts:

EμXig(X)=Eνφig(∇φ)=Eν∑jHij(∂jg)(∇φ)=Eμ∑jτij(X)∂jg(X).\mathbb E_\mu X_i g(X) =\mathbb E_\nu\varphi_i g(\nabla\varphi) =\mathbb E_\nu\sum_jH_{ij}(\partial_jg)(\nabla\varphi) =\mathbb E_\mu\sum_j\tau_{ij}(X)\partial_jg(X).

If a centered full-dimensional log-concave λ\lambda has a symmetric Stein kernel τλ∈L2(λ)\tau_\lambda\in L^2(\lambda), Cauchy–Schwarz gives for such gg

∣E(v⋅Z)g(Z)∣≤(E∣τλ(Z)v∣2)1/2(E∣∇g(Z)∣2)1/2.\left|\mathbb E (v\cdot Z)g(Z)\right| \le\left(\mathbb E|\tau_\lambda(Z)v|^2\right)^{1/2} \left(\mathbb E|\nabla g(Z)|^2\right)^{1/2}.

Proposition 27 extends this to every test in the dual norm: both sides of the Stein identity converge in an H1H^1 approximation by Cauchy–Schwarz, since v⋅Zv\cdot Z and τλv\tau_\lambda v are square integrable. Therefore ∥v⋅z∥H−1(λ)2≤E∣τλ(Z)v∣2\|v\cdot z\|_{H^{-1}(\lambda)}^2\le\mathbb E|\tau_\lambda(Z)v|^2.

Removing target regularity and matching the normalization

Transfer, hypotheses, and exclusions

The matrix statement matches source Theorem 2.5 on exactly the class of its Lemma A.1; the quadratic statement matches source Theorem 1.2 on the whole isotropic log-concave class. The preceding argument uses no regular moment potential for the limiting law. It does not assert convergence of Hessians under the approximation, only convergence of the polynomial moments needed for the quadratic variance.

HypothesisWhere it is used
Centering and covariance IIEH=I\mathbb EH=I, centered derivatives, and the final gradient normalization
Convex target and smooth convex VV near its closureRegular moment map, bounded coefficients, and positivity of AA
Bounded targetBounded ∇φ,H,LW\nabla\varphi,H,LW and bounded regular quadratic tests
Constant symmetric BBProduct differentiation and the tensor comparison; no derivatives of BB occur
Invertible MM at the intermediate stepFull-dimensional congruence and orthogonality of JJ; removed by an explicit limit
Log-concavity of the transported lawPublished H−1H^{-1} inequality and Sobolev density; isotropy is not needed there
Finite polynomial moments of the original lawGaussian conditioning and the fourth-moment passage

Fences respected. Neither imported node has a bounded_by edge. The program’s relevant boundaries are nevertheless explicit. Proposition 16.3 is respected because the conclusion is a fixed-matrix trace inequality, not EH2⪯cI\mathbb E H^2\preceq cI. Remark 25.4 is respected because the proof uses a full tensor comparison, not projection-only data. The two-tail scaling boundary Remark 25.3 is untouched: no unwhitened cut source is bounded here. The occupation, bootstrap and moving-competitor boundaries Remark 33.6, Remark 33.5, Remark 32.6 and Remark 31.6 are untouched because no stochastic-time or cut-dependent assertion is proved. No comparison between the trace-upgrade cluster members is asserted, and no open antecedent is discharged by assumption.

References
  1. Letwin, B. (2026). The KLS Constant is O(\log1/4 n).
  2. Klartag, B. (2014). Logarithmically-Concave Moment Measures I. In Geometric Aspects of Functional Analysis (Vol. 2116, pp. 231–260). Springer. 10.1007/978-3-319-09477-9_16
  3. Kolesnikov, A. V., & Milman, E. (2017). Brascamp–Lieb-Type Inequalities on Weighted Riemannian Manifolds with Boundary. The Journal of Geometric Analysis, 27(2), 1680–1702. 10.1007/s12220-016-9736-5
  4. Barthe, F., & Klartag, B. (2019). Spectral Gaps, Symmetries and Log-Concave Perturbations.