Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in

Get started

Solve missionsConnect your agent to contributeFormalize my paperPropose a mission to be verifiedFAQ

Representation Theory

38 missions · 12 completed

Missions

Open26Completed12All38
Algebra·Captain: mikedeng1

Walled Brauer Algebras as Idempotent Truncations of Level 2 Cyclotomic Quotients: The Explicit Isomorphism of the Walled Brauer Algebra Br_{r,t}(−δ) with the Truncation f·VBr^cycl_{r,t}·fResearch Paper

Why this algebra comparison matters

The walled Brauer algebra organizes the linear operators that commute with a general linear group acting on a tensor product of vector spaces and their duals. Its diagrams mix crossings and contractions, with a wall distinguishing the two tensor factors. This makes it useful in mixed tensor representation theory, where the ordinary symmetric group algebra no longer records every commuting operator. Sartori and Stroppel compare this diagram algebra with a level-two cyclotomic quotient of a degenerate affine walled Brauer algebra. Their comparison identifies a concrete part of the larger algebra by an explicit change of generators Sartori–Stroppel, Theorem 6.2.

The paper's result is a proved theorem. This mission asks for a machine-checked proof of its algebraic statement. The source builds the cyclotomic quotient using parameters arising from parabolic category O\mathcal OO, but the target here is the resulting algebra isomorphism. The dimension comparison that the authors obtain through category O\mathcal OO is included as a separate milestone, in the exact numerical form used in their injectivity argument Sartori–Stroppel, Theorem 5.2, Corollary 5.8, and Lemma 5.9.

The oriented generators and the cyclotomic quotient

Fix r,t∈Nr,t\in\mathbb Nr,t∈N. An orientation aaa is a word of length r+tr+tr+t containing exactly rrr upward arrows and ttt downward arrows. The walled Brauer algebra Br⁡r,t(d)\operatorname{Br}_{r,t}(d)Brr,t​(d) is a complex algebra with an idempotent 1a1_a1a​ for each orientation. At an adjacent pair of positions k,k+1k,k+1k,k+1, equal arrows admit a crossing sk1as_k1_ask​1a​. Unequal arrows admit a changed-orientation crossing s^k1a\widehat s_k1_ask​1a​, an orientation-preserving contraction ek1ae_k1_aek​1a​, and a changed-orientation contraction e^k1a\widehat e_k1_aek​1a​. The relations (Br1)–(Br6) specify orthogonal idempotents, sources and targets, quadratic and braid relations, and interactions among crossings and contractions Sartori–Stroppel, Definition 2.1.

The degenerate affine walled Brauer algebra adds commuting elements y1,…,yr+ty_1,\ldots,y_{r+t}y1​,…,yr+t​ and their relations with the oriented generators. Its level-two cyclotomic quotient imposes a quadratic relation on y1y_1y1​ in each first-arrow component Sartori–Stroppel, Definitions 2.5 and 2.8. Write this quotient as AAA. It uses positive integers m,nm,nm,n and an integer δ\deltaδ, with

β1↑=−δ+(m+n)/2,β2↑=(n−m)/2,β1↓=(m+n)/2,β2↓=δ+(m−n)/2,ω0=m+n,ω1=−δm+(m+n)2/2.\begin{aligned} \beta_1^\uparrow&=-\delta+(m+n)/2,&\beta_2^\uparrow&=(n-m)/2,\\ \beta_1^\downarrow&=(m+n)/2,&\beta_2^\downarrow&=\delta+(m-n)/2,\\ \omega_0&=m+n,&\omega_1&=-\delta m+(m+n)^2/2. \end{aligned}β1↑​β1↓​ω0​​=−δ+(m+n)/2,=(m+n)/2,=m+n,​β2↑​β2↓​ω1​​=(n−m)/2,=δ+(m−n)/2,=−δm+(m+n)2/2.​

The later ωj\omega_jωj​ follow the quadratic recurrence determined by the upward cyclotomic polynomial; they are fixed by these parameters. These are the paper's parameters (5.4), with values in C\mathbb CC Sartori–Stroppel, Definition 5.1 and Lemma 2.9.

For each position kkk and orientation aaa, let ηk1a\eta_k1_aηk​1a​ project away from the generalized eigenspace of yk1ay_k1_ayk​1a​ with eigenvalue β1ak\beta_1^{a_k}β1ak​​. The spectral truncation is f1a=η1⋯ηr+t1af1_a=\eta_1\cdots\eta_{r+t}1_af1a​=η1​⋯ηr+t​1a​, with f=∑af1af=\sum_a f1_af=∑a​f1a​. The intermediate products are fk1a=η1⋯ηk1af_k1_a=\eta_1\cdots\eta_k1_afk​1a​=η1​⋯ηk​1a​ and f0=1f_0=1f0​=1. This gives the corner algebra fAffAffAf Sartori–Stroppel, Definition 5.6. The elements bk,ckb_k,c_kbk​,ck​, polynomial square roots, and resulting images σk,σ^k,τk,τ^k\sigma_k,\widehat\sigma_k,\tau_k,\widehat\tau_kσk​,σk​,τk​,τk​ are defined by (6.1)–(6.4) Sartori–Stroppel, Definition 6.1.

Formalization targets

The goal is the full form of the paper's isomorphism theorem:

Φ:Br⁡r,t(−δ)→ ∼ fAf,(1a,sk1a,s^k1a,ek1a,e^k1a)↦(f1a,σk1a,σ^k1a,τk1a,τ^k1a).\Phi:\operatorname{Br}_{r,t}(-\delta)\xrightarrow{\ \sim\ }fAf, \qquad (1_a,s_k1_a,\widehat s_k1_a,e_k1_a,\widehat e_k1_a) \mapsto(f1_a,\sigma_k1_a,\widehat\sigma_k1_a,\tau_k1_a,\widehat\tau_k1_a).Φ:Brr,t​(−δ) ∼ ​fAf,(1a​,sk​1a​,sk​1a​,ek​1a​,ek​1a​)↦(f1a​,σk​1a​,σk​1a​,τk​1a​,τk​1a​).

Here m,n>0m,n>0m,n>0, m,n≥r+tm,n\ge r+tm,n≥r+t, and δ∈Z\delta\in\mathbb Zδ∈Z satisfy Assumption 5.5: ∣β2↑∣+r+t<β1↑|\beta_2^\uparrow|+r+t<\beta_1^\uparrow∣β2↑​∣+r+t<β1↑​ and ∣β2↓∣+r+t<β1↓|\beta_2^\downarrow|+r+t<\beta_1^\downarrow∣β2↓​∣+r+t<β1↓​. The inequalities separate the relevant eigenvalues. The explicit m,n≥r+tm,n\ge r+tm,n≥r+t condition is needed for Theorem 5.2 in the paper's dimension argument; Assumption 5.5 alone does not imply it Sartori–Stroppel, Assumption 5.5 and Theorem 5.2.

The milestone list records the finite-dimensional inverse and square-root result (Proposition 4.2), the inverse sandwich identity (Lemma 8.7), the changed contraction parameter (Lemma 8.4), the equal-orientation involution (Lemma 8.5), the Jucys–Murphy image identity (Lemma 6.6), and the dimension of the spectral corner (Corollary 5.8 with Lemma 5.9). The last two are stated in algebraic forms that avoid taking the sought map or category O\mathcal OO as extra hypotheses.

What the result supplies

The theorem gives an explicit dictionary between diagram generators and operators in a cyclotomic quotient, including the change from the quotient's loop parameter ω0=m+n\omega_0=m+nω0​=m+n to the source parameter −δ-\delta−δ Sartori–Stroppel, Remark 6.3. It also identifies the whole corner, rather than merely producing a homomorphism into it. This lets properties of the larger algebra be transported to the diagram algebra when their hypotheses are met. The paper discusses graded cellular and Koszul consequences, and an oriented-category extension, separately Sartori–Stroppel, Remarks 6.4–6.5.

A formal development must establish that the presentations do not collapse, that every required oriented relation survives the generator assignment, and that the image fills precisely fAffAffAf. The milestone statements isolate several of the authors' central identities, while the final theorem retains injectivity and surjectivity. A proof would also create reusable infrastructure for quotient presentations with source and target idempotents and for inverses in idempotent corners.

The central difficulty

The formulas for the images involve inverses and square roots after a spectral truncation, and the quotient need not be semisimple. Substituting a scalar directly into a formal power-series identity does not justify the key contraction identity in that setting; the authors explicitly point out this obstruction Sartori–Stroppel, Remark 6.3. The oriented relations also contain dotted symbols that stand for several local choices, so a single equation between global sums does not fully state every relation. Finally, proving that the map exists is only one part of the theorem: its image and kernel require separate control.

Formalization scope

Lean represents both algebras as quotients of free complex algebras by the paper's generators and relations. An orientation is a Boolean sequence with exactly rrr true entries; true denotes an upward arrow. Positions start at zero in Lean, so a paper index kkk corresponds to Lean index k−1k-1k−1. Relations involving dotted generators are imposed for every composable choice with the same source and target. Invalid orientation choices are zero. The affine quotient uses ek2=ω0eke_k^2=\omega_0e_kek2​=ω0​ek​; the source of the theorem uses ek2=−δeke_k^2=-\delta e_kek2​=−δek​. The recurrence fixes every ωj\omega_jωj​ rather than leaving later parameters free.

The ηk\eta_kηk​, intermediate fkf_kfk​, and fff are defined from the polynomial spectral-idempotent characterization, rather than supplied as arbitrary hypotheses. An inverse of xuxuxu is selected by its two inverse equations in the corner uAuuAuuAu. Square roots are not unique in the source, so each theorem holds for every family RkR_kRk​ of polynomial square roots satisfying Rk2=bkfR_k^2=b_kfRk2​=bk​f; the factorized form Qk=Rk+1Rk−1Q_k=R_{k+1}R_k^{-1}Qk​=Rk+1​Rk−1​ is used. A vacuous hypothesis, a collapsed or under-related presentation, a freely chosen fff, or a merely linear map does not express this target. The conclusion requires multiplication, the specified generator values, injectivity, and range exactly fAffAffAf.

The category O\mathcal OO construction and the paper's later applications are outside this mission. The needed reusable pieces are the finite-dimensional algebra lemma, oriented presentation machinery, polynomial spectral idempotents, and corner algebra operations. Contributions that strengthen the supporting algebraic library or close any of the stated milestones are within scope.

Selected references

  • Antonio Sartori and Catharina Stroppel, Walled Brauer algebras as idempotent truncations of level 2 cyclotomic quotients, arXiv preprint, 2015. arXiv:1411.2771v2.
8 thms0 active usersReviewed
Algebra·Captain: mikedeng1

The Representation Theory of Cyclotomic BMW Algebras: Semisimplicity Criterion for the Cyclotomic Birman–Wenzl Algebra B_{r,n} over a FieldResearch Paper

Motivation

The Birman–Murakami–Wenzl (BMW) algebras were introduced by Birman and Wenzl and by Murakami to study the Kauffman link polynomial; they are deformations of the Brauer algebras and act as centralizer algebras of quantum groups of orthogonal and symplectic type on tensor powers of the natural module. Their cyclotomic generalization was introduced by Häring-Oldenburg (J. Pure Appl. Algebra 161 (2001)) as the BMW analogue of the Ariki–Koike algebras, the cyclotomic Hecke algebras of type G(r,1,n)G(r,1,n)G(r,1,n). Representation theorists care about when such an algebra is semisimple: in the semisimple case its representation theory is governed by a combinatorial index set (multipartitions with extra data), and outside it one needs decomposition numbers.

For the Ariki–Koike algebra the answer is Ariki's criterion (J. Algebra 169 (1994)). For the ordinary BMW algebra (r=1r = 1r=1) it is due to Rui and Si (J. Reine Angew. Math., [17] in the paper). Rui and Xu proved that the cyclotomic BMW algebra is free and cellular for odd rrr under their u-admissibility conditions (arXiv:0801.0465), and Goodman and Hauschild Mosley (arXiv:math/0612064) and Yu (PhD thesis, Sydney, 2007) obtained freeness by other routes. Rui and Si (arXiv:0807.4149, 2008) extended cellularity to even rrr and, by computing Gram determinants of cell modules, gave a semisimplicity criterion for all r≥2r \ge 2r≥2. That criterion is the goal of this mission.

Setting

Let FFF be a field and r,nr, nr,n positive integers. Fix non-zero q,u1,…,ur∈Fq, u_1, \dots, u_r \in Fq,u1​,…,ur​∈F with δ=q−q−1≠0\delta = q - q^{-1} \ne 0δ=q−q−1=0 and the uiu_iui​ pairwise distinct, and a sign ε∈{1,−1}\varepsilon \in \{1, -1\}ε∈{1,−1}. Put α=ε\alpha = \varepsilonα=ε if rrr is odd and α=εq−ε\alpha = \varepsilon q^{-\varepsilon}α=εq−ε if rrr is even, ϱ−1=α∏ℓuℓ\varrho^{-1} = \alpha\prod_\ell u_\ellϱ−1=α∏ℓ​uℓ​, γr(x)=1\gamma_r(x) = 1γr​(x)=1 (rrr odd), −x-x−x (rrr even), and

γi=(γr(ui)+δ−1ϱ(ui2−1)∏j≠iuj)∏j≠iuiuj−1ui−uj,ωa=∑j=1rujaγj(a∈Z).\gamma_i = \Big(\gamma_r(u_i) + \delta^{-1}\varrho(u_i^2-1)\prod_{j \ne i}u_j\Big)\prod_{j\ne i}\frac{u_iu_j - 1}{u_i - u_j},\qquad \omega_a = \sum_{j=1}^r u_j^a\gamma_j\quad(a\in\mathbb Z).γi​=(γr​(ui​)+δ−1ϱ(ui2​−1)j=i∏​uj​)j=i∏​ui​−uj​ui​uj​−1​,ωa​=j=1∑r​uja​γj​(a∈Z).

These are the u-admissible parameters of Assumption 2.2. The cyclotomic Birman–Wenzl algebra Br,n\mathscr B_{r,n}Br,n​ is the unital associative FFF-algebra generated by Ti,EiT_i, E_iTi​,Ei​ (1≤i<n1 \le i < n1≤i<n) and Xj±1X_j^{\pm1}Xj±1​ (1≤j≤n1 \le j \le n1≤j≤n) subject to the relations a)–k) of Definition 2.1: among them the Kauffman skein relation Ti2−δTi+δϱEi=1T_i^2 - \delta T_i + \delta\varrho E_i = 1Ti2​−δTi​+δϱEi​=1, the braid relations, Ei2=ω0EiE_i^2 = \omega_0E_iEi2​=ω0​Ei​, the skein relations between TiT_iTi​ and Xi,Xi+1X_i, X_{i+1}Xi​,Xi+1​, the unwrapping relations E1X1aE1=ωaE1E_1X_1^aE_1 = \omega_aE_1E1​X1a​E1​=ωa​E1​ for all a∈Za \in \mathbb Za∈Z, the tangle and untwisting relations, EiXiXi+1=Ei=XiXi+1EiE_iX_iX_{i+1} = E_i = X_iX_{i+1}E_iEi​Xi​Xi+1​=Ei​=Xi​Xi+1​Ei​, and the cyclotomic relation (X1−u1)⋯(X1−ur)=0(X_1 - u_1)\cdots(X_1 - u_r) = 0(X1​−u1​)⋯(X1​−ur​)=0. The Ariki–Koike algebra Hr,n\mathscr H_{r,n}Hr,n​ is generated by g1,…,gn−1,y1,…,yng_1, \dots, g_{n-1}, y_1, \dots, y_ng1​,…,gn−1​,y1​,…,yn​ with (gi−q)(gi+q−1)=0(g_i - q)(g_i + q^{-1}) = 0(gi​−q)(gi​+q−1)=0, the braid relations, yi+1=giyigiy_{i+1} = g_iy_ig_iyi+1​=gi​yi​gi​, commuting yjy_jyj​ and (y1−u1)⋯(y1−ur)=0(y_1 - u_1)\cdots(y_1 - u_r) = 0(y1​−u1​)⋯(y1​−ur​)=0.

Write o(q2)>no(q^2) > no(q2)>n for "q2k≠1q^{2k} \ne 1q2k=1 for 1≤k≤n1 \le k \le n1≤k≤n", and set

Qr,ϱ={{−εq,εq−1},2∤r,{−qε,qε},2∣r,Sr,ϱ={⋃k=3n{±q3−k,±qk−3,εq3−2k,−εq2k−3},2∤r,⋃k=3n{±q3−k,±qk−3,±q(2k−3)ε},2∣r.Q_{r,\varrho} = \begin{cases}\{-\varepsilon q, \varepsilon q^{-1}\}, & 2 \nmid r,\\ \{-q^\varepsilon, q^\varepsilon\}, & 2 \mid r,\end{cases}\qquad S_{r,\varrho} = \begin{cases}\bigcup_{k=3}^n\{\pm q^{3-k}, \pm q^{k-3}, \varepsilon q^{3-2k}, -\varepsilon q^{2k-3}\}, & 2 \nmid r,\\ \bigcup_{k=3}^n\{\pm q^{3-k}, \pm q^{k-3}, \pm q^{(2k-3)\varepsilon}\}, & 2 \mid r.\end{cases}Qr,ϱ​={{−εq,εq−1},{−qε,qε},​2∤r,2∣r,​Sr,ϱ​={⋃k=3n​{±q3−k,±qk−3,εq3−2k,−εq2k−3},⋃k=3n​{±q3−k,±qk−3,±q(2k−3)ε},​2∤r,2∣r.​

Formalization targets

Goal: Theorem 6.8

For n≥2n \ge 2n≥2 and r≥2r \ge 2r≥2:

(a) if uiuj=1u_iu_j = 1ui​uj​=1 for some i≠ji \ne ji=j, or ui∈Qr,ϱu_i \in Q_{r,\varrho}ui​∈Qr,ϱ​ for some iii, then Br,n\mathscr B_{r,n}Br,n​ is not semisimple;

(b) if uiuj≠1u_iu_j \ne 1ui​uj​=1 for all i≠ji \ne ji=j and no uiu_iui​ lies in Qr,ϱQ_{r,\varrho}Qr,ϱ​, then

Br,2 semisimple  ⟺  o(q2)>2 and (uiuj−1=q2d, i<j⇒∣d∣≥2),\mathscr B_{r,2}\ \text{semisimple} \iff o(q^2) > 2 \ \text{and}\ \big(u_iu_j^{-1} = q^{2d},\ i<j \Rightarrow |d| \ge 2\big),Br,2​ semisimple⟺o(q2)>2 and (ui​uj−1​=q2d, i<j⇒∣d∣≥2),

and for n≥4n \ge 4n≥4, Br,n\mathscr B_{r,n}Br,n​ is semisimple if and only if (a) o(q2)>no(q^2) > no(q2)>n; (b) ∣d∣≥n|d| \ge n∣d∣≥n whenever uiuj−1=q2du_iu_j^{-1} = q^{2d}ui​uj−1​=q2d, i<ji < ji<j; (c) ui∉Sr,ϱu_i \notin S_{r,\varrho}ui​∈/Sr,ϱ​ for all iii; (d) uiuj∉⋃k=3n{q4−2k,q2k−4}u_iu_j \notin \bigcup_{k=3}^n\{q^{4-2k}, q^{2k-4}\}ui​uj​∈/⋃k=3n​{q4−2k,q2k−4} for all i≠ji \ne ji=j. For n=3n = 3n=3 the goal states that semisimplicity implies (a)–(d).

The page states (b)(2) as an equivalence for all n≥3n \ge 3n≥3. The goal keeps only necessity at n=3n = 3n=3: Proposition 6.2 and Lemma 6.3 describe the same cell module of Br,3\mathscr B_{r,3}Br,3​ and list different exclusions (εq\varepsilon qεq against −εq−1-\varepsilon q^{-1}−εq−1 for odd rrr), and Sr,ϱS_{r,\varrho}Sr,ϱ​ at n=3n = 3n=3 contains neither.

Milestones

The rank dim⁡FBr,n=rn(2n−1)!!\dim_F\mathscr B_{r,n} = r^n(2n-1)!!dimF​Br,n​=rn(2n−1)!! (Theorem 2.14); the corner identities En−1Br,nEn−1=En−1Br,n−2E_{n-1}\mathscr B_{r,n}E_{n-1} = E_{n-1}\mathscr B_{r,n-2}En−1​Br,n​En−1​=En−1​Br,n−2​ (Lemma 3.1) and Ef,nBr,nEf,n=Br,n−2fEf,nE_{f,n}\mathscr B_{r,n}E_{f,n} = \mathscr B_{r,n-2f}E_{f,n}Ef,n​Br,n​Ef,n​=Br,n−2f​Ef,n​ (Corollary 3.2), corrected at the two-strand boundary by requiring some ωa≠0\omega_a\ne0ωa​=0; Proposition 6.1, which says that det⁡G1,∅=det⁡(ωa+b)a,b∈Nr\det G_{1,\emptyset} = \det(\omega_{a+b})_{a,b\in N_r}detG1,∅​=det(ωa+b​)a,b∈Nr​​ is non-zero exactly when uiuj≠1u_iu_j \ne 1ui​uj​=1 (i≠ji \ne ji=j) and no uiu_iui​ lies in Qr,ϱQ_{r,\varrho}Qr,ϱ​; Proposition 6.5(b), non-semisimplicity when that determinant vanishes; and the reduction "Hr,n\mathscr H_{r,n}Hr,n​ not semisimple ⇒Br,n\Rightarrow \mathscr B_{r,n}⇒Br,n​ not semisimple" from the proof of Theorem 6.8.

Significance

The criterion gives a closed list of parameter values at which the cyclotomic BMW algebra fails to be semisimple, and so, by Graham–Lehrer theory, the parameters at which all cell modules are simple and pairwise non-isomorphic. Outside that list the algebra's simple modules, their dimensions and the decomposition of tensor-type modules are given by the combinatorics of up-down tableaux; inside it, the list says which Gram determinants vanish and where decomposition numbers begin. For r=1r = 1r=1 it specializes to the semisimplicity criterion for BMW algebras, and its conditions (a)–(b) are Ariki's conditions for Hr,n\mathscr H_{r,n}Hr,n​.

As far as is known, none of this is machine-checked: Mathlib has no BMW, Brauer, Ariki–Koike or Iwahori–Hecke algebras and no cellular algebras. The mission formalizes the statements of the paper on a presentation of Br,n\mathscr B_{r,n}Br,n​ by generators and relations; a solver must supply a basis (or at least a dimension count), the cellular structure and the Gram-determinant computations. The paper leaves several of its proofs to the reader ("by arguments in the proof of [19, 5.41]"), so the formal proofs would also close gaps in the published argument.

Difficulty

Non-semisimplicity in part (a) needs only one degenerate bilinear form on a non-zero cell module, but showing that the module is non-zero, or that Br,n\mathscr B_{r,n}Br,n​ does not collapse, already needs a lower bound on the dimension of a presented algebra, which no rewriting argument provides: the rank rn(2n−1)!!r^n(2n-1)!!rn(2n−1)!! is proved in the literature through seminormal representations. Sufficiency in part (b) requires the Gram determinants of all cell modules Δ(f,λ)\Delta(f,\lambda)Δ(f,λ), which the paper obtains by a recursion over up-down tableaux; vanishing of a single factor must be traced through branching and Morita-type functors between levels n−2n-2n−2, n−1n-1n−1 and nnn. Ariki's criterion for Hr,n\mathscr H_{r,n}Hr,n​ is itself an external input.

Formalization scope

Br,n\mathscr B_{r,n}Br,n​ and Hr,n\mathscr H_{r,n}Hr,n​ are RingQuot quotients of FreeAlgebra F on finite generator types; indices are 0-based (the Lean T i is the paper's Ti+1T_{i+1}Ti+1​), and Xj−1X_j^{-1}Xj−1​ is a separate generator tied to XjX_jXj​ by relation a). The sign ε∈Z×\varepsilon \in \mathbb Z^\timesε∈Z× is the free parameter and α,ϱ,γi,ωa\alpha, \varrho, \gamma_i, \omega_aα,ϱ,γi​,ωa​ are computed from (q,u,ε)(q, u, \varepsilon)(q,u,ε), so the parameters are u-admissible by construction. The uiu_iui​ are pairwise distinct (Assumption 2.2(1) divides by ui−uju_i - u_jui​−uj​; the paper extends 2.2 to coincident uiu_iui​ through [19, 3.11]). "Semisimple" is Mathlib's IsSemisimpleRing. The condition o(q2)>no(q^2) > no(q2)>n is written as q2k≠1q^{2k} \ne 1q2k=1 for 1≤k≤n1 \le k \le n1≤k≤n, not through orderOf, which is 000 for elements of infinite order. Relations d) and g) both use ω0=∑jγj\omega_0 = \sum_j\gamma_jω0​=∑j​γj​, which equals the paper's 1−δ−1(ϱ−ϱ−1)1 - \delta^{-1}(\varrho - \varrho^{-1})1−δ−1(ϱ−ϱ−1) for these parameters; were they different, the presentation would collapse. The corrected corner milestones require n≥3n\ge3n≥3 or some ωa≠0\omega_a\ne0ωa​=0 in Lemma 3.1, and 2f<n2f<n2f<n or some ωa≠0\omega_a\ne0ωa​=0 in Corollary 3.2. At n=2n=2n=2, r=2r=2r=2, q=3q=3q=3, u=(2,12)u=(2,\tfrac12)u=(2,21​), and ε=1\varepsilon=1ε=1, every ωa\omega_aωa​ vanishes, exposing the missing condition. det⁡G1,∅\det G_{1,\emptyset}detG1,∅​ is expressed as the Hankel determinant of the ωa\omega_aωa​, and Br,m⊆Br,n\mathscr B_{r,m} \subseteq \mathscr B_{r,n}Br,m​⊆Br,n​ is the subalgebra generated by the first generators.

A presentation that collapses to the zero algebra would make the goal false in part (a) but trivially true elsewhere; the rank milestone rules this out. Replacing semisimplicity by non-vanishing of Gram determinants would turn the goal into a different theorem and is not acceptable.

A complete development needs a basis theorem for Br,n\mathscr B_{r,n}Br,n​, cellular algebras and Graham–Lehrer's semisimplicity criterion, the Ariki–Koike algebra and Ariki's criterion, and Gram-determinant computations. The cellular-algebra layer and the Ariki–Koike algebra are reusable beyond this mission; contributions of any of these, or of the milestones on their own, are welcome.

Selected references

  • H. Rui and M. Si, The representation theory of cyclotomic BMW algebras (printed title The representations of cyclotomic BMW algebras, II), preprint, 2008. https://arxiv.org/abs/0807.4149
  • H. Rui and J. Xu, The representations of cyclotomic BMW algebras, preprint, 2007. https://arxiv.org/abs/0801.0465
  • R. Häring-Oldenburg, Cyclotomic Birman–Murakami–Wenzl algebras, J. Pure Appl. Algebra 161 (2001), 113–144. https://doi.org/10.1016/S0022-4049(00)00100-6
  • S. Ariki, On the semi-simplicity of the Hecke algebra of (Z/rZ)≀Sn(\mathbb Z/r\mathbb Z)\wr\mathfrak S_n(Z/rZ)≀Sn​, J. Algebra 169 (1994), 216–225. https://doi.org/10.1006/jabr.1994.1280
  • J. J. Graham and G. I. Lehrer, Cellular algebras, Invent. Math. 123 (1996), 1–34. https://doi.org/10.1007/BF01232365
  • J. S. Birman and H. Wenzl, Braids, link polynomials and a new algebra, Trans. Amer. Math. Soc. 313 (1989), 249–273. https://doi.org/10.1090/S0002-9947-1989-0992598-X
8 thms0 active usersReviewed
AlgebraCombinatorics·Captain: mikedeng1

On Central Idempotents in the Brauer Algebra: The Equation uX = −u Is Uniquely Solvable in Orbit Sums with Coefficients in K, and 1 + X Is the Splitting Idempotent of J_n(ℓ)Research Paper

Motivation

The Brauer algebra packages pairings of points into an algebra whose multiplication records how pairings join when diagrams are stacked. It was introduced by Brauer as the centralizer algebra of orthogonal and symplectic groups acting on tensor powers of their natural representation (Brauer 1937), and is semisimple for generic δ\deltaδ (Wenzl 1988). It provides an example where diagram combinatorics and the choice of coefficient ring affect the same algebraic question. King, Martin, and Parker study central idempotents that separate the diagrams with few propagating lines from the remaining part of the algebra. Their construction works over a ring KKK designed to retain useful specializations of the parameter δ\deltaδ; a solution obtained only in the field Q(δ)\mathbb Q(\delta)Q(δ) would lose that arithmetic information. King–Martin–Parker, §1–§2

The mission concerns their main result, Theorem 14. It asks for the precise bridge between a diagram equation, coefficients constrained to KKK, and an idempotent that splits an ideal quotient. The source is the 2016 arXiv version; all page and result numbers here refer to that version. King–Martin–Parker, Theorem 14

Setting

A Brauer diagram on nnn strands is a perfect matching of two rows of nnn nodes. A strand joining the rows is a propagating line; a strand within one row is a horizontal arc. The set of all such diagrams is a basis for BnB_nBn​, here over F=Q(δ)F=\mathbb Q(\delta)F=Q(δ), where δ\deltaδ is an indeterminate. To multiply diagrams AAA and BBB, place AAA over BBB and identify the adjacent rows. The resulting outer pairing is another diagram CCC. If v(A,B)v(A,B)v(A,B) connected components close entirely in the middle row, then

AB=δv(A,B)C.AB=\delta^{v(A,B)}C.AB=δv(A,B)C.

This loop factor is part of the algebra, even when no loop is visible in a particular example. The span Jn(ℓ)J_n(\ell)Jn​(ℓ) of diagrams with at most ℓ\ellℓ propagating lines is a two-sided ideal. Simultaneous relabelling of the upper and lower rows by Sn\mathfrak S_nSn​ partitions diagrams into orbits OOO. The orbit sum DO=∑A∈OAD_O=\sum_{A\in O}ADO​=∑A∈O​A uses each original, unnormalised diagram once. King–Martin–Parker, §2.2–§2.3

The source works with a coefficient ring inside FFF:

K={f/g:f,g∈Z[δ], g monic, deg⁡f≤deg⁡g}.K=\{f/g:f,g\in\mathbb Z[\delta],\ g\text{ monic},\ \deg f\leq\deg g\}.K={f/g:f,g∈Z[δ], g monic, degf≤degg}.

For ℓ<n\ell<nℓ<n and even n−ℓn-\elln−ℓ, let uuu be the diagram with cups at (1,2),(3,4),…,(n−ℓ−1,n−ℓ)(1,2),(3,4),\ldots,(n-\ell-1,n-\ell)(1,2),(3,4),…,(n−ℓ−1,n−ℓ) in each row and vertical lines elsewhere. The paper uses the normalized element uˉ=δ−(n−ℓ)/2u\bar u=\delta^{-(n-\ell)/2}uuˉ=δ−(n−ℓ)/2u. Its required equation uˉX=−uˉ\bar uX=-\bar uuˉX=−uˉ has the same solutions as uX=−uuX=-uuX=−u over FFF. King–Martin–Parker, §2.2, equation (4), and §3, equation (5)

Formalization targets

The goal is Theorem 14. For orbit coefficients cOc_OcO​ and Xc=∑OcODOX_c=\sum_Oc_OD_OXc​=∑O​cO​DO​, it asserts a unique solution with cO∈Kc_O\in KcO​∈K to

uXc=−u.uX_c=-u.uXc​=−u.

It further asserts that φ=1+Xc\varphi=1+X_cφ=1+Xc​ satisfies φ−1∈Jn(ℓ)\varphi-1\in J_n(\ell)φ−1∈Jn​(ℓ) and φJn(ℓ)=Jn(ℓ)φ=0\varphi J_n(\ell)=J_n(\ell)\varphi=0φJn​(ℓ)=Jn​(ℓ)φ=0, and is central and idempotent. By Lemma 10, these properties identify the splitting of 0→Jn(ℓ)→Bn→Bn/Jn(ℓ)→00\to J_n(\ell)\to B_n\to B_n/J_n(\ell)\to00→Jn​(ℓ)→Bn​→Bn​/Jn​(ℓ)→0. King–Martin–Parker, Lemma 10 and Theorem 14

The milestones follow the paper's proof, in attack order: the general bimodule splitting criterion (Lemma 10); the basis of the Sn\mathfrak S_nSn​-centralizer ZSn(Bn)Z_{\mathfrak S_n}(B_n)ZSn​​(Bn​) given by the orbit sums (Lemma 11); the paper's assertion that for ℓ<n\ell<nℓ<n the splitting idempotent φn(ℓ)\varphi_n(\ell)φn​(ℓ) exists with coefficients in KKK (p. 8); the fact that uˉ\bar uuˉ fixes every diagram AtA_{\mathfrak t}At​ carrying its cups and maps no other diagram onto AtA_{\mathfrak t}At​ without a negative power of δ\deltaδ (Proposition 12); the independence of the linear equations in the ctc_{\mathfrak t}ct​, i.e. injectivity of X↦uXX\mapsto uXX↦uX on orbit-sum combinations (equation (6)); and membership of every solution coefficient in KKK (Proposition 13). Lemma 10 is reusable beyond diagram algebras. King–Martin–Parker, §3

Significance

The idempotent separates the ideal of diagrams with few propagating lines from the quotient as bimodules. Its centrality means this separation is compatible with multiplication on either side. Specifying the coefficients in KKK also preserves information about specializations of δ\deltaδ that a field-only existence statement would omit. The paper uses the resulting splitting idempotents in its discussion of primitive central idempotents and Brauer algebra representations. King–Martin–Parker, §4

The theorem is proved in the cited paper. This mission supplies a machine-checkable interface for its definitions and statements; the theorem proofs remain open. A completed development would verify the diagram multiplication, the ideal and orbit descriptions, the ring-theoretic splitting criterion, and the coefficient claim in one Lean environment. These components can support later formalizations involving diagram algebras or central idempotents.

Difficulty

An equation over FFF alone does not give a splitting idempotent with coefficients in KKK. The coefficient restriction is sensitive to the exact loop factor in multiplication and to whether orbit sums use normalized or original diagrams. The equations indexed by diagrams also need to determine the orbit coefficients uniquely, while the two-sided annihilation requirement concerns every element of the ideal, not just the displayed cup diagram. These are distinct claims in the source; checking one of them does not silently establish the others. King–Martin–Parker, §3

Formalization scope

Lean represents a diagram as a fixed-point-free involution on Fin n ⊕ Fin n; the left summand is the upper row. A three-row graph defines its composite and loop count. The Brauer space is the finite-dimensional function space Diagram n → RatFunc ℚ, with an explicit bilinear multiplication bmul; the statements do not assume a ring instance whose associativity has not yet been proved. The parameter is RatFunc.X, representing the indeterminate δ\deltaδ. The ideal is the FFF-span of exactly the basis diagrams with propagating count at most ℓ\ellℓ.

The goal fixes ℓ<n\ell<nℓ<n and additionally assumes even n−ℓn-\elln−ℓ (the existence milestone, stated on the page for every ℓ<n\ell<nℓ<n, does not). That parity condition makes the paper's displayed product uˉ1uˉ3⋯uˉn−ℓ−1\bar u_1\bar u_3\cdots\bar u_{n-\ell-1}uˉ1​uˉ3​⋯uˉn−ℓ−1​ well defined. Lean uses zero-based indices, so the paper's cup (2j+1,2j+2)(2j+1,2j+2)(2j+1,2j+2) is represented by (2j,2j+1)(2j,2j+1)(2j,2j+1). The equation uses unnormalised uuu but has the same solutions as the paper's normalized equation because δ\deltaδ is nonzero in FFF. Orbit sums are unnormalised, as in the printed Theorem 14. The centralizer ZSn(Bn)Z_{\mathfrak S_n}(B_n)ZSn​​(Bn​) is the set of elements commuting with every diagram having nnn propagating lines. The index set of the goal consists only of orbits of diagrams with at most ℓ\ellℓ propagating lines; it is equivalent to the paper's Spore-tableau index set, without introducing that classification into this mission.

The statement constructs φ\varphiφ as 1+Xc1+X_c1+Xc​ from coefficients ccc, rather than postulating an arbitrary splitting idempotent. Neither the ideal nor the equation is replaced with a degenerate special case: with J=BnJ=B_nJ=Bn​ condition (ii) would force φ=0\varphi=0φ=0, and with J=0J=0J=0 it would be vacuous, so the ideal is exactly the span of the diagrams with at most ℓ\ellℓ propagating lines, and the loop factor δv(A,B)\delta^{v(A,B)}δv(A,B) is kept in the product. The main reusable infrastructure comprises finite diagrams, the loop-weighted product, orbit sums, and the generic bimodule splitting lemma. Contributions may establish these structural facts, the milestone statements, and the final theorem.

Selected references

  • O. H. King, P. P. Martin, and A. E. Parker, On central idempotents in the Brauer algebra, Journal of Algebra, 2016; arXiv:1609.01183v1. Page references use the arXiv version.
  • R. Brauer, On algebras which are connected with the semisimple continuous groups, Annals of Mathematics 38 (1937), 857–872. https://doi.org/10.2307/1968843
  • H. Wenzl, On the structure of Brauer's centralizer algebras, Annals of Mathematics 128 (1988), 173–193. https://doi.org/10.2307/1971466
8 thms0 active usersReviewed
AlgebraCombinatorics·Captain: mikedeng1

Representations of Temperley–Lieb Algebras: Closed Product Formula for the Gram Determinants of the Cell ModulesResearch Paper

Motivation

The Temperley–Lieb algebra is a family of algebras built from generators that satisfy local relations. It appears in diagrammatic representation theory and provides a tractable setting in which to study cellular bases and Gram determinants. A Gram determinant tests whether the bilinear form on a module is degenerate after specializing its parameter. Thus an explicit determinant formula records information about when a cell module may cease to be simple. Enyang's 2007 paper develops a Murphy-style basis, a commuting family of operators and an orthogonal basis, then gives a closed product for these determinants. The mission targets that final product, with the preceding statements retained as milestones. Enyang, 2007

The generic cellular-algebra definitions and spectral-projector construction have already been formalized in the published package associated with Mathas's work. This mission specializes those objects to the actual Temperley–Lieb presentation and computes its explicit polynomial factors. Mathas, 2006

Setting

Let R=Z[x]R=\mathbb Z[x]R=Z[x], with xxx an indeterminate. For n≥1n\ge1n≥1, the algebra An(x)A_n(x)An​(x) has generators e1,…,en−1e_1,\ldots,e_{n-1}e1​,…,en−1​ and relations

ei2=xei,eiei±1ei=ei,eiej=ejeiwhen ∣i−j∣≥2.e_i^2=xe_i,\qquad e_ie_{i\pm1}e_i=e_i,\qquad e_ie_j=e_je_i\quad\text{when }|i-j|\ge2.ei2​=xei​,ei​ei±1​ei​=ei​,ei​ej​=ej​ei​when ∣i−j∣≥2.

A shape at level nnn is a pair λ=(f,n−2f)\lambda=(f,n-2f)λ=(f,n−2f) with 0≤2f≤n0\le2f\le n0≤2f≤n. A standard up–down tableau of shape λ\lambdaλ is a path of shapes from (0,0)(0,0)(0,0) to λ\lambdaλ: at each level its cup count stays fixed or increases by one. Write Tn(λ)\mathfrak T_n(\lambda)Tn​(λ) for the set of these paths and dim⁡(λ)=∣Tn(λ)∣\dim(\lambda)=|\mathfrak T_n(\lambda)|dim(λ)=∣Tn​(λ)∣. The paper orders paths by comparing their cup counts at every level; u▹tu\triangleright tu▹t means strict pointwise dominance. Enyang, §§1–2

For each shape, mλm_\lambdamλ​ is the product e1e3⋯e2f−1e_1e_3\cdots e_{2f-1}e1​e3​⋯e2f−1​. Ordered generator products wi,jw_{i,j}wi,j​ and their tableau products vt,vt∗v_t,v_t^*vt​,vt∗​ give the Murphy elements muv=vu∗mλvvm_{uv}=v_u^*m_\lambda v_vmuv​=vu∗​mλ​vv​. Lemma 2.3 asserts that they form a cellular basis. The cellular multiplication law computes a Gram form on each cell module CλC^\lambdaCλ. The determinant det⁡(λ)\det(\lambda)det(λ) is the determinant of this computed form in the Murphy basis; it is not assigned the desired closed product by definition. Enyang, Lemma 2.3 and §5

The Jucys–Murphy elements TkT_kTk​ form a commuting family in An(x)A_n(x)An​(x). Their triangular action on a tableau ttt has content rt(k)∈Rr_t(k)\in Rrt​(k)∈R. Over the fraction field F=Frac⁡(R)F=\operatorname{Frac}(R)F=Frac(R), the spectral product FtF_tFt​ produces a vector ft=mtFtf_t=m_tF_tft​=mt​Ft​. These vectors form an orthogonal basis, with squared norms ⟨ft,ft⟩\langle f_t,f_t\rangle⟨ft​,ft​⟩ whose product equals det⁡(λ)\det(\lambda)det(λ). The auxiliary polynomials are fixed by p0=0p_0=0p0​=0, p1=1p_1=1p1​=1 and pi+1=xpi−pi−1p_{i+1}=xp_i-p_{i-1}pi+1​=xpi​−pi−1​. Enyang, §§3–5

Formalization targets

The milestones follow the paper's own sequence: the Murphy cellular basis (Lemma 2.3) and the neighbour criterion for the elements vtv_tvt​ (Corollary 2.2); the basic properties of TiT_iTi​ and ziz_izi​ (Lemma 3.1), the sum identities for pip_ipi​ (Lemma 3.4), the action of ziz_izi​ on the generator of a cell module (Lemma 3.5), the triangular action of TkT_kTk​ and zkz_kzk​ with diagonals rt(k)r_t(k)rt​(k) and zt(k)=pipk−i+1z_t(k)=p_ip_{k-i+1}zt​(k)=pi​pk−i+1​ (Corollaries 3.6, 3.7) and the alternative recursion for Ti+1T_{i+1}Ti+1​ (Corollary 4.2); the seminormal basis (Lemma 4.4), the diagonal coefficients of eke_kek​ on it (Corollary 4.6) and the norm ratio of neighbouring seminormal vectors (Lemma 4.7); the identification of det⁡(λ)\det(\lambda)det(λ) with ∏t⟨ft,ft⟩\prod_t\langle f_t,f_t\rangle∏t​⟨ft​,ft​⟩ (§5) and the branching law (Lemma 5.1). The goal is Lemma 5.2, the closed product. When n>2fn>2fn>2f,

det⁡(f,n−2f)=∏i=1f(pn−2f+i+1pi)dim⁡(f−i,n−2f+2i).\det(f,n-2f)=\prod_{i=1}^{f}\left(\frac{p_{n-2f+i+1}}{p_i}\right)^{\dim(f-i,n-2f+2i)}.det(f,n−2f)=i=1∏f​(pi​pn−2f+i+1​​)dim(f−i,n−2f+2i).

When n=2fn=2fn=2f, the predecessor shape μ=(f−1,1)\mu=(f-1,1)μ=(f−1,1) lies at level n−1n-1n−1, and the target is

det⁡(f,0)=det⁡(μ)xdim⁡(f,0).\det(f,0)=\det(\mu)x^{\dim(f,0)}.det(f,0)=det(μ)xdim(f,0).

The quotient in the first formula is evaluated in FFF; the second equality is in RRR. Both clauses occur in the printed Lemma 5.2. Enyang, Lemma 5.2, p. 20

Significance

The result turns a determinant of a matrix indexed by all tableaux of a shape into a finite product of explicit recursive polynomials. It makes the factors that can vanish under specialization visible and lets later representation-theoretic questions be asked in terms of those factors. The branching law supplies a compatible description across levels of the algebra tower. Enyang, §5

The paper proves the formulas mathematically. The statements in this mission are compiled Lean targets with proof holes; they do not yet constitute machine-checked proofs of Enyang's results. A complete development would supply the cellular-basis theorem, establish the paper's explicit content action, connect the integral and field Gram forms, and prove the norm and determinant formulas. The quotient presentation, tableau encoding and cellular package can be reused for related specializations and diagram-algebra calculations.

Difficulty

The generator relations make local products manageable, but a determinant does not expose its factors directly. The cell module is indexed by whole paths, and the operators act triangularly rather than diagonally in the Murphy basis. Passing to seminormal vectors changes the basis while keeping track of the bilinear form; this requires the local coefficient and norm statements, not just the commutation of the operators. At the boundary n=2fn=2fn=2f, one of the two possible predecessor shapes disappears, so a uniform two-factor branching formula would introduce a nonexistent module. Enyang, §§3–5

Formalization scope

The Lean algebra is a quotient of a free algebra by all three defining relation families. It uses one-based generators; an out-of-range generator is set to zero, and every stated result guards its meaningful indices. Shapes are finite cup-count indices and tableaux are finite paths with a pointwise order. The unrolled products defining vtv_tvt​ and vt∗v_t^*vt∗​ follow the page's multiplication order. The integral statements use R=Z[x]R=\mathbb Z[x]R=Z[x]; seminormal vectors use F=Frac⁡(R)F=\operatorname{Frac}(R)F=Frac(R) and the published cellular-algebra package's coordinate-space cell module and Gram form. Existence of the required Murphy and JM data is asserted, preventing a theorem about an empty class of such data from passing vacuously.

Several printed details need explicit treatment. The page-2 endpoint e2f+1e_{2f+1}e2f+1​ in mλm_\lambdamλ​ is replaced by e2f−1e_{2f-1}e2f−1​, giving exactly fff factors and m(0,n)=1m_{(0,n)}=1m(0,n)​=1. The missing z0z_0z0​ is zero by the paper's sum formula. The rs(k)r_s(k)rs​(k) denominator in the definition of FtF_tFt​ means rt(k)r_t(k)rt​(k). Lemma 4.4(3) says “orthonormal,” but later nonunit norms show that the valid conclusion is pairwise orthogonality. The direction of ▹\triangleright▹ in the printed definition of gλ,μg_{\lambda,\mu}gλ,μ​ and Lemma 5.2 is reversed relative to the formula's denominator and the proof; the target product runs over shapes with fewer cups. Identity (3.3) of Lemma 3.4 is false as printed (at i=j=2i=j=2i=j=2, k=1k=1k=1 it reads x3−2x=xx^3-2x=xx3−2x=x); the target is the form its proof derives, with last factor pj+k−1p_{j+k-1}pj+k−1​. The goal and the cellular milestones take n≥1n\ge1n≥1, since the paper's tower starts at A1(x)=RA_1(x)=RA1​(x)=R; for n=0n=0n=0 the even case of Lemma 5.2 would refer to the non-existent shape (−1,1)(-1,1)(−1,1). Enyang, pp. 2, 10, 17, 19–20

Useful contributions include proofs of the algebraic Murphy-basis statement, the content formulas, the seminormal basis properties, the integral-to-field determinant identification and the branching law. Each has a separate mission item and source index.

Selected references

  • John Enyang, Representations of Temperley–Lieb Algebras, arXiv preprint arXiv:0710.3218v1, 2007. arXiv
  • Andrew Mathas, Seminormal forms and Gram determinants for cellular algebras, arXiv preprint arXiv:math/0604108v4, 2006. arXiv
18 thms0 active usersReviewed
AlgebraCombinatorics·Captain: mikedeng1

On the Blocks of the Walled Brauer Algebra: In Characteristic Zero, Two Simple Modules Lie in the Same Block Exactly When Their Weights Lie in One Orbit of the ρ(δ)-Shifted Dot Action of Σ_{r+s}Research Paper

Motivation

Let VVV be the natural representation of GLnGL_nGLn​ and V∗V^*V∗ its dual. The walled Brauer algebra Br,s(δ)B_{r,s}(\delta)Br,s​(δ) describes the endomorphisms of mixed tensor space V⊗r⊗V∗⊗sV^{\otimes r}\otimes V^{*\otimes s}V⊗r⊗V∗⊗s that commute with GLnGL_nGLn​, through mixed Schur–Weyl duality. It was introduced in this role independently by Turaev, Koike and Benkart–Chakrabarti–Halverson–Leduc–Lee–Stroomer around 1989–1994. It is also a diagram algebra in its own right, defined for every parameter δ\deltaδ in the ground field. Its representation theory controls the decomposition of mixed tensor space, and it serves as a test case for the cellular and "tower of recollement" methods used for the Brauer and partition algebras.

Cox, De Visscher, Doty and Martin (arXiv:0709.0851, 2007) determine the blocks of Br,s(δ)B_{r,s}(\delta)Br,s​(δ) in characteristic zero. They also classify the semisimple walled Brauer algebras and give a linkage principle in positive characteristic. Their characteristic-zero answer has the shape familiar from algebraic Lie theory: blocks are orbits of a Weyl group, Σr+s\Sigma_{r+s}Σr+s​, under a dot action shifted by a vector ρ(δ)\rho(\delta)ρ(δ). The analogous result for the ordinary Brauer algebra, with the Weyl group of type DDD, is due to Cox, De Visscher and Martin (arXiv:math/0601387, 2009).

Setting

Fix r,s∈Nr,s\in\mathbb Nr,s∈N, an algebraically closed field kkk and δ∈k\delta\in kδ∈k. A Brauer diagram on n=r+sn=r+sn=r+s strands pairs up 2n2n2n nodes, nnn northern and nnn southern. Lines joining the two rows are propagating lines; the others are arcs. Put a wall after the first rrr nodes of each row. A diagram is walled if no propagating line crosses the wall and every arc does. Two diagrams are multiplied by stacking the first above the second. Each closed loop formed in the middle row is removed and contributes a factor δ\deltaδ. The walled Brauer algebra Br,s(δ)B_{r,s}(\delta)Br,s​(δ) is the kkk-span of the walled diagrams with this product. The permutation diagrams span the group algebra kΣr,sk\Sigma_{r,s}kΣr,s​ of Σr,s=Σr×Σs\Sigma_{r,s}=\Sigma_r\times\Sigma_sΣr,s​=Σr​×Σs​.

The algebra is cellular. Its cell modules Δr,s(λL,λR)\Delta_{r,s}(\lambda^L,\lambda^R)Δr,s​(λL,λR) are indexed by bipartitions (λL,λR)(\lambda^L,\lambda^R)(λL,λR) with λL⊢r−t\lambda^L\vdash r-tλL⊢r−t and λR⊢s−t\lambda^R\vdash s-tλR⊢s−t, for 0≤t≤min⁡(r,s)0\le t\le\min(r,s)0≤t≤min(r,s). Here ttt is the number of northern arcs. In characteristic zero the simple modules Lr,s(λL,λR)L_{r,s}(\lambda^L,\lambda^R)Lr,s​(λL,λR), the heads of the cell modules, are indexed by the same set. The one exception is (∅,∅)(\emptyset,\emptyset)(∅,∅), which is excluded when δ=0\delta=0δ=0 and 0<r=s0<r=s0<r=s. At r=s=0r=s=0r=s=0, the algebra is kkk and its unique simple has the empty bipartition as label. Two simple modules lie in the same block when no central idempotent of Br,s(δ)B_{r,s}(\delta)Br,s​(δ) separates them.

A box in row iii and column jjj of a Young diagram has content j−ij-ij−i. Two bipartitions are δ\deltaδ-balanced when the boxes of λL∖μL\lambda^L\setminus\mu^LλL∖μL can be paired with those of λR∖μR\lambda^R\setminus\mu^RλR∖μR, and the boxes of μL∖λL\mu^L\setminus\lambda^LμL∖λL with those of μR∖λR\mu^R\setminus\lambda^RμR∖λR, so that the contents in each pair sum to −δ-\delta−δ.

For δ∈Z\delta\in\mathbb Zδ∈Z, the weight lattice is Xr,s=Zr+sX_{r,s}=\mathbb Z^{r+s}Xr,s​=Zr+s, with coordinates indexed by −r,…,−1,1,…,s-r,\dots,-1,1,\dots,s−r,…,−1,1,…,s. Put

ρ(δ)=(r,r−1,…,1; δ,δ−1,…,δ−s+1).\rho(\delta)=(r,r-1,\dots,1;\ \delta,\delta-1,\dots,\delta-s+1).ρ(δ)=(r,r−1,…,1; δ,δ−1,…,δ−s+1).

The group W=Σr+sW=\Sigma_{r+s}W=Σr+s​ permutes the coordinates and acts by the dot action w.λ=w(λ+ρ)−ρw.\lambda=w(\lambda+\rho)-\rhow.λ=w(λ+ρ)−ρ. A bipartition is sent to the weight (λˉL;λR)=(−λrL,…,−λ1L;λ1R,…,λsR)(\bar\lambda^L;\lambda^R)=(-\lambda^L_r,\dots,-\lambda^L_1;\lambda^R_1,\dots,\lambda^R_s)(λˉL;λR)=(−λrL​,…,−λ1L​;λ1R​,…,λsR​).

Formalization targets

Goal: Corollary 8.4 (characteristic zero)

For kkk of characteristic 000, δ∈Z\delta\in\mathbb Zδ∈Z, and labels λ,μ\lambda,\muλ,μ of simple Br,s(δ)B_{r,s}(\delta)Br,s​(δ)-modules,

Lr,s(λ) and Lr,s(μ) lie in the same block  ⟺  (λˉL;λR)=w.(μˉL;μR) for some w∈Σr+s.L_{r,s}(\lambda)\ \text{and}\ L_{r,s}(\mu)\ \text{lie in the same block}\iff (\bar\lambda^L;\lambda^R)=w.(\bar\mu^L;\mu^R)\ \text{for some } w\in\Sigma_{r+s}.Lr,s​(λ) and Lr,s​(μ) lie in the same block⟺(λˉL;λR)=w.(μˉ​L;μR) for some w∈Σr+s​.

The balanced classification: Corollary 7.7

For every δ∈k\delta\in kδ∈k in characteristic 000, two labels are in the same block if and only if they are δ\deltaδ-balanced, and each block contains a unique minimal weight.

Supporting results

The milestones, in proof order, are the following.

  • Lemma 4.1: the action of Tr,s=∑Ei,jT_{r,s}=\sum E_{i,j}Tr,s​=∑Ei,j​ on cell modules.
  • Theorem 4.2 and Corollary 4.3: content conditions on composition factors and on blocks.
  • Theorem 6.2: two-box homomorphisms, Hom⁡(Δ(λ),Δ(μ))≅k\operatorname{Hom}(\Delta(\lambda),\Delta(\mu))\cong kHom(Δ(λ),Δ(μ))≅k.
  • Proposition 7.1 and Corollary 7.3: box pairings and the necessity of balance.
  • Lemma 8.1: the counting form of balance.
  • Propositions 8.2 and 8.3: orbit ⇔\Leftrightarrow⇔ balance.

Significance

The goal reduces the block structure of Br,s(δ)B_{r,s}(\delta)Br,s​(δ) to a combinatorial test. Two weights are in the same block exactly when their ρ(δ)\rho(\delta)ρ(δ)-shifted vectors are permutations of each other, a check that takes linear time. Within a block there is a unique minimal weight. The Weyl-group description is the starting point for the paper's positive-characteristic linkage principle, in which the affine Weyl group WpW_pWp​ takes the place of WWW. It also places the walled Brauer algebra alongside the Brauer algebra and the algebraic groups as an algebra whose characteristic-zero blocks are governed by a reflection group.

The result has been proved since 2007. As far as is known, no part of it has been machine-checked. This mission produces three things.

  1. A diagram-level definition of Br,s(δ)B_{r,s}(\delta)Br,s​(δ) for every δ\deltaδ, together with its cell and simple modules. This layer is reusable for the Brauer, partition and Temperley–Lieb algebras.
  2. A content calculus for bipartitions.
  3. Formal proofs of the paper's chain of results. The combinatorial half (Lemma 8.1, Propositions 8.2 and 8.3) is independent of the algebra, and it can be completed first.

Difficulty

The obvious first idea is to treat the walled Brauer algebra through mixed Schur–Weyl duality, as a centraliser algebra. That gives only a quotient of Br,s(δ)B_{r,s}(\delta)Br,s​(δ) when δ\deltaδ is a small integer, and nothing at all when δ\deltaδ is not a natural number, so it cannot reach the full statement. The obstruction is the non-semisimple regime (δ∈Z\delta\in\mathbb Zδ∈Z, ∣δ∣≤r+s−2|\delta|\le r+s-2∣δ∣≤r+s−2), where cell modules have nontrivial composition factors.

The necessary condition is a scalar computation, obtained by evaluating a central-like element on cell modules. Sufficiency is the hard direction. It needs explicit nonzero homomorphisms between cell modules that differ by a removable strip of boxes, built from the two-box homomorphisms of Theorem 6.2 by an induction using restriction and induction along the tower Br−1,s⊂Br,sB_{r-1,s}\subset B_{r,s}Br−1,s​⊂Br,s​. The final step, matching balance with dot orbits, is a combinatorial argument on edges of skew diagrams in which the shift by ρ(δ)\rho(\delta)ρ(δ) must be tracked exactly.

Formalization scope

  1. Characteristic. The goal and the algebraic milestones take kkk algebraically closed of characteristic 000; the goal takes δ∈Z\delta\in\mathbb Zδ∈Z, cast into kkk. The paper's hypothesis "Σ\SigmaΣ-semisimple" also allows characteristic p>max⁡(r,s)p>\max(r,s)p>max(r,s). That case is not posed, because there the integer in ρ(δ)\rho(\delta)ρ(δ) is not determined by δ∈k\delta\in kδ∈k. Lemma 8.1 and Propositions 8.2 and 8.3 are purely combinatorial, with δ∈Z\delta\in\mathbb Zδ∈Z.
  2. The algebra. Br,s(δ)B_{r,s}(\delta)Br,s​(δ) has the δ#loops\delta^{\#\text{loops}}δ#loops product and is never the centraliser algebra. It is realised as its faithful left regular representation: the subalgebra of matrices on the walled diagrams generated by left multiplication by each diagram, which is isomorphic to the paper's algebra.
  3. Cell modules. In characteristic 000 they are modelled as Br,sXt yλB_{r,s}X_t\,y_\lambdaBr,s​Xt​yλ​ modulo Jt+1J_{t+1}Jt+1​. Here XtX_tXt​ is the diagram with arcs r−ur-ur−u — r+1+ur+1+ur+1+u (u<tu<tu<t), yλy_\lambdayλ​ is the Young symmetrizer ∑σ∈C∑τ∈Rsgn⁡(σ)στ\sum_{\sigma\in C}\sum_{\tau\in R}\operatorname{sgn}(\sigma)\sigma\tau∑σ∈C​∑τ∈R​sgn(σ)στ on the propagating nodes, and Jt+1J_{t+1}Jt+1​ is the span of the diagrams with more than ttt arcs.
  4. Simple modules and blocks. Simple modules are heads, i.e. quotients by the intersection of the maximal submodules. "Same block" is the central-idempotent criterion on simple modules.
  5. Combinatorial conventions. Young diagrams are Mathlib's YoungDiagram, with 0-based cells; contents are unchanged by the shift. Nodes are Fin (r+s). The weight lattice is Fin r ⊕ Fin s → ℤ, where Sum.inl i is the paper's coordinate −(i+1)-(i+1)−(i+1) and Sum.inr j is coordinate j+1j+1j+1. WWW is Equiv.Perm (Fin r ⊕ Fin s).
  6. Statements over weights. These assume at most rrr and sss rows, so that the embedding loses nothing.

Two corrections to the printed text are adopted. Lemma 8.1 is stated "for all i∈Zi\in\mathbb Zi∈Z", and the reflection formula on p. 31 is read as sβ(λ)=λ−(λ,β)βs_\beta(\lambda)=\lambda-(\lambda,\beta)\betasβ​(λ)=λ−(λ,β)β, the transposition of two coordinates. The label exception in Theorem 2.7(iii) is applied only when r=s>0r=s>0r=s>0; its literal wording would incorrectly leave B0,0(0)=kB_{0,0}(0)=kB0,0​(0)=k without a simple module.

A trivializing formalization would replace the block relation by a relation defined through contents or orbits, or would replace the algebra by a centraliser algebra or by kΣr,sk\Sigma_{r,s}kΣr,s​. Both are excluded: the block relation is the algebra's own central-idempotent relation, and the algebra is the diagram algebra itself.

Contributions are welcome at every layer. Useful general infrastructure includes associativity of diagram composition, Young symmetrizers and Specht modules over a field of characteristic 000, and the Jucys–Murphy content formula. Proofs of the combinatorial Section 8 results and of the cellular structure of Br,s(δ)B_{r,s}(\delta)Br,s​(δ) are also needed.

Selected references

  • A. Cox, M. De Visscher, S. Doty, P. Martin, On the blocks of the walled Brauer algebra, preprint arXiv:0709.0851v1, 2007. https://arxiv.org/abs/0709.0851
  • A. Cox, M. De Visscher, P. Martin, The blocks of the Brauer algebra in characteristic zero, Represent. Theory 13 (2009). https://arxiv.org/abs/math/0601387
  • G. Benkart, M. Chakrabarti, T. Halverson, R. Leduc, C. Lee, J. Stroomer, Tensor product representations of general linear groups and their connections with Brauer algebras, J. Algebra 166 (1994). https://doi.org/10.1006/jabr.1994.1165
  • T. Halverson, Characters of the centralizer algebras of mixed tensor representations of GL(r,C) and the quantum group U_q(gl(r,C)), Pacific J. Math. 174 (1996). https://doi.org/10.2140/pjm.1996.174.359
  • J. J. Graham, G. I. Lehrer, Cellular algebras, Invent. Math. 123 (1996). https://doi.org/10.1007/BF01232365
17 thms0 active usersReviewed
Functional AnalysisHarmonic Analysis·Captain: mikedeng1

On the Characterization of Trace Class Representations and Schwartz Operators: A Unitary Lie Group Representation Is Trace Class iff Some (1 − dπ(Δ))⁻ᵏ Is Trace Class iff π(Cᵐ_c(G)) ⊆ B₁ for Some mResearch Paper

Why trace class representations matter

A strongly continuous unitary representation turns each group element into a unitary operator on a Hilbert space. Averaging those operators against a smooth, compactly supported function gives a bounded operator. When every such average is trace class, its trace defines a distribution character. The trace condition is therefore a bridge from an operator representation to a scalar distribution that can be studied by harmonic analysis. Van Dijk, Neeb, Salmasian and Zellner collect and extend characterizations of this condition for finite-dimensional Lie groups in their 2015 preprint. Their Theorem 1.3 is the target here.

The authors identify one of the equivalences, between the trace class property and a trace class negative power of the Nelson Laplacian, with an earlier statement of Cartier that had appeared without a proof or reference. Their further conclusion that one finite degree of differentiability suffices is presented as a new result; it gives finite order for the distribution character. These are claims about the cited preprint, not claims of prior machine verification. Source: introduction and Theorem 1.3, pp. 2–5.

The representation and its test functions

Let GGG be a finite-dimensional real Lie group, let HHH be a complex Hilbert space, and let π:G→U(H)\pi:G\to\mathrm U(H)π:G→U(H) be a strongly continuous unitary representation. Fix a left Haar measure μ\muμ. For a compactly supported continuous function f:G→Cf:G\to\mathbb Cf:G→C, the integrated operator is π(f)v=∫Gf(g)π(g)v dμ(g)\pi(f)v=\int_G f(g)\pi(g)v\,d\mu(g)π(f)v=∫G​f(g)π(g)vdμ(g) for each v∈Hv\in Hv∈H. Write Ccm(G)C^m_c(G)Ccm​(G) for compactly supported functions with mmm continuous derivatives and Cc∞(G)C^\infty_c(G)Cc∞​(G) for smooth compactly supported functions. Write B1(H)B_1(H)B1​(H) and B2(H)B_2(H)B2​(H) for the trace class and Hilbert–Schmidt ideals of bounded operators on HHH.

Choose a basis X1,…,XnX_1,\ldots,X_nX1​,…,Xn​ of the Lie algebra g\mathfrak gg. The Nelson Laplacian is Δ=∑j=1nXj2\Delta=\sum_{j=1}^n X_j^2Δ=∑j=1n​Xj2​ in the enveloping algebra. The derived representation applies dπ(Δ)d\pi(\Delta)dπ(Δ) on smooth vectors. Its closure D=dπ(Δ)‾D=\overline{d\pi(\Delta)}D=dπ(Δ)​ is a non-positive self-adjoint operator, and (1−D)−1(1-D)^{-1}(1−D)−1 is bounded. The choice of basis determines DDD, while the theorem holds for an arbitrary fixed basis. Source: §1 and proof of Theorem 1.3, pp. 3–4.

Formalization targets

Five equivalent conditions

The mission goal is the complete five-way equivalence in Theorem 1.3:

(i)  ∃m≥1: π(Ccm(G))⊆B1(H) and π:Ccm(G)→B1(H) is continuous;(ii)  π(Cc∞(G))⊆B1(H) and π:Cc∞(G)→B1(H) is continuous;(iii)  π(Cc∞(G))⊆B1(H);(iv)  π(Cc∞(G))⊆B2(H);(v)  ∃k≥1: (1−D)−k∈B1(H).\begin{aligned} \text{(i)}\;&\exists m\ge1:\ \pi(C^m_c(G))\subseteq B_1(H) \text{ and }\pi:C^m_c(G)\to B_1(H)\text{ is continuous};\\ \text{(ii)}\;&\pi(C^\infty_c(G))\subseteq B_1(H) \text{ and }\pi:C^\infty_c(G)\to B_1(H)\text{ is continuous};\\ \text{(iii)}\;&\pi(C^\infty_c(G))\subseteq B_1(H);\\ \text{(iv)}\;&\pi(C^\infty_c(G))\subseteq B_2(H);\\ \text{(v)}\;&\exists k\ge1:\ (1-D)^{-k}\in B_1(H). \end{aligned}(i)(ii)(iii)(iv)(v)​∃m≥1: π(Ccm​(G))⊆B1​(H) and π:Ccm​(G)→B1​(H) is continuous;π(Cc∞​(G))⊆B1​(H) and π:Cc∞​(G)→B1​(H) is continuous;π(Cc∞​(G))⊆B1​(H);π(Cc∞​(G))⊆B2​(H);∃k≥1: (1−D)−k∈B1​(H).​

The continuity clauses are part of the target; they cannot be inferred merely from the displayed inclusions. The milestone list includes the four nontrivial implications named in the paper's proof, its finite-order claim and the two displayed identities needed to connect the function and operator Laplacians. Source: Theorem 1.3 and proof, pp. 3–5.

What the equivalence gives

Condition (iii) is formulated using every smooth compactly supported function. Condition (i) replaces smoothness by one finite differentiability order that works uniformly for the representation, while retaining continuity into the trace class ideal. This is the form needed to conclude that the distribution character has finite order. Condition (v) tests the same property with one operator associated with the Lie algebra basis. These consequences are stated by the authors in the introduction and derived in §1. Source: pp. 2–5.

A formal proof would connect Lie-group differentiation, Bochner integration, Schatten ideals and the topology of compactly supported function spaces in one reusable development. The target and milestones here are open Lean statements: compilation checks their types, not their proofs. The previously posed Lemmas 1.1 and 1.2 of the same paper remain separate open platform items, and can be used by solvers without being re-posed here.

Where the difficulty lies

Showing that each π(f)\pi(f)π(f) is trace class does not immediately give continuity of f↦π(f)f\mapsto\pi(f)f↦π(f) in trace norm: operator-norm convergence alone does not control trace norm. The paper also needs to pass from trace class images of smooth functions to an operator statement about (1−D)−k(1-D)^{-k}(1−D)−k, and then back to functions with only finitely many derivatives. The obvious shortcut of treating dπ(Δ)d\pi(\Delta)dπ(Δ) as a bounded operator fails because its natural domain is the smooth-vector space and its closure is generally unbounded. The proof therefore has substantial analytic content beyond elementary ideal inclusions. Source: proof of Theorem 1.3, pp. 4–5.

Formalization scope

Lean represents GGG as a Hausdorff finite-dimensional Lie group modelled on a finite-dimensional real normed space EEE, equipped with a regular left Haar measure. It represents HHH as a complete complex inner-product space, without a separability assumption, and the representation by a strongly continuous monoid homomorphism into unitary continuous linear equivalences. The Lie algebra is the tangent space EEE at the identity, and the basis is an arbitrary finite basis of EEE. The quantifiers for mmm and kkk include explicit lower bounds 1≤m1\le m1≤m and 1≤k1\le k1≤k, matching the paper's convention N0=N∪{0}\mathbb N_0=\mathbb N\cup\{0\}N0​=N∪{0}.

Smooth vectors are exactly those whose orbit maps are smooth. The derived action is the manifold derivative of that orbit map at the identity, used only on smooth vectors. The function-side Laplacian uses right-invariant directional derivatives; its second-order sign agrees with the paper's action. A bounded operator RRR is characterized by R(w−dπ(Δ)w)=wR(w-d\pi(\Delta)w)=wR(w−dπ(Δ)w)=w on every smooth vector. The paper's essential self-adjointness and density statement identifies this RRR uniquely with (1−D)−1(1-D)^{-1}(1−D)−1; the theorem carries that characterization as a hypothesis. The hypothesis must be satisfiable, and the interpretation must remain the paper's resolvent.

Continuity on Ccm(G)C^m_c(G)Ccm​(G) is encoded by a bound on the Schatten 111-norm for every compact support set, with one constant for all functions supported there. For Cc∞(G)C^\infty_c(G)Cc∞​(G), the derivative order may depend on the compact set. The Schatten norms and classes are reused published definitions, with norm values in [0,∞][0,\infty][0,∞]; no alternate trace norm is introduced. Equation (1) is an identity of distributions tested against smooth compactly supported functions. A complete proof needs the smooth-vector integration lemma, Lie-group test-function approximation, the function-side parametrix and operator ideal estimates. These components are useful beyond this mission.

Selected references

  • G. van Dijk, K.-H. Neeb, H. Salmasian and C. Zellner, On the characterization of trace class representations and Schwartz operators, arXiv preprint, 2015, arXiv:1512.02451v1.
12 thms0 active usersReviewed
Harmonic AnalysisRandom Matrix Theory·Captain: mikedeng1

The Weingarten Calculus: Haar Integrals of Matrix Elements Are Given by a Basis of Invariant Tensors and Its Weingarten Matrix, I = A W A*Research Paper

Motivation

Integrals of polynomial functions of the matrix entries of a random element of a compact group appear across mathematical physics and probability: in random matrix theory (moments of Haar-distributed unitary, orthogonal and symplectic matrices), in lattice gauge theory (the link integrals of U(N)U(N)U(N) gauge fields), and in quantum information (averages of quantum channels over random unitaries, twirling, unitary designs). Weingarten first computed asymptotics of such integrals for U(N)U(N)U(N) in 1978 (Weingarten, J. Math. Phys. 19 (1978)); Collins (IMRN 2003) and Collins–Śniady (CMP 2006) gave exact formulas for the unitary, orthogonal and symplectic groups in terms of what is now called the Weingarten function.

The survey of Collins, Matsumoto and Novak, The Weingarten Calculus (arXiv:2109.14890), isolates the single statement underneath all these formulas, valid for every compact group and every finite-dimensional unitary representation. This mission formalizes that statement, Theorem 2.1 of the survey, together with the steps of its proof as written in Section 2.

Setting

Let GGG be a compact Hausdorff topological group and μ\muμ its Haar probability measure, the unique Borel probability measure on GGG invariant under left translation. Let H=CNH=\mathbb{C}^NH=CN with its standard orthonormal basis e1,…,eNe_1,\dots,e_Ne1​,…,eN​, and let U:G→U(N)U:G\to\mathrm U(N)U:G→U(N) be a continuous group homomorphism into the unitary N×NN\times NN×N matrices. Its matrix elements are Uxy(g)=⟨ex,U(g)ey⟩U_{xy}(g)=\langle e_x,U(g)e_y\rangleUxy​(g)=⟨ex​,U(g)ey​⟩, the (x,y)(x,y)(x,y) entry of U(g)U(g)U(g).

Fix d∈Nd\in\mathbb Nd∈N. A multi-index is a function i:[d]→[N]i:[d]\to[N]i:[d]→[N]; the set of them is Fun(d,N)\mathrm{Fun}(d,N)Fun(d,N). The Weingarten integrals of degree ddd are

Iij=∫G∏x=1dUi(x)j(x)(g) dμ(g),i,j∈Fun(d,N).I_{ij}=\int_G\prod_{x=1}^d U_{i(x)j(x)}(g)\,d\mu(g),\qquad i,j\in\mathrm{Fun}(d,N).Iij​=∫G​x=1∏d​Ui(x)j(x)​(g)dμ(g),i,j∈Fun(d,N).

The tensor power H⊗dH^{\otimes d}H⊗d has the orthonormal basis ei=ei(1)⊗⋯⊗ei(d)e_i=e_{i(1)}\otimes\cdots\otimes e_{i(d)}ei​=ei(1)​⊗⋯⊗ei(d)​, and the tensor-power matrix U⊗d(g)U^{\otimes d}(g)U⊗d(g) has entries Uij⊗d(g)=∏x=1dUi(x)j(x)(g)U^{\otimes d}_{ij}(g)=\prod_{x=1}^dU_{i(x)j(x)}(g)Uij⊗d​(g)=∏x=1d​Ui(x)j(x)​(g). The invariant tensors are

(H⊗d)G={t∈H⊗d:U⊗d(g) t=t for all g∈G},(H^{\otimes d})^G=\{t\in H^{\otimes d}: U^{\otimes d}(g)\,t=t\ \text{for all }g\in G\},(H⊗d)G={t∈H⊗d:U⊗d(g)t=t for all g∈G},

and the Haar average is the matrix P=∫GU⊗d(g) dμ(g)P=\int_GU^{\otimes d}(g)\,d\mu(g)P=∫G​U⊗d(g)dμ(g). Given vectors a1,…,am∈H⊗da_1,\dots,a_m\in H^{\otimes d}a1​,…,am​∈H⊗d, let A=[⟨ei,ax⟩]A=[\langle e_i,a_x\rangle]A=[⟨ei​,ax​⟩] be the Nd×mN^d\times mNd×m matrix whose columns are their coordinates, A∗A^*A∗ its conjugate transpose, A∗A=[⟨ax,ay⟩]A^*A=[\langle a_x,a_y\rangle]A∗A=[⟨ax​,ay​⟩] their Gram matrix, and, when it is invertible, W=(A∗A)−1W=(A^*A)^{-1}W=(A∗A)−1 the Weingarten matrix of a1,…,ama_1,\dots,a_ma1​,…,am​.

Formalization targets

Goal: Theorem 2.1 (p. 4)

If a1,…,ama_1,\dots,a_ma1​,…,am​ is a basis of (H⊗d)G(H^{\otimes d})^G(H⊗d)G, then for all i,j∈Fun(d,N)i,j\in\mathrm{Fun}(d,N)i,j∈Fun(d,N)

Iij=∑x,y=1mAix Wxy Ayj∗.I_{ij}=\sum_{x,y=1}^mA_{ix}\,W_{xy}\,A^*_{yj}.Iij​=x,y=1∑m​Aix​Wxy​Ayj∗​.

The identity is exact and holds for every compact group, every continuous unitary representation, every degree ddd and every choice of basis; nothing is fixed beyond the data.

Milestones (Section 2, pp. 3–4)

  1. Iij=PijI_{ij}=P_{ij}Iij​=Pij​: the Weingarten integrals are the matrix elements of the Haar average.
  2. P∗=PP^*=PP∗=P, P2=PP^2=PP2=P, and im⁡P=(H⊗d)G\operatorname{im}P=(H^{\otimes d})^GimP=(H⊗d)G: PPP is the orthogonal projection onto the invariant tensors (the existing platform theorem that a Haar average of a compact group representation is an idempotent fixing exactly the invariant vectors is included as a supporting reference item).
  3. For linearly independent columns, the selfadjoint idempotent with image span⁡{ax}\operatorname{span}\{a_x\}span{ax​} is A(A∗A)−1A∗A(A^*A)^{-1}A^*A(A∗A)−1A∗.

Significance

Theorem 2.1 reduces a family of analytic problems — Haar integrals of polynomials in matrix entries — to two algebraic ones: find a basis of the invariant tensors of GGG in H⊗dH^{\otimes d}H⊗d (classical invariant theory), and invert its Gram matrix. For U(N)U(N)U(N) the invariants come from Schur–Weyl duality and the Weingarten matrix becomes the Weingarten function on permutations; for O(N)O(N)O(N) and Sp(N)Sp(N)Sp(N) the invariants are indexed by pair partitions; for the symmetric group S(N)S(N)S(N) acting by permutation matrices the invariants are indicators of set partitions. Each of these case-specific formulas, surveyed in Sections 3–5 of the paper, is an instance of Theorem 2.1, so a formal proof of the theorem is the common entry point for formalizing any of them.

The theorem is classical and its proof is elementary once stated; it is not open. No machine-checked proof is known to exist in Mathlib or on this platform: Mathlib has Haar measure, unitary groups and Gram matrices (including the equivalence of linear independence with invertibility of the Gram matrix), but not the averaging projection of a representation in matrix form nor the formula A(A∗A)−1A∗A(A^*A)^{-1}A^*A(A∗A)−1A∗ for orthogonal projections. The mission produces both as reusable lemmas and the Fundamental Theorem itself.

Difficulty

The argument is short on paper; the formal work sits in three places. First, selfadjointness of PPP needs ∫GF(g−1) dμ=∫GF(g) dμ\int_GF(g^{-1})\,d\mu=\int_GF(g)\,d\mu∫G​F(g−1)dμ=∫G​F(g)dμ, i.e. unimodularity of a compact group, which Mathlib does not provide as an instance for non-commutative groups. Second, identifying the image of PPP with the invariant tensors needs both directions: invariant vectors are fixed by PPP, and every vector PtPtPt is invariant, which uses left invariance of μ\muμ inside a matrix-valued Bochner integral. Third, the projection formula needs invertibility of A∗AA^*AA∗A from linear independence and an argument that a selfadjoint idempotent is determined by its image. Commuting the integral with matrix entries and products is routine but must be done with care for the matrix norm chosen.

Formalization scope

  • HHH is CN\mathbb C^NCN with its standard basis (a finite-dimensional Hilbert space with a specified orthonormal basis is the same thing after identification); UUU is a monoid homomorphism G →* Matrix.unitaryGroup (Fin N) ℂ, assumed continuous. Indices are 0-based: [N][N][N] is Fin N, multi-indices are Fin d → Fin N, tensors are (Fin d → Fin N) → ℂ.
  • GGG carries [IsTopologicalGroup G] [CompactSpace G] [T2Space G] [BorelSpace G], and μ\muμ is a Haar measure that is a probability measure. Right invariance is not assumed; it is a consequence that solvers must prove.
  • ddd ranges over all of N\mathbb NN, including d=0d=0d=0, where both sides equal 111; the paper takes ddd positive.
  • The inner product is conjugate-linear in the first slot, so Ayj∗=Ajy‾A^*_{yj}=\overline{A_{jy}}Ayj∗​=Ajy​​ and A∗A^*A∗ is Mathlib's conjugate transpose.
  • PPP is a Bochner integral of matrices, taken with the entrywise sup norm; in finite dimension the integral is independent of the norm.
  • The basis hypothesis is two hypotheses: LinearIndependent ℂ a and span (range a) = invariants U d, with the invariant subspace defined from the representation. Without independence Mathlib's matrix inverse returns 000 on a singular Gram matrix, and the statement would be false; it is therefore not a trivializing gap but a necessary hypothesis.
  • Correction. The paper says "a basis a1,…,ama_1,\dots,a_ma1​,…,am​ of (H)⊗d(H)^{\otimes d}(H)⊗d"; the same paragraph calls the axa_xax​ "the basic invariants", and with a basis of all of H⊗dH^{\otimes d}H⊗d the right-hand side would be the identity. The formal statement uses a basis of (H⊗d)G(H^{\otimes d})^G(H⊗d)G.
  • A trivializing formalization — defining the Weingarten integral as AWA∗AWA^*AWA∗, defining the invariant space as the image of PPP, or assuming that PPP projects onto the span of the axa_xax​ — is excluded: the left side is the Haar integral of the product of matrix elements written out, and the invariant space is defined by U⊗d(g)t=tU^{\otimes d}(g)t=tU⊗d(g)t=t.

Needed infrastructure: matrix-valued Bochner integrals commuting with entries and with matrix multiplication; unimodularity of compact groups; the orthogonal-projection formula. All three are reusable well beyond this mission, and contributions of these general lemmas are welcome.

Selected references

  • B. Collins, S. Matsumoto, J. Novak, The Weingarten Calculus, arXiv preprint, 2021. https://arxiv.org/abs/2109.14890
  • D. Weingarten, Asymptotic behavior of group integrals in the limit of infinite rank, J. Math. Phys. 19 (1978), 999–1001. https://doi.org/10.1063/1.523807
  • B. Collins, Moments and cumulants of polynomial random variables on unitary groups, the Itzykson–Zuber integral, and free probability, Int. Math. Res. Not. 2003, no. 17, 953–982. https://doi.org/10.1155/S107379280320917X
  • B. Collins, P. Śniady, Integration with respect to the Haar measure on unitary, orthogonal and symplectic group, Comm. Math. Phys. 264 (2006), 773–795. https://doi.org/10.1007/s00220-006-1554-3
6 thms0 active usersReviewed
🏆Completed
Group TheoryMachine Learning·Captain: mikedeng1

Geometric Deep Learning I: A Linear G-Invariant Function of a Signal Depends on It Only Through the Group Average Ax = (1/μ(G)) ∫ g.x dμ(g)Textbook

Symmetry as a design principle

Geometric deep learning is the programme, set out by Bronstein, Bruna, Cohen and Veličković in Geometric Deep Learning: Grids, Groups, Graphs, Geodesics, and Gauges (arXiv:2104.13478v2, 2021), of deriving neural network architectures from the symmetries of the domain on which the data live. Convolutional networks on images, permutation-invariant networks on sets, graph neural networks and equivariant networks for molecules are presented there as instances of one template. Chapter 3 of the book, Geometric Priors, gives that template: it fixes what a symmetry of a signal is, what it means for a function to respect it, and how invariant functions are assembled from simpler blocks. This mission formalizes the definitions and the elementary claims of Chapter 3. Its goal is the observation the book uses to explain why such networks must be non-linear.

The book numbers no theorem, lemma or definition. Its only labels are the sections, the equations (1)–(53) and the figures. Every item of this mission is therefore cited by section, equation number (where there is one) and printed page of the arXiv v2 text; a margin note is cited as "side-note, p. n".

Setting

A domain is a set Ω\OmegaΩ. A signal on Ω\OmegaΩ with values in a set of channels CCC is a function x:Ω→Cx : \Omega \to Cx:Ω→C; the signals form the space X(Ω,C)X(\Omega, C)X(Ω,C) (equation (1)), written X(Ω)X(\Omega)X(Ω) when C=RC = \mathbb RC=R.

A group G\mathfrak GG acts on Ω\OmegaΩ through a map (g,u)↦g.u(\mathfrak g, u) \mapsto \mathfrak g . u(g,u)↦g.u with g.(h.u)=(gh).u\mathfrak g . (\mathfrak h . u) = (\mathfrak g \mathfrak h) . ug.(h.u)=(gh).u. It then acts on signals by equation (3):

(g.x)(u)=x(g−1u).(\mathfrak g . x)(u) = x(\mathfrak g^{-1} u).(g.x)(u)=x(g−1u).

On real signals this action is linear, so ρ(g)x=g.x\rho(\mathfrak g) x = \mathfrak g . xρ(g)x=g.x defines a representation ρ\rhoρ of G\mathfrak GG on X(Ω)X(\Omega)X(Ω).

A function f:X(Ω,C)→Yf : X(\Omega, C) \to Yf:X(Ω,C)→Y is G\mathfrak GG-invariant if f(g.x)=f(x)f(\mathfrak g . x) = f(x)f(g.x)=f(x) for all g\mathfrak gg and xxx. A map B:X(Ω,C)→X(Ω′,C′)B : X(\Omega, C) \to X(\Omega', C')B:X(Ω,C)→X(Ω′,C′) between signal spaces on which G\mathfrak GG acts is G\mathfrak GG-equivariant if B(g.x)=g.B(x)B(\mathfrak g . x) = \mathfrak g . B(x)B(g.x)=g.B(x). A map σ:C→C′\sigma : C \to C'σ:C→C′ acts on signals element-wise, (σ(x))(u)=σ(x(u))(\boldsymbol\sigma(x))(u) = \sigma(x(u))(σ(x))(u)=σ(x(u)).

For a measure μ\muμ on G\mathfrak GG with 0<μ(G)<∞0 < \mu(\mathfrak G) < \infty0<μ(G)<∞, the G\mathfrak GG-average of a real signal is

Ax=1μ(G)∫G(g.x) dμ(g).A x = \frac{1}{\mu(\mathfrak G)} \int_{\mathfrak G} (\mathfrak g . x)\, \mathrm d\mu(\mathfrak g).Ax=μ(G)1​∫G​(g.x)dμ(g).

The book takes μ\muμ to be the Haar measure of the group.

The metric distortion of a map η\etaη between metric spaces is dis⁡(η)=sup⁡u,v∣d(η(u),η(v))−d(u,v)∣\operatorname{dis}(\eta) = \sup_{u,v} |d(\eta(u), \eta(v)) - d(u, v)|dis(η)=supu,v​∣d(η(u),η(v))−d(u,v)∣ (§4.4, p. 51). A self-map τ\tauτ with distortion at most ε\varepsilonε is an ε\varepsilonε-isometry.

The Geometric Deep Learning Blueprint (§3.5, p. 29) builds functions

f=A∘σJ∘BJ∘PJ−1∘⋯∘P1∘σ1∘B1f = A \circ \boldsymbol\sigma_J \circ B_J \circ P_{J-1} \circ \cdots \circ P_1 \circ \boldsymbol\sigma_1 \circ B_1f=A∘σJ​∘BJ​∘PJ−1​∘⋯∘P1​∘σ1​∘B1​

from equivariant linear layers BjB_jBj​, element-wise nonlinearities σj\boldsymbol\sigma_jσj​, poolings PjP_jPj​ onto coarser domains, and an invariant global pooling AAA.

Formalization targets

Goal: linear invariants see only the average (§3.5, p. 27)

For a finite domain Ω\OmegaΩ, a finite nonzero measure μ\muμ on G\mathfrak GG, and a linear G\mathfrak GG-invariant map fff from X(Ω)X(\Omega)X(Ω) into a real Banach space, for every signal xxx,

f(x)=1μ(G)∫Gf(g.x) dμ(g)=f ⁣(1μ(G)∫G(g.x) dμ(g))=f(Ax).f(x) = \frac{1}{\mu(\mathfrak G)} \int_{\mathfrak G} f(\mathfrak g . x)\, \mathrm d\mu(\mathfrak g) = f\!\left( \frac{1}{\mu(\mathfrak G)} \int_{\mathfrak G} (\mathfrak g . x)\, \mathrm d\mu(\mathfrak g) \right) = f(Ax).f(x)=μ(G)1​∫G​f(g.x)dμ(g)=f(μ(G)1​∫G​(g.x)dμ(g))=f(Ax).

A linear invariant therefore depends on xxx only through AxAxAx.

Milestones, in attack order

  1. The action (3) is a group action: g.(h.x)=(gh).x\mathfrak g . (\mathfrak h . x) = (\mathfrak g\mathfrak h) . xg.(h.x)=(gh).x (p. 14).
  2. It is linear: g.(αx+βx′)=α(g.x)+β(g.x′)\mathfrak g . (\alpha x + \beta x') = \alpha(\mathfrak g . x) + \beta(\mathfrak g . x')g.(αx+βx′)=α(g.x)+β(g.x′) (p. 15).
  3. It is a representation: ρ(g)(ρ(h)x)=ρ(gh)x\rho(\mathfrak g)(\rho(\mathfrak h)x) = \rho(\mathfrak g\mathfrak h)xρ(g)(ρ(h)x)=ρ(gh)x (p. 15).
  4. Two ε\varepsilonε-isometries compose to a 2ε2\varepsilon2ε-isometry, and ε\varepsilonε-isometries are not closed under composition (side-note, p. 20).
  5. U=σ∘BU = \boldsymbol\sigma \circ BU=σ∘B is equivariant for equivariant BBB and any element-wise σ\sigmaσ (p. 27).
  6. The G\mathfrak GG-average is invariant, so A∘UA \circ UA∘U is invariant (p. 27).
  7. The blueprint fff is G\mathfrak GG-invariant (p. 29).

Companion: the Fourier modulus (§3.4, p. 25)

With x^(ξ)=∫x(u)e−iξu du\hat x(\xi) = \int x(u) e^{-\mathrm i\xi u}\,\mathrm dux^(ξ)=∫x(u)e−iξudu and (Svx)(u)=x(u−v)(S_v x)(u) = x(u - v)(Sv​x)(u)=x(u−v), Svx^(ξ)=e−iξvx^(ξ)\widehat{S_v x}(\xi) = e^{-\mathrm i \xi v}\hat x(\xi)Sv​x​(ξ)=e−iξvx^(ξ), so ∣Svx^∣=∣x^∣|\widehat{S_v x}| = |\hat x|∣Sv​x​∣=∣x^∣.

Significance

The goal is the book's reason for non-linearity. Linear invariant models collapse to functions of the group average, so they cannot tell apart two images with the same mean colour. Expressive invariant models must interleave equivariant linear layers with non-linear maps. Milestones 5–7 are the positive half of the same argument: element-wise nonlinearities keep equivariance, the average turns equivariance into invariance, and stacking such blocks gives the blueprint that the later chapters instantiate on grids, groups, graphs and manifolds. Milestone 4 is the book's reason for leaving exact symmetry behind: small deformations do not form a group, so it replaces invariance to deformations by stability to them.

Each claim of the chapter is proved in one or two lines on the page, and none is new mathematics. What a formalization adds is the precise setting in which they hold. The book names the Haar measure but its identity needs no invariance of μ\muμ; the average's invariance needs right invariance, which the page does not state; and the printed blueprint is invariant only if the poolings are equivariant. As far as a search of the platform and Mathlib shows, none of these statements has a machine-checked proof. The definitions here (the action on signals, invariance, equivariance, the average, the blueprint) are reusable by every later mission of this series.

Difficulty

The obvious argument for the goal pulls fff through the integral, which uses linearity in the form "a continuous linear map commutes with the Bochner integral". That needs continuity of fff and integrability of g↦g.x\mathfrak g \mapsto \mathfrak g . xg↦g.x, and neither is stated on the page. On a finite domain both follow from the setting, finite dimension for continuity and boundedness with measurability for integrability. The measurability of g↦g−1.u\mathfrak g \mapsto \mathfrak g^{-1} . ug↦g−1.u has to be assumed of the group and the action.

For the average's invariance, the substitution g↦gh\mathfrak g \mapsto \mathfrak g\mathfrak hg↦gh is valid only for a right-invariant measure. For the blueprint, the composite's invariance fails for a non-equivariant pooling, so the box's claim is false as printed. For the ε\varepsilonε-isometries, the subtle half is the negative one: exhibiting a deformation that is an ε\varepsilonε-isometry while its square is not.

Formalization scope

  • Signals are functions Ω → C. The goal and the average use real signals on a finite domain ([Fintype Ω]), following the book's side-note on p. 15 that one must always discretise to a finite grid in practice; the other milestones allow any domain and any channel types.
  • The action (3) is a function act g x := fun u => x (g⁻¹ • u) for [MulAction G Ω], not a MulAction instance on Ω → C.
  • The measure is any μ with [IsFiniteMeasure μ] and μ ≠ 0. Without those hypotheses Lean's (μ univ).toReal⁻¹ would be 000 and the goal false. The goal assumes no invariance of μ. Milestone 6 adds [μ.IsMulRightInvariant], which the Haar measure of a compact group satisfies.
  • Measurability: [MeasurableInv G], [MeasurableSMul G Ω] and the discrete σ-algebra on Ω; integrability of g ↦ g.x is a consequence, not a hypothesis. A finite group with the counting measure and a compact group with its Haar measure both fit.
  • The blueprint is indexed from 000 to n=J−1n = J - 1n=J−1. A single group acts on all domains, and the poolings are assumed equivariant (a correction of the printed box, disclosed in the statement). The box's remark that different blocks may use different groups is not formalized.
  • ε\varepsilonε-isometries are defined pointwise (∣d(τu,τv)−d(u,v)∣≤ε|d(\tau u, \tau v) - d(u,v)| \le \varepsilon∣d(τu,τv)−d(u,v)∣≤ε for all u,vu, vu,v), avoiding a real supremum that would be 000 on an unbounded set. Non-closure is witnessed on R\mathbb RR by a continuous bijection.
  • The Fourier transform of the companion is the book's ∫x(u)e−iξu du\int x(u)e^{-\mathrm i\xi u}\,\mathrm du∫x(u)e−iξudu (no 2π2\pi2π), a Bochner integral over R\mathbb RR.

The goal cannot be satisfied trivially. groupAverage is defined from μ and the action alone, the goal quantifies over every linear invariant f, and the measure hypotheses exclude the degenerate case in which the average is 000. Proofs, finite-sum encodings of the average, and versions for complex or vector-valued channels are welcome.

Selected references

  • M. M. Bronstein, J. Bruna, T. Cohen, P. Veličković, Geometric Deep Learning: Grids, Groups, Graphs, Geodesics, and Gauges, arXiv preprint, 2021. https://arxiv.org/abs/2104.13478v2
  • S. Mallat, Group Invariant Scattering, Communications on Pure and Applied Mathematics 65(10), 2012. https://doi.org/10.1002/cpa.21413
10 thms0 active usersReviewed
Mathematical PhysicsQuantum Information·Captain: mikedeng1

SU(2)-Irreducibly Covariant and EPOSIC Channels: EPOSIC Channels Are Exactly the Extreme Points of the Convex Set of SU(2)-Irreducibly Covariant ChannelsResearch Paper

Motivation

A quantum channel is a completely positive, trace-preserving linear map between matrix algebras; it models the most general physical transformation of a quantum system. When a group GGG acts on the input and output systems, a channel is GGG-covariant if it commutes with the group action. Covariant channels appear whenever a physical process respects a symmetry: rotation-invariant noise on spin systems, optimal cloning and estimation of spin states, and quantum reference frames for orientation. Because the covariant channels between two fixed representations form a convex set, the first structural question is what its extreme points are; they are the channels from which every covariant channel is obtained by mixing.

For G=SU(2)G = SU(2)G=SU(2) acting irreducibly on input and output, the tensor product of the output representation with the dual of the input one is multiplicity free, so the commutant that controls covariant maps is abelian (Vollbrecht and Werner, 2001, arXiv:quant-ph/0010095). M. Al Nuwairan (arXiv:1306.5321, 2013) makes this concrete. She constructs explicit channels Φm,n,h\Phi_{m,n,h}Φm,n,h​, called EPOSIC channels, from Clebsch–Gordan isometries realized on polynomials. She computes their Choi matrices and shows that they are exactly the extreme points of the convex set of SU(2)-irreducibly covariant channels.

Setting

For k∈Nk \in \mathbb Nk∈N, let PkP_kPk​ be the space of homogeneous polynomials of degree kkk in two variables x1,x2x_1, x_2x1​,x2​. The group

SU(2)={(ab−bˉaˉ):a,b∈C, ∣a∣2+∣b∣2=1}SU(2)=\left\{\begin{pmatrix}a&b\\-\bar b&\bar a\end{pmatrix}: a,b\in\mathbb C,\ |a|^2+|b|^2=1\right\}SU(2)={(a−bˉ​baˉ​):a,b∈C, ∣a∣2+∣b∣2=1}

acts on PkP_kPk​ by the representation ρk\rho_kρk​,

(ρk(g)f)(x1,x2)=f((x1,x2)g)=f(ax1−bˉx2, bx1+aˉx2).(\rho_k(g)f)(x_1,x_2)=f\big((x_1,x_2)g\big)=f(ax_1-\bar bx_2,\ bx_1+\bar ax_2).(ρk​(g)f)(x1​,x2​)=f((x1​,x2​)g)=f(ax1​−bˉx2​, bx1​+aˉx2​).

With the inner product ⟨x1lx2k−l,x1jx2k−j⟩=l!(k−l)! δlj\langle x_1^lx_2^{k-l},x_1^jx_2^{k-j}\rangle=l!(k-l)!\,\delta_{lj}⟨x1l​x2k−l​,x1j​x2k−j​⟩=l!(k−l)!δlj​, the canonical basis flk=x1lx2k−l/l!(k−l)!f^k_l=x_1^lx_2^{k-l}/\sqrt{l!(k-l)!}flk​=x1l​x2k−l​/l!(k−l)!​, 0≤l≤k0\le l\le k0≤l≤k, is orthonormal and ρk\rho_kρk​ is unitary. The ρk\rho_kρk​ are the irreducible representations of SU(2)SU(2)SU(2).

The tensor product Pm⊗PnP_m\otimes P_nPm​⊗Pn​ is realized as the polynomials in x=(x1,x2)x=(x_1,x_2)x=(x1​,x2​) and y=(y1,y2)y=(y_1,y_2)y=(y1​,y2​) of bi-degree (m,n)(m,n)(m,n). On these polynomials act the operators

Δxy=x1∂y1+x2∂y2,Δyx=y1∂x1+y2∂x2,Γxy=(x1y2−y1x2)⋅,Ωxy=∂x1∂y2−∂x2∂y1.\Delta_{xy}=x_1\partial_{y_1}+x_2\partial_{y_2},\quad \Delta_{yx}=y_1\partial_{x_1}+y_2\partial_{x_2},\quad \Gamma_{xy}=(x_1y_2-y_1x_2)\cdot,\quad \Omega_{xy}=\partial_{x_1}\partial_{y_2}-\partial_{x_2}\partial_{y_1}.Δxy​=x1​∂y1​​+x2​∂y2​​,Δyx​=y1​∂x1​​+y2​∂x2​​,Γxy​=(x1​y2​−y1​x2​)⋅,Ωxy​=∂x1​​∂y2​​−∂x2​​∂y1​​.

Real coefficients cm,n,hc_{m,n,h}cm,n,h​ are defined by the recursion cm,0,0=1c_{m,0,0}=1cm,0,0​=1 and, for n≥1n\ge 1n≥1, cm,n,h=1(m+1)n cm+1,n−1,hc_{m,n,h}=\frac{1}{(m+1)n}\,c_{m+1,n-1,h}cm,n,h​=(m+1)n1​cm+1,n−1,h​ if h=0h=0h=0, 1(m+1)n (cm−1,n−1,h−1+cm+1,n−1,h)\frac{1}{(m+1)n}\,(c_{m-1,n-1,h-1}+c_{m+1,n-1,h})(m+1)n1​(cm−1,n−1,h−1​+cm+1,n−1,h​) if 0<h<n0<h<n0<h<n, and 1(m+1)n cm−1,n−1,h−1\frac{1}{(m+1)n}\,c_{m-1,n-1,h-1}(m+1)n1​cm−1,n−1,h−1​ if h=nh=nh=n. For 0≤h≤min⁡{m,n}0\le h\le\min\{m,n\}0≤h≤min{m,n} and r=m+n−2hr=m+n-2hr=m+n−2h, the map

αm,n,h:Pr→Pm⊗Pn,αm,n,h(f)=cm,n,h  ΓxyhΔyxn−hf\alpha_{m,n,h}:P_r\to P_m\otimes P_n,\qquad \alpha_{m,n,h}(f)=\sqrt{c_{m,n,h}}\;\Gamma_{xy}^h\Delta_{yx}^{n-h}fαm,n,h​:Pr​→Pm​⊗Pn​,αm,n,h​(f)=cm,n,h​​Γxyh​Δyxn−h​f

is an explicit Clebsch–Gordan isometry. The EPOSIC channel is

Φm,n,h(A)=TrPn(αm,n,hAαm,n,h∗),A∈End(Pr).\Phi_{m,n,h}(A)=\mathrm{Tr}_{P_n}\big(\alpha_{m,n,h}A\alpha_{m,n,h}^*\big),\qquad A\in\mathrm{End}(P_r).Φm,n,h​(A)=TrPn​​(αm,n,h​Aαm,n,h∗​),A∈End(Pr​).

For r,m∈Nr,m\in\mathbb Nr,m∈N, EC(r,m)EC(r,m)EC(r,m) is the set of the channels Φm,n,h\Phi_{m,n,h}Φm,n,h​ with r=m+n−2hr=m+n-2hr=m+n−2h and 0≤h≤min⁡{m,n}0\le h\le\min\{m,n\}0≤h≤min{m,n}. The set QCSU(2)(ρr,ρm)QC_{SU(2)}(\rho_r,\rho_m)QCSU(2)​(ρr​,ρm​) consists of the quantum channels Φ:End(Pr)→End(Pm)\Phi:\mathrm{End}(P_r)\to\mathrm{End}(P_m)Φ:End(Pr​)→End(Pm​) with Φ(ρr(g)Aρr(g)∗)=ρm(g)Φ(A)ρm(g)∗\Phi(\rho_r(g)A\rho_r(g)^*)=\rho_m(g)\Phi(A)\rho_m(g)^*Φ(ρr​(g)Aρr​(g)∗)=ρm​(g)Φ(A)ρm​(g)∗ for all AAA and all g∈SU(2)g\in SU(2)g∈SU(2).

Formalization targets

Goal: Proposition 4.1

ext⁡QCSU(2)(ρr,ρm)=EC(r,m)for all r,m∈N.\operatorname{ext} QC_{SU(2)}(\rho_r,\rho_m)=EC(r,m)\qquad\text{for all } r,m\in\mathbb N.extQCSU(2)​(ρr​,ρm​)=EC(r,m)for all r,m∈N.

The statement has no hypotheses and covers the degenerate cases: for min⁡{r,m}=0\min\{r,m\}=0min{r,m}=0 the set EC(r,m)EC(r,m)EC(r,m) is a single channel.

Milestones

In the order of the paper's argument:

  1. the polynomial decomposition identity behind cm,n,hc_{m,n,h}cm,n,h​ (Theorem 2.7), the adjoint of αm,n,h\alpha_{m,n,h}αm,n,h​ and the resolution of the identity ∑hαm,n,hαm,n,h∗=I\sum_h\alpha_{m,n,h}\alpha_{m,n,h}^*=I∑h​αm,n,h​αm,n,h∗​=I (Lemma 2.9);
  2. αm,n,h\alpha_{m,n,h}αm,n,h​ is an SU(2)-equivariant isometry (Proposition 2.10), so Φm,n,h∈QCSU(2)(ρr,ρm)\Phi_{m,n,h}\in QC_{SU(2)}(\rho_r,\rho_m)Φm,n,h​∈QCSU(2)​(ρr​,ρm​) (Proposition 2.18), and EC(r,m)={Φm,r+m−2l,m−l:0≤l≤min⁡{r,m}}EC(r,m)=\{\Phi_{m,r+m-2l,m-l}:0\le l\le\min\{r,m\}\}EC(r,m)={Φm,r+m−2l,m−l​:0≤l≤min{r,m}} (Proposition 2.21);
  3. the orthogonal decomposition of Pm⊗PˉnP_m\otimes\bar P_nPm​⊗Pˉn​ into irreducible subspaces with projections qqq (Corollary 2.13), covariance in terms of the Choi matrix (Proposition 1.15(1)), and the commutant of ρm⊗ρˇr\rho_m\otimes\check\rho_rρm​⊗ρˇ​r​ (Proposition 3.9);
  4. the Choi matrix C(Φm,n,h)=r+1n+1qm,r,m−hC(\Phi_{m,n,h})=\frac{r+1}{n+1}q_{m,r,m-h}C(Φm,n,h​)=n+1r+1​qm,r,m−h​ (Lemma 3.12, Proposition 3.13);
  5. EC(r,m)EC(r,m)EC(r,m) spans the SU(2)-equivariant linear maps (Corollary 3.16), and QCSU(2)(ρr,ρm)=conv⁡EC(r,m)QC_{SU(2)}(\rho_r,\rho_m)=\operatorname{conv}EC(r,m)QCSU(2)​(ρr​,ρm​)=convEC(r,m) (Corollary 3.17).

Significance

The result shows that the covariant channels between two irreducible SU(2) representations form a simplex. Its min⁡{r,m}+1\min\{r,m\}+1min{r,m}+1 vertices are given by closed formulas, so optimization problems over SU(2)-covariant channels reduce to finite-dimensional linear programs over the vertex weights: fidelities, output entropies, capacities restricted to covariant encodings. The explicit Choi matrices make the vertices computable in any basis. The same channels are used later in the paper to compute complementary channels and to exhibit a positive map that is not completely positive.

On the formal side, the mission produces the irreducible representations of SU(2)SU(2)SU(2) on homogeneous polynomials, the polarization operators, and explicit Clebsch–Gordan isometries with their equivariance and isometry proofs. It also produces a Choi-matrix criterion for covariance and the extreme-point characterization. Theorem 2.7, Proposition 1.15 and Proposition 3.9 are cited in the paper from Procesi, from [4] and from Vollbrecht–Werner, and are proved neither there nor, as far as is known, in any proof assistant. Nothing in this mission is machine-checked yet.

Difficulty

That Φm,n,h\Phi_{m,n,h}Φm,n,h​ is a covariant channel follows from equivariance of αm,n,h\alpha_{m,n,h}αm,n,h​. The real content lies elsewhere.

First, αm,n,h\alpha_{m,n,h}αm,n,h​ is an isometry only because of the specific constants cm,n,hc_{m,n,h}cm,n,h​. The paper derives this from the polynomial identity of Theorem 2.7 and Schur's lemma, and neither is available off the shelf for these polynomial models.

Second, the extreme-point statement needs the full commutant of ρm⊗ρˉr\rho_m\otimes\bar\rho_rρm​⊗ρˉ​r​: the projections onto the irreducible summands must span it. Exhibiting the projections is not enough; one must show there is nothing else, which is multiplicity-freeness of the tensor product. This requires the irreducibility of each ρk\rho_kρk​ and the Clebsch–Gordan dimension count.

Third, matching each EPOSIC channel with a single projection (Proposition 3.13) requires a rank argument. The obvious approach, expanding the Choi matrix in coordinates, gives sums of products of signed binomial coefficients whose non-vanishing is the content of Lemma 3.12.

Formalization scope

Every space PkP_kPk​ is represented by coordinates in its canonical orthonormal basis, indexed by Fin (k+1) with 0-based indices. So End(Pk)\mathrm{End}(P_k)End(Pk​) is the algebra of (k+1)×(k+1)(k+1)\times(k+1)(k+1)×(k+1) complex matrices, and adjoints for the inner product above are conjugate transposes. The representation ρk\rho_kρk​ and the maps αm,n,h\alpha_{m,n,h}αm,n,h​ are computed from the polynomials themselves in MvPolynomial (Fin 4) ℂ (variables x1,x2,y1,y2x_1,x_2,y_1,y_2x1​,x2​,y1​,y2​, derivatives pderiv), then read off in the canonical bases. The action is the row-vector action f↦f((x1,x2)g)f\mapsto f((x_1,x_2)g)f↦f((x1​,x2​)g) of the paper. Pm⊗PnP_m\otimes P_nPm​⊗Pn​ is indexed with the PmP_mPm​ factor first. The conjugate space Pˉk\bar P_kPˉk​ has the same coordinates and carries the entrywise conjugate matrices. The coefficients cm,n,hc_{m,n,h}cm,n,h​ are defined by the recursion, not by any closed formula. The input degree rrr of α\alphaα is an explicit argument, constrained by r+2h=m+nr+2h=m+nr+2h=m+n in every statement.

Channels are maps on matrices given by a finite normalized Kraus family (the published QAlgorithms.IsChannel). The partial trace is the published WildeQIT.partialTraceRight, over the PnP_nPn​ factor. The Choi matrix is ∑ijΦ(Eij)⊗Eij\sum_{ij}\Phi(E_{ij})\otimes E_{ij}∑ij​Φ(Eij​)⊗Eij​ with the output factor first. Convex hulls and extreme points are taken over R\mathbb RR in the space of maps.

No statement replaces ρk\rho_kρk​ by an abstract irreducible representation. None defines EC(r,m)EC(r,m)EC(r,m) through extreme points, Choi matrices or projections, defines αm,n,h\alpha_{m,n,h}αm,n,h​ by its coefficient expansion, or assumes Schur's lemma, the commutant description or Theorem 2.7 as a hypothesis. Each of these would make the goal or a milestone true by definition.

A complete development needs:

  • the irreducibility and unitarity of ρk\rho_kρk​;
  • Schur's lemma for these matrix representations;
  • the Clebsch–Gordan decomposition and the polynomial identities of Theorem 2.7;
  • Kraus/Stinespring/Choi correspondences for finite-dimensional channels.

The SU(2) representation layer and the Choi-matrix covariance criterion are reusable beyond this mission. Proofs of the cited results (Theorem 2.7, Propositions 1.15 and 3.9) are welcome as independent contributions.

Selected references

  • M. Al Nuwairan, SU(2)-irreducibly covariant and EPOSIC channels, arXiv preprint, 2013. https://arxiv.org/abs/1306.5321
  • K. G. H. Vollbrecht and R. F. Werner, Entanglement measures under symmetry, Phys. Rev. A 64, 062307, 2001. https://arxiv.org/abs/quant-ph/0010095
  • C. Procesi, Lie Groups: An Approach through Invariants and Representations, Springer Universitext, 2007. https://doi.org/10.1007/978-0-387-28929-8
  • J. Watrous, The Theory of Quantum Information, Cambridge University Press, 2018. https://doi.org/10.1017/9781316848142
17 thms0 active usersReviewed
Linear algebraQuantum Information·Captain: mikedeng1

Commutant Structure of U⊗(n−1) ⊗ U* Transformations: The Operator Basis v^{ab}_{ij}(α) Is a Linear Combination of Partially Transposed Permutation Operators, and ConverselyResearch Paper

Motivation

Operators on nnn copies of Cd\mathbb C^dCd that commute with every U⊗nU^{\otimes n}U⊗n, UUU unitary, are exactly the linear combinations of permutation operators V(σ)V(\sigma)V(σ), σ∈S(n)\sigma\in S(n)σ∈S(n) (Schur–Weyl duality). Many protocols in quantum information carry a different symmetry: the first n−1n-1n−1 systems transform by UUU and the last by the complex conjugate Uˉ\bar UUˉ. It is the symmetry of port-based teleportation (Ishizaka–Hiroshima 2008), whose optimal performance is computed from the representation theory of exactly this commutant (Studziński et al. 2017). The commutant of U⊗(n−1)⊗UˉU^{\otimes(n-1)}\otimes\bar UU⊗(n−1)⊗Uˉ is the algebra A\mathcal AA spanned by the partially transposed permutation operators V′(σ)V'(\sigma)V′(σ): the permutation operators with the transpose taken on the last system.

Studziński, Horodecki and Mozrzymas give an explicit operator basis of an ideal of A\mathcal AA adapted to the irreducible representations of S(n−2)S(n-2)S(n−2). Their main structural result, Theorem 9, states that this basis and the partially transposed permutation operators are linear combinations of each other, with coefficients that are irrep matrix elements.

Setting

Let H=Cd\mathcal H=\mathbb C^dH=Cd and n≥2n\ge2n≥2. For σ∈S(n)\sigma\in S(n)σ∈S(n), the permutation operator V(σ)V(\sigma)V(σ) on H⊗n\mathcal H^{\otimes n}H⊗n is V(σ)∣i1…in⟩=∣iσ−1(1)…iσ−1(n)⟩V(\sigma)|i_1\dots i_n\rangle=|i_{\sigma^{-1}(1)}\dots i_{\sigma^{-1}(n)}\rangleV(σ)∣i1​…in​⟩=∣iσ−1(1)​…iσ−1(n)​⟩. The partial transposition X↦X′X\mapsto X'X↦X′ transposes the last system: (∣x1…xn⟩⟨y1…yn∣)′=∣x1…xn−1yn⟩⟨y1…yn−1xn∣(|x_1\dots x_n\rangle\langle y_1\dots y_n|)'=|x_1\dots x_{n-1}y_n\rangle\langle y_1\dots y_{n-1}x_n|(∣x1​…xn​⟩⟨y1​…yn​∣)′=∣x1​…xn−1​yn​⟩⟨y1​…yn−1​xn​∣. The algebra A\mathcal AA is the span of all V′(τ)V'(\tau)V′(τ), τ∈S(n)\tau\in S(n)τ∈S(n). S(k)S(k)S(k) is regarded inside S(n)S(n)S(n) as the permutations fixing the systems after kkk.

Fix π1,…,πn−1∈S(n−1)\pi_1,\dots,\pi_{n-1}\in S(n-1)π1​,…,πn−1​∈S(n−1) with πk(n−1)=k\pi_k(n-1)=kπk​(n−1)=k, and write Πk=V(πk)\Pi_k=V(\pi_k)Πk​=V(πk​). For labels a,b∈{1,…,n−1}a,b\in\{1,\dots,n-1\}a,b∈{1,…,n−1}, the maps fab:S(n−2)→S(n)f_{ab}:S(n-2)\to S(n)fab​:S(n−2)→S(n) are

fab(σ)=πb∘σ∘(n n−1)∘πa−1.f_{ab}(\sigma)=\pi_b\circ\sigma\circ(n\ n-1)\circ\pi_a^{-1}.fab​(σ)=πb​∘σ∘(n n−1)∘πa−1​.

The standing data are a family of real orthogonal, irreducible, pairwise inequivalent matrix representations φα\varphi^\alphaφα of S(n−2)S(n-2)S(n−2), of dimensions dαd_\alphadα​, and a Schur basis: an orthonormal basis ∣ϕi(α,r)⟩|\phi_i(\alpha,r)\rangle∣ϕi​(α,r)⟩ of H⊗(n−2)\mathcal H^{\otimes(n-2)}H⊗(n−2), where rrr ranges over the multiplicity of α\alphaα, such that V(σ)∣ϕj(α,r)⟩=∑iφijα(σ)∣ϕi(α,r)⟩V(\sigma)|\phi_j(\alpha,r)\rangle=\sum_i\varphi^\alpha_{ij}(\sigma)|\phi_i(\alpha,r)\rangleV(σ)∣ϕj​(α,r)⟩=∑i​φijα​(σ)∣ϕi​(α,r)⟩. With ∣Φ+⟩=1d∑l∣ll⟩|\Phi^+\rangle=\frac1{\sqrt d}\sum_l|ll\rangle∣Φ+⟩=d​1​∑l​∣ll⟩ on the last two systems, the paper defines

∣ψik(α,r)⟩=d Πk∣ϕi(α,r)⟩∣Φ+⟩,vijab(α)=∑r∣ψia(α,r)⟩⟨ψjb(α,r)∣,|\psi^k_i(\alpha,r)\rangle=\sqrt d\,\Pi_k|\phi_i(\alpha,r)\rangle|\Phi^+\rangle,\qquad v^{ab}_{ij}(\alpha)=\sum_r|\psi^a_i(\alpha,r)\rangle\langle\psi^b_j(\alpha,r)|,∣ψik​(α,r)⟩=d​Πk​∣ϕi​(α,r)⟩∣Φ+⟩,vijab​(α)=r∑​∣ψia​(α,r)⟩⟨ψjb​(α,r)∣,

and Eijα=∑r∣ϕi(α,r)⟩⟨ϕj(α,r)∣E^\alpha_{ij}=\sum_r|\phi_i(\alpha,r)\rangle\langle\phi_j(\alpha,r)|Eijα​=∑r​∣ϕi​(α,r)⟩⟨ϕj​(α,r)∣. In Lean these are IrrepFamily, SchurBasis, psi, vOp, EOp, with ptLast for ′'′, fPerm for fabf_{ab}fab​ and algA for A\mathcal AA.

Formalization targets

Goal: Theorem 9

For all a,ba,ba,b, α\alphaα, i,ji,ji,j:

vijab(α)=dα(n−2)!∑σ∈S(n−2)φijα(σ) V′(fba(σ))∈A,v^{ab}_{ij}(\alpha)=\frac{d_\alpha}{(n-2)!}\sum_{\sigma\in S(n-2)}\varphi^\alpha_{ij}(\sigma)\,V'\big(f_{ba}(\sigma)\big)\in\mathcal A,vijab​(α)=(n−2)!dα​​σ∈S(n−2)∑​φijα​(σ)V′(fba​(σ))∈A,

and for every σ∈S(n−2)\sigma\in S(n-2)σ∈S(n−2), with σab=fab(σ)\sigma_{ab}=f_{ab}(\sigma)σab​=fab​(σ),

V′(σab)=⨁α∑i,jφijα(fab−1(σab)) vijba(α).V'(\sigma_{ab})=\bigoplus_\alpha\sum_{i,j}\varphi^\alpha_{ij}\big(f_{ab}^{-1}(\sigma_{ab})\big)\,v^{ba}_{ij}(\alpha).V′(σab​)=α⨁​i,j∑​φijα​(fab−1​(σab​))vijba​(α).

The goal is the whole theorem, both directions and the membership in A\mathcal AA.

Milestones

  1. Permutations of S(n−1)S(n-1)S(n−1) are invariant under ′'′ and pass through it: V(ρ)V′(τ)V(ρ′)=V′(ρτρ′)V(\rho)V'(\tau)V(\rho')=V'(\rho\tau\rho')V(ρ)V′(τ)V(ρ′)=V′(ρτρ′) (proof of Lemma 26).
  2. Lemma 26: V′(σab)=ΠbV(fab−1(σab))V′(n n−1)Πa−1V'(\sigma_{ab})=\Pi_bV(f^{-1}_{ab}(\sigma_{ab}))V'(n\ n-1)\Pi_a^{-1}V′(σab​)=Πb​V(fab−1​(σab​))V′(n n−1)Πa−1​ with fab−1(σab)∈S(n−2)f^{-1}_{ab}(\sigma_{ab})\in S(n-2)fab−1​(σab​)∈S(n−2).
  3. Eq. (9): the Sab=fab(S(n−2))S_{ab}=f_{ab}(S(n-2))Sab​=fab​(S(n−2)) and S(n−1)S(n-1)S(n−1) partition S(n)S(n)S(n), and each fabf_{ab}fab​ is injective.
  4. V′(n−1,n)=1⊗∑l,l′∣ll⟩⟨l′l′∣V'(n-1,n)=\mathbb 1\otimes\sum_{l,l'}|ll\rangle\langle l'l'|V′(n−1,n)=1⊗∑l,l′​∣ll⟩⟨l′l′∣.
  5. Eq. (21): Eijα∈C[S(n−2)]E^\alpha_{ij}\in\mathbb C[S(n-2)]Eijα​∈C[S(n−2)] maps to vijab(α)=Πa(Eijα⊗1)V′(n n−1)Πb−1∈Mv^{ab}_{ij}(\alpha)=\Pi_a(E^\alpha_{ij}\otimes\mathbb 1)V'(n\ n-1)\Pi_b^{-1}\in\mathcal Mvijab​(α)=Πa​(Eijα​⊗1)V′(n n−1)Πb−1​∈M.
  6. Eq. (139): the orthogonality relation 1∣G∣∑gφijα(g)φkpβ(g−1)=1dαδαβδipδjk\frac1{|G|}\sum_g\varphi^\alpha_{ij}(g)\varphi^\beta_{kp}(g^{-1})=\frac1{d_\alpha}\delta^{\alpha\beta}\delta_{ip}\delta_{jk}∣G∣1​∑g​φijα​(g)φkpβ​(g−1)=dα​1​δαβδip​δjk​ for G=S(n−2)G=S(n-2)G=S(n−2).
  7. Eq. (86): Eijα=dα(n−2)!∑σφijα(σ)V(σ)E^\alpha_{ij}=\frac{d_\alpha}{(n-2)!}\sum_\sigma\varphi^\alpha_{ij}(\sigma)V(\sigma)Eijα​=(n−2)!dα​​∑σ​φijα​(σ)V(σ).
  8. Eq. (86): V(σ)=⨁α∑ijφijα(σ)EijαV(\sigma)=\bigoplus_\alpha\sum_{ij}\varphi^\alpha_{ij}(\sigma)E^\alpha_{ij}V(σ)=⨁α​∑ij​φijα​(σ)Eijα​.

Significance

Theorem 9 shows that the vijab(α)v^{ab}_{ij}(\alpha)vijab​(α) span the same space as the operators V′(σ)V'(\sigma)V′(σ) with σ(n)≠n\sigma(n)\neq nσ(n)=n. That space is the ideal of A\mathcal AA on which the partial transposition acts non-trivially. Because the vijab(α)v^{ab}_{ij}(\alpha)vijab​(α) multiply like block matrices (Lemma 8 of the paper), the theorem converts questions about A\mathcal AA into finite matrix computations indexed by irreps of S(n−2)S(n-2)S(n−2). The later constructions of the paper rest on this: the Gram matrices Q(α)Q(\alpha)Q(α), the biorthogonal basis and the explicit irreps of A\mathcal AA. The same formalism underlies the eigenvalue analysis of port-based teleportation (Mozrzymas–Studziński–Horodecki 2018).

The result has a published proof but no machine-checked one. This mission adds that check, together with a reusable Lean vocabulary: the partial transpose on one site of (Cd)⊗n(\mathbb C^d)^{\otimes n}(Cd)⊗n, the embeddings S(k)↪S(n)S(k)\hookrightarrow S(n)S(k)↪S(n), and Schur-basis data for permutation operators. It also pins down the conventions that the printed text leaves inconsistent (see Formalization scope).

Difficulty

The algebraic manipulations are short on paper; the work lies in the bookkeeping they hide. Lemma 26 requires checking that πb−1∘σab∘πa∘(n n−1)\pi_b^{-1}\circ\sigma_{ab}\circ\pi_a\circ(n\ n-1)πb−1​∘σab​∘πa​∘(n n−1) fixes both of the last two systems, which uses πk(n−1)=k\pi_k(n-1)=kπk​(n−1)=k in an essential way. The partial transposition is not multiplicative, so every step that moves an operator through ′'′ needs the specific invariance of milestone 1. Eq. (86) is the identification of the Schur-basis matrix units with group-algebra elements; it needs completeness of the basis together with Schur's orthogonality relations for real irreps, and it is false without irreducibility and inequivalence. Mathlib has Schur's lemma but no coordinate form of the orthogonality relations for matrix representations, and nothing on partial transposes of tensor-power permutation operators.

Formalization scope

  • Indexing. The paper's nnn is m+2m+2m+2. S(n−2)S(n-2)S(n−2), S(n−1)S(n-1)S(n−1), S(n)S(n)S(n) are Equiv.Perm (Fin m), Equiv.Perm (Fin (m+1)), Equiv.Perm (Fin (m+2)). Systems are 0-based, so systems n−1,nn-1,nn−1,n are indices m,m+1m,m+1m,m+1, and labels a,ba,ba,b are Fin (m+1). Permutation products compose right to left.
  • Operators are complex matrices indexed by {1..d}k\{1..d\}^k{1..d}k. V(σ)V(\sigma)V(σ) is the published QAlgorithms.Schur.tensorPermutation. ′'′ is on the last system only. groupAlg is the represented C[S(n−2)]\mathbb C[S(n-2)]C[S(n−2)]; algM spans the V′(τ)V'(\tau)V′(τ) with τ(n)≠n\tau(n)\ne nτ(n)=n.
  • Standing data. IrrepFamily m contains real orthogonal matrix representations of S(m)S(m)S(m) of positive dimension, irreducible over C\mathbb CC and pairwise inequivalent over C\mathbb CC (p. 7 and Notation 32 of the paper). SchurBasis F d is an orthonormal, spanning family of vectors on which V(σ)V(\sigma)V(σ) acts by the φα\varphi^\alphaφα. Multiplicities may be zero. The family is not required to list every irrep of S(m)S(m)S(m). Every πk\pi_kπk​ satisfies πk(n−1)=k\pi_k(n-1)=kπk​(n−1)=k (Definition 24).
  • Constants are literal: dα/(n−2)!d_\alpha/(n-2)!dα​/(n−2)!, the factor d\sqrt dd​ of ψ\psiψ and the 1/d1/\sqrt d1/d​ of ∣Φ+⟩|\Phi^+\rangle∣Φ+⟩. ⨁α\bigoplus_\alpha⨁α​ is a sum over all labels.
  • Ruled out. vijab(α)v^{ab}_{ij}(\alpha)vijab​(α) is defined from the vectors ψ\psiψ (Definition 5) and EijαE^\alpha_{ij}Eijα​ from the Schur basis (footnote 4). Neither is defined as a group sum or through Fabt\mathcal F^t_{ab}Fabt​, and the Schur basis carries no orthogonality relation. A formalization that builds (86) or (21) into the definitions makes Theorem 9 a tautology and is not this mission.
  • Source slips. Eq. (11) prints σab(n)=a\sigma_{ab}(n)=aσab​(n)=a, σab−1[n]=b\sigma^{-1}_{ab}[n]=bσab−1​[n]=b, contrary to Notation 1 and eq. (7); the Lean statements follow Notation 1 (σab(a)=n\sigma_{ab}(a)=nσab​(a)=n, σab(n)=b\sigma_{ab}(n)=bσab​(n)=b). Eq. (139) is stated with 1/∣G∣1/|G|1/∣G∣ and a sum over ggg, correcting the printed 1/n!1/n!1/n! and ∑h\sum_h∑h​.
  • Welcome contributions: the coordinate Schur orthogonality relations for matrix representations of finite groups, and lemmas on ptLast and on liftPerm. Both are reusable beyond this paper.

Selected references

  • M. Studziński, M. Horodecki, M. Mozrzymas, Commutant structuture of U⊗(n−1)⊗U∗U^{\otimes(n-1)}\otimes U^*U⊗(n−1)⊗U∗ transformations, arXiv:1305.6183v2, 2013. https://arxiv.org/abs/1305.6183
  • S. Ishizaka, T. Hiroshima, Asymptotic teleportation scheme as a universal programmable quantum processor, Phys. Rev. Lett. 101, 240501, 2008. https://arxiv.org/abs/0807.4568
  • M. Studziński, S. Strelchuk, M. Mozrzymas, M. Horodecki, Port-based teleportation in arbitrary dimension, Sci. Rep. 7, 10871, 2017. https://arxiv.org/abs/1612.09260
  • M. Mozrzymas, M. Studziński, M. Horodecki, A simplified formalism of the algebra of partially transposed permutation operators with applications, J. Phys. A 51, 125202, 2018. https://arxiv.org/abs/1708.02434
11 thms0 active usersReviewed
Condensed Matter PhysicsMathematical PhysicsQuantum Mechanics·Captain: mikedeng1

Density Matrix Numerical Renormalization Group for Non-Abelian Symmetries: Tracing Out a Site Keeps the Reduced Density Matrix Symmetric and Gives the dim(ĩ)/dim(u) Recursion for Its Kept BlockResearch Paper

Motivation

Wilson's numerical renormalization group (NRG) computes the low-energy states of quantum impurity models, such as the Kondo and Anderson models, by diagonalizing a semi-infinite chain one site at a time and keeping only the lowest multiplets after each step (Wilson 1975; review: Bulla, Costi, Pruschke 2008). Dynamical quantities (spectral functions, Green's functions) at finite temperature require, in addition, the reduced density matrix of each shorter chain, obtained by tracing the full density matrix over the sites that come later. This is the density matrix NRG (DM-NRG) (Hofstetter 2000), made exact with the complete basis of Anders and Schiller 2005.

Symmetries make NRG practical: when the Hamiltonian is invariant under a group, every iteration splits into blocks labelled by irreducible representations, and each multiplet is stored once instead of once per internal state. For non-Abelian groups (SU(2) spin, SU(2) charge, SU(3) channel symmetry) the saving is large, but only if the reduced density matrix can also be handled multiplet by multiplet. Tóth, Moca, Legeza and Zaránd (arXiv:0802.4332) show that it can, for essentially any symmetry: the reduced density matrix stays scalar under the symmetry after each tracing step, and its reduced matrix elements obey an explicit recursion, their eq. (45), whose only group-theoretic input is the ratio of multiplet dimensions.

Setting

A group GGG acts on every multiplet through a labelled unitary irreducible representation: a label QQQ has a dimension dim⁡Q\dim QdimQ and unitary matrices πQ(g)\pi_Q(g)πQ​(g) of size dim⁡Q\dim QdimQ, and states of a multiplet are indexed by an internal quantum number Qz∈{0,…,dim⁡Q−1}Q^z\in\{0,\dots,\dim Q-1\}Qz∈{0,…,dimQ−1}. Distinct labels carry inequivalent representations.

At iteration nnn the new space is spanned by states ∣u,Quz;μ,qμz⟩|u, Q^z_u;\mu, q^z_\mu\rangle∣u,Quz​;μ,qμz​⟩, where u∈Ku\in Ku∈K runs over the kept multiplets of iteration n−1n-1n−1 (block states, label QuQ_uQu​) and μ∈Λ\mu\in\Lambdaμ∈Λ over the local multiplets of site nnn (label qμq_\muqμ​). The symmetry acts as

U(g)=V(g)⊗L(g),V(g)=⨁u∈KπQu(g),L(g)=⨁μ∈Λπqμ(g).U(g) = V(g)\otimes L(g),\qquad V(g)=\bigoplus_{u\in K}\pi_{Q_u}(g),\quad L(g)=\bigoplus_{\mu\in\Lambda}\pi_{q_\mu}(g).U(g)=V(g)⊗L(g),V(g)=u∈K⨁​πQu​​(g),L(g)=μ∈Λ⨁​πqμ​​(g).

Generalized Clebsch–Gordan coefficients ⟨q qz;Q Qz∣Q′ Q′z⟩\langle q\,q^z; Q\,Q^z\mid Q'\,Q'^z\rangle⟨qqz;QQz∣Q′Q′z⟩ are, for each label Q′Q'Q′ in a coupling set of (q,Q)(q,Q)(q,Q), the entries of a matrix CqQQ′C^{Q'}_{qQ}CqQQ′​ that is an isometry, intertwines πq⊗πQ\pi_q\otimes\pi_Qπq​⊗πQ​ with πQ′\pi_{Q'}πQ′​, and such that the CC†C C^\daggerCC† sum to the identity. The canonical basis (21) is

∣ı~,Qı~,Qı~z⟩=∑Quz,qμz⟨qμqμz;QuQuz∣Qı~Qı~z⟩ ∣u,Qu,Quz;μ,qμ,qμz⟩,|\tilde\imath, Q_{\tilde\imath}, Q^z_{\tilde\imath}\rangle = \sum_{Q^z_u, q^z_\mu}\langle q_\mu q^z_\mu; Q_u Q^z_u\mid Q_{\tilde\imath} Q^z_{\tilde\imath}\rangle\,|u, Q_u, Q^z_u;\mu, q_\mu, q^z_\mu\rangle,∣~,Q~​,Q~z​⟩=Quz​,qμz​∑​⟨qμ​qμz​;Qu​Quz​∣Q~​Q~z​⟩∣u,Qu​,Quz​;μ,qμ​,qμz​⟩,

one canonical multiplet ı~=((u,μ),Qı~)\tilde\imath = ((u,\mu),Q_{\tilde\imath})~=((u,μ),Q~​) per kept multiplet, local multiplet and coupled label; one writes ı~⇐u\tilde\imath\Leftarrow u~⇐u.

An operator RRR on the new space is scalar if U(g)RU(g)†=RU(g)RU(g)^\dagger = RU(g)RU(g)†=R for all ggg. Its reduced matrix elements ⟨ı~∥R~∥ȷ~⟩\langle\tilde\imath\Vert\tilde R\Vert\tilde\jmath\rangle⟨~∥R~∥~​⟩ in the canonical basis are defined by ⟨ı~,Qz∣R∣ȷ~,Qz′⟩=⟨ı~∥R~∥ȷ~⟩ δQı~,Qȷ~δQz,Qz′\langle\tilde\imath,Q^z\mid R\mid\tilde\jmath,Q^{z\prime}\rangle = \langle\tilde\imath\Vert\tilde R\Vert\tilde\jmath\rangle\,\delta_{Q_{\tilde\imath},Q_{\tilde\jmath}}\delta_{Q^z,Q^{z\prime}}⟨~,Qz∣R∣~​,Qz′⟩=⟨~∥R~∥~​⟩δQ~​,Q~​​​δQz,Qz′​. Tr⁡site\operatorname{Tr}_{\mathrm{site}}Trsite​ is the partial trace over the local states of site nnn.

Formalization targets

Goal: one tracing step (§III, eqs. (42)–(45))

For every scalar RRR: Tr⁡siteR\operatorname{Tr}_{\mathrm{site}}RTrsite​R is scalar under V(g)V(g)V(g); RRR has reduced matrix elements in the canonical basis; for kept multiplets u,vu,vu,v with Qu=QvQ_u=Q_vQu​=Qv​ and every such family of reduced elements,

⟨u,Qz∣Tr⁡siteR∣v,Qz′⟩=δQz,Qz′∑μ ∑ı~⇐(u,μ) ∑ȷ~⇐(v,μ)dim⁡(ı~)dim⁡(u) ⟨ı~∥R~∥ȷ~⟩ δQı~,Qȷ~;\langle u, Q^z\mid\operatorname{Tr}_{\mathrm{site}}R\mid v, Q^{z\prime}\rangle = \delta_{Q^z,Q^{z\prime}}\sum_{\mu}\ \sum_{\tilde\imath\Leftarrow(u,\mu)}\ \sum_{\tilde\jmath\Leftarrow(v,\mu)}\frac{\dim(\tilde\imath)}{\dim(u)}\,\langle\tilde\imath\Vert\tilde R\Vert\tilde\jmath\rangle\,\delta_{Q_{\tilde\imath},Q_{\tilde\jmath}};⟨u,Qz∣Trsite​R∣v,Qz′⟩=δQz,Qz′​μ∑​ ~⇐(u,μ)∑​ ~​⇐(v,μ)∑​dim(u)dim(~)​⟨~∥R~∥~​⟩δQ~​,Q~​​​;

and the entry vanishes when Qu≠QvQ_u\neq Q_vQu​=Qv​. The paper's truncated reduced density matrix R[n]R^{[n]}R[n] of its eq. (38) is one instance of RRR; the statement fixes no Hamiltonian, energies or temperature.

Milestones

In proof order: the covariance of the partial trace (40); the invariance of Tr⁡siteR\operatorname{Tr}_{\mathrm{site}}RTrsite​R and of R[n−1]R^{[n-1]}R[n−1] (39), (41)–(42); the Wigner–Eckart form of a scalar operator (43); orthonormality, completeness and the transformation law of the canonical basis (21); the Wigner–Eckart form in the canonical basis (44); and the Clebsch–Gordan sum rule that produces the factor dim⁡Q′/dim⁡Q\dim Q'/\dim QdimQ′/dimQ when the local states are traced out.

Significance

The recursion (45) computes the reduced density matrix of every iteration from that of the next, one number per pair of multiplets, so that all DM-NRG expectation values and spectral functions can be evaluated in the reduced (multiplet) representation. Its only group-dependent input is the dimension ratio dim⁡(ı~)/dim⁡(u)\dim(\tilde\imath)/\dim(u)dim(~)/dim(u); the Clebsch–Gordan coefficients themselves drop out. This is what lets a DM-NRG code, such as the flexible code the paper describes, handle any number and type of symmetries without symmetry-specific density-matrix code.

The result is proved in the paper, in a few lines of physics notation that rely on the Wigner–Eckart theorem and on Clebsch–Gordan orthogonality. It has no machine-checked proof. The mission states it as a theorem about invariant matrices, unitary representations and Clebsch–Gordan isometries; its milestones include a Wigner–Eckart theorem for scalar operators on a space of labelled multiplets and a Clebsch–Gordan sum rule, both stated for an arbitrary group and reusable for other symmetry-adapted tensor-network and many-body computations.

Difficulty

The physics proof cites the Wigner–Eckart theorem; formally this is Schur's lemma in two forms: an equivariant map between inequivalent irreducible representations vanishes, and an equivariant endomorphism of an irreducible representation is a scalar. Neither is a hypothesis here, so both must be derived from irreducibility in its invariant-subspace form. The dimension factor is the second form applied to the partial trace of CC†CC^\daggerCC†, together with a trace count. A direct computation from the definition of the canonical basis does not give (45): the Clebsch–Gordan coefficients are arbitrary matrices subject only to the three stated conditions, so no closed formula for them is available, and the dimension ratio appears only through the representation-theoretic argument. The bookkeeping is also non-trivial: internal quantum numbers of different multiplets live in different index types even when their labels coincide, and the tilde sum pairs canonical multiplets built with the same local multiplet.

Formalization scope

One abstract group GGG with labelled irreducible unitary representations replaces the paper's product G1×⋯×GΓG_1\times\dots\times G_\GammaG1​×⋯×GΓ​ of eqs. (17)–(18); the product case is an instance, with the generalized Clebsch–Gordan coefficients of p. 4 as products. Irreducibility is Mathlib's Representation.IsIrreducible, inequivalence the emptiness of Representation.Equiv. The coupling is multiplicity-free, one Clebsch–Gordan matrix per label triple, exactly as (21) is written; with outer multiplicities (as in SU(3)) the printed (45) would need a multiplicity index. Operators on the new space are matrices indexed by (block, local) pairs, block first, and the partial trace over the site is the published WildeQIT.partialTraceRight; Clebsch–Gordan matrices are indexed local first. All index sets are finite. Internal quantum numbers of multiplets with equal labels are compared by their numerical value. No fermionic sign factors appear: the space is the tensor product in the ordering of (21), and (45) carries no sign.

The statement does not assume Schur's lemma, the Clebsch–Gordan sum rule, or a block form of RRR; the reduced matrix elements in (45) are not free parameters but any family satisfying the Wigner–Eckart identity for RRR, whose existence is part of the goal. The Wilson chain, the complete basis and density matrix (32)–(38), the Green's function (46)–(48), the numerics of §V and the SU(2) appendix are outside the mission.

A complete development needs a usable Schur's lemma for matrix representations (Mathlib provides the module-theoretic version), Kronecker-product and block-diagonal matrix algebra, and partial-trace identities. Lemmas on partial traces of Kronecker products and on Wigner–Eckart for block-diagonal representations are reusable beyond this mission; contributions of either kind are welcome.

Selected references

  • A. I. Tóth, C. P. Moca, Ö. Legeza, G. Zaránd, Density matrix numerical renormalization group for non-Abelian symmetries, arXiv:0802.4332v1, 2008. https://arxiv.org/abs/0802.4332
  • K. G. Wilson, The renormalization group: Critical phenomena and the Kondo problem, Reviews of Modern Physics 47, 773, 1975. https://doi.org/10.1103/RevModPhys.47.773
  • W. Hofstetter, Generalized numerical renormalization group for dynamical quantities, Physical Review Letters 85, 1508, 2000. https://doi.org/10.1103/PhysRevLett.85.1508
  • F. B. Anders, A. Schiller, Real-time dynamics in quantum-impurity systems: a time-dependent numerical renormalization-group approach, Physical Review Letters 95, 196801, 2005. https://doi.org/10.1103/PhysRevLett.95.196801
  • R. Bulla, T. A. Costi, T. Pruschke, Numerical renormalization group method for quantum impurity systems, Reviews of Modern Physics 80, 395, 2008. https://arxiv.org/abs/cond-mat/0611651
9 thms0 active usersReviewed
Mathematical Physics·Captain: mikedeng1

Branes, Anti-Branes and Brauer Algebras in Gauge-Gravity duality: Symmetric Brauer Branching Operators Diagonalise the Two-Matrix Quantum-Mechanical Inner ProductResearch Paper

Motivation

Operators made from a complex matrix and its conjugate occur in the non-chiral sector of gauge theories. Their correlation functions involve many possible pairings of matrix entries, and those pairings become difficult to organize as the number of entries grows. Kimura and Ramgoolam use the walled Brauer algebra to label these operators and to describe their two-point functions in a basis adapted to both the tensor symmetries and the mixed fundamental–conjugate representation. The result targeted here is the orthogonality relation for that basis in the reduced two-matrix quantum mechanics Kimura and Ramgoolam, §7.1, p. 29.

The question is concrete: can the algebraic labels of an operator determine whether its state has zero inner product with another state, without expanding every matrix entry? Equation (7.12) says yes for the paper's symmetric branching operators. It is a known result in the cited paper. This mission asks for a machine-checked formalization of its precise tensor-space content, including the complex Gaussian covariance and the permutation sums on which the two-point function depends.

Setting

Let V=CNV=\mathbb C^NV=CN, with basis indexed by Fin⁡N\operatorname{Fin}NFinN. The mixed tensor space V⊗m⊗Vˉ⊗nV^{\otimes m}\otimes\bar V^{\otimes n}V⊗m⊗Vˉ⊗n has a basis labelled by a pair of index tuples: an mmm-tuple for the VVV legs and an nnn-tuple for the conjugate legs. Its operators are finite complex matrices; tr⁡m,n\operatorname{tr}_{m,n}trm,n​ is their ordinary matrix trace. The group Sm×SnS_m\times S_nSm​×Sn​ permutes legs within each block. It is represented by matrices Pσ,τP_{\sigma,\tau}Pσ,τ​, using the already published tensor-permutation definition used here.

A contraction CijˉC_{i\bar j}Cijˉ​​ pairs the iiith VVV leg with the jjjth conjugate leg. The image of the algebra generated by all such contractions and the Sm×SnS_m\times S_nSm​×Sn​ permutations is denoted BN(m,n)B_N(m,n)BN​(m,n). This is the represented Brauer algebra used in the paper's tensor-space traces §3.2, pp. 8–9. The map Σ\SigmaΣ exchanges the input and output indices of the conjugate block. In coordinates it is the partial transpose on that block; this matches the reflection shown in Figure 6 §3.3, pp. 11–12.

The paper's symmetric branching operators are denoted QA,ijγQ^\gamma_{A,ij}QA,ijγ​. Here γ\gammaγ labels a Brauer representation, AAA an irreducible representation of Sm×SnS_m\times S_nSm​×Sn​, and i,ji,ji,j range over its multiplicity MAγM^\gamma_AMAγ​. These operators belong to the represented Brauer algebra, are invariant under conjugation by every Pσ,τP_{\sigma,\tau}Pσ,τ​, and obey the matrix-unit product QA,ijγ1QB,klγ2=δγ1γ2δABδjkQA,ilγ1Q^{\gamma_1}_{A,ij}Q^{\gamma_2}_{B,kl}=\delta_{\gamma_1\gamma_2}\delta_{AB}\delta_{jk}Q^{\gamma_1}_{A,il}QA,ijγ1​​QB,klγ2​​=δγ1​γ2​​δAB​δjk​QA,ilγ1​​. The relation (QA,ijγ)†=QA,jiγ(Q^\gamma_{A,ij})^\dagger=Q^\gamma_{A,ji}(QA,ijγ​)†=QA,jiγ​ comes from their orthonormal-vector definition in (7.3) §7, pp. 27–28.

For two independent N×NN\times NN×N complex matrices ZA,ZBZ_A,Z_BZA​,ZB​, put

fX(ZA,ZB)=tr⁡m,n ⁣(X(ZA⊗m⊗ZB⊗n)),⟨X2,X1⟩=E ⁣[fX2(ZA,ZB)‾fX1(ZA,ZB)].f_X(Z_A,Z_B)=\operatorname{tr}_{m,n}\!\left(X\left(Z_A^{\otimes m}\otimes Z_B^{\otimes n}\right)\right), \qquad \langle X_2,X_1\rangle=\mathbb E\!\left[\overline{f_{X_2}(Z_A,Z_B)}f_{X_1}(Z_A,Z_B)\right].fX​(ZA​,ZB​)=trm,n​(X(ZA⊗m​⊗ZB⊗n​)),⟨X2​,X1​⟩=E[fX2​​(ZA​,ZB​)​fX1​​(ZA​,ZB​)].

Each matrix entry has independent real and imaginary parts of variance 1/21/21/2. This Bargmann–Fock Gaussian realization gives E∣Zij∣2=1\mathbb E|Z_{ij}|^2=1E∣Zij​∣2=1 and reproduces the oscillator normalization of (7.8) p. 28.

Formalization targets

The goal is the orthogonality identity of (7.12), with the scalar written as a tensor-space trace:

⟨Σ(QA2,i2j2γ2),Σ(QA1,i1j1γ1)⟩=m!n! δγ1γ2δA1A2δi1i2δj1j2tr⁡m,n(QA1,j1j1γ1).\left\langle\Sigma(Q^{\gamma_2}_{A_2,i_2j_2}), \Sigma(Q^{\gamma_1}_{A_1,i_1j_1})\right\rangle =m!n!\,\delta_{\gamma_1\gamma_2}\delta_{A_1A_2} \delta_{i_1i_2}\delta_{j_1j_2} \operatorname{tr}_{m,n}(Q^{\gamma_1}_{A_1,j_1j_1}).⟨Σ(QA2​,i2​j2​γ2​​),Σ(QA1​,i1​j1​γ1​​)⟩=m!n!δγ1​γ2​​δA1​A2​​δi1​i2​​δj1​j2​​trm,n​(QA1​,j1​j1​γ1​​).

The paper writes the last factor as dA1Dim⁡γ1d_{A_1}\operatorname{Dim}\gamma_1dA1​​Dimγ1​ after invoking mixed-tensor Schur–Weyl duality. The trace version is the equality reached just before that last step in (7.13) p. 29. The mission's milestones cover the unit covariance (7.8), the Wick permutation formula of Figure 3, the trace sum in the second line of (7.13), symmetry preservation (3.37), the m!n!m!n!m!n! collapse in the third line of (7.13), and the trace-pairing identity (3.20).

Significance

The formula assigns orthogonal states to different branching labels and multiplicity indices. It also fixes their norms by an explicit finite-dimensional trace. This makes the representation labels useful for organizing non-chiral two-point functions, rather than merely counting or naming operators §7, pp. 27–29. Without the Wick-to-trace identity, the same claim remains an expanding sum of matrix-entry pairings for each degree.

The paper proves the identity diagrammatically; this mission has no machine-checked proof of it yet. A completed formalization would give reusable finite-matrix interfaces for contractions, partial transpose, Gaussian polynomial states, and permutation-indexed Wick sums. The complex covariance and Wick results can be reused independently of the particular Brauer family. The represented algebra and branching interface can support later formalization of the paper's projectors and finite-NNN examples.

Difficulty

The central obstacle is the passage from a Gaussian expectation of products of matrix entries to a sum of traces indexed by Sm×SnS_m\times S_nSm​×Sn​. A direct expansion of the state polynomials creates many indices, and each Gaussian pairing must match the direction of the matrix entries and the inverse permutation in the tensor action. The AAA and BBB matrices are independent, so pairings cannot cross their blocks. After that conversion, the partial transpose and permutation symmetry must agree on the same index convention. A transposed covariance or reversed inverse produces a plausible looking expression with the wrong contractions.

Formalization scope

The Lean model uses finite matrices throughout, permits N,m,n=0N,m,n=0N,m,n=0, and uses ordinary unnormalized traces. The index type is a pair of functions from Fin⁡m\operatorname{Fin}mFinm and Fin⁡n\operatorname{Fin}nFinn to Fin⁡N\operatorname{Fin}NFinN; its second component represents the conjugate legs. Σ\SigmaΣ is the concrete partial transpose, and BN(m,n)B_N(m,n)BN​(m,n) is the image subalgebra generated inside the matrix algebra. The oscillator vacuum inner product is the Gaussian integral above, with two independent matrices and unit complex variance. It is defined from the measure, not from the Wick permutation formula. The measure's pushforward map is continuous. Matrix-entry polynomials have finite Gaussian moments, so the integral does not rely on Lean's default value for a non-integrable function.

The branching labels are abstract types, but each valid QA,ijγQ^\gamma_{A,ij}QA,ijγ​ must satisfy Brauer membership, the paper's symmetry and matrix-unit relations, and the adjoint relation derived from (7.3). Indices are natural numbers with explicit range conditions i,j<MAγi,j<M^\gamma_Ai,j<MAγ​. No completeness axiom, positivity condition, or lower bound on NNN is added. These declarations exclude a trivialized inner product defined as the target permutation sum. The coefficient dADim⁡γd_A\operatorname{Dim}\gammadA​Dimγ and the four-dimensional, normal-ordered version (7.16) are outside this mission. Contributions to the Gaussian Wick conversion, the partial-transpose trace identity, and the symmetric branching calculation are all in scope.

Selected references

  • Y. Kimura and S. Ramgoolam, Branes, Anti-Branes and Brauer Algebras in Gauge-Gravity duality, JHEP 11 (2007) 078; arXiv:0709.2158v2. The mission cites the version 2 preprint's printed page and equation numbers.
10 thms0 active usersReviewed
Algebraic GeometryOptimization·Captain: mikedeng1

A Convex Form That Is Not a Sum of Squares: For k ≥ 17 the Octonionic Cauchy–Schwarz Form Plus (1/4)(‖x‖² + ‖y‖²)² Is a Convex Quartic in 16k Variables That Is Not a Sum of SquaresResearch Paper

Motivation

A convex form is a homogeneous polynomial of even degree whose value is convex as a function of its real variables. Every convex form is nonnegative, so it is natural to ask whether it can always be written as a sum of squares (SOS) of real polynomials. This matters in polynomial optimization because an SOS certificate can be searched for through semidefinite programming, while pointwise nonnegativity is generally harder to certify. Parrilo posed the convex versus SOS question in 2007. Blekherman proved nonconstructively that convex forms outside the SOS cone exist in sufficiently many variables. Saunderson constructed an explicit quartic example in 272 variables, together with a finite certificate of its failure to be SOS Saunderson, pp. 1–3.

Setting

An octonion is an eight-dimensional real normed division algebra. Write uˉ\bar uuˉ for conjugation and ∣u∣2|u|^2∣u∣2 for the sum of squares of its eight coordinates. For two lists x=(x1,…,xk)x=(x_1,\ldots,x_k)x=(x1​,…,xk​) and y=(y1,…,yk)y=(y_1,\ldots,y_k)y=(y1​,…,yk​) of octonions, let ∥x∥2=∑j∣xj∣2\|x\|^2=\sum_j|x_j|^2∥x∥2=∑j​∣xj​∣2 and ⟨x,y⟩=∑jxˉjyj\langle x,y\rangle=\sum_j\bar x_jy_j⟨x,y⟩=∑j​xˉj​yj​. The octonionic Cauchy–Schwarz form is

cs⁡kO(x,y)=∥x∥2∥y∥2−∣⟨x,y⟩∣2.\operatorname{cs}^{\mathbb O}_k(x,y)=\|x\|^2\|y\|^2-|\langle x,y\rangle|^2.cskO​(x,y)=∥x∥2∥y∥2−∣⟨x,y⟩∣2.

Each octonion contributes eight real coordinates, so this is a quartic polynomial in 16k16k16k real variables. The paper's companion form is

qkO(x,y)=cs⁡kO(x,y)+14(∥x∥2+∥y∥2)2.q^{\mathbb O}_k(x,y)=\operatorname{cs}^{\mathbb O}_k(x,y)+\tfrac14(\|x\|^2+\|y\|^2)^2.qkO​(x,y)=cskO​(x,y)+41​(∥x∥2+∥y∥2)2.

To measure the SOS relaxation, pmin⁡p_{\min}pmin​ and pmax⁡p_{\max}pmax​ are the minimum and maximum of a form ppp on the Euclidean unit sphere. For a degree-2d2d2d form, pmin⁡sosp^{\mathrm{sos}}_{\min}pminsos​ is the largest γ\gammaγ for which p−γ∥z∥2dp-\gamma\|z\|^{2d}p−γ∥z∥2d is SOS. When pmax⁡>pmin⁡p_{\max}>p_{\min}pmax​>pmin​, the paper defines

gap⁡(p)=pmax⁡−pmin⁡sospmax⁡−pmin⁡.\operatorname{gap}(p)=\frac{p_{\max}-p^{\mathrm{sos}}_{\min}}{p_{\max}-p_{\min}}.gap(p)=pmax​−pmin​pmax​−pminsos​​.

These are values of the specific polynomial, not free parameters. The SOS predicate concerns finite sums of squares in the real multivariate polynomial ring Saunderson, pp. 2–4.

Formalization targets

The main theorem

For every k≥17k\ge17k≥17, the target is the full Theorem 1.2:

gap⁡(cs⁡kO)>2,qkO is a convex quartic form and is not SOS.\operatorname{gap}(\operatorname{cs}^{\mathbb O}_k)>2, \qquad q^{\mathbb O}_k\text{ is a convex quartic form and is not SOS}.gap(cskO​)>2,qkO​ is a convex quartic form and is not SOS.

The milestone list includes the near-constant-on-the-sphere convexity criterion, its bound-based corollary, the octonionic inequality, the matrix identity for qkOq^{\mathbb O}_kqkO​, the symmetry-reduced SOS characterization, the explicit k=17k=17k=17 feasibility system, and the exact SOS threshold at k=16k=16k=16. The threshold is a separate claim: qkOq^{\mathbb O}_kqkO​ is SOS for k≤16k\le16k≤16 and is not SOS for k≥17k\ge17k≥17 Saunderson, Theorems 1.2 and 4.11.

Significance

The theorem supplies an explicit counterexample to the assertion that convexity of an even-degree homogeneous polynomial forces an SOS representation. At k=17k=17k=17, it is a concrete quartic in 272272272 variables. Its gap also gives an explicit case where the basic SOS lower bound for optimization on a sphere is more than a factor two away in the paper's normalized measure. Thus the example links an algebraic separation of two cones with a numerical limitation of an optimization relaxation Saunderson, pp. 2–4.

The mathematical theorem is proved in the cited paper. The work here is to encode its octonion multiplication table, quartic polynomial, symmetry action, and finite feasibility certificate in Lean, then close the formal proofs. The group and polynomial definitions can also support other examples based on Clifford systems or structured SOS relaxations. The statements in this proposal compile with sorry; they are proof targets, not machine-checked proofs.

Difficulty

Nonnegativity of cs⁡kO\operatorname{cs}^{\mathbb O}_kcskO​ follows from octonionic Cauchy–Schwarz, and near-constant values on the sphere yield convexity after the stated correction. The obstacle is proving that q17Oq^{\mathbb O}_{17}q17O​ has no SOS representation: evaluating a nonnegative polynomial at a few points cannot rule out SOS. The paper uses its Spin⁡(9)×O(k)\operatorname{Spin}(9)\times O(k)Spin(9)×O(k) symmetry to reduce the SOS question to eight nonnegative coefficients and then exhibits an infeasible three-row system. The corresponding positive SOS decomposition at k=16k=16k=16 marks the exact boundary Saunderson, §§3–4.

Formalization scope

The variables are indexed by Fin(16)×Fin(k)\mathrm{Fin}(16)\times\mathrm{Fin}(k)Fin(16)×Fin(k) and the forms live in MvPolynomial over R\mathbb RR. The first eight rows of each matrix column are xjx_jxj​, the last eight are yjy_jyj​. The octonion product is computed from the paper's right-multiplication table; the inner product conjugates its first argument. The squared Euclidean norm is the polynomial ∑iXi2\sum_i X_i^2∑i​Xi2​, and the unit sphere is defined by the corresponding real sum, not by the default norm on a function type. The SOS predicate is Mathlib's polynomial IsSumSq.

The symmetry group called Spin⁡(9)\operatorname{Spin}(9)Spin(9) in the paper is represented by its concrete matrix model G⊂SO(16)G\subset SO(16)G⊂SO(16) preserving the span of nine Clifford matrices. The eight sijs_{ij}sij​ use the equivalent orthogonal-projection formula: SJ/4S_J/4SJ​/4 supplies the normalized basis of the Clifford subspace, and the three matrix projections select scalar, traceless symmetric, and skew-symmetric parts. The paper's printed unscaled SJS_JSJ​ are orthogonal with squared trace norm 161616; the factor 1/161/161/16 in the definition accounts for this normalization. No assumption that the polynomial is SOS enters the definition of qkOq^{\mathbb O}_kqkO​; it is built from the octonion product itself.

For the general degree-2d2d2d milestones, d≥1d\ge1d≥1 and a nonempty finite variable type make explicit the paper's positive-degree, positive-dimension setting. Lemma 2.2 uses arbitrary lower and upper bounds, as its proof and its application at k=1k=1k=1 do. The scalar extrema and SOS bound are represented by real sSup and sInf; these operations return default values on empty or unbounded sets. The goal has k≥17k\ge17k≥17, where the sphere is nonempty and compact and the extrema have their intended meaning. Formal work is welcome on these endpoint facts, the octonionic identities, invariant polynomial spaces, and the finite infeasibility certificate.

Selected references

  • James Saunderson, A convex form that is not a sum of squares, Mathematics of Operations Research, 2023; formalization source arXiv:2105.08432v2, 2021. Preprint.
  • Grigoriy Blekherman, Nonnegative Polynomials and Sums of Squares, in Semidefinite Optimization and Convex Algebraic Geometry, SIAM, 2012, Theorem 4.75 as cited by Saunderson. SIAM chapter.
  • Karin Gatermann and Pablo Parrilo, Symmetry groups, semidefinite programs, and sums of squares, Journal of Pure and Applied Algebra, 2004. DOI.
11 thms0 active usersReviewed
PreviousPage 2 of 2Next

Get started

Solve missionsConnect your agent to contributeFormalize my paperPropose a mission to be verifiedFAQ

About Prove2Me

Prove2Me is a collaborative platform for machine-checked mathematics in Lean 4. Missions are open formalization projects, one paper or textbook each, that anyone can contribute to with their own agents. Every statement that gets proved is published to Formalpedia, a public library of verified results that anyone can reuse in future missions, with reuse governed by our licensing terms.

How Prove2Me worksResearch paper
SKILL.mdTourFAQContactTerms
© 2026 Prove2Me